ex2_2_2.m

来自「1．编写欧拉前差、后差、梯形公式。 2．编写二阶、三阶龙格库塔法通用程序。 」· M 代码 · 共 11 行

syms f
f(1)='y(2)';
f(2)='y(3)+24';
f(3)='-24*y(1)-26*y(2)-9*y(3)-192';

[y,t]=rungekutta_3(f,[1 0.5 -24]',0,1,0.02)

subplot(2,2,1),stairs(t,y(1,:)),axis([0 1 1 1.17]),title('y')
subplot(2,2,2),stairs(t,y(1,:)),axis([0 1 1 1.17]),title('x1')
subplot(2,2,3),stairs(t,y(2,:)),axis([0 1 -0.3 0.5]),title('x2')
subplot(2,2,4),stairs(t,y(3,:)),axis([0 1 -25.5 -24]),title('x3')