[t,y]=ode45(@fode3,[0,5],[3,-5]) plot(t,y(:,1),'r-') hold on % plot(t,y(:,2)) -- % This above the derivative of sol. Don't really need to plot % True solution -- really this time ezplot('-exp(-t)*sin(2*t)+(11/4)*exp(-t)*cos(2*t)+(1/4)*exp(-t)',[0,5])