import yaml
import os
import mpmath

path = 'data/P-adic_numbers/P-adic_Gamma_function_at_integers/'

prec10 = 50 #relative precision in base 10

p_range = prime_range(30)
k_range = [-40..40] 

numbers = {}
for p in p_range:
	print("p:",p)
	numbers_p = {}

	prec_p = ceil(prec10*log(10,p))
	Q_p = Qp(p, prec=prec_p, print_mode='val-unit')

	for k in k_range:
		number = Q_p(k).gamma()
	
		number_str = str(number)
		numbers_p[str(k)] = number_str

	numbers[str(p)] = numbers_p

filename = os.path.join(path, 'numbers.yaml')
yaml.dump(numbers, stream = open(filename, 'w'), sort_keys = False)