def lpol_fima(d, n=20):
 from scipy.special import gammaln
 j = np.arange(n)
 return np.exp(((gammaln((d + j)) - gammaln((j + 1))) - gammaln(d)))

def lpol_fima ( d , n = 20 ) : from scipy . special import gammaln j = np . arange ( n ) return np . exp ( ( ( gammaln ( ( d + j ) ) - gammaln ( ( j + 1 ) ) ) - gammaln ( d ) ) )
ma representation of fractional integration .