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

2   
-(-mean + log(z))   
——————-   
2   
___ 2*std   
\/ 2 *e   
————————–   
____
2*\/ pi *std*z 

# Import sympy and LogNormal
from sympy.stats import LogNormal, density
from sympy import Symbol, pprint
  
z = 2.1
mean = 7.6
std = 4
  
# Using sympy.stats.LogNormal() method
X = LogNormal("x", mean, std)
gfg = density(X)(z)
  
pprint(gfg)

0.0136890249307238*\/ 2   
————————   
____
\/ pi