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

2   
-lambda*(-mean + z)   
——————–   
____ 2   
___ _______ / 1 2*mean *z   
\/ 2 *\/ lambda * / — *e   
/ 3   
\/ z   
———————————————–   

____

2*\/ pi 

# Import sympy and Wald
from sympy.stats import Wald, density
from sympy import Symbol, pprint
  
z = 0.86
mean = 6
lambda = 2.35
  
# Using sympy.stats.Wald() method
X = Wald("x", mean, lambda)
gfg = density(X)(z)
  
pprint(gfg)

0.498668646362573   
—————–   

____

\/ pi