%
% HBaa, 16.09.2004
%
% Algorithmus zur Bestimmung von Pi
%

% Grafik
hold on
axis square

% Initialisierung
loops = 10000;

radius = 1.0;

kreis = 0;

% Ablauf
for i = 1:loops

%  zwei Zufallszahlen im Intervall ]0 ... radius[
    x1 = rand * radius;
    x2 = rand * radius;
    
%  Abfrage "Treffer innerhalb Kreisbogen" ?
    if (radius-x1)^2 + x2^2 <= radius^2
        kreis = kreis + 1;
        plot(x1,x2)
    end
end % i

% Ergebnis- / Ausgabe
pi_est = 4.0 * kreis / loops

% Fehler absolut
err_abs = pi - pi_est

% Fehler relativ
err_rel = 1.0 - pi_est / pi

%
% *** ENDE ***