Задача про модули

[Tattoquardas]

Написать модуль, который содержит подпрограммы приблизительного подсчета интеграла (x^2 - 1)3/2dx (от 0 до 1) методом прямоугольников и Симпсона. подпрограммы должны содержать как один из параметров функцию типа function ff(x:real):real.

[IT-man]

симпсон

{код Delphi}
type
  TDoubleFunc = function(X: Double): Double;

function Integral(A, B, Eps: Double; F: TDoubleFunc): Double;

  function InternalCalc(A, B: Double; F: TDoubleFunc; N: Integer): Double;
  var
    x, dx: Double;
    i: Integer;
  begin
    dx := (B - A) / N;
    Result := 0;
    x := A;
    for i := 1 to N do
    begin
      Result := Result + dx * (F(x) + 4 * F(x + dx / 2) + F(x + dx)) / 6;
      x := x + dx;
    end;
  end;

var
  N: Integer;
  Prev: Double;
begin
  Result := InternalCalc(A, B, F, 4);
  N := 4;
  repeat
    Prev := Result;
    N := N shl 1;
    Result := InternalCalc(A, B, F, N);
  until (Result = 0) or (Abs((Result - Prev) / Result) < Eps);
end;

(с) Dimka Maslov

методом прямоугольников

каких ? левых? правых?

{код Delphi}
function LeftPramInt(a, b: Double; eps: Double; IntF: TDoubleFunc): Double;
var
  //S - площадь на предыдущей итерации,
  //step - "толщина" прямоугольника
  //gran - передвигаемая от a до b граница
  //n - число прямоугольников, удваивается на каждой итерации
  S, step, gran: Double;
  n: integer;
begin
  //Сначала приближение одного прямоугольника
  step := b - a;
  Result := IntF(a) * step;
  n := 1;
  repeat
    S := Result;
    n := n * 2;
    step := (b - a) / n;
    Gran := a;
    Result := 0;
    //Ниже - просто вычисляем площади новых прямоугольников
    while gran < b do
    begin
      Result := Result + IntF(gran) * step;
      gran := gran + step;
    end;
  until abs(S - Result) <= eps;
end;

function RightPramInt(a, b: Double; eps: Double; IntF: TDoubleFunc): Double;
var
  //S - площадь на предыдущей итерации,
  //step - "толщина" прямоугольника
  //gran - передвигаемая от a до b граница
  //n - число прямоугольников, удваивается на каждой итерации
  S, step, gran: Double;
  n: integer;
begin
  //Сначала приближение одного прямоугольника
  step := b - a;
  Result := IntF(b) * step;
  n := 1;
  repeat
    S := Result;
    n := n * 2;
    step := (b - a) / n;
    Gran := b;
    Result := 0;
    //Ниже - просто вычисляем площади новых прямоугольников
    while a < gran do
    begin
      Result := Result + IntF(gran) * step;
      gran := gran - step;
    end;
  until abs(S - Result) <= eps;
end;

(c) Алексей Глеб

[Tattoquardas]

одним из методов прямоугольников.
программа пишет, имя типа Double неизвестное... что делать?

[IT-man]

имя типа Double неизвестное... что делать

Заменить везде Double на Real