function [P_in, v_otn, Q] = ...
    P_forsun(d_pr, v_pr, L_k, d_k, v_otn_min, Q_min)
% Определяет давление воды перед форсункой P_in [кгс/см2], относительную
% скорость потока v_otn [м/с] и расход воды Q [л/с] при диаметре проката
% d_pr [мм], скорости прокатки v_pr [м/с], длине и диаметре камеры
% охлаждения L_k [м] и d_k [мм].
% Необходимо создать относительную скорость движения потока v_otn_min [м/с]
% и обеспечить прохождение требуемого количества воды Q_min [л/с] при
% настройке форсунки на полное заполнение.

Lambla_k = 0.015; % Коэффициент сопротивления по длине камеры охлаждения
Lambla_r = 0.075; % Коэффициент сопротивления по длине раската

% Миллиметры диаметров переводим в метры
d_k = d_k./1000;
d_pr = d_pr./1000;
% Литры переводим в м3
Q_min = Q_min./1000;

% Фактическая площадь сечения потока
S_pot =pi./4.*(d_k.^2 - d_pr.^2);

% Скорость потока абсолютная и относительная для обеспечения требуемого
% расхода воды
v_pot = Q_min./S_pot;
v_otn = v_pot - v_pr;

% Относительная скорость потока обеспечивающая требуемый расход воды Q и
% заданную в условиях v_otn_min
v_otn = max(v_otn, v_otn_min);
% Скорость потока и расход воды для такой относительной скорости
v_pot = v_pr + v_otn;
Q = v_pot.*S_pot;

% Гидравлический радиус камеры охлаждения и раската
R_gidr_k = S_pot./pi./d_k;
R_gidr_r = S_pot./pi./d_pr;

% Общий коэффициент потерь от трения охлаждающей воды о стенки камеры
% охлаждения и о поверхность раската
Xi_k = Lambla_k.*L_k./4./R_gidr_k;
Xi_r = (v_otn./v_pot).^2.*Lambla_r.*L_k./4./R_gidr_r;

% При полном заполнении камеры охлаждения коэффициент потерь напора на
% формирование потока в форсунке составит
Xi_f = (1 - 1./sqrt(1 + Xi_k + Xi_r)).^2;

% Требуемый напор и давление воды перед форсункой
h_in = v_pot.^2./2./9.8.*(1 + Xi_f + Xi_k + Xi_r);
P_in = h_in./10;

% Для красивости вновь переводим м3 в литры
Q = Q.*1000;