利用python求积分

python的numpy库集成了很多的函数。利用其中的函数可以很方便的解决一些数学问题。本篇介绍如何使用python的numpy来求解积分。代码如下：

# -*- coding: utf-8 -*-
import numpy as np
from scipy.integrate import quad,dblquad,nquad

def main():

    print quad(lambda  x:np.exp(-x),0,np.inf)
    '''求积分，np.inf代表正无穷。
    结果第一个数值代表运算结果，第二个数值代表误差
    '''
    print dblquad(lambda t,x:np.exp(-x*t)/t**3,0,np.inf,lambda x:1,lambda x:np.inf)
    '''
    求二重积分 然后给t,x赋积分区间
    lambda是匿名函数
    '''


if __name__ == "__main__":
    main()

结果如下：

(1.0000000000000002, 5.842607038578007e-11)
(0.3333333333366853, 1.3888461883425516e-08)