import numpy as np
import matplotlib as mpl
import matplotlib.pyplot as plt
from scipy.special import gamma #导入伽马函数
from scipy.special import factorial
mpl.rcParams['axes.unicode_minus'] = False
mpl.rcParams['font.sans-serif'] = 'SimHei'
if __name__ == '__main__': #之后的语句作为模块被调用的时候,语句之后的代码不执行;直接使用的时候,语句之后的代码执行。通常,此语句用于模块测试中使用。
N = 5
x = np.linspace(0, N, 50)
y = gamma(x+1)
plt.figure(facecolor='w')
plt.plot(x, y, 'r-', x, y, 'mo', lw=2, ms=7) #lw表示线宽,ms表示点的大小
z = np.arange(0, N+1)
f = factorial(z, exact=True) # 求阶乘,当exact=True时,返回整数
print(f)
plt.plot(z, f, 'go', markersize=9)
plt.grid(b=True)
plt.xlim(-0.1,N+0.1)
plt.ylim(0.5, np.max(y)*1.05)
plt.xlabel('X', fontsize=15)
plt.ylabel('Gamma(X) - 阶乘', fontsize=15)
plt.title('阶乘和Gamma函数', fontsize=16)
plt.show()
伽马函数与阶乘关系的python实现
最新推荐文章于 2024-05-13 18:21:36 发布