def heigth_max_suction(Pa, rho_F_20, g, P_satur_vapor_feed, w_liq_real_enter_feed, hydraulic_losses_suct_feed, heigth_cavitation_feed):
    """
    Calculates the maximum theoretical suction height.
    Parameters
    ----------
    Pa : float
        The atmosphere pressure [Pa]
    rho_F_20 : float
        The density of feed, [kg/m**3]
    P_satur_vapor : float
        The  pressure of saturated vapor, [Pa]
    w_liq_real_enter_feed : float
        The speed of suction line , [m/s]
    hydraulic_losses_suct_feed : float
        The hydraulic losses of suction line, [m]      
    heigth_cavitation_feed : float
        The losses due to cavitation, [m]  
    Returns
    -------
    heigth_max_suction : float
        The maximum theoretical suction height, [m]
    References
    ----------
    &&&&
    """   
    return ((Pa/(rho_F_20 * g) - ((P_satur_vapor_feed)/(rho_F_20 * g) + ((w_liq_real_enter_feed / (2 * g))) + hydraulic_losses_suct_feed + heigth_cavitation_feed)))