python使用integrate函数计算球体体积,具体方法如下:
from scipy import integrate
import numpy as np
import time
start=time.clock()
def half_circle(x):
return (1-x**2)**0.5
def half_sphere(x,y):
return (1-x**2-y**2)**0.5
volume,error=integrate.dblquad(half_sphere,-1,1,
lambda x: -half_circle(x),
lambda x: half_circle(x))
end=time.clock()
print('{}\n{}\n{}'.format(volume,error,np.pi*4/3/2))
print('time={} Sconds'.format(end-start))