MATLAB程序ea703m a= input(输入分母系数向量a=[a1,a2,…=); n=length(a)-1 YO= input(输入初始条件向量Y=y0,Dyo, D2y0,…= p=roots(a), V=roto(vander(p)); C= VYO dt=input('dt=); tf=input(=) t=0 dt: tf; y=zeros(1, length(t)) for k=1 n y= y+c(k)*exp(p(k)*t) end plot(t, y), gridMATLAB程序ea703.m a=input('输入分母系数向量a=[a1，a2，...]= '); n=length(a)-1; Y0=input(‘输入初始条件向量 Y0=[y0，Dy0， D2y0，...]= '); p=roots(a);V=rot90(vander(p));c= V\Y0'; dt=input('dt='); tf=input('tf= ') t=0:dt:tf; y=zeros(1，length(t)); for k=1:n y= y+c(k)*exp(p(k)*t);end plot(t，y)，grid