expfj.asv

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

function expfj(A)
syms x t
[n,m]=size(A);
y=eig(A);
z=zeros(1,n);
for i=1:n;
    for j=1:n;
        for k=i-1:-1:1;
            if y(k)==y(i);
               z(k)=z(k)+1;
           end
       end
       c=diff(x^(j+z(k)),x,z(k));
       c=inline(c);
       B(i,j)=c(y(i));
       a=diff(exp(x*t))
       a=inline(a);
       b(i)=a(y(i),t);
   end
end
gauss()