diff --git a/EstimatesForAccordionLattice.m b/EstimatesForAccordionLattice.m
index 33dfbaa..ca7a3d2 100644
--- a/EstimatesForAccordionLattice.m
+++ b/EstimatesForAccordionLattice.m
@@ -173,20 +173,18 @@ TwoPhotonRecoilEnergy           = (2*PlanckConstantReduced*2*pi/Wavelength)^2 /
 TrapDepthsInUnitsOfRecoilEnergy = TrapDepth ./ TwoPhotonRecoilEnergy;
 RabiOscillationFrequency        = (1/PlanckConstantReduced) .* (sqrt(TrapDepth.^2/2 + TwoPhotonRecoilEnergy^2));
 
-PulseDurations                  = linspace(1E-6, 150E-6, 1000);
-
 alpha = TwoPhotonRecoilEnergy / PlanckConstantReduced;
 beta  = TrapDepth / PlanckConstantReduced;
 C     = beta^2 / ((2*beta^2)  + (4*alpha^2));
 
-RabiOscillations          = C .* (sin(0.5 .* PulseDurations .* (RabiOscillationFrequency)));
+PulseDurations                  = linspace(1E-6, 150E-6, 1000);
+PopulationInFirstOrders         = C .* (sin(0.5 .* PulseDurations .* (RabiOscillationFrequency)));
 
 figure(6);
 set(gcf,'Position',[100 100 950 750])
-plot(PulseDurations .* 1E6, RabiOscillations, LineWidth=2.0, DisplayName=['\bf Power = ' num2str(Power) ' W / Trap depth = ' num2str(round(TrapDepthsInUnitsOfRecoilEnergy, 1)) ' E_r'])
+plot(PulseDurations .* 1E6, PopulationInFirstOrders, LineWidth=2.0, DisplayName=['\bf Power = ' num2str(Power) ' W / Trap depth = ' num2str(round(TrapDepthsInUnitsOfRecoilEnergy, 1)) ' E_r'])
 xlabel('Pulse duration (µs)', FontSize=16)
 ylabel('Fraction of atoms in first order', FontSize=16)
 title('\bf Expected Rabi oscillation', FontSize=16)
 grid on
 legend(FontSize=16)
-