fun = @(x) [x(1)^2 + x(2)^2 - 4; % Prima equazione x(1) - x(2) - 1]; % Seconda equazione x0 = [1; 1]; % Punto iniziale sol = fsolve(fun, x0)
