diff --git a/.idea/workspace.xml b/.idea/workspace.xml
index 7e3bd25..e9d8af3 100644
--- a/.idea/workspace.xml
+++ b/.idea/workspace.xml
@@ -2,15 +2,8 @@
 <project version="4">
   <component name="ChangeListManager">
     <list default="true" id="f4439dc0-6756-4612-8f7d-596d8949f300" name="Default Changelist" comment="">
-      <change afterPath="$PROJECT_DIR$/data/csvs/may_pitzer_params.csv" afterDir="false" />
-      <change afterPath="$PROJECT_DIR$/data/csvs/zeroes_removed_PC88A_HCL_NdPrCeLaDySmY.csv" afterDir="false" />
-      <change afterPath="$PROJECT_DIR$/data/xmls/PC88A_HCL_NdPrCeLaDySmY_w_pitzer.xml" afterDir="false" />
-      <change afterPath="$PROJECT_DIR$/docs/Examples/iterative_fitter.py" afterDir="false" />
-      <change afterPath="$PROJECT_DIR$/llepe/objectives.py" afterDir="false" />
-      <change afterPath="$PROJECT_DIR$/llepe/optimizers.py" afterDir="false" />
       <change beforePath="$PROJECT_DIR$/.idea/workspace.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/workspace.xml" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/data/xmls/twophase.xml" beforeDir="false" afterPath="$PROJECT_DIR$/data/xmls/twophase.xml" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/llepe/__init__.py" beforeDir="false" afterPath="$PROJECT_DIR$/llepe/__init__.py" afterDir="false" />
+      <change beforePath="$PROJECT_DIR$/llepe/optimizers.py" beforeDir="false" afterPath="$PROJECT_DIR$/llepe/optimizers.py" afterDir="false" />
     </list>
     <option name="SHOW_DIALOG" value="false" />
     <option name="HIGHLIGHT_CONFLICTS" value="true" />
@@ -269,7 +262,14 @@
       <option name="project" value="LOCAL" />
       <updated>1594739749122</updated>
     </task>
-    <option name="localTasksCounter" value="14" />
+    <task id="LOCAL-00014" summary="added example file to do iterative optimization">
+      <created>1594767749286</created>
+      <option name="number" value="00014" />
+      <option name="presentableId" value="LOCAL-00014" />
+      <option name="project" value="LOCAL" />
+      <updated>1594767749286</updated>
+    </task>
+    <option name="localTasksCounter" value="15" />
     <servers />
   </component>
   <component name="Vcs.Log.Tabs.Properties">
@@ -293,7 +293,8 @@
     <MESSAGE value="Edited parity plot to allow color to represent 3rd dimension. Still need to improve colorbar axis name." />
     <MESSAGE value="updated to LLEPE" />
     <MESSAGE value="fixed setup file by changing LLEPE to llepe" />
-    <option name="LAST_COMMIT_MESSAGE" value="fixed setup file by changing LLEPE to llepe" />
+    <MESSAGE value="added example file to do iterative optimization" />
+    <option name="LAST_COMMIT_MESSAGE" value="added example file to do iterative optimization" />
   </component>
   <component name="WindowStateProjectService">
     <state x="-1330" y="212" key="#Inspections" timestamp="1590787654691">
@@ -313,11 +314,11 @@
       <screen x="-1920" y="2" width="1920" height="1040" />
     </state>
     <state x="-1213" y="379" key="ANALYSIS_DLG_com.intellij.analysis.BaseAnalysisAction$1/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1590787657711" />
-    <state x="169" y="115" key="CommitChangelistDialog2" timestamp="1594739747678">
+    <state x="169" y="115" key="CommitChangelistDialog2" timestamp="1594767747695">
       <screen x="0" y="0" width="1920" height="1040" />
     </state>
     <state x="-1364" y="117" key="CommitChangelistDialog2/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1592321237444" />
-    <state x="169" y="115" key="CommitChangelistDialog2/0.0.1920.1040@0.0.1920.1040" timestamp="1594739747678" />
+    <state x="169" y="115" key="CommitChangelistDialog2/0.0.1920.1040@0.0.1920.1040" timestamp="1594767747695" />
     <state x="-1828" y="94" width="1736" height="856" key="DiffContextDialog" timestamp="1591048879404">
       <screen x="-1920" y="2" width="1920" height="1040" />
     </state>
@@ -367,11 +368,11 @@
       <screen x="-1920" y="2" width="1920" height="1040" />
     </state>
     <state x="-1245" y="438" width="437" height="168" key="VCS.ChangelistChooser/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1592321346068" />
-    <state x="552" y="254" key="Vcs.Push.Dialog.v2" timestamp="1594739751769">
+    <state x="552" y="254" key="Vcs.Push.Dialog.v2" timestamp="1594767758886">
       <screen x="0" y="0" width="1920" height="1040" />
     </state>
     <state x="-1368" y="256" key="Vcs.Push.Dialog.v2/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1592321242949" />
-    <state x="552" y="254" key="Vcs.Push.Dialog.v2/0.0.1920.1040@0.0.1920.1040" timestamp="1594739751769" />
+    <state x="552" y="254" key="Vcs.Push.Dialog.v2/0.0.1920.1040@0.0.1920.1040" timestamp="1594767758886" />
     <state x="263" y="182" width="1283" height="717" key="com.intellij.ide.util.TipDialog" timestamp="1594694030135">
       <screen x="0" y="0" width="1920" height="1040" />
     </state>
diff --git a/llepe/optimizers.py b/llepe/optimizers.py
index 221a069..49d3646 100644
--- a/llepe/optimizers.py
+++ b/llepe/optimizers.py
@@ -1,6 +1,6 @@
 import scipy.optimize as scipy_opt
 from scipy.optimize import minimize
-import skopt
+# import skopt
 
 def dual_anneal_optimizer(objective, x_guess):
     bounds = [(1e-1, 1e1)] * len(x_guess)
@@ -21,19 +21,19 @@ def diff_evo_optimizer(objective, x_guess):
     return est_parameters, res.fun
 
 
-def forest_lbfgsb_optimizer(objective, x_guess):
-    x_guess = list(x_guess)
-    bounds = [(1e-1, 1e1)]*len(x_guess)
-    bounds[1] = (1e-1, 2)
-    res = skopt.forest_minimize(objective,
-                                bounds,
-                                random_state=1,
-                                acq_func='LCB',
-                                n_random_starts=30,
-                                x0=x_guess,
-                                xi=1e-4)
-    x_guess = res.x
-    optimizer_kwargs = {"method": 'l-bfgs-b',
-                        "bounds": bounds}
-    res = minimize(objective, x_guess, **optimizer_kwargs)
-    return res.x, res.fun
+# def forest_lbfgsb_optimizer(objective, x_guess):
+#     x_guess = list(x_guess)
+#     bounds = [(1e-1, 1e1)]*len(x_guess)
+#     bounds[1] = (1e-1, 2)
+#     res = skopt.forest_minimize(objective,
+#                                 bounds,
+#                                 random_state=1,
+#                                 acq_func='LCB',
+#                                 n_random_starts=30,
+#                                 x0=x_guess,
+#                                 xi=1e-4)
+#     x_guess = res.x
+#     optimizer_kwargs = {"method": 'l-bfgs-b',
+#                         "bounds": bounds}
+#     res = minimize(objective, x_guess, **optimizer_kwargs)
+#     return res.x, res.fun