Update toy_notebook_en.ipynb

module2/exo1/toy_notebook_en.ipynb

@@ -22,4 +22,26 @@
  "nbformat": 4,
  "nbformat_minor": 2
 }
+from math import * print (pi)
+import numpy as np
+np.random.seed(seed=42)
+N = 10000
+x = np.random.uniform(size=N, low=0, high=1)
+theta = np.random.uniform(size=N, low=0, high=pi/2)
+2/(sum((x+np.sin(theta))>1)/N)
+
+ %matplotlib inline
+import matplotlib.pyplot as plt
+np.random.seed(seed=42)
+N = 1000
+x = np.random.uniform(size=N, low=0, high=1)
+y = np.random.uniform(size=N, low=0, high=1)
+1
+accept = (x*x+y*y) <= 1
+reject = np.logical_not(accept)
+fig, ax = plt.subplots(1)
+ax.scatter(x[accept], y[accept], c='b', alpha=0.2, edgecolor=None)
+ax.scatter(x[reject], y[reject], c='r', alpha=0.2, edgecolor=None)
+ax.set_aspect('equal')
+