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

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

Определенный интеграл I <= ò f(x)dx по методу Монте-Карло

по формуле I = (1/n)* å (f(xi))/(g(xi)) ,где

распределения Увспомогательной случайной величины X, причем ò g(x)dx = 1 ,

В программе

Программа написана на языке TURBO PASCAL 7.0

Program pmk;

Uses crt;

Var k,p,s,g,x,Integral : real;

BEGIN

randomize;

readln(a);

readln(b);

readln(n);

for i:= 1 to n do begin {проведем

промежутка [0;1]}

delay(1); {задержка,чтобы произвольные значения не повторялись<}

end;{конец испытаний<}

Integral:=(1/n)*k*s ;

readln;

END.

Требуется ввести промежуток интегрирования и количество испытаний, интегрируемая функция уже задана в программе(но ее можно поменять).

3 3

ò(x+1)dx = 6а ;а ò (x*x)dx = 9; (По методу Ньютона-Лейбница).

Blog