t = 0:0.2:100;
y = 0.4985*exp(-0.6615*t)+(0.00784264*cos(5*t)-0.00126003*sin(5*t))+exp(-0.0693*t).*(-0.506344*cos(0.2088*t)+1.441432*sin(0.2088*t));
plot(t,y);