Составить подпрограмму для вычисления определенного интеграла методом трапеций

uses crt;
var a,b,S,h,integ: real;
    i,n: integer;
 
function f(x: real): real;
begin
 f:=2*x+1; //сюда интеграл.
end;
 
begin
 clrscr;
 write('a/b/n: '); //нижняя/верхняя границы и число разбиений (обычно 10 хватает)
 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.