[t,y] = euler (function, [t0 tf], y0, dt)
Integration of a given function using Euler's method.
[t,y]=euler('vdp1',[0,10],[2,0],0.1)