Obliczanie całki metodą Monte Carlo:
a := 0.2
b := 1.2
n:=50
n
przykładowa funkcja w przedziale (a ,b): z := a, a + A.. b y(z) := sin(2 z) Max_y := 1
Losowanie punktów w przedziale całkowania: i := 0.. 2999 p. := a + md(b - a) q. := md(Max_y)
Ocena czy "dobry strzał" to znaczy czy punkt znajduje się pod krzywą (1) czy nie (0): d°biy_stizal. : = if^q. < y ^, 1, O^J wszystlde : = length( dobiy_stizal) wszystkie = 3 • 10^
trafione := ^dobry.staal. ttafione = 2.506-103 i
Przybliżona wartość całki - wg met. Monte Carlo:
MC h&one .Maxjr.^b _ ^ MC = 0.335
wszystlde
Całka wyliczona konwencjonalnie wynosi zaś:
rb
całka :=
y(x)dx
całka = 0.329
• a