Skip to content
Snippets Groups Projects
Commit 7ce5646f authored by Hanna Schilbert's avatar Hanna Schilbert
Browse files

Upload New File

parent b8c68ccf
No related branches found
No related tags found
No related merge requests found
%% Cell type:markdown id: tags:
# Python course 2023 - Exercises C
%% Cell type:markdown id: tags:
## Part1 - file handling
%% Cell type:markdown id: tags:
---
1.1) Count number of sequences (number of headers) in AtCol0_Exons.fasta!
%% Cell type:code id: tags:
```
from google.colab import drive
drive.mount('/content/drive')
```
%% Output
Mounted at /content/drive
%% Cell type:code id: tags:
```
datei = open("/content/drive/MyDrive/Python_Kurs_2022/Python_course_2021_data/AtCol0_Exons.fasta", "r")
lines = datei.readlines()
datei.close()
```
%% Cell type:code id: tags:
```
```
%% Cell type:markdown id: tags:
---
1.2) Count number of sequence lines!
%% Cell type:code id: tags:
```
```
%% Cell type:markdown id: tags:
---
1.3) Count number of characters in document! (How many per line?)
%% Cell type:code id: tags:
```
```
%% Cell type:markdown id: tags:
---
1.4) How long are all contained sequences combined?
%% Cell type:code id: tags:
```
```
%% Cell type:markdown id: tags:
---
1.5) Calculate the average sequence length in this file!
%% Cell type:code id: tags:
```
```
%% Cell type:markdown id: tags:
---
#Bonus tasks
%% Cell type:markdown id: tags:
![image.png](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAACQ0AAAA4CAYAAABdEUznAAAgAElEQVR4Xu2dX8guVfX4p0ARo4vEA3UpBUcQE9QOHPHPTaCCSCiVkqkogiLHpChJQQ+mKRZCnUQDTRAkb0LDm0SyA4XQH43I0CSIMrAupC4UUoN+v/XEfr/7nTN7rz+z53lnnufzQpBnZvZe67PWXmvtmfXMfOC///+v4w8CEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQ2BoCH6BpaGtsjaIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhBYEaBpCEeAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIDAlhGgaWjLDI66EIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQoGkIH4AABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgMCWEaBpaMsMjroQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCgaQgfgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAwJYRoGloywyOuhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEKBpCB+AAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIDAlhGgaWjLDI66EIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQoGkIH4AABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgMCWEaBpaMsMjroQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCgaQgfgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAwJYRoGloywyOuhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACENi4pqFf/OIX3Xnnndd95jOf6Z5++mksDAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQj0CGxc09C///3vbv/+/d3f//737uc//3l34MCBRRj9P//5T/foo492r7zySvfAAw90J5544iLktgj5xhtvdN/61re6I0eOrE4//vjjux/96EfdJZdcYrl8tudsis2WumZaOca26l9alxdddNFgLNpWTq38bBPGwQc2wYrr02GT/WWTdVufhzDTNhKY09rZlDp+yI/mwnkuckTX2tLlj+i9jn17ies28o7YiGt2E8Bv2nkELNuxlJE2nedS9dvk+q+tB69vtKX60roILd1n52zfddS9Vj+ZMyerDpwHAQiUCUTX+JQ5YMqxNV+I8tDG5XgbApM0DUnSvfnmm7tnn322yRt/xIF/+tOfdnfffXf34osv7mh+wQUXdLfeeutqjvzvBz/4QXf99dd3d911V3f48OEqqT//+c/dWWed1b3zzjvVJqP0BiNpeKk1I0UdPsnxz3/+0yR3G/NPP4rY7rLLLlv5gvydc845KxtabDO9dONm2BSbRX12HL35XL2N+tfW5TXXXLOKif1YtI2c5uOl85AEH5iHHZYixSb7yybrthT/Qs5lEpjT2tmUOn7IE+bCeS5yRFfL0uX36r2ufTtNQ17LcH6NwLat0ym9AZZt6W46z6Xqt8n1X1sPXt9oS/WlPqGvfOUr3YMPPngMOHmO9PnPf7677bbbutNOO80Nduk+O1f7rqvutRp8rpys8nMeBCBQJxBd41PmgCnH1vwhykMbl+NtCDRvGnruuee6K6+8cvXAWf7GfiZMGpCuvvrq7ujRo6vxpNg6++yzu1dffXU1x5lnnrlq4snfzJMc/pRTTjnmWB9bctC//vWv1UaW1Igk19caXlJz0Yc+9KFO5Ni3b5/JUqLnueeeu3pDUv8tPGnusSxNgjQ+KcqjsRih4TTuNZuFJpzoIk2PdQdpTZ4ahjHXlsZdt/4Tmdk1bG1dlvx6Gzm5oG7ByfjAFhi5oYpL9hct1yxZt4YmZqhGBDR/azTNLIaZ09qZoo5fly21eebCeS5yRJ1/6fJ79V7Xvp2mIa9ltvv8ucU7TZ4lW2vbYt5YW2m+sHSem6rfFPXfWF/a9uuXvlaS/VLTUHpuJf/+7rvvdi+99NKOiSM/oJ7CZ7X13dIn52rfddW9VpZz5WSVn/MgsO0EtLgaXeNT5IBkqynHnorHtvvZuvRv1jSUv9JPCiRpmpGmnjGNLnlDz+mnn9798Ic/3NWVLV3Bf/jDH7ozzjhjF6/ULSwNTJZPlKXCbqgBSQbudx+XzpNz01hj9O4bX1tk63KWyDypCKsxi4y7jmuWzD3no+kRTVpRG2jy1MYdc21p3HXrH+XW8rrIutxGTi2Zb8JY+MAmWHF9OizZX7Rcs2Td1ucBzGQloPmbdZwlnLfpa2ddttTmmQvnucgRXRtLl9+rd2R/4J1DzqdpKEJte6+ZW7zT5FmypbYt5o21leYLS+e56fqNtT/XtyOw9LWSSKTnQf3GIHmmJG8ZkrcQaV+vaEe1PpK2vlvKMVf7rqvutbKcKyer/JwHgW0noMXVbVvj8Fj2imjWNJSKoy9+8YurYkgcQwqiMc0zybm8b+0Rk5SKtSFzad3F+au65PqSPHlzkTQ4XXHFFU28Q1tkTSaZaJAkO01DEwE2DKv5z7qTliZPTaUx15bGXbf+BpNNfkpkXW4jp8kNsbAJ8IGFGWyPxV2yv2i5Zsm67bFbMP0AAc3fNgnapq+dddlSm2cunOciR3QNLV1+r96R/YF3DjmfpqEIte29Zm7xTpNnyZbatpg31laaLyyd56brN9b+XN+OwNLXSiJRew6VPy8a85ysFXVtfbeap1b3tZwjMta66l6rbJuyDqz6ch4ENo2AFle3bY3DY9ke3qxp6JFHHuk++tGPrpqE5G/sG3fSQpLPhkUacDwdw9pcaaxLLrmk+8AHPtA9++yzgzKl5qL33ntv16fJ0vhvvfXW6t9ffvnlnU+4pWaaD37wg93+/ftXnydLb0dKi2vIxYa606UIffTRR7vbb7995/Nwp5566uq/pZnL+vfmm292P/7xj1dj5a/RvOCCC7qHHnrI9A3emuw1nXMZtRuKiaecl3/CTt5KJbIfOHBgUGXR78iRI534bPqMnnC69tpru0OHDnVPPfVUd/311w9em3PXgn2yR86x9i3jvp949eoLbPWfvh4f+9jHuhtuuKGTN3XJX42nx1es8gyB91zrkUnbwMjav/TSS7t33nnnmLeWedebnH/nnXce43eXX355d88996jL06uXl6O2Lq3+Pjb+jFkHXkb9ud5+++3uuuuu2/kcZj/m5Z/flLV81VVXrWJJ/nnMnLvXR0pOMFYvbyzJ3xwoMomuEsc++9nPHpOnao7rlbs2lnf9RNiX9JZPrl522WWrmJi/vVBbE9pxj4xj1oVw1fJe7sMeuWo2m5P+Yrs77rhjV02Ty57iXzT/e/Jmq3Uu43jXRX9u7xq11LPJl6zM5XzNV6z14F7G8Mga9dQ2e23rvcojQ1zzGrW0z7H46tDeK68LPfuMddnSOs+Y2j7FlhZ7ynXL4Y1pKSZG6x5PvrT4pMTDsWs96eTZh1r27aVaN9UYnnsHWkzP74doGyQvL4/NNP+gPuy62r0mr202+V7GmHVSu4eY1mXEr4fWllaPeeJKNJ/mcnlqyX7ss9yTsMbluftmC9a5jkurNUR2j68M+bmleaD08Muz/iL7Bct+21Oz5v5Syrdank7ztdp7afWi9lxjjA2Gns/Uap6aPbQfrw89J7PEob3as7S6b1HLLd5awZuHhnTw1L3e+Sz2LPmXtg60Z10eWbWXJ2gx0frMd8p77t59XyT+ahzEv7TGiNwHo/5uqWnYtxza9WzGkxfG3ludum5qGUP7MXGKsafmEan5tXsKaa9W6xdIsdsbeyL+5fVfi37ec5o1DfUntiaQksClBhyrgp7r847voe/LJl2keenXv/518Q1KpUalPDHJQysJ9ukvvbXowx/+8DEPY2Xjc/fdd3evvvrqTnPLOeecs7r0uOOO29XAI3NcfPHFq4fd8oBXGmDef//97plnnllda33TT/+tSlKcn3TSSd2LL764mtf6Ks2a7DKmLMJS4ZvYWDYqOU9hk+SUMYaazfr6yTX5N4bF/gcPHjRxrwXW3B4iS5+j/NsTTzyxq5mr5CcWvYbWhdV/cj3uvffe7pvf/ObKZ7R5vb5ilWeMLl6ZZK6SHfOx+nHBu97SHNIEKb5w/vnnr9SUprz//ve/q0aE2sY0opeXo7Yurf4+Jv7k9pAbEZ71HWFkWXOiz89+9rPu/vvvXzWMyn/LZi01VJZ+JeT1ES0X5g2OWky26OWJkSmuXnjhhd2vfvWrwSa6vvwRe5QYeNdPhP1Qbkh6S9OtPFCU/27VNOSVsbVN+3nv8OHDK/xeuUo2q8XWyDxj9E81nKxdadI84YQTdmojkUVy3cc//vFVQ6fY11J3RfPmFOt8nXnFUs/u27dv54cDFuatmoZK+WJdMTzio9a6yBsDW8TjiD5p3lo89eSRfk3w8MMPd1/60pdWNarULP/4xz929kdDdZr8GKNfS8iYtb1XtA5Zly2988jDIE9tH4nP1jwwtRzRuiPqr958aYmf6X5ANK4nW3j3oZZ9e2mfEuGu7fGtTUPe2Oi1mfCkPozVh17bePet+b2waE3m9V1r/B3SxTtXLRfleUxqrohfa3uuoTXojSvRfJpk89TvKT55WVjisjDu/1l9YV2+OZb1uvVrXWt4fWUoB6U1OvRjQZG39Myglc9p9161Neu9d1a7t1fbx1v2B5G9Vy7P17/+9e673/2uqd6P+FJ03dfq3X7cGnquJOdoTUOe50XRdW+NXyV9x+Szfm7x1gqRPOStL9J98XSPIj1nk3HGPNfp5+4huaI5o78OLLLm7IeenSVfHfrSiucrK1PFjDF+6ImX0dxQs691n+fNLzIn+5bDK/Redt5YtO66Sds7W31KiztDzzsiY1vzTKuYN/Y5ZGntDD03GRN7rCy9/mupTyLnzLZpyNLNWVNYK4L715bmS8kwvWFA3iJ03nnnDTbhlLrK8+Aj80oDjzSMnHbaad2f/vSn7hOf+ET1V9Za12qesOUB9pNPPrnTgKA1RPU5iPPL2zNuvPHGXQ0t/Tmefvppk7/V7KjZyBIUE8/UdCFypm8F9xulNE7yFphPfvKTO5+V07iX5MvnkeYCKcDE1vKXuorvu+++YxqwhvzEopdmCKseErw8PKO+oslT00e7NiJT7aaBPJga+h61vHVEGkis602TW7NhRC8Lx6FmQm3d9Td62rpKrEob576c0XUQYdSfK7enNImce+65XVoXImeug+SEiy66aLCRsiWTFnrldrbGyP418rD2lltuWZnL0jwakbvks571E2GfX+PRO5rDIjJG14U2V573tHOja7lFzIjqX/oFlXajQ/O5qDwlH4+sF03GKfKKpZ6NMI+upTnF8DE+odlSO74uW1tq0mg8LelQ8zm5JrHp5yWLr2r1jrcuzuUpNRSPtWXipI0T9cep8oCXZUSOSByN+mtEPotPanbV1nraa6Z87dmH5v5r/bGRXBPhrq09a9OQh1fEZlH/iOa0iIxTrfWx9aHHNppfa2NFGUR81xLnh/SJzGWJGRGfqfGey/2tSC0ZYWFhbLmnUvsRkTQwe3NgK3+x1G17oV/EViU5I74y5OfaffOhB8cRPaLxylMnW+wezVNT7b1qXGr1fisb9J/PaHmpdLz2pqFcx/y+oiUOaXWTN8ZEc1nruk/L7znn6HMWS4wbqnuj81ns2Xo9R2W1PsPsNxWleNj/ysqQXlPFjHXlyUhu4L52vdE2UhdpsX4u+xYtpkVrgFIO0Oaz5LIpx9bki/LQ/MH77EIbL/evSOzROJRynfWZs8XOkXNm3zQkgKTB4tZbb1294SB/20Lts1u541k+b1bqHk3/fsopp6x+ga69FrL0q5wUFIe6dMVwtaJdcy7ttYKeT7XVnCgyztRNQ/0FlIrYs84665i3YWic+rpr3Es2S/OUHqz3g1FqwMp91qOXtvCteshGo9W8NV/R5LEU9ZFvQGtvAktrVxrH0nqVT1U99thju0TS/GhoHu01tZoN57AuNX8vxTZv3JhiHWi2T77fb4bMX6lYax7rH4v4SMQHrHrljaRjYmTp4axXdq9PeNZPhL12jfZQ2vtab22+IT7RdaHNldtOO9drt5YxI6p/st1Q06LlmOUhQKu8qd0M79/M8qyLVms0t0Mp5lu4lt4K411Lfb/Yyxge9VHLTdu9trXHx7U44s0jOVfthm7pBrx377UEW2q1dFQHzX7RPODdY7SWoyS3Nk/JX7XrtDxe8skWaz26D81jkadpKLJH0R5+WZuGPLwiNtOuoT78n/Wn3vtOFe8ivmvJ2eusuTQfjcbs/hqMxpVoLorUkhEWlrp2zL2pqP4R35xirqnWXsRWJSYRX9EekA3lwCEWET1a2yk6Xu35Q+35xFR7L02P0oP7sTYo1WPeOJ7OL9Ul8oPEm2++efXD0/6cljik1U3eOnuKXFaqB2r+FKnjvM9ZLDF8aM23yHte/9LWQem5ZVTW0vO6NF56A1rpfrflucxUMWOdebL2XLN2rC9jxN89z1m0eMi+Zfp9y1R1UykHeHyqtGamHHsqHpqve/c/2njWmqA0r8dOmixe3ayyD50326ahBDQXeuhVgKUkpf2KvA+jdH4yRpqnVKz2m4vyV3hrBXmtiLIUc9oi9HQB15wpMs6UTUPeh0veT+ZpXLXAWiughrhofqId1zbS2sNPL8+or2hcLUW9pTjtj1Py35xr/hmqoYYhy3ocmifp3P/F75jgna5d17rUbq6UbOKVT/Nz7fgQU4vt89cwpjG0ZFyKKZqPe5mU/CSqV9SWEfYee2gxzLJ+Iuy1a7SbNN7Yqc03ZFeN/Zh8lLhH5KrFsKiftdQ/ciPZEuej9ojE/NI634u8oumds/M0amnjRtfgOmN4VAeLv83R1tH1rXHqrxHL+ZGaurb30uasHdfiaCtbavNEddDG9dYvc5FDi6PePVKEk8YijwWWeqeUTyz73dL+3HPj25LPvHWqhVE+r2c9RWymXRPNTdH42bI+svjJmPrQYxvNl6J2SON6/Uquq8U6TR5Nn8g9iqF9qqV+aBWzLf4SycXaWvDUkppdIutHs6U2p+Z72nHPPlobSzs+NNdU+mnjevw2stcrsSg9kC89K4joodlBO+6tk1vnKU2+6N5LGzePd3mDxxQ20NZ97Xj+qbyzzz57dWr+eZOhhheL7q3taMkfEQ7eus9TK0TzUE2PWt0bnc9iz5JM2rWl41FZSzFPuNx0003d888/v/oCycknn7x6iUJ6zul5AK7pFI0ZNbt6/VCrF725oSSbx98jsU27pnUc0eaL1F1jfHwO+xaNibYevDby+JQ37rQYeyoe2rieOlK4WGKoJSdq95ws93ha62aRu3TO7JuGpKh69NFHd30qKy9w5P8PvUlIe4XcEJAhJ0n/ls8xZMCaUbWgILKMuQGdZBTnO3jw4DGqvf/++90zzzxzzJt3NMd57bXXuqNHj3a//e1vV6fKOI8//vgxnfHRIkzj4g2WtUQf8QdtoQ7JZ51nqHMwykOzY0SPfExNLjnX4yuaPBZ/sjQNWWVK+uWfoar9wjay3vpzSOEtn9FLn63TbJgft+pl4dji82QRHjXZNH/Tjnv8URtL28CUfLk1k8TLantNL60QL31KThu3ZFer3FoRm9Zobf1E2Gub3pa5KC9IPTlbYz8mHyXuEXaRtRyZJ6K/yFZaw3lMHqohtTwVlccS663rZS/yiqZ3lLk2bnQNrjOGR3UQZlZ/s8RAzce8/uVtiozGUy3+1948UnurS+26qF+N2bO1WrdWn4naz5OfInmgtmeL5slIrRz11ynymMg/1j/G7EPzWBR909DY+KLF0b6veXhFbBb1D00PrQ73rL/IXFY/GVMfemyj5a2p4l0kZlhytqbP2HWSjx/xa2/MtvpLy/tbkfo9wkJbP5ot5+Sbmi7a8SFdp9IvYquSLSK+UmKh/UC4/ymeiB6aHbTjpTzorfe0eaI1cnTvpclT2mNOYQNt3deOD/34Xc6XH8Bfe+213aFDh3YaL7RaOJ8nao8xexYLh7H5zForjMlDNT1KTUNj5rP4cnTfO+ae31COHmqITP/2wgsvdK+//np34403ds8999yqaejAgQPVZ5dDemk8ojEjUsNpsrTKDZp9Lfd0IrGNfcv/yEfYWWORFhenqpu0HGDxKc0v+3VECyZT8YjYuKS/Nd4PXe/NgRY7tdRN81ft+GybhjTH0oyqHR8C00+iH/nIRzr5Dl6eIOW6oc6xoeaiFkWgjFFjkeupGdv6akTR98orr9z5FFx/XOs4ueye5gSNWyTRR/xB88ExBVvLjlvN7hE9LBsUOSfiK5o8lqK+1jTklSm347333tt99atfXYmgNSNq3PvrRHzwy1/+cnfkyJGdSy1dpulkr14Wjp51qfm7l4e3YNHiQsQftVgS2cDMISZremkPK1o1DbX0Wcv6ibJf5+YqKmPEpp68F5Wrtu5bxoyI/iJbrpc0oV9++eXdCSecsGqkls/dlh6MankqKk+NV2S9WNZFyzk1vaPMtXG1TXLpRvk6Y3hUB63OT/abm63XlUc0rsIn8qBSrov61dgb8GNtafEZjVvL+OzNA9oeY0w+8sbRSP6PyqfZpMVat+Z97Vdw3qYhL/fo2hvyNct6itos4h+12FLby0Rl1PxKW+ulmj/3R7kfJp9U0f5a733TfFPWZF7ftcTfEifvXJptoz7jjdlj4oqmg+UBnaV+j7LQ5NN8fk6+qemiHR/SdQr9orYq2SKy17PUcdrbbKJ6aHbQjvc5aOdH8230usjey5I3h+r99LzGm580Ztq6rx3X6oahay3yRO1h8XXLD3L7crfMZ946rla3eN8uMbZpKPJcx5uDPfumCJt+nE8+IzFQPv2ejqfnI0nnU045Zdfbh0p6af4djRkyX0s/1OLQkJxajmTf8j8C+f5A84el71s0n4jor/mmJYZG487YsafgEa2/LHWkti9OY3hjj1xnYdlatzH1jFw7+6Yhyxs/hoqcHPTQw39LsbZv377urLPO6vrJMC1y6RCTsWX+/fv3d2+99daqoUiu0xK7p9jXFlmkMC05TkqEcvyOO+7obr/99p1OeG8BJmOUijAt6NWOR4Ks9WZHzkXjPiaZRR5waHqXbBrRw+K/UV/R5KkFNe3aiEx9rq+88kp3/fXXd6VvJY9db+KL8vY0WVvy0NrShBfRy8JxbNOQzDGWhycW1uJChJG2pqIbmJZMptCrpLcmt8Yrt2VEbktBo60fTYehObRrWt+k0eaz1CfyKx8tTnvzXkSums2iftZKfxkn+aG8neqNN95YvT1R/mq/9strF+8na5LsnrWSyxmtvbR1McQ0skYtekWYa+NG1+A6Y3hUB4u/5fabi62j61vj5K0J8vWT1zSWeaJ+1eoGfMSWiY9WD2v6R+1nydNabrIcj+SjSEzT5mnJSbNJn23EP6x5f2gfmsciT9NQhHt07dX8T+Ol2Zr6sEw3wm5s3sqvnyreRXzXm7OTHpG5LDFjrG0sOXdMXNF0KB2P1JIRFpp8Ws6bk29qumjHh3SdSr+IrUq2iPhKjUX/MzSSD1PjZP+5QkQPzQ7accuatdRZ2jzRPB3de2nylOr9KWygrfva8ankidqj1Z4l13mqfFar48bkoZq9xjYNRZ7r1OTR1sHQ8bFs+s/5/vjHP3bnnXdelx6eJx3TfbHEzPpwXdMpGjOm8MOarP1jtdxgiSHsW/6+8/YqLWdZfTyNE4nDrfYtU9VN2joS+TWfKvnllGNPxWOsjfv2TnWeJa5FYo/Hv1rqZolFtXNm2zRUupmVK1MDmZw+NfZcccUVJlb5mJ/+9Kd3JcvS3Om86MMlGXdMMZdkjnSI96HUmC65aUj09HKKBjfLPEPFqRaoteMlB4/qkcaL3DCv+Yomj6WoL/l6xH9rxd/QTXOLfS3BxhOjInpZOLZsGmoRf7RYWDseYaStqegGppWP5HFrqHhp/a1mTW6NlzU/R3JJ359L60fTYWhdaNe0vkmjzTcko8Zei9OW9RmRK3LTIzJPRP+0uey/LdISn7U8FZGnNm8kftX8xFL7RubU9I4y18aNrsF1xvCoDmJHzd/maOtozNE4lWJ97TNjQzHFMk/Ur8bs2cbaMl2v+Yymf9R+lviZn7NOOSIxTctHLTlpLEpsPfsF63639JCk9mOfknwR7tG1Z/E/6kO9qdzqJ4m3tk4sdsn3cJYaJR9z6njn2WNFc3bLdTK037LU+RY7jYl5re5vRWvJiJ9G4/K6crHHNzVdtOND/jH12hvrt1FfqbHoP5y85pprVj8k7n+azBtHE1/NDtpxb50czbfR66J7L4veQ/64F+s+uqfX6ry5fVJ5HXWfZ19ksbW3jq2dH53P4stRX2ido/PaLPmf/Ij6pptuOuZzZOklCffff3/3ve99b7Dho2bPlm+GzuPvOvNkXsvVcoOl3urnhH5tbPG//jzaNdG4PsbvWuUsTbcpamPvHtyyR9DiQ9RGQ/dfrPstTaYxY8+9jozsdSN7uta5zhNjxpw726YhbXFqx/u/EOj/Kr8ELS8aPvWpT3Xf//73Bz9TlBfD8m3aW265ZacT1xsUhxJ1Lq9W+FgarCxOonVvRh701mXYVugAABDTSURBVGTX3gaVru2/8UULaKXjabzSm11Ev9/97nerN0flwb70a8vSPJo9Spyjemm21fwnMu8YX9Hkqelj9aexhWMeX/o3MzT7avYYKmZq3axjWFvinHyb+MQTT9w51VuktOShxcLS8SgjzfejNz1aMZlKr2jsKsXkvp9F5fasndJGMcJeu6akdzSHafNFNtst8l5ErprNon7WSv/cPt/+9rdXjcPWPy1PabFDO57L0Xq91DYxad7onJpeUebRtaTJs84YrslSO675W8lv99LW0fVtzSNJ53499uSTTw7WLHJTJK+lNHvU6g3t2r2yZe4Hms9EdVhXHujbt38j2StHNKZp85T8Vbsuksdq+cmy1tP1mmw1Vppftar3vHsNa+5O5w3x0rgMzaFdQ33otUzsDbWaX0biXTRmWO4RTb1O8vE1H/VaKJrbW97fitaSERaa72j85uSbmi7a8SFdp9BP5onYaki+qK9oLHK9073+oQaniB7a3NrxPgft/NLxue29avdf85q9X+9PYQNt3beq17RaOJ8nWje13rNEc6fmp32mkTpOk80b4zTfiua9mv9onFrn6KG6+eWXX+5eeOGFXV9REXtIo9Dzzz/fyRu8Tz75ZNOnyWr77f6+pfTcbajpQLN16x/XDslayw3eGBLxd/Ytu7/wkz/P1dauxz6ePbhljxBd49p1lhha0nvKsedeRyYm1n6BSy65ZOftk57nziX2rda+x6c95+5p01Ae6K+77rruscce2yV7gnf66aevGndOO+201XH5xMTNN9+8+rZ6KbFEGlxk7HSdfDZI/moNJvKLg3Re6VNGliSpnaMFvJyj8HjiiSd2WMnYclw+h/S3v/2tu+eee6r+kXdv5jfmcy6Wzyn1F17JTmk+GfOXv/zlSm6R9+GHH141Yslfq6ah2oOH9D1CmfPw4cOreTXuWuEuG5q+PUS32267rXvwwQeP8S0tUGvHS4aN6pHGK80b9RVNnpqDatdGZNIKb5EnfxWxd73J+VdeeWV36aWXdpdffvkxn/uTGKJ9QjGiV41jLXFrft1/sOPloSUozc9b+qM2V/SBc0smEdtremk2ltglN8ZSDhB97rzzzu6+++4bjMlDNo3IPTSOd/1E2Jdyg0XvSA4bI2PpFzpemwrrft6LyFVbz5ablNaaJerTyT59Oc8555zVP0lcPnTo0K4mhDH5X8ubJV7e9eJdF63WqGYHmSfKPLKWNHnWGcM1WWrHa7XNXG3tjTmWeDrkp3l8luP5fmJMTS1jaTp4460WO1rYMjHS6uGoP64rD2ixMiKHN47mPuCteyLyaTZp5R/5mvHsQ4WHdmNv6lyiMcrn9/IaYzOvf+S50HOPY4yM3njluS/ilctrG21vOFW8i8QMLc63qvNqOaq/DtIr9a01daR2j8YVbU1r+3pP/e71UyvjGq85+WaU9V7oF7GVtrY8vqKxyn+ALM8iXnrppcH7dRE9tLm14309tfNrx+e094rW+1PYQBjn4w49q9L80fKJE60WzudY956lFhciuXNIfm+tEM1DNV0s98VbP9eJ5GDNV8aySbnswgsv7F599dXVfvvpp5/eETVxkuca8uzC85Y4LUZF79e08kPLOuvzF59If9qznHSe19/HxDb2Lf/Zaeyw1MZe2+xF3dQihpbknnLsqerkyPqwxN7+2pFr+s9NvLHH618R3aI1S42JHGvWNCRJJG/6efHFF3fmTg9njjvuuO6hhx7aaWjpJ7bSGy9SQD711FNXjRZStMtffhOmr2jkple/OJT/LjW75AapnSfHtCSpndMvoIXnb37zm12vA+w3Owmrk046aZXwU2OTJbHn48gY559/fveXv/xltUjkv1977bViI9WQs2l26MstuiXfkQ7mf/3rX6u5xTfS25c0ntqDGXlImJjIfIlRqTkp+V+fe22eIXvkvjvUZDZGL2vwk/M8etR8M+orFn/WklnJJhGZLBvpvr08660fK4T/u+++uxPH5rYuIxtSDw8tKUXXQWvbi5zRDYxc24rJFHppMVK+YZ3+LDF5yKYRuYfGiayfCPtk6yG9v/a1r3VPPfVUN/RAJpLDIv4RXRfJj615L8JOi9cWblrNEtVfrrv44ou7o0ePrmoYqYvy+JtqSdF7377/+7WKlqei8pRYeddLZF305/bOqdWqafwo88ha0uywzhiuyVI7XvM32YekB4KphltHDTFGn1o8LdX2ljjyhS98oXv88cdXp+a5Sf67f8NOk7/mz9q1e2nLfJ3t37+/i+xRtLW8rjzQWo5ITMtrPW/d4+Wk+VWLuJ50iOxD5Vpt396q3ovsNfpzR3h5bab5B/Xh8L2miG1KOSCPE63jXTRmaDXi1OukVs/JMa2mtrC21u5T3d9qVUtqLLS4XGM1N9/UdNGOD+mq+bo2pvd+aeT+dcRXNLmFRf4jiNqPZr15RZtbO963k3a+xwaW+z7afNG9VxpXPoF0xx13dLfffrup3peTWtugv7ZLz4SG1oz3rRRaLdyv+/uxWbNHdM9S+1JHJHcOyRGpFaL1bSmWa3VvZD7NJtEcrPlKRNa+f6U6q9/01h/b2iijyZzX2Z43DfXXvfXZpWYb7bgnN7BvGX4WrTGuHZf8MsV97UgsKq3jqeqmVjG0Vuvl+aUVk6l4RHK/VtNb/cubAyMso3VNanLt99doupeON2saKv2yOJ+4/+Dd0gkl58gbDh555JGdRg9JCPIquKFfh6f5IsVaujYVDvLftQ5x63laULQkUnm70tVXX7168CV/klB/8pOf7HrIJazkjULyv9RYJcwPHjzY3XrrrTuf3dKcRea64YYbVs068idjXHXVVd0DDzywmjd9yzR/wBYtwuS6/nzya4677767S6/9atk0NDRf0u/IkSPHvG2gxl2z65tvvtnJmLnv1ubSxtOO1+w6Ro/avFFfsfhzSR/tWq9MNf3yGNVv7vGstyFfkDgmm2F5gGb58+pVG9PyiwrvhtTDoyab5uct/VGbK3rTI+nXionX9ppe2nEpUOQXVSnfaDG5tlZb5JLI+omwL+l90UUXdfJwtvQrbm8Oi/iHZjPteMmHhvJehJ1185Gf55lH0690PMW6oVouf2vlUPPmVHmz1XqJrIv+3K1ji4w/lnkeM7TYo/nFOmO4Jot2vOZvc7S1pk+rPNKf5+233+6+8Y1v7OSnCy64YNcPUpKPa/LV9l7atdrxqW2ZdJwyRnni81Q1nYzrlcMb05LsUX/1yKf5jcjSYq0nnbz70Dx+ex6Q1fbWpXsHJRYWRrm/RXh5bKb5B/Xh/+4TDd1rithm0+9leOOTZz1E/DpSu3vjiqbDFPW7h4Umn/UeTek+rTb+Ou9raLJE7r1pY2rHPbYqyRfZd2hyyVx5A7w0sdTe2O/RQ5tbO97noJ2vHffex9DGi+690rjpWcPvf/97U70/xX2VVH+mH21s6puGUv2mPWcaWnut8lmkVvDmoVoc15qGSrX5mOc6U+6bxrBJz1BLP3D/3Oc+t/rai+fLI7X9duIQjRkt9x+e+wf5w3zLD8DZt+x+Fq3lEO34VPe1I7FonXVTiUsLuaccW8szmr214576y1rTDz3D6j838ebAiJ08uln6ayz6989p1jQUmXyqa3JYng7YqeRhXAhAAAIQgAAElklAK1SXqdXmS92/+TjU6KzdpNh8Sm01hHlbnoxme1srnCAAAQjsBQHqw72gzpybToBactMt3E4/fKUdS0aCAAQgAAGdwJJr/yXLrluGMyAAgdYENrJpKHV+vvfee6tXVVreiNMaLONBAAIQgAAEILB8AmyulmlDS0NQ+nWZ91dCyyQyvdQwn57xts1A/N02i6MvBJZDgPi0HFsh6XIIUEsux1Z7LSm+stcWYH4IQAAC20XAknfmSoR9y1wtg1wQmCeBjWwaqr2idJ5mQCoIQAACEIAABOZIgM3VHK2iy5QayN95551Ovul74MCBXRfJK0XPPffcTr77y1spdZ6WM2BuocQ5HgLEXw8tzoUABNZJgPi0TtrMtS0EqCW3xdLj9cRXxjNkBAhAAAIQsBHIv2pz1113dYcPH7ZdOJOz2LfMxBCIAYGFENi4piGC4EI8DzEhAAEIQAACCyBAXbEAIxVETN9ml8Onnnpqd/7556/OfP3117ujR4+u/v91113XPfbYY8tVcmaSw3xmBlm4OMTfhRsQ8SGwwQSITxtsXFTbUwLUknuKf1GT4yuLMhfCQgACEFgsAa1Rde6KsW+Zu4WQDwLzIrBxTUPpVXF8amJejoY0EIAABCAAgSUSYHO1RKv9n8wvvPBC98ADD3TPPffczj8ef/zx3cGDB7vvfOc73RlnnLFsBWcoPcxnaJSFikT8XajhEBsCW0CA+LQFRkbFPSNALbln6Bc3Mb6yOJMhMAQgAIHFEUhftTnzzDNXbzI/8cQTF6UD+5ZFmQthIbDnBDauaWjPiSIABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEZk6ApqGZGwjxIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JoATUOtiTIeBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEZk6ApqGZGwjxIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JoATUOtiTIeBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEZk6ApqGZGwjxIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JoATUOtiTIeBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEZk6ApqGZGwjxIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JoATUOtiTIeBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEZk6ApqGZGwjxIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JoATUOtiTIeBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEZk6ApqGZGwjxIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JoATUOtiTIeBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEZk6ApqGZGwjxIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JoATUOtiTIeBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEZk6ApqGZGwjxIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JoATUOtiTIeBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEZk6ApqGZGwjxIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JrA/wP0jwJhmonrBQAAAABJRU5ErkJggg==)
%% Cell type:code id: tags:
```
```
%% Cell type:markdown id: tags:
![image.png](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABdsAAAA0CAYAAAByxuFVAAAfC0lEQVR4Xu2dW+htVdXAp1AdERGMo09REUgYkQ+mib740FOSUaIHzeON8AKlouSFEi1NkiC8oIKJN0o4IF4eDMJLkCiiPqSCEIooCl4RRJIi4fsa+2P+v3nWmXOOMeaae++19v/3fynPXpcxfmPMMcYca6659vuf//4F/iAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEGgmsB/N9mZ2nAgBCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQWBCg2Y4jQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAARGEqDZPhIgp0MAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEaLbjAxCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIACBkQRoto8EyOkQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAZrt+AAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAYCQBmu0jAXI6BCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQIBmOz4AAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEBhJgGb7SICcDgEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCg2Y4PQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAARGEqDZPhIgp0MAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEaLbjAxCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIACBkQRoto8EyOkQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAgdk22//4xz+GCy+8MPzhD38IJ510EpaEAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCKyNwGyb7f/617/C0UcfHb71rW8FabxP+e/dd98NBx10UDjggAOmLKYq2xNPPBF+8pOfhI8//jh8/vOfDzfffHPYtWuXet7UDpiTPd58881w/PHHh+OOO27yfj7GzpuuZ2ns5Hxx01mM8ZO5nRtt+ZWvfCX8+c9/nkwOwMfm5knIC4G+BKYam/pqOc2rrTv+bkotP03rrkeqnjZdt39aCc5pLmfVKT1uLnZo0W0q50zJh0r2xg+m4i3bU45W/5vS2IqWW1fd27XZ/tlnn4Xrr78+3HjjjeF73/ueuTko5+3evTv85S9/qXqyNKv//ve/h507dy6O++Uvf7lY2S5NjCOPPLJ67rPPPruQSZrzpaZHXC1vOcaqX7zvUPa5DdnooJ9++mk48cQTw44dO8KXvvSlcMUVV8xKlbnZozXIzcoo/xV2k/UsjR15iCJxZBgbNpnF3PxyrLzrSuya3PiYRojfIbDZBKYamzab+v9pt874uym1/HbwE6uOvW26Tv+06jy3uZxVL5rtLaTazpmaD9Fsb7MjZy2XQEs+mNrY2phm+9tvvx0uuOCC8PTTTy90sjaj5djYpH/mmWeyHvPvf/87vPjii/s0paIxL7vsMrXpG1fCv/fee9nmfNrwrzXGpcF/2223mbevKTX54/2efPJJ08OC5Q4l/epz2bZH42p56KLT6HeEJm9LkPNKp8mQu17LOTW5VqGnl0uv40tjp+SLm8yiF9O5XGeqDa1N9rHesWkuvoac/Qlssi+tKzZtMtPogZqO64y/c6nl+4/mzb2i16ZT9k+rlaY2l7PK7TlunXHCI2fLsZoPtlzTe87UfIhmu9eCHL8KAi1xaGpjK3JaV93bZWV7TPSytchpp50W7r33XlezXXOWeP1hAz820L/4xS+aXtGPjfJccz5eSx4ayF9uL/h4zEcffbTXCntN/lU0K1tk8JzjfcjguXbPY6eQwD36aPK2BDnP/eVYTYZV+O8q9PRy6XW8d+xsMoteTOdynXUldo3PJvtYSzzTePH79iSwyb60rti0yUzjKNF0XGf89dYj23Pkz0trr02n7J/zIr9cadcZJ5arWdu8c9kyrfv6NNvXbQHunyOwSXFoXXXv6GZ7bEAfe+yx4fe//3146aWXFo12z8r2mnvXigKtYBheNz5pyckWf5N94J977rms/D2N5JV93SHAW8ytS965cdXkXUWQ02Sg2T7Om71jZxU2H6cRZ1sJ9MwZ1ntajttkH2uJZxZmHLP9CGyyL60rNm0yU5rt2y9GTEFjb42pjcFNrg+mYC+rDJtsB80HrYw26Tia7Ztkzc3RZZPi0Lrq3tHNdnGn119/PXzta19beFatod3ietqrCJ4iI0L+z3/+s8/KdFk9f+mll4Zbb701/PznPw+5D9rFFfbpyvh4zfPOOy+ceuqp4ac//eliKx3Z9/2BBx4I//znP/f5wGWUOcdjuKL+b3/7W/jd7363tT2PvD1w5ZVXhnPPPdf8sT35SMEdd9wR7r777sXHTeVPHo786le/Uve6jzoPZY1b7bz22mvFhyslp06ZnX766eGWW25ZvEkgf9/4xjfCTTfdVJRLPgJ03XXXLbYVkr/DDjtswV14yPcCZIufGlctaAhv+fCrbO8T/4SV2ODwww/f69Jj9JALWfwglfeee+4Jv/jFL0ysrDa3yDDkaT3HKoNcv2QXGf9i39x49IyN3LGnnHJKuOGGG6rjyKPDkJM2dnKxocYiXt+j91CmVp/1cNDuIWNdxtJw6zGxhzywzX1IeozOPWUvxSaZONx3333h2muv3fqA9M9+9rPwox/9KJxwwglZ/83FKY+stXxai5PCd9lxJfXjUm6Ub6+06Ds2B3jHkZbj4zdkhvbw6raKGGXRRXw5zTNDveJihVqt15r7lxkbtLg0HNvWPNdqN+Hq8RGv/NFuY2OTxWfkXpYYbWE61q9qdbjEIk/N2WJbi46t8dfKOZcbrPVILV633H/Pnj2LbT5l7hHnL+ecc06Q+7zzzjt7vZXcYvsx8dzjC3KfWu558MEHw4UXXhhK25rm5o+lHG6dh2g2zeWmqfpnGq8s+Sc3Z9AYl363xK5hTeP1nZytqdP2pdLaA7HmqfSOWnPbmy/G1vSaPMcdd5z6LURrzrL6fNRpu8TxdP7pYWTl3hIH0hzXqz/lqbVr9Yr04C655JKwa9euvVTTfLllbHnsMbburc2vW37r0mxPb9yz2R6ffMqHU3Pbush9teSaylZ6khr//amnngrPP//8wnGGe6mXzo0O9d3vfjc8/vjji4LyiCOOCK+++uqioZ9rqInDPPTQQ+Hhhx9eHH/SSSeFAw88MHzuc58LUoTGpm66PY98lFSOiefUPuKa6vzCCy+Ek08+eXEfaRrLwHjrrbcW+knhq31cVs7/05/+tNBFGtzxGvvvv3+46qqrqm8yaBPuyEw+uirbD2lyxSJR5B7yEP845JBDVK61ZntauEabyMOE+B2CoQ8ObW/VIw2cmh/Eexx66KHh/fffD5Z7eGxu9cXUpyzneGRIi9i0mIixJOennrGRXif6jfjahx9+WN1+yqvDMABrY+eDDz7Y50FciUW8tkfvXEJo8Vkvh9I94liSYubyyy8PV1999SIeSUyJ8SUX18bo3Et2T2xKY1m0gSVee2UtJXwtTkpsW3ZcSf24lBvfeOMNd27SdLPkgLRuyOWSoa20HJ9raHhtuaoYpekiNYa83ScPwmKuj2MzbhMo9pQHSC2NsXXGBm/ss+S5VruJD3p9xCt/jA/pmGmJTZrPiP9bY7SHae4NVGtNWarDPTVnq20tOrbEX2/caq1HSvHaY2fN92RRzcEHHxw++eSTLs12q//l8pKlppbztNzz7W9/e1HP5RaGeFbxeuYhWo1ZWriwjHnHWP+U89PtXLX8k5vLleJDeu3hd9tW4Tu5Wk3zJ+q0/+/vLKtO0xqCnnyRs3Frjh821Wt9i1z/LeWVm+t6fH4Y+9IaYhPjuIw7byxrrRUseWUoj5wzpj/VMg/P1SvSS9mxY0e4//77F7yGD5l7j61ePuuZk5fm2C3/Pulmu7aqXRT2NvdzzfmY4GWSLb/XjhnuDx8dSor7XEOl5HBa8VUqGuS8a665xvyRVuEjK6IvuuiivVZml/bBLzlR6Q2CMRNuYSYTKll1H4vCklzxPhLc5YFD2uCQyY00WeQhhZXrMJmlX05+7LHH9mIlqw/kgcXww7mp7a16DPla5fWy8thckyHnE9o5Xr8bjpP43zIRGj4Q8owNTc5a0PTq4B07WjIqFVzDyZwnJrT4rJdD6R7pw9Nhci59wNpj6xz/XrKnxdewCZTGDmlMxtgkDUs5Vv7X0mz3ylrSV+6pxclV+0FOf6++vXOAdRxpOX6s360yRmm65Gqfknw9c/8qYkOLz/fcxnDoJ17/b5G/R2zSfMYbozV/H+tXtTrcWkdpMmqTLe38Flt6ObfWI9pcxho3W3xvjO2tcrWwt+QeWTiwe/fukFscZp2ntsxD0oZNaVHaKuYdPfzTk39y9XNt3OVs4JW5xXeo03ZuIWiZq7eMa0u9LUJpczBrvijF2tYc39Js13JOlNHr89stjksfycPIyn1MHGjJC8uYh+fGQ8k/eo4tjz3SvvCwZ+edk2t1nuf3STfbax80jUpai5haoInXiE9mck5Suk88Nrc1TS2YawO0tj1OyYk9hvdeYxnN9twKkNKgsm4XpHEt6a1dP/d7SdbU7jkdUztZ5fWwKvlBSXdNhpZmu1eGVLbf/OY3ixVCuUa7NrEY6tiimzaOeo0dLRkNC64eMaGHzw5jqfWhQK1hnf6WThh76Jyzp8beM95qearGW/MzjXPufC2OaQW3J3Zp8mm50RsjrLppY97rU616WP1Ok9fqJ+lxmn+X6hVhI+Nv+IAzNkHSsdmzMbaK2NAS+5bZbPf6f4v8PWKT5v/e8aT5+xi/0upwa1zXZNTGpHb+GFvmGqqemqRkr952bvG9Ftt7/a83+9QXSouFrLlLO670u3be0F+n6p+e/FPy+dLb7jlGq/Ad6rR8tNR8sNU2pRxgrce0GniZNb1Wt9W2kdF4Rv29XLdjHPcwsnIfEwe0+O7tT9UW2+Zq/ZrPpwtl0tpE82VrLebt+aTH57Z06zF+tfov9/tkm+1a4ReV8Tbb4wrKjz76aGvf9mGCzx2Tc8A0KJeCYEuDM3Xe+LpIajxJJvLqhvXprZwre4PKqktZdSF/8RrWD9kuo9meYxbZp28QeIKZdmzOHto5aSMgt1+/VY/cANTuXZtE5Vil97DaXJOhRe54jlWGqGfcLkf885FHHtln7/6WsRFjhLxaJ/uCyzcZvvzlL5vjpVWH0gW1ya1lFUOL3p7iUo7t5U81nx0+2MxNVOMY66WzJ/55x5s2dloSe6u/abKkrL16tsQVawPIoq9HN0uDVPKgNbda9ciNN4tuct6qYpSmi6fZ0dIYW2dsaPF5ze/G2m0OsanGrSVGW5m2bCPjrcNreW+MbTUdvb7YwrlXPSLX8d7/C1/4wmKV93A7zihTKS96Y4pXruH3SmRel/6NnX/k5o65f2utqT2N5FqBO1X/9OQfrZmTNl1zNliF77TaeThOWuaY1tpDqwk8cznNr1IeU6jTNB9q4T6mDmyRJ72flrO8Pr8d47iXkeQUjfuYOGAZU7m80KvWFtm1GJF7INPiy7X8a52ztfqsuSHUeOBkm+21p2m54GJtGsu5afPrBz/4waIoHH6oJ3eM7OmebhNgccKxzfaaXS3Ndhmolo/NaP6zzma7tVhNJwVakZ8mUcv1c08CvROmloDbcg+vzS3BfCi7do5Xhqjnd77znbDffvstHgrlnkqmidAzNoYf1pBtNuIHQkrX8erQa3JbeyAUH5aV7qXFhFX4U68k32rrYUHviX9ePto48DTbx/qbJY6NmcR55dMKNM/1PLpZJ3HW+KHpUYrrHr+Ta6wiRmm6eF7j9zbGtHqp9UFcDzuWHjRq47vVbrFW8fjIumKTtdlutYPGtLdfaX5Xe8jcMiZb69DIT5twWjn3qkdSfaz1R+uE12v7lhrBO448uWc4v5QHutYFYZb7jF05H31CG4NeRi12yPmnJ/9oc+t0O5+cDVpk9nLJ6WixM3XaZ1tbMlnjnVbb5K7TqyHYaw7ZIs/w3rWc5fX57RjHvYziFsjeWsEaByzHeftTnlpbq5/kd0+z3xtDvfZo9dlanOnx2ySb7aX9e3MKWwuZ9Nz0nLhlRWkPb2n4nXrqqdkPGVqcUCsIck1hrRDyGD42yaUJd/vtt2/tRe5NTOtstnt4aMfWGpmlBr03mMjx2irhZRW9abFvtbnGLOdv2jlev0vtcuedd259pG/4yrR2X21syGoP+UbAjTfeqH4g2KtDr8ltq49qumvxquSzXg69m+21canp3FP2llVvnma7V9ah7p6x4S14WuKKlmM8+np0szTbPT6l6ZHzQY9uw/OXGaM0XeLv6QfV4wdSS98q6LUC2TMB8PjD2MZFbYVuajuP3XqPp2XGJkuz3TOeNNt5G65aXtN+t9RqXttqOnrjr3Y9LQ+mv3vftJNzvffXji/lRa/ttfvkuCyb/VAHbQsA61ygNA9JY0nvPdutK3tb7FCzjSX/WOrMmJtyNmiR2es7LXOn9JyW+3lrD60m8FzPw3QKdVpLc9uSL6INPexqeUqzUc7PcjnriCOOqL5x5J1TbGIc9/iwlfuRRx65z6HW+1iO8zS7RRBPra3VT2n+sez84I1pFv1TuNrxnjm5p67Sjp1ks93z8c7Sa3U1xdNgeemll4Yzzzxzsco1fvU3dS7ZV+iqq65afCCzZa/ElmZ7S/HkTere4K0V6J79l7yDzcPDOtA8+2GnT9Yse1KJvNaE3Cpv6R6167X6ote3WmQYypa+ljXcP9g6eanFgNK2UJaJT++xY9lGxjMGanp7x14PW6byeJP8GFv3ll17GFHLD2O+3eDxNyuvdftBi22sumkx1Xqd6Lce/nJOi265MbuMGGVtnEq98+ijjwZpesj2W1deeWU499xztz5kPqyRJEbHFT613ywNktwbTbk6r6cdx6xsH9pOs1urj3jHrJYzrJMOzf+9dtDGZ02u0m+ajC3sWsakpX7QJq8teUab6KW/a7V8aSser51rx2t29MwnvHK1+ILnHulKxFjvpNti1mxVu09pHqKNc2/9vi7/TOOCJ//UHgjIdrE1G3js2solx996X6+vttQe1ppgOA8rnWfVTcsD1uu01mk1e3q5e8eYt8mv5TfPXNfLdTvGcS+j1lrBep+WvNCz1rbEiOGHwb0+XuppteS3Fp/11FAtx6682a6tWtd+Hyppddb0vDTIy3YSr7zyyj4fBLMcoyVf7fea7LFAKG0LkX7Vu2T4UjKTf7/mmmvCbbfdFqzb75RkLdmr9tXflkQWJ7A5Hs8880w4+OCDt1btWwZa6U0GaRg89thjW9cStk888cTiYcvw3i16eAsv7z1abd4yjkrntMiQ07Nkc8/YEP8866yzwtVXX72XTTV9W3QojUPv5LZkc4/eJVlW4U89k/wYnVts6OUjnEsyxmtJw1Lb3qdF1lqBp8VJr54t8lkKtOGK2Fpu6pUDvD7lneh4Wa0yRlnsftBBB4WHH3447Ny5s1pT9s793gdxPe1Yezheq31acotA9fqIVkuW5O8RmzT/99pBmzz19qsWdq1jMh0wLXVobcLZwtlTb/a2c4vvtdjey8USA4fNcU/uEebx+KOPPjo899xz2S0Rc7aJunjmIdp48talLeOlVgfF+2tz1mh7a/7R/DWOP1lMd++992ZtsArfoU47YAuBZ67e0zbe+VBLjBjeo2eO13w95g1LPeLluh3juIfRmFrBmlda8kLPeXg6l5WeoewOEBfXePtjLWPLY4/WObm391ydHGV+7NJsv+uuu8LLL7+8uPxbb721+CjOYYcdFo499tjFv8kqhYsvvnjx/1Oj5VYC1pzPGswsEOJ95NhSI8RyjBYEa7/H68vKsdNOO23BTrasiSvsY7EgMgpLYRo/bCr/ZnllML3HiSeeGHbs2LG1Yk2uMbbZnhZ6oke8h3zAVeR9//33F/ZPV761DLZ0VUe8z4EHHrhoCkgzK2VR41q7d8o7fjzvtddeC08//XR2y5EWPWqFV84PWu7RYnPNF71ye2Uo6RltMvRT69iIAVQe/sQxFG2qNUC9OpTiTq9mezrWWmPCKvypZ5Ifq7PXhi18hrEpxnLJg/IwVx4EfvLJJ3vFQG08tcZqa5xs0bMnS9Hfez2rbsNrj82tWo4fa8tVxihNFylqpQaRnJr+SdPoq1/9arjooouyDy175H5vs90bG1p8vuZLJ5xwwtZ2Z97c0uL/LfL3iE2az3jt4BmfPfwqnXdYt8UYMybjuGmtQy0Pflrzf2qr0tu8pZXtXjvXfK80N2iZT3jl6jGOJD+X5h+pv8U3g4argku14lAXyzykZtPafabon578o8WlmE+EgcSRkg2sc4mWOFLib61lWnzVW1dpHL3Xs+rmyQOWeKfpkbNF6ZwW7mPrwJp/WXTz5CyPz2/XOG5l5OE+9BHPWOnZn/LW2tH/vv71r4d//OMf4dNPP93qWco8V/6s9UTr2LLaQ2RpqXu13nQtl1p+G91sT5Uq3TBtasUVz++9915xNfnwdYTSdSMc7RX9WpCVYuiSSy5ZfDR1+PfCCy8sVjTLMblXm7Xka/n917/+9WL/aPmTJ0WPPPJISPd3kqdG1113XXjxxRcXx0jRcNRRR4Xf/va34Zvf/KbFxmHPnj3hiiuu2JpAy8RQrnnGGWcsHjRIMtX+tNdY7rvvvnDttdfu9dq5XP/4448Pw1UirYNNfC29j8gsRbs0B4avuJe4akkrx/uHP/zhPk0GzbbWbWQi9xZ5a/dosbnmizkfqZ3jkaFkF9FRxqA88JBvDuzatWtLDOvYkJhzyy23LJJBHEOnnHJKuOGGG/ba+iCnn0eH0hjq2WyXe1j1rsU9a9MhXsPDoXezfazOvWSvjTeJTddff/1WLJc4LduPnX/++dkPcJd8xSNr6RqWONkagz3yabFW5PdcT4636KbF1Pi7dRxZ9BgbO1YVozRdYqySh0RSZ8S/uG/7sGkxtEfccqYl93snAD3tqOXrUp4bYzev/7eO2bGxSfMZrx0s47OnX7XWamNtK/ftXdeNzYVyvrceGcY2a9yM8XqYF2+++eZwzDHHLOYGubmb1/Ze/xszjqzzj3RerC3qyOUOzzykZtNSjaCNwVZGY/3Tk3+0uCQx/eyzzw7ST9BsYPXpMVyGtrDUMq3389RVGkdvnppTnbbsZruXXYs8qV95cpbV57dzHLcy8nBviQO1HNfSn/LW2tEvzzvvvHD66aeHCy64YNGTkT+ZK9x000179S1rNVdrTPPmFm/dW+tNa3nU8vvoZrvlJss6puYwy7on14UABCAAAQhAAAKbRCDWU6W33Vq+j7NJfNAFAhDoQ2DMQqk+Eiz/KtpWhcuXYF53IP/My15IC4HtEMexMgR6EJh1s12KGVmt6nlFrwc0rgEBCEAAAhCAAAQ2hUBsdvz4xz9evAmU/qV76Vu2r9sUJugBAQj0J7DpTZr0I6nyVpD2/Yv+hOd3RfLP/GyGxNubwKbH8e1tXbTvSWC2zXbLq089QXEtCEAAAhCAAAQgsIkE0r0n4/dDRE/5Tsyjjz662CJOex1/E7mgEwQg0JfApjdp4ltA1m9i9aU7z6uRf+ZpN6TevgQ2PY5vX8uieW8Cs222x2KGVVa9XYLrQQACEIAABCCw3Qi8++674Y477gh//etft74TIwzkA6my/+33v/999TsX240Z+kIAAj4Cm9ykifu1y4fjeOva5xfkHx8vjobAOglschxfJ1fuvXkEZtts3zxToBEEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQjMlQDN9rlaDrkhAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEBgMgRotk/GFAgCAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACcyVAs32ulkNuCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQmAwBmu2TMQWCQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAwFwJ0Gyfq+WQGwIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCZDgGb7ZEyBIBCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIDBXAjTb52o55IYABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAIHJEKDZPhlTIAgEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQjMlQDN9rlaDrkhAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEBgMgRotk/GFAgCAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACcyVAs32ulkNuCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQmAwBmu2TMQWCQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAwFwJ0Gyfq+WQGwIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCZDgGb7ZEyBIBCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIDBXAjTb52o55IYABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAIHJEKDZPhlTIAgEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQjMlcD/Aq/jswf4sjnzAAAAAElFTkSuQmCC)
%% Cell type:code id: tags:
```
```
%% Cell type:markdown id: tags:
![image.png](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABYYAAAA0CAYAAAA+LOytAAAeDElEQVR4Xu2dbchmU/fAjydvCRkZn4TUfCDRhKHxxQcfZJoRGjKMtzQoZkQYMQ1mTKTkJRTCjJcaKTOK0ngpIjHFUL6QkilmaEoTKer5/9f1tO/Z95m9z15rnXO9nPv+3fX0PM9c5+y99m+97LXX2WefA/77/38VfxCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIDBrCBxAYXjW6JqBQgACEIAABCAAAQhAAAIQgAAEIAABCEAAAhAYEKAwjCFAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEJhlBCgMzzKFM1wIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhSGsQEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAArOMAIXhWaZwhgsBCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAgMIwNgABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAYJYRoDA8yxTOcCEAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIUBjGBiAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIzDICFIZnmcIZLgQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCsPYAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIACBWUaAwvAsUzjDhQAEIAABCEAAAhCAAAQgAAEIQAACEIAABCBAYRgbgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAwywhQGJ5lCme4EIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAR6Wxh+9dVXq5UrV1bPP/98demll6JJCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQUBLobWH477//rhYsWFCddtpplRSJJ/nv119/rY488sjqsMMOm2Qxi7J98MEH1Q033FD98ccf1UEHHVQ9+eST1eWXX168b9Iu6JM+fvrpp+q8886rzj333Im38zZ6nunjzPlOyhZnOos2dtK3e4MuTzjhhOrdd9/t/RzQN/7jlLdvfoytjtNa6HuSCeAb49NO3+Lo+EjRcx8J9Gk9KnxHHQs9fEYVM0bVT1u79jBs22cX9/eFbxjrqH2jC8apNjotDP/777/Vhg0bqscff7y68MILzYUsuX/Tpk3VK6+8Uu3YsWMgrxQgL7744mrVqlXVySefPG0M991332DHsCy4zzjjjEZGn3/++UAmKSTnFuhhF7LmGu34Qr9SFP7666+rY445Zli6HGq7weD/+uuvasmSJdUhhxxSHXfccdXq1auH2m/XjfdNH30LjF59zeRx5nxHCv4SR+qxYSaz8NpHX++bKYlCX/mPU+6++TG2Ok5roe9JJoBvjE87fYuj4yNFz30j0Lf16KgLw14+o4oZo+qnjV17Gbbps6t7+8A3HutMyRM6Kwzv3Lmzuvnmm6tPP/10wElbOA1QpSi8fPny6r333hsUg5ctWzb4acuWLYMdqqnCajD4u+66q1igDDuMd+3alSwkx/03FXGlGP3MM8+oj7DIFaRDfx9++KGqsN2Vo3nb6cvRHSWumgcEXkae+0ryjiIwlmRIjctzTxOfUYzTo58u7sn5Ts4WZzKLLnj2qY2Zkij0ifmkyNo3P8ZWJ8Vy+itH13nBpJAYp2/MVKb1tV9uLdS3ODopNosck0+gq/XoKGPEKGOht37SVcwoce2qn2Faalc2NkwZc233gS+F4Yz2QuEjFHQ3btxoLgw3GW8oxtYLwKHYe/TRR6te0821I8MKbUmBW/5SZxeHa/bs2dN6928p4IzDCZv6tBbExyV/37iW5B1FYCzJQGG4nTVbfWcUOm83Iu7WEhhlEq2VietGQ6BvfoytjsYuZnIvnlyiDzzG6RszlSmF4T5YPjL2gcAoY8Q4Y+GoY0aJa99yvD7YcqrQ2pejNCfBN7rQcesdw6FYunDhwuqxxx6rvvnmm0FR2LpjOBSXU7t/Q9G43mbJaeuAcu3IdeE3Obf4iy++SMrfpdKtsneh7DZtWItbbfpqc2/fuJbkHcXEU5KBwnAbi6wqq++MQuftRsTdWgJdzhnaPrluMgj0zY+x1cmwmz5L4ckl+jDecfrGTGU66iJPH+wMGSHgITDKGDHOWDjqmFHi2rccz2Nb47ynb3wnwTe60FfrwrAI8eOPP1YnnXTSQJ6m4muTwE33NRWNLUWXoLR//vlnvx2/0scdd9xRPf3009Wdd95ZpT4WlJIjtHnjjTdWV1xxRXXLLbcMjtOQc4rffPPN6s8//9zv42FB5hSP+k7ljz/+uHr00UenjuiQXdn33HNPtWLFCvWHjOTg8eeee6566aWXBsdyyJ8U8h944IHi2cxhzHVZw3EbP/zwQ/ZBQM5JYmZXXXVV9dRTTw12aMvfKaecUj3xxBNZueQjXuvXr586g3revHkD7sJDzreWYz6auJYCjfCWj+rJa23hT1iJDupnXLcZh7StsYNY3pdffrm69957Vay0OtfIUOepvUcrg7Sf04vEBdFvyh8tvpG69rLLLqseeeSRRj+yjKHOqeQ7qdjQxCK0bxl3XSavzVo4lPoQXxdfqh8/JPqQh4upj3S2GXOXsudiUzgff926dVMf57z11lurSy65pFq0aFHSflNxyiJrbj7V2rqV6ebNmwdHNoWPj8o8dP3111cy9/3yyy/T3tppms9LyZNWrpKdeeeR2P60sqR04Y3d0palX6vNeG1VK1MpJ0p9Z8E6Bq2Nt/UxzViEZzwv1/sMGxo8PlGy8WHG0lLfdf/S5gVe3VlsxCp70JnXN2Kda2xG6+Mapm3tqmndIrHdkqNbdasZ36jiKPPp/z50bdW5cGtal3U1n0o/lhiQ87F63q+1WcvaMNhSiUvTelQ7Vo0PedYObWKhpp6T2jxYr7Gk+GjG2yZmBFbD7kebT+XikjbvSDFsqqlJf7nftTJ759/celV2DFvqH9r5NfSn9bW2eYKHX2p+Hva3yjopDMfK9BaG46Mcbr/99kr+IxOKBNalS5dmPxpXMvBYttzTn/Dvn3zySfXll18O+q6fd5W7NzjA+eefX73//vuDRfPpp59eff/994Pic6r4I8bx1ltvTZ2ffOmll1aHH354deCBBw4W2qEAGR/RIR98k2vCmctNH8iLx7x9+/YBP5FLCpxSSP35558H45Mic+nDfXL/a6+9NhiLfBAwtHHooYdWa9asadwhXioMB2byQTs5U7okVwjUInedhyyQ5s6dW+TaNBHHhbygEyl8h3Oz60X7uu6144gn6ZIdhD6OPfbYavfu3ZWmD4vOtbYY25TmHosM0nZKLyGWpOzU4htxO8FuxNZ+//33xiNorGOoT2ol3/ntt9/2e2iUYxHatow7lVB4bNbKIddH8CWJ63fffXe1du3aQTySmBLiSyqutRlzV7JbYlMcy4IONPHaKmtKv1pbtzKNY288PnmQN2fOnGrv3r2dFIYtcnlsWZiV5hGJ/XFynJpvNPr0xG5rvx6byemyyVbb6KWeE9WTWusYtDae8g9rX6X8TnIyectMHnKF3CjEsnC0muQ58nCoTQGvniuNIpZa/UuTF3h111ZvmrypHhu8cbxkM2L/Wn+yME0VWrQ5eG7dYsnRPbrVjG8UcTRXfNGOSavPetEpxIiQXzCfHpZTRWWNARrdaa6J52T535q1oTbPyK1HLWPV+JA1r2gbC5uO3Qz5RyqHkt9knR3qEik+mvF6Y4Z1nevtxxov6k4R18tKeUeKYW5ekH5y3+OyyGzNHXJOPwq+Fl8rxe9x5dDZoOn8YWIKwyJ/vINMJsz58+cPjnVo2klmLUSnCsnBESRQye9N19TPMw6GK4lVKtDlAr/2FYX6Tkm57/7771d/AE/4yJOWVatWTdvxGsaoPfIjtzO7zWJHmEn/sps5PFXOyRX6kcRJiuPx4lImCikKS0Fdy7V+Zk1oX+TYtm3bNFbh4UT9o4Sx7rXjqPupVl4rK4vOSzKkYkvpHqvd1f0k/H9Z0NUfXuQmtZRvlORsipvWMeTayvlOLjaU/r1NTPDYrJVDro+gC/nAqPzFxwblkhGLrlP8u5I9TqzrMTOOHVIUCrFJ5jO5Vv5bU0i0ymqNJeF6K1PP+NrMC1r79tqy6KQ0j1gZpWzPI5+1X6vNeHRplamUE9VZWcbQJp5Lv5a+5PrSWFK5Yk7GNj5Rn/9HEUs99tuknza6a6M3bW7m8Y2S3zetCbRxrsStrV01yajNO0syNuVZpXs9dmiNWZ4cl/l03xvCo5hPLTGgZFOiO801Yc6QGGJZG2rXq7lc3zJWzVis/tBFLEyte+KCZn3DT+qbUW3rJ9r4lYtPJRsZV2yy5B0phk3jSs0nVvvxcCnNpdp53Cqr1dc8vmGVqZR3Ouu96tsmqjBcLw7L/5fgIa/2L168OPmKsbUwnFJQaCMUKlKOlOsnXJs6niJeWNQLkaWA03RERi5YqrUeLXi0h3oPozCcOh4g50DaI0NKXHPsSu2nfs/JGus9NcZYT1p5LaxyduCdZNskzfV7czLE//7QQw8NdtKmisLSnsU3SnwtPlNfCLT1HQ0LSQDCn2XcJRsYpj01+UXTA6nwW7wzv4sxNyUcdR16fDrImDofv6k9re1p473W1q1MPePzFCuscrXRVerjsrE+rLI02ZjF17rotynv8OjSKlMpJ2pj91ob1/ZRiumlsdR3OIX2UrHM4xPjjKUe/xpWYbg0n01SHNfYjBx5lopB1kW8cGljV6V1izZ2tfHL0r0eO7TGrDY5rrUvTwz26NgqV9ecu55PLTGgZFPSluaa0noj93tpPVmadyxj1YzFagseG63LnMr1gx1fc8011caNG6dtDmkqSFrrJx5b9sQATz9WXaTksuQdufVD7m37lHxWmT1cyKHPHWxKjfP33PxszW+t109UYTgYqjx9lPMLjzjiiKnzdVNPJEtJUQpG6hWHupOlrkkl+k0LsFLg1yTQsrMuvLYSj0WM5fXXX1ftQgv3yRkqspst7NYLbYxzx3CqsJZ6aqidwDUTpCfxlnZTQbSpYJMaR5uJR8sq7kOrcwvf0L72Hq0MgWU4MkPsc+vWrfudNR3vlNL6Rkg45CGTvH0gZ4gff/zx6lipHUOuwS52DHvG3TTRDtOemvyi/hAulrHuY12NWfrQ6tDq0yU/8BSGtbKm9FuydSvTgw8+uFq+fPl+RyvV57d6AcG6kLXKJbt3utZVPba1nX+98nn61diMx1a71ktT0NWMIc75vPG8q3gg7VgWaFafKOWXw46lVvvV5F+l+FSalLU2YpXd4xs5WZv69vhTSbau7apkd7nc1qvb0vi8uvTE0bpOS2Oy6pP59H/HJrZZz2pjQEl3mrmkZJuptaHmnnr+lNtkoh1rl/UEr43WfSeV+0qOL99ykiMUb7rppip+A7tp85W3MOxZ58TjKOlyXLHJkneUNiDFxcdU7csa4zy5uWcubaoVeWK/xte09hCvh7rmV8qRuvh9YgrD8dOl+GgBgRqOTki98mTdMRyS+fDE/qKLLhosfOsfz4mDVLhGziCOXxUuJVBNv2sDeZOSNa8nSz+aD6OUjClX3Oo6KU05e9N5RXW5tY4bTxia9lNPP62TQoqxR97QTi5Bt+q8JINHbqsMgeXZZ59dHXDAAYMHGKkdmHGQtfhG/dB3edAUPtyTa8c6hlw7XReGLePOJWuWhMnKYRiFYe+Yu5Td89DKUhi2yppj0mTrVv/xLhKs84JVLk/yqYnzcXErPETNcS7Nv9b5wcPAYjOlOJ+yVY9MTePOzSXWHMUTz4NuLX2VxmJ5pdPqE6X80lsY1sZSq/3GvlP/Rkfcp0d3FjsvcRt2HNcWhrV6KPlt13bl4RfG4tWt5uGjNmfxxKwmXTCf7jt3NC7kjXo+tcYA0anGHpuu0YyxvjbU3BPsLRcrrGPtsp7gzflKa3H59lKotbz99tuDjxyGmspRRx01+K1eY2ni02XMaFoDdtlPV7HJkneUGEqeG95eSc0lHpk9uUNKB9Z2PLJafK00F48rhy7V76y/T0RhOHe+ZL0AtmvXrv3OG/UUhuN7wmvruTNnpTglXwWUV9tTiUlp4VByylQCXTI+i5JDUUoWsc8+++zU2bkluet9jLMwbOFRujY17tI9wqJPO4atOteMvzTpn3HGGdMuscoQ6+WFF16Y+qBP/XVLj6yxYOHro48//njx44vWMeT8ssvCcNOCWxMXrBOttGnl0HVhuM2Yu5R92AUFq6wlfadsPSTnWqYlf8sVvq3FilI/XSaNpbF7ZBmlfHFfFpspjaspqS0xi2Xy5haeHMUSz7uOZXEBLf64b/j4XO4scs9HwlK5p6cw3JUemx5KNy2kRzUXW+c5j2/k4q+mMGzRQ0k2a6wtFX5Lv2vehrP4ZWl8XeuyNG/mfmc+XTBth2dJb4Gj9rqSXixzXb0tjT169Zt7003j4znbto5VUxjWyCPcSvqybHaI2YQ6ihSEV69ePbWmljXemWeeafoQt1ZG7cOknO113U+pvZIPhN+DDjR5h2Y9FnKS1FrVI7M1ZnvmUs96LNWPxddKLMaVQ2vtRnvdRBWG9+zZs9+O3FKgyp2T0gQgNih5rUHOu6kXoGIFr1mzplq6dKn6TDDNIqlkYNoziprG2dSHd/HWxEk+EhY+Ihcnl/XXjD1BQ8ujxLU0EafOfYufQsW/e8ZR15dXXmnHGhg9Dyk8E6bH7uqyxa+B1T8+p7WFJt/IHQ2jSWa79p3c+Yj1f+9i3Fab7UKXsR4sxYy4kFM6Ezal665lzy2INedvtTlr3GpvdRZ1W7fakWd8TYuH3G9Wuay2bLEnqywp+xumfKW8yDPPjUIvlvnPYveleO7lpS3ySX74zjvvVLJQkyMu5Ci0FStWqPIhb6407Fjqsd9SPpPyk5Lu+hTHY13mXg+3xpYSU0+sLfmWR/ce3Zb8ssTTMydrF8e565hP963ztLasvc6z3ijZctxmKdbItRb95taG2vGmZPfEu1KM0MoTWHlyvqYcSPLfeh0lHvsFF1xQrVy5MluHmSlHSVhyUI0vWPKOpiK51N5CPhG/FaCxh65y367asdi6x9c8vmGRqTTntZ2/NPePvDCc2h0cB9fU6+PBYKXoWD/KwQq8nojIK+XffffdfjuRY4PJXaNRYNOk1SR7GHPudVV5/WXu3LlTO4AthREZWzieo+0Zw7nd3jt37qykbfnv+hg8SWeYrFM8Pvvss2rOnDlTLDSOm9sh3vTl2S7GkdKTR15px1IYLunc40e5e3LBtkmGlE3kdG7xDWF07bXXVmvXrp3mK6XxesaQC7i5vnJ+kPt3y7hzslh9z8OhqQ9rMaPNmLuWvXS2Yi4+SLGodOyAR9a6jrW2bmWauz7oOTU+z7xglctqy/GirzSPWGXpKqm19OuxGY8uLTJpcqKYlXUMWhvvKh/S2NiRRx5ZbdmypTrmmGMac26PT4wzlmrGbllEenVntZGSDQ4zjpf6lt+t/lQqJnRtV6Ux5I578+RZmsKDxw49jJlPf8q+GZvzGe26rK0+rDFAE2sWLVqkWhvE9Ydt27ZNW0d88MEHg01j9XxCy8VSGG6zhrPy9+QJTXOuHBEh9QDZ5BNqOEGn8tv8+fOrr776ynRUZ5fr5qaJu+t+rLpIxaUFCxZU2ryj9OAkjC/1QcDQt1VmT8wmh95XzyvprHTKgqb423RNJ4XhF198sfr2228H/fz888+DD9XMmzevWrhw4eDf5GnRbbfdNvjf8WIy3gUWDE+ukR0Xy5YtG1wfXsuTf5NzgeWDU+Gv9JSsaeAhcMs1uUW75pqSApt+D+2H8Qo7ed0ijDE4rMgoLIVp+Gic/JtmF13cx5IlS6pDDjlkameLtNG2MCxtBDllHKEP+TieyLt79+6B/uPdxJ6gET88CP3IRwxkQSaFiZhFE1dNoV7GFD5s9sMPP1Sffvpp8tgBzzhSNumVt5SsxfoIu5lyOi/ZolVuq93lWAbbqtup1jcCI3lAEXwo6LRUrLOOIRdvuioMx77mjQkem7Vy6LKY0XbMXcpeen1azuqKY7nMg/JQUR5a7d27d1oMLPmTJ1ZbbF3rPyJnPfbK3BzP86kY75kXrLr22LJlHrEwSunTI5+VgdW+m3TZZKsWFqWcqM7KMgaLjXfhY6WxSN4qOZvkIPGfLNxOPPHEatWqVckHkl3kStaHbFbb8tpvLpeQQoxwGcVc7JHd6xt1OyvZjFUPcn0pP+syB5f+rPy69Mswv4S1kFWWoA9LzErFCsuYLH0xn/rWs13PE//5z3/U8SjWr2ZtqM0zcrZtGWuwXW2M0KwduoqFcexKrT1jrk1HLKV2u3a9bs6t34bRjyVepOSy5B2l+ahee6u/oeuJp96YbZlLm440svC1+prXNywylXSWq6PmbNj6760LwzGkXOdxASbsJE2dFyy/PfXUU9Ubb7wxlWhLAn3xxRfvl2DHyUvpNd2mRZsk9LfffvuguFn/2759++BpoFyT2slcSqA0vz/44IOVnHcqf7JbdevWrVV8Xqs8kVy/fn21Y8eOwTXC46yzzqoefvjh6tRTT1Xpe/PmzYMzfcLiRQpk0ubVV189KIqLY5T+mp6aiQ1s2rSpWrdu3bRXKaV9OZu5vqvEGzTq/YjMUnyWBLL+2maOa8nhUrxL9tf2HCMZh0fepsDo0XnJFlM20nSPRYacXmSM4oNSnJczsi+//PIpMbS+EeKKPDwIPnTZZZdVjzzyyLRXfVPjs4wh50NdFoalD+24m+Ke1WYtHLouDLcdc1eyN/mbxKYNGzZMxXKJ0/LqnHx5OfVx05ytWGRNtWGxdYsdpcb35JNPVuecc84gxqfmYOu8EMajlWsU84hWli59zWrvVpvx2qqWRWmObRtnLTbetq/SWEJsl6K65GXhL97QEC+0rD4xzljq9a+mfKaN7ix27pXd6xuxnZVsxhrnwvVNuVaXdiX9efi10e0wcmBrHGU+te8YFmaWdZl2DukiP9LYo+aaJn/NrQ21XJr8zBLvNDHC6g9dxMI4lqTqKHFRMrXRrRRLu1435+xuGP208QVL3lFiKOua6667bvDB99JmKa3MnvmDHHpfPa+ks6Y6as6GLf/eujBs6azra5t2THTdF+1BAAIQgAAEILCPQEhgPA9n4QiBPhMI+WfurauwEyW3qaDPY0d2CECgewLMp90zpUUIzCQC5B0zSZuTOZZeF4blqYk8ZcptfZ9M5EgFAQhAAAIQ6D8BFrL91yEj8BEIC7Qrr7xy8KZb/BefB6k58ssnAXdBAAIziQDz6UzSJmOBQPcEyDu6Z0qL0wn0tjBc2mqNoiEAAQhAAAIQGB4BFrLDY0vLk00gPoM0/i6GfAcinOlfejVzskeIdBCAwCgJMJ+OkjZ9QaB/BMg7+qezvknc28JweE2P3Rh9MznkhQAEIACBmUCAhexM0CJj8BL49ddfq+eee6766KOPpr4DIW3JR9bk3L7FixcXz7H39s19EIDAzCLAfDqz9MloIDAMAuQdw6BKm4FAbwvDqBACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAR8BCgM+7hxFwQgAAEIQAACEIAABCAAAQhAAAIQgAAEIACB3hKgMNxb1SE4BCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEfAQrDPm7cBQEIQAACEIAABCAAAQhAAAIQgAAEIAABCECgtwQoDPdWdQgOAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMBHgMKwjxt3QQACEIAABCAAAQhAAAIQgAAEIAABCEAAAhDoLQEKw71VHYJDAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEPARoDDs48ZdEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAR6S4DCcG9Vh+AQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABHwEKAz7uHEXBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAIHeEqAw3FvVITgEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAR8BCsM+btwFAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQKC3BCgM91Z1CA4BCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAwEeAwrCPG3dBAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEOgtAQrDvVUdgkMAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQ8BGgMOzjxl0QgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABHpLgMJwb1WH4BCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEfAT+D9MUi9pPSD3FAAAAAElFTkSuQmCC)
%% Cell type:code id: tags:
```
```
%% Cell type:markdown id: tags:
![image.png](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAByYAAAA0CAYAAADPGF+9AAAgAElEQVR4Xu2dbchlVdmAt1FNhEiG+qfoCySKKMKaQv/0ox+RNFGDDn71RagDNcr0pZRoaZIIoSMVTFGN1YAQpEGC2BQYRphQUxBEEojBqMWASFNQ8L7vfXrXtJ79rLXuj73PPmef53r+ZHP23ute1/2x7rXuvdc643/+76/jDwIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgMASCZxBYXKJdHk0BCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCCwIEBhEkOAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASWToDC5NIR0wAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIEBhEhuAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASWToDC5NIR0wAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIEBhEhuAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASWToDC5NIR0wAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIEBhEhuAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASWToDC5NIR0wAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIEBhEhuAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASWToDC5NIR0wAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIEBhEhuAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASWToDC5NIR0wAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIDDbwuT3v//97sCBA903v/nNbu/evWgSAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhBYYwKzLUz+85//7Hbv3t29+c1v7qRIuc5/Tz/9dHfWWWd1L33pS9dZTFW2Y8eOdR//+Me75557rnvRi17UHTp0qNu3b59637pdMCd9PPnkk9273vWu7qKLLlp7Ox+i503vZ813Sra46SyG2Mnc7k26fPWrX909+OCDsx8D5sZ/J8g7t3iBT+wEq6SPyySwah9aRcxZRZvL1GHr2Zve11XOwVbZtujc2/6qfV3zAW9/tOdpv2+6b2j9X7ff0ce6aQR5IFAmMHWsHqqHVYx9EUZTxMAp2hiqr3R/hOFYbUefMye+0sdV+EaUrfe+UQuT//73v7vbb7+9u+uuu7r3vve97kKK3H/vvfd23/ve97rjx48v+nL++ed3Bw8eLBbAvvCFLyy+mJQF3wsuuKDZ91/96lcLmaSQWVsgTl9hWq6x9i+1K0XJ3/72t90555zj1dFaXJ+c4NSpU92ePXu6Xbt2da985Su7G264YS3kswoxN33MLVha9dC/bpP7WfMdKThLHOnHhk1mEbWPud63ycnDXHWyaXLPLV7gE5tmgfRnagKr9qFVxJxVtDm1XlN7m9zXVc7BVtm26DbS/qp9veUDkf4M9alN9o2hbFZxP/pYBXXahICPwCpitU/C7VdPPfZFGU0RA6doY6i+ojnOGO0OfcZc+PbnCJv40cNohcm//OUv3f79+7tHH310wc1auEuQpSh51VVXdQ899NDia7zLL7+8+9e//tUdPXp0cckVV1zR3XPPPVtsLwWRz372s2qBLH1h+cwzzxQLmXn7rSKiFEO//vWvm7eQrRVEU3s/+9nPTIXVoU439P65bJ2rcbUUqIey8tyvyTtFsNRkKPUnck+LyxT99OhlzGtrvlOzxU1mMSbXOTzLk1iP7VNz4IOMwwnMLV54fGI4HZ6wUwlscjxdtQ+tIuaM1eY62IUmw1h9XUffX+UcbJVt54t2rZef+zpbpa9rdrqK9Y0pfEPr9zr61bJk0lhMoY8p1iCWxY/nQmAdCIw19mnxYMy+Tj32RcezMWKgxnWMNsbUTe1ZY9nZFLLmbcyFb5J5at+YUh+jFCbTwnsqKB45csRdmEzPkILm4cOHT295JwVP+bdSQTEVG1/+8pebtslLRcVSITM9S9qTv9LZlemakydPDv76UQtCUxqBpS1vQdbyzGVcMzeumrxTBEtNhikmBVP0cxn2Znmm13c2mYWF1yZd40keIn64SazoS4zA3OKFxydiRLgLAl23yfF01T60ipgzVpvrYBeaDGP1lTgwfwKr9HXNTmt0o/dZtDWFbyxTfksf1+kajcUU+phiDWKdmCMLBNaVgBYPxpR7lWNf3g+tz2PEwCnaGFM3m/asMXQ4JZN18Y1l9HlwYTIV6y688MLuq1/9ave73/1uUUj0fDGpOWQqWvYLitp9fWCpkl+SLf0m51Y+9thjRfnHNASv7MtQvueZ3uKK59ljXjs3rpq8UwRLTYYpJgVT9HNMO/M8y+s7m8zCw20TrvWMGRE/3ARG9GEYgbnFC49PDCPD3TuZwCbH01X70CpizlhtroNdaDKM1ded7P+b0vdV+rpmpzXG0fssOpvCN5Ypv6WP63SNxmIKfUyxBrFOzJEFAutKQIsHY8q9yrEv74fW5zFi4BRtjKmbTXvWGDqcksm6+MYy+jy4MClC/fnPf+5e97rXLeRrFf9qHdC+RGx9GuxZ9E+KlC1i++c9SvHzU5/6VPe1r32t+8xnPtOV9u0tFUjTM6+55prusssu6z7xiU8strOVrVp++MMfdn//+987OUvuoosuOn3mZpK5xKP/peYjjzzS3Xnnnae3yJWvUm+88cbu6quvPv1VqWYYchCtfIX6ne98p3vuuecWl0sh+Ytf/KJ6Nmfqc7+NtN3tE088US1E1xwnZ3bllVcutuiVfsvfG9/4xu7uu++uynXs2LHutttu23IGqXAXHnK+qWyz2+KqBR/hfejQoU622E1/wkp08IY3vGHLo4f0Qx5ksYNc3u9+97vd5z//eRMrq84tMvR5Wu+xyiDPr+lFfF/0W/JHj2+Urr300ku7O+64o+lHnj70OWm+U4oNLRbp+Z5+92WK2qyHg9aG+Lr4Un/7b9GHvNwisaX/N6TPY8pei02SVMr5yLfeeusixkqc/uQnP9l98IMf7C6++OKi/eZ91HzqFa94RTXO5rsNyH/nf7WXejxxrja+WH3Kq7v77rtvsTV74ijj3cc+9rFOxtgTJ05s2R2hlW9oiZtVLs2eo+NVbudWWUq6iI4R8ixPux4/kmdHfcIqk5Z7yXneVhvtc/X2dYp2LP0V5nmO0O9Xeikv4jeaHywzrmtt931Qi6d79+5doInqTe712IhX/qS3qA/19d7Km1McivQnn9fkbWrtDbG/UptW2a124Y2NHjuyyDBVTG/lh5/+9Ke3xBLJaWRutG/fvm3DkJV/K79txTeZm3/4wx8uvjDsyX9a8z+rDq3jkzZW93O1qK9b5fHGIK+dpv5o9/3jH//oDhw40NWO4KnlrIlnyzfOP//87uDBg0Ubtfq0Jv8ycnHPGs+y9F2yV42FjKOriFWaXMvQkXVdiPnSk4s1T5mvjbm+F50vWeP6FHMAqyyR+bmWc7XGPuv4rfldyqutsXasPLeVTw7NDyx9HhIDhcGy27COGbXYlebzlrllyc60Mb21TmUZG725jZabedbavbZu9bWhvmHVuWVdwRq3Wvbj+W2UwmTe4JDCZO38x9YCo2bwuWy1NxLSv//iF7/ofv3rXy8S2/7Zj7V7k2zvfve7u5/+9KeLxdS3vOUt3Z/+9KdF8bNUfBAl/+hHP+ruv//+xfUSzM8888zuhS984WIBNhXA8i1y9+zZs7gm3WM9o+Lxxx/vLrnkkkU7UmCTxP2pp55a9E8mmg8++GCzOCn3/+AHP1j05fjx46ef8ZKXvKS76aabml/IaoXJxOzUqVOLM0U1uVLwFrn7PGRR7Nxzz1W5tgbnvJCUdCKF13Ruar9o3Ne9tR/JJi12kNo477zzumeffbaztOHRuUWGfkCx3OORobZYkWJJyU49vpE/J9mN2Nrf/va35hbQ3j70OWm+89e//nXbSwuthRv5zdPv1uDr8T0vh5pfJF+SCcbnPve57uabb17EI4lLKb6U4tqQPo8luyc25bEs6UCL15pPvfa1r+3ka/7SNuIpLpbakN8kbuVx3hvnSnZk9Smv7vIYn3OUF1bOPvvs7vnnnx+lMOmRKxrntfEqTeI8srR82jNGeGOJ14/yyVba4j/5UMsnPCy03Cu9NJXnC8uI+1ZfKI0PntxM66/khxIj5IWPlOuluJp0IHFfXpQYUhjqjx1TxHWvD2rxVHLsqN5Ej15/8MqfbKUWD63jSt8PS3mzxKFof0pFQkvcG2J//TY9slvswhsbvXZkkSGS93vltuaHYjO7du3qjh49urilX1Dy8G/lt634lhgPzX9q8z+rDj3jU4uv5jfWHNIjjzcGeew07492n+T/UrgovWyqfTWS20+e72g26vENTf5l5OKlNY2hcyGvvkv2qrGQcXQVsUqTaxk6Ej6WdSHmS/8pTHrWGKx5iseP5VprXB+am1vasVyT989jb5acqzb2ecZvze+ia9dD89zWx02e9ZESI0ufozEw2d0y2/DkCKXYJf+WH3WnzS1LDFs1nPTsfv3HI/cYY11Uhx45Pb42ZA7okUlbV4iuo9RsyfLva1GYTMnoQw89VHyLLn01VfrS0VsILRUyk2NIIii/t67pn2eZlCqFv1IiWRsMtAS85shy3y233LL4MrB0DmZf6cJHqv/XXXfdli/+Wm+RlAyn9mXqkAUGYdY/U7QmV2pHFqilOCtfQ+RBXYqSMihaufYnhun5Mml6+OGHt7CSN5FkATF9JZraznVv7UefrVVeLyuPzjUZSvag3eO1u76fpP8vhdh+8dzjG5qcrSDp7UPtWTXfqcUG7d/7E3pPTIjYrJdDrY08zvcXuGrJiUfXtQmhxxeifMT/+/EhnY8s/6sVJkV2zVZLdpQnjP0CfukM5kic88asdL1Xd7ls+Y4GLY5Dxh+rHw2xCW288jIq2XhEPm+73hgQ0aVXplbupfnSWHF/qnZEXi3XLOWtNfmG+E0/F5kirkdsvKWbIXpLC15zGFesebPXv2s5iqe92pEftTjQatOjC033njikPasWZ7T7IvbukbsmV6vdWkwfy3a0+DZG/lOyIU0X0XymNVa35p/W3Mer74hNaWwi6xv5eNFfw7Cs50Rs1MtK6/eYtujNQ73Xe9YOInP9iF159RGRawwdMV/6b965TLuz5g0eu9F8eMo5gFWWiL152ZXGnlXlTymXHmP9ZAx/j4xn/bnZOq0Fe/yl5Q+euaU3xyqN+V65I2NQv7+RZ3jl9ObKc15HsRQga9esRWEyD07y3+krGilE/uQnPzm9/Wh/4Te/z3qmZcmQkvLTm6Alx6olzOnaUtE0D1j9wUAbqGrFjNYzPYZQC8De4sqQBa7SG5M1R2/xyGXWuNb6rT2/9HtN1lxHpT5G5PWw0hYdvLYYmRR4Zcj18uUvf3nxxl2pKCnP9fiGZg8en+kvUNS2Mes/c6zCpKffGv9l2lPLLyzbauSLFWP02bNAFPHpJGNpa6rW8/pyabZaYpfir2xzduTIkS0v9pRicyTOeeVM13t1F+EYGX+8cg2xCe3lIa8sLVv2+PQY7bZykYguvTK1ci/Nl8aK+1O1k7Ou5Zqlr6PlvhQ3cluM+M0q43rEB5dZmNTG135uEJE/4kPW3MNq/7V8OZpHp3aH2J8174ouMnniUNT/tfuG2EtpzLHO91rttgpKkRyr5iO1+DZG/uNdNMv75bELb3yI+LpXnohNWe3UO6eszQO0/FSbW9ds1MtK6/cYtuiVyXt9RN+RuX6kHW9fInItW0faWkSSWbMl5kvtNZ1oDLZyt+ZCtfmOpR3LNRZ7KvmNJWa25mresUrri9e3I2NfSeYx/D2aM0ZioHddJdKGVxc1W/DMLWsMa7tbem269PwIm35fI88Yi682n/KsLXplmnodxRJv16YwKcLKlxHXX3/9tvP95NzG/fv3F7f+sLxhl4MoffLdd7rSNaXFHUuwjwS6PLlOW0fkfZCJm2yrY/kKJ90n+xrLm5jyVar8pWdYC7o1Yx97gaH0lY82EOZstGujE9NSUG1N9kv9iCTXQ9qw6lxjFpHba3epn2l7HrHPBx54YNs2wxHfSDYqX5TJOYZyhuyrXvUqS3xcXGPl2BrUS184exb9Iv0uyTOFPbXa6L8EksvY97Gx+uzRoZeP5jutZGeMxFSYyflLsjX4tdde2+Vf9PdjtiaryFNLHvuyaj7l1d2LX/zi7qqrrtq2hXpqt8bRO/545ZKXoca2idSniCxj+PSQdi2xULOzki4jMmmL7ZqNagOApa/yjKna0frrmTx6/UbLdZcd170+KPJqdjhUb+s+rmj9L9m/1eajeXRqc2z78+jCUrCWeZJ17hWxI003XnuPxE/vWCLXtxYRh9iOFl/y3/MXcDz5T6sNTYfL5Gu1hbzfEXm8NmWJoZH1DXluaY2ltR1fbquWcVDmW2kxL8IqohOPLXplSvmxJy5F9F2KCVYWpRdGWus5nr4MkSsaL7R+M1/67zpivhbptTsL5zwWeexGi+tTzgE0WSwcamsk/aPHImO7dfweM3+Kzvtb/Yv6eys/0HTjtfkh8WzqOJvyvv6RQHn8y1+G09Y18xfPSmO+d2yMrNF4/WPscczia1abW2ZeOMb6hhZjS7+vVWGy1gHLQre1wJZPrsSZ3v/+9y8WRE+cOLHl7Kp8UTddI2dQ5tuseCZTnjcKc8dsKdVSmJRnWQ6s1YxnlYVJ64QlOpGyPL/0Ns46DkaJgUfnWgCMDKJeu0ss3/GOd3RnnHHGooBeeksk6hv9w3tli0Xx/7Qnfq2PHo41H6r5jjaA5zEj2u++TBGbjeqylEBFF7CjcXBM2SMvTYxZmOz7qZxlnMauH//4x90111zTpTHqZS972eK3fMyKxrka+5ZPee01OkHxLnB75YokvRbO/Yl21L61HERLqK3tevxIG0+0wqRVJm2BUp4zVdyfoh2tv57tdrx+o9nZsuN6dNxqvewQtY9IjuWVP+JDfb+xxqGx+uNpb0z788QmbY4QGR8idmTVb2QRyho/vYszcn0txnjyZE/em8s4NP/RYtiY+YzkDVa+VluoLUBZ9e2NQZqvtHhqfZJ7+/Mi60vm3nEw4tOa/ENt0StTXphcpr5Lz9ZYeO3K2/eaL2lyDdWRZTzzHIk0ZnxhvlT3gv666BS5uTUHaMkSsTfLPYlU6yUSz/htLUxa4lTUji1xyrs+MmQ888ZAi/wXXHDBlsu8bYwVZ1t5X2l+pb2sJGu7qZBZGvMjcnvZWHOzdB3rKI90d955Z/foo48ukFjWz1v+r/02i8Jk6ysOazKbg8jvSdtG1s58kOLIZZddtthasjRJ1JJkzUlLb7poCZem1Pz3lPzLYP2Nb3zjdCFGk7vfxioLkx4e2rXRN73n9MWkV+cas8gg6pUh18u3vvWtbvfu3YsvqPvbUkVkzeWXN1UOHz7c3XXXXV3/TL6azS/LdzwLNEP7rSWn8nvtK98huhS/KcXeUtF5yNuAtZg4puyrLkz2E8Q0LklB8oYbbtiyZePb3va2bWOWxYasX0xqPpUmBZY3OeVZmmy1Aq93gVtrZ8yEVet7RJYp5YvmEVq/WoVJjVkukyeHmSruL7Mdrb/pdzlT6sILL1wciSAv0h0/fnxxlncei71+05q0y2+RwuRYuq6NW5odajGsvzAQ9QeN3bLGFU//xxgnPe2NaX9e2S0Lax7bjNiRxsq70KI9zzpn1GJM6YtJL39P3tuXO89PvPmP5oeprVIM9+YzNd6R+ecU42U0hkbWNxKbfgwYa0vCvo1GfMNyzxBbtDw/tyHv9ZqtW3d08uTm3pcoojE2cbEwWbaOmC9tj3TLGrss+m6Nc8vMzb05QHSMGbJGUtOLd/weM3/SdFqb99f0PMTfWzHTKqc1BpbkH7sN7XnWnDDnYplbtvzfMuZH5PbGHNZR/ktAy/m9sc1jV7Vr174wmRKoZ555ZvFFY3+xIJIY5EmZbIMn53P1CyB5QLzpppu6Sy65ZNs1WvI3JNBZk/SWEbQc3GOM0ob21ZfnfKtIELHy0IKaNjiXzmfJ3+CwfLIuvKyJf1TeWhsRnWsyeAfRiAx9veTbXvT93moLLd+obc1smfSM7Tu1s3b6/z5Gv72+N4Yucz14FrBbcccy+I0tu1a4bZ3tpJ01a1kAyMcTeV5/XMp1+573vKc7cODAtjGrZUO1OGdhLdf0fcprr5a98fscWxOX2m9eubw+47FbryxjJNQe+TS7jIynU+ilZrNTxf1ltNOywzzWSa6azmaXl29uvPHG7uqrr+7yLw8ifmOZcFpeOPHan5ZnRxfVSzai6S3qD2PHEOuCjSW+jDlOWtrrj2OS31ls01vUiRZNrH1YVn4Ztfcx5LbEmPyN9zFtx5JXt+bllvzH0kbS69B8xjNWR3Ifr74jMUibG0Z9LJ+/njx58vSLLflRBDX/8tpoZLzR+r3sXLzU9yn07Z3rrzJWrVpHzJfKHhqJM1bbtl63rLF5yPhhyTcj83Mrk3XNnyJjnzY2RNdHouNZxOb7fdDiWaQNq220/CU6t2wVabUx3yt3hE2/z5FneOSM5MoR3/DIpI3f0XWUlj1Zfpu8MKkVGnOh84G/tlWrVwn9hQX5JPUPf/jDtqJnbkS1ayxKbRl7S/a0gF/brlW2BDj33HPVrShLn1pL32655ZZOzmGwboFbk7WmT/naTZ4t/9vvQyQApEliiccvf/nL7uyzzz7NwuLMtS9kZYHk4Ycf3sL12LFji8L0GP3wTji8rGrBT9N5xI9q90RkKPWzpnOPb4h9fuQjH+luvvnmLTrV+hvpQy3g1tqq6bb2755+a4mb9c2uCIcxF7CH9Hls2WsL4DUZEwd5y8yy7bZlASX1SbZolfgqi7ppe/H8t7e+9a3db37zm21bjydZPXGub0tWn/LqLsIxMv545fLGYOFlHa+8sngWO/PFv/6Cn6fdiB9FdOmRScu9rDZqnRzWxs+p2rH0V3YZOOuss7r777+/O+ecc5r5f8RvVhnXIz7YiqdRvfXnD/lLU60cKyJ/xIf6SrfEIfmy1jtPqPXH0p5snz+W/UVikzbOeuLQEDuKzFPGiulafig5i+QXsrtIKhqX5kER/t68t7QmEM1/Sm1bdeixC41vbf5Zm2OWckivPJEYpPlKdH0j8UnxQsauxx57rHh8R20B0Wqjcr+XldbvfByI2qJXJu/1UX0ve40iqg+vXGPqiPnSe6tz2P5aZMTurHmDxwescX2KOcDFF19sWoeKzM+t7DyFySHrhR4dtWJBZP1k6PpIdDyL2Lw3nkXa8OqiJFPK1a1zS614n/Je+SDsyJEjxTHfK3eETS2vsK6NesexSK4cmQOOzc6yfu6p6zUXJv7/x1EKk9/+9re73//+94tHPvXUU51s0SCTXdlOSv7k7YXrr79+8d95sMm/MEnGunfv3u7MM8/s5IDUo0ePLu6pLepqbxi0AKRg3nq+5RrNCVu/p+fLW+2XX375gp1sTyMM8sRY/jttzZVzKX2h0+9z3saePXu6Xbt2nX6bXq4dWpjM5ZR+pDZEd2IDzz777EL/+VvRkSCSF6lTO2InsggnE5ScRYurZeCRPiU7fOKJJxb7Kpe2/Yz0o2STUXlrhZKIzjVb9MrtlaHGMsWFvp2mf9d8IzGSAnnyoaRTrVjk7UMt3oxVmBwjJkRs1sthzAXsoX0eU3btyxz5kiCP5TIOykst8tLE888/vyUG1mzF4of5uNTyC+1lHmuc68vq8Smrn0ob/RifxsSUT5TGksj447WpiM94xisPo5LdROTzMvD6UUuXLZ/wsLD0e9lx3+MLQ3MzLdeUCYnkj5IP5X+y6Pua17ymu+6664ov54yRt3m/hPfan0XXpa9tavFUForSdvHevEBk9/pDRP6oD+W6t8ahsfpjbS8St7ViaG7H6Yvh2hxHG2etcWgs/+/P/yL24vWp1lhy3nnnLeZup06dOj03lXFY/vrzzrFsR4tvSd4h+U+tMGmNBVa7qOV3tT5Gfd0jT9SmovNTzceEUb4epB2tkZime17/+td3f/zjH002GvENi/xDbDEi0xT6Ltlu1AZaO0d5+sJ86T9roXnsZb7UXhcda2zW1k0t7bzgBS8w55u5X1jWIa051zrmT9q837t+kufmpdwvZ1taH4mu10fH1taccF1yQpHRM7fU8rg0V5TntsZ8z/gwBv/oMzxyenPldckLtXWUWl2vNm5r/z64MJmDqzWWFwDSl3T9rVmlY7Ktapr8yLNkwWD//v2dLCK0JlGWbfL69+cgDx48uNiqtP/3+OOPL76Uk4We0hZVeWJdqrJbfv/Sl760OO9O/uStrAceeGDLdrXyluptt922OCcoOfLb3/727itf+Ur3pje9SdPv4vf77rtvcQZZWrASrvLMD33oQ4uirziL9td6u1ds4N577+1uvfXWRRtp+zB5vpzN2V8wigaAfjsisxQ/ZTGuv1VZjasWNEu8P/CBD2xb1NN0a93KNXGPyNtqI6JzzRZLNtK6xyNDTS/SR/FBKQ7LGan79u07LYbVNyTm3HPPPYukPvnQpZde2t1xxx1bthAr9c/Th5oPjVmYlDas/W7FTM9bQd4YMnZhcmifPTocEptuv/3207FcYqBstXrttdcuvkY5ceKEqTApfdX8MB+7+uNSnvS1XlzxxLmSHXl8ymOvEuP7HA8dOtS9853vXIwlpbHeO/6k/ljlGmIT+bjYGq+ssozp016/8viRPLukS4tPWFloY7nHRvtcPX2dqh2tv2mckUHMBcMAAAWHSURBVMm85IjpL50z2Z8Eev1mlXE96oOteDpEb2OOia08LupDuT1b82aPzbf0YW1vTPvzyK7l3N7xYYgdjZ33e2N6P+7lepU5prxQnObkElckp5AvXofEy5rtaPEttTkk/6m14dGhdXzyjtVRX7fKs4wYqulMy2Xz9SPtRdG+/uVs9SuvvHKxRiTzQ/kTG7377ru3Hfnj9WlrjBhii1GZptB3yXbXLVatWke1dSHmS08u5mreNQZr3uAZ4zxxfciYZmnHck0rJrTszcKuFatXmT8NmSu2xtjSur22PhIdz4aMrf0+rGOc9cwtNYYy3/noRz/ayYv82pg/5Vg3RIdWOb3zxiG+YZVJ05clbtXqeiX/tPzb4MKkpZFlXdN6S3tZbfJcCEAAAhCAAASmI5CSp8hLSNNJSUsQWA2BlAvXvpJOb2rWXrBbjdS0CgEIrBMBbZFinWRFls0goB2rsRm9pBcQmI4A86XpWNMSBDaZAHPLTdbuevZt1oVJSWjlDc78nJf1xIxUEIAABCAAAQhECDDRjlDjnp1CIE0er7jiisXuAPlffl6MZfv/ncKMfkIAAlsJUJjEIqYkkL4WP3ny5Lbz0KeUg7YgsEkEmC9tkjbpCwRWR4C55erY79SWZ1uYZAK1U02WfkMAAhCAwE4iwER7J2mbvnoJ5OfcpHNu5RlyJnk6b0/bNsfbJtdDAAKbRYB59Wbpc917k77kr33pv+7yIx8E1pEA86V11AoyQWB+BJhbzk9nc5d4toXJlNDyBvjcTRD5IQABCEAAAnUCTLSxDgi0CTz99NPd4cOHu5///OenzySXO3bv3r040+N973ufeqYyjCEAgZ1LgMLkztX91D1P50vKGabsejU1fdrbZALMlzZZu/QNAtMSYG45Le+d3tpsC5M7XXH0HwIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAJzIkBhck7aQlYIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIzJQAhcmZKg6xIQABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIDAnAhQm56QtZIUABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIDATAlQmJyp4hAbAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAnMiQGFyTtpCVghAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQjMlACFyZkqDrEhAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgMCcCFCbnpC1khQAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgMBMCVCYnKniEBsCEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACcyJAYXJO2kJWCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCMyUAIXJmSoOsSEAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAwJwIUJuekLWSFAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAwEwJUJicqeIQGwIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAJzIkBhck7aQlYIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIzJQAhcmZKg6xIQABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIDAnAhQm56QtZIUABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIDATAlQmJyp4hAbAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAnMi8L9sNo73oWQnggAAAABJRU5ErkJggg==)
%% Cell type:code id: tags:
```
```
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment