diff --git a/MOT Capture Process Simulation/+Plotting/plotFreeMolecularFluxVsTemp.m b/MOT Capture Process Simulation/+Plotting/plotFreeMolecularFluxVsTemp.m
index b2edac2..1900384 100644
--- a/MOT Capture Process Simulation/+Plotting/plotFreeMolecularFluxVsTemp.m	
+++ b/MOT Capture Process Simulation/+Plotting/plotFreeMolecularFluxVsTemp.m	
@@ -17,21 +17,22 @@ function plotFreeMolecularFluxVsTemp(obj, Temperature)
     
     for i=1:length(Temperature)
         beta = linspace(0.01,0.5,200);
-        L    = 2*obj.NozzleRadius./beta;
         obj.OvenTemperature = Temperature(i);
-        flux = zeros(1,length(L));
-        for j=1:length(L)
-            obj.NozzleLength = L(j);
-            flux(j) = obj.calculateFreeMolecularRegimeFlux();
+        flux = zeros(1,length(beta));
+        for j=1:length(beta)
+            obj.Beta = beta(j);
+            [ReducedClausingFactor, ~] = obj.calculateReducedClausingFactor();
+            flux(j)  = ReducedClausingFactor * obj.calculateFreeMolecularRegimeFlux();
         end
         plot(beta, flux, 'DisplayName', sprintf('T = %.1f ℃', Temperature(i)), 'Linewidth', 1.5)
     end
     set(gca,'yscale','log')
 
     obj.reinitializeSimulator();
-    
+    [ReducedClausingFactor, ~] = obj.calculateReducedClausingFactor();
+
     xline(obj.Beta, 'k--','Linewidth', 0.5);
-    fmf = obj.calculateFreeMolecularRegimeFlux();
+    fmf = ReducedClausingFactor * obj.calculateFreeMolecularRegimeFlux();
     yline(fmf, 'k--', 'Linewidth', 1.5);
     textstring = [sprintf('%1.e',fmf) ' atoms/s for ' '$$ \beta $$ = ' num2str(obj.Beta, '%.2f')  sprintf(' @ %.2f K', obj.OvenTemperatureinKelvin)];
     txt = text((obj.Beta + 0.05*obj.Beta), (max(fmf) + 0.2*fmf), textstring, 'Interpreter','latex');