From 20ebd93c6a1c86c080d198fd6fe925b9a520bb6f Mon Sep 17 00:00:00 2001
From: Yifan Wang <yifan.wang@aei.mpg.de>
Date: Tue, 17 Jan 2023 13:55:58 +0100
Subject: [PATCH] add two pulse
---
notebook.ipynb | 124 +++++++++++++++++++++++++++++++++++++++++++++++--
1 file changed, 120 insertions(+), 4 deletions(-)
diff --git a/notebook.ipynb b/notebook.ipynb
index 0ba2b67..a9022d0 100644
--- a/notebook.ipynb
+++ b/notebook.ipynb
@@ -1003,7 +1003,7 @@
},
{
"cell_type": "code",
- "execution_count": 82,
+ "execution_count": 17,
"metadata": {},
"outputs": [],
"source": [
@@ -1029,18 +1029,134 @@
" sigma is the unknown width of the pulse\n",
" '''\n",
" s=signal_twopulse(A,sigma,x_1,x_2)\n",
- " template = signal(A,sigma,x_1+offset1,x_2+offset2)/A\n",
+ " template = signal_twopulse(A,sigma,x_1+offset1,x_2+offset2)/A\n",
" \n",
" stat_h0 = []\n",
" stat_h1 = []\n",
" for i in range(1000):\n",
" n = np.random.normal(0,1,100)\n",
" d = n+s\n",
- " stat_h0.append(t_unknownA(n,template))\n",
- " stat_h1.append(t_unknownA(d,template))\n",
+ " stat_h0.append(t_unknownAsigma(n,template))\n",
+ " stat_h1.append(t_unknownAsigma(d,template))\n",
" return stat_h0,stat_h1"
]
},
+ {
+ "cell_type": "code",
+ "execution_count": 20,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnAAAAGbCAYAAABEeK1/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAa+UlEQVR4nO3dTW9c14Em4PfEHzI6o4TieGAEDXhsGvA+lHovA+Q/kKL5AybXmYUIzR8QqE3WpH+BLe961RCDcN8Wt8EEcKUnQDcGA5siLKVhybDPLOpSKtFVVBVFVvGwngcQ6Lr3VN0T4oR8eT5LrTUAALTjF7OuAAAAkxHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhxw4ZVSNkspa7OuB8BpKfaBAy66UsqjJL1a69VzUJelJBtJFpMsJzlI8lWSjVrrwSzrBrTjzVlXAOAslVJWkiwkWS6lLNVaezOuy2oGwlopZSHJ/SSPSik3a61fzqp+QDsMoQIX3c30e7yS5MasKtEFtZu11pd62mqtB7XW1SS9JPe7HjqAYwlwwEV3kGS7++/1GdbjTpLNY+5vHPkKMJIAB1xY3cKFz7ser50kS6WU5RlVZznJw2Oev9d9XZlSfYCGCXDARXa11noYjLa6r7PshVtIcm3Evf1pVgRomwAHXEjdnLOvD18PLA743WxqlJvpB8rtEfcP577NbJEF0A4BDrio1pIcXdH5ZZKFbjXoVHWLFfaOKXKr+7p1TBmAJAIccHF9NGTLkPMwjDrKWvp71dlGBHglAQ64cLqtOB4evV5r3Ul/Veqx24mUUm6XUuop/FsYs75b6c+PWz3B/1xgDglwwEW0fsxcs+0kKaWMDHG11nu11nIK/155skI3nLuWZHWWmwwDbXGUFnDhlFK+zujFAEvdv51uA92Z6Xro/prkU0OnwCQEOOBCOTw667hA1AW8pSRXZnn+aCnlYZK7whswKUOowEUzznmih4sZ1s66MqOUUu4n2RLegJMQ4IB5dDg/7tawm2e9iKGUspnkwbB5erPY4gRojwAHXBjdwoT7ryrXDZvuJVkednj8WS5i6I73+vaYRRZWogKvJMABF8mtbquQcUx9T7iB+Xn3jik21tYjwHyziAG4ELqera2MuTChO1T+YZKDWuuVKdTvcG+6UQFzIf2FFVuvCHgAeXPWFQB4HV2v1v286Ln6aynl2G05SikP8uJQ+YVSyqP0tx25eYZ7sR3W8dhNhOMsVGAMeuAAABpjDhwAQGMEOACAxghwAACNEeAAABrT9CrUd999t37wwQdn/py///3v+eUvf3nmz2F+aWNMg3bGWdPGTtfDhw+/qbX+t2H3mg5wH3zwQb766qszf87u7m6uX79+5s9hfmljTIN2xlnTxk5XKeX/jLpnCBUAoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0JimN/IFAE7u6dOn2d/fz+PHj/Pjjz++9uf9+te/zp///OdTqNnF88Ybb+Ty5ctZXFzMpUuXXvvzBDgAmENPnz7N3/72t1y5ciUffPBB3nrrrZRSXuszHz9+nMuXL59SDS+OWmt++OGHfPfdd/nb3/6W999//7VDnCFUAJhD+/v7uXLlSt599928/fbbrx3eGK2Ukrfffjvvvvturly5kv39/df+TAEOAObQ48eP86tf/WrW1Zg7v/rVr/L48ePX/hwBDgDm0I8//pi33npr1tWYO2+99dapzDcU4ABgThk2nb7T+p4LcAAAjbEKdQz/7/HT/OHBX56//v3qxzOsDQAw7/TAAQA0Rg8cAPAzgyNP43r27Gnefvv1N6kdx7yPhumBAwDmzsbGRq5evZpSStbX17Ozs/P83vb2dlZXV1NKyc2bN7O9vT3Dmg6nBw4AmDubm5vZ3t7OxsZGtra2Xrq3traWa9eu5erVq/nss8+ysLDw0v29vb3s7OxkaWkp+/v7WVpaysrKyjSrL8ABAIyr1+vl7t27uX///vNrN2/ezOLiYpaXl6dWD0OoAABj2tzczPr6+kvX7ty5k42NjanWQ4ADABjTF198kaWlpZeuLS0tvTSHbhoEOACAMfR6vRwcHPwswB3Okdvb25taXcyBAwDm1sHBQe7du/ez619//fXQssfZ398/tXq9igAHAMythYWF3L59+2fX9/b2zuX2IYcMoQIATOBVPXHTIMABAIzhcK7b0aHSw0C3uLg4tboIcAAAY1haWsrCwsLPeuAOA5194AAAzqGVlZX0er2XrvV6vamfxCDAAQCMaXNzM3fv3n3p2tbWVjY3N6daD6tQAYCf+f3qxxO/5/Hjx7l8+fIZ1Ob0bWxsZGdnJwcHB1lfX8/Nmzef96Jtb28/Pyrr008/za1bt3Ljxo0k/WHUw3NUl5aW0uv1sr6+PtXh00SAAwDm0HE9Zmtra1lbWxt5f9rDpcMYQgUAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGvDnrCgAA59Cf7k78lrefPU3evnQGlRnikzvTec45pQcOAJg7GxsbuXr1akopWV9fz87OzvN729vbWV1dTSklN2/ezPb29gxrOpwABwDMnc3Nzayvr2dhYSFbW1tZWVl5fm9tbS2bm5tJks8++yxra2tDP2NjYyN7e3tTqe9RhlABAMbU6/WyubmZhYWF5z11szBRgCulrCX5aODSt7XWeyPKLSQ56L72aq1fHvOZY5UFAJilpaWlbG1tJUm+/HJ2cWXsAFdKuZ/kX2utG4PXSilrtdbtgWubORLsSimbpZSlo2FvkrIAAPSNFeC6oLU/JFStJPl8oNxSkrVa65XBQrXWjVLKo1LKdq31YNKyAAC88MoA1wWt23l56DRJcjR8JVlPsnO0XKeXZC3JvROUBQA4dQcHB7l37+dx4+uvv55BbcY3Tg/cevrz0npjlF3J8aFsNS9C2SRlAQBO3cLCQm7fvv2z63t7e+dy+5BD42wjspJkL0lKKQullBtdr9wwy0m+HXGvl+TaCcsCANAZJ8AtJ9kvpdzIi16zhW4Bw6ggN8rCGZUFAJgb465CXUxyUGs9HPLcK6VsJHlYSvmw1npQSjkMXK9ceDBJ2SHvXUt/flzee++97O7uTvoRE3vrp6f5x+//+vz17u5/nPkzmS9PnjyZSltmvmlnDPr1r3+dx48fj7z/9rOnE39m/anm6QnedxLPjqn7uL7//vskGfp9+Pvf//783htvvDH0/T/99FP+8z//89jv46jnvu7/F48NcANBa3kgvCVJaq29UkovyZ0kGz978xnptizZTpJr167V69evn/kzv/jnf8m/v/Ph89e/u/7xmT+T+bK7u5tptGXmm3bGoD//+c+5fPny6AInONP06bOnuTSls1AvHVf3Mb3zzjtJMvT78Mtf/vL5vVHfp1/84hf5h3/4h+O/jyOe+9vf/nbC2h559nE3B7bxGHVORC/JjSNlXzn0OUlZAABeNs4Q6kGS/RH39pOMOw/u8LSF0y4LAJy2T+5M/JZnjx+fSs/YNGxsbGRnZycHBwdZX1/PzZs3n5+Hur29nfv37ydJPv3009y6dSs3btxI0t925O7du+n1eun1etnY2MjKykpWV1dfOk/1rI0T4Hrpz4EbZTBo7WTIfnGdxSRfnbAsAMCpOTysfpi1tbWRB9gvLCwc+95pGWcV6ufpr0QdZikvB629jO6RW0py/4RlAQDojBPgtpMsjdgyZCXJYAzdypD927rFEMtJvjhhWQAAOq8McN2Cg430A9dzpZStJNuDq1O70xq2u7NTB20m2Rg823SSsgAAvDDWPnC11nullIMutB162G3pcbTsRillrZRyO/35cR8leVBr/fJ1ygIA0DfuRr4ZFtamXRYAgPHmwAEAcI4IcAAwp2qts67C3Dmt77kABwBz6I033sgPP/ww62rMnR9++GHk2aqTEOAAYA5dvnw533333ayrMXe+++67ic9OHUaAA4A5tLi4mEePHuWbb77Js2fPDKeeoVprnj17lm+++SaPHj3K4uJxB1yNZ+xVqADAxXHp0qW8//772d/fz7/927/lxx9/fO3P/P777/POO++cQu0unjfeeCOXL1/O+++/n0uXLr325wlwADCnLl26lN/85jf5zW9+cyqft7u7m9/+9ren8lkczxAqAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0Jg3T/rGUsqDWuvqiHtrSRaSHHRfe7XWL1+3LAAAJwxwpZQbSVZG3NtM8m2t9d7gtVLK0uC1ScsCANA3cYArpSwkuTXi3lKStVrrlcHrtdaNUsqjUsp2rfVg0rIAALxwkjlwa0k+H3FvPcnOiHu97r0nKQsAQGeiAFdKWU4/dI3qGVtJP3wN00uyesKyAAB0Ju2BW6m17h1zfznJtyPu9ZJcO2FZAAA6Ywe4brXo9ms+b+GMygIAzI2xFjF0Cw72j1tU0C1uSEYPr56o7JD3rqWbH/fee+9ld3d30o+Y2Fs/Pc0/fv/X5693d//jzJ/JfHny5MlU2jLzTTvjrGlj0zPuKtQb52Vbj1rrdrqewGvXrtXr16+f+TO/+Od/yb+/8+Hz17+7/vGZP5P5sru7m2m0ZeabdsZZ08am55VDqN2eb6/cWHegd+6VQ5+TlAUA4GXHBrhuqHOx1jpqtegkDk9bOO2yAABz5VVDqGtJ/qmUcvXI9aUkKaVsda/v11p30t9i5KMRn7WY5KuB15OUBQCgc2yAGzXv7fAorVrr+pFbe+lvDzLMUpKtE5YFAKBz4sPsR9hK8vDoxW4odjnJFycsCwBA5yRHaSX9Ic7B7UCSJN1cue3ukPpBm0k2BrchmaQsAAAvTNQD1w2d3kr/GKwk+WMppZfk08PA1R1Gv1ZKuZ3+QoSPkjyotf5sJeskZQEA6JsowHXBapwtRcY+sWGSsgAAnHwIFQCAGRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0Jg3Z10BgFPxp7uzrsGrfXJn1jUALgg9cAAAjRHgAAAaI8ABADTGHDiAaZn1PL0nHx5fB3P0oBl64AAAGiPAAQA0RoADAGiMAAcA0JixFzGUUm4n+a9JlrtLW7XWL0eUXUuykOSg+9o7jbIAAIwZ4Eopm0nu1loPutfLSf5YSrlVa705pOy3tdZ7g9dKKUuD1yYtCwBA3ysDXNdDtnUY3pKk1rpXSrmZ5EEp5cZhj1kpZSnJWq31yuBn1Fo3SimPSinbAyFw7LIAALwwzhy4j2qtvaMXa6076Q97rg9cXk+yM+JzeknWTlgWAIDOOAHudinl/oh7XyW5NvB6Jf3wNUwvyeoJywIA0BknwI3qJUuSxST7A6+Xk3w7omwvL4e9ScoCANB55Ry4WutxPWHLSbYneN7CGZUFAJgbJz4LtVvckCSb3evDwPXKhQeTlB3x3LUkee+997K7uzvpR0zsrZ+e5h+//+vz17u7/3Hmz2S+PHnyZCpt+UJ78uGsa3DuPfnpUnaP+z5pg7wmP8um50QBrgtgm0luDlvgcJZqrdvpev2uXbtWr1+/fubP/OKf/yX//s6LH3q/u/7xmT+T+bK7u5tptOULbdYHxTdg98mHuf5f/jq6wPX/Mb3KcCH5WTY9Jz2J4X6SjcENdwe2/Hjl0OckZQEAeNnEAa7bfPd+1xM2icPTFk67LADAXJloCLU7TuvrY8LbTpKPRtxbTH/bkZOUBWbNECXAuTF2D1wp5UaSg6PhrZSyMvByL8nSiI9YSn/o9SRlAQDojBXgDkPaiJ63wW1GtjJk/7Zu0cNyki9OWBYAgM44Z6Eup3/s1eddL9yhxfTnqj1fiFBr7ZVStkspm7XWjYGym+kvejg4SVkAAF4YZw7cH9MPaTdG3B8MX4eH0a918+UO0p/n9mBwxepJygIA0DfOSQxXJv3QSVaonmA1KwDAXDvpPnAAAMyIAAcA0BgBDgCgMSc+zH6e/eHBX156/ftVZ6MCANOjBw4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjXGYPQB9f7o76xq82id3Zl0DOBf0wAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAY96cdQUAYGx/ujvrGhzvkzuzrgFzQg8cAEBjBDgAgMYIcAAAjRHgAAAaYxEDnAeP/+/5n5wNwLkhwJ2CPzz4y0uvf7/68YxqAgDMA0OoAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMa8OesKAMCF8ae7s67Bq31yZ9Y14BQIcGfgDw/+8tLr369+PKOaAAAXkSFUAIDGCHAAAI0R4AAAGiPAAQA0xiKGKbCoAQA4TQIc8+HcL+3/cNYVAKAhhlABABoz8x64UspakoUkB93XXq31y9nWCgDg/JppgCulbCb5ttZ6b/BaKWVp8BoAAC/MLMCVUpaSrNVarwxer7VulFIelVK2a60HM6oeAMC5NcseuPUkOyPu9ZKsJbmQvXBWpQIAr2OWAW4lxwe41VzQAHeUQAcATGKWAW45yecj7h32wM094Q6AuXLut31K8smdWddg9qtQj7Ew6wrMytHQdtw9gQ6AiZxlQHryYRsB7AKYSYArpRyGs4kXKXTbjhz2zj0ppfzvU6vYaO8m+WYKz5nY/5x1BTgt57aNcaFoZ5y1OWlj/2taD/rvo26c5x64oWqt20m2p/nMUspXtdZr03wm80UbYxq0M86aNjY9MzmJYWB7kLkdJgUAOKnzepTW4ckMAAAcMcsAt5PkoxH3FpN8NcW6vMpUh2yZS9oY06Cdcda0sSkptdbZPLh/jNZyrXV1yL2HSba6+W4AAAyYZYBbSvLw6FFa3QrVR0muOEoLAODnZhbgkue9cKm1bgxc20rytcPsAQCGm2mAS57v63a4aOGjJP9aa/1yppXqHKnbQpLeeakbbSml3E7yX9M/gSTpTxEY2pa0O05LKeXBsGkq3T3tjBPp2s7gHPZvh3W6aGNna+YB7rzqegdfapTDrsGrdO3m7uGUgFLKcpI/Jtmptd4cUla747WVUm4kuV9rLUPuaWecSCnlfvodLfeOXHswOG9dGzt7AtwQo+bndfceJfnQ/DzG0f0FulNr7R25vpLkQZKbh3+Raneclm4u8WdJbhwNcNoZJ9UFsIVa6/qR64+SfOpn2XSd133gZm09/W1OhunlxVFe8CofHQ1vSVJr3Ul/WGHwB6F2x2lZS/L5iHvaGRPrQtntJJtH79VarxwZGtXGpkCAG24l/UY2TC/J0DklMMTtbnhhmK+SDB45o93x2roh+sM/EIbRzjiJ9fTnsI1qO4O0sSkQ4IZbTvLtiHu9vPxLF44z6q/QpL9h9f7Aa+2O07BSa9075r52xkmsJNlL+kP0pZQbXa/cMNrYFAhwJ+MMV8ZSa109ulBhwGFPybi0O47Vzbl83Q3QtTOGWU6y3y2OWUn/Z9dCKeX+MUFuFG3sFAhwR3STfxNnsXKGul+0STefRLvjdXW/RPePmxyunfGaFpMc1Fq/rLUedD29G0keHrYtbWx6BDiYsu4H3Gb6K1DHmU8C47hhjy3OwkAoW+4WYD3X/QzrJbkz9YrNOQHuiIG/XnXxclbuJ9kY/GWr3fE6umGtV4Y37YyTGGg3o+ZW9pLcOFJWGztjAtzkDneVhol1+yjdH9zwckzaHUN1vSOLp9Sbq50xykFeXnQ1aD/JuPPgtLFT8uasK3BO7eTlY0IGLaa//QNMpDtO6+tjwpt2x0msJfmnUsrVI9eXkufnSyf9Pxx2op1xMr3028cog6FMG5sCAW64vbw4s/KopSRbI+7BUN0Q18HR8FZKWRmYU6LdMbFRxxIdrhY8umt+tDNO5vO8vPH4oKW8HMq0sSlwlNYQo44B6YYqHiW54hgQxtUdm7UwbIJ5KWWz1rrR/bd2x6kZdRaqdsZJDLSPn50uU0qpSVYP/xjVxqZDgBuhm6uUw1+u3bWt9IfAHMTLWLpd8e/k58caLaY/F+SjwR4S7Y7T0m1Vs5Uhvyy1M06imwayWmtdHbi2lSRDzkfVxs6YAHeM7gfg4YTLj5L8q2X6TKI7uPm41VgbR3+YaXe8jq7n7Vb6m60upD+c1Uv/sPGDgXLaGRPr2s3gfMuHo+b1amNnS4ADAGiMbUQAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA05v8DCj/7aLdmw7kAAAAASUVORK5CYII=\n",
+ "text/plain": [
+ "<Figure size 735.243x454.405 with 1 Axes>"
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnAAAAGbCAYAAABEeK1/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAaS0lEQVR4nO3dwW9bx70v8N88x3aAPqe06osgKF5eQgPZV/YF3ts5gPQfWPX7B0KtuzHh9w8Y8qZrMX9BbO+6Kqwg2t1FbW2LW8BsX4B2cRHLQuwCsYN03oKHDi2TEimShxrx8wECmeeMyMlo4nw1c2Ym5ZwDAIBy/LdFVwAAgMkIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOOPNSSlsppdai6wEwK8k+cMBZl1J6HhHdnPO1RddlUErpZkSs55w3F10XoCxG4IAzLaW0FhGNiFhNKTVPQX2aKaWbKaUHEfFg0fUByiTAAWfdRkS0qz/fXGRFqhG3rYhoRsQXi6wLUDYBDjjrDiKiU/15oVOVOeeHOeeNnPO9nPPBIusClE2AA86sauHCV1VY2omIZkppdcHVApiaAAecZddyznvVn7errxYMAMUT4IAzKaXUiIin/dc554fVH3+7mBoBzI4AB5xVrYh4eOjaw4hoVCtTAYolwAFn1dWcc/fQNdOowJkgwAFnTrXf25PD13POO9FblXrkdiIppdsppTyDfxpz+lcElpwAB5xFmznnzoh7nYg3e7INVW3zkWbwj61CgLl4b9EVAJiDm0dsF9I/jWEz3n1GDqAIAhxwplQLFNoDq06HlXkaEWsppYZRMqBEplCBs2bjqPBW6S9maM27MgDzIMABy6j/fNytYTctYgBOOwEOODOqhQkPjitXTZvuRcRqtWL18H2LGIBTTYADzpJb1VYh4zgte8KtLPjzgQIJcMCZUB1cf3OCacvH1deFPAc3UE/TrMDEBDigaCmltZTS8/h5RO2vR+3xVn3Po4j4unrZSCk9Tyk9GTadOoe6PqpWwT6vLq8NfP6jeX4+cHaknPOi6wAAwASMwAEAFEaAAwAojAAHAFAYAQ4AoDBFn4V65cqV/Mknn8z9c/75z3/GL37xi7l/Dj3au17au17auz7aul7ae/aePHnyXc7534bdKzrAffLJJ/H48ePjC05pd3c3bty4MffPoUd710t710t710db10t7z15K6f+NumcKFQCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQmKI38gUATu7Vq1exv78fL168iJ9++mmq9/rlL38Zf/7zn2dUs7Pn3LlzcenSpVhZWYmLFy9O/X4CHAAsoVevXsW3334bly9fjk8++STOnz8fKaUTv9+LFy/i0qVLM6zh2ZFzjh9//DG+//77+Pbbb+Pjjz+eOsSZQgWAJbS/vx+XL1+OK1euxIULF6YKbxwtpRQXLlyIK1euxOXLl2N/f3/q9xTgAGAJvXjxIj744INFV2PpfPDBB/HixYup30eAA4Al9NNPP8X58+cXXY2lc/78+amfN4wQ4ABgaZk2rd+s2lyAAwAojFWoY/ivF6/i94/+cmy5361/VkNtAIBlZwQOAKAwRuAAgHeMM/M06PXrV3HhwvQb1I5r2We9jMABAEun3W7HtWvXIqUUm5ubsbOz8+Zep9OJ9fX1SCnFxsZGdDqdBdZ0OCNwAMDS2draik6nE+12O7a3t9+612q14vr163Ht2rX48ssvo9FovHV/b28vdnZ2otlsxv7+fjSbzVhbW6uz+gIcAMC4ut1u3L17Nx48ePDm2sbGRqysrMTq6mpt9TCFCgAwpq2trdjc3Hzr2p07d6LdbtdaDwEOAGBM9+/fj2az+da1ZrP51jN0dRDgAADG0O124+Dg4J0A139Gbm9vr7a6eAYOAFhaBwcHce/evXeuP336dGjZo+zv78+sXscR4ACApdVoNOL27dvvXN/b2zuV24f0mUIFAJjAcSNxdRDgAADG0H/W7fBUaT/Qrays1FYXAQ4AYAzNZjMajcY7I3D9QGcfOACAU2htbS263e5b17rdbu0nMQhwAABj2trairt37751bXt7O7a2tmqth1WoAMA7frf+2UTlX7x4EZcuXZpTbWav3W7Hzs5OHBwcxObmZmxsbLwZRet0Om+Oyvriiy/i1q1bcfPmzYjoTaP2z1FtNpvR7XZjc3Oz1unTCAEOAFhCR42YtVqtaLVaI+/XPV06jClUAIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKMx7i64AAHAKfXN3ouIXXr+KuHBxTpUZ4vM79X3WKWQEDgBYOu12O65duxYppdjc3IydnZ039zqdTqyvr0dKKTY2NqLT6SywpsMJcADA0tna2orNzc1oNBqxvb0da2trb+61Wq3Y2tqKiIgvv/wyWq3W0Pdot9uxt7dXS30PM4UKADCmbrcbW1tb0Wg03ozULYIABwAwpmazGdvb2xER8fDhw4XVwxQqAEBhBDgAgMKYQgUAltbBwUHcu3fvnetPnz5dQG3GJ8ABAEur0WjE7du337m+t7d3KrcP6TOFCgBQGAEOAKAwAhwAQGEEOACAwghwAACFsQoVAHjX53cmKv76xYu4eOnSnCoze+12O3Z2duLg4CA2NzdjY2PjzXmonU4nHjx4EBERX3zxRdy6dStu3rwZEb1tR+7evRvdbje63W602+1YW1uL9fX1t85TnTcBDgBYOv3D6odptVojD7BvNBpHfm9dTKECABRGgAMAKIwABwBQGAEOAKAwAhwAQGEmWoWaUmpFxNWBS89yzvdGlGtExEH1tZtzfnjEe45VFgCACQJcSulBRPwp59wevJZSauWcOwPXtuJQsEspbaWUmofD3iRlAYDZyjlHSmnR1VgqOeeZvM9YAa4KWvtDQtVaRHw1UK4ZEa2c8+XBQjnndkrpeUqpk3M+mLQsADBb586dix9//DEuXLiw6KoslR9//DHOnTs39fsc+wxcFbRuR8Q7u9blnC8fmu7cjIidEW/VjYjWCcsCADN06dKl+P777xddjaXz/fffx6UZnFgxziKGzeg9l9Ydo+xa9MLXMN2IWD9hWQBghlZWVuL58+fx3XffxevXr2c2tce7cs7x+vXr+O677+L58+exsrIy9XuOM4W6FhF7EREppUb/9YhAtxoDU6qHHB5Vm6QsADBDFy9ejI8//jj29/fjb3/7W/z0009Tvd8PP/wQ77///oxqd/acO3cuLl26FB9//HFcvHhx6vcbJ8CtRsTjlNLN6vVORDSrRQ3tMUfm+hpzKgsATOjixYvx0UcfxUcffTT1e+3u7sZvfvObGdSKcYy7CnUlIg5yzv1n1vZSSu2IeJJS+jTnfFCNzkX0tgM50iRlh3xvK6rRuQ8//DB2d3cnfYuJnf/Xq/j1D389ttzu7j/mXpdl8PLly1p+rvRo73pp7/po63pp73odGeAGgtbqQHiLiIicczel1I2IOxHRfueb56TasqQTEXH9+vV848aNuX/m/T/8Mf7+/qfHlvvtjc/mXpdlsLu7G3X8XOnR3vXS3vXR1vXS3vU6chHDwDYeeyOKdCPi5qGyx059TlIWAIC3jbMK9SAi9kfc24+I5pif1T9tYdZlAQCWyjgBrhu9Z+BGGQxaO/H2UVuDViLi8QnLAgBQGSfAfRW9lajDNOPtoLUXo0fkmhHx4IRlAQCojBPgOtHbNmRY2FqLt09o2I6I64cLVYshViPi/gnLAgBQOTbAVQsO2tELXG+klLYjojO4OrXaE65TnZ06aCt6e8YdnKQsAAA/G2sfuJzzvZTSQRXa+p5UW3ocLttOKbVSSrej93zc1Yh4dOjM1InLAgDQM+5GvjEsrNVdFgCA8Z6BAwDgFBHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMK8d9JvTCk9yjmvj7jXiohGRBxUX7s554fTlgUA4IQBLqV0MyLWRtzbiohnOed7g9dSSs3Ba5OWBQCgZ+IAl1JqRMStEfeaEdHKOV8evJ5zbqeUnqeUOjnng0nLAgDws5M8A9eKiK9G3NuMiJ0R97rV956kLAAAlYkCXEppNXqha9TI2Fr0wtcw3YhYP2FZAAAqk47AreWc9464vxoRz0bc60bE9ROWBQCgMnaAq1aLdqb8vMacygIALI2xFjFUCw72j1pUUC1uiBg9vXqiskO+txXV83Effvhh7O7uTvoWEzv/r1fx6x/+emy53d1/zL0uy+Dly5e1/Fzp0d710t710db10t71GncV6s3Tsq1HzrkT1Ujg9evX840bN+b+mff/8Mf4+/ufHlvutzc+m3tdlsHu7m7U8XOlR3vXS3vXR1vXS3vX69gp1GrPt2M31h0YnTt26nOSsgAAvO3IAFdNda7knEetFp1E/7SFWZcFAFgqx02htiLi31NK1w5db0ZEpJS2q9cPcs470dti5OqI91qJiMcDrycpCwBA5cgAN+q5t/5RWjnnzUO39qK3PcgwzYjYPmFZAAAqJz7MfoTtiHhy+GI1FbsaEfdPWBYAgMpJjtKK6E1xDm4HEhER1bNyneqQ+kFbEdEe3IZkkrIAAPxsohG4aur0VvSOwYqI+Dql1I2IL/qBqzqMvpVSuh29hQhXI+JRzvmdlayTlAUAoGeiAFcFq3G2FBn7xIZJygIAcPIpVAAAFkSAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKMx7i64AzMU3dxddg9n5/M6iawDAKWMEDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYd5bdAWAY3xzd77v//LT+X9G3+d36vkcgDPOCBwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCjH2UVkrpdkT8KiJWq0vbOeeHI8q2IqIREQfV1+4sygIAMGaASyltRcTdnPNB9Xo1Ir5OKd3KOW8MKfss53xv8FpKqTl4bdKyAAD0HBvgqhGy7X54i4jIOe+llDYi4lFK6WZ/xCyl1IyIVs758uB75JzbKaXnKaXOQAgcuywAAD8b5xm4qznn7uGLOeed6E17bg5c3oyInRHv042I1gnLAgBQGSfA3U4pPRhx73FEXB94vRa98DVMNyLWT1gWAIDKOAFu1ChZRMRKROwPvF6NiGcjynbj7bA3SVkAACrHPgOXcz5qJGw1IjoTfF5jTmUBAJbG2NuIHFYtboiI2Kpe9wPXsQsPJik74nNbEREffvhh7O7uTvoWEzv/r1fx6x/+emy53d1/zL0uy+Dly5fT/1xffjqTuiyDl/+6GLt1tVcN/72edjPp34xFW9dLe9frRAGuCmBbEbExbIHDPOWcO1GN+l2/fj3fuHFj7p95/w9/jL+/f/z/4H5747O512UZ7O7uxtQ/12/uzqQuy2D35adx478f/wvKTNz4P/V8zik2k/7NWLR1vbR3vU56EsODiGgPbrg7sOXHsVOfk5QFAOBtEwe4avPdB9VI2CT6py3MuiwAwFKZKMBVx2k9PSK87UTE1RH3VqK37chJygIAUBk7wKWUbkbEweHwllJaG3i5FxHNEW/RjN7U60nKAgBQGSvA9UPaiJG3wW1GtmPI/m3VoofViLh/wrIAAFTGOQt1NXrHXn1VjcL1rUTvWbU3CxFyzt2UUieltJVzbg+U3YreooeDk5QFAOBn42wj8nX0QtrNEfcHw1f/MPpW9bzcQfSec3s0uGL1JGUBAOgZ5ySGy5O+6SQrVE+wmhUAYKmddB84AAAWRIADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKMw4Z6GyLL65u+ga9Lz89PTUhdk6Sz/Xz+8sugbAEjMCBwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFMY+cDP0+0d/Gbvs79Y/m2NNAICzzAgcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAw7y26AgBF+ubuyb7v5acn/955+fzOomsATMgIHABAYQQ4AIDCCHAAAIXxDNyC/P7RX8Yq97v1z+ZcEwCgNEbgAAAKYwRuCv/r2878P+SbX83/MwCAohiBAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCWIV6yv1H99nYZf9304pVAFgGRuAAAAojwAEAFEaAAwAojGfgAJbdN3cXXYPZ+vzOomsAc2cEDgCgMAIcAEBhBDgAgMIIcAAAhbGIAYCzpb8o4+Wn5S/QsCCDERYe4FJKrYhoRMRB9bWbc3642FqVadxTG5zYAABlW2iASyltRcSznPO9wWsppebgNQAAfrawAJdSakZEK+d8efB6zrmdUnqeUurknA8WVD0AgFNrkSNwmxGxM+JeNyJaEWEUbg7GnWqdhGlZAKjPIlehrkUvqA3TjYj1GusCAFCMRY7ArUbEVyPu9UfgWGKTjBQaAQTOpJJW0Z6FVb/jOgWrgxe+CvUIjUVXgPHNclr2n7/8H/Ef/zX7aV4AOCsWEuBSSv1wNvEihWrbkf7o3MuU0n/OrGKjXYmI72r4HHq0d720d720d320db2WqL3/b10f9D9H3TjNI3BD5Zw7EdGp8zNTSo9zztfr/Mxlpr3rpb3rpb3ro63rpb3rtZBFDAPbg5gmBQCY0Gk9C7V/MgMAAIcsMsDtRMTVEfdWIuJxjXU5Tq1Ttmjvmmnvemnv+mjremnvGqWc82I+uHeM1mrO+Z393lJKTyJiu3reDQCAAYsMcM2IeHL4KK1qherziLjsKC0AgHctLMBFvBmFi5xze+DadkQ8dZg9AMBwCw1wEW/2desvWrgaEX/KOT9caKUqh+rWiIjuaanbWZJSehC90ze2c87danR2LSKu5Zw3F1u7slUj2lsR8eiovquvz8Y47a2/z05K6XZE/Cp6J/tE9Np0VLvr41Mat7318XosPMCdVtXo4LPBkcBh15he9czj6qHLO8Oej2Q81Uj2SkT8KXqBYuOI/7Hp61OasL319xmo+ujd/qM2KaXViPg6em25MaSsPj6FCdtbH6+BADfEqOfzqnvPI+JTz+fNTvUXw6OIaFaXHuec9xZYpTMlpZRjRKDQ12fvqPau7uvvU6pG03Zyzt1D19ei17Zv2l8fn94k7V1d18drUNxJDDXZjN42J8N0o3eUl9/aZudZznlUezNf+nr99PfpXR22S0HOeSeldBC9ft0PFPr49CZp7wh9vBandSPfRVuL3n/Yw3QjwjAwZ4W+ToluV89ZDfM4IgaPc9LHpzdJe1MTAW641Yh4NuJeN3TWuUgpNVJKa9WUB/XQ1xdEf5/KUaM7KxGxP/BaH5/eJO39hj4+XwLcyTjDdbZ+VT1jcT16v801UkqPqhV9LJafwezp71PKOa8ffnB+wGocHTgO0+7HOEF76+M18AzcIQMdzEOt9bo/8CDxXvUQ7JMYfdwaU9LXF0p/n4MqNET0VgLr43N2uL0P0cfnzAgcC5dzbh9eBVY9ALsy8BcEnAn6+3wM7MG3cXi1JLN3VHvr4/UQ4A4Z6HSGehevGxGjhu2Zkr5+6ujv03kQEe3B7Sz08bl6p73HoI/PkAA3uf5O3sxfN37eR4j66ev10t9PqJqeezBsq4tj6OMnMEV76+MzJMANtxOj5+lXovdQJjOQUnrSPxOXhdDXa6S/z151vNPTI8KEPj5Dx7W3Pl4fAW64vRj9W0IzekPHzEYjIp6OuNeM3s+C+dHX66W/z1BK6WZEHBwOE9UJAX36+IyM2d76eE0EuOG2Y8jeQNVDm6sRcb/2Gp1d28N+k6v2DVqN3s+C+dHX66W/z0g/NIwYCRrcnFcfn4FJ2lsfr4ezUEfoDwHnnNsD17ajN3Ts2JUZqf4SvTPYztX1RxGxd/g6k6vO5tw8YspDX5+ho9pbf5+N6iD1OxHx1aFbK9EbAbqac94cKK+PT2GS9tbH6yPAHaFa7tx/yPVqRPxpwhU3jKH/H3z1shG9vxS2naV3ctVzKv8evd94m9Hrw4+j9xDxO0v89fXpTNLe+vv0qkPoj1pZ2j4czPTxk5u0vfXxeghwAACF8QwcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwvx/HkgjrisWxs4AAAAASUVORK5CYII=\n",
+ "text/plain": [
+ "<Figure size 735.243x454.405 with 1 Axes>"
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnAAAAGbCAYAAABEeK1/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAczklEQVR4nO3dT2+cx4En4F+t4iiYgWyK9sLIYbxOC8h5hlKAyU3BUt9AWs0XMHkOsBDh/QICddkz6etcbOk6JzGwbhkgEjG3YAdw76xPi8Ga5lgKEMMxag/9Um5R3VR3s8nuMp8HMGi+b7G7VCS7f6y/pdYaAADa8Z8WXQEAAKYjwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4oFmllO1Sysai6wFw3op94IBWlVK+TtKvtV5fdF2SpJRyL8m1JF8keTdJaq1bc3jc20lu1Vo3T/tYwI/DTxZdAYBZlFLWk6wkWSul9Gqt/QXX53GSh8Mhq5SyXkp5NkvALKX0kqwluZvkdpLduVUWaJ4hVKBVd5Ic9W7dXmRFSinbSVJrfSVk1Vr3kjwtpTyc8vFuJ9lO0kvy0bzqCfx4GEIFmtSFpvtJjoZRry2wLjXJnVrroxH31pI8S3K11np4isffNYQKHNEDBzSnW7jwaReI9pL0uqC0iLoc9f7tjylyNLRrsQUwNwIc0KLrtdajwLTTfVxU79StJBk3B2+o1+3WudUI+NET4ICmlFJWMljlmSQZGrb8b4upUXoTlOlPWA5gIgIc0JqNJMfnmj1KstKtTD1vq3MuB/BGAhzQmmsjhisXOYy6kmSSxQkrZ10R4OIQ4IBmdHujPTt+vduu4zBv2E6klHKvlFLn8N9wGJukZ22m1acA4whwQEs2j++1NmQ3eWVV6GtqrQ9qrWUO/00byPS+AXPlJAagJbdP2C7kaJHAZl6fI3eWDjLZAgW9cMDcCHBAE7oFClujNssdKvNFkvVSysqsm+bOYJLnWc0g6AHMhSFUoBUjTzo45mgxw3lumvs0ebm9yTgrGb/RL8DUBDjgx+RoftzdUTfPaBHD0aKKkcOoQ2Ufz+efCCDAAQ3oFia88UD4bth0P8lat2L1+P2zWMTwWffxxphqHe1NtzfFPxngRAIc0IK73VYhkzjXPeG6MPfghOe7m8FB9K8dtTXl+a02AgZeEuCApdYdXH/7DXPMhj3tPp7bPLha61bysq4vdT2Ha0m2jn9NKeVekmellJ3j946VO/p324oEeEmAA5ZSKWW9lPJ1fuhR+98n7fHWfc3jJL/rPl0ppXxdSnk2ajh13mqt17vn3Onm2m0n+VWt9dqYFbH7GaxgfW1j4u7f/rhbVft1d3l96N9jPh1ccKXWuug6AAAwBT1wAACNEeAAABojwAEANEaAAwBoTNNnob733nv1ww8/PPPn+dOf/pS//uu/PvPnuUi06fxp0/nSnvOnTedPm87fMrXps2fP/l+t9T+Putd0gPvwww/z9OnTNxc8pSdPnuTmzZtn/jwXiTadP206X9pz/rTp/GnT+VumNi2l/J9x9wyhAgA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAa0/RGvgDA7L799tscHBzk+fPn+f777xddnaXwzjvv5I9//OPcH/fSpUu5cuVKVldXc/ny5VM/ngAHABfQt99+my+//DJXr17Nhx9+mLfeeiullEVXa+GeP3+eK1euzPUxa6357rvv8s033+TLL7/MBx98cOoQZwgVAC6gg4ODXL16Ne+9915++tOfCm9nqJSSn/70p3nvvfdy9erVHBwcnPoxBTgAuICeP3+et99+e9HVuHDefvvtPH/+/NSPI8ABwAX0/fff56233lp0NS6ct956ay7zDQU4ALigDJuev3m1uQAHANAYq1An8fz/Jp/fX3Qt5uc3Hy+6BgDAKeiBAwBojB44AOA1//Pxvy66Cif67a1fLroKC6UHDgC4cLa2tnL9+vWUUrK5uZm9vb2X93Z3d3Pr1q2UUnLnzp3s7u4usKaj6YEDAC6c7e3t7O7uZmtrKzs7O6/c29jYyI0bN3L9+vV88sknWVlZeeX+/v5+9vb20uv1cnBwkF6vl/X19fOsvgAHADCpfr+f+/fv5+HDhy+v3blzJ6urq1lbWzu3ehhCBQCY0Pb2djY3N1+59vHHH2dra+tc6yHAAQBM6LPPPkuv13vlWq/Xe2UO3XkQ4AAAJtDv93N4ePhagDuaI7e/v39udTEHDgC4sA4PD/PgwYOXn3/77be5fPlyvvjii5FlT3JwcDD3+o0jwAEAF9bKykru3bv38vPnz5/nypUr2d/fX8rtQ44YQgUAmMKbeuLOgwAHADCBo7lux4dKjwLd6urqudVFgAMAmECv18vKysprPXBHgc4+cAAAS2h9fT39fv+Va/1+/9xPYhDgAAAmtL29nfv3779ybWdnJ9vb2+daD6tQAYDX/PbWLxddhTO1tbWVvb29HB4eZnNzM3fu3HnZi7a7u/vyqKyPPvood+/eze3bt5MMhlGPzlHt9Xrp9/vZ3Nw81+HTRIADAC6gk3rMNjY2srGxMfb+eQ+XjiLATeBP3/4lv//3rxZdjbF+3Xt30VUAAM6ROXAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIzD7AGA131+f9E1ONlvPl50DRZKDxwAcOFsbW3l+vXrKaVkc3Mze3t7L+/t7u7m1q1bKaXkzp072d3dXWBNRxPgAIALZ3t7O5ubm1lZWcnOzk7W19df3tvY2Mj29naS5JNPPsnGxsbIx9ja2sr+/v651Pc4Q6gAABPq9/vZ3t7OysrKy566RRDgAAAm1Ov1srOzkyR59OjRwuphCBUAoDECHABAYwyhAgAX1uHhYR48ePDy82+//TaXL1/OF198scBavZkABwBcWCsrK7l3797Lz58/f54rV65kf39/KbcPOWIIFQCgMVP1wJVSNpJcG7r0Va31wZhyK0kOu4/9WuvIpRrTlAUAYIoAV0p5mOQPtdat4WullI1a6+7Qte0cC3allO1SSu942JumLAAAAxMFuC5oHYwIVetJPh0q10uyUWu9Olyo1rpVSvm6lLJbaz2ctiwAAD944xy4LmjdS7J9/F6t9eqx4c7NJHvHy3X6STZmLAsAQGeSHrjNDOal9Scou56TQ9mtJA9mKAsAnKfffLzoGpypra2t7O3t5fDwMJubm7lz587L81B3d3fz8OHDJMlHH32Uu3fv5vbt20kG247cv38//X4//X4/W1tbWV9fz61bt145T/WsTRLg1pPsJ0kpZeXo8zGBbi1DQ6rHHO9Vm6YsAMDcHB1WP8rGxsbYA+xXVlZO/NrzMsk2ImtJDkopt/NDr9lKt4ChN+XzrZxRWQCAC2PSVairSQ5rrUdDnvullK0kz0opv6i1Hna9c8lgO5ATTVN2xNdupOude//99/PkyZNpH2Jqf7n0Vzl452/P/Hlm9eTFlPsxn0ObvcmLFy/O5Xt3kWjT+dKe86dN5+80bfrOO+/k+fPn863Qj8D3339/5u3y5z//+dS/Cye+8w8FrbWh8JYkqbX2Syn9JB8n2Xrti89It2XJbpLcuHGj3rx588yf858e/WNW/+Nfzvx5ZvXr3rvTfcHNfzibikzhyZMnOY/v3UWiTedLe86fNp2/07TpH//4x1y5cmW+FfoRODqJ4Sz97Gc/y9/93d+d6jFOHEId2sZjf0yRfpLbx8q+cehzmrIAALxqkjlwh0kOxtw7SDLpPLij0xbmXRYA4EKZJMD1M5gDN85w0NrLq0dtDVtN8nTGsgAAdCYJcJ9msBJ1lF5eDVr7Gd8j10vycMayAMCc1VoXXYULZ15tPkmA203SG7NlyHpePaFhJ8mN44W6xRBrST6bsSwAMEeXLl3Kd999t+hqXDjfffddLl26dOrHeWOA6xYcbGUQuF4qpewk2R1endpt7rvbnZ06bDvJ1vDZptOUBQDm68qVK/nmm28WXY0L55tvvpnLKteJNhCrtT4opRx2oe3Is25Lj+Nlt0opG6WUexnMj7uW5PGxM1OnLgsAzM/q6mq+/PLLJMnbb7+dt956K6WUBdfqx6nWmu+++y7ffPNNvv7663zwwQenfsyJd4AdFdbOuywAMB+XL1/OBx98kIODg/zbv/1bvv/++0VXaSn8+c9/zs9+9rO5P+6lS5dy5cqVfPDBB7l8+fKpH2/KLfwBgB+Ly5cv5+c//3l+/vOfL7oqS+PJkyen3mT3PEyyiAEAgCUiwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBozE9m/cJSyuNa660x9zaSrCQ57D72a62PTlsWAIAZA1wp5XaS9TH3tpN8VWt9MHytlNIbvjZtWQAABqYOcKWUlSR3x9zrJdmotV4dvl5r3SqlfF1K2a21Hk5bFgCAH8wyB24jyadj7m0m2Rtzr9997SxlAQDoTBXgSilrGYSucT1j6xmEr1H6SW7NWBYAgM60PXDrtdb9E+6vJflqzL1+khszlgUAoDNxgOtWi+6e8vlWzqgsAMCFMdEihm7BwcFJiwq6xQ3J+OHVmcqO+NqNdPPj3n///Tx58mTah5jaXy79VQ7e+dszf55ZPXkx5VqUc2izN3nx4sW5fO8uEm06X9pz/rTp/GnT+WulTSd957+9LNt61Fp30/UE3rhxo968efPMn/OfHv1jVv/jX878eWb16967033BzX84m4pM4cmTJzmP791Fok3nS3vOnzadP206f6206RuHULs93964se5Q79wbhz6nKQsAwKtODHDdUOdqrXXcatFpHJ22MO+yAAAXypuGUDeS/KqUcv3Y9V6SlFJ2us8f1lr3Mthi5NqYx1pN8nTo82nKAgDQOTHAjZv3dnSUVq1189it/Qy2Bxmll2RnxrIAAHRmPsx+jJ0kz45f7IZi15J8NmNZAAA6sxyllQyGOIe3A0mSdHPldrtD6odtJ9ka3oZkmrIAAPxgqh64buj0bgbHYCXJ70op/SQfHQWu7jD6jVLKvQwWIlxL8rjW+tpK1mnKAgAwMFWA64LVJFuKTHxiwzRlAQCYfQgVAIAFEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBozE8mLVhKuZfk3SRr3aWdWuujMWU3kqwkOew+9udRFgCACQNcKWU7yf1a62H3+VqS35VS7tZa74wo+1Wt9cHwtVJKb/jatGUBABh4Y4Dresh2jsJbktRa90spd5I8LqXcPuoxK6X0kmzUWq8OP0atdauU8nUpZXcoBE5cFgCAH0wyB+5arbV//GKtdS+DYc/NocubSfbGPE4/ycaMZQEA6EwS4O6VUh6Oufc0yY2hz9czCF+j9JPcmrEsAACdSQLcuF6yJFlNcjD0+VqSr8aU7efVsDdNWQAAOm+cA1drPaknbC3J7hTPt3JGZQEALoyJtxE5rlvckCTb3edHgeuNCw+mKTvmeTeS5P3338+TJ0+mfYip/eXSX+Xgnb898+eZ1ZMXU34bz6HN3uTFixfn8r27SLTpfGnP+dOm86dN56+VNp0pwHUBbDvJnVELHM5SrXU3Xa/fjRs36s2bN8/8Of/p0T9m9T/+5cyfZ1a/7r073Rfc/IezqcgUnjx5kvP43l0k2nS+tOf8adP506bz10qbznoSw8MkW8Mb7g5t+fHGoc9pygIA8KqpA1y3+e7DridsGkenLcy7LADAhTJVgOuO0/rihPC2l+TamHurGWw7MktZAAA6Ewe4UsrtJIfHw1spZX3o0/0kvTEP0ctg6HWWsgAAdCYKcEchbUzP2/A2IzsZsX9bt+hhLclnM5YFAKAzyVmoaxkce/Vp1wt3ZDWDuWovFyLUWvullN1SynatdWuo7HYGix4OZykLAMAPJtlG5HcZhLTbY+4Ph6+jw+g3uvlyhxnMc3s8vGJ1lrIAAAxMchLD1WkfdJoVqjOsZgUAuNBm3QcOAIAFEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYybZyJcl9/v+V1OV/+e//OsZ1WS039765bk+HwD82OmBAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADTmJ4uuAOfv77/cPd8n/Pzd16+9+EXy+f35PP5vPp7P4wBAI/TAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANCYnyy6AnBqn99fdA3m5zcfL7oGADRAgOPM/b7/1WvX/vTO3+T3//769UX4de/dRVcBAKZiCBUAoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQmJ8sugLAkM/vz/61L35xuq8/C7/5eNE1APhR0gMHANAYPXBceL/vf7XoKpzo1713F10FAJaMHjgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGNs5AtLbtKNhv/0zt/k9/9+/psSn7jR8LId7TWNUUeTORoMWBJ64AAAGiPAAQA0ZuFDqKWUjSQrSQ67j/1a66PF1goAYHktNMCVUraTfFVrfTB8rZTSG74GMKtJ5xAeN2pO4T//5V/nUaXX/PbWL8/kcYEfr4UFuFJKL8lGrfXq8PVa61Yp5etSym6t9XBB1QMmNGtAYgmc9yKTUQtD5skiEy6QRfbAbSbZG3Ovn2QjiV44YGn8/Ze7Z/PAn5+wkhdghEUGuPWcHOBuRYADLoBl78U8casYYCEWuQp1Lcm4V61+khvnWBcAgGYsfBXqCVYWXQEA5tdDeOabTff/+9k99imcaQ/mWc8rvIgmadMlmG+5kABXSjkKZ1MvUui2HdnoPn1RSvlfc6vYeO8l+X/n8DwXiTadP206X9pz/rTp/GnT+ZugTf/HuVQkyX8Zd2OZe+BGqrXuJjmjmcSjlVKe1loN6c6RNp0/bTpf2nP+tOn8adP5a6VNFzIHbmh7EMOkAABTWtajtI5OZgAA4JhFBri9JNfG3FtN8vQc6/Im5zpke0Fo0/nTpvOlPedPm86fNp2/Jtq01FoX88SDY7TWaq23Rtx7lmSnm+8GAMCQRQa4XpJnx4/S6laofp3kqqO0AABet7AAl7zshUutdWvo2k6SLxxmDwAw2kIDXPJyX7ejRQvXkvyh1vpooZXqHKvbSpL+stStNV3P6naSx9pwfkop95K8m8HJJslg6oH2nVH3c7qRwWvRYZJeBifD3DciMD+llMejps/wZqWUhxn8TO7UWvvdaNZ6kuu11s3F1q5t3Xv+8Nz8r5a5M2nhAW5Zdb2Dr3zzRl3jZF2P6mqSP2QQ4O4IGPPR/Ty+DBallLUkv0uyV2u9s9DKNagLbx8Pjwh01x9mEJCvC3GnV0q5neRhrbUsui4t6uaIrx27vCcQn073e/6HY+/5DzPodFjK+fgC3Ajj5ud1975O8gsv5NMrpdQIcHPR/aW4V2vtH7u+nuRxtPPUut7Mu0n+6/DvdxeMnyXZXNYX8lZ0IfmTJLcFuNl0f7g9zqB3OEme1lr3F1il5nVtunK8B7N7v/9oWV9LmzuJ4ZxsZrDNySj9DIZY9MKxSNdGhYla614p5TCDn+GlfNFZYv388KbI2dhI8mmS24uuSMO+qrWOe39iSl2Hzb2M2NZsVCfOMlnWjXwXbT2DF/NR+kl0VbNo97ru/VGeJln6Y2CWTa31Ua111Or3o7b0pnkKXU/mXmzSznLZzGB++7j3/KUlwI22luSrMff68ebI4p0UJlaTHJxXRS6ArSS7Lb7AL5l1Q33zU0pZKaWsdz1IzG49yX7ysk1vt9KmAtxsnOHKQtVab52wUOGop4NTKKX0ul7OR1b3nU43Z9P8wfl4t2vPGxn0tq+UUh538wuZ3lqSg25xzXoGr50rpZSHyx7kzIE7ZuiXQDc/zele2JPBil9mMLSVyLsZrJ7W83YK3ZvggYVfc/XZUHvud5Pwn2X88ZScbDXJ4dDcwv1SylaSZ6WUpV20qAcOfiSG9tq7Y7hvdrXWw1rrg1rrVrelwK1SyuNF16tht5d1FV+Lup/Lw2PX9pKsDv0BxwSGOmzWji8M6V5D+0k+PveKTUiAO2boF0N3NK15mGTLm+XcbSVZPzo5hsl1w1J+Hs9HP4n9H6cw3Is5pkg/S7xiWoCb3tHJDLA0unDx0D5l89e9yO9nMKzKhLrejVW9wefGNjizOcz4RV8HWeI2NQdutL2Mn0uwmsHEUVgK3Qa0Xwhvp1NK+SKDBQtbI24fJFkrpaws63yYJbSR5FellOvHrveSl6e0JIM/PCy6mUB3CsPemJ9RZtPP4H19nKX9fRfgRtvP60eVHOkl2RlzD85VN0R1eDy8lVLWvSlOrustOukv7V4G7by0L+bLZtyRg0er/azsnclKki/G3Otl/FAg432awV5wo/SyxB02hlBH28mIvd66F/m1JJ+de43gmO7YrIzpebPZ9BS6YLY7qmdjKNz5vWfRdkb9vncrfdeic2EWu0l6Y7YMWc8Sr+gX4Ebo5mzsjpi0vJ3BJHF/hc/upK5qJtTtar/Z/f/tof82uiFVi3Cmtz00rDfskyT7eozmZjV5ZQUgkxv1vpQMgtsDve7T697Pt3Is/HavBbvL3KYOsz9BtyT7aNHCtSR/sMJvOl2Y+FUGfx32MmjLpxnMOxCGZ9QdsnzSG+DWuCEsxutCxfC2Ab0Mfu+15Sl1Q6d3M+jVWMlguK+fwWHhXgcmdOxndCWDQLyzzEGjBd37/fB8zWfLPq9YgAMAaIwhVACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI35//VHNvjdOzXdAAAAAElFTkSuQmCC\n",
+ "text/plain": [
+ "<Figure size 735.243x454.405 with 1 Axes>"
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnAAAAGbCAYAAABEeK1/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAcf0lEQVR4nO3dT2jd14Ev8O95buIyxY2s5BFKwC+jQNaJ7IEuZuHy5NXsBrveDWQRaVfo42ERGGZrZPqma6kMmcVsEntXZmWVatdCbTO78lJyOy+EWQyJrIlTiJOG8xb3J+da1pXvlSXde3I/HwiOfr+j3+9YB8lfnb+l1hoAANrx3yZdAQAAxiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDhgqpVS1kopy5OuB8A0KfaBA6ZZKeV+kl6t9fyk65IkpZRrSV5L8mGSF5Ok1ro6ieeVUi4nuVRrXTns+4E2fWfSFQAYppSylGQuyWIpZaHW2ptwfW4nuTkYmEopS6WUu4cJmId5XillIclikqtJLifZGPe9QPsMoQLT7EqS3d6oy5OsSCllLUlqrY8FplrrZpI7pZSbx/28rsdtLclCkrfH+gsA3yqGUIGp1YWc60l2h1Ffm2BdapIrtdZb+9xbTHI3ydla685JPa97xoYhVJg9euCAqdQtXHivCzCbSRa6YDOJuuz2/t0bUmR3aHekxRZH/Txg9ghwwLQ6X2vdDTjr3Z+T6mm6lCTD5uAN9JJdmtDzgBkjwAFTp5Qyl/6qzCTJwDDjjydToyyMUKY3YrnjeB4wYwQ4YBotJ9k7N+xWkrluZepJm5/ycsCMEeCAafTaPsOLkxxGnUsyyuKEuQk9D5gxAhwwVbp9zu7uvd5tr7GTp2wnUkq5VkqpR/DfYHgapSdspNWnx/Q8YMYIcMC0Wdm7N9qAjeSxVZxPqLXeqLWWI/hv3AB11L1let+AoZzEAEybywdsF7I7qX8lT86RO07bGW1Bwaih76ifB8wYAQ6YGt0ChdX9NrcdKPNhkqVSytwheskOa5T3zKcfzCbxPGDGGEIFpsm+JxPssbuY4SQ3ub2TPNreZJi5DN+Y97ifB8wYAQ5oze78uKv73TymRQy7iyr2HfYcKHt7xL/DUT8PmDECHDAVuoUJTz0Qvhs2vZdksVuxuvf+cSxieL/788KQau3uTbc54l/3qJ8HzBgBDpgWV7utQkZxonvCdWHuxgHvu5r+ofJPHI2134KMZ3nePmz2CzNIgAMmrju4/vJT5oQNutP9eWLz4Gqtq8mjuj7S9RwuJlnd+zmllGtJ7pZS1vfeO8zz9pTb/VrZbgRmkFWowMR0q05v5psQ8sdSyttPWYV6O98MPc6VUu6nf27olRF7rA6t1nq+m2O3nv5ZrS92118b8in30l9x+sTGxId53u4q3fTnzu0OHy8NfA22a62XDvWXA5pSaq2TrgMAAGMwhAoA0BgBDgCgMQIcAEBjBDgAgMY0vQr1pZdeqq+++uqxv+dPf/pTvve97x37exiPdplO2mU6aZfppF2m07S0y927dz+ptf73/e41HeBeffXV3Llz5+kFn9HW1lYuXrx47O9hPNplOmmX6aRdppN2mU7T0i6llP837J4hVACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAY5reyBcAOLyHDx9me3s7Dx48yNdffz3p6kyNF154Ib///e+P9JmnTp3KmTNnMj8/n9OnTz/z8wQ4AJhBDx8+zEcffZSzZ8/m1VdfzXPPPZdSyqSrNRUePHiQM2fOHNnzaq356quv8tlnn+Wjjz7KuXPnnjnEGUIFgBm0vb2ds2fP5qWXXsrzzz8vvB2jUkqef/75vPTSSzl79my2t7ef+ZkCHADMoAcPHuT73//+pKsxc77//e/nwYMHz/wcAQ4AZtDXX3+d5557btLVmDnPPffckcw3FOAAYEYZNj15R/U1F+AAABpjFeoI/vTJx/nNP/3vSVfjmf323HKS5KeXXp9wTQCAZ6EHDgCgMXrgAIAn/Pz2B5OuwoFmfTRJDxwAMHNWV1dz/vz5lFKysrKSzc3NR/fefffdXLp0KaWUXLlyJRsbGxOs6f70wAEAM2dtbS0bGxtZXV3N+vr6Y/feeuut/PVf/3XOnz+fX/ziF5mbm3vs/r1797K5uZmFhYVsb29nYWEhS0tLJ1l9AQ4AYFS9Xi/Xr1/PzZs3H127cuVK5ufns7i4eGL1MIQKADCitbW1rKysPHbtnXfeyerq6onWQ4ADABjR+++/n4WFhceuLSwsPDaH7iQIcAAAI+j1etnZ2XkiwO3Okbt3796J1cUcOABgZu3s7OTGjRuPXXv48GE+/vjjfcseZHt7+0jrdhABDgCYWXNzc7l27dpj1x48eJA//OEPU7l9yC5DqAAAY3haT9xJEOAAAEawO9dt71DpbqCbn58/sboIcAAAI1hYWMjc3NwTPXC7gc4+cAAAU2hpaSm9Xu+xa71e78RPYhDgAABGtLa2luvXrz92bX19PWtraydaD6tQAYAn/PTS65OuwrFaXV3N5uZmdnZ2srKykitXrjzqRXv33Xfzy1/+Mkny9ttv5+rVq7l8+XKS/jDq7jmqCwsL6fV6WVlZOdHh00SAAwBm0EE9Zm+99VZ+8pOfDL1/0sOl+zGECgDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBozFgb+ZZSlpO8NnDp01rrjSHl5pLsdH/2aq23DnjmSGUBABgjwJVSbib5Xa11dfBaKWW51roxcG0te4JdKWWtlLKwN+yNUxYAgL6RAlwXtLb3CVVLSd4bKLeQZLnWenawUK11tZRyv5SyUWvdGbcsAADfeOocuC5oXUvyxKFhtdaze4Y7V5JsDnlUL8nyIcsCANAZZRHDSvrz0nojlF1KP3ztp5fk0iHLAgDQGWUIdSnJvSQppcztfjwk0C1mYEh1j729auOUBQBO0q+vT7oGB/vRO5OuwUSN0gO3mGS7lHI5/fC2mWSuW8CwMOb75o6pLADAyFZXV3P+/PmUUrKyspLNzW9mdb377ru5dOlSSim5cuVKNjY2DnjSZIy6CnU+yU6tdfdvd6+UsprkbinlL2utO13vXNLfDuRA45Td53OX0/XOvfzyy9na2hr3EWP786m/yPYLbxz7e47bK1/8MUmytfUfE67J0fj8889PpP0Zj3aZTtplOk2yXV544YU8ePBg6P3nv3x4grUZ35cH1H0Uf//3f59XXnklvV4vP/vZz5Lk0dfj7/7u7/Lmm29mc3Mz//iP/5i5ubl9v1b/8A//kL/927/NG2+MlxG++OKLZ273AwPcQNBaHAhvSZJaa6+U0kvyTpLVJz75mHRblmwkyYULF+rFixeP/Z3/eutfMv9f/3bs7zluH5zrj0r/+OLrE67J0dja2spJtD/j0S7TSbtMp0m2y+9///ucOXNmeIHnT59cZQ7h9EF1H9F3v/vdJHni6/DgwYN873vfe3Rv8H6v18va2lrm5ubyz//8z/mbv/mbg7+OQ9775ptvPlPdDwxwXc9a0s2B20cvyeUkqwNlnzr0OU5ZAIBpsbCwkPX19STJrVuTO3dglDlwO0m2h9zbTjLqPLjd0xaOuiwAwEwZJcD10p8DN8xg0NrM40dtDZpPcueQZQEA6IyyiOG99PeC289CHg9a99JftTqs7PohywIAHLmdnZ3cuPH4QVMPHz7Mxx9/PKEajWaUALeRZPd80r17vy3l8Q1315Pc3fuAbjHEYpL3D1kWAODIzc3N5dq1a49de/DgQf7whz9M5fYhu546hNqdR7qaPT1ipZT1JBuDq1O7gLfRnZ06aC3dQofDlAUA4Bsj7QNXa71RStnpQtuuu92WHnvLrpZSlksp19KfH/daktt7zkwduywAAH2jbuSb/cLaSZcFAGC0VagAAEwRAQ4AoDEjD6ECADPkR+9MugbHanV1NZubm9nZ2cnKykquXLmSpaWlJP3D7H/5y18mSd5+++1cvXo1ly9fTtLfduT69evp9Xrp9XpZXV3N0tJSLl269OjzT4IABwDMnLW1vZtgfOOtt97KT37yk33vzc3NHfi5J8UQKgBAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcAM6rWOukqzJyj+poLcAAwg06dOpWvvvpq0tWYOV999VVOnTr1zM8R4ABgBp05cyafffbZpKsxcz777LOcOXPmmZ8jwAHADJqfn8/9+/fzySef5MsvvzSceoxqrfnyyy/zySef5P79+5mfn3/mZzpKCwBm0OnTp3Pu3Llsb2/n3//93/P1119PukpT44svvsh3v/vdI33mqVOncubMmZw7dy6nT59+5ucJcAAwo06fPp0f/OAH+cEPfjDpqkyVra2tvPnmm5OuxoEMoQIANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADTmO4f9xFLK7VrrpSH3lpPMJdnp/uzVWm89a1kAAA4Z4Eopl5MsDbm3luTTWuuNwWullIXBa+OWBQCgb+wAV0qZS3J1yL2FJMu11rOD12utq6WU+6WUjVrrzrhlAQD4xmHmwC0neW/IvZUkm0Pu9brPPUxZAAA6YwW4Uspi+qFrWM/YUvrhaz+9JJcOWRYAgM64PXBLtdZ7B9xfTPLpkHu9JBcOWRYAgM7IAa5bLbrxjO+bO6ayAAAzY6RFDN2Cg+2DFhV0ixuS4cOrhyq7z+cup5sf9/LLL2dra2vcR4ztz6f+ItsvvHHs7zlur3zxxyTJ1tZ/TLgmR+Pzzz8/kfZnPNplOmmX6aRdplML7TLqKtTL07KtR611I11P4IULF+rFixeP/Z3/eutfMv9f/3bs7zluH5zrrwv58cXXJ1yTo7G1tZWTaH/Go12mk3aZTtplOrXQLk8dQu32fHvqxroDvXNPHfocpywAAI87MMB1Q53ztdZhq0XHsXvawlGXBQCYKU8bQl1O8lellPN7ri8kSSllvfv4Zq11M/0tRl4b8qz5JHcGPh6nLAAAnQMD3LB5b7tHadVaV/bcupf+9iD7WUiyfsiyAAB0Dn2Y/RDrSe7uvdgNxS4mef+QZQEA6BzmKK2kP8Q5uB1IkqSbK7fRHVI/aC3J6uA2JOOUBQDgG2P1wHVDp1fTPwYrSX5VSukleXs3cHWH0S+XUq6lvxDhtSS3a61PrGQdpywAAH1jBbguWI2ypcjIJzaMUxYAgMMPoQIAMCECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMd8ZtWAp5VqSF5MsdpfWa623hpRdTjKXZKf7s3cUZQEAGDHAlVLWklyvte50Hy8m+VUp5Wqt9co+ZT+ttd4YvFZKWRi8Nm5ZAAD6nhrguh6y9d3wliS11nullCtJbpdSLu/2mJVSFpIs11rPDj6j1rpaSrlfStkYCIEjlwUA4BujzIF7rdba23ux1rqZ/rDnysDllSSbQ57TS7J8yLIAAHRGCXDXSik3h9y7k+TCwMdL6Yev/fSSXDpkWQAAOqMEuGG9ZEkyn2R74OPFJJ8OKdvL42FvnLIAAHSeOgeu1npQT9hiko0x3jd3TGUBAGbGyNuI7NUtbkiSte7j3cD11IUH45Qd8t7lJHn55ZeztbU17iPG9udTf5HtF9449vcct1e++GOSZGvrPyZck6Px+eefn0j7Mx7tMp20y3TSLtOphXY5VIDrAthakiv7LXA4TrXWjXS9fhcuXKgXL1489nf+661/yfx//duxv+e4fXCun7l/fPH1CdfkaGxtbeUk2p/xaJfppF2mk3aZTi20y2FPYriZZHVww92BLT+eOvQ5TlkAAB43doDrNt+92fWEjWP3tIWjLgsAMFPGCnDdcVofHhDeNpO8NuTefPrbjhymLAAAnZEDXCnlcpKdveGtlLI08OG9JAtDHrGQ/tDrYcoCANAZKcDthrQhPW+D24ysZ5/927pFD4tJ3j9kWQAAOqOchbqY/rFX73W9cLvm05+r9mghQq21V0rZKKWs1VpXB8qupb/oYecwZQEA+MYo24j8Kv2QdnnI/cHwtXsY/XI3X24n/XlutwdXrB6mLAAAfaOcxHB23IeOs0L1EKtZAQBm2mH3gQMAYEIEOACAxghwAACNOfRh9rTr57c/mHQVjsQrXzzMz29/kJ9e+nac7QoAo9IDBwDQGD1wM+SHH317Fvz+9tzypKsAABOjBw4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACN+c6kKwCH8cOPNrL9whv54X/eTn794qSrc3g/emfSNQCgQXrgAAAaI8ABADTGECrN+03v00lX4dB+++cPHvv4p5den1BNAGiJHjgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANMZh9jBFfn77g6cXmnKvfPHw0d/jp5den3BtAL6dBDiYoB9+tDHpKhyZ355bnnQVAGaGIVQAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxjiJATgSu6dKbL/wRn74n7f7F3/94gRrdEg/emfSNQB4KgEOODa/6X066SqM7bd/Hn4erbNdgWlhCBUAoDECHABAYyY+hFpKWU4yl2Sn+7NXa7012VoBAEyviQa4Uspakk9rrTcGr5VSFgavAZyU3cUY+2ptUYYFGfCtNbEAV0pZSLJcaz07eL3WulpKuV9K2ai17kyoegBPaG1RxuCCjFe+eJif3358gYZFGdCuSfbArSTZHHKvl2Q5iV44gGOyN9C1TiBllkwywC3l4AB3KQIcwKENDgc/tj/ft1ULQ9yGtTkikwxwi0neG3JvtwcOAEbSwhD33n0G9xvaThroTfz19bGKT2vb/Pbc/lFjWLvsmob2mfgq1APMTboCAHCU9i6SGdYz+pt/Oqka0aqJBLhSym44G3uRQrftyG5k/ryU8n+PrGLDvZTkkxN4D+PRLtNJu0wn7TKdtMtE/Z9hNw5sl/91LHXZ1/8YdmOae+D2VWvdSHLAOv+jV0q5U2u9cJLv5Om0y3TSLtNJu0wn7TKdWmiXiZzEMLA9iGFSAIAxTetRWrsnMwAAsMckA9xmkteG3JtPcucE6/I0Jzpky8i0y3TSLtNJu0wn7TKdpr5dSq11Mi/uH6O1WGu9tM+9u0nWu/luAAAMmGSAW0hyd+9RWt0K1ftJzjpKCwDgSRMLcMmjXrjUWlcHrq0n+dBh9gAA+5togEse7eu2u2jhtSS/q7XemmilOnvqNpekNy11m1VdD+1aktvaYrqUUq4leTH9U1aS/jQIbTRB3ffLcvo/W3eSLKR/0s11IxzTpZRye78pRZycUsrN9L8/1mutvW6kcCnJ+VrrymRr96SJB7hp1fUOfjrYE7jfNU5G1zM7n+R36Qe4K8LB9Oi+Nx6FglLKYpJfJdmstV6ZaOVmVBfe3hkc4eiu30w/ZJ8X4qZDKeVykpu11jLpusyybv794p7Lm9MarAW4fQybn9fdu5/kL/3gm5xSSo0ANzW6nurNWmtvz/WlJLejrSai6xG9muR/Dv686sL13SQrFopNXhe0f5HksgA3Wd0vorfT76lOkju11nsTrNKBpnUfuElbSX+bk/308s1RXkDy2t7wliS11s30h+2mbuhhRvTyzT9ETK/lJO9NuhIk6Y+wbdZaN7r/pja8JQLcMEvp//DbTy/JVHanwoRc64bl9nMnyVQfR/NtVWu9VWvdbzX/bnsM+yWVE9L1hu7+ogNjEeD2t5jk0yH3evEPEgw6KAjMJ9k+qYowktUkG/v1mnLilqa9l2cWlVLmSilL3XSqqSXAHY4zXKFTa710wEKF3R4GJqyUstD1lN6axhV1s6abO2oO4nR5sWuXC+mPHsyVUm538xSnzncmXYFpM9BQurThGXQ/CJP+qmEmZGArkRfTX8Wt523Cup6dbYvhptL7A+1yr1vYcDfDj/6cGD1wwJEb2K/viqG6yaq17tRab9RaV7stkC6VUm5Pul4z7rKV2dOn+x7Z2XNtM8n8wC+kU0OA22Og8aayyxQacTPJqn+kptJqkqXdk3A4Wd2eb74v2tJLMnX7WQpw49s9mQHYRxcMbtpjbDp1v6Tei+2QTlzXMz2vV7o5U7kljzlw+9vM8PHu+fQnNwJ7dJvHfii8TV4p5cP0Fyys7nN7O8liKWXOPKwTtZzkr0op5/dcX0genTiT9H8BsvjnBHWnMGwO+X6ZSgLc/u7lyeM0di0kWR9yD2ZWNzS0sze8lVKW/GN0srqenoN6DBbSbyvh7QQNO4ax+95Zsjp4ouaSfDjk3kL6uWCqGELd33r22eut+6G4mOT9E68RTLHu2KwM6Xmz8fUJ64LZxn69CQPhzs8x+Mb6fj+/uhXDi5nCjhsBbh/d/ISNfSb5rqU/MdtvrZM3P+kK0NftJr/S/f/lgf+WuyFVC4ImY21gSG7QL5Lc09szVeaTx7ax4uTt929+0g9uN6ZxFMFh9gfolg3vLlp4LcnvrKqbjC4I/FX6vwktpN8md9KfXCpUT1Ap5X4ODmmrw4aOOF5dIHhn4NJC+j/HtMcU6IZOr6Z/fONc+sN0vSRv+5l28vZ8v8ylH6zXpzG8JQIcAEBzDKECADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBozP8H0EdEgT/aMNQAAAAASUVORK5CYII=\n",
+ "text/plain": [
+ "<Figure size 735.243x454.405 with 1 Axes>"
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "for A in [2,1,0.1,0.01]:\n",
+ " stat_h0,stat_h1 = dist_unknownAsigma_twopulse(A=A,sigma=5,x_1=30,x_2=60)\n",
+ " plt.figure()\n",
+ " plt.hist(stat_h0,alpha=0.5,label='H0')\n",
+ " plt.hist(stat_h1,alpha=0.5,label='H1')\n",
+ " plt.legend()\n",
+ " plt.title('A='+str(A))"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 21,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnkAAAGbCAYAAAC4bYa1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAd1klEQVR4nO3dT2+cx4En4F+t/w0wq4SivfBhsLM2BeQDUBpg96ZgyG8gxfsFTJ7nIkL7BQTqMmcyX2Bt+ZZTIAbhbQ6xeB3sAO7BBshgMbBkwlKAkY2k9tBvKy2qm+xuNrup4vMAgcT3LfZbKZaln6reqiq11gAA0Jb/tOwKAAAwf0IeAECDhDwAgAYJeQAADRLyAAAa9O6yKwAwi1LKRpK7Sb7pLn2Y5HGSr5Psdtc/TPK7WutXS6kkwBIVW6gAb5tSyuMkvVrr9tC19SS/6b78+1rrUSnlSZK1Wuv1ZdQTYJmM5AFvlVLKXpLVWuvmiVsbSVaSHNRaj7prX6Q/mnfRdVrp6tS76GeNef69JDfyl9HL1Fp3llEX4PIwkge8VUopNcl2rXV/zL2dWuvDBddpI8nKMqaFu1HNR8Pt0dVnt9Z6c9H1AS4PCy+At0YpZa377WkjZseLqMsJa2cXmb9Sym6SnAy8tdaDJF+XUh4to17A5SDkAW+jZzPeuyh3l/DMJLmXZG/Mvb0kd7qpZOAKEvIAzqEbXdxYwnPvdL89GlNkMNq5tYDqAJeQkAcwo26U7PGSHr+ZJOMWe9Raj4fLAVePhRfAUg3eK0vyNP0VosdJHgyFlMGo1XaS1STr6Y9SnQw3o+4dDVaZllJOjmitJDlIfxHH9vCNLrztdnUZ1OvJicUNW+lP046r087QKt+56xZcbNRayyllvkmSWuuNi6oHcHkJecBSdNOcj9IPQwdD19eT/DLJ5ydDUvc93yTZHP6eofsjV952267sjfi8vfSD0o2ha+tdvTaHR8lG7c03VP5JkruLXF3b7QG4PkHIW7VPIFxNpmuBZXmcfvB6Lax1QWwnyW/msWigC4a3Ro2qnQxsnd+kv/3IyZHC7SRb3fYkl8FKJltJbOEFXFFCHrBw3ea9a6P2uktebQHyLP0p0/NaS386dZzhUcTd9Pe7e6NeXegbBNCplFLulVLqHP43HNhWJ3j0MraTAS4JIQ9Yhu0MhasxDjKflaFfJ/3pzTGjcMNBciun78HXS3Jr2grUWh/WWssc/jdtaDOKB1eYY82AZVjL+K0/Bo6T/nTreY4Lq7Uel1I2058eflxKSfoB8nGS/cFnd6NkK0mejVikMTAocxk8y2SbMBvNgytKyAMWqluokEy+afFaTh9dO1Ot9aCUcj39kbrN9Pe120hyv5Ryt5seHgSmg3HTyEnGXV+GScLbapazOTRwCZiuBRZtENgmeadsuPxMSikrpZSVWutxN2262a1I3ew++3F3fzCyOGm9znrunbNLnctgGvq0kcWVnD1iCjRKyAMWqnuv7DhnTzWudeXPFfLSf4fuFyPqcVBrvdnVZTA925ugXpP6u8FvLmjhxZPu15H1HSq7rM2agSUT8oBl2E+yfsYo1J0k89p37rRTH74c+v1XOX0l7qhNlc+cNr2ghReDeo9bCDJYZHLWAhegUUIesHDdKRTDI2ivGZrq/HyGjx8VHO+cEihX001pdvXqDZ3CcbJebyy8GBppnNcI4ES6wPcw/ZXKo3yWoYUlwNUj5AHLcjPJ9sl317qFGbvpnzhxcpRs5cSvw983uPbhiGcdZ8See92zVk5syLyZfig8Wa+VJPdrrQ9HfP5+hsJWt1XLFyPKzdW4I9u6uq9nhj39gHY41gxYqm7U7GleD2cnz65dSf8kirX85aSHZ+mOHuuO+BrcS/rv1h3VWu92gWtwTu397v7TwWePCW2Deg2v7H06rmxXfq8r/zjJ8SkrdOeu21z6RvpHvn2Y/CUAAleXkAcA0CDTtQAADRLyAAAaJOQBADRIyAMAaNBbfXbtRx99VD/55JMLf84f//jH/PVf//WFP+eq086Lo60XQzsvhnZeHG29GIN2fvLkybe11v8y6+e81SHvk08+yddff33hzzk8PMzt27cv/DlXnXZeHG29GNp5MbTz4mjrxRi0cynl/57nc0zXAgA0SMgDAGiQkAcA0CAhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAa9FZvhgwAzO7ly5d59uxZnj9/nj/96U9nlv/pT3+af/7nf15Azdrzzjvv5Nq1a1ldXc0HH3ywkGcKeQBwBb18+TK///3vc/369XzyySd57733Uko59XueP3+ea9euLaiG7ai15scff8z333+f3//+9/nbv/3bhQQ907UAcAU9e/Ys169fz0cffZT333//zIDH7Eopef/99/PRRx/l+vXrefbs2UKeK+QBwBX0/Pnz/OQnP1l2Na6cn/zkJ3n+/PlCniXkAcAV9Kc//Snvvffesqtx5bz33nsTvf84D0IeAFxRpmgXb5FtLuQBADTI6toJ/Pvzl/nHx/9yZrl/2PzZAmoDAHA2I3kAAA0ykgcAvGHUDNYPP7zM++8vZiPfs5g9O5uRPADgytnZ2cnNmzdTSsn29nYODg5e3dvf38/m5mZKKbl792729/eXWNPZGckDAK6c3d3d7O/vZ2dnJ3t7e6/d29rayq1bt3Lz5s388pe/zMrKymv3j46OcnBwkLW1tTx79ixra2vZ2NhYZPUnIuQBAEyo1+vlwYMHefTo0atrd+/ezerqatbX15dYszeZrgUAmNDu7m62t7dfu3b//v3s7OwsqUbjCXkAABP68ssvs7a29tq1tbW1197puyyEPACACfR6vRwfH78R8gbv7B0dHS2jWmN5Jw8AuLKOj4/z8OHDN65/8803I8ue5tmzZ3Or1zwIeQDAlbWyspJ79+69cf3o6Oit3TplwHQtAMAUzhrRuyyEPACACQzevTs5LTsIfaurqwuv02mEPACACaytrWVlZeWNkbxB6LNPHgDAW2pjYyO9Xu+1a71e71KeeCHkAQBMaHd3Nw8ePHjt2t7eXnZ3d5dUo/GsrgUA3vAPmz9749rz589z7dq1JdRm/nZ2dnJwcJDj4+Nsb2/n7t27r0bj9vf3Xx1b9vnnn+ezzz7LnTt3kvSnbAfn3q6traXX62V7e/vSTdUmQh4AcAWdNvK2tbWVra2tsfcv49TsKKZrAQAaJOQBADRIyAMAaJCQBwDQICEPAKBBQh4AQIOEPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgIQ8AoEHvLrsCAMAl9NsHb1x6/4eXyfsfLKEyI/z8/rJrcOkZyQMArpydnZ3cvHkzpZRsb2/n4ODg1b39/f1sbm6mlJK7d+9mf39/iTWdnZAHAFw5u7u72d7ezsrKSvb29rKxsfHq3tbWVnZ3d5Mkv/zlL7O1tTXyM3Z2dnJ0dLSQ+s7CdC0AwIR6vV52d3ezsrLyasTvshLyAAAmtLa2lr29vSTJV199teTanM50LQBAg4Q8AIAGma4FAK6s4+PjPHz48I3r33zzzRJqM19CHgBwZa2srOTevXtvXD86Onprt04ZMF0LANAgIQ8AoEFCHgBAg4Q8AIAGCXkAAA2yuhYAeNPP779x6Yfnz/PBtWtLqMz87ezs5ODgIMfHx9ne3s7du3dfnV+7v7+fR48eJUk+//zzfPbZZ7lz506S/pYrDx48SK/XS6/Xy87OTjY2NrK5ufna+beXgZAHAFw5u7u7Y+9tbW1la2tr5L2VlZVTv/cyMV0LANAgIQ8AoEFTTdeWUraS3Bi69LTW+sZZIF25lSTH3a+9WutXp3zmRGUBAJjMxCGvlPIoye9qrTvD10opW7XW/aFruzkR/kopu6WUtZOBcJqyAABMbqKQ14WxZyOC10aSL4bKrSXZqrVeHy5Ua90ppXxXStmvtR5PWxYAgOmc+U5eF8buJXljKUmt9fqJqdXtJAdjPqqXZGvGsgAATGGShRfb6b8n15ug7Eb6AW2UXpLNGcsCAHNWa112Fa6cRbb5JCFvI8lRkpRSVkopd7rRvVHWkzwdc6+X5NaMZQGAOXrnnXfy448/LrsaV86PP/6Yd955ZyHPmiTkrSd5Vkq5k37gO0iy0i26GBf2xlm5oLIAwBSuXbuW77//ftnVuHK+//77XFvQqSGT7pO3muS41vpVrfW41nqUZCfJk1LKStIf5evKnrlYYpqyAMD8ra6u5rvvvsu3336bH374wdTtBaq15ocffsi3336b7777Lqurqwt57qmra4fC2Hqt9bVFErXWXimll+R++oFvIbp99baS5OOPP87h4eGFP/O9P7/M3/zHv55Z7vDw3y68Li178eLFQn6eaOtF0c6LoZ1nV0rJ06dP8/7776eUcmb5WutE5XjTIOi9fPkyf/jDH04tO68+fWrIq7Uedz/MozFFeknuJNkZKnvmNOs0ZUd8736S/SS5detWvX379rQfMbUvf/Xr/OGvPj2z3C9u/+zC69Kyw8PDLOLnibZeFO28GNp5cbT1YsyrnSeZrj1O8mzMvWdJJn0vb3CqxbzLAgBwwiQhr5f+O3njDIexg7x+7Nmw1SRfz1gWAIApTBLyvkh/he0oa3k9jB1l/MjeWpJHM5YFAGAKk4S8/SRrY7ZL2cjrJ2HsZcT+dt0CjvUkX85YFgCAKZwZ8rrzY3fSD2WvlFL2kuwPr7rtTsXY7866HbabbnHGLGUBAJjOqatrB2qtD0spx12wG3jSrXQ9WXanlLJVSrmX/vt6N5I8PnHG7dRlAQCY3EQhL3m1dclSywIAMJlJT7wAAOAtIuQBADRIyAMAaJCQBwDQICEPAKBBQh4AQIOEPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgIQ8AoEFCHgBAg4Q8AIAGCXkAAA0S8gAAGiTkAQA0SMgDAGiQkAcA0CAhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAaJOQBADRIyAMAaJCQBwDQICEPAKBBQh4AQIOEPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgIQ8AoEFCHgBAg4Q8AIAGCXkAAA0S8gAAGiTkAQA0SMgDAGiQkAcA0CAhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAaJOQBADRIyAMAaJCQBwDQICEPAKBBQh4AQIOEPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgIQ8AoEFCHgBAg4Q8AIAGCXkAAA0S8gAAGiTkAQA0SMgDAGiQkAcA0CAhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAa9O6s31hKeVxr3RxzbyvJSpLj7tderfWr85YFAGAyM4W8UsqdJBtj7u0meVprfTh8rZSyNnxt2rIAAExu6pBXSllJ8tmYe2tJtmqt14ev11p3SinflVL2a63H05YFAGA6s7yTt5XkizH3tpMcjLnX6753lrIAAExhqpBXSllPP5iNG2HbSD+gjdJLsjljWQAApjDtSN5GrfXolPvrSZ6OuddLcmvGsgAATGHikNetgt0/5/NWLqgsAABDJgp53SKJZ6cthOgWZCTjp3JnKgsAwPQmXV1757JsadKNKG4lyccff5zDw8MLf+Z7f36Zv/mPfz2z3OHhv114XVr24sWLhfw80daLop0XQzsvjrZejHm185khr9sT78zNiWutx6WUZIJp1mnKjvje/XTTxrdu3aq3b9+e9iOm9uWvfp0//NWnZ5b7xe2fXXhdWnZ4eJhF/DzR1ouinRdDOy+Otl6MebXzqdO13bTqaq113CrYaQxOtZh3WQAATjhrJG8ryd+VUm6euL6WJKWUve7rR7XWg/S3V7kx5rNWk3w99PU0ZQEAmMKpIW/ce3iDY81qrdsnbh2lvzXKKGtJ9mYsCwDAFGY6u/YUe0menLzYTfuuJ/lyxrIAAExhlmPNkv506vBWKEmS7t29/VLK7onyu0l2hrdgmaYsAADTmWokr5um/Sz9I8mS5DellF6SzwehrNa6U0rZKqXcS3/xxI0kj2utb6zQnaYsAACTmyrkdeFrku1UJj4ZY5qyAABMZtbpWgAALjEhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAaJOQBADRIyAMAaJCQBwDQICEPAKBBQh4AQIOEPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgIQ8AoEFCHgBAg4Q8AIAGCXkAAA0S8gAAGiTkAQA0SMgDAGiQkAcA0CAhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAaJOQBADRIyAMAaJCQBwDQICEPAKBBQh4AQIOEPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgIQ8AoEHvLrsCsDS/fbDsGszHz+8vuwYAXEJG8gAAGiTkAQA0SMgDAGiQkAcA0CAhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAaJOQBADRIyAMAaJCQBwDQICEPAKBBQh4AQIOEPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgIQ8AoEFCHgBAg4Q8AIAGCXkAAA0S8gAAGiTkAQA0SMgDAGiQkAcA0CAhDwCgQUIeAECD3p20YCnlXpIPk6x3l/ZqrV+NKbuVZCXJcfdrbx5lAQCYzEQhr5Sym+RBrfW4+3o9yW9KKZ/VWu+OKPu01vpw+FopZW342rRlAQCY3Jkhrxtp2xsEvCSptR6VUu4meVxKuTMYeSulrCXZqrVeH/6MWutOKeW7Usr+UFCcuCwAANOZ5J28G7XW3smLtdaD9KdYt4cubyc5GPM5vSRbM5YFAGAKk0zX3uumT++OuPd1kltDX2/k9OC2meThDGWBcX77YPrvefHpbN93kX5+f9k1AGjKJCN544JYkqwmeTb09XqSp2PK9vJ6IJymLAAAUzhzJK/WunnK7fUk+1M8b+WCygIAMGTmffK6BRlJstt9PQhlZy6WmKYsAADTm3ifvGFdSNtNcnfUooyL1IXLrST5+OOPc3h4eOHPfO/PL/M3//GvZ5Y7PPy3C6/LpfD8/13Ix7748wc5/NX/vpDPHu3TBT7rcnnx5w9y+OKS/f9fwH/Li/bixYuF/Bl11WnnxdHWizGvdp4p5CV5lGRneNPiWutxKSWZYJp1mrIjvnc/3RTxrVu36u3bt6f9iKl9+atf5w9/dfZfiL+4/bMLr8ulcEEv7B+++DS3//PZYZrzu5Rtfft/LrsGc3d4eJhF/Bl11WnnxdHWizGvdp56urbbwPhRF7amMTjVYt5lAQA4YaqQ1x1t9s0pAe8gyY0x91bT33JllrIAAExh4pBXSrmT5PhkwCulbAx9eZRkbcxHrKU/zTtLWQAApjBRyBsEuTEjeMNbrOxlxP523UKN9SRfzlgWAIApTHJ27Xr6R5B90Y3mDaym/+7cq8UTtdZeKWW/lLJba90ZKrub/kKN41nKAgAwnUlW1/4m/SB3Z8z94YCWWutOKWWre3/vOP337h4Pr8SdpSwAAJOb5MSL69N+6DQrb2dYpQsAwBlmPvECAIDLS8gDAGiQkAcA0CAhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAaJOQBADRIyAMAaJCQBwDQICEPAKBBQh4AQIPeXXYFWvKPj/9lonL/sPmzC64JAHDVGckDAGiQkAcA0CAhDwCgQd7JAy6H3z5Ydg3m4+f3l10DgCRG8gAAmiTkAQA0SMgDAGiQkAcA0CAhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAaJOQBADRIyAMAaJCQBwDQICEPAKBBQh4AQIOEPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgIQ8AoEFCHgBAg4Q8AIAGCXkAAA0S8gAAGiTkAQA0SMgDAGiQkAcA0KB3l12Bq+S//36//5vffrjcigAAzTOSBwDQICEPAKBBQh4AQIOEPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgJ14swT/1nk5c9n+sOR0DAJiekTwAgAYZyQOYp98++MvvX3z6+tdvk5/fX3YNgHMykgcA0CAhDwCgQUIeAECDhDwAgAYJeQAADRLyAAAaZAuVS27SjZNtmgwADDOSBwDQICEPAKBBQh4AQIO8kwfAm96m49hOOz7O8WxcYUbyAAAaJOQBADTIdG0jbLUCAAwzkgcA0KClj+SVUraSrCQ57n7t1Vq/Wm6tAADebksNeaWU3SRPa60Ph6+VUtaGrwEAMJ2lhbxSylqSrVrr9eHrtdadUsp3pZT9WuvxkqoHQAvepq1gTmMrGGawzJG87SQHY+71kmwlMZq3JBZyAMDbbZkhbyOnh7zNCHlzN2l4AwDebssMeetJvhhzbzCSBwBclmnn004XuUrekunzpa+uPcXKsivAfJ01ivjHn/7X/NO/Tz7SaKoYAMZbSsgrpQwC3NQLK7otVwajfC9KKf9nbhUb76Mk3y7gOVeddl4cbb0Y2nkxtPPiaOskyf+66AcM2vm/nedDLvNI3ki11v0k+4t8Zinl61rrrUU+8yrSzoujrRdDOy+Gdl4cbb0Y82rnpZx4MbQ1iilZAIALcFmPNRucgAEAwAyWGfIOktwYc281ydcLrMtZFjo9fIVp58XR1ouhnRdDOy+Otl6MubRzqbXO43Omf3D/SLP1WuvmiHtPkux1798BADClZYa8tSRPTh5r1q28/S7JdceaAQDMZmkhL3k1mpda687Qtb0k39RanXYBADCjpYa85NW+d4OFFjeS/K7W+tVSK9U5UbeVJL3LUrdWlFIepX/CyV6ttdeN8G4kuVlr3V5u7d5O3Wj4bpLHp/VX/ft8Jmln/Xs+Sin3knyY/klJSb89x7W5fn0Ok7a1vj277s+OrfQzz3GStfTb8sGoGczz9Omlh7zLqhtlfDo8ojjqGufTvX+5fuLywah3NTldNwq+muR36YePu6f8Rah/z2jKdta/z6nrl6/+8iulrCf5TfrteHdEWf16RlO2tb49gy7g3R+eweyuP0q/PW8OB73z9mkhb4Rx7wt2975L8qn3Beej66yP0/+XTJJ8XWs9WmKVmlBKqRkTPvTv+Tmtnbv7+vc5dCMYB7XW3onrG+m366u216/PZ5q27q7r2zPoRko/S/L3J8LcepInSbYHi07n0affuhMvFmQ7/S1eRumlP8zqX4Xz8bTWOq6tuRj69+Lo3+dzY9QuC7XWg1LKcfp9eRA89OvzmaatE317Vr38JRif5dx9+rJuhrxsG+k34Ci9JIajeZvp37wt7nXTWKN8nWT42Cf9+nymaWtmVGv9qtY6aveQQfsOh7pz92khb7T1JE/H3OtFZ5+7UspKKWWjG57mYunfC6Z/z+y0kaLVJM+Gvtavz2eatn5F356bnST7J6bLz92nhbzZOHN3fj7s3gW5lf6/FldKKY+7l1NZDm0/P/r3OdRaN0++8D9kPacHk5O0+SlmaGt9ew5KKWvdCOpXM6xKPrOtvZN3wlAH9YLu4nw5NHR91L3Q+yTjj71jRvr3Uujfc9aFi6S/ulm/vkAn2/oEfXtGQ9uofJj+Sv2TC17m0qeN5LFUtdadk+8mdC/zrg794QJvJf17/ob2J7x7ciUo83VaW+vb51NrPa61Puza8WGSzVLK43k/R8g7YajTGnJerl6ScVMHzEj/vjT079k9SrIzvJ2Hfn1h3mjrCejbs9lJsjF0Ethc+rSQN73BrtNcrGmWmTM/+vdi6N8z6P4CfDRqq48z6NdTOkdb69sz6ELdUfpTuJOYqE8LeaMdZPw7Bavpv2TKOZVSngz+1cJC6d8LoH/PV7eJ7DenhA79ek7Oamt9e3allG9Oabtn6S9gGYzenbtPC3mjHWX8v0TW0h/C5vxWknwz5t5a+j8H5k//Xgz9e05KKXeSHJ8MHd1pDAP69RxM2Nb69gy68HbaKOda+m3/ajHLKeUn6tNC3mh7GbH/TPcDWk/y5cJr1Ka9Uf9S7PZbWk//58D86d+LoX/PwSBcjBlVGt4MVr8+p2naWt+eXhfe9k+eW5u8FgCH++m5+7Sza8cYevlxZ+jaXvpD2I7GmYNTDmp+nORo1H8ITKY7U3X7lOkW/XsOTmtn/fv8uvM87yf54sSt1fRHk24M7y2mX89umrbWt2fXBeGdk3vidXvlrdVab564fq4+LeSdolsGPni58UaS3025yogzDP6w6L5cSf8PlD1nIk6ve4/m79L/F95a+v326/RfhH5juwP9ezbTtLP+fT7dIeynrS7cOfkXnX49m2nbWt+e3Ym2S/p/jvxuXGg7T58W8gAAGuSdPACABgl5AAANEvIAABok5AEANEjIAwBokJAHANAgIQ8AoEFCHgBAg4Q8AIAGCXkAAA36/8geJe1uzsM3AAAAAElFTkSuQmCC\n",
+ "text/plain": [
+ "<Figure size 735.243x454.405 with 1 Axes>"
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnAAAAGbCAYAAABEeK1/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAdkklEQVR4nO3dQW8cx90n4H9FlmXEK3tEe2EYwXrtEeC7KQXY9yZjyW8gxvsFTJ5z0UD7BQTqkjNp7Aew5FtOgWiEtxwi8WpsAE2yBpLDC4siLDmwZNi1h+mRRsOZ4Qw5M83iPA8gUOwuTle6WfFPVV1VKeccAACU41d1VwAAgMkIcAAAhRHgAAAKI8ABABRGgAMAKMxrdVcA4DhSSisRsRYRD6tD70TEvYi4HxGb1fF3IuKvOeevaqkkwIwky4gApUkp3YuIds55o+fYckR8XX37P3POeymlBxHRzDlfqqOeALOiBw4oSkppKyKWcs6rfadWIqIRETs5573q2JfR6YWbdZ0aVZ3as77WkGtvRkQzIg6icw8OIqJVR32A+dADBxQlpZQjYiPnvD3kXCvnfHvOdVqJiMa8h2qrXseN6PxvPug5fiM6oW7N8DGcTSYxAMVIKTWrv47qWToYcW5WmkcXmYnNnPNGb3iLiKgCbDsi7lY9dMAZI8ABJdo/5rlZWZv3Bavet5XqPb9BusPIK3OqEjBHAhzACVS9gnWEpIPqz/0xygFnjEkMAMdUDU/eq+Pa1QSFUbNrl6uvRwU8oEAmMQC1SiltVn99FBGXo9NjdKvvpfzr0XlZfyk6waQdh9+DG3RuL+fcqj5jva98IyJ2ojMhYqP3RM/MzoOeej3onThRfd7aiDq1embDzlU1qeJe1DChA5gPAQ6oRTX0eDc6IWOn5/hyRHwREZ/3B6DqZx5GxGrvz/ScHzhDtVp6ZGvA521FxErO+XLf9e9W12j3HD+09lxP+QdxSmZ8VuHzQUR81Q2vwNljCBWoy73ozKJ8JYhVC/C2IuLrlNJH/TMsJ1WFvqv9wau61kZK6WHf4a9j8BpqGxHxMKV0d1B4rFPV47YaL4dN1+rq/QPmwyQGYO6qdcqag9Zyi4ioAtJ+dIYxT6oZL4PNIL29f5vRWc/tUL2qQLcXERP3aqWUbqSU8hT+DFwSJOe8k3NuVYsbb0bEZs/QNHAGCXBAHTaiJzgNsRMR/e+tHcf9iIiU0oOqp6pfb9BZj9FrzLUj4uqkFcg53845pyn8ObI3sgpzqxFxfcQSI0DhDKECdWjGy3XKhjmI6AyBnmRLqJzzQUppNTpDtvdSShGdcHgvIra7n131bjUiYn/AhIeubpkStKKzkO/WoOFjoGwCHDBX1Uv/EeMvuNuM0b1iR8o576SULkWnh201Ouu2rUTEzZTSWjVk291NYWfY0G5EDDt+GnUD8u+i0+MJnCGGUIF564axpQnLH0tKqZFSauScD6qhzNWcc4pOkGtHp1eu0fPS/7j1Ouq616fxOSM+/15K6eGI9+K6961hOy04ewQ4YK6q97gO4uj9Q5tV+RMFuOi8s/a7AfXYyTlfqerSHTJtj1Gvcf22+5cZTWJYqeo6sL49+8bGOO/OAWUR4IA6bEfE8hE9Q9cjYlrrqq2OOHen5+9fxegZq4MWBB5nYsEsJjHsRGetumHvEjZ7ygFnjAAHzF21wGxvz9creoYfPz/Gxw8KhddHhMWlqN4Xq+rVHrYER89Ehxd6egin1XM3rq0YHTZbfV+BM0SAA+pyJSI2+t8VqyY5bEand6m/d6vR97X357rH3hlwrYMYsKZcda1G38K83SU4+uvViIibQ7am2o6eiQLVciVfDig3Nd1dH1JKW/3htAqgK2FBXzizbKUF1KoKG4/i1eDVvxdqIzo7JDSjE94OojOLdTXn3K7WO+uei+i8y7aXc16rwlR339Ob1flH3c8etldoVa/eGbCPRu0rWm3L1YzO8iQHI2ayTlX1rlsrXk6+aESnzptTeH8QOKUEOACAwhhCBQAojAAHAFAYAQ4AoDACHABAYYreC/Xdd9/NH3744cyv88MPP8Sbb7458+swmPtfP8+gXu5/vdz/+i3qM3jw4MF3Oef/Ouhc0QHuww8/jPv378/8Oru7u3Ht2rWZX4fB3P/6eQb1cv/r5f7Xb1GfQUrp/w07ZwgVAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFCYohfyBQCO79mzZ7G/vx9PnjyJn3/+ue7qDPX222/HN998U3c1TuTcuXNx8eLFWFpaigsXLpz48wQ4AFhAz549i2+//TYuXboUH374YZw/fz5SSnVXa6AnT57ExYsX667GseWc46efforvv/8+vv322/jggw9OHOIMoQLAAtrf349Lly7Fu+++G6+//vqpDW9nQUopXn/99Xj33Xfj0qVLsb+/f+LPFOAAYAE9efIk3nrrrbqrsXDeeuutePLkyYk/R4ADgAX0888/x/nz5+uuxsI5f/78VN43FOAAYEEZNp2/ad1zAQ4AoDBmoY7hP588iz/c+9uR5X6/+vEcagMALDo9cAAAhdEDBwAcMs7I07w8f/4sXn/91XXTFn3USw8cALBwWq1WXLlyJVJKsbGxETs7Oy/ObW9vx+rqaqSUYm1tLba3t2us6WB64ACAhbO5uRnb29vRarVia2vrlXPr6+tx9erVuHLlSnzxxRfRaDReOb+3txc7OzvRbDZjf38/ms1mrKyszLP6AhwAwLja7XbcunUr7t69++LY2tpaLC0txfLy8tzqYQgVAGBMm5ubsbGx8cqxmzdvRqvVmms9BDgAgDHduXMnms3mK8eazeYr79DNgwAHADCGdrsdBwcHhwJc9x25vb29udXFO3AAwMI6ODiI27dvHzr+8OHDgWVH2d/fn1q9jiLAAQALq9FoxI0bNw4d39vbO5XLh3QZQgUAmMBRPXHzIMABAIyh+65b/1BpN9AtLS3NrS4CHADAGJrNZjQajUM9cN1AZx04AIBTaGVlJdrt9ivH2u323HdiEOAAAMa0ubkZt27deuXY1tZWbG5uzrUeZqECAIf8fvXjuqvwwpMnT+LixYtT/cxWqxU7OztxcHAQGxsbsba29qIXbXt7+8VWWZ9//nl89tlncf369YjoDKN291FtNpvRbrdjY2NjrsOnEQIcALCARvWYra+vx/r6+tDz8x4uHcQQKgBAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABTmtborAACcQn++VXcNXnj9+bOI1y+8evDTm/VU5pTQAwcALJxWqxVXrlyJlFJsbGzEzs7Oi3Pb29uxuroaKaVYW1uL7e3tGms6mAAHACyczc3N2NjYiEajEVtbW7GysvLi3Pr6emxubkZExBdffBHr6+sDP6PVasXe3t5c6tvPECoAwJja7XZsbm5Go9F40VNXh4kCXEppPSIu9xx6lHO+PaRcIyIOqq/tnPNXIz5zrLIAAHVqNpuxtbUVERFffVVfXBk7wKWU7kbEX3POrd5jKaX1nPN2z7HN6At2KaXNlFKzP+xNUhYAgI6xAlwVtPYHhKqViPiyp1wzItZzzpd6C+WcWymlxyml7ZzzwaRlAQB46cgAVwWtG/Hq0GlERPSHr4jYiIid/nKVdkSsR8TtY5QFAJi6g4ODuH37cNx4+PBhDbUZ3zg9cBvReS+tPUbZlRgdylbjZSibpCwAwNQ1Go24cePGoeN7e3uncvmQrnGWEVmJiL2IiJRSI6V0veqVG2Q5Ih4NOdeOiKvHLAsAQGWcALccEfsppevxstesUU1gGBbkhmnMqCwAwMIYdxbqUkQc5Jy7Q557KaVWRDxIKX2Ucz5IKXUD15ETDyYpO+Bn16Pzfly89957sbu7O+lHTOz8L8/iNz/+/chyu7v/mnldFtHTp0/n8pwZzjOol/tfr7N6/99+++148uTJ0POvP382x9qMln/J8ayvPs9H1H1cP/74Y0TEwPvwww8/vDh37ty5gT//yy+/xL///e+R93HYdU/6OzUywPUEreWe8BYRETnndkqpHRE3I6J16IdnpFqyZDsi4urVq/natWszv+adP/4p/vnGR0eW+921j2del0W0u7sb83jODOcZ1Mv9r9dZvf/ffPNNXLx4cXiB/r1Ha/Ts+bO40FefC6PqPqY33ngjImLgfXjzzTdfnBt2n371q1/Fr3/969H3cch1P/nkkwlr23ftUSd7lvEYtk9EOyKu95U9cuhzkrIAALxqnCHUg4jYH3JuPyLGfQ+uu9vCtMsCANP26c26a/DC8ydPptLj1qvVasXOzk4cHBzExsZGrK2tvdgPdXt7O+7evRsREZ9//nl89tlncf369YjoLDty69ataLfb0W63o9VqxcrKSqyurr6yn+qsjRPg2tF5B26Y3qC1EwPWi6ssRcT9Y5YFAJia7mb1g6yvrw/dwL7RaIz82XkZZxbql9GZiTpIM14NWnsxvEeuGRF3j1kWAIDKOAFuOyKaQ5YMWYmI3hi6FQPWb6smQyxHxJ1jlgUAoHJkgKsmHLSiE7heSCltRcR27+zUareG7Wrv1F6bEdHq3dt0krIAALw01jpwOefbKaWDKrR1PaiW9Ogv20opraeUbkTn/bjLEXEv5/zVScoCANAx7kK+MSiszbssAADjvQMHAMApIsABwILKOdddhYUzrXsuwAHAAjp37lz89NNPdVdj4fz0009D91adhAAHAAvo4sWL8f3339ddjYXz/fffT7x36iACHAAsoKWlpXj8+HF899138fz5c8OpM5RzjufPn8d3330Xjx8/jqWlURtcjWfsWagAwNlx4cKF+OCDD2J/fz/+8Y9/xM8//1x3lYb68ccf44033qi7Gidy7ty5uHjxYnzwwQdx4cKFE3+eAAcAC+rChQvx/vvvx/vvv193VUba3d2NTz75pO5qnCqGUAEACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCvHbcH0wp3cs5rw45tx4RjYg4qL62c85fnbQsAADHDHAppesRsTLk3GZEPMo53+49llJq9h6btCwAAB0TB7iUUiMiPhtyrhkR6znnS73Hc86tlNLjlNJ2zvlg0rIAALx0nHfg1iPiyyHnNiJiZ8i5dvWzxykLAEBlogCXUlqOTuga1jO2Ep3wNUg7IlaPWRYAgMqkPXArOee9EeeXI+LRkHPtiLh6zLIAAFTGDnDVbNHtE16vMaOyAAALY6xJDNWEg/1RkwqqyQ0Rw4dXj1V2wM+uR/V+3HvvvRe7u7uTfsTEzv/yLH7z49+PLLe7+6+Z12URPX36dC7PmeE8g3q5//Vy/+vnGRw27izU66dlWY+c83ZUPYFXr17N165dm/k17/zxT/HPNz46stzvrn0887osot3d3ZjHc2Y4z6Be7n+93P/6eQaHHTmEWq35duTCuj29c0cOfU5SFgCAV40McNVQ51LOedhs0Ul0d1uYdlkAgIVy1BDqekT8NqV0pe94MyIipbRVfX8357wTnSVGLg/5rKWIuN/z/SRlAQCojAxww957626llXPe6Du1F53lQQZpRsTWMcsCAFA59mb2Q2xFxIP+g9VQ7HJE3DlmWQAAKsfZSiuiM8TZuxxIRERU78ptV5vU99qMiFbvMiSTlAUA4KWJeuCqodPPorMNVkTE1ymldkR83g1c1Wb06ymlG9GZiHA5Iu7lnA/NZJ2kLAAAHRMFuCpYjbOkyNg7NkxSFgCA4w+hAgBQEwEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIV5re4KwMz9+VbdNZiOT2/WXQMATgk9cAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFCY1+quADCmP9+q79pPP5ru9T+9Ob3PAlhAeuAAAAojwAEAFGbsIdSU0o2IeCcilqtDWznnr4aUXY+IRkQcVF/b0ygLAMCYAS6ltBkRt3LOB9X3yxHxdUrps5zz2oCyj3LOt3uPpZSavccmLQsAQMeRAa7qIdvqhreIiJzzXkppLSLupZSud3vMUkrNiFjPOV/q/Yyccyul9DiltN0TAscuCwDAS+O8A3c559zuP5hz3onOsOdGz+GNiNgZ8jntiFg/ZlkAACrjBLgbKaW7Q87dj4irPd+vRCd8DdKOiNVjlgUAoDJOgBvWSxYRsRQR+z3fL0fEoyFl2/Fq2JukLAAAlSPfgcs5j+oJW46I7Qmu15hRWQCAhXHsnRiqyQ0REZvV993AdeTEg0nKDrnuekTEe++9F7u7u5N+xMTO//IsfvPj348st7v7r5nXZRE9ffr0ZM/56UdTq8uievrLhdid5n2cQ7s9S07cBjgR979+nsFhxwpwVQDbjIi1QRMcZinnvB1Vr9/Vq1fztWvXZn7NO3/8U/zzjaP/4/W7ax/PvC6LaHd3N070nOvcguqM2H36UVz7L0f/I2Zs1/7X9D5rAZy4DXAi7n/9PIPDjrsTw92IaPUuuNuz5MeRQ5+TlAUA4FUTB7hq8d27VU/YJLq7LUy7LADAQpkowFXbaT0cEd52IuLykHNL0Vl25DhlAQCojB3gUkrXI+KgP7yllFZ6vt2LiOaQj2hGZ+j1OGUBAKiMFeC6IW1Iz1vvMiNbMWD9tmrSw3JE3DlmWQAAKuPshbocnW2vvqx64bqWovOu2ouJCDnndkppO6W0mXNu9ZTdjM6kh4PjlAUA4KVxlhH5Ojoh7fqQ873hq7sZ/Xr1vtxBdN5zu9c7Y/U4ZQEA6BhnJ4ZLk37oJDNUjzGbFQBgoR13HTgAAGoiwAEAFEaAAwAozLE3s+eMO037hz796HTVBwBqpgcOAKAwAhwAQGEMoQLzd5aGxD+9WXcNgAWkBw4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAUxjIiU/SHe38bq9zvVz+ecU0AgLNMDxwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACiMAAcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACFEeAAAAojwAEAFEaAAwAojAAHAFAYAQ4AoDACHABAYQQ4AIDCCHAAAIUR4AAACiPAAQAURoADACjMa3VXAKBof741+2s8/Wg+1/n05uyvAUyFAFeDP9z725Flfr/68RxqAgCUyBAqAEBhBDgAgMIIcAAAhRHgAAAKI8ABABTGLNQp+h/fbk/ts/7yf8Yr9x/Nd6Z2TQCgDHrgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEsIwJAx59v1V2D6fn0Zt01gJnSAwcAUBgBDgCgMAIcAEBhBDgAgMIIcAAAhRHgAAAKI8ABABTGOnAAnD3TXNPu6Uf1rZFnPTuG0AMHAFAYAQ4AoDCGUAv3l/ajscr9R/OdGdcEAJiX2gNcSmk9IhoRcVB9beecv6q3VmePoAcAZ0etAS6ltBkRj3LOt3uPpZSavccAAHiptgCXUmpGxHrO+VLv8ZxzK6X0OKW0nXM+qKl6AACnVp2TGDYiYmfIuXZErM+xLgAAxagzwK1EJ6gN0o6I1TnWBQCgGHW+A7ccEV8OOacHribjTnYY1ziTIo665g9v/7f4y38+MsECoGQnWQy5zsWUBzkFCyzXPgt1hEbdFeDkph0IARbKaQotnCq1BLiUUjecTTxJoVp2pNs79zSl9H+nVrHh3o2I7+ZwHQZz/+vnGdTL/a+X+1+/U/YM/ve8LvTfh504zT1wA+WctyNie57XTCndzzlfnec1ecn9r59nUC/3v17uf/08g8NqmcTQszyIYVIAgAmd1r1QuzszAADQp84AtxMRl4ecW4qI+3Osy1HmOmTLIe5//TyDern/9XL/6+cZ9Ek553ou3NlGaznnfGi9t5TSg4jYqt53AwCgR50BrhkRD/q30qpmqD6OiEu20gIAOKy2ABfxohcucs6tnmNbEfHQZvYAAIPVGuAiXqzr1p20cDki/ppz/qrWSlX66taIiPZpqdtZl1K6G50dObZyzu2qx3YlIq7knDfqrd3ZUvV6b0bEvVG/39rDbIxz/7WH2Uop3YiId6KzQ1BE5z4PexbawQyM+wy0hZdqD3CnVdU7+Ki3J3DQMWajeg9yue/wzqB3Jjmeqrd7KSL+Gp0AsTbiP1raw5RNeP+1hxmpfo9vdV/ZSSktR8TX0bm/awPKagdTNuEz0BYqAtwAw97Pq849joiPvJ83W1WDvhcRzerQ/ZzzXo1VOtNSSjmGBAjtYfZG3f/qvPYwA1Vv2k7Oud13fCU69/vFM9EOZmOSZ1Ad1xYqxe3EMCcb0VnmZJB2dLby8q+t2XqUcx72DJgv7aF+2sNsXB602kHOeSeldBCd3/1ueNAOZmOSZxChLbxwWhfyrdtKdBrkIO2IWLiuWhaa9sBZdaN6p2qQ+xHRu3WTdjAbkzwDeghwgy1HxKMh59rhF2puUkqNlNJKNXxBPbSHU0J7mLpRPTlLEbHf8712MBuTPIMXtAUB7rjs4Tp771TvRlyNzr/CGimle9WMPU4Xz2T2tIcZyDmv9r8k32M5RoeLfp7FMRzjGWgLFe/A9en5JfAyav3u9LwUvFe9vPoghm/BxpRpD6eK9jAnVUCI6MwO1g5q0P8M+mgLoQeOUyrn3Oqf0VW9uLrU07BhIWgP89OzLt9a/8xI5mPUM9AWXhLg+vT8Yixcd2wh2hExrLudKdMeTj3tYfruRkSrd+kK7WDuDj2DMSxcWxDgJtddgZt6tOPl+j/UT3uol/YwRdVQ3N1By1ocQTuYkhM8g4VrCwLcYDsxfCx9KTovTjIjKaUH3X1yORW0hxppD/NRbeX0cERw0A5m7KhnoC28SoAbbC+GJ/lmdLp3mZ1GRDwccq4ZnefD/GgP9dIeZiyldD0iDvqDQ7UbQJd2MENjPgNtoYcAN9hWDFjTp3qxcjki7sy9Rotla9C/wKr1fpaj83yYH+2hXtrDDHUDwpBen97FebWDGZnkGWgLL9kLdYhuN23OudVzbCs63bu2S5mh6v8Qb/be++r4vYjY6z/OyVV7cW6MGLrQHmZo1P3XHman2jT9ZkR82XdqKTq9PZdzzhs95bWDKZvkGWgLrxLgRqimJHdfTr0cEX+dcFYMx9RtqNW3jeg05i174E1P9b7Jb6PzL9dmdH7P70fnZeBDU/W1h+ma5P5rD7NRbUI/amZpqz+YaQfTNekz0BZeEuAAAArjHTgAgMIIcAAAhRHgAAAKI8ABABRGgAMAKIwABwBQGAEOAKAwAhwAQGEEOACAwghwAACF+f8q/7DL5I8OkwAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ "<Figure size 735.243x454.405 with 1 Axes>"
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnAAAAGbCAYAAABEeK1/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAegklEQVR4nO3dz25cV2In4N8Zt60GMjIo2gOjEcDTpoDeh1aAzE5GyDcQ2/MCJtfZiNC8gEBtsib9BLa8yyoQDXM3i7aY9QRQdaaBZDGwJMJSgLYN95lF3ZJKpSqyimKxeMzvAwy57j1175H7tPjT+VtqrQEAoB3/ZdEVAABgNgIcAEBjBDgAgMYIcAAAjRHgAAAa86tFVwBgklLKWpKNJI+6S+8leZDk2yQ73fX3kvyh1vrVQioJsADFNiLARVRKeZCkV2vdGrq2muTr7uPf11oPSykPk6zUWq8top4Ai6AHDrhwSim7SZZrresjt9aSLCXZr7Uedte+SL8Xbt51Wurq1Jv3u06ox60k68PB9piyt5Ncz8ueytRat+dbQ+A8CHDARbSZ5LWAUmu9V0rZSX8Y9cW1c6rTjfTD47kHuFLKSpLVJJ8muZVkb4rvPEhyf6QHc62U8rDW+vHcKgucC4sYgAulCyvJ8UHp6DzqMmLl5CJnr+tx2+ne/9mU39lJklrrK0Gv1rqf5NtSyv2zridwvvTAARfVk1Pem5eNJLvn/dJuccaLBRqllGm+djv9+o6zm+RhKWWp1rqIIAycAT1wACfoegXXFl2PaXQ9dklyOKHIoGdz8xyqA8yJAAdwjG7xwoMTC14c60kyabHFUK/b6AIRoCGGUIG5G8zJSvI4/VWRR0nuDg/hdT1HW0mWu0v3SymjIWRwb6eUMpicfzhYWVlKGe1VWkqyn2RrdNVmF8x2uroM6vVweN5Y97yNvJxzN/zeJNkeWg17UUwzV683ZTnggrIPHDA33dDj/fSDzv7Q9dUknyf5bDQAdd95lP5WGfsZUUqp6QeyvZHru0l2xzxvN8larfX6yPvvd+/oDV1/be+5ofIPk2wsesPg7ve/N2kbkW5fvNVa68TJcqWUR+lviWLvPGiUIVRgnh6kH6peCWJdyNpO8nXXE/ZGutB3Y1xv2ISg83WSnTHDjFtJNrsTIFq1lOlW6b7xf3dgcQQ4YC66TWRXRnvKBrpQ9yT9Ycw3NdgnbZLh3r+dJEvj6tUFukG4nEkp5XYppZ7BP28arJZPLrKQbViAMyTAAfOylaHgNMF+zmY15LdJf/hwQu/ZcEjczPF7zPXS37R3JrXWe7XWcgb/nEe40vsGjbOIAZiXlUzeymLgKOkPgb7JEVW11qNSynr6Q7YPur3S9rvPe4Nnd71bS0mejFnwMDAo06onmW6Bgl44aJgAB5y5btJ/Mv2Guyt5wyOqaq37pZRr6fewrae/b9takjullI1uyHYQbPYnDe1mimOqLrhpgtlyFrMZMnBGDKEC8zAIY9PMxxoufyqllKXByQLdUOZ6twpzvXv2g+7+oEdw2nqd9N5bJ5c6d4Ph5ON6EZdycu8ocIEJcMCZ6+ZxHeXkobyVrvybHhB/I8nvx9Rjvzu4/Sgv59qd5R5ofzv4lwu0iOFh9+vY3+PQ81vanBgYIcAB87KXZPWEQHIrQ+d8vqHjThb4cujfv8rxK1bHbQh84rDkBVrEMPi9TlqIMVjkcdICE+ACE+CAuehORxju+XrF0PDjZ6d4/LhQeOuYsLicbsiwq1dv6HSI0Xq9tohhqIfwopxeMHEIuAuA99JfBTzOpxla2AG0SYAD5unjJFujc8W6RQ476Z+EMNrjtDTy6/D3BtfeG/Ouo4zZU65719LIZsLr6Qe+0XotJblTa7035vl7GQpF3XYlX4wpNzdDv/9jh1knHS3W/X5Xc4p97oCLxVFawNx1vV2P82rwGj0LdSn9ExJW8vI0gSfpjrvqjoga3Ev6c9kOa60bXZganHt6p7v/ePDsCYFsUK/hFbCPJ5Xtyu925R8kOTpmJeuZ6X5v2917h3sAj9Kv95Na69jh424z5evpH032XvIy3AFtE+AAABpjCBUAoDECHABAYwQ4AIDGCHAAAI1p+izU999/v/72t7+d+3v+8z//M3/1V3819/fAJNogF4F2yKJdtjb48OHD72qt/23cvaYD3G9/+9t8++23c3/PwcFBbt68Off3wCTaIBeBdsiiXbY2WEr5v5PuGUIFAGjM1D1w3YaQ7+XlJpqPjtkcczMvN+JcStKrtY4973CWsgAATBngut3Hd4bPziulrJVSHozuAD7YcX043JVSdkopK6OBb5ayAAD0nRjgumNcHowefFxr3S+lrJdS1gZnDJZSVpJs1lqvjZTdLqU8LaXsDY7OmaUsAAAvTTMHbj394c1x/pD+wcgDW+mfRThOL8nmKcsCANCZJsA9TrIz4d56ksOhz2t5eSj0qF5X/jRlAQDoTBPg9pKslFIelVJe9LYN/n0wfNpZTT/wjdNLcuOUZQEA6Jw4B67WelRK+fskXyd5WEq5l/7QaWqtWzO+b+nkIqcqCwBwaUy1CrXWelhK+SjJH5PcTr+HbGO4TCllELhOXHgwS9kx391MNz/ugw8+yMHBwayPmNnz58/P5T0wiTbIRaAd/vKUUnLlypW88847KaUsujonevfdd/Mv//Ivi67GqdRa8+OPP+aHH35IrfWNnzftNiKrST5N8lH6c9c+T783buO892yrte6lP6ybGzdu1PPYkfmy7fzMxaMNchFoh78sP/zwQ/70pz/l2rVreffdd/P2229f+BD37NmzXL16ddHVmFmtNT/99FO+//77PH36NB9++GGuXLnyRs88cQ5cF97u1Fq3a61HXWD7KP0VpPdLKbe6yg16004c+pylLABw9p48eZJr167l/fffb6YHrlWllLzzzjt5//33c+3atTx58uSNnznNIobPk3w2fKELcuvp94R9PuW7BqctnHVZAGBGz549y7vvvrvoalw67777bp49e/bGzzk2wA3mqk3aUHewiGFoTtt+kusTHrecZPjk+VnKAgBn6Oeff87bb7+96GpcOm+//XZ+/vnnN37OWRxm3xsKeIdJViaUW0lyf+jzLGUBgDNm2PT8ndV/82MDXBfMlrpjr8ZVYimvbsa7mzH7t3XlVpN8ecqyAAB0plmFup7+YoWNkcPsV9I/oeHF/Lhaa6+UsldK2am1bg89YyfJ9vBQ7CxlF+3/Pfsh//jgX6cu/w/rv5tjbQCAy26ajXx73Ua+O1233yBYPa61bowpv11K2Syl3O7KXk/yYNx2I7OUBQCgb9qNfI/SP3x+Kt1ebWdeFgA4H7OMPJ2XH3/8Ie+8098/7bKPdp3FIgYAgKZsb2/n448/TiklW1tb2d9/ebT73t5e1tfXU0rJxsZG9vYuXl/TVD1wAAC/JDs7O9nb28v29nZ2d3dfube5uZkbN27k448/zueff56lpVfPHTg8PMz+/n5WVlby5MmTrKysZG1t7TyrL8ABAEyr1+vl7t27uX//5W5nGxsbWV5ezurq6rnVwxAqAMCUdnZ2srX16rKAO3fuZHt7e8I35kOAAwCY0pdffpmVlVe3x11ZWXllDt15EOAAAKbQ6/VydHT0WoAbzJE7PDw8t7qYAwcAXFpHR0e5d+/ea9cfPXo0tuxxnjx5cmb1OokABwBcWktLS7l9+/Zr1w8PDy/k9iEDhlABAGZwUk/ceRDgAACmMJjrNjpUOgh0y8vL51YXAQ4AYAorKytZWlp6rQduEOjsAwcAcAGtra2l1+u9cq3X6537SQwCHADAlHZ2dnL37t1Xru3u7mZnZ+dc62EVKgDwmn9Y/92iq/CaZ8+e5erVq2fyrO3t7ezv7+fo6ChbW1vZ2Nh40Yu2t7f34qiszz77LJ9++mlu3bqVpD+MOjhHdWVlJb1eL1tbW+c6fJoIcADAJXRcj9nm5mY2Nzcn3j/v4dJxDKECADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAY3616AoAABfQN3cXXYPXvPPjD8k7V/ofPrmz2MosmB44AODS2d7ezscff5xSSra2trK/v//i3t7eXtbX11NKycbGRvb29hZY0/EEOADg0tnZ2cnW1laWlpayu7ubtbW1F/c2Nzezs7OTJPn888+zubk59hnb29s5PDw8l/qOMoQKADClXq+XnZ2dLC0tveipWwQBDgBgSisrK9nd3U2SfPXVVwurhyFUAIDGCHAAAI0xhAoAXFpHR0e5d+/ea9cfPXq0gNpMT4ADAC6tpaWl3L59+7Xrh4eHF3L7kAFDqAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGKtQAYDXfXJn0TV4zY/PnuXK1atn8qzt7e3s7+/n6OgoW1tb2djYeHEe6t7eXu7fv58k+eyzz/Lpp5/m1q1bSfrbjty9eze9Xi+9Xi/b29tZW1vL+vr6K+epzpsABwBcOoPD6sfZ3NyceID90tLSsd89L4ZQAQAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOAC4pGqti67CpXNW/80FOAC4hN5666389NNPi67GpfPTTz/lrbfeeuPnCHAAcAldvXo133///aKrcel8//33uXoGp0kIcABwCS0vL+fp06f57rvv8uOPPxpOnaNaa3788cd89913efr0aZaXl9/4mY7SAoBL6MqVK/nwww/z5MmT/Nu//Vt+/vnnRVfpRH/+85/z61//etHVOJW33norV69ezYcffpgrV6688fMEOAC4pK5cuZLf/OY3+c1vfrPoqkzl4OAgf/M3f7PoalwIhlABABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0ZqajtEopm0muD116XGu9N6HcUpKj7tderfWrY545VVkAAGYIcKWU+0n+UGvdHr5WStmste4NXdvJSLArpeyUUlZGw94sZQEA6JsqwHVB68mYULWW5IuhcitJNmut14YL1Vq3SylPSyl7tdajWcsCAPDSiXPguqB1O8nO6L1a67WR4c6tJPsTHtVLsnnKsgAAdKZZxLCV/ry03hRl19IPX+P0kqyfsiwAAJ1pAtxaksMkKaUslVJudb1y46wmeTzhXi/JjVOWBQCgM02AW03ypJRyK/0wt59kqVvAMCnITbI0p7IAAJfGtKtQl5Mc1VoHc9YOSynbSR6WUj6qtR6VUgaB68SFB7OUHfPdzXTz4z744IMcHBzM+oiZvf2XH/LXf/7j1OUPDv5jjrXhMnr+/Pm5tHU4jnbIommDLx0b4IaC1upQeEuS1Fp7pZRekjtJtl/78px0W5bsJcmNGzfqzZs35/7OL//pn/Pvv/5o6vK/v/m7OdaGy+jg4CDn0dbhONohi6YNvnTsEOrQNh6HE4r0ktwaKXvi0OcsZQEAeNU0c+COkjyZcO9JkmnnwQ1OWzjrsgAAl8o0Aa6X/hy4SYaD1n5ePWpr2HKSb09ZFgCAzjQB7ov0V6KOs5JXg9ZhJvfIrSS5f8qyAAB0pglwe0lWJmwZspZXT2jYzZj927rFEKtJvjxlWQAAOicGuG7BwXb6geuFUspukr3h1andaQ173dmpw3aSbA+fbTpLWQAAXppqH7ha671SylEX2gYedlt6jJbdLqVsllJupz8/7nqSByNnps5cFgCAvmk38s24sHbeZQEAmG4OHAAAF4gABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMb867RdLKQ9qresT7m0mWUpy1P3aq7V+9aZlAQA4ZYArpdxKsjbh3k6Sx7XWe8PXSikrw9dmLQsAQN/MAa6UspTk0wn3VpJs1lqvDV+vtW6XUp6WUvZqrUezlgUA4KXTzIHbTPLFhHtbSfYn3Ot13z1NWQAAOjMFuFLKavqha1LP2Fr64WucXpL1U5YFAKAzaw/cWq318Jj7q0keT7jXS3LjlGUBAOhMHeC61aJ7b/i+pTmVBQC4NKZaxNAtOHhy3KKCbnFDMnl49VRlx3x3M938uA8++CAHBwezPmJmb//lh/z1n/84dfmDg/+YY224jJ4/f34ubR2Oox2yaNrgS9OuQr11Ubb1qLXupesJvHHjRr158+bc3/nlP/1z/v3XH01d/vc3fzfH2nAZHRwc5DzaOhxHO2TRtMGXThxC7fZ8O3Fj3aHeuROHPmcpCwDAq44NcN1Q53KtddJq0VkMTls467IAAJfKSUOom0n+tpTy8cj1lSQppex2n+/XWvfT32Lk+oRnLSf5dujzLGUBAOgcG+AmzXsbHKVVa90auXWY/vYg46wk2T1lWQAAOqc+zH6C3SQPRy92Q7GrSb48ZVkAADqnOUor6Q9xDm8HkiTp5srtdYfUD9tJsj28DcksZQEAeGmmHrhu6PTT9I/BSpKvSym9JJ8NAld3GP1mKeV2+gsRrid5UGt9bSXrLGUBAOibKcB1wWqaLUWmPrFhlrIAAJx+CBUAgAUR4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACN+dWiKwBz8c3dRdfg7HxyZ9E1AOCC0QMHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANOZXi64AF8Q3dxddAwBgSlMHuFLK7STvJVntLu3WWr+aUHYzyVKSo+7X3lmUBQBgygBXStlJcrfWetR9Xk3ydSnl01rrxpiyj2ut94avlVJWhq/NWhYAgL4TA1zXQ7Y7CG9JUms9LKVsJHlQSrk16DErpawk2ay1Xht+Rq11u5TytJSyNxQCpy4LAMBL0/TAXa+17o1erLXul1KOkmwlGQx5biXZn/CcXpLNJPdOURYur2/uJs8/+mXMU/zkzqJrAPCLMM0q1NullPsT7n2b5MbQ57X0w9c4vSTrpywLAEBnmgA3qZcsSZaTPBn6vJrk8YSyvbwa9mYpCwBA58Qh1FrrcT1hq0leG149xtKcygIAXBqn3geuW9yQJDvd50HgOnHhwSxlJ7x3M0k++OCDHBwczPqImb39lx/y13/+49TlDw7+Y461mZPnHy26Bhzj+V+u5OCX8L/ROfz/lfl5/vz5ufyZC5Nogy+dKsB1AWwnyUatddI8trnoFlTsJcmNGzfqzZs35/7OL//pn/Pvv57+h+fvb/5ujrWZk1/CBPlfsIPnH+Xmf53+LxEX1s3/uega8AYODg5yHn/mwiTa4EunPUrrfpLt4Q13h7b8OHHoc5ayAAC8auYA122+e3/c1iInGJy2cNZlAQAulZkCXHec1qNjwtt+kusT7i2nv+3IacoCANCZOsCVUm4lORoNb6WUtaGPh0lWJjxiJf2h19OUBQCgM1WAG4S0CT1vw9uM7GbM/m3doofVJF+esiwAAJ1pzkJdTf/Yqy+6XriB5fTnqr1YiFBr7ZVS9kopO7XW7aGyO+kvejg6TVkAAF6aZhuRr9MPabcm3B8OX4PD6De7+XJH6c9zezC8YvU0ZQEA6JvmJIZrsz50lhWqp1jNCgBwqZ12HzgAABZEgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAY3616Ar8Ev3jg3+duuw/rP9ujjUBAH6J9MABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojKO0gPPzzd1F1+DsfHJn0TUALjE9cAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANAYAQ4AoDECHABAYwQ4AIDGCHAAAI0R4AAAGiPAAQA0RoADAGiMAAcA0BgBDgCgMQIcAEBjBDgAgMb8atEVaNnf/WnvzR/yzXtv/gwA4FLRAwcA0BgBDgCgMQIcAEBjBDgAgMZYxABwGt/cXXQNzs4ndxZdA2BGeuAAABqjB27B/nfv8Uzl/8eKbUcA4LLTAwcA0BgBDgCgMQIcAEBjBDgAgMYIcAAAjRHgAAAaI8ABADRGgAMAaIwABwDQGAEOAKAxCz9Kq5SymWQpyVH3a6/W+tViawUAcHEtNMCVUnaSPK613hu+VkpZGb7GS85OBQAWFuBKKStJNmut14av11q3SylPSyl7tdajBVUPAODCWuQcuK0k+xPu9ZJsnmNdAACascgh1LUcH+DWkxhGPUeGZ+GS+ubudOWefzR92UX55M6iawDnYpEBbjXJFxPu6YE7I7OGMoCmXfSAOSuBlAkWvgr1GEuLrgDHm2c4nLV3T+8hwAV3FuH6ovQCX4BgvZAAV0oZhLOZFyl0244Meueel1L+z5lVbLL3k3x3Du+BSbRBLgLt8Nz9r0VX4KK5IG3w3P53+e+TblzkHrixaq17SfbO852llG9rrTfO850wTBvkItAOWTRt8KWFrEId2h7EMCkAwIwu6lFag5MZAAAYscgAt5/k+oR7y0m+Pce6nORch2xhDG2Qi0A7ZNG0wU6ptS7mxf1jtFZrretj7j1MstvNdwMAYMgiA9xKkoejR2l1K1SfJrnmKC0AgNctLMAlL3rhUmvdHrq2m+SRw+wBAMZbaIBLXuzrNli0cD3JH2qtXy20Up2Rui0l6V2UuvHL0vU87yR5cFwb0yaZp1LK7STvpX9STtKfyjK2fWmLnLXuz8HN9LPAUZKV9E9mujtuRO6yt8GFB7iLqusdfDzcEzjuGryJrsd5Ockf0g9wG8f8wNQmmZuuLb34QVlKWU3ydZL9WuvGmLLaImemC293hkfkuuv30/8LxcfDIU4bFODGmjQ/r7v3NMlH5udx1kopNRMCnDbJPHU9Gfu11t7I9bUkDzLULrVF5qHr/f00yd+PBLXVJA+TbA0WNmqDfRd1H7hF20p/m5Nxenl5lBecF22Sebo+Gt6SpNa6n/7w1NbQZW2ReeilP2Q6DW0wAtwka+k3gnF6SV7b+gTmTJtknm53Q1XjfJtk+OgibZEzV2v9qtY6bveJQdsbDmzaYAS4SVaTPJ5wr5dX/zCD86BNMk+TejOS/hzNJ0OftUXO03aSvZEeYm0wDR5mf0E4w5WLRpvk1MZtqD5kNbPtfq8t8sa6eW47Sb4aXdgwhUvRBgW4Ed1KmMRZrFwQ2iSL0i1uSPo/SLVF5m5oK5H30l+dP7qwRhvsCHAAvGZob8KNcQscYB66OXDDW4PsllK2TuglvpTMgRsxNIHyUnTBcvFpkyzI/STbw9vaaIsswHaStaGTm7TBjgA3u8Guz3BRaJOcqe6H5f3Bvlsz0BY5U11gO8z0W4NcmjYowI23n/5RHuMsp7+sHs6TNsm56DZUfXRMeNMWOXOllEeDXrYxniRZGpr/pg1GgJvkMJM3FFxJf2gBzpM2ydyVUm4lORoNb92JDAPaImeqC2bHbeK7kn67HPSsaYNxlNZYk47p6BrZ0yTjNhuEN3Kao7S0Sc5KF9KWJrS/ncFWDtoi81BK2a21bo25PmhXe4P72mCfADfB0ITJ7aFru+kPLVyKg3I5X12A25o0dKVNMi/deZN3knwxcms5/TlF14d/uGqLnLUulG2PhrjuhJCVWuvHI9cvfRsU4I7R7YE0mBB5Pckfxv3tFE6rm2/0t+lvlrqSflv7Nv29j7ZH/xapTTIP3QHgx63q2x79oagtcta6HrQ7Q5dW0m9XYwPZZW+DAhwAQGMsYgAAaIwABwDQGAEOAKAxAhwAQGMEOACAxghwAACNEeAAABojwAEANEaAAwBojAAHANCY/w+fkb4qHzhn4AAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ "<Figure size 735.243x454.405 with 1 Axes>"
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "for offset in [0,3,10]:\n",
+ " stat_h0,stat_h1 = dist_unknownAsigma_twopulse(A=1,sigma=5,x_1=30,x_2=60+offset)\n",
+ " plt.figure()\n",
+ " plt.hist(stat_h0,alpha=0.5,label='H0')\n",
+ " plt.hist(stat_h1,alpha=0.5,label='H1')\n",
+ " plt.legend()\n",
+ " plt.title('offset='+str(offset))"
+ ]
+ },
{
"cell_type": "code",
"execution_count": null,
--
GitLab