X1=[0 0.2 0.4 0.6 0.8 1.0 1.01];
Y1=[0 1.5 4   7   9.3 10  50];

X2=[0 0.2 0.4  0.6 0.8 1.0 1.2];
Y2=[0 0.1 0.75 2   3.9 6   9];

X3=[0 0.2  0.4  0.6 0.8 1.0 1.2];
Y3=[0 0.05 0.2  0.5 1   1.5 2.2];

X=zeros(12,1);
Y=zeros(12,1);
for i = 1:1:12
[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(i);
X(i)=a1+a4;
Y(i)=a1*a4+a2;
A3=b1*(a2+a1*a4)+b2*(b4+a1*b7-b6)-a6*(a4*b6+a2*b7);
rezultat_proverki=b1*(a1*a4+a2)/(A3-b2*b4);

[ MUn, sigma, Tn ] = ABUparams( i );
fprintf('режим %d\n', i);
fprintf('sigma = %f\n\n', sigma);

%MUn = 0.8;

a2_2 = a2+ (sigma*a3*a4)/(57.3*b4-sigma*a7);
a4_2 = a4*57.3*b4/(57.3*b4-sigma*a7);
a1_2 = a1 + (57.3*a3*b4*MUn)/(57.3*b4-a7*sigma);

if (i==4 || i==9 || i==12)
    a4=a4_2;
    a2=a2_2;
    %a1=a1_2;
    %MUn = 2*0.4*(sqrt(a1*a4_2+a2_2)-a1-a4_2)/a3;
end


if (i==1 || i==4 || i==6 || i==7 || i==8 || i==9 || i==10 || i==11 || i==12)
    %a4=a4_2;
    %a2=a2_2;
    a1=a1_2;
end

%a4=a4_2;
%a2=a2_2;
%a1=a1_2;

X(i)=a1+a4;
Y(i)=a1*a4+a2;
end

p=plot(X1,Y1, X2, Y2, X3,Y3,X,Y,'x','LineWidth',2)
%p=plot(X1,Y1, X2, Y2, X3,Y3,X,Y,'x','LineWidth',2, {0,1,0,0,10})
%p=plot(X1,Y1, X2, Y2, X3,Y3,X,Y,'x','LineWidth',2,...
%                'MarkerEdgeColor','k',...
%                'MarkerFaceColor','g',...
%                'MarkerSize',10);
grid on;
set(p,'Color','red','LineWidth',2);
%set(p,'Color','red','LineWidth',2);
%axis([0.5 5 10]);
%ylim([0, 11]);

%p=plot(X1,Y1, X2, Y2, X3,Y3,X,Y,'x','LineWidth',2)