Corrected spherical bessel definitions

charged_shells/functions.py

@@ -7,7 +7,7 @@ def sph_bessel_i(n, x, **kwargs):
 
 
 def sph_bessel_k(n, x, **kwargs):
-    return sps.spherical_kn(n, x, **kwargs)
+    return 2 / np.pi * sps.spherical_kn(n, x, **kwargs)
 
 
 def sph_harm(l, m, theta, phi, **kwargs):
@@ -15,7 +15,7 @@ def sph_harm(l, m, theta, phi, **kwargs):
 
 
 def interaction_coef_C(l, p, x):
-    return x * sps.iv(l + 1 / 2, x) * sps.iv(p + 1 / 2, x)
+    return np.pi / 2 * x * sps.iv(l + 1 / 2, x) * sps.iv(p + 1 / 2, x)
 
 
 def coefficient_Cpm(l, x):