用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的web.py框架实现类似Django的ORM查询的教程
May 02 Python
Python中的日期时间处理详解
Nov 17 Python
Python实现曲线拟合操作示例【基于numpy,scipy,matplotlib库】
Jul 12 Python
Python实现按逗号分隔列表的方法
Oct 23 Python
Python3.5基础之变量、数据结构、条件和循环语句、break与continue语句实例详解
Apr 26 Python
六行python代码的爱心曲线详解
May 17 Python
Django框架封装外部函数示例
May 28 Python
Python flask框架post接口调用示例
Jul 03 Python
Python代码注释规范代码实例解析
Aug 14 Python
Python基于mediainfo批量重命名图片文件
Dec 29 Python
python批量更改目录名/文件名的方法
Apr 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/01/12 PHP
php打包压缩文件之ZipArchive方法用法分析
2016/04/30 PHP
PHP在innodb引擎下快速代建全文搜索功能简明教程【基于xunsearch】
2016/10/14 PHP
PHP实现上传多文件示例代码
2017/02/20 PHP
thinkphp3.2实现在线留言提交验证码功能
2017/07/19 PHP
PHP回调函数与匿名函数实例详解
2017/08/16 PHP
php实现对文件压缩简单的方法
2019/09/29 PHP
Aster vs KG BO3 第一场2.18
2021/03/10 DOTA
从阿里妈妈发现的几个不错的表单验证函数
2007/09/21 Javascript
JQuery 学习技巧总结
2010/05/21 Javascript
基于datagrid框架的查询
2013/04/08 Javascript
jquery实现弹出窗口效果的实例代码
2013/11/28 Javascript
JS弹出窗口插件zDialog简单用法示例
2016/06/12 Javascript
JavaScript编写一个简易购物车功能
2016/09/17 Javascript
canvas 弹幕效果(实例分享)
2017/01/11 Javascript
nodejs socket实现的服务端和客户端功能示例
2017/06/02 NodeJs
通过命令行创建vue项目的方法
2017/07/20 Javascript
nuxt+axios解决前后端分离SSR的示例代码
2017/10/24 Javascript
在Vue中使用echarts的方法
2018/02/05 Javascript
React中如何引入Angular组件详解
2018/08/09 Javascript
js神秘的电报密码 哈弗曼编码实现
2019/09/10 Javascript
vue+iview分页组件的封装
2020/11/17 Vue.js
深入解析Python设计模式编程中建造者模式的使用
2016/03/02 Python
Scrapy的简单使用教程
2017/10/24 Python
python交换两个变量的值方法
2019/01/12 Python
pyinstaller打包opencv和numpy程序运行错误解决
2019/08/16 Python
使用卷积神经网络(CNN)做人脸识别的示例代码
2020/03/27 Python
python openCV自制绘画板
2020/10/27 Python
python中用Scrapy实现定时爬虫的实例讲解
2021/01/18 Python
html5 canvas 简单画板实现代码
2012/01/05 HTML / CSS
Desigual德国官网:在线购买原创服装
2018/03/27 全球购物
英国领先的野生鸟类食品供应商:GardenBird
2018/08/09 全球购物
TripAdvisor日本:全球领先的旅游网站
2019/02/14 全球购物
应届生求职推荐信
2013/10/28 职场文书
黄石寨导游词
2015/02/05 职场文书
粗暴解决CUDA out of memory的问题
2021/05/22 Python