Approximative Bestimmung von π

© Prof. Dr. Jürgen Schnack, Herbstakademie 2007

Random uses the Wolfram rule 30 cellular automaton generator for integers. 
It uses a Marsaglia-Zaman subtract-with-borrow generator for real numbers.

Wie oft wollen wir würfeln?

Ng=10000;

Funktion

ViertelKreis[x_]=(Sqrt[1-(x^2)]);

Jetzt wird Ng mal gewürfelt

(Schauen Sie doch mal in der Hilfe nach, was Random macht.)

PiApprox=0;

For[n=1,n<Ng+1,n++,
x=Random[];
y=Random[];
If[(x^2)+(y^2)<1,PiApprox=PiApprox+1]
];

PiApprox=PiApprox*(4/Ng);

Ergebnis

Print["Zahl der Würfe = ",Ng];
Print["Approximatives Ergebnis für Pi = ",N[PiApprox]];