F(p, l)=\left(  \frac{1}{\rho^2} - \frac{1}{\rho_0^2}   \right)  
+ \left(  \frac{f}{\rho^2} + \frac{f_0}{\rho_0^2}   \right)  
\cdot \frac{l}{ 2 \, d}  - (2/j_m^2) \,  \int_p^{p_0} \frac{dp}{\rho} - (2/j_m^2) \, g \, \Delta z(l) = 0

\dot m =  
A \cdot \sqrt{ 2 \cdot \frac
{
g \, \Delta z + \int_p^{p_0} \frac{dp}{\rho}
}
{
\left( \frac{1}{\rho^2} - \frac{1}{\rho_0^2} \right) 
+ \left( \frac{f}{\rho^2} + \frac{f_0}{\rho_0^2} \right) 
\cdot \frac{l}{ 2 \, d}
}
}

q_0 =  
\frac{A}{\rho_s} \cdot \sqrt{ 2 \cdot \frac
{
g \, \Delta z + \int_p^{p_0} \frac{dp}{\rho}
}
{
\left( \frac{1}{\rho^2} - \frac{1}{\rho_0^2} \right) 
+ \left( \frac{f}{\rho^2} + \frac{f_0}{\rho_0^2} \right) 
\cdot \frac{l}{ 2 \, d}
}
}

j_m =  \sqrt{ 2 \cdot \frac
{
g \, \Delta z + \int_p^{p_0} \frac{dp}{\rho}
}
{
\left(  \frac{1}{\rho^2} - \frac{1}{\rho_0^2}   \right)  
+ \left(  \frac{f}{\rho^2} + \frac{f_0}{\rho_0^2}   \right)  
\cdot \frac{l}{ 2 \, d}
}
}