diff --git a/flip/__init__.py b/flip/__init__.py
index 466843f..9a6f6a7 100644
--- a/flip/__init__.py
+++ b/flip/__init__.py
@@ -6,6 +6,8 @@
 from . import likelihood
 from . import model_evaluation
 from . import utils
+from . import power_spectra
+
 
 __version__ = "1.0.0"
 __flip_dir_path__ =  os.path.dirname(__file__)
diff --git a/flip/fitter.py b/flip/fitter.py
index b4511e6..e6c482d 100644
--- a/flip/fitter.py
+++ b/flip/fitter.py
@@ -280,6 +280,7 @@ def run(
         migrad=True,
         hesse=False,
         minos=False,
+        n_iter=1
     ):
         """
         The run function is the main function of the class. It takes in a number of
@@ -297,7 +298,9 @@ def run(
 
         """
         if migrad:
-            log.add(self.minuit.migrad())
+            for i in range(n_iter):
+                log.add(f'Iteration {i}/{n_iter}\n')
+                log.add(self.minuit.migrad())
         if hesse:
             log.add(self.minuit.hesse())
         if minos: