jia1.asv

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

function X=jia1(A,t)
[n,m]=size(A);
%a=eig(A);
syms a
f=det(A-a*eye(n));
dsolve
Q=eye(n);
syms s
a(1)=round(a(1));
r=exp(a(1)*t);
S=r*Q;
for i=2:n;
    r1=inline(r);
    a(i)=round(a(i));
    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;