text: # Take 10,000 bootstrap replicates of the mean: bs_replicates
  bs_replicates = draw_bs_reps(rainfall, np.mean, size=10000)
  
  # Compute and print SEM
  sem = np.std(rainfall) / np.sqrt(len(rainfall))
  print(sem)
  
  # Compute and print standard deviation of bootstrap replicates
  bs_std = np.std(bs_replicates)
  print(bs_std)
  
  # Make a histogram of the results
  _ = plt.hist(bs_replicates, bins=50, normed=True)
  _ = plt.xlabel('mean annual rainfall (mm)')
  _ = plt.ylabel('PDF')
  
  # Show the plot
  plt.show()
  