F1 = sqrt( h1 )/ Rv1 ;
F12 = (1/ Rv12 )*sign( h1 - h2 )*sqrt(abs( h1 - h2 )) ;
f2 = (1/ Rv2 )*sqrt( h2 ) ;
f23 =( 1/ Rv23 )*sign( h2 - h3 )*sqrt(abs( h2 - h3 )) ;
F3 = (1/ Rv3 )*sqrt( h3 ) ;
d_h1 = Fin - F12 - F1 ;
d_h2 = F12 - f23 - f2 ;
d_h3 = f23 - F3 ;