用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模拟登录百度代码分享(获取百度贴吧等级)
Dec 27 Python
跟老齐学Python之眼花缭乱的运算符
Sep 14 Python
Python使用xlwt模块操作Excel的方法详解
Mar 27 Python
解决Python3 被PHP程序调用执行返回乱码的问题
Feb 16 Python
Python基本数据结构之字典类型dict用法分析
Jun 08 Python
python替换字符串中的子串图文步骤
Jun 19 Python
对django layer弹窗组件的使用详解
Aug 31 Python
python破解bilibili滑动验证码登录功能
Sep 11 Python
windows10环境下用anaconda和VScode配置的图文教程
Mar 30 Python
Python尾递归优化实现代码及原理详解
Oct 09 Python
python 实现aes256加密
Nov 27 Python
Python tkinter实现日期选择器
Feb 22 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代码把全角数字转为半角数字
2007/12/10 PHP
解决ThinkPHP关闭调试模式时报错的问题汇总
2015/04/22 PHP
PHP hex2bin()函数用法讲解
2019/02/25 PHP
Laravel jwt 多表(多用户端)验证隔离的实现
2019/12/18 PHP
页面中body onload 和 window.onload 冲突的问题的解决
2009/07/01 Javascript
JS遮罩层效果 兼容ie firefox jQuery遮罩层
2010/07/26 Javascript
使用jQuery轻松实现Ajax的实例代码
2010/08/16 Javascript
jQuery.ajax 用户登录验证代码
2010/10/29 Javascript
javascript full screen 全屏显示页面元素的方法
2013/09/27 Javascript
javascript获取鼠标点击元素对象(示例代码)
2013/12/20 Javascript
jQuery中:last选择器用法实例
2014/12/30 Javascript
一看就懂:jsonp详解
2015/06/01 Javascript
分享网页检测摇一摇实例代码
2016/01/14 Javascript
JavaScript的MVVM库Vue.js入门学习笔记
2016/05/03 Javascript
深入分析javascript中的错误处理机制
2016/07/17 Javascript
javascript获取指定区间范围随机数的方法
2017/09/08 Javascript
JS中的算法与数据结构之栈(Stack)实例详解
2019/08/20 Javascript
js获取本日、本周、本月的时间代码
2020/02/01 Javascript
ant-design-vue 时间选择器赋值默认时间的操作
2020/10/27 Javascript
python 获取et和excel的版本号
2009/04/09 Python
Python命令行参数解析模块optparse使用实例
2015/04/13 Python
python迭代器与生成器详解
2016/03/10 Python
关于python的list相关知识(推荐)
2017/08/30 Python
Python numpy生成矩阵、串联矩阵代码分享
2017/12/04 Python
使用Python向DataFrame中指定位置添加一列或多列的方法
2019/01/29 Python
Python学习笔记之变量、自定义函数用法示例
2019/05/28 Python
Python range与enumerate函数区别解析
2020/02/28 Python
在Python IDLE 下调用anaconda中的库教程
2020/03/09 Python
python读取hdfs上的parquet文件方式
2020/06/06 Python
Python‘==‘ 及 ‘is‘相关原理解析
2020/09/05 Python
Nike英国官网:Nike.com (UK)
2017/02/13 全球购物
英国体育器材进口商店:UK Sport Imports
2017/03/14 全球购物
英语演讲稿3分钟
2014/04/29 职场文书
基层党建工作汇报材料
2014/08/15 职场文书
党员剖析材料范文
2014/09/30 职场文书
漫画《催眠麦克风-Dawn Of Divisions》第二卷PV公开
2022/04/05 日漫