clc
clear
f=inline('x-y')
x0=0;
y0=2;
xx=1;
n=4;
i=1;
h=(xx-x0)/n;
disp('-n-----x-------------y--------pendiente')
while i<=n
x1=x0+h;
m=f(x0,y0);
y1=y0+m*h;
fprintf('\n%2i %10.5f %10.5f %10.5f \n ',i,x1,y1,m)
y0=y1;
x0=x1;
i=i+1;
end
fprintf('\n y(1)= %2f\n ',y1)