def METHOD_NAME(fieldstype=3, n=10, L=12, options=None):
    quadtree = quadmesh(n)
    NN = quadtree.number_of_nodes()
    print('The number of mesh:', NN)
    mesh = quadtree.to_pmesh()
    obj = SCFTVEMModel(mesh, options=options)
    mu = obj.init_value(fieldstype=fieldstype)  # get initial value
    problem = {'objective': obj, 'x0': mu, 'quadtree': quadtree}
    return problem