Untitled

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)