\frac{\partial}{\partial x^\alpha}\left(\sqrt{-g}\ g^{\alpha\beta} \frac{\partial \phi}{\partial x^\beta} \right) \right]+ \frac{m^2c^2}{\hbar^2}\phi = 0