数学建模习题2.2

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import quad

def fun(t, x):
return np.exp(-t) * (t ** (x - 1))

x = np.linspace(0, 10, 100) # x 的范围
y = [quad(fun, 0, np.inf, args=i)[0] for i in x] # 计算积分

plt.plot(x, y)
plt.xlabel('x')
plt.ylabel('$ y = \int_0^{\infty} e^{-t} \cdot t^{x-1} dt $')
plt.grid(True)
plt.show()

print("学号：2023310143005")