Ошибка: Floating point division by zero

program Project1;
 
function fA(x: real; j,k:integer) : real;
begin
  result := (1+x*x/8)*(9/8-x)*(j*ln(x-1/8))*(9/8-x)*(k*ln(x-1/8))+
            ((j*ln(x-1/8))*(-1)+(9/8-x)*j*((j-1)*ln(x-1/8)))*
            ((k*ln(x-1/8))*(-1)+(9/8-x)*k*(k-1*ln(x-1/8)));
end;
function fB(x,a,b: real; k:integer) : real;
begin
  result := (9/64-(-14/8))/(9/8-1/8)*((k*ln(x-1/8))*(-1)+(9/8-x)*k*(k-1*ln(x-1/8)))+
  (9/8-x)*(k*ln(x-1/8))*((1+(x*x)/8)*(9/64+(-3.5-9/64)*(x-1/8)/(9/8-1/8))-(1+x)/x*x*x+ln(7));
end;
 
var
    h,a,b:real;
    prov:real;
    i,k,j:integer;
    x,y:array[0..3] of real;
    mA:array[1..3,1..3] of real;
    mB:array[1..3] of real;
Begin
    a:=1/8; b:=9/8;
    h:=(b-a)/3;
    begin
      writeln('Расспишем точки интегрирования');
      x[0]:=a;
      for I := 1 to 2 do
        x[i]:= x[0]+i*h;
      x[3]:=b;
      for I := 0 to 3 do
        writeln('x[',i,']: ',x[i]:6:4);
      readln;
    end;
//============================
    begin
      writeln('Найдем A и B');
      writeln('Матрица A:');
      for k := 1 to 3 do
        begin
        for j := 1 to 3 do
          begin
          for I := 0 to 3 do
            y[i]:= fA(x[i],j,k);
            mA[k,j]:=h/2*(y[0]+2*(y[1]+y[2])+y[3]);
            write(mA[k,j]);
          end;
        writeln;
        end;
    end;
 
end.

if x = 0 then Exit;