import decimal import numpy as np from decimal import * getcontext().prec = 1000 def my_factorial (number): facto = Decimal ('1') for n in range (1, round (number) + 1, 1): facto = Decimal (facto) * Decimal (n) return facto bern = [] nmax = 2000 # n = 0 bern.append ((Decimal ('1'))) # n = 1 bern.append ((- Decimal ('1') / Decimal ('2'))) for n in range (2, nmax + 10, 1): temp = Decimal ('0') for i in range (0, round (n), 1): temp = Decimal (temp) - my_factorial (Decimal (n)) * Decimal (bern[i]) / (my_factorial (Decimal (n) + Decimal ('1') - Decimal (i)) * my_factorial (Decimal (i))) bern.append ((Decimal (temp))) for i in range (2, nmax + 4, 2): temp = Decimal (bern[i + 2]) / Decimal (bern[i]) if temp < 0: temp = Decimal (temp) * Decimal ('-1') temp = Decimal ('2') * Decimal (temp) ** Decimal ('0.5') / Decimal (i) print ('{:22.20f}'.format (temp))