def eval_poly_vs_fct(poly, function, test_values):
    """ Compute the maximum absolute difference between <poly> and
        <function> evaluated at each point in <test_values> vector """
    diff = max(abs(poly(v) - function(v)) for v in test_values)
    return diff