YHFT-DX 源代码

S = dsolve('Dx = y,Dy = -x'); disp([blanks(12),'x',blanks(21),'y']),disp([S.x,S.y]) y = dsolve('x*D2y-2*Dy = x^2','y(1) = 0,y(5)…

S=dsolve('Dx=y,Dy=-x'); disp([blanks(12),'x',blanks(21),'y']),disp([S.x,S.y])

function dx=apolloeq(t,x) mu=1/82.45; mu1=1-mu; r1=sqrt((x(1)+mu)^2+x(3)^2); r2=sqrt((x(1)-mu1)^2+x(3)^2); …

function dx=exc7ffo2(t,x,z) xlag1=z(:,1); xlag2=z(:,2); xlag3=z(:,3); dx=[x(2); x(3); x(4); exp(-t.^2)-4*xlag1(4)-…

function dx=c6ximp(t,x) A=[sin(x(1)) cos(x(2)); -cos(x(2)) sin(x(1))]; B=[1-x(1); -x(2)]; dx=inv(A)*B;

S = dsolve('Dx = y,Dy = -x'); disp([blanks(12),'x',blanks(21),'y']),disp([S.x,S.y]) y = dsolve('x*D2y-2*Dy = x^2','y(1) = 0,y(5)…

S = dsolve('Dx = y,Dy = -x'); disp([blanks(12),'x',blanks(21),'y']),disp([S.x,S.y]) y = dsolve('x*D2y-2*Dy = x^2','y(1) = 0,y(5)…

function dx=c2ximp(t,x) A=[sin(x(1)) cos(x(2)); -cos(x(2)) sin(x(1))]; B=[1-x(1); -x(2)]; dx=inv(A)*B;

function dx=apolloeq(t,x) mu=1/82.45; mu1=1-mu; r1=sqrt((x(1)+mu)^2+x(3)^2); r2=sqrt((x(1)-mu1)^2+x(3)^2); dx=[x(2); 2*x…

function dx=c3ximp(t,x) A=[sin(x(1)) cos(x(2)); -cos(x(2)) sin(x(1))]; B=[1-x(1); -x(2)]; dx=inv(A)*B;