j-jith / orthopoly Goto Github PK

Example:

def pdf(x):
    return math.exp(-x)

pp = OrthoPoly(pdf, intlims=(0,1))
pp.gen_poly(n)

# a new class field 'betas' has been added to store the beta recurrence coefficients for demonstration
print(pp.betas)

Output:

[1.0,
0.07932640579220766,
0.06711669572587618,
0.06433876970545459,
0.06350596398737984,
0.06313659912953457,
0.06293951687236084,
0.06282180581953442,
0.0627458612867561,
0.06269478868364445,
0.06283407650853606,
0.06089373421606997,
0.21595648578197688,
0.669882586586399,
-1.3444323914957315,
4.962392290085074,
0.4906144952293512,
0.20512208494961423,
-686.035257250946,
15925.79931665891,
56662.12826844703]

The last elements demonstrate instability.