diff --git a/Examples/fit55_SpecularIntro/FitWithUncertainties.py b/Examples/fit55_SpecularIntro/FitWithUncertainties.py
index 6ff15c09a5b8f140cad5f519eef68b468acc793b..1980a03a1eabacd00eca656520657f6705e63244 100755
--- a/Examples/fit55_SpecularIntro/FitWithUncertainties.py
+++ b/Examples/fit55_SpecularIntro/FitWithUncertainties.py
@@ -12,12 +12,12 @@ Besides we will set Chi squared with L1-normalization as the objective metric
 and use genetic algorithm as the minimizer.
 """
 
-import bornagain as ba
-import ba_fitmonitor as bafim
-from bornagain import ba_plot
+import os
 import numpy as np
 from matplotlib import pyplot as plt
-import os
+import bornagain as ba
+from bornagain import ba_plot
+import ba_fitmonitor as bafim
 
 
 def get_sample(params):
@@ -53,7 +53,7 @@ def get_sample(params):
     substrate_layer = ba.Layer(m_substrate)
     multi_layer = ba.MultiLayer()
     multi_layer.addLayer(vacuum_layer)
-    for i in range(n_repetitions):
+    for _ in range(n_repetitions):
         multi_layer.addLayer(ti_layer)
         multi_layer.addLayer(ni_layer)
     multi_layer.addLayer(substrate_layer)
diff --git a/Examples/fit55_SpecularIntro/PolarizedSpinAsymmetryFit.py b/Examples/fit55_SpecularIntro/PolarizedSpinAsymmetryFit.py
index 31aff477bc8185842cec82a1d27063095298d885..ad5d18c445dd420f8f3108fee2251d70211ec219 100755
--- a/Examples/fit55_SpecularIntro/PolarizedSpinAsymmetryFit.py
+++ b/Examples/fit55_SpecularIntro/PolarizedSpinAsymmetryFit.py
@@ -9,12 +9,12 @@ For simplicity, here we only reproduce the first part of that
 demonstration without the magnetically dead layer.
 """
 
-import bornagain as ba
-from bornagain import deg, angstrom
-
-import numpy, os
+import os
 from sys import argv
+import numpy
 import matplotlib.pyplot as plt
+import bornagain as ba
+from bornagain import deg, angstrom
 
 # q-range on which the simulation and fitting are to be performed
 qmin = 0.05997
diff --git a/Examples/scatter2d/CorrelatedRoughness.py b/Examples/scatter2d/CorrelatedRoughness.py
index 3c6e7a0d3ecce192d5a512be4504af496cfc11b5..1a1382cf6fce408a19cb5178337f31d1179bdfd9 100755
--- a/Examples/scatter2d/CorrelatedRoughness.py
+++ b/Examples/scatter2d/CorrelatedRoughness.py
@@ -31,7 +31,7 @@ def get_sample():
     my_sample.addLayer(l_ambience)
 
     n_repetitions = 5
-    for i in range(n_repetitions):
+    for _ in range(n_repetitions):
         my_sample.addLayerWithTopRoughness(l_part_a, roughness)
         my_sample.addLayerWithTopRoughness(l_part_b, roughness)
 
diff --git a/Examples/scatter2d/MesoCrystal.py b/Examples/scatter2d/MesoCrystal.py
index 12d1a655c7be3109b8344c98c34a33d6a0fd5664..5edd036f9e6c81aba999dc30b4f05e9f306013ec 100755
--- a/Examples/scatter2d/MesoCrystal.py
+++ b/Examples/scatter2d/MesoCrystal.py
@@ -3,7 +3,7 @@
 Cylindrical mesocrystal on a substrate
 """
 import bornagain as ba
-from bornagain import deg, nm, R3
+from bornagain import deg, nm
 
 
 def get_sample():
diff --git a/Examples/specular/AlternatingLayers.py b/Examples/specular/AlternatingLayers.py
index ed9de46be6aee4aeb665b774fae30eee40c35823..95597dea773c225d0e6a93b84d0a4389fb9d735c 100755
--- a/Examples/specular/AlternatingLayers.py
+++ b/Examples/specular/AlternatingLayers.py
@@ -23,7 +23,7 @@ def get_sample():
     # Define sample
     sample = ba.MultiLayer()
     sample.addLayer(ambient_layer)
-    for i in range(10):
+    for _ in range(10):
         sample.addLayer(ti_layer)
         sample.addLayer(ni_layer)
     sample.addLayer(substrate_layer)
diff --git a/Examples/specular/BeamAngularDivergence.py b/Examples/specular/BeamAngularDivergence.py
index 098b17af940b40de9103ca14507417fee2f93f70..a67ae771468df1c02eb4ac8ce845f9f73a210205 100755
--- a/Examples/specular/BeamAngularDivergence.py
+++ b/Examples/specular/BeamAngularDivergence.py
@@ -45,7 +45,7 @@ def get_sample():
     substrate_layer = ba.Layer(m_substrate)
     multi_layer = ba.MultiLayer()
     multi_layer.addLayer(vacuum_layer)
-    for i in range(n_repetitions):
+    for _ in range(n_repetitions):
         multi_layer.addLayer(ti_layer)
         multi_layer.addLayer(ni_layer)
     multi_layer.addLayer(substrate_layer)
diff --git a/Examples/specular/BeamFullDivergence.py b/Examples/specular/BeamFullDivergence.py
index ff9cc38a9df9de2794b49930f135aa5ca71adefa..6281c40374938b15ddf0c7e57a7e4e24a7d95b52 100755
--- a/Examples/specular/BeamFullDivergence.py
+++ b/Examples/specular/BeamFullDivergence.py
@@ -44,7 +44,7 @@ def get_sample():
     substrate_layer = ba.Layer(m_substrate)
     multi_layer = ba.MultiLayer()
     multi_layer.addLayer(vacuum_layer)
-    for i in range(n_repetitions):
+    for _ in range(n_repetitions):
         multi_layer.addLayer(ti_layer)
         multi_layer.addLayer(ni_layer)
     multi_layer.addLayer(substrate_layer)
diff --git a/Examples/specular/TOFRWithResolution.py b/Examples/specular/TOFRWithResolution.py
index 1763e488fb608ffa4d221384909c0129d5aa1742..275eb2fa801d296edba7b097d92dcd7670648b15 100755
--- a/Examples/specular/TOFRWithResolution.py
+++ b/Examples/specular/TOFRWithResolution.py
@@ -35,7 +35,7 @@ def get_sample():
     # creating multilayer
     multi_layer = ba.MultiLayer()
     multi_layer.addLayer(ambient_layer)
-    for i in range(10):
+    for _ in range(10):
         multi_layer.addLayer(ti_layer)
         multi_layer.addLayer(ni_layer)
     multi_layer.addLayer(substrate_layer)
diff --git a/Examples/specular/TimeOfFlightReflectometry.py b/Examples/specular/TimeOfFlightReflectometry.py
index 72eb026a6f79fe695ecf6293c16bd8f81ce63611..31bed9fd273046b22d6cfb8e6dcb51964307d7de 100755
--- a/Examples/specular/TimeOfFlightReflectometry.py
+++ b/Examples/specular/TimeOfFlightReflectometry.py
@@ -33,7 +33,7 @@ def get_sample():
     # creating multilayer
     multi_layer = ba.MultiLayer()
     multi_layer.addLayer(ambient_layer)
-    for i in range(10):
+    for _ in range(10):
         multi_layer.addLayer(ti_layer)
         multi_layer.addLayer(ni_layer)
     multi_layer.addLayer(substrate_layer)
diff --git a/Examples/varia/FootprintCorrection.py b/Examples/varia/FootprintCorrection.py
index 56249203ae78615064e616abf56ff86768211322..e0a4f2e5a369ac6e94efa3bb93d4fa03ce06d15a 100755
--- a/Examples/varia/FootprintCorrection.py
+++ b/Examples/varia/FootprintCorrection.py
@@ -29,7 +29,7 @@ def get_sample():
     # creating multilayer
     multi_layer = ba.MultiLayer()
     multi_layer.addLayer(ambient_layer)
-    for i in range(10):
+    for _ in range(10):
         multi_layer.addLayer(ti_layer)
         multi_layer.addLayer(ni_layer)
     multi_layer.addLayer(substrate_layer)
diff --git a/Examples/varia/MaterialProfile.py b/Examples/varia/MaterialProfile.py
index 79fdf65e4aff081580800108c9b858e1624dc040..d1c079e807f07c55a2b263e69634586c406a6b0a 100755
--- a/Examples/varia/MaterialProfile.py
+++ b/Examples/varia/MaterialProfile.py
@@ -32,7 +32,7 @@ def get_sample():
     # creating multilayer
     multi_layer = ba.MultiLayer()
     multi_layer.addLayer(ambient_layer)
-    for i in range(4):
+    for _ in range(4):
         multi_layer.addLayerWithTopRoughness(ti_layer, roughness)
         multi_layer.addLayerWithTopRoughness(ni_layer, roughness)
     multi_layer.addLayer(substrate_layer)
diff --git a/Examples/varia/SpecularSimulationWithRoughness.py b/Examples/varia/SpecularSimulationWithRoughness.py
index 64bded655e2588a5a7f37edaf83bac382fdf3a89..c5483f7dafbb56af1d05eff520c44fdd61e0b775 100755
--- a/Examples/varia/SpecularSimulationWithRoughness.py
+++ b/Examples/varia/SpecularSimulationWithRoughness.py
@@ -30,7 +30,7 @@ def get_sample():
     # Sample
     multi_layer = ba.MultiLayer()
     multi_layer.addLayer(ambient_layer)
-    for i in range(10):
+    for _ in range(10):
         multi_layer.addLayerWithTopRoughness(ti_layer, roughness)
         multi_layer.addLayerWithTopRoughness(ni_layer, roughness)
     multi_layer.addLayerWithTopRoughness(substrate_layer, roughness)