function roz = predic (fn,x0,y0,dt,z,xjest);    
for i=1:z
yg=y0+dt*(fn(x0,y0)) %funkcja y*
y=y0+0.5*dt*(fn(x0,y0)+fn(x0,yg))
y0=y
if xjest ~= 0 %sprawdzanie czy zmienna x wystepuje w podanym wzorze funkcji
x0=x0+dt
end
end;
fprintf('\n\n\nWartość różniczki = %0.4f',y) %zaokragla do 4 miejsc po przecinku