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

uses crt;
var a,b,S,h,integ: real;
nbsp;nbsp;nbsp;nbsp;i,n: integer;
nbsp;
function f(x: real): real;
begin
nbsp;f:=2*x+1; //сюда интеграл.
end;
nbsp;
begin
nbsp;clrscr;
nbsp;write(a/b/n: ); //нижняя/верхняя границы и число разбиений (обычно 10 хватает)
nbsp;readln(a,b,n);
nbsp;h:=(b-a)/n;
nbsp;for i:=1 to n-1 do
nbsp;nbsp;S:=s+f(a + h*i);
nbsp;integ:=h*((f(a)+f(b))/2+S);
nbsp;writeln(S=,integ);
nbsp;readln;
end.