def METHOD_NAME(self, p=1, n=1, plot=True):
    pde = XYData()
    mesh = pde.init_mesh(n=n)
    node = mesh.node
    cell = mesh.entity("cell")
    name = 'RobinBCTest.mat'
    space = LagrangeFiniteElementSpace(mesh, p=p)
    A = space.stiff_matrix()
    F = space.source_vector(pde.source)