From df2c01ef3704da364f6a5bad2b5f4dc02cffd84b Mon Sep 17 00:00:00 2001
From: Joshua Coles <josh@coles.to>
Date: Thu, 2 Mar 2023 00:08:58 +0000
Subject: [PATCH] Check stds for sp aren't all equal, compare to mva

---
 notebooks/stick-probability.ipynb | 46 +++++++++++++++++++++++++++++--
 1 file changed, 43 insertions(+), 3 deletions(-)

diff --git a/notebooks/stick-probability.ipynb b/notebooks/stick-probability.ipynb
index b7cf035..55a812e 100644
--- a/notebooks/stick-probability.ipynb
+++ b/notebooks/stick-probability.ipynb
@@ -4300,9 +4300,49 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
-   "outputs": [],
-   "source": [],
+   "execution_count": 74,
+   "outputs": [
+    {
+     "data": {
+      "text/plain": "(0.03445962437314559,\n 0.03483188535565984,\n 0.03297533970742226,\n 0.03436640348443668,\n 0.03337523830983404,\n 0.03260085335147516,\n 0.03906515517195682,\n 0.03254700858568894,\n 0.039006238984560744,\n 0.0345566212847466,\n 0.043316132964867936,\n 0.03544901612316204,\n 0.03283492691646963,\n 0.03451787755427769,\n 0.03159104220855551,\n 0.030857756004029485,\n 0.03462365543867886,\n 0.03132245597904412,\n 0.035128459151237926,\n 0.03079112194136338)"
+     },
+     "execution_count": 74,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "stds"
+   ],
+   "metadata": {
+    "collapsed": false
+   }
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 82,
+   "outputs": [
+    {
+     "data": {
+      "text/plain": "(1.0, 1.7343792431761986)"
+     },
+     "execution_count": 82,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "ps[7], means[7]"
+   ],
+   "metadata": {
+    "collapsed": false
+   }
+  },
+  {
+   "cell_type": "markdown",
+   "source": [
+    "Observing the result from mva $1.7311567309540636$ we obtain an answer to 3sf which isn't bad, although it may be worth seeing if we can push the maxParticles higher to achieve better outcomes"
+   ],
    "metadata": {
     "collapsed": false
    }