diff --git a/qpms/qpms_p.py b/qpms/qpms_p.py
index e55b50c..ecd4ceb 100644
--- a/qpms/qpms_p.py
+++ b/qpms/qpms_p.py
@@ -1,6 +1,7 @@
 import numpy as np
 from qpms_c import *
 ň = np.newaxis
+import scipy
 from scipy.constants import epsilon_0 as ε_0, c, pi as π, e, hbar as ℏ, mu_0 as μ_0
 eV = e
 from scipy.special import lpmn, lpmv, sph_jn, sph_yn, poch