# Import sympy and Maxwell
from sympy.stats import Maxwell, density
from sympy import Symbol, pprint
  
z = Symbol("z")
a = Symbol("a", positive = True)
  
# Using sympy.stats.Maxwell() method
X = Maxwell("x", a)
gfg = density(X)(z)
  
pprint(gfg)

2   
-z   
—-   
2   
___ 2 2*a   
\/ 2 *z *e   
————–   
____ 3   
\/ pi *a 

# Import sympy and Maxwell
from sympy.stats import Maxwell, density
from sympy import Symbol, pprint
  
z = 1.2
a = 3
  
# Using sympy.stats.Maxwell() method
X = Maxwell("x", a)
gfg = density(X)(z)
  
pprint(gfg)

0.0492328718072872*\/ 2   
————————   

____

\/ pi