In de onderstaande Monte-Carlosimulatie wordt de inhoud van een bol met straal $r$ benanderd.
Het experiment dat wordt herhaald is het kiezen van een willekeurig punt in de ruimte, binnen een 2r X 2r X 2r kubus waarvan het midden (zwaartepunt) op (0, 0, 0) ligt. Dan test je of het gekozen punt in de bol valt, of niet.
Vul zelf de uiteindelijke schatting van de inhoud van de bol in.
Let op: numpy.random.rand(n) geeft lijstje met $n$ randoms in $\big[0, 1\big\rangle$ terug.
In [ ]:
import numpy as np
r = 1
n = 25000
aantal_succes = 0
for i in range(n):
x, y, z = np.random.rand(3) * 2*r - r
if sqrt(x**2 + y**2 + z**2) < r:
aantal_succes += 1
print("Schatting inhoud bol:", VUL HIER JE BEREKENING VAN DE SCHATTING IN)
In [ ]: