Метод Симпсона

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

function f(x,c,d: real): real;

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.