Решение системы нелинейных уравнений с 3 переменными (или как понять что решения не существует) - Matlab - Ответ 14060123

function main
clear, clc
clear functions             % clear persistent variables
% x(1)->a. x(2)->b. x(3)->c
 
    function y = F(x)
        persistent iter x00
        if isempty(iter), iter = 1; end
        if isempty(x00), x00 = x; end
        y(1) = x(3)*x(1)*5.451 + x(3)*x(2)*5.662 + x(3)*(1 -x(1)- x(2))*6.135 +x(1)* (1 - x(3)) * 5.449 + x(2) * (1 - x(3)) * 5.641 + (1 - x(3)) * (1 - x(1) -x(2)) * 6.094 - 6.094;
        y(2) = x(3) * x(1) * 4.7 + x(3) * x(2) * 5.2 + x(3) * (1 - x(1) - x(2)) * 4.88 + x(1) * (1 - x(3)) * 5.81 + x(2) * (1 - x(3)) * 6.86 + (1 - x(3)) * (1 - x(1) - x(2)) * 6.7 - 6.7;
        y(3) = x(3) * x(1) * 3.6 + x(3) * x(2) * 2.9 + x(3) * (1 - x(1) - x(2)) * 2.2 + x(1) * (1 - x(3)) * 2.78 + x(2) * (1 - x(3)) * 1.35 + (1 - x(3)) * (1 - x(1) - x(2)) * 0.72 - x(1) * x(2) * (x(3) * 0.22 + (1 - x(3)) * 0.22) - x(1) * (1 - x(1)- x(2)) * (x(3) * 1.2 + (1 - x(3)) * 1.2) -x(2) * (1 - x(1) - x(2)) * (x(3)* 0.72 + (1 - x(3)) * 1.21) - x(3) * (1 - x(3)) * (x(1) * 0.04 +x(2) * 0.026 + (1 -x(1) - x(2)) * 0.17);
        
        plot([iter iter+1],[x; x00]','.-')
        iter = iter + 1;
        x00 = x;
    end
 
   x0 = [1 0 0]; % начальное предположение
   % Option to display output
options = optimoptions('fsolve','Display','iter');
figure(1)
hold on
X = fsolve(@F,x0,options);
hold off
a = X(1);
b = X(2);
c = X(3);
end