procedure bernoulli_diffuzio;
begin
  (i,j):=veletlenhely;
  if a[i]='a' then dec(abd);
  if a[i]='b' then inc(adb);
  csereab(i,j);
end;

procedure sikbeli_diff_modell;
{"itt vonzo-ero van a molekulak kozott,
de ez mommentan nem a tavasznak koszonheto" - abonyita}
{-> a vonzoero a tavolsag negyzetevel forditottan ara'nyos}
{-> nem lephet molekula oda, ahol kevesebb szomszadja lenne}
begin
  (i,j):=veletlenhely;  (n,m)-es tablaban
  (k,l):=veletlenszomszed(i,j);
  if T(i,j)=1 and T(k,l)=0 and szomszedszam(i,j)<=(k,l)
       then csere(t(i,j),t(k,l));
end;