# -*- 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)