psi 函数计算

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 展开式

• C++ 版本， http://people.sc.fsu.edu/~jburkardt/cpp_src/asa103/asa103.cpp

向量化

import numpy as np
psi = np.frompyfunc(r8_psi,1,1)