"Roger Stafford" wrote in message <email@example.com>... > "Hemming" wrote in message <firstname.lastname@example.org>... > > Im trying to extract a scattering angle for a photon using the Klein-Nishina scattering angle distribution (KN in the code) and for this I need the CDF (of KN) to be able to use the Monte Carlo method when that is achieved. All i've managed so far is to plot the PDF between 0 degrees and Pi to see that it looks alright, and that it has that "peanut shape". Ive tried to use the built in CDF function but it seems very slow. > - - - - - - - - - - > Do I understand correctly that this is a probability density function with respect to a solid angle of scattering? Wouldn't that need to be known in obtaining a cumulative distribution? ....... - - - - - - - - - - One more thought about your question. If it is solid angle that the probability density is taken with respect to, then its cumulative distribution in terms of scatter angle is easy to find. However if it is this scattering you wish to simulate in a Monte Carlo process using matlab's 'rand' generator, you would need the inverse of this CDF function and that may not be so easy to calculate. There are some functions in matlab that may able to help in this. There is also the alternative of a random procedure using rejection to achieve the proper simulation distribution. Anyway please give us the details of what it is you wish to do.