[x,y]=forward_euler(@pendulum_ode,[0,20],[0;3.5],800);
plot(x,y(1,:),'k','linewidth',10)
axis('off')
print -dpng large_euler.png

