ex2_4.m

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

syms T;
A=[0 1 0;0 0 1;-24 -26 -9];
B=[0;24;-192];
[ft fmt]=sszhuanyi(A,B,T);
ftd=double(subs(ft,'T',1))
fmtd=double(subs(fmt,'T',1))