function [x1, n] = newtonsys (f, J, x0, eps, nmax)
% Az eljárás a newton módszer segítségével határozza meg egy többváltozós 
% nemlineáris egyenletrendszer megoldását
%
% Bemenő paraméterek:
%   f - többváltozós egyenletrendszer
%   J - az egyenletrendszer Jacobi mátrixa
%   x0 - kezdő érték
%   eps - leállási feltétel, közelítés küszöbértéke
%   nmax - leállási feltétel, maximális iterácó szám
% Kimenet: x1 - a megoldás helye, n - iterációk száma
    dx = J(x0)\-f(x0); % első iteráció
    x1 = x0 + dx;      % első iteráció
    n = 1;
    while norm(x1-x0)>eps && n<=nmax
        x0 = x1;
        dx = J(x0)\-f(x0);
        x1 = x0 + dx;
        n = n + 1;
    end;
end