Метод Монте-Карло

[Виктория Кусяк]

помогите пожалуйста, у меня задача методом статистических испытаний (метод Монте-Карло) вычислять определенные интегралы и находить погрешность вычислений, в Pascal оганизировать программу в Delphy))
Вот код на Паскале, только без погрешности((
Вычислить определенный интеграл.

program MonteCarlo;
uses
 crt;
const k=100;
Var a,b,c,d,ng,vg,x,y,s,integral : real;
 n,i,j : integer;
 integr:array[1..k]of real;
Function f(x,y:real):real;
 Begin
 f:=Sqrt(x+y);
 end;
Function nm(x:real):real;
 Begin
 nm:=3*x;
 end;
Function vm(x:real):real;
 Begin
 vm:=8*x;
 end;
BEGIN
clrscr; 
writeln('Vvedit znachennya granyts integruvannya ');
 write('a='); readln(a);
 write('b='); readln(b);
 writeln('Vvedit chyslo vyprobuvan:');
 readln(n);
 c:=nm(a);
 d:=vm(b);
 randomize;
 for j:=1 to k do
 begin
 s:=0; integral:=0;
 For i:=1 to n do
 begin
 x:=a+(b-a)*random;
 y:=c+(d-c)*random;
 ng:=nm(x);
 vg:=vm(x);
 If (y <= vg) and(y >= ng) then s:=s + f(x, y);
 end;
 integr[j]:=(b-a)*(d-c)*s/n;
 writeln(integr[j]:10:4);
 end;
 for j:=1 to k do
 Integral:=integral+ integr[j];
 writeln('Userednenyy integral=',(integral/k):10:4);
 readln;
END.

________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE]
(это кнопочка на панели форматирования с решёточкой #)
Не забывайте об этом!

Модератор.