-- Maxwell-Boltzmann velocity distribution

fMB = (m / (2 pi k T))^(3/2) exp(-m v^2 / (2 k T))
fMB

-- Maxwell speed distribution (integrate over theta, then phi)

fM = defint(fMB v^2 sin(theta), theta, 0, pi, phi, 0, 2 pi)
fM

"Verify integral"

check(fM == 4 pi v^2 fMB)

"ok"
