Scipy积分

数值积分

一重积分(SciPy.integrate.quad):

例子为求解单位半圆的面积:

from scipy import integrate

def half_circle(x):
    return (1-x**2)**0.5

pi_half, err = integrate.quad(half_circle, -1, 1)

结果为:

In [4]: pi_half*2
Out[4]: 3.1415926535897967
In [5]: err
Out[5]: 1.0002354500215915e-09

https://docs.scipy.org/doc/scipy/reference/tutorial/integrate.html


  转载请注明: 石锅拌饭 Scipy积分

 上一篇
解决numpy 0d arrays error 解决numpy 0d arrays error
使用numpy array时,可能会遇到 “iteration over a 0-d array” 的错误。 例1:import numpy as np a = 1 a = np.asarray(a) for i in a: prin
2018-12-04
下一篇 
LeetCode 4 LeetCode 4
问题描述:给定两个排好序的数组,大小分别为$m$和$n$,找到两个数组的中位数(median),要求时间复杂度不超过$O(log(m+n))$. Example 1: nums1 = [1, 3] nums2 = [2] The medi
2018-12-03
  目录