香港云服务器python怎么用函数算体积

已关闭留言

香港云服务器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))