%% Высоты каждого режима
[ H ] = H_vector();
[ M ] = M_vector();

%% Результаты из материалов Окоемова
% -------------------------------------------
% ( корректированные - ХЗ как посчитанные, неправильно, если использовать этот метод)
% Окоемов так и не показал каким образом получили эти числа
% качественно все сходится, в числах - нет!

% R - расчетные
% K - корректированные
% A - аппроксимированные
MU_R = [0.335, 0.033, 0, 0.737, 0.077, 0.354, 0.306, 0.319, 0.898, 0.735, 0.598, 1.995];
MU_K = [0.098, 0.033, 0, 0.241, 0.033, 0.144, 0.058, 0.048, 0.348, 0.243, 0.209, 0.566];
MU_A = [0.341, 0.341, 0.341, 0.341, 0.341, 0.341, 0.341, 0.341, 0.341, 0.341, 0.341, 0.566];

I_R = [1.534, 0.527, 0.806, 2.761, 1.179, 1.406, 1.588, 1.698, 3.047, 2.744, 2.250, 6.429];
I_K = [0.676, 0.527, 0.670, 0.838, 0.777, 0.660, 0.683, 0.698, 0.870, 0.826, 0.769, 1.087];
I_A = [0.838, 0.838, 0.838, 0.838, 0.838, 0.838, 0.838, 0.838, 0.838, 0.838, 0.838, 1.087];

NU_R = I_R;
NU_K = [0.527, 0.527, 0.527, 0.527, 0.527, 0.527, 0.527, 0.527, 0.527, 0.527, 0.527, 0.527];
NU_A = NU_K;


%для курса
Kurs_muMin = [0.1916, 0.1267, 0.1049, 0.3614, 0.2093, 0.3343, 0.2913, 0.3667, 0.5929, 0.4652, 0.5778, 0.9672];
Kurs_muMax = [0.8510, 1.9060, 1.0990, 2.0632, 1.3954, 1.5712, 1.5405, 1.9202, 2.4561, 2.0597, 2.6959, 3.9575];
Kurs_muMaxH = [0.9060, 0.9060, 0.9060, 2.0632, 1.3954, 1.5712, 1.5405, 1.9202, 2.4561, 2.0597, 2.6959, 3.9575];
Kurs_muAP = [0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 1.0];

% для АП крена
%[ MU, I, NU ] = getCorrectedVectors( )


% для АП курса:
t=2;
%[ Kurs_MU_R, Kurs_I_R ] = getKursRVectors( 2 ); % Tmin
%[ Kurs_MU, Kurs_I ] = getKursRVectors( 5 ); % Tmax

%% собственно строим графики тут
% ранжирование по высоте

% АП курса
% расчетные значения mu при Tmin и Tmax
%p=plot(H,Kurs_muMin,'Color','blue');
%hold on
%p=plot(H,Kurs_muMaxH,'Color','red');
%hold on
%p=plot(H,Kurs_muAP,'Color','green');
%set(p,'Color','green','LineWidth',2);
%grid on;

IEpsi1 = [3.3, 5.8, 3.8, 4.3, 4.9, 5.5, 5.9, 6.0, 4.4, 5.7, 6.0, 4.4];
IEpsi2 = [6.6, 11.6, 7.6, 8.6, 9.8, 11, 11.8, 12, 8.8, 11.4, 12, 8.8];
IEpsiAP = [6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6];
% ie psi
p=plot(H,IEpsi1,'Color','blue');
hold on
p=plot(H,IEpsi2,'Color','red');
hold on
p=plot(H,IEpsiAP,'Color','green');
set(p,'Color','green','LineWidth',2);
grid on;



[ MU1, I1, NU1 ] = getKrenRVectors( 2 )
[ MU2, I2, NU2 ] = getKrenRVectors( 5 )
% расчетные значения при Tmin.
%plot(H,MU1,'Color','blue');
%hold on
% 
%plot(H,MU2,'Color','red');
%hold on

%plot(H,MU_K,'Color','green');
%hold on

%p=plot(H,MU_A,'Color','green');
%grid on;


%p=plot(H,I_R,'Color','blue');
%hold on
%p=plot(H,I,'Color','red');
%hold on
%p2=plot(H,MU_A,'Color','black');
%grid on;
%set(p2,'LineWidth',2);



%set(p,'Color','red','LineWidth',2);

%hg