function [ IE, TF1 ] = getKursRieTFvectors( MUe, Ie )
% Возвращает 2 вектора расчетных значений Ie-пси, Тф1

%MU=zeros(12,1);
IE=zeros(12,1);
TF1=zeros(12,1);

%t=2;
%[ MUe, Ie ] = getKursRVectors( t );
%asd
for k = 1:1:12
   fprintf('режим: %d\n', k);
   
   [ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(k);
   [ MUn, sigma, Tn ] = ABUparams( k );
   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);
   
   i=k;
   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
    %a2=a2_2;
    %a4=a4_2;

   %sigma = 57.3*b4*(4-(a1*a4+a2))/(4*a7+a3*a4-a2*a7);
%sigma = 57.3*b4*(4 - a1*a2 -a2 )/(4*a7 + a3*a4 - a2*a7);

    %MUn = 2*0.4*(sqrt(a1*a4+a2)-a1-a4)/a3;
    %Tn = a3/(a3*a4-a2*a7);
   
   
   
   A0 = 1;
   A1 = a1 + a4 + b1 + MUe(k)*b3; %
   A2 = b1*(a1+a4) + a2 + a1*a4 + b2*b7 + MUe(k)*b3*(a1+a4) + Ie(k)*b3; %
   A3 = b1 * (a1*a4 + a2) +b2*(a1*b1 + b4) + MUe(k)*b3*(a2+a1*a4) + Ie(k)*b3*(a1+a4); %
   %A4 = a1*b2*b4 + Ie(k)*b3*(a2+a2*a4);
   A4 = a1*b2*b4 + Ie(k)*b3*(a2+a1*a4);
   
   A0_2 = Tn*A0;
   A1_2 = A0 + Tn*A1 + MUn*Tn*a3;
   A2_2 = A1 + Tn*A2 + MUn*Tn*(a3*(a4+b1+MUe(k)*b3 ) - a2*a7 );
   A3_2 = A2 + Tn*A3 + MUn*Tn*(a3*(a4*b2+b2*b7) + Ie(k)*a3*b3 -MUe(k)*b3*(a2*a7+a3*a4) - a2*a7*b1 );
   A4_2 = A3 + Tn*A4 + MUn*Tn*(a3*b2*b4 - Ie(k)*b3*(a2*a7+a3*a4)); %
   A5_2 = A4;
   
   syms x;
   S = solve(A0_2*x^5 + A1_2*x^4 + A2_2*x^3 + A3_2*x^2 + A4_2*x + A5_2);
   fprintf('roots:\n');
   
   %S = solve(A0*x^4 + A1*x^3 + A2*x^2 + A3*x + A4);
   %1/Tn
   %S
   %min(S.x)
   
   one_div_lambda = 9999999;
   lambda_i = 9999999; % минимальный по абс. величине вещественный корень
   % минимальные по модулю комплексно-сопряженные корни уравнения
   lambda_j = 999999;
   lambda_j_2 = 9999999;
   noRealRoots = true;
   
   for idx = 1:numel(S)
       element = S(idx);
       element = element(1);
       if isreal(element)% && (abs(element)>abs(lambda_i))
           %if 
           noRealRoots = false;
           fprintf('корень: %f\n', double(element));
           
           % минимальный лямбда_i
           if (double(abs(lambda_i))>double(abs(element)))
              lambda_i=element;
              %element
              
           end
           
           if (one_div_lambda>abs(1/double(element)))
               one_div_lambda=abs(1/double(element));
           end
           
           % 
       else % complex
           fprintf('корень: %s\n', char(element));
            if (double(abs(element))<double(abs(lambda_j)))
              lambda_j=element;
           end
       end

       % solve(A0*x^4 + A1*x^3 + A2*x^2 + A3*x + A4)
        
   end
    
   %lambda_i
   %lambda_j
   %double(abs(lambda_i))
   lambda_j_2 = (lambda_j)';
   
   %lambda_i
   one_div_lambda = abs(1/lambda_i);
   T1 = double((one_div_lambda)/(sqrt(abs(lambda_j)^2)));
   %T1 = 0.4;
   %A5_2
   if (double(abs(lambda_i))<double(abs(sqrt(lambda_j*lambda_j_2))))
       IE(k) = 0.25*A5_2  /(Ie(k)*b3*b4*a2*T1);
       
   else
       IE(k) = 0.125*A5_2 /(Ie(k)*b3*b4*a2*T1);
       %IE(k) = 0;
   end
   
   fprintf('Ie=%f\n', Ie(k));
   fprintf('T1=%f\n', T1);
   
   if (noRealRoots)
       %S
       %lambda_j
   else
       fprintf('1/Tn = %f', 1/Tn);
   end
    
   fprintf('\n--------\n');
    
    
   %[mu_korr, i_korr, nu_korr] = get_korr(k, 2);
   %IE(k) = (9.48-b1*t)/(b3*t);
   TF1(k) = b7/b4;
   %NU(k) = nu_korr;
end
end