diff --git a/examples/coulomb_blockade/pentacene/hybr.py b/examples/coulomb_blockade/pentacene/hybr.py
index cfdb510..eae0d20 100755
--- a/examples/coulomb_blockade/pentacene/hybr.py
+++ b/examples/coulomb_blockade/pentacene/hybr.py
@@ -45,7 +45,7 @@ def hybridize_orbitals(
     gfp = ProjectedGreenFunction(gf, los_indices)
     hyb = Hybridization(gfp)
 
-    energies = np.arange(E_min, E_max + E_step / 2.0, E_step).round(7)
+    energies = np.linspace(E_min, E_max, int((E_max - E_min) / E_step) + 1)
 
     no = len(los_indices)
     gd = GridDesc(energies, no, complex)
diff --git a/examples/coulomb_blockade/pentacene/trans.py b/examples/coulomb_blockade/pentacene/trans.py
index 3cc8ac5..b8ef15c 100755
--- a/examples/coulomb_blockade/pentacene/trans.py
+++ b/examples/coulomb_blockade/pentacene/trans.py
@@ -45,7 +45,7 @@ def compute_transmission(
         eta=eta,
     )
 
-    energies = np.arange(E_min, E_max + E_step / 2.0, E_step).round(7)
+    energies = np.linspace(E_min, E_max, int((E_max - E_min) / E_step) + 1)
 
     i1 = los_indices - leads_nao
     s1 = hs_list_ii[1][1]