def METHOD_NAME(year, city_code, expected_revenu_disponible):
    simulation = simulations.Simulation(
        period=periods.period(year), tax_benefit_system=tax_benefit_system
    )
    famille = simulation.populations["famille"]
    famille.count = 3
    famille.roles_count = 2
    famille.step_size = 1
    individu = simulation.populations["individu"]
    individu.count = 6
    individu.step_size = 2
    simulation.get_or_new_holder("city_code").array = np.array(
        [city_code, city_code, city_code]
    )
    famille.members_entity_id = np.array([0, 0, 1, 1, 2, 2])
    simulation.get_or_new_holder("salaire_brut").array = np.array(
        [0.0, 0.0, 50000.0, 0.0, 100000.0, 0.0]
    )
    revenu_disponible = simulation.calculate("revenu_disponible")
    assert_near(
        revenu_disponible, expected_revenu_disponible, absolute_error_margin=0.005
    )