第8章-随机数应用-蒙特卡洛方法-课堂练习I

第1关：计算不规则区域 (Viviani 体) 体积

任务描述
本关任务：编写一个能通过使用蒙特卡洛方法计算不规则区域 (Viviani 体) 体积的小程序。

相关知识
为了完成本关任务，你需要掌握：1.什么是不规则区域 (Viviani 体) 体积，2.蒙特卡洛方法如何使用。

编程要求
根据提示，在右侧编辑器补充代码，计算并输出不规则区域 (Viviani 体) 的体积。

from numpy import random
from math import *

B = 2; N = 100000; M = 0

X = random.uniform(-B, B, size=N)
Y = random.uniform(-B, B, size=N)
Z = random.uniform(-B, B, size=N)

# please write your code below from begin to end.
# -------------------- begin --------------------
for x,y,z in zip(X,Y,Z):
    if x**2+y**2+z**2<=4 and (x-1)**2+y**2<=1:
        M+=1
volume=M/N*64



# -------------------- begin --------------------

 
print('%.f'%volume)