y:=k1*cos(sqrt(2)*t)+k2*sin(sqrt(2)*t)+2/(2-omega^2)*cos(omega*t);yp:=diff(y,t);y0:=subs(t=0,y);yp0:=subs(t=0,yp);solve({y0=0,yp0=0},{k1,k2});y:=subs(k1=2/(omega^2-2),y);y:=subs(k2=0,y);yh:=2*cos(sqrt(2)*t)/(omega^2-2);yp:=2*cos(omega*t)/(2-omega^2);yh1:=subs(omega=.5,yh);yp1:=subs(omega=.5,yp);with(plots);plot1:=plot(yh1,t=0..60,color=BLACK):plot2:=plot(yp1,t=0..60):display(plot1,plot2);y1:=subs(omega=0.5,y);plot(y1,t=0..60);yh2:=subs(omega=1.2,yh);yp2:=subs(omega=1.2,yp);plot3:=plot(yh2,t=0..60,color=BLACK):plot4:=plot(yp2,t=0..60):display(plot3,plot4);y2:=subs(omega=1.2,y);plot(y2,t=0..60);yh3:=subs(omega=1.3,yh);yp3:=subs(omega=1.3,yp);plot5:=plot(yh3,t=0..60,color=BLACK):plot6:=plot(yp3,t=0..60):display(plot5,plot6);y3:=subs(omega=1.3,y);plot(y3,t=0..60);yr:=1/(2*sqrt(2))*t*sin(sqrt(2)*t);plot(yr,t=0..60);