Gestern (13.3. oder 0313in US-amerikanischer Schreibweise) war Pi-Day. Aus diesem Anlaß gab es bei Heise auf der Webseite einen Artikel zur Berechnung der Zahl π. Es wurde auch gleich das Programm geschrieben in Python dazu mitgeliefert. Die Zahl π wurde darin nach dem Monte-Carlo-Verfahren bestimmt. Die Beschreibung fand ich toll und dachte mir, das kannst Du auch. Später kannst Du dann dein eigenes Programm mit dem von Heise vergleichen. Ich bin spazieren gegangen und habe mir dabei überlegt, wie ich das mache und habe es nach der Rückkehr gleich umgesetzt.
Das Programm war schnell geschrieben und kam zu dem erwünschten Ergebnis. Danach habe ich das Programm von Heise aber gar nicht ans Laufen bekommen, da bei mir auf dem Rechner anscheinend eine Library fehlt, die Library zur Darstellung der berechneten Punkte auf einem Grafen. Leider stand nirgendwo, wo ich diese her bekomme. Mein Programm ist jetzt also bei weitem nicht so elegant, wie das Programm von Heise, dafür funktioniert es.
import random
import math
def Abstand(a,b):
c=pow(a*a +b*b, 0.5)
return c
def Zufallszahl():
z=random.random()
return(z)
Zaehler=1
Treffer=0
while abs(Treffer/Zaehler*4-3.1416)>.0001:
x=Zufallszahl()
y=Zufallszahl()
d=Abstand(x, y)
if d<1:
Treffer=Treffer+1
Zaehler=Zaehler+1
print(Zaehler, Treffer, Treffer/Zaehler*4)