[tout,yout]=eul('quadratic',[0,1],0.1,0.01,50)
plot(tout,yout)