Как вычислить этот интеграл с помощью scipy/numpy? ⇐ Python

[Anonymous]

У меня есть следующий вопрос: введите описание изображения здесь
оно написано плохо, но для этих значений альфа и бета (0,32, 12) мне нужно найти минимальное k, при котором интеграл от -inf до inf это выражение меньше числа "пи".
но этот код, который я попробовал, не сработал, поскольку в примере я получил значение 4, а не 15, поэтому я не знаю, как продолжить.
вот что я получил с помощью ChatGPT: когда последний отпечаток показывает положительное значение, это означает, что я выполнил свое условие, и это происходит при k = 2. >
import numpy as np
from scipy.integrate import quad

def a(i):
return (12 * i - 11) ** 0.32

def sinc(t):
if t == 0:
return 1.0
else:
return np.sin(t) / t

def p(x, k):
product = 1.0
for i in range(1, k + 1):
product *= sinc(x / a(i))
return product

def P(k):
integral, _ = quad(lambda x: p(x, k), -np.inf, np.inf)
return integral

for value in range(1,5):
print(f"P({value}) =", P(value))
print(f'pi - P({value}) =', np.pi - P(value))

### solution was 2


Подробнее здесь: https://stackoverflow.com/questions/786 ... cipy-numpy