diff --git a/Quasi2DBogoliubovSpectrum.m b/Quasi2DBogoliubovSpectrum.m
index e21d5a3..0a24720 100644
--- a/Quasi2DBogoliubovSpectrum.m
+++ b/Quasi2DBogoliubovSpectrum.m
@@ -30,16 +30,16 @@ DyMagneticMoment             = 9.93*9.274009994E-24;
 
 %% Bogoliubov excitation spectrum for quasi-2D dipolar gas with QF correction
 AtomNumber        = 1E5;                                                                             % Total atom number in the system
-wz                = 2*pi*72.4;                                                                       % Trap frequency in the tight confinement direction
-lz                = sqrt(PlanckConstantReduced/(Dy164Mass*wz));                                      % Defining a harmonic oscillator length
-as                = 102.515*BohrRadius;                                                              % Scattering length
-Trapsize          = 7.5815;                                                                          % Trap is assumed to be a box of finite extent , given here in units of the harmonic oscillator length
+wz                = 2 * pi * 72.4;                                                                   % Trap frequency in the tight confinement direction
+lz                = sqrt(PlanckConstantReduced/(Dy164Mass * wz));                                    % Defining a harmonic oscillator length
+as                = 102.515 * BohrRadius;                                                            % Scattering length
+Trapsize          = 7.5815 * lz;                                                                     % Trap is assumed to be a box of finite extent , given here in units of the harmonic oscillator length
 alpha             = 0;                                                                               % Polar angle of dipole moment
 phi               = 0;                                                                               % Azimuthal angle of momentum vector
-MeanWidth         = 5.7304888515*lz;                                                                 % Mean width of Gaussian ansatz
+MeanWidth         = 5.7304888515 * lz;                                                               % Mean width of Gaussian ansatz
 k                 = linspace(0, 3e6, 1000);                                                          % Vector of magnitudes of k vector
 
-AtomNumberDensity = AtomNumber / (Trapsize * lz)^2;                                                  % Areal density of atoms
+AtomNumberDensity = AtomNumber / Trapsize^2;                                                         % Areal density of atoms
 add               = VacuumPermeability*DyMagneticMoment^2*Dy164Mass/(12*pi*PlanckConstantReduced^2); % Dipole length
 eps_dd            = add/as;                                                                          % Relative interaction strength
 gs                = 4 * pi * PlanckConstantReduced^2/Dy164Mass * as;                                 % Contact interaction strength