clc;clearsyms t n kA=[-7,-16,4;6,13,-2;12,16,1];% 对角化[~,r]=eig(A)exp(A*t)f=symsum(k/2^k,1,n);%前25项和subs(f,n,25)%极限limit(f,n,inf)
