plot the exp function

I'd like to draw graph like this with using plot, not ezplot.

I've tried several expression but not worked.

Answer 1

[Birdman]

syms y(t)
eq=diff(y,2)+13*diff(y)+40*y==4;
Dy=diff(y);
Y(t)=dsolve(eq,[y(0)==0 Dy(0)==0])
t=0:0.001:0.9;
plot(t,Y(t))