function [x,nit] = fikspunkt(g,x0,tol)
% [x,nit]=fixpunkt(g,x0,tol) 
% loeser en ikke-lineaer ligning paa formen 
%                  x = g(x). 
% Inn:   g er et funkjsonshaandak til funksjonen g (@g).
%        x0 er en startverdi for iterasjonene
%        tol er en feiltoleranse.
% Ut:    x er loesningen (hvis den er funnet)
%        nit er antall iterasjoner. 

x = x0;
feil = 10;
MaxIter = 100;
nit = 0;             % Antall iterasjoner
while (feil > tol)
    xny = g(x);
    feil = abs(xny-x);
    x = xny;
    nit = nit + 1;
    if nit >= MaxIter
        warning('For mange iterasjoner');
        break;
    end
end