[sebastain]

xmax=sqrt(x0^2+(dotx0/omega0)^2)+min([0.5*abs(f0)*Dt/(m*omega0) f0/omega0^2]);
omegad=omega0*sqrt(1-zeta^2)
上面的语句写成的公式
万分感谢