# This script file is a solution to exercise 3 in Eriksen og Lauritzen (2004).
# Answers are held in  a, b og c

import scipy
from scipy.stats import norm
x = norm.rvs(30, .05, (1000, 4))
y = x.min(axis = 1)
z = x.max(axis = 1)
a = sum(y > 30)
b = sum(z < 29.9)
c = 1000 - sum(z < 30.15)
print(a, b, c)

# a=number of samples where the minimum is greater than 30
# b=number of samples where the maximum is less than 29.9
# c=number of samples where the maximum is less than 30.15