def METHOD_NAME():
    mesh = IntervalMesh(10, 0, 1)
    c = Constant((3.0, 4.0))
    J = assemble(c**2*dx(domain=mesh))
    Jhat = ReducedFunctional(J, Control(c))
    deriv = Jhat.derivative()
    assert(taylor_test(Jhat, c, Constant((0.3, 0.5))) > 1.9)