diff --git a/src/pynlo/interactions/FourWaveMixing/SSFM.py b/src/pynlo/interactions/FourWaveMixing/SSFM.py
index 91cc905..091a04e 100755
--- a/src/pynlo/interactions/FourWaveMixing/SSFM.py
+++ b/src/pynlo/interactions/FourWaveMixing/SSFM.py
@@ -524,7 +524,7 @@ def propagate(self, pulse_in, fiber, n_steps, output_power=None, reload_fiber_ea
         
 
         for i in range(n_steps):                        
-            print ("Step:", i, "Distance remaining:", fiber.length * (1 - np.float(i)/n_steps) )
+            print ("Step:", i, "Distance remaining:", fiber.length * (1 - float(i)/n_steps) )
             
             if reload_fiber_each_step:
                 self.load_fiber_parameters(pulse_in, fiber, output_power, z=i*delta_z) 
diff --git a/src/pynlo/media/fibers/calculators.py b/src/pynlo/media/fibers/calculators.py
index ba1260d..fad816e 100644
--- a/src/pynlo/media/fibers/calculators.py
+++ b/src/pynlo/media/fibers/calculators.py
@@ -24,7 +24,7 @@
 
 
 import numpy as np
-from scipy.misc import factorial
+from math import factorial
 from scipy import constants
 import matplotlib.pyplot as plt
 
diff --git a/src/pynlo/media/fibers/fiber.py b/src/pynlo/media/fibers/fiber.py
index 1471be8..f3060c4 100755
--- a/src/pynlo/media/fibers/fiber.py
+++ b/src/pynlo/media/fibers/fiber.py
@@ -24,7 +24,7 @@
 import numpy as np
 import scipy.interpolate
 from pynlo.media.fibers.calculators import DTabulationToBetas
-from scipy.misc import factorial
+from math import factorial
 from scipy import constants
 from scipy.optimize import minimize
 from pynlo.util.pynlo_ffts import IFFT_t