% 12 режим
%0,566	1,087	0,527
mu=0.566;
i=1.087;
nu=0.527;
b1 = 0.62;
b3 = 4.2;



ie_psi=6;

TP_r = [0.9026 0.3722  0.0667 1.6469 0.2304 1.3227 0.6239 0.6037 2.8329 1.5541 1.5319 4.6747];
A1_r = [6.98   19.4    14.85  3.93   10.09  7.15 6.99 6.7 3.32 3.92 4.80 2.61];
A2_r = [29.93  117.21  166.82 15.95  21.07  658.58 65.32 191.3 222.08 14.8 42.71 66.51 7.09];
A3_r = [64.05  730.63  1453.3 21.07  658.58 65.32 191.3 222.08 14.8 42.71 66.51 7.09];
A4_r = [84.12  1128.36 1601.25 28.01 791.44 104.51 261.46 294.71 24.06 68.48 101.83 13.8];

A12_r = [26.83  75.77   64.64   15.50  44.87 27.64 30.59 28.72 13.40 17.22 19.18 11.26];
A22_r = [142.22 786.86  985.51  66     505.58 138.8 225.63 223.97 52.99 116.17 103.35 38.06];
A32_r = [315.56 3338.77 7779.21 110.06 3408.25 322.33 903.33 1023.49 80.91 439.76 294.06 46.26];
A42_r = [297.85 3638.43 5700.60 103.1  2868.46 365.14 940.3 1071.16 86.36 272.05 366.03 48.27];
Tn_r = [5.4787 3.0259 1.6111 9.3684 2.1146 9.7592 4.2424 3.3368 20.4196 8.9422 6.4477 18.7062];

k=1;
A=zeros(12,1);
for k = 1:1:12
    if (k==6 || k==7 || k==10 || k==11 || k==12)
        %continue;
    end
    [ MUe, Ie ] = getKursRVectors( k );
    i=1;
    [ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(k);
    [ MUn, sigma, Tn ] = ABUparams( k );
    Tn=3;
    
    MUn=0.6;
    sigma = 10;
    if k==12
        MUn=1;
        sigma = 18;
    end
    
    %MUe(k) = 1;
    %Ie(k) = 1;
    %Ie(k) = 1;
    
    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;

    A0=A0_2;
    A1=A1_2;A2=A2_2;A3=A3_2;A4=A4_2;A5=A5_2;
    Tpsi = TP_r(k);
    A1 = A12_r(k);
    A2 = A22_r(k);
    A3 = A32_r(k);
    A4 = A42_r(k);
    A5 = A4_r(k);
    
    ie_psi=6;
    %Tpsi=b7/b4;
    a2_2 = a2+ (sigma*a3*a4)/(57.3*b4-sigma*a7);
    if k==1
        %continue
        ie_psi=3.5;
    end
    if (k==4 || k==9 || k==12)
        if (k==4 || k==9)
            ie_psi=5;
        end
        %i=1;
        %continue;
        if (k==4)
            i=0.6;
            %continue;
        end
        a2=a2_2;
    end
    
    w1 = get_tf_kurs_zamkn( Tn, Tpsi, ie_psi, i, a2, b3, b4, b7, A0, A1, A2, A3, A4, A5 );
    %A(k) = w1;

    step(w1)
    hold on
    %w2 = get_tf_kren_zamkn( mu, i, nu, b1, 2 )
    %step(w2)
end
axis([0 40 0 1.5])
grid on;