Dear Hamza Alhasan

These links might be useful, have a look:

1. https://stackoverflow.com/questions/63300833/maxwellian-distribution-in-python-scipy

2. https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.maxwell.html

3. https://scipython.com/blog/the-maxwellboltzmann-distribution-in-two-dimensions/

4. https://www.geeksforgeeks.org/python-maxwell-distribution-in-statistics/

Kind Regards

Qamar Ul Islam

may this link usefull for you:

https://stackoverflow.com/questions/63300833/maxwellian-distribution-in-python-scipy/63311101#63311101

Thank you... Useful links. I found the code:

from scipy.stats import maxwell

maxwell.rvs(loc=0, scale=np.sqrt(kB*T/m), size=Number_of_particles)