Determination of Non-adiabatic Scattering Wave Functions in a Born-Oppenheimer Model