def METHOD_NAME(branch_pit, der_lambda, p_init_i_abs, p_init_i1_abs,
                                    height_difference):
    v_init_abs = np.abs(branch_pit[:, VINIT])
    v_init2 = v_init_abs * branch_pit[:, VINIT]
    lambda_term = np.divide(branch_pit[:, LENGTH] * branch_pit[:, LAMBDA], branch_pit[:, D]) \
                  + branch_pit[:, LC]
    const_p_term = np.divide(branch_pit[:, RHO], P_CONVERSION * 2)
    df_dv = const_p_term * (2 * v_init_abs * lambda_term + der_lambda
                            * np.divide(branch_pit[:, LENGTH], branch_pit[:, D]) * v_init2)
    load_vec = p_init_i_abs - p_init_i1_abs + branch_pit[:, PL] \
               + const_p_term * (GRAVITATION_CONSTANT * 2 * height_difference
                                 - v_init2 * lambda_term)
    mass_flow_dv = branch_pit[:, RHO] * branch_pit[:, AREA]
    df_dv_nodes = mass_flow_dv
    load_vec_nodes = mass_flow_dv * branch_pit[:, VINIT]
    df_dp = np.ones_like(der_lambda) * (-1)
    df_dp1 = np.ones_like(der_lambda)
    return load_vec, load_vec_nodes, df_dv, df_dv_nodes, df_dp, df_dp1