Объясните, плз, строчку кодо

[Дмитрий Огнев]

программа вычисляет интеграл по методу трапеции
Вопрос состоит в следующем: Формула в моем коде (отмеченная " * ") считает конечный результат как-то иначе, по сравнению со строчкой кода из интернета ( ** ), объясните что не так

мой код:

program qq;
var a,b,n,i,x:real;
    S,h: double;
begin
a:=0;
b:=1;
S:=0;
write('Введите количество делений отрезка: ');
readln(n);
h:=(b-a)/n;
i:=0;
 while i<=n do
  begin
    x:=a + h*i;
    S:=S+( 1/(1+x*x*x) );
    i:=i+1;
  end;
S:=h*( ( 1/(1+a*a*a) ) + ( 1/(1+b*b*b) ) )/2+S; *
writeln('S= ',S);
readln;
end.

код из интернета:
uses crt;
var a,b,S,h,integ: real;
    i,n: integer;
function f(x: real): real;
begin
 f:=1/(1+x*x*x);
end;
begin
 clrscr;
 write('a/b/n: ');
 readln(a,b,n);
 h:=(b-a)/n;
 for i:=1 to n-1 do
  S:=s+f(a + h*i);
 integ:=h*((f(a)+f(b))/2+s); **
 writeln('S=',integ);
 readln;
end.

[p51x]

Есть ли разница между h*((a+b) + c) и h*(a+b) + c?

[Дмитрий Огнев]

Согласен, тут ошибка. Но я это иправил и все равно что-то не то
h*(((1/(1+a*a*a))+(1/(1+b*b*b)))/2+s) вот так правильно, но ответ все равно не получается, тк если n=5, то результат 1,13..., а нужно 0,8...

[p51x]

Давайте дальше.
Тут идет суммирование по внутренним точкам

 for i:=1 to n-1 do
  S:=s+f(a + h*i);

а тут включая граничные

i:=0;
 while i<=n do
  begin
    x:=a + h*i;
    S:=S+( 1/(1+x*x*x) );

[Дмитрий Огнев]

Спасибо огромное, очень помогли