% Section 4.5 Exercise 26
% y'' +  4*y = 4*cos(2*t)
% with the 
% particular solution:
% yp=t.*sin(2*t)
% and the 
% homogeneous solution:
% yh = c1*cos(2*t) + c2*sin(2*t) 
t = linspace(-1,10,100);
c1=1;c2=0.001;
yh = c1*cos(2*t) + c2*sin(2*t) ;
yp = t.*sin(2*t);
y = yh + yp;

figure(26)
plot(t,yh,'.-.r',t,yp,'-ob',t,y,'-g','Linewidth',2)
legend('homogeneous','particular','solution')
title('Again - notice the transient !!!')
