Yes, judging from a quick look at the code (and not being familiar with the GiNaC library used) this seems to square the wave function instead of squaring the a

Yes, judging from a quick look at the code (and not being familiar with the GiNaC library used) this seems to square the wave function instead of squaring the absolute value. [1] https://github.com/c0nrad/hydrogen/blob/master/hydrogen.cpp#L205 https://github.com/c0nrad/hydrogen/blob/master/hydrogen.cpp#...