def lpol_fiar(d, n=20):
 from scipy.special import gammaln
 j = np.arange(n)
 ar = (- np.exp(((gammaln(((- d) + j)) - gammaln((j + 1))) - gammaln((- d)))))
 ar[0] = 1
 return ar

def lpol_fiar ( d , n = 20 ) : from scipy . special import gammaln j = np . arange ( n ) ar = ( - np . exp ( ( ( gammaln ( ( ( - d ) + j ) ) - gammaln ( ( j + 1 ) ) ) - gammaln ( ( - d ) ) ) ) ) ar [ 0 ] = 1 return ar
ar representation of fractional integration .