Решить интеграл методом трапеций.

[Gy9vin]

Добрый день. Не получается решить интеграл методом трапеций.

Вот код программы.

uses crt;
var a,b,S,h,integ: real;
    i,n: integer;
 
function f(x: real): real;
begin
 f:=1/x;
end;
 
begin
 clrscr;
 write('a/b: ');
 readln(a,b);
 n:=10;
 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.

[Sciv]

Перед циклом добавьте:

S:=0;

а то у Вас начальное значение суммы неинициализировано никак