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