[lmikle]

Общий алгоритм:
На интервале a-b случайным образом генерируем много точек (x,y). Далее проверяем, где они находятся относительно функции (выше или ниже). Соотв. складываем в 2 счетчика. В конце считаем площадь прямоугольника и умножаем на кол-во точек под функцией и делим на общее кол-во точек. Код примерно такой:

Код:

function F(x : Double) : Double;
begin
  // твоя функция
end;

function GetIntegralMonteCarlo(nProb, a, b : Integer) : Double;
var
  I : Integer;
  nIn : Integer;
  yMax  : Double;
  x, y : Integer;
begin
  yMax := Max(F(a),F(b));
  nIn := 0;
  nOut := 0;
  For I := 1 To nProb Do
    Begin
      x := Random(b-a) + a;
      y := Random(yMax);
      If y <= F(x)  Then Inc(nIn);
    End;
  Result := (b-a) * yMax * nIn / nProb;
end;