tmp_sys.m

求解常微分方程的工具

function Y=tmp_sys(t,X); 
Y=zeros(4,1); 
x1=X(1); 
x2=X(2); 
x3=X(3); 
x4=X(4); 
Y(1)=-x2-x3; 
Y(2)=x1+0.25*x2+x4; 
Y(3)=3+x3*x1; 
Y(4)=-0.5*x3+0.05*x4;