def METHOD_NAME(f_eval, eval_points, target_points, p):
    alpha_s_eval = eval_points
    beta_s = target_points
    U_dec = gen_Lagrange_coeffs(beta_s, alpha_s_eval, p)
    # print U_dec
    f_recon = np.mod((U_dec).dot(f_eval), p)
    # print f_recon
    return f_recon