Przejdź do treści

Metoda Monte Carlo

Opis problemu

Obliczanie wartości liczby PI

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
function MonteCarloPi(pointsCount)
    pointsInCircleCount = 0
    centerX = 1
    centerY = 1
    radius = 1

    for i in 1:pointsCount
        x = rand() * 2
        y = rand() * 2
        distance = (x - centerX)^2 + (y - centerY)^2

        if distance <= radius^2
            pointsInCircleCount += 1
        end
    end

    return (4 * pointsInCircleCount) / pointsCount
end

println(MonteCarloPi(100000000))