Update pi_monte_carlo_estimation.py

Author: cclauss

maths/pi_monte_carlo_estimation.py

@@ -26,7 +26,7 @@ def random_unit_square(cls, ran: random.Random):
 
 def estimate_pi(number_of_simulations: int, seed: int | None = None) -> float:
     """
-    Generates an estimate of the mathematical constant PI.
+    Generate an estimate of the mathematical constant PI.
     See https://en.wikipedia.org/wiki/Monte_Carlo_method#Overview
 
     The estimate is generated by Monte Carlo simulations. Let U be uniformly drawn from
@@ -46,7 +46,7 @@ def estimate_pi(number_of_simulations: int, seed: int | None = None) -> float:
             circle, which is called m.
         3. An estimate of P[U in unit circle] is m/n
 
-    'seed' provides seed for the number generator - if None, no seed is used.
+    'seed' provides a seed for the number generator; if None, no seed is used.
 
     >>> estimate_pi(100, 1)
     3.2
@@ -74,3 +74,6 @@ def estimate_pi(number_of_simulations: int, seed: int | None = None) -> float:
     import doctest
 
     doctest.testmod()
+    prompt = "Please enter the desired number of Monte Carlo simulations: "
+    my_pi = estimate_pi(int(input(prompt).strip()))
+    print(f"An estimate of PI is {my_pi} with an error of {abs(my_pi - pi)}")