problem1.m

来自「平时练习时的一些仿真程序」· M 代码 · 共 14 行

x=[0;0;0];               
A=[-8,1,0;-19,0,1;-12,0,0];
B=[0;4;10];
C=[1 0 0];
h=0.05;
t=0.5;
for n=h:h:t
    k1=A*x+B;
    k2=A*[x+h*k1/2]+B;
    k3=A*[x+h*k2/2]+B;
    k4=A*[x+h*k3]+B;
    x=x+h*(k1+2*k2+2*k3+k4)/6;
end
y=C*x