Метод Симпсона - Turbo Pascal - Ответ 436620

Program  Kurs;
uses crt;
{Функция, площадь которой нужно вычислить}
function f(x,c,d: real; n:integer): real;
begin
case n of
1: f:= x*x+3;
2: f:=c*d+x;
3: f:=d*cos(x)+c;
end;
end;
var  a,b,c,d,h,x,integ :real;
     n,k,I :integer;
begin
clrscr;
writeln ('vviberite funkciu');
writeln (' 1: f := x*x+3');
writeln (' 2: f:=c*d+x');
writeln (' 3: d*cos(x)+c');
readln(k);
write('Введите начало интегрирования a=');
readln(a);
write('Введите конец интегрирования b=');
readln(b);
writeln ('vvedite c,d=');
readln(c,d);
repeat
write('Введите количество разбиений интервала (четное число) n=');
readln(n);
if (n mod 2)<>0 then writeln('Число n было введено нечетное, оно было заменено на n=',n);
until n mod 2=0;
h:=(b-a)/n;
integ:=F(a,c,d,k)+F(b,c,d,k)+4*F(a+h,c,d,k);
for i:=1 to (n div 2)-1 do
 begin
  x:=a+2*h*i;
  integ:=integ+2*F(x,c,d,k)+4*F(x+h,c,d,k);
 end;
integ:=h*integ/3;
writeln;
writeln('Результат по методу Симпсона равен ',integ:0:2);
readln
end.