Hi;
I'm trying to do a Monte Carlo electron scattering problem. Briefly what I have to do is start with a number of electrons (say 100), test the probability that each one scatters into one state or another (by testing a random number) and then place the electrons in the proper states. I have to do this for as many time intervals (of 1 picosecond) as it takes to get the initial number to get to zero.
So, No is the initial number of electrons, Nf is one scattering state based on the conditions in the code, and Ng is the other state. Here's the code I have so far:
tau1 = 30*(10^(-12));
tau2 = 60*(10^(-12));
deltatau = 1*(10^(-12));
No = 100;
Nf = 0;
Ng = 0;
Ni = No;
P1 = (deltat/tau1);
P2 = (deltat/tau2);
Psum = P1+P2;
r = rand(10,1)
for i=1
if (r(i)>0 & r(i)<P1)
Ni(i) = Ni-1
Nf(i) = Nf+1
end
if (r(i)>P1 & r(i)<Psum)
Ni(i) = Ni-1
Ng(i) = Ng+1
end
end
The thing I'm struggling with is how to loop through the values (hence why I only have i=1 in the for loop). So, if the first random number meets the condition, then I need to add 1 to either Nf or Ng, subtract 1 from Ni (so I'll have 99 Ni's left) and then go back to the beginning, look at the next random number and if that one meets the condition, add another value to Nf or Ng and subtract another value from Ni, until the value of Ni reaches zero.
It's a pretty convoluted problem, but it's basically a for loop within a for loop (I think) and I just can't set it up (my Matlab skills aren't that good). Any help is very much appreciated! Thanks!
