function xdot=tut2_eq(t,x)
global a
xdot(1)=x(2);
xdot(2)=4*exp(-t)*cos(2*t)-a*x(2)-5*x(1);
xdot=xdot';
