Page 311
A=[-0.0361  0.0124  0.000035; 0.0111 -0.0286 0.0; 0.0039 0 -0.000035];
f=[49.3; 0; 0]; % column vector
t=0:2:366;   % in steps of 2
% make an m-file, At.m, with
%    function Y=At(A,t)
%    Y=t*A;
AinvF=inv(A)*f;
% make up the seq. of sol. vectors
u=[]; s=size(t);
for i=1:s(2)
   u=[u,-AinvF+expm(At(A,t(i)))*AinvF];
end
plot(t,u)

                Page 312
[evect,eval]=eig(A)

leadlim=-AinvF