Метод средних прямоугольников

[Олеся Хусаинова]

Всем привет!Не мог понять в чем дело,помогите,пожалуйста!
Вычислить определенный интеграл методом средних прямоугольников на отрезке [0.15-0.75] от функции 1/(x^2(x^2+z^2)),z - положительный параметр

вот,что у меня получилось
var A, B, h, x, S,z : real;
N, I : integer;
function F(X: real): real;

begin
f := 1/(x*x*(x*x+z*z)) end;
begin
write ('Input A, B, N ');
readln (A, B, N);
begin
S:=0;
h:=(B-A)/N;
x:=h/2;
for i:=1 to N-1 do begin
S:=S+f(x);
x:=x+h
end;

Writeln ('Sum = ',S:10:7);
readln
end.

[BDA]

Проверьте правильность счета:

var
    A, B, h, x, S, z : Extended;
    N, i : integer;

function F(X, z : Extended) : Extended;
var
    d : Extended;
begin
    d := x * x * x * x + x * x * z * z;
    f := 1 / d;
end;

BEGIN
    write('Input A, B, N : ');
    readln(A, B, N);
    S := 0;
    z := 1;
    h := (B - A) / N;
    x := A + h / 2;
    for i := 1 to N do
    begin
        S := S + f(x, z) * h;
        x := x + h;
    end;
    Writeln('Sum = ', S:10:7);
    readln;
END.