切比雪夫求积分法（附全过程代码）

前言

  这是一篇数值计算算法的开胃菜，主要介绍切比雪夫求积分法的使用；如果你遇到编程和数学的问题，强烈建议找到你身边的算法工程师或者计算机、电子、通信、数学专业的朋友，往往能够事半功倍。

求解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]