s10fig24.m

Industrial Mathematics

% Figure 10.24
clear
a=1;
b=2;
%
t=0:0.1:5;
y=1/(a*b) +  exp(-a*t)/(a*(a-b)) + exp(-b*t)/(b*(b-a));
plot(t,y,'k-')