Mathematik für Physiker - Numerische Physik: Modellierung

5.5. NUMERISCHE INTEGRATION IN MATLAB 191
Abbildung 5.9: Bestimmung von π mit
Monte Carlo
§ 749 Das Verfahren lässt sich auch in MATLAB anwenden. Dazu benötigen wir Zufallszahlen,
die mit der Funktion rand erzeugen. Wird rand ohne Argument auf gerufen, so erzeugt
es im Intervall (0,1) gleich verteilte Zufallszahlen. 5 Werden Argumente an rand übergeben,
so kann rand z.B. Matrizen aus Zufallszahlen erzeugen
>> rand(3) ←↪
ans =
0.2028 0.2722 0.7468
0.1987 0.1988 0.4451
0.6038 0.0153 0.9318
oder einen Vektor
>> rand(1,3) ←↪
ans =
0.7095 0.4289 0.3046
rand
§ 750 Für die Bestimmung der Fläche des Kreises benötigen wir die folgende Sequenz:
n=1000; z=0;
x=2*rand(1,n)-1; y=2*rand(1,n)-1;
r= sqrt(x.*x + y.*y);
for i=1:n
if r(i)