scipy5.01：积分应用：椭球体积的数值积分计算

积分应用：椭球体积的数值积分计算

计算方法：

代码：

"""
scipy1.1：椭球体积的数值积分计算
马玉华   2019.11.6
"""
import numpy as np
import scipy
from scipy import integrate

def main():
    a = float(input('请输入椭球参数a的值：'))
    b = float(input('请输入椭球参数b的值：'))
    c = float(input('请输入椭球参数c的值：'))
    def f(x,y):
        return 2*c*np.sqrt(1-(x**2)/(a**2)-(y**2)/(b**2))
    def bound_y(x):
        bound = b*np.sqrt(1-(x**2)/(a**2))
        return [-bound,bound]
    def bound_x():
        return [-a,a]
    result = integrate.nquad(f,[bound_y,bound_x])
    print('椭球的体积为{}'.format(result[0]),',误差值为{}'.format(result[1]))


i