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

Program integral_simpson;
Uses crt;
var
a,b,s: real;
n:integer;
function f(x:real):real;
begin
f:=(2*sin(x)*cos(x)+sin(x));
end;
procedure simpson(a, b: real; n: integer; var s: real); 
var h, h2, x : real; {h - шаг интегрирования, h2 - полушаг интегрирования, x - текущее значение переменной интегрирования} 
i : integer; 
begin h := (b - a)/n; {определение шага интегрирования} 
h2 := h/2;
 x := a; 
 s := (f(a)+f(b))/2+2*f(a + h2); 
 for i := 1 to n-1 do 
 begin 
 x := x + h;
 s := s + f(x) + f(x + h2);
 end;
 s := s*(h/3);
 end;
 
begin
clrscr;
write ('Введите границу промежутка интегрирования a = ');
read (a);
write ('Введите границу промежутка интегрирования b = ');
read (b);
write ('Введите число разбиений n = ');
read (n);
simpson (a,b,n,s);
Writeln ('Значение интеграла = ', s:3:4);
end.