runge kutta cuarto orden

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