[RowdyScot]

Arg...well, I've started out by deriving the Schroedinger equation in order to get energy. Solving through all of that is how I get the value in my above post of E = (n + 1/2)*hbar*omega. However, I'm not sure how I should be getting probability. Wait...maybe I've just realized something. If I derived the Schroedinger equation to get E, I should just be coding for ∫[(n + 1/2)*hbar*omega], right? How do I go about setting that up with limits from -x to x, which I can just say are -1 and 1, respectively? The whole integration in Python still throws me a bit.