Лаба горит!

2.2.3. Формула Симпсона
:{
Формула Симпсона для приближенного вычисления определенного интеграла имеет вид

S = h/3*(( y0 + y2*k ) + 4*( y1 + y3 + ... + y2*k-1 ) + 2*( y2 + y4 + ... + y2*k-2 )),

где n - четно ( 2*k ),
h = ( b-a )/n,
xj = a+ j*h ( j = 0, 1, 2, ..., n ),
yj = f( xj ).
------------------------------------------------
для функций
1.f(x) = (x2-5) *(10-x); [0, 10].
2.f(x) = exp(2*x) * sin(2*x2 + 1); [0, 1].
3.f(x) = arccos(exp(-(3*x)(1/3) ); [0.2, 0.56].
через Процедурные типы данных

пример работающей проги:
program Tables;

type
Func = function(X,Y: integer): integer;

function Add(X,Y: integer): integer; far;
begin
Add := X + Y;
end;

function Multiply(X,Y: integer): integer; far;
begin
Multiply := X*Y;
end;

function Funny(X,Y: integer): integer; far;
begin
Funny := (X+Y) * (X-Y);
end;

procedure PrintTable(W,H: integer; Operation: Func);
var
X,Y : integer;
begin
for Y := 1 to H do
begin
for X := 1 to W do Write(Operation(X,Y):5);
Writeln;
end;
Writeln;
end;

begin
PrintTable(10,10,Add);
PrintTable(10,10,Multiply);
PrintTable(10,10,Funny);
end.
:(