Идеальный газ в распределенных коробчатых вероятностях

Question

У меня есть 2D идеальный газ в симуляции разделенной коробки. Все частицы газа изначально находятся на левой стороне коробки. Затем они начинают постоянно перемещаться между сторонами с течением времени.

Это мой код:

N=20; %total number of particles
nstep = 100;
n = randi(nstep,1);

n(1) = N; % initial conditions- number of particles on left side
for i = 2:nstep_5
r = rand(1,1);
if (r<n(i-1)/N)
n(i) = n(i-1) - 1; % Move atom from left to right
else
n(i) = n(i-1) + 1; % Move atom from right to left
end
end
time=(1:nstep)

Мне нужно, чтобы они вели себя, так как есть вероятность 60% двигаться слева направо и вероятность 40% для перемещения справа налево.

Я пытаюсь вставить вероятностные функции в этот код. Есть идеи, пожалуйста?

Answer 1

вы не проверяете состояние, когда частицы уже находятся на левой стороне и снова попали в go влево (то есть, чтобы остаться там, где они есть? Или сталкиваются со стеной, а go вправо? Нет Чисто). Что касается условия if в вашем l oop, вы можете написать его следующим образом:

if rand(1)<=0.6
   % move left to right
else
    % move right to left
end