Python的剧情的积分与离散和连续参数

Question

我有以下功能：

def sbeta(beta,gamma,y): 
    k = 2/(np.pi)**2. 
    return k * np.sqrt(1 - (np.sqrt(1-y**2.)*np.sin(beta)*np.cos(gamma) - y*np.cos(beta))**2.) 

其中的β为常数，y是在-1和1之间定义为：

beta = 23.4 
y = numpy.linspace(-1, 1, 100) 

欲绘制该函数对于从0到2pi评估的伽玛的积分：

def integral(beta,gamma,y): 
    for i in range(len(y)): 
     I = integrate.quad(sbeta, 0., 2*np.pi, args=(beta, y[i])) 
     print(I) 
     plt.plot(y[i],I[0]) 

gamma = np.linspace(0., 2*np.pi, 10) 
integral(beta,gamma,y) 
plt.show() 

此时没有错误，但我不认为这是正确的。我想计算gamma作为从0到2pi的变量的积分，但y是一个离散数组。如何从-1和1计算并绘制10 y值？我应该使用scipy.integrate.cumtrapz吗？

Answer 1

这是我尝试了答案：

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


def f(gamma, y): 
    k = 2/(np.pi)**2. 
    beta = 23.4*(np.pi/180.) 
    return k * np.sqrt(1 - (np.sqrt(1-y**2.)*np.sin(beta)*np.cos(gamma) - y*np.cos(beta))**2.) 


y = np.linspace(-1., 1.,10) 
low = 0. 
high = 2*np.pi 

Ivals = [] 
for i in range(len(y)): 
    I = integrate.quad(f, low, high, args=(y[i])) 
    plt.scatter(y[i],I[0]) 
    Ivals.append(I) 
plt.show() 

这显示了以下情节：

谁能确认这是否是正确的吗？