Monte_carlo_pi.py

# Heitetään pisteitä satunnaisesti neliölle
# Osa pisteistä osuu neliön sisällä olevalle neljännesympyrälle
# Lasketaan tuloksista Piin arvo
# Juhani Kaukoranta 25.3.2019
from numba import njit
import random
import time
@njit
def monte_carlo_pi(n):
    accum = 0
    for i in range(n):
        x = random.random()
        y = random.random()
        if (x**2 + y**2) < 1.0 :
            accum += 1
    return 4.0*accum/n
pisteita = int(input("Kuinka monta pistettä arvotaan "))
time0 = time.perf_counter()
print("Pii = ",monte_carlo_pi(pisteita))
time1 = time.perf_counter()
print("Aikaa kului ",time1-time0," sekuntia")