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)…

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;

//********************** Start of Code in Segment: 1 ************** :0001.0100 E9FD2E jmp 3000 :0001.0103 …

sx=24, sy=20 dx=16, dy=16 sx=24, sy=20 dx=16, dy=16 8888888888888888888888888888888888 88 8 88…

sx=24, sy=20 dx=16, dy=16 sx=24, sy=20 dx=16, dy=16 8888888888888888888888888888888888 88 8 88…

function dx=plantmodel(t,x,flag,para) dx=zeros(2,1); u=para; dx(1)=x(2); dx(2)=-50*x(2)+400*u;

function dx=PlantModel(t,x,flag,para) dx=zeros(2,1); u=para; dx(1)=x(2); dx(2)=-50*x(2)+400*u;

function dx=PlantModel(t,x,flag,para) dx=zeros(2,1); u=para; dx(1)=x(2); dx(2)=-50*x(2)+400*u;