Add files via upload

This commit is contained in:
Peter Norvig authored and GitHub committed 2019-09-27 22:12:01 -07:00
1 parent af96a79962
commit fe11dbbe54
1 file changed
+24 -7
+24 -7
View File
@@ -234,7 +234,7 @@
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAykAAAKSCAYAAADWEJRoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsnXt0nFW5/7+TTJI2adq0aUjSlg4BQqVc9YBoIREBWy4pBDWCIiJKUBDOQdCjB8PtELmJB5dwFCkiChQhWlJaLsKPc4SmgrZcLBfF9BAmaTKZJulcMjOdmczl90eckMvkTebN+zx79vB81mItTSaf7P2Z3aGb92bzer1JCIIgCIIgCIIgZAl5qgcgCIIgCIIgCIIwHtmkCIIgCIIgCIKQVcgmRRAEQRAEQRCErEI2KYIgCIIgCIIgZBWySREEQRAEQRAEIauQTYogCIIgCIIgCFmFbFIEQRAErTnqqKNw1FFHzfr1TqcTZWVluOyyywhHJQiCIMwF2aQIgiB8iNi8eTO++93v4owzzsCBBx6IsrIyfO1rX5v29X19fbjvvvvwhS98AccccwwqKyvhcDhw5plnYuPGjUgmpz5q67LLLkNZWZnhP2effTblNAVBEATNsasegCAIgsDHj370I7z11ltYsGABli1bhuHhYcPX33ffffjJT36CFStW4KSTTsKyZcvgcrmwdetWXH755Xjuuefw4IMPTviZs846CytXrkzre+yxx/D+++/jM5/5jFVTwpNPPmmZSxAEQcgObPLEeUEQhA8PL730EpYvX46DDz4YHR0dWL9+PT772c/igQceSPv6J598EmVlZairq4PNZhv7usvlwmmnnYbe3l785je/mdWREa/Xi8MPPxzxeBx/+9vfUF5ebtm8MsHpdOKYY47BF7/4Rfz85z9XMgZBEATBGDndSxAE4UNEfX09DjnkkAkbDiPOPvts1NfXT3l9dXX12Gli27Ztm5Xrsccew/79+7F+/fpZb1C2bds2dv1IZ2cnLrroIhxyyCFYvHgxdu3aBWD6a1KGh4dx7bXXYvXq1aisrMTxxx+Pe+65J+0painC4TDuvvtufOpTn8Ly5cuxbNkynHzyyXjggQfS/twzzzyDc845Bx/5yEdwwAEHYNWqVVi3bh1+/OMfz2p+giAIQnrkdC9BEATBFAUFBQAAu312/yr59a9/DQC46KKLMv5dXV1dOPXUU7Fq1Sqcd955GB4eRnFx8bSvj0QiOOecc/Daa69h9erVaGpqgt/vx5133ont27en/Znh4WE0Njbi1VdfxdFHH40vfelLAIAXXngBV199NXbs2DHhyMuDDz6Iq666CgcccADWrVuHiooKDA0N4d1338WvfvUrXHPNNRnPUxAEQRhFNimCIAhCxsRiMfz2t78FAJx22mkzvv4vf/kL3nnnHRx66KGor6/P+Pe98soruOaaa3DdddfN6vX33HMPXnvtNZx55pl4+OGHkZc3euLAt7/9bZx88slpf+baa6/Fq6++ihtvvBFXXXXV2NcjkQguvPBCPProozj77LNxxhlnAAB+9atfobCwEB0dHTjggAMmuIaGhjKeoyAIgvABcrqXIAiCkDEtLS145513sG7dOpx66qkzvj51cb2ZoygAcMABB+Df//3fZ/36Rx55BDabDTfddNPYBgUAVq5ciW984xtTXu/xePDoo4/i6KOPnrBBAYCioiJcf/31AEZPWRuP3W4fO6I0HlXX2wiCIOQKciRFEARByIj/+q//wr333ovDDz8cv/jFL2Z8vc/nQ3t7OwoLC8dOocqUI488EkVFRbN67fDwMN577z1UVVWhtrZ2yvdPPPHEKV979dVXEYvFkJeXh1tvvXXK92OxGADgH//4x9jXmpqa0NLSghNOOAHnnnsu1qxZgxNOOAFVVVWznZYgCIIwDbJJEQRBEGbNnXfeidbWVhxxxBHYvHkzysrKZvyZxx9/HKFQCJ/73OdMH2GYfDqVEX6/HwBQUVExa9e+ffsAAG+88QbeeOONad2BQGDsf19xxRWoqKjAAw88gPvvv39sw3b88cfj+uuvR11d3azHLAiCIExETvcSBEEQZkVraytaW1tx7LHHYuvWrVi6dOmsfi51wfxXv/pV0797tncjA4CFCxcCAAYGBtJ+f+/evdP+zKWXXgqv1zvtP6k7iqU477zz8Ic//AFdXV3YtGkTmpub8de//hVNTU3YvXv3rMcsCIIgTEQ2KYIgCMKMXHvttbjzzjtx/PHHY/PmzVi8ePGsfm7nzp146623cOihh7IdWSgtLcXBBx8Mt9uddqOQ7u5exx13HPLy8vDyyy+b+p0LFy7EKaecgh/96Ee44oorEA6H8fzzz5tyCYIgCLJJEQRBEAxIJpO4+uqr8bOf/Qxr1qzBE088gUWLFs365+d6wbxZLrjgAiSTSVx//fVIJBJjX+/u7k57Hc3SpUtx3nnn4c0338Stt946dg3KeHp7eydck/Liiy+mfXaK2+0GAMNbJAuCIAjGyDUpgiAIHyK2bt2Kp556CsAHpz3t3LkTl1122dhrxj8L5Pbbb8cDDzyAwsJCHHfccbj77runOFeuXIkLLrhgytf9fj+eeOIJFBUVmb5g3ixXXHEFnnrqKTz99NOoq6vDaaedNjaeT37yk3jmmWem/Mwdd9yB9957D7fffjsee+wxrFmzBpWVlWNHZHbs2IEf/vCHOOywwwAAF154IUpKSnDcccdh5cqVsNlsePXVV/Hyyy+jpqYGjY2NrHMWBEHIJWSTIgiC8CHizTffxKOPPjrha93d3eju7h77/+M3KU6nEwAQjUbx05/+NK3zxBNPTLtJaWtrQzAYnNMF82YpKipCe3s7brvtNjzxxBO49957sXLlSlxzzTVYv3592k1KaWkptm7dioceeghtbW3YunUrwuEwKioq4HA4cMMNN+Dcc88de/2NN96I//mf/8Gbb76JF154AXa7HStWrMD3vvc9fOMb38joiJMgCIIwEZvX6516rFoQBEEQBEEQBEERck2KIAiCIAiCIAhZhWxSBEEQBEEQBEHIKmSTIgiCIAiCIAhCViGbFEEQBEEQBEEQsgrZpAiCIAiCIAiCkFXIJkUQBEEQBEEQhKxCNimCIAiCIAiCIGQVskkxwO12a+kGgM7OTjK3zl0o/ZTNAX27y1pX45e1zu+m9kt3fjcgnzEq3LLW1fizba3LJsUAyickcz992Up07iLdc8tNjc7rUbrzuzn8lOjaXZqr8Ut3fjeHnwoz487JTYrb7UYsFkM0GsXAwAAAwOfzIRAIAAD6+/sRj8cRiUQwODgIAPB6vQgGgwAAl8uFRCKB/fv3Y2hoCADg8XgQCoUAAH19fQCAUCgEj8cDABgaGkI4HEYikYDL5QIABINBeL1eAMDg4CAikQji8Tj6+/uRSCQQCATg8/kAAAMDA4hGo4jFYmO7Tb/fD7/fb+mcwuEw2ZwAZP2c9u3bp+2cEokE2fsUiUTI5uR2u0nXXgqK9ykcDmv75ymRSJD9eUpBNafh4eGM3qds+twbHBzMus+92c5pZGSE7HNvZGSEbE6pMVH9eUqRLe9TJnMaHh4mW3uJRIJsTimy8e8RM81paGiI9O8Re/fuJZtT6udn+z5lMqdoNDplTjNh83q9yRlf9SFlYGAAFRUV2rmB0UN2tbW1JG6du1D6KZsD+naXta7GL2ud303tl+78bkA+Y1S4Za2r8WfbWs/JIylWQbnIKN3U6NxFuueWmxqd16N053dz+CnRtbs0V+OX7vxuDj8VZsYtmxQDZnMoKhvd1OjcRbrnlpsandejdOd3c/gp0bW7NFfjl+78bg4/FWbGLZsUA/Lz87V0U6NzF+meW25qdF6P0p3fzeGnRNfu0lyNX7rzuzn8VJgZt2xSDFiwYIGWbmp07iLdc8tNjc7rUbrzuzn8lOjaXZqr8Ut3fjeHnwoz45ZNigGpOxjo5qZG5y7SPbfc1Oi8HqU7v5vDT4mu3aW5Gr9053dz+KkwM27ZpBig84VVlOjcRbrnlpsandejdOd3c/gp0bW7NFfjl+78bg4/FXLhvMXEYjEt3dTo3EW655abGp3Xo3Tnd3P4KdG1uzRX45fu/G4OPxVmxi2bFANSDxTTzU2Nzl2ke265qdF5PUp3fjeHnxJdu0tzNX7pzu/m8FNhZtyySTFg6dKlWrqp0bmLdM8tNzU6r0fpzu/m8FOia3dprsYv3fndHH4qzIxbNikGeL1eLd3U6NxFuueWmxqd16N053dz+CnRtbs0V+OX7vxuDj8VZsYtmxQDCgoKtHRTo3MX6Z5bbmp0Xo/Snd/N4adE1+7SXI1fuvO7OfxUmBm3bFIMKCkp0dJNjc5dpHtuuanReT1Kd343h58SXbtLczV+6c7v5vBTYWbcskkxwOVyaemmRucu0j233NTovB6lO7+bw0+Jrt2luRq/dOd3c/ipMDNu2aQYUFlZqaWbGp27SPfcclOj83qU7vxuDj8lunaX5mr80p3fzeGnwsy4ZZNiQDQa1dJNjc5dpHtuuanReT1Kd343h58SXbtLczV+6c7v5vBTYWbcskkxIBgMaummRucu0j233NTovB6lO7+bw0+Jrt2luRq/dOd3c/ipMDPunNykuN1uxGIxRKNRDAwMAAB8Ph8CgQAAoL+/H/F4HJFIBIODgwBGb42WCuhyuZBIJFBSUoKhoSEAgMfjQSgUAgD09fUBAEKhEDweDwBgaGgI4XAYiURi7Ly7YDA4dsu1wcFBRCIRxONx9Pf3o7y8HIFAAD6fDwAwMDCAaDSKWCwGt9sNAPD7/fD7/ZbOKRwOz3lO+fn5aecEYM5zstvtpHMCkNH7ZMWcrHqfysvLLVl76eZUWlpKNqdoNGrZ2ks3pxQU79OCBQtI/zzNmzePbO2Vl5eTfUaksGLtpZtT6i4wFH+e8vLyLP0snzynRCKhzWf55DmVlZWLine truncated
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAykAAAKSCAYAAADWEJRoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsnXt0nFW5/7+TTJI2adq0aUjSlg4BQqVc9YBoIREBWy4pBDWCIiJKUBDOQdCjB8PtELmJB5dwFCkiChQhWlJaLsKPc4SmgrZcLBfF9BAmaTKZJulcMjOdmczl90eckMvkTebN+zx79vB81mItTSaf7P2Z3aGb92bzer1JCIIgCIIgCIIgZAl5qgcgCIIgCIIgCIIwHtmkCIIgCIIgCIKQVcgmRRAEQRAEQRCErEI2KYIgCIIgCIIgZBWySREEQRAEQRAEIauQTYogCIIgCIIgCFmFbFIEQRAErTnqqKNw1FFHzfr1TqcTZWVluOyyywhHJQiCIMwF2aQIgiB8iNi8eTO++93v4owzzsCBBx6IsrIyfO1rX5v29X19fbjvvvvwhS98AccccwwqKyvhcDhw5plnYuPGjUgmpz5q67LLLkNZWZnhP2effTblNAVBEATNsasegCAIgsDHj370I7z11ltYsGABli1bhuHhYcPX33ffffjJT36CFStW4KSTTsKyZcvgcrmwdetWXH755Xjuuefw4IMPTviZs846CytXrkzre+yxx/D+++/jM5/5jFVTwpNPPmmZSxAEQcgObPLEeUEQhA8PL730EpYvX46DDz4YHR0dWL9+PT772c/igQceSPv6J598EmVlZairq4PNZhv7usvlwmmnnYbe3l785je/mdWREa/Xi8MPPxzxeBx/+9vfUF5ebtm8MsHpdOKYY47BF7/4Rfz85z9XMgZBEATBGDndSxAE4UNEfX09DjnkkAkbDiPOPvts1NfXT3l9dXX12Gli27Ztm5Xrsccew/79+7F+/fpZb1C2bds2dv1IZ2cnLrroIhxyyCFYvHgxdu3aBWD6a1KGh4dx7bXXYvXq1aisrMTxxx+Pe+65J+0painC4TDuvvtufOpTn8Ly5cuxbNkynHzyyXjggQfS/twzzzyDc845Bx/5yEdwwAEHYNWqVVi3bh1+/OMfz2p+giAIQnrkdC9BEATBFAUFBQAAu312/yr59a9/DQC46KKLMv5dXV1dOPXUU7Fq1Sqcd955GB4eRnFx8bSvj0QiOOecc/Daa69h9erVaGpqgt/vx5133ont27en/Znh4WE0Njbi1VdfxdFHH40vfelLAIAXXngBV199NXbs2DHhyMuDDz6Iq666CgcccADWrVuHiooKDA0N4d1338WvfvUrXHPNNRnPUxAEQRhFNimCIAhCxsRiMfz2t78FAJx22mkzvv4vf/kL3nnnHRx66KGor6/P+Pe98soruOaaa3DdddfN6vX33HMPXnvtNZx55pl4+OGHkZc3euLAt7/9bZx88slpf+baa6/Fq6++ihtvvBFXXXXV2NcjkQguvPBCPProozj77LNxxhlnAAB+9atfobCwEB0dHTjggAMmuIaGhjKeoyAIgvABcrqXIAiCkDEtLS145513sG7dOpx66qkzvj51cb2ZoygAcMABB+Df//3fZ/36Rx55BDabDTfddNPYBgUAVq5ciW984xtTXu/xePDoo4/i6KOPnrBBAYCioiJcf/31AEZPWRuP3W4fO6I0HlXX2wiCIOQKciRFEARByIj/+q//wr333ovDDz8cv/jFL2Z8vc/nQ3t7OwoLC8dOocqUI488EkVFRbN67fDwMN577z1UVVWhtrZ2yvdPPPHEKV979dVXEYvFkJeXh1tvvXXK92OxGADgH//4x9jXmpqa0NLSghNOOAHnnnsu1qxZgxNOOAFVVVWznZYgCIIwDbJJEQRBEGbNnXfeidbWVhxxxBHYvHkzysrKZvyZxx9/HKFQCJ/73OdMH2GYfDqVEX6/HwBQUVExa9e+ffsAAG+88QbeeOONad2BQGDsf19xxRWoqKjAAw88gPvvv39sw3b88cfj+uuvR11d3azHLAiCIExETvcSBEEQZkVraytaW1tx7LHHYuvWrVi6dOmsfi51wfxXv/pV0797tncjA4CFCxcCAAYGBtJ+f+/evdP+zKWXXgqv1zvtP6k7iqU477zz8Ic//AFdXV3YtGkTmpub8de//hVNTU3YvXv3rMcsCIIgTEQ2KYIgCMKMXHvttbjzzjtx/PHHY/PmzVi8ePGsfm7nzp146623cOihh7IdWSgtLcXBBx8Mt9uddqOQ7u5exx13HPLy8vDyyy+b+p0LFy7EKaecgh/96Ee44oorEA6H8fzzz5tyCYIgCLJJEQRBEAxIJpO4+uqr8bOf/Qxr1qzBE088gUWLFs365+d6wbxZLrjgAiSTSVx//fVIJBJjX+/u7k57Hc3SpUtx3nnn4c0338Stt946dg3KeHp7eydck/Liiy+mfXaK2+0GAMNbJAuCIAjGyDUpgiAIHyK2bt2Kp556CsAHpz3t3LkTl1122dhrxj8L5Pbbb8cDDzyAwsJCHHfccbj77runOFeuXIkLLrhgytf9fj+eeOIJFBUVmb5g3ixXXHEFnnrqKTz99NOoq6vDaaedNjaeT37yk3jmmWem/Mwdd9yB9957D7fffjsee+wxrFmzBpWVlWNHZHbs2IEf/vCHOOywwwAAF154IUpKSnDcccdh5cqVsNlsePXVV/Hyyy+jpqYGjY2NrHMWBEHIJWSTIgiC8CHizTffxKOPPjrha93d3eju7h77/+M3KU6nEwAQjUbx05/+NK3zxBNPTLtJaWtrQzAYnNMF82YpKipCe3s7brvtNjzxxBO49957sXLlSlxzzTVYv3592k1KaWkptm7dioceeghtbW3YunUrwuEwKioq4HA4cMMNN+Dcc88de/2NN96I//mf/8Gbb76JF154AXa7HStWrMD3vvc9fOMb38joiJMgCIIwEZvX6516rFoQBEEQBEEQBEERck2KIAiCIAiCIAhZhWxSBEEQBEEQBEHIKmSTIgiCIAiCIAhCViGbFEEQBEEQBEEQsgrZpAiCIAiCIAiCkFXIJkUQBEEQBEEQhKxCNimCIAiCIAiCIGQVskkxwO12a+kGgM7OTjK3zl0o/ZTNAX27y1pX45e1zu+m9kt3fjcgnzEq3LLW1fizba3LJsUAyickcz992Up07iLdc8tNjc7rUbrzuzn8lOjaXZqr8Ut3fjeHnwoz487JTYrb7UYsFkM0GsXAwAAAwOfzIRAIAAD6+/sRj8cRiUQwODgIAPB6vQgGgwAAl8uFRCKB/fv3Y2hoCADg8XgQCoUAAH19fQCAUCgEj8cDABgaGkI4HEYikYDL5QIABINBeL1eAMDg4CAikQji8Tj6+/uRSCQQCATg8/kAAAMDA4hGo4jFYmO7Tb/fD7/fb+mcwuEw2ZwAZP2c9u3bp+2cEokE2fsUiUTI5uR2u0nXXgqK9ykcDmv75ymRSJD9eUpBNafh4eGM3qds+twbHBzMus+92c5pZGSE7HNvZGSEbE6pMVH9eUqRLe9TJnMaHh4mW3uJRIJsTimy8e8RM81paGiI9O8Re/fuJZtT6udn+z5lMqdoNDplTjNh83q9yRlf9SFlYGAAFRUV2rmB0UN2tbW1JG6du1D6KZsD+naXta7GL2ud303tl+78bkA+Y1S4Za2r8WfbWs/JIylWQbnIKN3U6NxFuueWmxqd16N053dz+CnRtbs0V+OX7vxuDj8VZsYtmxQDZnMoKhvd1OjcRbrnlpsandejdOd3c/gp0bW7NFfjl+78bg4/FWbGLZsUA/Lz87V0U6NzF+meW25qdF6P0p3fzeGnRNfu0lyNX7rzuzn8VJgZt2xSDFiwYIGWbmp07iLdc8tNjc7rUbrzuzn8lOjaXZqr8Ut3fjeHnwoz45ZNigGpOxjo5qZG5y7SPbfc1Oi8HqU7v5vDT4mu3aW5Gr9053dz+KkwM27ZpBig84VVlOjcRbrnlpsandejdOd3c/gp0bW7NFfjl+78bg4/FXLhvMXEYjEt3dTo3EW655abGp3Xo3Tnd3P4KdG1uzRX45fu/G4OPxVmxi2bFANSDxTTzU2Nzl2ke265qdF5PUp3fjeHnxJdu0tzNX7pzu/m8FNhZtyySTFg6dKlWrqp0bmLdM8tNzU6r0fpzu/m8FOia3dprsYv3fndHH4qzIxbNikGeL1eLd3U6NxFuueWmxqd16N053dz+CnRtbs0V+OX7vxuDj8VZsYtmxQDCgoKtHRTo3MX6Z5bbmp0Xo/Snd/N4adE1+7SXI1fuvO7OfxUmBm3bFIMKCkp0dJNjc5dpHtuuanReT1Kd343h58SXbtLczV+6c7v5vBTYWbcskkxwOVyaemmRucu0j233NTovB6lO7+bw0+Jrt2luRq/dOd3c/ipMDNu2aQYUFlZqaWbGp27SPfcclOj83qU7vxuDj8lunaX5mr80p3fzeGnwsy4ZZNiQDQa1dJNjc5dpHtuuanReT1Kd343h58SXbtLczV+6c7v5vBTYWbcskkxIBgMaummRucu0j233NTovB6lO7+bw0+Jrt2luRq/dOd3c/ipMDPunNykuN1uxGIxRKNRDAwMAAB8Ph8CgQAAoL+/H/F4HJFIBIODgwBGb42WCuhyuZBIJFBSUoKhoSEAgMfjQSgUAgD09fUBAEKhEDweDwBgaGgI4XAYiURi7Ly7YDA4dsu1wcFBRCIRxONx9Pf3o7y8HIFAAD6fDwAwMDCAaDSKWCwGt9sNAPD7/fD7/ZbOKRwOz3lO+fn5aecEYM5zstvtpHMCkNH7ZMWcrHqfysvLLVl76eZUWlpKNqdoNGrZ2ks3pxQU79OCBQtI/zzNmzePbO2Vl5eTfUaksGLtpZtT6i4wFH+e8vLyLP0snzynRCKhzWf55DmVlZWLine truncated
"text/plain": [
"<Figure size 864x720 with 1 Axes>"
]
@@ -249,6 +249,17 @@
"show(pct2, mph2, xlabel='Segment Grade (percent)') "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"So I can cruise at 17 mph on a 1% grade, but only a bit over 6 mph at 6% grade, and around 5 mph on 8% grades.\n",
"\n",
"# VAM\n",
"\n",
"The number of vertical feet climbed per hour (regardless of horizontal distance travelled) is known as [VAM](https://en.wikipedia.org/wiki/VAM_(bicycling)), which stands for *velocità ascensionale media* (when you have Campagnolo components) or *mean ascent velocity* (on SRAM) or 平均上昇率 (on Shimano). Legend has it that at grades between about 6% to 12%, most of your power is going into lifting against gravity, so your VAM should be about constant no matter what the grade in that range. (I guess that below that grade, power is spent on wind and rolling resistance, and above it, power goes largely to cursing *sotto voce*, as they say in Italian.) Here's a plot of VAM versus grades (in percent):"
]
},
{
"cell_type": "code",
"execution_count": 9,
@@ -256,7 +267,7 @@
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAzoAAAKSCAYAAAAaia6DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsnX98U/W9/1+lhSLlRxmU/phQi3ZMUcfdJl7xC1cFdSpgvdqLypxzu1Vx3rn5c7iIMruryH64+QMdTt3m5lxRq6KO66Y4cU5EHygTx+3VkgJNQhuahLQkIUm/f7BEWtKStHl/0tc57+fj0YeQnLxy3s/zzrFvTs45eT6frxuKoiiKoiiKoigWYliuV0BRFEVRFEVRFCXb6KCjKIqiKIqiKIrl0EFHURRFURRFURTLoYOOoiiKoiiKoiiWQwcdRVEURVEURVEshw46iqIoiqIoiqJYDh10FEVRFAXACSecgBNOOCHt5Z1OJ4qLi7FkyRLBtVIURVEGig46iqIoSsY899xzuOmmm3DOOedg8uTJKC4uxje+8Y0+l29tbcUvfvEL/Md//Ae+8IUvoLS0FJWVlTj33HPxu9/9Dt3dh97SbcmSJSguLu73Z+HChZJlKoqiKMQU5HoFFEVRFD5WrlyJv//97xg9ejQqKiqwd+/efpf/xS9+gXvvvRdHHnkk/t//+3+oqKiAy+XC2rVrcc011+B//ud/8Pjjj/d4zXnnnYcpU6akzHvqqaewfft2nHnmmdkqCc8//3zWshRFUZTck+fz+Q79ZzRFURRF6Ye//OUv+OxnP4upU6diw4YNWLBgAf793/8djz76aMrln3/+eRQXF2P27NnIy8tLPu5yuTBv3jzs2rULv/71r9M6QuPz+XDsscciFovho48+woQJE7JWVyY4nU584QtfwCWXXIJVq1blZB0URVGUvtGvrimKoigZM2fOHBx99NE9hpb+WLhwIebMmXPI8uXl5cmvvL3xxhtpZT311FPYt28fFixYkPaQ88YbbyTPp2lqasLll1+Oo48+GuPHj8cHH3wAoO9zdPbu3Ytbb70Vxx13HEpLS3HSSSfh/vvvT/l1uwShUAj33Xcf/u3f/g2f/exnUVFRgdNOOw2PPvpoyte9/PLLOP/88/H5z38ekyZNwrRp03D22Wfjxz/+cVr1KYqiKIeiX11TFEVRcsrw4cMBAAUF6f0v6Ve/+hUA4PLLL8/4vZqbmzF37lxMmzYNixYtwt69ezFq1Kg+lw+Hwzj//PPx3nvv4bjjjkNtbS0CgQB+9KMf4c0330z5mr1796KmpgbvvvsuTjzxRFx66aUAgD//+c+4/vrr8c477/Q4AvT444/jO9/5DiZNmoSzzz4bJSUl8Hq92LZtGx577DHccMMNGdepKIqi6KCjKIqi5JBoNIrf//73AIB58+YddvmNGzdi69atOOaYYzBnzpyM3+9vf/sbbrjhBtx2221pLX///ffjvffew7nnnosnnngCw4Yd+CLEd7/7XZx22mkpX3Prrbfi3XffxR133IHvfOc7ycfD4TAuu+wyPPnkk1i4cCHOOeccAMBjjz2GESNGYMOGDZg0aVKPLK/Xm3GNiqIoygH0q2uKoihKznA4HNi6dSvOPvtszJ0797DLJy5YMJCjOQAwadIk3HzzzWkv/9vf/hZ5eXlYvnx5csgBgClTpuCqq646ZPmOjg48+eSTOPHEE3sMOQBQWFiIZcuWATjw9buDKSgoSB7ZOphcnX+kKIpiBfSIjqIoipITfvKTn+Chhx7Csccei4cffviwy/v9fjQ2NmLEiBHJr4NlyvHHH4/CwsK0lt27dy8++eQTlJWVobq6+pDnTz311EMee/fddxGNRjFs2DDcddddhzwfjUYBAP/7v/+bfKy2thYOhwMnn3wyLrjgAsyaNQsnn3wyysrK0i1LURRFSYEOOoqiKIpxfvSjH6G+vh7Tp0/Hc889h+Li4sO+5g9/+AO6urpw4YUXDvhIR++vhvVHIBAAAJSUlKSdtWfPHgDA5s2bsXnz5j6zg8Fg8s/XXnstSkpK8Oijj+KRRx5JDn0nnXQSli1bhtmzZ6e9zoqiKMqn6FfXFEVRFKPU19ejvr4eM2bMwNq1azFx4sS0Xpe4CMHXv/71Ab93uleJA4CxY8cCANra2lI+v3v37j5fc+WVV8Ln8/X5k7jSW4JFixZh3bp1aG5uxjPPPIO6ujq8//77qK2txf/93/+lvc6KoijKp+igoyiKohjj1ltvxY9+9COcdNJJeO655zB+/Pi0Xrdp0yb8/e9/xzHHHGPsCMeYMWMwdepUeDyelMNGqquuffnLX8awYcPw1ltvDeg9x44dizPOOAMrV67Etddei1AohFdeeWVAWYqiKHZHBx1FURRFnO7ublx//fV48MEHMWvWLDz77LMYN25c2q8f7EUIBsrixYvR3d2NZcuWIR6PJx9vaWlJeV7RxIkTsWjRImzZsgV33XVX8pycg9m1a1ePc3Ref/31lPfW8Xg8ANDv5a8VRVGUvtFzdBRFUZSMWbt2LV588UUAn36Fa9OmTViyZElymYPvFbNixQo8+uijGDFiBL785S/jvvvuOyRzypQpWLx48SGPBwIBPPvssygsLBzwRQgGyrXXXosXX3wRL730EmbPno158+Yl1+eUU07Byy+/fMhr7rnnHnzyySdYsWIFnnrqKcyaNQulpaXJI0PvvPMOfvjDH+Jzn/scAOCyyy5DUVERvvzlL2PKlCnIy8vDu+++i7feegtVVVWoqakxWrOiKIpV0EFHURRFyZgtW7bgySef7PFYS0sLWlpakn8/eNBxOp0AgEgkgp///OcpM0899dSUg05DQwM6OzsHdRGCgVJYWIjGxkbcfffdePbZZ/HQQw9hypQpuOGGG7BgwYKUg86YMWOwdu1a/OY3v0FDQwPWrl2LUCiEkpISVFZW4vbbb8cFF1yQXP6OO+7Aq6++ii1btuDPf/4zCgoKcOSRR+KWW27BVVddldGRL0VRFOVT8nw+36HHyxVFURRFURRFUYjRc3QURVEURVEURbEcOugoiqIoiqIoimI5dNBRFEVRFEVRFMVy6KCjKIqiKIqiKIrl0EFHURRFURRFURTLoYOOoiiKoiiKoiiWQwcdRVEURVEURVEshw46wng8Htr8pqYmsWxmL9LrzuqdeZtKOgfUe1+w9rp0vu5jzGdL5+s+xnw2wNvr0vl28q6DjjDSd/E2fZfwbMHshdU5wO1FvZvPNpEvBbMXVucAtxf1bq1saZj70U7eddARJh6PU+dLweyF1TnA7UW9m882kS8FsxdW5wC3F/VurWxpmPvRTt5zNuisXr0as2bNwuTJkzF58mSceeaZWLduXfL5JUuWoLi4uMfPvHnzemSEw2HcdNNNmDp1KioqKnDxxRdj165dPZbZsWMHFi1ahIqKCkydOhU333wzIpGIkRoBwO/3U+dLweyF1TnA7UW9m882kS8FsxdW5wC3F/VurWxpmPvRTt4LhNbjsFRUVGD58uU4+uijEY/H8eSTT2Lx4sVYv349jj/+eADAaaedhocffjj5mhEjRvTIWLp0KV566SX88pe/xPjx4/H9738fixYtwuuvv478/HzEYjEsWrQI48ePx0svvYSOjg4sWbIE3d3dWLlypZE6S0pKqPOlYPbC6hzg9qLezWebyJeC2Qurc4Dbi3q3VrY0zP1oJ+85O6Jz3nnn4cwzz8TUqVNxzDHH4LbbbsPo0aPxzjvvJJcpLCxEaWlp8mf8+PHJ5/x+P37zm9/gBz/4AU4//XTMmDEDDz/8MD788EOsX78eAPDqq6/io48+wsMPP4wZM2bg9NNPx/Lly/HrX/8agUDASJ3ME78kzF5YnQPcXtS7+WwT+VIwe2F1DnB7Ue/WypaGuR/t5H1InKMTi8Xw9NNPo7OzEzNnzkw+/tZbb+GYY47Bl770JXz7299GW1tb8rnNmzdj//79OOOMM5KPHXnkkZg2bRrefvttAMDGjRsxbdo0HHnkkcll5s6di3A4jM2bNxuoDMjPz6fOl4LZC6tzgNuLejefbSJfCmYvrM4Bbi/q3VrZ0jD3o5285+yrawDw4Ycf4qyzzkIoFEJRURGeeOIJTJ8+HQAwb948LFiwAJWVlWhpaUF9fT0WLlyI9evXo7CwELt370Z+fv4hV18oKSnB7t27AQC7d+8+5BDXhAkTkJ+fn1xGmtGjR1PnS8HshdU5wO1FvZvPNpEvBbMXVucAtxf1bq1saZj70U7eczroVFdX44033kAgEMBzzz2HJUuWYO3atTjuuONw4YUXJpebPn06ZsyYgRNOOAHr1q3DwoULxddN+pr3LKiH3KDezaPOc4N6zw3q3TzqPDeo99xg0nt1dXWfz+V00BkxYgSmTp0KAJgxYwbee+89PPjgg7j//vsPWba8vBwVFRX45JNPAACTJk1CLBaD1+vFxIkTk8u1tbXhlFNOSS6T+BpbAq/Xi1gshkmTJvW7bv1Jy4RYLCZ6iFAyv6mpKWseesPsRXrdWb0zb1NJ54B67wvWXpfO132M+WzpfN3HmM8GeHtdOt9O3ofEOToJ4vF4n5d+9nq9cLlcKC0tBXBgMBo+fDhee+215DK7du3Ctm3bcPLJJwMAZs6ciW3btvW45PRrr72GwsJCzJgxQ7CST4lGo9T5UjB7YXUOcHtR7+azTeRLweyF1TnA7UW9WytbGuZ+tJP3nA06d9xxB/7617/C6XTiww8/xPLly7FhwwbU1tYiGAzC4XBg48aNcDqdeOONN3DxxRejpKQE8+fPBwCMGzcOl112GW6//XasX78e77//Pq666ipMnz4dp512GgDgjDPOwLHHHourr74a77//PtavX49ly5bha1/7GsaOHWukzr1791LnS8HshdU5wO1FvZvPNpEvBbMXVucAtxf1bq1saZj70U7ec/bVNY/HgyuvvBKLine truncated
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAzoAAAKSCAYAAAAaia6DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsnX98U/W9/1+lhSLlRxmU/phQi3ZMUcfdJl7xC1cFdSpgvdqLypxzu1Vx3rn5c7iIMruryH64+QMdTt3m5lxRq6KO66Y4cU5EHygTx+3VkgJNQhuahLQkIUm/f7BEWtKStHl/0tc57+fj0YeQnLxy3s/zzrFvTs45eT6frxuKoiiKoiiKoigWYliuV0BRFEVRFEVRFCXb6KCjKIqiKIqiKIrl0EFHURRFURRFURTLoYOOoiiKoiiKoiiWQwcdRVEURVEURVEshw46iqIoiqIoiqJYDh10FEVRFAXACSecgBNOOCHt5Z1OJ4qLi7FkyRLBtVIURVEGig46iqIoSsY899xzuOmmm3DOOedg8uTJKC4uxje+8Y0+l29tbcUvfvEL/Md//Ae+8IUvoLS0FJWVlTj33HPxu9/9Dt3dh97SbcmSJSguLu73Z+HChZJlKoqiKMQU5HoFFEVRFD5WrlyJv//97xg9ejQqKiqwd+/efpf/xS9+gXvvvRdHHnkk/t//+3+oqKiAy+XC2rVrcc011+B//ud/8Pjjj/d4zXnnnYcpU6akzHvqqaewfft2nHnmmdkqCc8//3zWshRFUZTck+fz+Q79ZzRFURRF6Ye//OUv+OxnP4upU6diw4YNWLBgAf793/8djz76aMrln3/+eRQXF2P27NnIy8tLPu5yuTBv3jzs2rULv/71r9M6QuPz+XDsscciFovho48+woQJE7JWVyY4nU584QtfwCWXXIJVq1blZB0URVGUvtGvrimKoigZM2fOHBx99NE9hpb+WLhwIebMmXPI8uXl5cmvvL3xxhtpZT311FPYt28fFixYkPaQ88YbbyTPp2lqasLll1+Oo48+GuPHj8cHH3wAoO9zdPbu3Ytbb70Vxx13HEpLS3HSSSfh/vvvT/l1uwShUAj33Xcf/u3f/g2f/exnUVFRgdNOOw2PPvpoyte9/PLLOP/88/H5z38ekyZNwrRp03D22Wfjxz/+cVr1KYqiKIeiX11TFEVRcsrw4cMBAAUF6f0v6Ve/+hUA4PLLL8/4vZqbmzF37lxMmzYNixYtwt69ezFq1Kg+lw+Hwzj//PPx3nvv4bjjjkNtbS0CgQB+9KMf4c0330z5mr1796KmpgbvvvsuTjzxRFx66aUAgD//+c+4/vrr8c477/Q4AvT444/jO9/5DiZNmoSzzz4bJSUl8Hq92LZtGx577DHccMMNGdepKIqi6KCjKIqi5JBoNIrf//73AIB58+YddvmNGzdi69atOOaYYzBnzpyM3+9vf/sbbrjhBtx2221pLX///ffjvffew7nnnosnnngCw4Yd+CLEd7/7XZx22mkpX3Prrbfi3XffxR133IHvfOc7ycfD4TAuu+wyPPnkk1i4cCHOOeccAMBjjz2GESNGYMOGDZg0aVKPLK/Xm3GNiqIoygH0q2uKoihKznA4HNi6dSvOPvtszJ0797DLJy5YMJCjOQAwadIk3HzzzWkv/9vf/hZ5eXlYvnx5csgBgClTpuCqq646ZPmOjg48+eSTOPHEE3sMOQBQWFiIZcuWATjw9buDKSgoSB7ZOphcnX+kKIpiBfSIjqIoipITfvKTn+Chhx7Csccei4cffviwy/v9fjQ2NmLEiBHJr4NlyvHHH4/CwsK0lt27dy8++eQTlJWVobq6+pDnTz311EMee/fddxGNRjFs2DDcddddhzwfjUYBAP/7v/+bfKy2thYOhwMnn3wyLrjgAsyaNQsnn3wyysrK0i1LURRFSYEOOoqiKIpxfvSjH6G+vh7Tp0/Hc889h+Li4sO+5g9/+AO6urpw4YUXDvhIR++vhvVHIBAAAJSUlKSdtWfPHgDA5s2bsXnz5j6zg8Fg8s/XXnstSkpK8Oijj+KRRx5JDn0nnXQSli1bhtmzZ6e9zoqiKMqn6FfXFEVRFKPU19ejvr4eM2bMwNq1azFx4sS0Xpe4CMHXv/71Ab93uleJA4CxY8cCANra2lI+v3v37j5fc+WVV8Ln8/X5k7jSW4JFixZh3bp1aG5uxjPPPIO6ujq8//77qK2txf/93/+lvc6KoijKp+igoyiKohjj1ltvxY9+9COcdNJJeO655zB+/Pi0Xrdp0yb8/e9/xzHHHGPsCMeYMWMwdepUeDyelMNGqquuffnLX8awYcPw1ltvDeg9x44dizPOOAMrV67Etddei1AohFdeeWVAWYqiKHZHBx1FURRFnO7ublx//fV48MEHMWvWLDz77LMYN25c2q8f7EUIBsrixYvR3d2NZcuWIR6PJx9vaWlJeV7RxIkTsWjRImzZsgV33XVX8pycg9m1a1ePc3Ref/31lPfW8Xg8ANDv5a8VRVGUvtFzdBRFUZSMWbt2LV588UUAn36Fa9OmTViyZElymYPvFbNixQo8+uijGDFiBL785S/jvvvuOyRzypQpWLx48SGPBwIBPPvssygsLBzwRQgGyrXXXosXX3wRL730EmbPno158+Yl1+eUU07Byy+/fMhr7rnnHnzyySdYsWIFnnrqKcyaNQulpaXJI0PvvPMOfvjDH+Jzn/scAOCyyy5DUVERvvzlL2PKlCnIy8vDu+++i7feegtVVVWoqakxWrOiKIpV0EFHURRFyZgtW7bgySef7PFYS0sLWlpakn8/eNBxOp0AgEgkgp///OcpM0899dSUg05DQwM6OzsHdRGCgVJYWIjGxkbcfffdePbZZ/HQQw9hypQpuOGGG7BgwYKUg86YMWOwdu1a/OY3v0FDQwPWrl2LUCiEkpISVFZW4vbbb8cFF1yQXP6OO+7Aq6++ii1btuDPf/4zCgoKcOSRR+KWW27BVVddldGRL0VRFOVT8nw+36HHyxVFURRFURRFUYjRc3QURVEURVEURbEcOugoiqIoiqIoimI5dNBRFEVRFEVRFMVy6KCjKIqiKIqiKIrl0EFHURRFURRFURTLoYOOoiiKoiiKoiiWQwcdRVEURVEURVEshw46wng8Htr8pqYmsWxmL9LrzuqdeZtKOgfUe1+w9rp0vu5jzGdL5+s+xnw2wNvr0vl28q6DjjDSd/E2fZfwbMHshdU5wO1FvZvPNpEvBbMXVucAtxf1bq1saZj70U7eddARJh6PU+dLweyF1TnA7UW9m882kS8FsxdW5wC3F/VurWxpmPvRTt5zNuisXr0as2bNwuTJkzF58mSceeaZWLduXfL5JUuWoLi4uMfPvHnzemSEw2HcdNNNmDp1KioqKnDxxRdj165dPZbZsWMHFi1ahIqKCkydOhU333wzIpGIkRoBwO/3U+dLweyF1TnA7UW9m882kS8FsxdW5wC3F/VurWxpmPvRTt4LhNbjsFRUVGD58uU4+uijEY/H8eSTT2Lx4sVYv349jj/+eADAaaedhocffjj5mhEjRvTIWLp0KV566SX88pe/xPjx4/H9738fixYtwuuvv478/HzEYjEsWrQI48ePx0svvYSOjg4sWbIE3d3dWLlypZE6S0pKqPOlYPbC6hzg9qLezWebyJeC2Qurc4Dbi3q3VrY0zP1oJ+85O6Jz3nnn4cwzz8TUqVNxzDHH4LbbbsPo0aPxzjvvJJcpLCxEaWlp8mf8+PHJ5/x+P37zm9/gBz/4AU4//XTMmDEDDz/8MD788EOsX78eAPDqq6/io48+wsMPP4wZM2bg9NNPx/Lly/HrX/8agUDASJ3ME78kzF5YnQPcXtS7+WwT+VIwe2F1DnB7Ue/WypaGuR/t5H1InKMTi8Xw9NNPo7OzEzNnzkw+/tZbb+GYY47Bl770JXz7299GW1tb8rnNmzdj//79OOOMM5KPHXnkkZg2bRrefvttAMDGjRsxbdo0HHnkkcll5s6di3A4jM2bNxuoDMjPz6fOl4LZC6tzgNuLejefbSJfCmYvrM4Bbi/q3VrZ0jD3o5285+yrawDw4Ycf4qyzzkIoFEJRURGeeOIJTJ8+HQAwb948LFiwAJWVlWhpaUF9fT0WLlyI9evXo7CwELt370Z+fv4hV18oKSnB7t27AQC7d+8+5BDXhAkTkJ+fn1xGmtGjR1PnS8HshdU5wO1FvZvPNpEvBbMXVucAtxf1bq1saZj70U7eczroVFdX44033kAgEMBzzz2HJUuWYO3atTjuuONw4YUXJpebPn06ZsyYgRNOOAHr1q3DwoULxddN+pr3LKiH3KDezaPOc4N6zw3q3TzqPDeo99xg0nt1dXWfz+V00BkxYgSmTp0KAJgxYwbee+89PPjgg7j//vsPWba8vBwVFRX45JNPAACTJk1CLBaD1+vFxIkTk8u1tbXhlFNOSS6T+BpbAq/Xi1gshkmTJvW7bv1Jy4RYLCZ6iFAyv6mpKWseesPsRXrdWb0zb1NJ54B67wvWXpfO132M+WzpfN3HmM8GeHtdOt9O3ofEOToJ4vF4n5d+9nq9cLlcKC0tBXBgMBo+fDhee+215DK7du3Ctm3bcPLJJwMAZs6ciW3btvW45PRrr72GwsJCzJgxQ7CST4lGo9T5UjB7YXUOcHtR7+azTeRLweyF1TnA7UW9WytbGuZ+tJP3nA06d9xxB/7617/C6XTiww8/xPLly7FhwwbU1tYiGAzC4XBg48aNcDqdeOONN3DxxRejpKQE8+fPBwCMGzcOl112GW6//XasX78e77//Pq666ipMnz4dp512GgDgjDPOwLHHHourr74a77//PtavX49ly5bha1/7GsaOHWukzr1791LnS8HshdU5wO1FvZvPNpEvBbMXVucAtxf1bq1saZj70U7ec/bVNY/HgyuvvBKLine truncated
"text/plain": [
"<Figure size 864x720 with 1 Axes>"
]
@@ -288,7 +299,7 @@
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAzoAAAKSCAYAAAAaia6DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsvXt4VNW9//8OAQIESGgIuQiEoBQVsZxW8VusFgXvgPEoxUvVak+seDzWO1UjaDs9ilhra73iraf2tIjVKKi1toIVa0X0oV5QmnOMwyXJJBkymWTCzDAz+f3Bb+Yk5sKePfNZe6/F+/U8Pk8bxtd81mtWYhYzsycnEAh0gxBCCCGEEEIMYojTAxBCCCGEEEJItuFBhxBCCCGEEGIcPOgQQgghhBBCjIMHHUIIIYQQQohx8KBDCCGEEEIIMQ4edAghhBBCCCHGwYMOIYQQAmDmzJmYOXOm5dt7vV4UFhZi6dKlglMRQgixCw86hBBC0ubFF1/ETTfdhDPOOAOTJk1CYWEhLr/88gFv39DQgMceewzf+c538LWvfQ0lJSWoqKjAmWeeif/+7/9Gd3ffj3RbunQpCgsLB/1n0aJFksskhBCiMUOdHoAQQoh+rFq1Ch9//DFGjx6N8vJydHR0DHr7xx57DPfffz8mTpyIb33rWygvL0djYyPWr1+Pq666Cn/605/w9NNP9/p3zjrrLEyePLlf35o1a/DFF1/glFNOydaS8NJLL2XNRQghxHlyAoFA379GI4QQQgbhr3/9Kw455BBMnToVmzZtwsKFC/Gv//qvePLJJ/u9/UsvvYTCwkKccMIJyMnJSX29sbER8+fPx+7du/Ff//Vflp6hCQQCOOKIIxCPx/Hpp5+iqKgoa+tKB6/Xi6997Wu44IIL8PDDDzsyAyGEkIHhS9cIIYSkzYknnohDDz2016FlMBYtWoQTTzyxz+3LyspSL3l76623LLnWrFmDvXv3YuHChZYPOW+99Vbq/TR1dXW49NJLceihh2LcuHH48MMPAQz8Hp2Ojg7ceuutOPLII1FSUoJjjz0Wv/rVr/p9uV2ScDiMBx54AN/+9rdxyCGHoLy8HHPnzsWTTz7Z77/36quv4uyzz8bhhx+OCRMmYPr06TjttNPws5/9zNL6CCGE9IUvXSOEEOIow4YNAwAMHWrtP0m//vWvAQCXXnpp2vdVX1+PefPmYfr06ViyZAk6OjowatSoAW8fiURw9tln44MPPsCRRx6JxYsXIxgM4t5778Xbb7/d77/T0dGBqqoqvP/++zj66KNx4YUXAgD+8pe/4Prrr8d7773X6xmgp59+Gtdeey0mTJiA0047DcXFxfD7/di+fTueeuop3HDDDWmvkxBCCA86hBBCHCQWi+H3v/89AGD+/PkHvP3mzZuxbds2HHbYYTjxxBPTvr+///3vuOGGG3D77bdbuv2vfvUrfPDBBzjzzDPxzDPPYMiQ/S+EuO666zB37tx+/51bb70V77//Pu644w5ce+21qa9HIhFcfPHF+N3vfodFixbhjDPOAAA89dRTGD58ODZt2oQJEyb0cvn9/rTXSAghZD986RohhBDHqKmpwbZt23Daaadh3rx5B7x98oIFdp7NAYAJEybg5ptvtnz73/72t8jJycGdd96ZOuQAwOTJk/GDH/ygz+3b2trwu9/9DkcffXSvQw4A5OXlYfny5QD2v/yuJ0OHDk09s9UTp95/RAghJsBndAghhDjCfffdh0ceeQRHHHEEHn300QPevr29HbW1tRg+fHjq5WDpctRRRyEvL8/SbTs6OvD555+jtLQU06ZN6/Pnxx9/fJ+vvf/++4jFYhgyZAjuuuuuPn8ei8UAAP/85z9TX1u8eDFqampw3HHH4ZxzzsGcOXNw3HHHobS01OqyCCGE9AMPOoQQQpRz7733wuPxYMaMGXjxxRdRWFh4wH/n2WefRVdXF84991zbz3R8+aVhgxEMBgEAxcXFll179uwBAGzduhVbt24d0N3Z2Zn631dffTWKi4vx5JNP4vHHH08d+o499lgsX74cJ5xwguWZCSGE/B986RohhBCleDweeDwezJo1C+vXr8f48eMt/XvJixB873vfs33fVq8SBwBjx44FALS0tPT7583NzQP+O1dccQUCgcCA/ySv9JZkyZIleO2111BfX4/nn38e1dXV+Mc//oHFixfjf/7nfyzPTAgh5P/gQYcQQogybr31Vtx777049thj8eKLL2LcuHGW/r0tW7bg448/xmGHHabsGY4xY8Zg6tSp8Pl8/R42+rvq2jHHHIMhQ4bgnXfesXWfY8eOxcknn4xVq1bh6quvRjgcxuuvv27LRQghBzs86BBCCBGnu7sb119/PR566CHMmTMHL7zwAgoKCiz/+5lehMAuF110Ebq7u7F8+XIkEonU13fs2NHv+4rGjx+PJUuW4KOPPsJdd92Vek9OT3bv3t3rPTpvvvlmv5+t4/P5AGDQy18TQggZGL5HhxBCSNqsX78eL7/8MoD/ewnXli1bsHTp0tRten5WzMqVK/Hkk09i+PDhOOaYY/DAAw/0cU6ePBkXXXRRn68Hg0G88MILyMvLs30RArtcffXVePnll/HKK6/ghBNOwPz581PzfPOb38Srr77a59+555578Pnnn2PlypVYs2YN5syZg5KSktQzQ++99x5++tOf4qtf/SoA4OKLL0Z+fj6OOeYYTJ48GTk5OXj//ffxzjvvoLKyElVVVUrXTAghpsCDDiGEkLT56KOP8Lvf/a7X13bs2IEdO3ak/n/Pg47X6wUARKNR/PKXv+zXefzxx/d70Fm7di1CoVBGFyGwS15eHmpra3H33XfjhRdewCOPPILJkyfjhhtuwMKFC/s96IwZMwbr16/Hb37zG6xduxbr169HOBxGcXExKioqsGLFCpxzzjmp299xxx1444038NFHH+Evf/kLhg4diokTJ2LZsmX4wQ9+kNYzX4QQQv6PnEAg0Pf5ckIIIYQQQgjRGL5HhxBCCCGEEGIcPOgQQgghhBBCjIMHHUIIIYQQQohx8KBDCCGEEEIIMQ4edAghhBBCCCHGwYMOIYQQQgghxDh40CGEEEIIIYQYBw86FvH5fGLuuro6Mbfk3JJuaT+bq3Wzt3o/m6t1s7d6P5urdevaW9rP5u5286BjEdWfxp0tJOeWbsLm5rgl0bkJm5vjlkTnJmxujlsSnfchm7vbzYOORRKJhNMj2EJybukmbG6OWxKdm7C5OW5JdG7C5ua4JdF5H7K5u9086Fikvb3d6RFsITm3dBM2N8ctic5N2NwctyQ6N2Fzc9yS6LwP2dzdbscOOqtXr8acOXMwadIkTJo0Caeccgpee+211J8vXboUhYWFvf6ZP39+L0ckEsFNN92EqVOnory8HOeffz52797d6zY7d+7EkiVLUF5ejqlTp+Lmm29GNBpNe97i4mJ7C3UYybmlm7C5OW5JdG7C5ua4JdG5CZub45ZE533I5u52O3bQKS8vx5133ok333wTGzZswIknnoiLLroIH3/8ceo2c+fOxfbt21P/rF27tpfjlltuwbp16/DEE0/glVdeQUdHB5YsWYJ4PA4AiMfjWLJkCTo7O/HKK6/giSeewEsvvYTbbrst7Xl5YlfrVuGXQtfm7K3WrcIvha7N2VutW4VfCl2bs7d6P5u72z1UaI4DctZZZ/X6/7fffjueeOIJvPfeezjqqKMAAHl5eSgpKen3329vb8dvfvMbPPjggzjppJMAAI8++ihmzpyJjRs3Yt68eXjjjTfw6aef4qOPPsLEiRMBAHfeeSeuueYa3H777Rg7dqzleXNzc+0s03Ek55ZuwubmuCXRuQmbm+OWROcmbG6OWxKd9yGbu9vtivfoxONx/OEPf0AoFMLs2bNTX3/nnXdw2GGH4Rvf+AauueYatLS0pP5s69at2LdvH04++eTU1yZOnIjp06fj3XffBQBs3rwZ06dPTx1yAGDevHmIRCLYunVrWjOOHj3a7vIcRXJu6SZsbo5bEp2bsLk5bkl0bsLm5rgl0Xkfsrm73Y49owMAn3zyCU499VSEw2Hk5+fjmWeewYwZMwAA8+fPx8KFC1FRUYEdO3bA4/Fg0aJF2LhxI/Ly8tDc3Izc3Nw+l5krLi5Gc3MzAKC5ubnPa/mKioqQm5ubus1ASF673A33R9hcNeytHjZXC3urh83Vwt7qYfMDM23atAH/zNGDzrRp0/DWW28hGAzixRdfxNKlS7F+/XoceeSROPfcc1O3mzFjBmbNmoWZM2fitddew6JFi5TM1pN4PC72VFxdXd2gD1ImSM4t6Zb2s7laN3ur97O5Wjd7q/ezuVq3rr2l/WzubrejL10bPnw4pk6dilmzZmHFihWYOXMmHnrooX5vW1ZWhvLycnz++ecAgAkTJiAej8Pv9/e6XUtLCyZMmJC6Tc+XuwGA3+9HPB5P3cYqsVgsrdu7Bcm5pZuwuTluSXRuwubmuCXRuQmbm+OWROd9yObudrviPTpJEonEgJd+9vv9aGxsTF2cYNasWRg2bBg2bNiQus3u3buxfft2HHfccQCA2bNnY/v27b0uOb1hwwbk5eVh1qxZac3W0dGR7nJcgeTc0k3Y3By3JDo3YXNz3JLo3ITNzXFLovM+ZHN3ux176dodd9yBU089FYcccgg6Ozvx3HPPYdOmTXj22WfR2dmJu+++G4sWLUJJSQl27NiBH//4xyguLsaCBQsAAAUFBbj44ouxYsUKFBcXY9y4cbjtttswY8YMzJ07FwBw8skn44gjjsCVV14Jj8eDtrY2LF++HJdccklaV1wDgPHjx2c7gRIk55ZuwubmuCXRuQmbm+OWROcmbG6OWxKd9yGbu9vt2DM6Pp8PV1xxBY499licffbZ+OCDD/Dcc8/hlFNOQW5uLrZt24YLL7wQxxxzDJYuXYrDDjsMf/rTnzBmzJiU46677sJLine truncated
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAzoAAAKSCAYAAAAaia6DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsvXt4VNW9//8OAQIESGgIuQiEoBQVsZxW8VusFgXvgPEoxUvVak+seDzWO1UjaDs9ilhra73iraf2tIjVKKi1toIVa0X0oV5QmnOMwyXJJBkymWTCzDAz+f3Bb+Yk5sKePfNZe6/F+/U8Pk8bxtd81mtWYhYzsycnEAh0gxBCCCGEEEIMYojTAxBCCCGEEEJItuFBhxBCCCGEEGIcPOgQQgghhBBCjIMHHUIIIYQQQohx8KBDCCGEEEIIMQ4edAghhBBCCCHGwYMOIYQQAmDmzJmYOXOm5dt7vV4UFhZi6dKlglMRQgixCw86hBBC0ubFF1/ETTfdhDPOOAOTJk1CYWEhLr/88gFv39DQgMceewzf+c538LWvfQ0lJSWoqKjAmWeeif/+7/9Gd3ffj3RbunQpCgsLB/1n0aJFksskhBCiMUOdHoAQQoh+rFq1Ch9//DFGjx6N8vJydHR0DHr7xx57DPfffz8mTpyIb33rWygvL0djYyPWr1+Pq666Cn/605/w9NNP9/p3zjrrLEyePLlf35o1a/DFF1/glFNOydaS8NJLL2XNRQghxHlyAoFA379GI4QQQgbhr3/9Kw455BBMnToVmzZtwsKFC/Gv//qvePLJJ/u9/UsvvYTCwkKccMIJyMnJSX29sbER8+fPx+7du/Ff//Vflp6hCQQCOOKIIxCPx/Hpp5+iqKgoa+tKB6/Xi6997Wu44IIL8PDDDzsyAyGEkIHhS9cIIYSkzYknnohDDz2016FlMBYtWoQTTzyxz+3LyspSL3l76623LLnWrFmDvXv3YuHChZYPOW+99Vbq/TR1dXW49NJLceihh2LcuHH48MMPAQz8Hp2Ojg7ceuutOPLII1FSUoJjjz0Wv/rVr/p9uV2ScDiMBx54AN/+9rdxyCGHoLy8HHPnzsWTTz7Z77/36quv4uyzz8bhhx+OCRMmYPr06TjttNPws5/9zNL6CCGE9IUvXSOEEOIow4YNAwAMHWrtP0m//vWvAQCXXnpp2vdVX1+PefPmYfr06ViyZAk6OjowatSoAW8fiURw9tln44MPPsCRRx6JxYsXIxgM4t5778Xbb7/d77/T0dGBqqoqvP/++zj66KNx4YUXAgD+8pe/4Prrr8d7773X6xmgp59+Gtdeey0mTJiA0047DcXFxfD7/di+fTueeuop3HDDDWmvkxBCCA86hBBCHCQWi+H3v/89AGD+/PkHvP3mzZuxbds2HHbYYTjxxBPTvr+///3vuOGGG3D77bdbuv2vfvUrfPDBBzjzzDPxzDPPYMiQ/S+EuO666zB37tx+/51bb70V77//Pu644w5ce+21qa9HIhFcfPHF+N3vfodFixbhjDPOAAA89dRTGD58ODZt2oQJEyb0cvn9/rTXSAghZD986RohhBDHqKmpwbZt23Daaadh3rx5B7x98oIFdp7NAYAJEybg5ptvtnz73/72t8jJycGdd96ZOuQAwOTJk/GDH/ygz+3b2trwu9/9DkcffXSvQw4A5OXlYfny5QD2v/yuJ0OHDk09s9UTp95/RAghJsBndAghhDjCfffdh0ceeQRHHHEEHn300QPevr29HbW1tRg+fHjq5WDpctRRRyEvL8/SbTs6OvD555+jtLQU06ZN6/Pnxx9/fJ+vvf/++4jFYhgyZAjuuuuuPn8ei8UAAP/85z9TX1u8eDFqampw3HHH4ZxzzsGcOXNw3HHHobS01OqyCCGE9AMPOoQQQpRz7733wuPxYMaMGXjxxRdRWFh4wH/n2WefRVdXF84991zbz3R8+aVhgxEMBgEAxcXFll179uwBAGzduhVbt24d0N3Z2Zn631dffTWKi4vx5JNP4vHHH08d+o499lgsX74cJ5xwguWZCSGE/B986RohhBCleDweeDwezJo1C+vXr8f48eMt/XvJixB873vfs33fVq8SBwBjx44FALS0tPT7583NzQP+O1dccQUCgcCA/ySv9JZkyZIleO2111BfX4/nn38e1dXV+Mc//oHFixfjf/7nfyzPTAgh5P/gQYcQQogybr31Vtx777049thj8eKLL2LcuHGW/r0tW7bg448/xmGHHabsGY4xY8Zg6tSp8Pl8/R42+rvq2jHHHIMhQ4bgnXfesXWfY8eOxcknn4xVq1bh6quvRjgcxuuvv27LRQghBzs86BBCCBGnu7sb119/PR566CHMmTMHL7zwAgoKCiz/+5lehMAuF110Ebq7u7F8+XIkEonU13fs2NHv+4rGjx+PJUuW4KOPPsJdd92Vek9OT3bv3t3rPTpvvvlmv5+t4/P5AGDQy18TQggZGL5HhxBCSNqsX78eL7/8MoD/ewnXli1bsHTp0tRten5WzMqVK/Hkk09i+PDhOOaYY/DAAw/0cU6ePBkXXXRRn68Hg0G88MILyMvLs30RArtcffXVePnll/HKK6/ghBNOwPz581PzfPOb38Srr77a59+555578Pnnn2PlypVYs2YN5syZg5KSktQzQ++99x5++tOf4qtf/SoA4OKLL0Z+fj6OOeYYTJ48GTk5OXj//ffxzjvvoLKyElVVVUrXTAghpsCDDiGEkLT56KOP8Lvf/a7X13bs2IEdO3ak/n/Pg47X6wUARKNR/PKXv+zXefzxx/d70Fm7di1CoVBGFyGwS15eHmpra3H33XfjhRdewCOPPILJkyfjhhtuwMKFC/s96IwZMwbr16/Hb37zG6xduxbr169HOBxGcXExKioqsGLFCpxzzjmp299xxx1444038NFHH+Evf/kLhg4diokTJ2LZsmX4wQ9+kNYzX4QQQv6PnEAg0Pf5ckIIIYQQQgjRGL5HhxBCCCGEEGIcPOgQQgghhBBCjIMHHUIIIYQQQohx8KBDCCGEEEIIMQ4edAghhBBCCCHGwYMOIYQQQgghxDh40CGEEEIIIYQYBw86FvH5fGLuuro6Mbfk3JJuaT+bq3Wzt3o/m6t1s7d6P5urdevaW9rP5u5286BjEdWfxp0tJOeWbsLm5rgl0bkJm5vjlkTnJmxujlsSnfchm7vbzYOORRKJhNMj2EJybukmbG6OWxKdm7C5OW5JdG7C5ua4JdF5H7K5u9086Fikvb3d6RFsITm3dBM2N8ctic5N2NwctyQ6N2Fzc9yS6LwP2dzdbscOOqtXr8acOXMwadIkTJo0Caeccgpee+211J8vXboUhYWFvf6ZP39+L0ckEsFNN92EqVOnory8HOeffz52797d6zY7d+7EkiVLUF5ejqlTp+Lmm29GNBpNe97i4mJ7C3UYybmlm7C5OW5JdG7C5ua4JdG5CZub45ZE533I5u52O3bQKS8vx5133ok333wTGzZswIknnoiLLroIH3/8ceo2c+fOxfbt21P/rF27tpfjlltuwbp16/DEE0/glVdeQUdHB5YsWYJ4PA4AiMfjWLJkCTo7O/HKK6/giSeewEsvvYTbbrst7Xl5YlfrVuGXQtfm7K3WrcIvha7N2VutW4VfCl2bs7d6P5u72z1UaI4DctZZZ/X6/7fffjueeOIJvPfeezjqqKMAAHl5eSgpKen3329vb8dvfvMbPPjggzjppJMAAI8++ihmzpyJjRs3Yt68eXjjjTfw6aef4qOPPsLEiRMBAHfeeSeuueYa3H777Rg7dqzleXNzc+0s03Ek55ZuwubmuCXRuQmbm+OWROcmbG6OWxKd9yGbu9vtivfoxONx/OEPf0AoFMLs2bNTX3/nnXdw2GGH4Rvf+AauueYatLS0pP5s69at2LdvH04++eTU1yZOnIjp06fj3XffBQBs3rwZ06dPTx1yAGDevHmIRCLYunVrWjOOHj3a7vIcRXJu6SZsbo5bEp2bsLk5bkl0bsLm5rgl0Xkfsrm73Y49owMAn3zyCU499VSEw2Hk5+fjmWeewYwZMwAA8+fPx8KFC1FRUYEdO3bA4/Fg0aJF2LhxI/Ly8tDc3Izc3Nw+l5krLi5Gc3MzAKC5ubnPa/mKioqQm5ubus1ASF673A33R9hcNeytHjZXC3urh83Vwt7qYfMDM23atAH/zNGDzrRp0/DWW28hGAzixRdfxNKlS7F+/XoceeSROPfcc1O3mzFjBmbNmoWZM2fitddew6JFi5TM1pN4PC72VFxdXd2gD1ImSM4t6Zb2s7laN3ur97O5Wjd7q/ezuVq3rr2l/WzubrejL10bPnw4pk6dilmzZmHFihWYOXMmHnrooX5vW1ZWhvLycnz++ecAgAkTJiAej8Pv9/e6XUtLCyZMmJC6Tc+XuwGA3+9HPB5P3cYqsVgsrdu7Bcm5pZuwuTluSXRuwubmuCXRuQmbm+OWROd9yObudrviPTpJEonEgJd+9vv9aGxsTF2cYNasWRg2bBg2bNiQus3u3buxfft2HHfccQCA2bNnY/v27b0uOb1hwwbk5eVh1qxZac3W0dGR7nJcgeTc0k3Y3By3JDo3YXNz3JLo3ITNzXFLovM+ZHN3ux176dodd9yBU089FYcccgg6Ozvx3HPPYdOmTXj22WfR2dmJu+++G4sWLUJJSQl27NiBH//4xyguLsaCBQsAAAUFBbj44ouxYsUKFBcXY9y4cbjtttswY8YMzJ07FwBw8skn44gjjsCVV14Jj8eDtrY2LF++HJdccklaV1wDgPHjx2c7gRIk55ZuwubmuCXRuQmbm+OWROcmbG6OWxKd9yGbu9vt2DM6Pp8PV1xxBY499licffbZ+OCDD/Dcc8/hlFNOQW5uLrZt24YLL7wQxxxzDJYuXYrDDjsMf/rTnzBmzJiU46677sJLine truncated
"text/plain": [
"<Figure size 864x720 with 1 Axes>"
]
@@ -321,7 +332,7 @@
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAykAAAKSCAYAAADWEJRoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsvX14XFW5v/+ZZJK0SdMkTdMkfRsChkIV8CDI10IiAlLAtA1gBEVEj0SFU85BQOXwjaCHCIgofgUFCQfx1IIYhZQi4OsPaRGOUCjlTQw2TNomM03SeclkOu/z+yNMaJLJzszOfp41a3ju6/K6zpnM3Fnr3qvTLmb23jav15uEIAiCIAiCIAhCjlCgegCCIAiCIAiCIAiHIpsUQRAEQRAEQRByCtmkCIIgCIIgCIKQU8gmRRAEQRAEQRCEnEI2KYIgCIIgCIIg5BSySREEQRAEQRAEIaeQTYogCIKgNccccwyOOeaYjJ/vdDpRWVmJyy67jHBUgiAIwlyQTYogCMJ7iC1btuBrX/sazj77bKxYsQKVlZX413/91xmfPzAwgHvuuQef+tSncNxxx6G2thYOhwPnnHMOHnjgASST02+1ddlll6GystLwf+vXr6ecpiAIgqA5dtUDEARBEPj43ve+h1dffRULFizA0qVLMTo6avj8e+65Bz/84Q+xfPlynHLKKVi6dCkGBwfx2GOP4fLLL8fvf/973H///ZNe84lPfAIrV65M63vooYfw9ttv4+Mf/7hVU8Kjjz5qmUsQBEHIDWxyx3lBEIT3Dk8//TSWLVuGww8/HNu3b8e6detw3nnn4b777kv7/EcffRSVlZVoamqCzWabeHxwcBBnnHEG9u3bh//5n//J6JMRr9eLo48+GvF4HG+88Qaqq6stm1c2OJ1OHHfccfj0pz+Nu+66S8kYBEEQBGPk616CIAjvIZqbm3HEEUdM2nAYsX79ejQ3N097fn19/cTXxLZt25aR66GHHsLBgwexbt26jDco27Ztmzh/pLe3F5dccgmOOOIIVFVVYdeuXQBmPidldHQU1113HVavXo3a2lqceOKJuPPOO9N+RS1FKBTCHXfcgY9+9KNYtmwZli5dilNPPRX33Xdf2tc98cQT2LBhA4466igsWbIEq1atwtq1a/H9738/o/kJgiAI6ZGvewmCIAimKCoqAgDY7Zn9VfLzn/8cAHDJJZdk/bv6+vpw+umnY9WqVbjgggswOjqK0tLSGZ8fDoexYcMGvPjii1i9ejXa2trg9/tx22234Zlnnkn7mtHRUbS2tmLHjh049thj8ZnPfAYA8Kc//QlXXXUVnn/++UmfvNx///248sorsWTJEqxduxY1NTUYGRnBm2++iZ/97Ge4+uqrs56nIAiCMI5sUgRBEISsicVi+OUvfwkAOOOMM2Z9/t/+9je8/vrreN/73ofm5uasf99zzz2Hq6++Gt/85jczev6dd96JF198Eeeccw5+8YtfoKBg/IsDX/3qV3Hqqaemfc11112HHTt24Fvf+hauvPLKicfD4TAuvvhiPPjgg1i/fj3OPvtsAMDPfvYzFBcXY/v27ViyZMkk18jISNZzFARBEN5Fvu4lCIIgZE1HRwdef/11rF27Fqeffvqsz0+dXG/mUxQAWLJkCb7+9a9n/PzNmzfDZrPh29/+9sQGBQBWrlyJL3/5y9Oe7/F48OCDD+LYY4+dtEEBgJKSElx//fUAxr+ydih2u33iE6VDUXW+jSAIQr4gn6QIgiAIWfGDH/wAd999N44++mj89Kc/nfX5Pp8PPT09KC4unvgKVbZ84AMfQElJSUbPHR0dxe7du1FXV4fGxsZpPz/55JOnPbZjxw7EYjEUFBTg5ptvnvbzWCwGAPjHP/4x8VhbWxs6Ojpw0kkn4dxzz8WaNWtw0kknoa6uLtNpCYIgCDMgmxRBEAQhY2677TZ0dnbi/e9/P7Zs2YLKyspZX/OrX/0KwWAQ559/vulPGKZ+ncoIv98PAKipqcnYdeDAAQDAzp07sXPnzhndgUBg4v/euHEjampqcN999+Hee++d2LCdeOKJuP7669HU1JTxmAVBEITJyNe9BEEQhIzo7OxEZ2cnPvjBD+Kxxx7D4sWLM3pd6oT5z3/+86Z/d6ZXIwOAhQsXAgCGhobS/nz//v0zvuZLX/oSvF7vjP9LXVEsxQUXXIDf/e536Ovrw8MPP4z29na8/PLLaGtrw1tvvZXxmAVBEITJyCZFEARBmJXrrrsOt912G0488URs2bIFVVVVGb3uhRdewKuvvor3ve99bJ8slJeX4/DDD4fb7U67UUh3da8TTjgBBQUFePbZZ039zoULF+K0007D9773PWzcuBGhUAh/+MMfTLkEQRAE2aQIgiAIBiSTSVx11VX4yU9+gjVr1uCRRx5BRUVFxq+f6wnzZrnooouQTCZx/fXXI5FITDze39+f9jyaxYsX44ILLsArr7yCm2++eeIclEPZt2/fpHNS/vKXv6S9d4rb7QYAw0skC4IgCMbIOSmCIAjvIR577DH89re/BfDu155eeOEFXHbZZRPPOfReIN/97ndx3333obi4GCeccALuuOOOac6VK1fioosumva43+/HI488gpKSEtMnzJtl48aN+O1vf4vHH38cTU1NOOOMMybG85GPfARPPPHEtNfceuut2L17N7773e/ioYcewpo1a1BbWzvxiczzzz+P73znOzjyyCMBABdffDHKyspwwgknYOXKlbDZbNixYweeffZZNDQ0oLW1lXXOgiAI+YRsUgRBEN5DvPLKK3jwwQcnPdbf34/+/v6J///QTYrT6QQARCIR/OhHP0rrPPnkk9NuUrq7uzE2NjanE+bNUlJSgp6eHtxyyy145JFHcPfdd2PlypW4+uqrsW7durSblPLycjz22GPYtGkTuru78dhjjyEUCqGmpgYOhwM33HADzj333Innf+tb38Kf//xnvPLKK/jTn/4Eu92O5cuX4xvf+Aa+/OUvZ/WJkyAIgjAZm9frnf5ZtSAIgiAIgiAIgiLknBRBEARBEARBEHIK2aQIgiAIgiAIgpBTyCZFEARBEARBEIScQjYpgiAIgiAIgiDkFLJJEQRBEARBEAQhp5BNiiAIgiAIgiAIOYVsUgRBEARBEARByClkk2KA2+3W0g0Avb29ZG6du1D6KZsD+naXta7GL2ud303tl+78bkDeY1S4Za2r8efaWpdNigGUd0jmvvuylejcRbrnl5sandejdOd3c/gp0bW7NFfjl+78bg4/FWbGLZsUAxKJhJZuanTuIt3zy02NzutRuvO7OfyU6NpdmqvxS3d+N4efCjPjzstNitvtRiwWQyQSwdDQEADA5/MhEAgAAFwuF+LxOMLhMIaHhwEAXq8XY2NjAIDBwUEkEgmMjIxgZGQEAODxeBAMBgEAAwMDAIBgMAiPxwMAGBkZQSgUQiKRwODgIABgbGwMXq8XADA8PIxwOIx4PA6XyzUxHp/PBwAYGhpCJBJBLBab+EjM7/fD7/dbOqdQKDTnOe3fvz/tnADMeU779+8nnZPb7c7qOFkxJ6uOk8/ns2TtpZvTgQMHyOY0MDBg2dpLN6cUFMfpwIEDpH+ehoeHydaez+cje49IYcXaSzen1DGh+PPkdrstfS+fOqfBwUFt3sunzsnj8ZC973k8HrI5pX5O9V6eguI4jYyMWPpePnVOQ0NDZGvP5/NZ/m+j1JxSUP2dm/q/Kf48uVwu0n9HDAwMkL1HpF6f6XHKZk4ej2fanGbD5vV6k7M+S9CO3t5eNDY2qh7Gewpprgbpzo80V4N0V4N050eaqyHXuuflJylWkckuLxfd1OjcRbrnl5sandejdOd3c/gp0bW7NFfjl+78bg4/FWbGLZsUAwoLC7V0U6NzF+meX25qdF6P0p3fzeGnRNfu0lyNX7rzuzn8VJgZt2xSDFiwYIGWbmp07iLd88tNjc7rUbrzuzn8lOjaXZqr8Ut3fjeHnwoz45ZNigGpk4N0c1Ojcxfpnl9uanRej9Kd383hp0TX7tJcjV+687s5/FSYGbdsUgyoqanR0k2Nzl2ke365qdF5PUp3fjeHnxJdu0tzNX7pzu/m8FNhZtyySTEgFotp6aZG5y7SPb/c1Oi8HqU7v5vDT4mu3aW5Gr9053dz+KkwM27ZpBgwOjqqpZsanbtI9/xyU6PzepTu/G4OPyW6dpfmavzSnd/N4afCzLhlk2LA4sWLtXRTo3MX6Z5fbmp0Xo/Snd/N4adE1+7SXI1fuvO7OfxUmBm3bFIMSN1pUzc3NTp3ke755aZG5/Uo3fndHH5KdO0uzdX4pTu/m8NPhZlxyybFgKKiIi3d1OjcRbrnl5sandejdOd3c/gp0bW7NFfjl+78bg4/FWbGLZsUA8rKyrR0U6NzF+meX25qdF6P0p3fzeGnRNfu0lyNX7rzuzn8VJgZt2xSDBgcHNTSTY3OXaR7frmp0Xk9Snd+N4efEl27S3M1funO7+bwU2Fm3LJJMaC2tlZLNzU6d5Hu+eWmRuf1KN353Rx+SnTtLs3V+KU7v5vDT4WZccsmxYBIJKKlmxqdu0j3/HJTo/N6lO78bg4/Jbp2l+Zq/NKd383hp8LMuGWTYsDY2JiWbmp07iLd88tNjc7rUbrzuzn8lOjaXZqr8Ut3fjeHnwoz487LTYrb7UYsFkMkEsHQ0BAAwOfzIRAIAABcLhfi8TjC4TCGh4cBjF8aLRVwcHAQiUQCZWVlGBkZAQB4PB4Eg0EAwMDAAAAgGAzC4/EAAEZGRhAKhZBIJCa+dzc2NjZxybXh4WGEw2HE43G4XC5UV1cjEAjA5/MBAIaGhhCJRBCLxeB2uwEAfr8ffr/f0jmFQqE5z6mwsDDtnADMeU52u510TgCyOk5WzMmq41RdXW3J2ks3p/LycrI5RSIRy9ZeujmloDhOCxYsIP3zNG/ePLK1V11dTfYekcKLine truncated
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAykAAAKSCAYAAADWEJRoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsvXt4XGW5v/+ZZJK0SdMkTdMkPQ0BQ6EKuBHkayERATmZtgGMoIjolqiwy96IqGy+EXQTARHFr6AgYSPuWhCjkFIUPP6QFmEL5VBOYpAwaZvMNEnnkJnpnOf3R5jQJJOVmZX1PO+8w3Nfl9e192Tmzvve6+20LzNrLZvX601BEARBEARBEAQhTyhSPQBBEARBEARBEISDkU2KIAiCIAiCIAh5hWxSBEEQBEEQBEHIK2STIgiCIAiCIAhCXiGbFEEQBEEQBEEQ8grZpAiCIAiCIAiCkFfIJkUQBEHQmqOOOgpHHXVU1s93Op2orq7GpZdeSjgqQRAEYT7IJkUQBOFdxNatW/HVr34VZ511FlatWoXq6mr867/+66zPHxoawl133YVPfOITOOaYY1BfXw+Hw4Gzzz4b9913H1KpmbfauvTSS1FdXW34vw0bNlBOUxAEQdAcu+oBCIIgCHx897vfxcsvv4xFixZh+fLlGB8fN3z+XXfdhR/84AdYuXIlTjrpJCxfvhzDw8N45JFHcNlll+H3v/897r333imv+djHPobVq1dn9D3wwAN466238NGPftSqKeHhhx+2zCUIgiDkBza547wgCMK7hyeeeAIrVqzAoYceih07dmD9+vU499xzcc8992R8/sMPP4zq6mq0tLTAZrNNPj48PIzTTjsNe/fuxf/8z/9k9cmI1+vFkUceiUQigddeew21tbWWzSsXnE4njjnmGHzyk5/EHXfcoWQMgiAIgjHydS9BEIR3Ea2trTjssMOmbDiM2LBhA1pbW2c8v7GxcfJrYtu3b8/K9cADD+DAgQNYv3591huU7du3T54/0t/fj4svvhiHHXYYampqsGvXLgCzn5MyPj6Oa665BmvXrkV9fT2OP/543H777Rm/opYmHA7jtttuw4c//GGsWLECy5cvx8knn4x77rkn4+seffRRbNy4EUcccQSWLVuGNWvW4IwzzsD3vve9rOYnCIIgZEa+7iUIgiCYoqSkBABgt2f3V8nPfvYzAMDFF1+c8+8aGBjAqaeeijVr1uD888/H+Pg4ysvLZ31+JBLBxo0b8dxzz2Ht2rXo6OiA3+/HLbfcgieffDLja8bHx9He3o6dO3fi6KOPxqc+9SkAwJ/+9CdceeWVeOaZZ6Z88nLvvffiiiuuwLJly3DGGWegrq4OY2NjeP311/HTn/4UX/nKV3KepyAIgjCBbFIEQRCEnInH4/jFL34BADjttNPmfP7f/vY3vPrqq3jPe96D1tbWnH/f008/ja985Sv4xje+kdXzb7/9djz33HM4++yz8fOf/xxFRRNfHPjyl7+Mk08+OeNrrrnmGuzcuRPf/OY3ccUVV0w+HolEcNFFF+H+++/Hhg0bcNZZZwEAfvrTn6K0tBQ7duzAsmXLprjGxsZynqMgCILwDvJ1L0EQBCFnurq68Oqrr+KMM87AqaeeOufz0yfXm/kUBQCWLVuGr33ta1k/f8uWLbDZbPjWt741uUEBgNWrV+OLX/zijOd7PB7cf//9OProo6dsUACgrKwM1157LYCJr6wdjN1un/xE6WBUnW8jCIJQKMgnKYIgCEJOfP/738edd96JI488Ej/5yU/mfL7P50NfXx9KS0snv0KVK+973/tQVlaW1XPHx8fx5ptvoqGhAc3NzTN+fuKJJ854bOfOnYjH4ygqKsKNN9444+fxeBwA8I9//GPysY6ODnR1deGEE07AOeecg3Xr1uGEE05AQ0NDttMSBEEQZkE2KYIgCELW3HLLLeju7sZ73/tebN26FdXV1XO+5pe//CVCoRDOO+88058wTP86lRF+vx8AUFdXl7Vr//79AIAXXngBL7zwwqzuQCAw+X9v2rQJdXV1uOeee3D33XdPbtiOP/54XHvttWhpacl6zIIgCMJU5OtegiAIQlZ0d3eju7sb73//+/HII49g6dKlWb0ufcL8Zz/7WdO/O9urkQHA4sWLAQAjIyMZf75v375ZX/OFL3wBXq931v+lryiW5vzzz8fvfvc7DAwM4MEHH0RnZydefPFFdHR04I033sh6zIIgCMJUZJMiCIIgzMk111yDW265Bccffzy2bt2KmpqarF737LPP4uWXX8Z73vMetk8WKisrceihh8LtdmfcKGS6utdxxx2HoqIiPPXUU6Z+5+LFi3HKKafgu9/9LjZt2oRwOIw//OEPplyCIAiCbFIEQRAEA1KpFK688kr8+Mc/xrp16/DQQw+hqqoq69fP94R5s1x44YVIpVK49tprkUwmJx8fHBzMeB7N0qVLcf755+Oll17CjTfeOHkOysHs3bt3yjkpf/nLXzLeO8XtdgOA4SWSBUEQBGPknBRBEIR3EY888gh+85vfAHjna0/PPvssLr300snnHHwvkO985zu45557UFpaiuOOOw633XbbDOfq1atx4YUXznjc7/fjoYceQllZmekT5s2yadMm/OY3v8Fvf/tbtLS04LTTTpscz4c+9CE8+uijM15z8803480338R3vvMdPPDAA1i3bh3q6+snP5F55pln8O1vfxuHH344AOCiiy5CRUUFjjvuOKxevRo2mw07d+7EU089haamJrS3t7POWRAEoZCQTYogCMK7iJdeegn333//lMcGBwcxODg4+f8fvElxOp0AgGg0ih/+8IcZnSeeeGLGTUpvby+CweC8Tpg3S1lZGfr6+nDTTTfhoYcewp133onVq1fjK1/5CtavX59xk1JZWYlHHnkEmzdvRm9vLx555BGEw2HU1dXB4XDguuuuwznnnDP5/G9+85v485//jJdeegl/+tOfYLfbsXLlSnz961/HF7/4xZw+cRIEQRCmYvN6vTM/qxYEQRAEQRAEQVCEnJMiCIIgCIIgCEJeIZsUQRAEQRAEQRDyCtmkCIIgCIIgCIKQV8gmRRAEQRAEQRCEvEI2KYIgCIIgCIIg5BWySREEQRAEQRAEIa+QTYogCIIgCIIgCHmFbFIMcLvdWroBoL+/n8ytcxdKP2VzQN/ustbV+GWt87up/dKd3w3Ie4wKt6x1Nf58W+uySTGA8g7J3HdfthKdu0j3wnJTo/N6lO78bg4/Jbp2l+Zq/NKd383hp8LMuGWTYkAymdTSTY3OXaR7Ybmp0Xk9Snd+N4efEl27S3M1funO7+bwU2Fm3AW5SXG73YjH44hGoxgZGQEA+Hw+BAIBAIDL5UIikUAkEsHo6CgAwOv1IhgMAgCGh4eRTCYxNjaGsbExAIDH40EoFAIADA0NAQBCoRA8Hg8AYGxsDOFwGMlkEsPDwwCAYDAIr9cLABgdHUUkEkEikYDL5Zocj8/nAwCMjIwgGo0iHo9PfiTm9/vh9/stnVM4HJ73nPbt25dxTgDmPad9+/aRzsntdud0nKyYk1XHyefzWbL2Ms1p//79ZHMaGhqybO1lmlMaiuO0f/9+0j9Po6OjZGvP5/ORvUeksWLtZZpT+phQ/Hlyu92WvpdPn9Pw8LA27+XT5+TxeMje9zweD9mc0j+nei9PQ3GcxsbGLH0vnz6nkZERsrXn8/ks/7dRek5pqP7OTf/fFH+eXC4X6b8jhoaGyN4j0q/P9jjlMiePxzNjTnNh83q9qTmfJWhHf38/mpubVQ/jXYU0V4N050eaq0G6q0G68yPN1ZBv3QvykxSryGaXl49uanTuIt0Ly02NzutRuvO7OfyU6NpdmqvxS3d+N4efCjPjlk2KAcXFxVq6qdG5i3QvLDc1Oq9H6c7v5vBTomt3aa7GL9353Rx+KsyMWzYpBixatEhLNzU6d5HuheWmRuf1KN353Rx+SnTtLs3V+KU7v5vDT4WZccsmxYD0yUG6uanRuYt0Lyw3NTqvR+nO7+bwU6Jrd2muxi/d+d0cfirMjFs2KQbU1dVp6aZG5y7SvbDc1Oi8HqU7v5vDT4mu3aW5Gr9053dz+KkwM27ZpBgQj8e1dFOjcxfpXlhuanRej9Kd383hp0TX7tJcjV+687s5/FSYGbdsUgwYHx/X0k2Nzl2ke2G5qdF5PUp3fjeHnxJdu0tzNX7pzu/m8FNhZtyySTFg6dKlWrqp0bmLdC8sNzU6r0fpzu/m8FOia3dprsYv3fndHH4qzIxbNikGpO+0qZubGp27SPfCclOj83qU7vxuDj8lunaX5mr80p3fzeGnwsy4ZZNiQElJiZZuanTuIt0Ly02NzutRuvO7OfyU6NpdmqvxS3d+N4efCjPjlk2KARUVFVq6qdG5i3QvLDc1Oq9H6c7v5vBTomt3aa7GL9353Rx+KsyMWzYpBgwPD2vppkbnLtK9sNzU6LwepTu/m8NPia7dpbkav3Tnd3P4qTAzbtmkGFBfX6+lmxqdu0j3wnJTo/N6lO78bg4/Jbp2l+Zq/NKd383hp8LMuGWTYkA0GtXSTY3OXaR7Ybmp0Xk9Snd+N4efEl27S3M1funO7+bwU2Fm3LJJMSAYDGrppkbnLtK9sNzU6LwepTu/m8NPia7dpbkav3Tnd3P4qTAz7oLcpLjdbsTjcUSjUYyMjAAAfD4fAoEAAMDlciGRSCASiWB0dBTAxKXR0gGHh4eRTCZRUVGBsbExAIDH40EoFAIADA0NAQBCoRA8Hg8AYGxsDOFwGMlkcvJ7d8FgcPKSa6Ojo4hEIkgkEnC5XKitrUUgEIDP5wMAjIyMIBqNIh6Pw+12AwD8fj/8fr+lcwqHw/OeU3FxccY5AZj3nOx2O+mcAOR0nKyYk1XHqba21pK1l2lOlZWVZHOKRqOWrb1Mc0pDcZwWLVpE+udpwYIFZGuvtraW7D0Line truncated
"text/plain": [
"<Figure size 864x720 with 1 Axes>"
]
@@ -332,9 +343,8 @@
],
"source": [
"points = [2, 6, 7]\n",
"grid()\n",
"show(pct2, mph2, xlabel='Grade (percent)') \n",
"plt.plot(points, points, 'gD-', lw=1);"
"plt.plot(points, points, 'gD-', lw=2);"
]
},
{
@@ -390,7 +400,7 @@
"cell_type": "markdown",
"metadata": {},
"source": [
"My **Eddington Number is 54**: I've done at least 54 miles on at least 54 days. This is a bit above [the mean for Strava users](https://swinny.net/Cycling/-4687-Calculate-your-Eddington-Number). \n",
"My **Eddington Number is 54**: I've done at least 54 miles on at least 54 days. This is a bit above [the mean for Strava users](https://swinny.net/Cycling/-4687-Calculate-your-Eddington-Number). Eddington was about my age when he died with a number of **77**, and his available roads, bicycles, and navigation aids were not nearly as nice as mine, so good for him.\n",
"\n",
"If I wanted a bigger number, I could use kilometers:"
]
@@ -464,6 +474,13 @@
"Egaps(mls, range(54, 67))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"OK, if I do one 65 mile ride every 3 weeks, I could have my Eddington number match my age before I turn 65."
]
},
{
"cell_type": "code",
"execution_count": 15,