Cosmetic changes only.
This commit is contained in:
parent
b3a6f61d75
commit
48463a9d0c
@ -50,20 +50,20 @@ function plotPhaseSpaceWithAccelerationField(obj, MaximumVelocity, NumberOfBins,
|
||||
%-------------------------------------------------------------------------
|
||||
|
||||
Y = linspace(0,MaximumVelocity,N);
|
||||
Trajectories = zeros(length(Y),int64(T/tau),9);
|
||||
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
|
||||
|
Loading…
Reference in New Issue
Block a user