% Octave Programm zur Bestimmung des Maschinen-Epsilons:

x=2;
while 1-x/2 < 1
  x=x/2;
end
x   % Unser Ergebnis
eps % Eingebauter Befehl fuer das Maschinenepsilon

% Bemerkung:
%
% x und eps koenne sich u.U. um einen Faktor 2 unterscheiden.
%
% Unsere Berechnung erfolgt nach der "akademischen" Definition des
% Maschinen-epsilons aus der Vorlesung.  Man weiss, dass die IEEE
% Floting-Point Zahlen 2^(-52) as Maschinenepsilon haben, weil die
% Mantisse aus 52 Bits (ohne das Vorzeichen) besteht.  Dieses Epsilon
% wird auch von Octave's "eps" Befehl so ausgegeben.
%
% Warum der Unterschied?  Die Hardware (meine jedenfalls) schneidet
% das Ergebnis einer Addition nicht ab, sondern rundet.  Dadurch kann
% das Ergebnis noch von Stellen abhaengen, die scheinbar nicht mehr
% signifikant sind.
