From 55ac7649cfddf63cf41d16e6d9a2211d363706f6 Mon Sep 17 00:00:00 2001 From: d932cd8408420be151b222f5aab7b293 Date: Thu, 15 Dec 2022 21:50:40 +0000 Subject: [PATCH] Update toy_notebook_en.ipynb --- module2/exo1/toy_notebook_en.ipynb | 22 ++++++++++++++++++++++ 1 file changed, 22 insertions(+) diff --git a/module2/exo1/toy_notebook_en.ipynb b/module2/exo1/toy_notebook_en.ipynb index 0bbbe37..fac4d46 100644 --- a/module2/exo1/toy_notebook_en.ipynb +++ b/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') + -- 2.18.1