function tau=tau_cylinder_1(T_0, T_1, d)
 % Определяет время охлаждения tau [сек] круглого проката диаметром d [мм]
 % высокотурбулизированным потоком воды с температуры конца прокатки
 % T_0 [град. C] до температуры конца ускоренного охлаждения T_1 [град. C]

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

% Теплосодержание проката в начале и в конце ускоренного охлаждения
i_0 = spline(i_ot_T(1,:), i_ot_T(2,:), T_0);
i_1 = spline(i_ot_T(1,:), i_ot_T(2,:), T_1);

% Теплосодержание стали, соответствующее температуре поверхности раската
% в процессе ускоренного охлаждения
T_p = 40; % Принимаем температуру охлаждающей воды равной 40 ?C
i_p = spline(i_ot_T(1,:), i_ot_T(2,:), T_p);

% Относительное теплосодержание составит
Tetta_i = (i_1 - i_p)./(i_0 - i_p);

% Числу Фурье, соответствующее такому относительному теплосодержанию
Fo = spline(Tetta_ot_Fo(2,:), Tetta_ot_Fo(1,:), Tetta_i);

% Интегральная величина коэффициента температуропроводности для исходных
% теплосодержаний 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=(A_1-A_0)./(i_1-i_0);

% Время охлаждения составит
tau = Fo.*(d/2).^2./a_pr;