YHFT-DX 源代码

[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z','t'); x=simple(x) y=simple(y) …

function dx=shier(t,x) dx=zeros(2,1); dx(1)=x(1)*(1-0.1*x(2)); dx(2)=x(2)*(-0.5+0.02*x(1));

function dx=myfun02(x,y) dx=zeros(2,1); dx(1)=y(1)*(1-0.1*y(2)); dx(2)=y(2)*(-0.5+0.02*y(1));

function dx=c6exdde(t,x,z) xlag1=z(:,1); xlag2=z(:,2); dx=[1-3*x(1)-xlag1(2)-0.2*xlag2(1)^3-xlag2(1); x(3); …

function dx=exc7ffo1(t,x,flag,m) dx=[x(2); -m*(x(1)^2-1)*x(2)-x(1)];

function dx=exc7eqa(t,x) r=sqrt(x(1)^2+x(3)^2); dx=[x(2); -x(1)/r^3; x(4); -x(3)/r^3];

function dx=shier(t,x) dx=zeros(2,1); dx(1)=x(1)*(1-0.1*x(2)); dx(2)=x(2)*(-0.5+0.02*x(1));

[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z','t'); x=simple(x) y=simple(y) …

function dx=shier1(t,x) dx=zeros(2,1); dx(1)=x(1)*(0.7-0.1*x(2)); dx(2)=x(2)*(-0.8+0.02*x(1));

function dx=c3xdae2(t,x,lam) dx=[x(2); -lam*x(1); x(4); -lam*x(3)-9.81; x(1)*x(1)+x(3)*x(3)-1; x(1)*x(2)+x(3)*x(4); x(1)*(-l…