scipy 特殊函数

概率统计的伽马函数

import scipy.special as Sprint S.gamma(4)6.0print S.gamma(0.5)1.77245385091print S.gamma(1+1j) #复数也可以(0.498015668118-0.154949828302j)print S.gamma(1000)inf

数太大使用对数伽马函数

S.gammaln(1000)Out[6]: 5905.2204232091808

ｌｏｇ（１＋ｘ）函数ｌｏｇ１ｐ（），精度较浮点数高

print 1+1e-201.0print np.log(1+1e-20)0.0print S.log1p(1e-20)1e-20

ellipj(u,m)计算雅可比椭圆函数

可以传入数组，支持广播运算

result=S.sllipj(u[:,None],m[None:])