METODO DE EULER

clear all
disp('METODO DE EULER')
clc
syms x
syms y
f=inline(input('ingrese la derivada:','s'));
x=input('ingrese el valor de x:');
y=input('ingrese el valor de y:');
h=input('ingrese el valor de h:');
n=input('ingrese numero de iteraciones:');
clc     
disp('x(n)      y(n)      y´(n)     hy´(n)');
for i=1:n
   y1=feval(f,x,y);
   hy1=h*y1;
   fprintf('n%0.1f     %0.4f      %0.4f       %0.4f  ',x,y,y1,hy1);
   y=y+hy1;
   x=x+h;
end