import numpy as np

V13 = -0.2886751346
V15 = -0.4082482906
V16 = 0.06763826297
V17 = 0.8633800241
V36 = -0.02153282025
V37 = 0.9571849376
V46 = 0.1961805388
V47 = -0.9805677930


count = 0
trials = 30000000
for i in range (1, trials):
  Z1 = np.random.normal(0, 1, 1)
  Z2 = np.random.normal(0, 1, 1)
  Z3 = np.random.normal(0, 1, 1)
  Z4 = np.random.normal(0, 1, 1)
  Z5 = np.random.normal(0, 1, 1)
  Z6 = np.random.normal(0, 1, 1)
  Z7 = np.random.normal(0, 1, 1)
  if((Z7 > (-V13*Z3-V15*Z5-V16*Z6)/V17) and (Z7 > (V13*Z3+V15*Z5-V16*Z6)/V17) and (Z7 > (-V13*Z5-V36*Z6)/V37) and (Z7 < (V46*Z6)/(-V47)) and (Z7 > (V13*Z5-V36*Z6)/V37) and (Z7 > (-V13*Z4+V15*Z5-V16*Z6)/V17) and (Z7 > (V13*Z4+V15*Z5-V16*Z6)/V17)): 
    count = count + 1
print(count)
print(count/trials)