Say I want to randomly sample the Maxwell Boltzmann Distribution in python: This

Question

Say I want to randomly sample the Maxwell Boltzmann Distribution in python:

This function is going from 0 to infinity. In python, I would usually generate a random point from a function like:

def f(x):

return somenumber

xval = numpy.random.uniform(xmin,xmax)

randomyval = f(xval)

But how would I do this for a function with an infinite domain? I can't just choose an arbitrarily large xmax. I think I need to use importance sampling, but I'm not sure how. Can someone show me in python how to randomly sample a distribution with infinte range, where I'm still able to tweak the standard deviation if I like?

I've tried using scipy.stats.maxwell, but it doesn't give you an option to change the standard deviation, only scale and shape. And I'm not sure how they relate.

Explanation / Answer

Since there is a maxInt in Python ( 231-1 ) or (263 -1) depending on 32 and 64 bit systems , or in any computer languages , you cannot choose numbers beyond that and hence domain cannot be infinite.

The best you can do that choose numbers between minInt and maxInt .

*** Hope it helps , if not , we can follow this up in the comments ***

In Python 3 you can use the following code

Related Questions

Navigate

• Browse (All)
• Browse S
• Subjects

Integrity-first tutoring: explanations and feedback only — we do not complete graded work. Learn more.