Definition at line 123 of file numerov_solve.py. 00123 : 00124 "generate an eigenfunction with given potential with emin<(2m/hbar**2)E<emax, assuming the potential well has a single allowed region bounded on both sides" 00125 e=zbrent(lambda x: single_well_numerov_match(V0-x, dx, parity, overlapsize=overlapsize)[1], emin, emax, 1e-15, trace=0) 00126 psi, slope_error=single_well_numerov_match(V0-e, dx, parity, overlapsize) 00127 return e, psi 00128 def find_spectrum(V0, dx, levels, ground_state_estimate, state_separation_estimate=None, trace=0, subdivide=4):
|