7楼: Originally posted by zlt11 at 2014-12-08 14:58:29
function df = odefunction(t,f,fp)
%隐式微分方程组
%   记 y(t)=f(1);z(t)=f(2);y'=fp(1);z'=fp(2)
x=@(t) 0.01*sin(2*pi*t);
dx=@(t) 0.01*2*pi*cos(2*pi*t);
alpha=963;c0=53;k0=14;c1=930;k1=5.4;gamma=2 ...