human
I wrote the below code

    def f(x, y):
        return ((x ** 2) + (y ** 2))
    actual = discretize_model(f, ((-5), 6), ((-5), 6))
    (y, x) = (np.indices((11, 11)) - 5)
    desired = ((x ** 2) + (y ** 2))
    assert_allclose(actual, desired)

What's a good function header?
 Here's a good function header?
assistant
#Here's a good function header.

def test_discretize_callable_2d():