function [x,y]=Euler1(x0,y0,h,n,a,b)
x0=1;
y0=1;
h=0.1;
n=10;
a=1;
b=2;
for i=1:n
i
x=1+h*i
if (x>a)&(x<b)
yp=y0+h*(1/x0.^2-y0/x0);
yc=y0+h*(1/x.^2-yp/x);
y=(yp+yc)/2
else
x
y
end
x0=x;
y0=yp;
plot(x,y,'+')
hold on
end