Hi every body,(adsbygoogle = window.adsbygoogle || []).push({});

i've a simple question in the programation for a certain methode called Monte Carlo method for Smoluchowski equation:

i programe this method in maple but the problem is that the positions of particules z[l]is not changing ( it most be changed ): can any one tell me where is the problem?

this is the method:

//calculate natural random numbers

test:=proc(N)

local a,i,L;

L:=[];

a:=rand(1..N):

for i from 1 to N do L:=[op(L),a()] od;

L;

end:

\\z is the position of particules

z:=test(10);

pi:=test(10);

\\ calculate uniform random numbers U[0,1]

r:=stats[random, uniform](10):for i from 1 to 10

do

printf("r[%d]=%f\n",i,r);

od;

\\ position of particules at time t+dt=t+1/20

for i from 1 to 10

do

if

r>1/(20*z) or r=1/(20*z)

then

z(t+dt)=z:

end if;

if

r<1/(20*z)

then

z(t+dt)=z+z[pi]:

end if;

od:

thankx verry much

**Physics Forums | Science Articles, Homework Help, Discussion**

Dismiss Notice

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Monte Carlo method for Smoluchowski equation

**Physics Forums | Science Articles, Homework Help, Discussion**