使用Python实现定积分、二重定积分(简单实用)
1.概述
最近项目需要使用程序实现数学微积分,最初想用java实现,后来发现可用文档太少,实现比较麻烦,后来尝试使用python实现,代码量较少,主要有sympy与scipy两种实现方式,本文主要记录scipy的实现方式。
2.内容
2.1 所求函数
2.2 python代码
# 引入需要的包
import scipy.integrate
from numpy import exp
from math import sqrt
import math
# 创建表达式
f = lambda x,y : exp(x**2-y**2)
# 计算二重积分:(p:积分值,err:误差)
# 这里注意积分区间的顺序
# 第二重积分的区间参数要以函数的形式传入
p,err= scipy.integrate.dblquad(f, 0, 2, lambda g : 0, lambda h : 1)
print(p)
2.3 注意问题
1. exp尽量使用numpy的exp
2. 注意积分区间参数的顺序
3. 第二重积分的区间参数要以函数的形式传入