diff --git a/Bogoliubov-deGennes-Solver/solveBogoliubovdeGennesIn2D.m b/Bogoliubov-deGennes-Solver/solveBogoliubovdeGennesIn2D.m
index cc26b1e..84a3552 100644
--- a/Bogoliubov-deGennes-Solver/solveBogoliubovdeGennesIn2D.m
+++ b/Bogoliubov-deGennes-Solver/solveBogoliubovdeGennesIn2D.m
@@ -1,11 +1,11 @@
-function [evals, modes] = solveBogoliubovdeGennesIn2D(psi, Params, VDk, VParams, muchem)
+function [evals, modes] = solveBogoliubovdeGennesIn2D(psi, Params, VDk, VParams, Transf, muchem)
 
 wz_tilde     = Params.wz / Params.w0;
 gs           = Params.gs; 
 gdd          = Params.gdd;
 gammaQF      = Params.gammaQF; 
 
-KEop         = 0.5*KX.^2 + 0.5*KY.^2;
+KEop         = 0.5*(Transf.KX.^2+Transf.KY.^2);
 Ez           = (0.25*VParams.sigma^2) + (0.25*wz_tilde^2*VParams.sigma^2);
 muchem_tilde = muchem - Ez;
 
@@ -31,7 +31,7 @@ X            = @(w,psi) (psi.*real(ifft(VDk.*fft(psi.*w)))) + (3/2)*(gQF_pf_2D*g
 BdG          = @(g) muHC(muHC(g) + (2.*X(g)));
 
 syssize      = size(psi);
-opts.v0      = psi(:,:);
+opts.v0      = psi(:);
 opts.tol     = 1e-16;
 opts.disp    = 1;
 opts.issym   = 0;