Just for fun I'm doing a graphical representation of a QM problem out of Merzbacher. Basically we have a regular old potential barrier, with incoming particle of energy E>Vo. I am showing that as you increase the energy, the transmission is increased, i.e. the wavefunctions in all three areas match. Anyway..I'm not sure how to go about matching the wave functions at the boundary, I know all I have to do is change the phase of each wave fn but that's a pain, and I want this to be as automated as possible. Please see pdf of my plot, code can be shown if requested. Any ideas?