function T_uo = T_uo_ot_tau_cylinder_1(T_0, d, tau)
 % Определяет температуру конца ускоренного охлаждения T_uo [град. C]
 % круглого проката диаметром d [мм] с начальной температуры T_0 [град. C]
 % при длительности охлаждения tau [сек].

% При вычислениях будем использовать сплайн-интерполяцию трех зависимостей
load Cylinder_1_M_Spline.mat i_ot_T Tetta_ot_Fo A_ot_i

% Поскольку величина приведенного коэффициента температуропроводности a_pr
% зависит от конечной температуры, будем использовать метод последователь-
% ного приближения. Начинаем расчеты с величины приведенного коэффициента
% температуропроводности равного 5 мм.кв./с
a_pr_1=5;
a_pr=0.1;

while any(abs(a_pr_1-a_pr)./a_pr.*100 >= 0.1)
    % Расчеты производятся до тех пор, пока ошибка в определении a_pr
    % не будет ниже 0.1 %
    a_pr = a_pr_1;
    % Число Фурье при принятом нами a_pr
    Fo=a_pr.*tau./(d/2).^2;
    
    % Относительное теплосодержание в конце ускоренного охлаждения
    % как функция, зависящая от времени охлаждения (числа Фурье)
    Tetta_i = spline(Tetta_ot_Fo(1,:), Tetta_ot_Fo(2,:), Fo);
    
    % Теплосодержание проката в начале ускоренного охлаждения
    i_0 = spline(i_ot_T(1,:), i_ot_T(2,:), T_0);
    
    % Теплосодержание стали, соответствующее температуре поверхности
    % раската в процессе ускоренного охлаждения
    T_p = 40; % Принимаем температуру охлаждающей воды равной 40 град. C
    i_p = spline(i_ot_T(1,:), i_ot_T(2,:), T_p);
    
    % Теплосодержание раската в конце ускоренного охлаждения
    i_1 = Tetta_i.*(i_0 - i_p) + i_p;
    
    % Интегральная величина коэффициента температуропроводности для
    % теплосодержаний i_0 и i_1 составит
    A_0 = spline(A_ot_i(1,:), A_ot_i(2,:), i_0);
    A_1 = spline(A_ot_i(1,:), A_ot_i(2,:), i_1);
    
    % Приведенный коэффициент температуропроводности
    a_pr_1=(A_1-A_0)./(i_1-i_0);
end

% Температура конца ускоренного охлаждения как функция от теплосодержания
% стали
T_uo = spline(i_ot_T(2,:), i_ot_T(1,:), i_1);