<jupyter_start><jupyter_text>Testing the Bessel Basis<jupyter_code>import dvr_1d
%matplotlib inline<jupyter_output><empty_output><jupyter_text>With a Morse potential<jupyter_code>d = dvr_1d.BesselDVR(npts=100, R=20, dim=3, lam=0)
d.morse_test(xmax=20., ymax=0.5)<jupyter_output>Testing 1-D DVR with a Morse potential
The first 5 energies are:
[-1.75751018 -0.76051377 -0.21815647 -0.00674606  0.05559304]<jupyter_text>With a Woods-Saxon potential<jupyter_code>d = dvr_1d.BesselDVR(npts=100, R=10, dim=3, lam=0)
d.woods_saxon_test()<jupyter_output>Testing 1-D DVR with a Woods-Saxon potential
The first 5 energies are:
[-43.55090359 -36.58134219 -29.25606017 -21.92813792 -14.94875837]
StarcoderdataJupyterStructuredCleanDedup