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)