function x = fixpunkt(phi,x0)
%FIXPUNKT: Demo-Programm zur Fixpunkt-Iteration
% Das Verfahren iteriert gemäss der Vorschrift
%       x_i+1 = phi(x_i)
% bis Änderungen unter eine Toleranzschwelle sinken
% oder maximale Iterationszahl überschritten wird.
% 
% Eingabe   phi... Funktion, deren Fixpunkt gesucht ist
%                  (mit Funktions-Henkel @)
%           x0 ... Startwert oder -vektor
% Ausgabe   x .... Bei Konvergenz: Fixpunkt, 
%                  sonst: NaN
% Beispiel  fixpunkt(@cos,1) 
%
%------------------------------------------NMI,SS09-18 C.Brand
itmax = 100;                % maximale Iterationszahl
errlim = 1.e-9;             % Fehlerschranke
for i=1:itmax
    x = phi(x0);            % Funktionsauswertung  
    if norm(x-x0)<errlim    % Abbruchkriterium: 2-Norm des Fehlers     
        return
    end
    x0 = x;
end
x = NaN;
end
