hello
can anyone please guide me to good resources about this subject. i understand the theory but i cant understand where can i find an algorithm to base my work on it, also i am wondering if i can use the probability to control the bacteria movements and how, i need to do a program with the same idea but much simpler

hello
can anyone please guide me to good resources about this subject. i understand the theory but i cant understand where can i find an algorithm to base my work on it, also i am wondering if i can use the probability to control the bacteria movements and how, i need to do a program with the same idea but much simpler

Much better now. Problems I noted previously are repaired. I also did not see a truly n-d normal generator after another check, only the 2-d case. So this is more general.

Why would you use the central limit theorem to generate approximate normal variates, when randn does so much more accurately, as well as doing so much more quickly? I'm sorry, but that is just a downright ridiculous thing to do.
Looking at the help, it is pretty good. There is a reference, there are readable descriptions of the parameters, etc. I found error checks for the errors one might expect to see.
One issue is the H1 line. There was none. The H1 line is the very FIRST line of comment in the help block. The H1 line is used by lookfor, as well in creating a contents file when necessary.
The funny thing is, the author has written a very serviceable H1 line, but left it as the second comment line, with a blank comment line preceding it. This is silly, since it fails to serve its proper utility.
There are other issues, more grievous ones, I am afraid. The code computes the cholesky factor of the covariance matrix, but for no good reason, since then it next computers the eigenvalue decomposition. Since either of these decompositions will serve to transform a set of independent normal variates into the desired structure, the code does extra work for no valid reason.
Finally, there are at least a few other multivariate normal deviate generation codes on the file exchange. I don't see anything that this improves upon. I've given it as high as a 2 rating only because the help was clearly written.

Comment only