import numpy as np from scipy.special import gamma x = np.array([0.5, 1.0, 5.0, 0.0, -1.0])