前言
这是一篇数值计算算法的开胃菜,主要介绍切比雪夫求积分法的使用;如果你遇到编程和数学的问题,强烈建议找到你身边的算法工程师或者计算机、电子、通信、数学专业的朋友,往往能够事半功倍。
求解Gamma函数
Gamma函数对特殊值
以下资料来自于维基百科
切比雪夫多项式展开
通过求解方程组的方式解出c[i]的值,即可通过c[I]求解出Gamma函数的值
#写于2019年3月13日
import numpy as np
Gam = np.array([2.363271801207,-3.544907701811,1.772453850906,1,\
0.886226925453,1,1.329340388179,2,3.323350970448,6])
a = [-1.5,-0.5,0.5,1,1.5,2,2.5,3,3.5,4]
n = len(a)
T = np.ones((n,n))
for i in range(1,n):
for j in range(n):
if 1 == i:
T[j][i] = a[j]