chap5_section53.m

来自「wonderful tutorial of matlab include man」· M 代码 · 共 12 行

A=[0,0,1,0;0,0,0,1;-6,4,-0.2,0.1;2,-2,0.05,-0.05];

T=0.02;
F0=eye(4,4);

for j=1:500
    F=F0+(A^j)*(T^j)/prod(1:j);
    F0=F;
end
F