如前所述，我需要用 Dirichlet BC 数值求解 Bunimovich 体育场形无限势阱中的薛定谔方程（这在分析上是不可能的）。为此，我需要定义系统的离散哈密顿矩阵。由于在任何坐标系中都没有很好地定义形状，我不知道该怎么做。

通常（就像我对方形井所做的那样），我会定义井的大小，然后使用 MATLAB eig 函数为每个维度求解 1D TISE。可以类似解决吗？

任何帮助将不胜感激。