用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中的引用和拷贝浅析
Nov 22 Python
python实现将html表格转换成CSV文件的方法
Jun 28 Python
利用 python 对目录下的文件进行过滤删除
Dec 27 Python
mvc框架打造笔记之wsgi协议的优缺点以及接口实现
Aug 01 Python
Python基本socket通信控制操作示例
Jan 30 Python
使用python实现ftp的文件读写方法
Jul 02 Python
Python通过cv2读取多个USB摄像头
Aug 28 Python
Python实现Keras搭建神经网络训练分类模型教程
Jun 12 Python
什么是python的函数体
Jun 19 Python
python2和python3哪个使用率高
Jun 23 Python
pycharm 如何取消连按两下shift出现的全局搜索
Jan 15 Python
python中使用asyncio实现异步IO实例分析
Feb 26 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后门代码解析
2014/07/05 PHP
微信利用PHP创建自定义菜单的方法
2016/08/01 PHP
php curl上传、下载、https登陆实现代码
2017/07/23 PHP
快速解决PHP调用Word组件DCOM权限的问题
2017/12/27 PHP
thinkPHP3.2.3结合Laypage实现的分页功能示例
2018/05/28 PHP
centos7上编译安装php7以php-fpm方式连接apache
2018/11/08 PHP
php使用gearman进行任务分发操作实例详解
2020/02/26 PHP
JavaScript 对象的属性和方法4种不同的类型
2010/03/19 Javascript
JQuery插件Style定制化方法的分析与比较
2012/05/03 Javascript
jQuery后代选择器用法实例
2014/12/23 Javascript
javascript跨域原因以及解决方案分享
2015/04/08 Javascript
JavaScript知识点总结之如何提高性能
2016/01/15 Javascript
Bootstrap模态框调用功能实现方法
2016/09/19 Javascript
vue-quill-editor富文本编辑器简单使用方法
2018/09/21 Javascript
详解Angular模板引用变量及其作用域
2018/11/23 Javascript
JQuery样式与属性设置方法分析
2019/12/07 jQuery
详解Vue数据驱动原理
2020/11/17 Javascript
[07:08]2014DOTA2西雅图国际邀请赛 小组赛7月11日TOPPLAY
2014/07/11 DOTA
[04:45]上海特级锦标赛主赛事第三日TOP10
2016/03/05 DOTA
python传递参数方式小结
2015/04/17 Python
python中引用与复制用法实例分析
2015/06/04 Python
python实现mysql的读写分离及负载均衡
2018/02/04 Python
Python 字符串操作(string替换、删除、截取、复制、连接、比较、查找、包含、大小写转换、分割等)
2018/03/19 Python
解决python有时候import不了当前的包问题
2019/08/28 Python
pytorch标签转onehot形式实例
2020/01/02 Python
python通过opencv实现图片裁剪原理解析
2020/01/19 Python
python查找特定名称文件并按序号、文件名分行打印输出的方法
2020/04/24 Python
CSS3中的clip-path使用攻略
2015/08/03 HTML / CSS
HTML5样式控制示例代码
2013/11/27 HTML / CSS
伦敦高级内衣品牌:Agent Provocateur(大内密探)
2016/08/23 全球购物
Martinelli官方商店:西班牙皮鞋和高跟鞋品牌
2019/07/30 全球购物
黄河的主人教学反思
2014/02/07 职场文书
大学四年个人的自我评价
2014/02/26 职场文书
2019年新郎保证书3篇
2019/10/17 职场文书
TensorFlow中tf.batch_matmul()的用法
2021/06/02 Python
ipad隐藏软件app图标方法
2022/04/19 数码科技