diff --git a/Bogoliubov-deGennes-Solver/solveBogoliubovdeGennesIn2D.m b/Bogoliubov-deGennes-Solver/solveBogoliubovdeGennesIn2D.m
index 11e1757..8458d8d 100644
--- a/Bogoliubov-deGennes-Solver/solveBogoliubovdeGennesIn2D.m
+++ b/Bogoliubov-deGennes-Solver/solveBogoliubovdeGennesIn2D.m
@@ -56,7 +56,7 @@ end
 % Obtain f from g
 f = zeros(size(g));
 for ii = 1:Neigs
-    f(:,ii) = muHC(g(:,ii)) + (2.*X(g(:,ii)));
+    f(:,ii) = (1/evals(ii)) * (muHC(g(:,ii)) + (2.*X(g(:,ii))));
 end
 
 % Obtain u and v from f and g