eatfjs2.m

来自「为配套王能超版本数值分析简明教程编写,方便学生教师演示!」· M 代码 · 共 16 行

function X=eAtfjs2(A,t,e)
[n,m]=size(A);
a=eig(A);
Q=eye(n);
syms s
a(1)=fix(a(1)*10^e)/10^e;
r=exp(a(1)*t);
S=r*Q;
for i=2:n;
    r1=inline(r);
    a(i)=fix(a(i)*10^e)/10^e;
    r=int(exp(a(i)*(t-s))*r1(s),s,0,t);
    Q=(A-a(i-1)*eye(n))*Q;
    S=S+r*Q;
end
X=S;