function lsg = veuler(f,n,x,h,y)

lsg = y + h*f(n,x+0.5*h,y+0.5*h*f(n,x,y));