Вычислить интеграл методом Монте-Карло

##
function f(x: real): real := 5 - exp(-2 * x);
function i(x: real): real := 5 * x + exp(-2 * x) / 2;
 
  randomize;
  var (xl, xc, xr, hn, hp, eps) := (-1.0, -ln(5) / 2, 1.0, f(-1), f(1), 0.00001);
  var (wn, wp, x, y, s) := (xc - xl, xr - xc, 0.0, 0.0, 0.0);
  var cn, cnf, cp, cpf: Int64;
  (cn, cnf, cp, cpf) := (1, 1, 1, 1);
  loop Int64(1 / sqr(eps)) do
    begin
      x := random(xl, xr);
      if x < xc then
        begin
          y := random(hn, 0.0);
          cn += 1; 
          if y >= f(x) then cnf += 1
        end
      else
        begin
          y := random(0.0, hp);
          cp += 1;
          if y <= f(x) then cpf += 1
        end
    end;
  s := cpf / cp * wp * hp + cnf / cn * wn * hn;
  $'Метод Монте-Карло: {s}'.PrintLn;
  $'  Точное значение: {i(xr) - i(xl)}'.PrintLn

##
function f(x: real): real := 5 - exp(-2 * x);
function i(x: real): real := 5 * x + exp(-2 * x) / 2;
 
  randomize;
  var (xl, xc, xr, hn, hp, eps) := (-1.0, -ln(5) / 2, 1.0, f(-1), f(1), 0.00001);
  var (wn, wp, x, y, s) := (xc - xl, xr - xc, 0.0, 0.0, 0.0);
  var cn, cnf, cp, cpf: Int64;
  (cn, cnf, cp, cpf) := (1, 1, 1, 1);
  var n := Int64(1 / sqr(eps));
  var m := Int64(n / hp * hn / wp * wn);
  n -= m;
  loop n do
    begin
      x := random(xc, xr);
      y := random(0.0, hp);
      cp += 1;
      if y <= f(x) then cpf += 1
    end;
  loop m do
    begin
      x := random(xl, xc);
      y := random(hn, 0.0);
      cn += 1; 
      if y >= f(x) then cnf += 1
    end;
  s := cpf / cp * wp * hp + cnf / cn * wn * hn;
  $'Метод Монте-Карло: {s}'.PrintLn;
  $'  Точное значение: {i(xr) - i(xl)}'.PrintLn

## uses NumLibABC;
var o := new StopWatch;
o.Start;
var I := new Quanc8(x -> 5 - exp(-2 * x), -1, 1, 1e-5, 1e-5);
var r := I.Value;
o.Stop;
r.Println;           // (6.37313959215278,1.85660046042368E-12,0,33)
// Значение 6.37313959215278
// абсолютная погрешеость 1.85660046042368E-12
// индикатор надежности 0 (надёжное решение)
// потребовалось вычислить функцию в 33 точках
Println(o.Elapsed);  // 00:00:00.0021802 - время вычисления примерно 2 мс