Solving Schrodinger's Equation Numerically in a Bunimovich Stadium

I need to solve, as mentioned, Schrodinger's equation in a Bunimovich stadium-shaped infinite potential well with Dirichlet BC Numerically (this isn't possible analytically). In order to do so, I need to define the discrete Hamiltonian matrix of the system. Since the shape isn't defined well in any coordinates system, I have no idea how to do so.

Normally (as I did with a square well), I'd define the size of the well and then solve 1D TISE for each dimension using MATLAB eig function. Can it be solved similarly?

Any help will be greatly appreciated.

• Can you describe the equations and the form of the potential you want to solve for? This may allow people to understand the question who are not already familiar with the Bunimovich Stadium. – Wolfgang Bangerth Apr 28 '16 at 14:32
• I don't see how is different from this question. You can use Finite Differences or Finite Methods for that purpose. BTW, the drawing you are referring to is really confusing. – nicoguaro Apr 28 '16 at 19:50
• Not different at all. In the approaches I mentioned (and in the codes) the potentials are turn into a discrete representation. – nicoguaro Apr 28 '16 at 21:06
• You still discretize the 2D domain, and the row and column numbers correspond to your discretization process. It depends if you use Finite Differences or Finite Elements but in essence is the same. – nicoguaro Apr 28 '16 at 23:15
• That's right. You replace your domain by a set of triangles or quadrilaterals, and then you just do the regular finite element approach on it (or finite differences, if you're so inclined). – Wolfgang Bangerth Apr 29 '16 at 2:45