# Import sympy and betanoncentral
from sympy.stats import BetaNoncentral, density
from sympy import Symbol, pprint

alpha = Symbol("alpha", positive = True)
beta = Symbol("beta", positive = True)
lambda = Symbol("lamda", nonnegative = True)
z = Symbol("z")

# Using sympy.stats.BetaNoncentral() method
X = BetaNoncentral("x", alpha, beta, lamda)
gfg = density(X)(z)

pprint(gfg, use_unicode = False)

oo _____ \ ` \ -lamda \ k ——- \ k + alpha – 1 /lamda\ beta – 1 2 ) z *|—–| *(1 – z) *e / \ 2 / / ———————————————— / B(k + alpha, beta)*k! /____, k = 0 

# Import sympy and betanoncentral
from sympy.stats import BetaNoncentral, density
from sympy import Symbol, pprint

alpha = 4
beta = 5
lambda = 1

# Using sympy.stats.BetaNoncentral() method
X = BetaNoncentral("x", alpha, beta, lamda)
gfg = density(X)(2)

pprint(gfg, use_unicode = False)

oo ____ \ ` \ -k k + 3 -1/2 \ 2 *2 *e / —————- / B(k + 4, 5)*k! /___, k = 0