Программа

{ Вычисление интеграла
  методом Симпсона }
TYPE Func=function(x:real):real; { Тип функция }
VAR Left,Right:real;       { Левая и правая границы интервала интегрирования }
    Break:integer;         { Разбиение }
FUNCTION F(X:real):real;far; { Вычисление интегрируемой функции }
  Begin
    F:=x
  End;
FUNCTION Integral(Left,Right:real;Break:integer;F:Func):real;
                                              { Функция вычисления интеграла }
  Var Sum1,Sum2,                              { Промежуточные суммы }
      X,                                      { Промежуточная абсцисса }
      LengthStep,                             { Длина шага разбиения }
      LeftY,RightX,             { Левая и правая границы шага интегрирования }
      RightY:real;              { Значение интегрируемой функции
                                  на правой границе шага интегрирования }
      i:integer;                { Счетчик }
  Begin
    LengthStep:=(Right-Left)/Break/2; { Вычисление длины шага интегрирования }
    X:=Left+LengthStep; { Начальное определение промежуточной абсциссы }
    Sum1:=0;     { Обнуление промежуточной суммы 1 }
    Sum2:=F(X);  { Начальное определение промежуточной суммы 2 }
    for i:=2 to Break do { Цикл набора промежуточных сумм }
      begin
        X:=X+LengthStep;   { Переопределение абсциссы }
        Sum1:=Sum1+F(X);   { Набор промежуточной суммы 1 }
        X:=X+LengthStep;   { Переопределение абсциссы }
        Sum2:=Sum2+F(X);   { Набор промежуточной суммы 2 }
      end;
    Integral:=LengthStep/3*(F(Left)+F(Right)+2*Sum1+4*Sum2)
                                        { Определение значения интеграла }
  End;
BEGIN
{ Ввод данных: }
  Write('Введите левую границу интервала интегрирования ');
  ReadLn(Left);
  Write('Введите правую границу интервала интегрирования ');
  ReadLn(Right);
  Write('Введите разбиение интервала интегрирования ');
  ReadLn(Break);

  WriteLn('=',Integral(Left,Right,Break,F)) { Вычисление и вывод результата }
END.