def METHOD_NAME():
    mesh = UnitCubeMesh(1, 1, 1)
    fs = FunctionSpace(mesh, 'CG', 1)
    f = Function(fs)
    out = Function(fs)
    u = TrialFunction(fs)
    v = TestFunction(fs)
    return inner(u, v) * dx, inner(f, v) * dx, out