%disp('a'' = 2(c1+c5-c4) / c3;');
%a1=
%disp('b'' = [(c1+c4+c5)^2 - 4(c1*c4 + c2)] / c3^2;');

C=[ 0.82  11  8.8  1.2  0.3   7.3
    0.3   3   4.5  0.3  0.08  7.5];

IV=[ 2.1  1.8  1.05  1.75  1.6 ];
MUV=[3.0 2.6 1.5 2.5 2.25];


rezim=1;

c1=0.5;     % 0.1 .. 1.8
c2=4;       % -0.05 .. 8.0
c3=2;       % 0.3 .. 4.0
c4=0.1;     % 0.08 .. 2
c5=0.3;     % 0 .. 0.7
c6=8;       % 1.3 .. 16
c7=0.1712;  % 0.1712
c8=0.4;     % 0.005 .. 0.8
c9=0.1;     % 0.014 .. 0.25

%Tp=sqrt(1/(C(1,2)+C(1,4)*C(1,1)));
%kzeta=(C(1,1)+C(1,4)+C(1,5)) / 2 / Tp;
Tp=0.001;
kzeta=0.7;

c1=C(rezim,1);c2=C(rezim,2);
c3=C(rezim,3);c4=C(rezim,4);
c5=C(rezim,5);c6=C(rezim,6);
iv=IV(rezim);
muv=MUV(rezim);

walfa=sqrt(c1*c4+c2);
Ta=1/walfa;
Tv=1/c4;

ke=6;
A2ee=ke/(0.95/Tv);
A2e=A2ee+ke*Tv;

A0e=A2ee*Tv*0.9; %A0e<
%A2ee/A0e
%1/Tv
A1e=2*sqrt(A0e*A2ee);