python中使用scipy.integrate求积分(二重、三重)
#1.1
"""
from sympy import *
#integrate(函数,(变量,下限, 上限))
x = symbols('x')
print(integrate(2*x**2+1, (x, 0, 1)))
"""
from scipy.integrate import quad
def integrand(x, a, b):
return a*x**2 + b
a = 2
b = 1
I = quad(integrand, 0, 1, args=(a,b))
print("1.1的I=",I)
print("*********************************************************")
#1.2
from scipy.integrate import quad
from numpy import sin, cos, pi
def integrand(w):
return sin(w)*cos(b)**3
I = quad(integrand, 0, pi/2)
print("1.2的I=",I)
print("**********************************************************")
#1.3
import numpy as np
from scipy.integrate import quad
def integrand(t):
return t*np.exp(-(t*t)/2)
I = quad(integra