用python求一重积分和二重积分的例子


Posted in Python onDecember 06, 2019

首先是对一元函数求积分,使用Scipy下的integrate函数:

from scipy import integrate
def g(x):
  return (1-x**2)**0.5
 
#用integrate下的quad函数可以同时求出积分结果和误差
res,err=integrate.quad(g,-1,1) #-1和1表示积分上下限,如果是正无穷用np.inf
print(res,err)

得到的结果如下,前者是积分结果,后者是误差

1.5707963267948986 1.0002356720661965e-09

如果是二重积分,比如是关于t和x的积分:

import numpy as np
from scipy.integrate import dblquad #dblquad用于二重积分
def main():
  print(dblquad(lambda t,x:np.sin(t)*np.exp(-x*t)/t**5,0.5,0.8,lambda x:0.2,lambda x:0.7))
 
if __name__ == "__main__":
  main()

在上面的括号中,被积函数是sint*exp(-xt)/t^5, 其中t的积分上下限是0.5和0.8,x的积分上下限是0.2和0.7,结果如下:

(10.041536868759145, 1.3645801372419182e-11)

以上这篇用python求一重积分和二重积分的例子就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python发布模块的步骤分享
Feb 21 Python
使用Python脚本将文字转换为图片的实例分享
Aug 29 Python
对Python3中的print函数以及与python2的对比分析
May 02 Python
带你认识Django
Jan 15 Python
浅谈pandas筛选出表中满足另一个表所有条件的数据方法
Feb 08 Python
详解python 利用echarts画地图(热力图)(世界地图,省市地图,区县地图)
Aug 06 Python
python rsync服务器之间文件夹同步脚本
Aug 29 Python
Python基础之函数原理与应用实例详解
Jan 03 Python
如何验证python安装成功
Jul 06 Python
python字典通过值反查键的实现(简洁写法)
Sep 30 Python
Python中常见的反爬机制及其破解方法总结
Jun 10 Python
Python基本数据类型之字符串str
Jul 21 Python
解决Numpy中sum函数求和结果维度的问题
Dec 06 #Python
numpy按列连接两个维数不同的数组方式
Dec 06 #Python
使用Python实现分别输出每个数组
Dec 06 #Python
Python 获取numpy.array索引值的实例
Dec 06 #Python
python的json中方法及jsonpath模块用法分析
Dec 06 #Python
python输出数组中指定元素的所有索引示例
Dec 06 #Python
关于numpy.where()函数 返回值的解释
Dec 06 #Python
You might like
PHP中几种常见的超时处理全面总结
2012/09/11 PHP
PHP实现的简单日历类
2014/11/29 PHP
学习php设计模式 php实现原型模式(prototype)
2015/12/07 PHP
使用PHPUnit进行单元测试并生成代码覆盖率报告的方法
2019/03/08 PHP
使用jQuery UI的tooltip函数修饰title属性的气泡悬浮框
2013/06/24 Javascript
JS实现多物体缓冲运动实例代码
2013/11/29 Javascript
jquerydom对象的事件隐藏显示和对象数组示例
2013/12/10 Javascript
JavaScript实现两个Table固定表头根据页面大小自行调整
2014/01/03 Javascript
jquery操作select大全
2014/04/25 Javascript
从零学jquery之如何使用回调函数
2014/05/16 Javascript
nodejs中的fiber(纤程)库详解
2015/03/24 NodeJs
JavaScript中实现键值对应的字典与哈希表结构的示例
2016/06/12 Javascript
BootStrap智能表单实战系列(四)表单布局介绍
2016/06/13 Javascript
JS 在数组指定位置插入/删除数据的方法
2017/01/12 Javascript
bootstrap按钮插件(Button)使用方法解析
2017/01/13 Javascript
使用 NodeJS+Express 开发服务端的简单介绍
2017/04/07 NodeJs
JS中的Replace()传入函数时的用法详解
2017/09/11 Javascript
使用Bootstrap和Vue实现用户信息的编辑删除功能
2017/10/25 Javascript
JQuery实现ajax请求的示例和注意事项
2018/12/10 jQuery
vue props 单项数据流实例分享
2020/02/16 Javascript
[46:37]LGD vs TNC 2019国际邀请赛小组赛 BO2 第二场 8.15
2019/08/16 DOTA
在Python中操作时间之tzset()方法的使用教程
2015/05/22 Python
利用Python暴力破解zip文件口令的方法详解
2017/12/21 Python
Python3导入自定义模块的三种方法详解
2018/04/13 Python
Python实现获取nginx服务器ip及流量统计信息功能示例
2018/05/18 Python
Tensorflow 实现修改张量特定元素的值方法
2018/07/30 Python
Python 基于FIR实现Hilbert滤波器求信号包络详解
2020/02/26 Python
python switch 实现多分支选择功能
2020/12/21 Python
css3 transform过渡抖动问题解决
2020/10/23 HTML / CSS
英国男女豪华配饰和礼品网站:Black.co.uk
2020/02/28 全球购物
迪士尼法国在线商店:shopDisney FR
2020/12/03 全球购物
用人单位终止解除劳动合同证明书
2014/10/06 职场文书
2015年学校保卫部工作总结
2015/05/11 职场文书
入党积极分子党支部意见
2015/06/02 职场文书
《曾国藩家书》读后感——读家书,立家风
2019/08/21 职场文书
Python作用域和名称空间的详细介绍
2022/04/13 Python