Check stds for sp aren't all equal, compare to mva

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
    }