clear all

wzor=input('Podaj wzor funkcji=');

f=inline(wzor);

a=input('Podaj x0:');

y(1)=input('Podaj y0');

h=input('Podaj krok h:');

b=input('Podaj ostatni argument:');

x=a:h:b;

for k=2:length(x)

K1=h*f(x(k-1),y(k-1));

K2=h*f(x(k-1)+h/2,y(k-1)+1/2*K1);

K3=h*f(x(k-1)+h/2,y(k-1)+1/2*K2);

K4=h*f(x(k-1)+h,y(k-1)+K3);

y(k)=y(k-1)+1/6*(K1+2*K2+2*K3+K4);

end

plot(x,y,'ro')

hold on

disp('Wpisz wzor rozwiazania dokladnego')

dokl=input('y=');

roz=inline(dokl);

y(length(x))

ezplot(roz,[a-h,b+h])

---