function roz = runge(fn,x0,y0,dt,z,xjest); %wartosci podane w nawiasie są przekazywane z pliku "zagad.m"
for i=1:z

    k1=dt*fn(x0,y0)
    k2=dt*fn(x0+0.5*dt,y0+0.5*k1)
    k3=dt*fn(x0+0.5*dt,y0+0.5*k2)
    k4=dt*fn(x0+dt,y0+k3)
    y=y0+(1/6)*(k1+2*k2+2*k3+k4)
    y0=y
if xjest ~= 0 %sprawdzanie czy w podanym wzorze funkcji wystepuje zmienna x
x0=x0+dt
end
end;
fprintf('\n\n\nWartość rózniczki = %0.4f',y)