def METHOD_NAME():
    mesh = UnitSquareMesh(10, 10)
    V = VectorFunctionSpace(mesh, "CG", 1)
    f = Function(V)
    c = Constant(2.0, domain=mesh)
    x, y = SpatialCoordinate(mesh)
    g = interpolate(as_vector([sin(y)+x, cos(x)*y]), V)
    f.assign(c * g)
    J = assemble(inner(f, f) ** 2  * dx)
    rf = ReducedFunctional(J, Control(c))
    h = Constant(0.1)
    r = taylor_to_dict(rf, c, h)
    assert min(r["R0"]["Rate"]) > 0.9
    assert min(r["R1"]["Rate"]) > 1.9
    assert min(r["R2"]["Rate"]) > 2.9