Вычислить интеграл используя квадратурную формулу Симпсона.

program funkcia;
var
   h,a,b,s,s1,x,int:real;
   i,m:integer;
   
 
begin
   f :=(x+1)*(1-x*x*x) (1,9);
end;
 
begin
    a:=1;
    b:=9;
    read(m)

const
  a = 1;{нижний предел}
  b = 9;{верхний предел}
  e = 0.0001;{точность}
 
function F(x: real): real;
begin
  F := (x + 1) * power((1 - x * x * x), 1 / 5);
end;
 
var
  n, i: integer;
  h, k, s1, s2: real;
 
begin
  write('Количество разбиений n=');
  readln(n);//начальное число
  s1 := 0;//итоговый интеграл
  h := (b - a) / n;//шаг
  s2 := h * F(a);//текущий интеграл
  repeat
    s1 := s2;//запомним текущий
    h := (b - a) / n;
    k := a;//левый край
    for i := 1 to n - 1 do
    begin
      k := k + h;//наращиваем
      if i mod 2 = 0 then 
        s2 := s2 + 2 * f(k)//считаем новый интеграл
      else 
        s2 := s2 + 4 * f(k);
    end;
    s2 := (s2 + f(a) + f(b)) * h / 3;//конечное значение
    n := n * 2;//если точность не достигнута удваиваем количество разбиений
  until abs(s1 - s2) <= e;
  write('Интеграл=', s2:0:6);
end.