function Loese

    globals;

    ep0  = ep;
    ep   = ep0*2^rmax;

    for relax = 1:rmax
        ep  = ep/2;
        for it = 1:itmax
            switch Methode
                case 1         % Abstieg
                    KriegF;
                    du = -F;
                case 2         % Implizit
                    KriegD;
                    du = (I + mue*D/(h*h)) \ v - u;
                case 3         % Newton
                    KriegF;
                    KriegDF;
                    du = -DF \ F;
            end
            Liniensuche;
            du = alf*du;
            u  = u + du;
            if (norm(du) <= (tol*norm(u)))
                break;
            end
        end
        Ergebnis;
        disp(['relax = ',num2str(relax),': ep = ',num2str(ep),', Iterationen bei Austieg = ',num2str(it)]);
    end

end