def test_orbital_radial():
r = np.linspace(0, 10, 1000)
f = np.exp(- r)
o = sisl.SphericalOrbital(2, (r, f))
sisl.plot(o)
fig = plt.figure()
sisl.plot(o, axes=fig.gca())
plt.close('all')
def test_orbital_harmonics():
r = np.linspace(0, 10, 1000)
f = np.exp(-r)
o = sisl.SphericalOrbital(2, (r, f))
sisl.plot(o, harmonics=True)
plt.close('all')