test

parent 272867e7
...@@ -70,7 +70,7 @@ ...@@ -70,7 +70,7 @@
"## 1.3 Using a surface fraction argument\n", "## 1.3 Using a surface fraction argument\n",
"\n", "\n",
"A method that is easier to understand and does not make use of the sin function is based on the\n", "A method that is easier to understand and does not make use of the sin function is based on the\n",
"fact that if *X ∼ U(0, 1)* and *Y ∼ U(0, 1)* then P[X^2 + Y^2 ≤ 1] = &#960;/4 (see [\"Monte Carlo method\"\n", "fact that if *X ∼ U(0, 1)* and *Y ∼ U(0, 1)* then P[X<sup>2</sup> + Y<sup>2</sup> ≤ 1] = &#960;/4 (see [\"Monte Carlo method\"\n",
"on Wikipedia](https://en.wikipedia.org/wiki/Monte_Carlo_method)). The following code uses this approach:" "on Wikipedia](https://en.wikipedia.org/wiki/Monte_Carlo_method)). The following code uses this approach:"
] ]
}, },
...@@ -108,6 +108,16 @@ ...@@ -108,6 +108,16 @@
"ax.set_aspect('equal')\n" "ax.set_aspect('equal')\n"
] ]
}, },
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"It is then straightforward to obtain a (not really good) approximation to π by counting how\n",
"many times, on average, X<sup>2</sup> + Y<sup>2</sup> is smaller than 1:"
]
},
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 4, "execution_count": 4,
...@@ -125,25 +135,7 @@ ...@@ -125,25 +135,7 @@
} }
], ],
"source": [ "source": [
"4*np.mean(accept)\n" "4*np.mean(accept)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"ename": "SyntaxError",
"evalue": "invalid syntax (<ipython-input-6-2b522dbd7382>, line 1)",
"output_type": "error",
"traceback": [
"\u001b[0;36m File \u001b[0;32m\"<ipython-input-6-2b522dbd7382>\"\u001b[0;36m, line \u001b[0;32m1\u001b[0m\n\u001b[0;31m 4*np.mean(accept.5f)\u001b[0m\n\u001b[0m ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
]
}
],
"source": [
"\n"
] ]
}, },
{ {
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment