function dxdt=myproblemm(t,x)
dx1dt=x(2);
dx2dt=x(3);
dx3dt=-3*x(3)-5*x(2)-7*x(1)+6;
dxdt=[dx1dt; dx2dt; dx3dt];