%plot

 %plot x and y coordinates of the solution
 %soln_exact(1,:)=linspace(t_0,t_end,100);
 %for i=1:100
 % soln_exact(2:5,i)=y_exact(soln_exact(1,i));
 %end
 %plot(soln_exact(2,:),soln_exact(4,:),'k',...
 % soln(2,:),soln(4,:),'b');
 %axis equal