numpy积分计算结果错误

Question

我正在尝试计算 sin(x)*sin(2x) 在 - pi,pi 范围内的积分，但我找不到 0。我尝试使用 math.randiand(math.pi) 但仍然错误。

import scipy as sp
import numpy as np

def integrand(x):
  return math.sin(x)*math.sin(2*x)

def int_test(a,b):
  return sp.integrate.quad(integrand,-math.pi,math.pi)

Answer 1

我不得不尝试如何导入集成库：

import scipy.integrate as integrate
import numpy as np
import math

def integrand(x):
  return math.sin(x)*math.sin(2*x)

def int_test(a, b):
  return integrate.quad(integrand, a, b)

x = int_test(-math.pi, math.pi)
print(x[0])

→ python3 test.py 
1.5156111913114e-16

这非常接近于 0。

请注意，您会得到一个元组作为响应。 您可以在这里阅读返回值是什么： https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.quad.ZFC35FDC70D5FC69D2693883A822EZ5A