ea622a.m

来自「用MATLAB程序解决高阶问题；通过20多个实例来说明线性代数在后续课和工程中的」· M 代码 · 共 12 行

% 《工程线性代数（MATLAB版）》第六章例题6.22程序ea622a
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 物料混合问题
% 
clear
A=[-1/10,1/40;1/10,-1/10]
for t=0:100;
Y(:,t+1)=expm(A*t)*[60;0];
end
plot([0:100],Y(1,:),[0:100],Y(2,:),':')
grid on
set(gcf,'color','w')