diff --git a/MOT Capture Process Simulation/+Plotting/plotPhaseSpaceWithAccelerationField.m b/MOT Capture Process Simulation/+Plotting/plotPhaseSpaceWithAccelerationField.m index 7edb2c7..c68faae 100644 --- a/MOT Capture Process Simulation/+Plotting/plotPhaseSpaceWithAccelerationField.m +++ b/MOT Capture Process Simulation/+Plotting/plotPhaseSpaceWithAccelerationField.m @@ -49,21 +49,21 @@ function plotPhaseSpaceWithAccelerationField(obj, MaximumVelocity, NumberOfBins, %------------------------------------------------------------------------- - Y = linspace(0,MaximumVelocity,N); - Trajectories = zeros(length(Y),int64(T/tau),9); + Y = linspace(0,MaximumVelocity,N); + ParticleDynamicalQuantities = zeros(length(Y),int64(T/tau),6); for i=1:length(Y) x =-L/2; vx = Y(i)*cos(Theta); vz = Y(i)*sin(Theta); r = [x,0,z]; v = [vx,0,vz]; - [Trajectories(i,:,:),~] = obj.solver(r, v); + ParticleDynamicalQuantities(i,:,:) = obj.solver(r, v); end hold on for i=1:length(Y) - plot(Trajectories(i,:,1),Trajectories(i,:,4),'w','linewidth',1.3) + plot(ParticleDynamicalQuantities(i,:,1),ParticleDynamicalQuantities(i,:,4),'w','linewidth',1.3) end hold off