def estimate_fuse_passengers(fus_nb, FLOORS_NB, PASS_PER_TOILET, cabin_area,\
                             MASS_PASS, pass_density):
    """ The function evaluates the number of passengers members on board in
        case of an unconventional aircraft with fuselage.

    Source : passengers density is defined averaging the ratio between
    the maximum passengers allowed and the cabin area on different
    conventional aircraft.

    Args:
        fus_nb          (int): Number of fuselages.
        FLOORS_NB       (int): Number of floors.
        PASS_PER_TOILET (int): Number of passengers per toilet.
        cabin_area    (float): Cabin Area [m^2].
        MASS_PASS     (float): Passenger mass [kg].
        pass_density  (float): Base passengers density [pass/m^2].

    Returns:
        pass_nb       (int): Number of passengers.
        toilet_nb     (int): Number of toilet.
        mass_pass_tot (int): Total passengers mass [kg].

    """

    MF = FLOORS_NB * 0.5 + 0.5
    pass_nb = 0

    for i in range(0, fus_nb):
        pass_nb += cabin_area[i-1] * pass_density

    pass_nb = int(round(pass_nb * MF,0))
    mass_pass_tot = round(pass_nb * MASS_PASS,3)
    toilet_nb = round(pass_nb/PASS_PER_TOILET,0)

    return(pass_nb, toilet_nb, mass_pass_tot)