用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命令行参数解析模块getopt使用实例
Apr 13 Python
Python操作列表之List.insert()方法的使用
May 20 Python
Python函数any()和all()的用法及区别介绍
Sep 14 Python
Python使用numpy模块实现矩阵和列表的连接操作方法
Jun 26 Python
Django之创建引擎索引报错及解决详解
Jul 17 Python
Django Rest framework频率原理与限制
Jul 26 Python
PyQt5通信机制 信号与槽详解
Aug 07 Python
windows下python安装pip方法详解
Feb 10 Python
Python描述数据结构学习之哈夫曼树篇
Sep 07 Python
如何向scrapy中的spider传递参数的几种方法
Nov 18 Python
一文带你了解Python 四种常见基础爬虫方法介绍
Dec 04 Python
用Python自动清理系统垃圾的实现
Jan 18 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表单提交问题的解决方法
2011/04/12 PHP
浅析iis7.5安装配置php环境
2015/05/10 PHP
nicejforms——美化表单不用愁
2007/02/20 Javascript
jQuery插件实现表格隔行换色且感应鼠标高亮行变色
2013/09/22 Javascript
JS操作Cookie写入和读取实例代码
2013/10/20 Javascript
Jquery 获取指定标签的对象及属性的设置与移除
2014/05/29 Javascript
在AngularJS中使用jQuery的zTree插件的方法
2016/04/21 Javascript
JavaScript基于DOM操作实现简单的数学运算功能示例
2017/01/16 Javascript
jQuery Pagination分页插件_动力节点Java学院整理
2017/07/17 jQuery
JS常见DOM节点操作示例【创建 ,插入,删除,复制,查找】
2018/05/14 Javascript
教你使用vue-cli快速构建的小说阅读器
2019/05/13 Javascript
javascript实现简易聊天室
2019/07/12 Javascript
element 中 el-menu 组件的无限极循环思路代码详解
2020/04/26 Javascript
JS事件循环机制event loop宏任务微任务原理解析
2020/08/04 Javascript
[51:30]OG vs LGD 2018国际邀请赛淘汰赛BO3 第二场 8.26
2018/08/30 DOTA
[53:50]CHAOS vs Mineski 2019国际邀请赛小组赛 BO2 第一场 8.16
2019/08/18 DOTA
python使用xauth方式登录饭否网然后发消息
2014/04/11 Python
Python解释执行原理分析
2014/08/22 Python
python实现搜索指定目录下文件及文件内搜索指定关键词的方法
2015/06/28 Python
Python(TensorFlow框架)实现手写数字识别系统的方法
2018/05/29 Python
使用PyQtGraph绘制精美的股票行情K线图的示例代码
2019/03/14 Python
Python实现的删除重复文件或图片功能示例【去重】
2019/04/23 Python
TensorFlow keras卷积神经网络 添加L2正则化方式
2020/05/22 Python
python爬虫如何解决图片验证码
2021/02/14 Python
recorder.js 基于Html5录音功能的实现
2020/05/26 HTML / CSS
车间班组长岗位职责
2013/11/13 职场文书
电气工程和自动化自荐信范文
2013/12/25 职场文书
文字自荐书范文
2014/02/10 职场文书
国窖1573广告词
2014/03/21 职场文书
詹天佑教学反思
2014/04/30 职场文书
机关作风整顿个人整改措施2014
2014/09/17 职场文书
2014迎国庆演讲稿
2014/09/19 职场文书
单位个人查摆问题及整改措施
2014/10/28 职场文书
2015新年寄语(一句话)
2014/12/08 职场文书
2014年扶贫帮困工作总结
2014/12/09 职场文书
2019年二手房买卖合同范本
2019/10/14 职场文书