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

% 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;

% параметры АБУ, тут это не надо
for k = 1:1:12
[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(k);
%MU_K(k)=(18-b1*t_reg_max)/(b3*t_reg_max);
sigma = 57.3*b4*(4-(a1*a4+a2))/(4*a7+a3*a4-a2*a7);
a2_2 = a2+ (sigma*a3*a4)/(57.3*b4-sigma*a7);
a4_2 = a4*57.3*b4/(57.3*b4-sigma*a7);
mun = 2*0.4*(sqrt(a1*a4+a2)-a1-a4)/a3;
%fprintf('%d\n', k);
%fprintf('x=%0.1f\n',a1+a4_2);
%fprintf('y=%0.1f\n',a1*a4_2+a2_2);
%fprintf('y=%0.1f\n',a1*a4_2+a2_2);
%fprintf('mun=%0.2f\n',mun);
%fprintf('\n');
end