function [x k] = NewtonSys(x0,LGSSolv,Func,TOL,ItMax,flag)

x = x0;
[f Df] = Func(x);

k=0;
while norm(f)>=TOL & k<ItMax
    
    z = LGSSolv(Df,-f');
    x = x+z;
    
    [f Df] = Func(x);
    k = k+1;
    if flag 
        fprintf('Res.: %.2e nach %5d Iterationen\n',norm(f),k)
    end
end

end

