Как оформить функцию Pascal

[Тетрадь]

const
  n=10000;
var
  k1,k2,k3,k4: real;
  x,y: array [0..n] of real;
  h: real;
  i: integer;

function f(x,y: real): real;
begin
  f:=(3-x-sqr(y))/(2-x*sqr(y));
end;
begin
  h:=1;
  y[0]:=1;


  for i:=0 to n do
    x[i] :=i*h;
    for i:=0 to n-1 do begin
      k1:=f(x[i],y[i]);
      k2:=f(x[i]+h/2,y[i]+h*k1/2);
      k3:=f(x[i]+h/2,y[i]+h*k2/2);
      k4:=f(x[i]+h,y[i]+h*k3);
      y[i+1]:=y[i]+h/6*(k1+2*k2+2*k3+k4);
    end;
    for i := 0 to n do
      writeln('x[',i:2,'] = ',x[i]:0:1,'   y[',i,'] = ',y[i]:0:4);
end.

Нужно оформить эту часть программы в виде функции

for i:=0 to n do
    x[i] :=i*h;
    for i:=0 to n-1 do begin
      k1:=f(x[i],y[i]);
      k2:=f(x[i]+h/2,y[i]+h*k1/2);
      k3:=f(x[i]+h/2,y[i]+h*k2/2);
      k4:=f(x[i]+h,y[i]+h*k3);
      y[i+1]:=y[i]+h/6*(k1+2*k2+2*k3+k4);
    end;

Но так, чтобы она была независима от фун f

[Аватар]

type TMyFunction = function(x,y: real): real;

//процедура
procedure Abc(ff: TMyFunction);
var a: Real;
begin
  a:=ff(1.5,0.0125);
...
end;

//обращение
Abc(f);