def METHOD_NAME():
    m = UnitSquareMesh(2, 2)
    V = FunctionSpace(m, 'CG', 1)
    METHOD_NAME = Function(V, name="f")
    x = SpatialCoordinate(m)
    METHOD_NAME.interpolate(x[0]*x[1])
    return METHOD_NAME