function x = fixpunkt(phi,x0)
itmax = 100;errlim = 1.e-9;for i=1:itmaxx = phi(x0);if norm(x-x0)<errlimreturnendx0 = x;endx = NaN;end