scipy.special.psi
odps中不支持 scipy.special.psi，需要改写
基于 chebyshev_polynomial
- https://people.sc.fsu.edu/~jburkardt/m_src/chebyshev_polynomial/r8_psi.m
- https://people.sc.fsu.edu/~jburkardt/py_src/polpak/r8_psi.py
- java 函数 http://jdistlib.sourceforge.net/ 可以下到源代码
基于Stirling 展开式
向量化
import numpy as np
psi = np.frompyfunc(r8_psi,1,1)