import numpy as np
   
n = 10000  # number of trials e.g. 10, 100, 1000, 10000 

# H=head, T=tail
sample_space = ['H', 'T'] 
event = 0 #Head

for k in range (0, n):
    sample_k = np.random.choice(sample_space)
    
    if ( sample_k == 'H'): 
        event = event + 1
              
print('Head occured:', f'{event/n :.4f}')