% euler_demo

f='-y^2';
x0=1;
y0=1;
h=0.5;
xn=5;
u='1./x';
ee_demo(f,x0,xn,y0,h,u)
pause
ie_demo(f,x0,xn,y0,h,u)
pause
lm_demo(f,x0,xn,y0,h,u)
pause
clf
h=0.25;
ee_demo(f,x0,xn,y0,h,u)
pause
ie_demo(f,x0,xn,y0,h,u)
pause
lm_demo(f,x0,xn,y0,h,u)
pause
clf
h=0.1;
ee_demo(f,x0,xn,y0,h,u)
pause
ie_demo(f,x0,xn,y0,h,u)
pause
lm_demo(f,x0,xn,y0,h,u)
pause
clf
f='-2*y';
x0=0;
y0=1;
h=0.5;
xn=5;
u='exp(-2*x)';
ee_demo(f,x0,xn,y0,h,u)
pause
ie_demo(f,x0,xn,y0,h,u)
pause
lm_demo(f,x0,xn,y0,h,u)
pause
clf
h=0.25;
ee_demo(f,x0,xn,y0,h,u)
pause
ie_demo(f,x0,xn,y0,h,u)
pause
lm_demo(f,x0,xn,y0,h,u)
pause
clf
h=1;
ee_demo(f,x0,xn,y0,h,u)
pause
ie_demo(f,x0,xn,y0,h,u)
pause
lm_demo(f,x0,xn,y0,h,u)
pause
clf
h=1.25;
xn=10;
ee_demo(f,x0,xn,y0,h,u)
pause
ie_demo(f,x0,xn,y0,h,u)
pause
lm_demo(f,x0,xn,y0,h,u)
