{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "e3876bdf", "metadata": {}, "outputs": [], "source": [ "import sympy as sp" ] }, { "attachments": { "dipole-gauss.png": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAIAAAC6s0uzAAAABGdBTUEAALGPC/xhBQAAAAFzUkdCAK7OHOkAAAADc0JJVAgICNvhT+AAAAAJb0ZGcwAAAAAAAAAAANoqts4AAAAHdElNRQfmBwsGMhnkfoktAAAAI3RFWHRTb2Z0d2FyZQBQT1YtUmF5IDMuNy4wLjgudW5vZmZpY2lhbMtjldoAAABWdEVYdENvbW1lbnQAUmVuZGVyIERhdGU6IDIwMjItMDctMTEgMDY6NTA6MjVaClBsYXRmb3JtOiB4ODZfNjQtcGMtbGludXgtZ251CkNvbXBpbGVyOiBnKysKcjwBAAAAIABJREFUeJzs3XecJGd97/tPVafJcSduzkG7K6EsIQkJGQSSEGCBbGNsfC+2D9g+XPvidG2fl891OPj6+hzucTjwOtiYZHIQlgAh0JGEclppk7TaHCfs7sxO7OlU9dw/aqq3eqZDzSi0Wvt9v/aFhuqnqp+q7qpfPxlERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERM6LVDsDIrJYd8EmeKna2RCRRbGrnQERWay3wnXVzoOILFa02hkQkcXaDDPVzoOILJYCsEhtug6aobna2RCRxVIVtEht2gQxiMMnqp0TEVkUlYBFatPl/u27ocoZEZHFUQAWqU0rwQJgdZUzIiKLoypokRp0B8TBAgsa4LZq50dEFk4BWKQGbYM6/+84XF3NvIjI4qgKWqQGrYaIXwUdgZ4qZ0dEFkElYJFaswWWBGaxi0A/bKtmjkRkERSARWrNrdDoF38BC+rglmrmSEQWQQFYpNb0QCxw79oQhe5q5khEFkEBWKSmbIP+QAMwYEEElsKmauZLRBZKAVikpmyBpnnLmEWgQTNyiNQYBWCRmrIVEmD7g4C9fzbUw2XVzpuILIQCsEjtWA/LIBqof/ZYEIXlsLE6+RKRRVAAFqkdDVAPEb8E7PFKwF4t9NJq5k5EFkQBWKR23O4H4Pkl4AjUazSwSC1RABapHf3+ACSrsATsFYJjWphBpJYoAIvUiM3QBNESd603GrgNrnm98yUii6MALFIjbgzUP2cASEEK0jAJWb8ZeHl1cykiYWkxBpEaMQRPQcr/vw40ATDt/5BuBGde87CIvFHpZhWpNaugDRrhBsjBvTADbfBCtTMmIguhErBIrTkKQBdcBzaMw0B1MyQii6E2YJHalAYDFjRWOycisigKwCK1qR1scOFAtXMiIouiACxSm475f7RXMxcismgKwCK1qQ8AA+eqnBERWRwFYJHaNAjoDhapYbp9RWrTGnDB1QpIIrVKAVikNsX8YfzZKmdERBZHAVikltmQrnYeRGRRFIBFalPOHwecqHZORGRRFIBFatMhsNULWqSGKQCL1KaNYABoqHJGRGRxFIBFatOkXwU9Ue2ciMiiKACL1KZG//ZVG7BIbVIAFqlNB8AB4GyVMyIii6MALFKbOsACA73VzomILIoCsEht6gTA6CYWqVW6d0Vqkzf6SAFYpGbp3hWpTWfBBgty1c6JiCyKArBIber1hyHNVDsnIrIoCsAitcmbBstRJyyRWqUALFKbBgCwIVPljIjI4kSrnQERWSw78L8iUmsUgEVqUxsYcCFZ7ZyIyKLox7NIbaoDB2wYrHZORGRRFIBFalMXWGDBhmrnREQWRQFYpDalwQIH9lc7JyKyKArAIrXJW4XQgsYqZ0REFkcBWKQ2RQCwoa/KGRGRxVEAFqlNM2ABGsogUqsUgEVqUx0Ajr8qg4jUGgVgkdqUAAMRaKt2TkRkURSARWrTEbDAQEO1cyIii6IALFKb+gEwcKrKGRGRxVEAFqlNA2DAQGu1cyIii6IALFLLvFpoEalBCsAitanXH4akm1ikNuneFalNLWDA0lSUIrVKAVikNu33S8DrqpwREVkcBWCR2rQMAKP1gEVqlQKwSG3yhv9GYKDKGRGRxVEAFqlNcbDAhfXVzomILIoCsEhtSoIDBuLVzomILIoCsEhtmvI7Ye2tckZEZHEUgEVqUwdYYMGaaudERBZFAVikNqX8ObBiVc6IiCyO1vIWqU0JPwDXVTkjIrI4CsAitSntV2ANVTkjIrI4qoIWqU3T/h/D1cyFiCyaArBIbWoBwIKuKmdERBZHAVikNh0CFww0VTsnIrIoCsAitanHHwccqXJGRGRxFIBFalOjvxxhpto5EZFFUQAWqU2HIeLPRikiNUjDkERqlgsRaOLLX2bgOIMDTIzz+a9UO1ciEo5V7QyISCgf+Qhf/KL352a4kf46GtqIOuyb+c53jjQmjnU0n01EjoydYdce9u1n334eeKq6WRaRchSARd7o/uJPuewK2rr4wheu/9znboEOaICoPxm0Cw4kb33HyK/c9cT29T+wHZwMUxO8sIe7f8J9T1b7BESkGAVgkTeoay/jtz7Cli00tJB01959z3v+4q+3QCPUQQxsvw+HCy5kIQXjf/jxfav779u0an93G5kk01O8dIDf+Jsqn4uIzKcALPKGs2U1H/8Q115JQzMplxcPX/LPX77zwUeWQbMffaNg+/evARdyfgyehBMXb/jSn/7xwZ4c7c1EDaOj3PsIf/ONKp+XiASpE5bIG8sffJh3XU9/P1mLHTv5+r3c88BtsAKaoR7iEIOIX/8MGDDgQBQiXsl42foPrdn21Wzy4M4X2JCgo5E7b2DLMr54Pw/srvIJiohHY/hF3ihuvoQ/+zVuvY7mZg6c4Gv38if/lf1Hfge2QItf+Rz3S8BRP+JG/Oro/D/rwIHEyqX12zY/29/LmMvhg3TXs7KbqzcwMcGLp6p9qiKiACzyBvG+q/jkL7JlNakMP3uO//4Vjp9j+PRGuBXaoRESfvTNB2Av3OYDcL5FyQImJ+3pyaNrl492d1LfzGM7abXoauSajfQ08fC+Kp6riIACsMgbwYffxid+nv52hsf4yv38xb9yxx38/iepS9z53I7V0AJ1fgDOl3rzZV+r8B9epfTAoP3ok87efTs3r6annRXL2HGU5CjdTWzp5y3LueeFKp+1yAVOAVikyn797Xz0FnpaefEE/3g3D+/lz3+fu95PawtPPXf7czt6/eJvvuwbKSz1BkMvfpOw1y3LnBg48LV/n1izlLV9rOzmpSH2H2JVByvbuWwFR4cYnqrmuYtcyBSARarpYzfzsXfR1sTuY3ziMzx7iE/9AddfhW3z7AuX/Ke/ug7aCpt+I/Oqnb1/3oyU+ejrQA5y8A3g/sdoq2fjCtb3Y2L85AnW97Cmk43dfPv56l4AkQuXArBI1fzG9Xz8nTTU8ex+PvkvnBrnU7/LO67Fsvn+j/nd/7QJLoamwNCjSKkW30AA9qKvA1lIw/7L3zIyMMSjOxkb49INLO/ERPjBE1zUy/J2Luvn++oXLVINCsAi1fGb1/MfbqYxxu7j/NI/MJnmj36F224gGuWJHXzyU8BK2A5NgaFHZWqe8aOvGxgWPPOtL+x7y7aTd/8QYO9RcNi2ku5WcvCDJ9jez8p2rljG3Xuqcg1ELmhaDUmkCq5Zya/fQHOMvSf54N8D/O4H+eVbaIjx7z/hY//ZS5UurFvOM4UVzsE/CCQ24DY1tF25lW9/ZvaFf/4Rn/4mbpZrNnH9VXz+UTI5Ll/O3976Wp2piJSiErBIFXzuw/S3cWyEP/kmp6f4w7u46+3UxXliN7/36Xyqk3ALNPqdn8uXfYOVz/n657Gmuh9dsnG0p53VffzkcYDdx8ikuHgVq5eQi/DITi7qY1UbccPTGh8s8jpSABZ5vX35w2zu58wEn32QBw9y22X89ntpqueRnfzXrzIyOZvs0osZHH47NBZ2vyoTfee0/s7Aued2faOnjS0rWdPP9CS7DgK8cJSZGa5cS38be4bITLG2g/5mjpzm+OT8/IrIa0IBWOR19fGruG0bruGbz/K5JwE+/VGWdnHgBL/8N+ej79I+PvNpLr24/76fLg8M/7UKu1+Zwujr+j2fM5CCaXgRnnn4Ga67mGVL2Lycf713ds9dx7lkGas7Wd/Dd59jVQvLW+hp5uAwp1Ov9zURuTCpDVjk9XPdMn7tKupsHtvP3z0I8E8fZVUX4+N85vsFKf+fP6Wvjb6efTANachA1o+vOX+UkRP4I/9SFjKQhkl41jvan3yak6foaee+//f8W3zpIQ4M0BLnA1fz3Rc5l+SyPt6/5XW7GCIXOpWARV4///NO+lt5cZBf/TrAb9zIB6/BMXz+fr762Plkt9/E+9+B4/DTh08+9sxa6Awsf2QFCr7Bmud82derfJ6EA/At74Dj05wY4OpNtDcQt3h6P8CJMdIpLltBfytZl51HWN/Jpk7cLDtOv/7XRuSCoxKwyOvkH29jRSvj03zlmdktH7iSBDy9j3+6vyDlnTfTmmB4gL/9J+B7MABTMANpvzRc9F/aj75Tt9127C//8u7gMR/aydO7qY9w6+VsWza78Xs7+ZeHSc5w5SqOpXjhFAmbX93OJe2v/eUQueCpBCzyevi9K7ljMw588Rm+sBPgz9/D1es4l+Qz93PozPmUv/tL3HYdyST/+BVePAxMwACshkSg+OuVel2/y1XWXwzYK/sOfelLX9q47uDgKfYfPH/kmRmu20xLHZ2N3OdPBL3jBNesoKeJvmbu+jq3rmZFC5s7+cbLr9/FEbkwKQCLvOZuW81HL6e9gZeH+eT9ALdt4aM3EI/w8Et85qGCxH/3OzTV8dAz/Ld/y287C8dglX/DBns7Z/0+zylIwjichM+9552H+5bQlOD7Pzp/5GNniMG16+lvY+As+4dnt586y03r6Wvhoi6+9ixX9tOWoNHmyWFE5LWjKmiR19wHt9DTwOEzfPqR2S2/9TY6Euwf4JPfLEj5m7fRGmd6gh89OucYh7797S98/vOPXH/dETgLYzAOEzAB4zAGp+EoPAp/DscefQR3mvXL5+bks/dzYpiWej5w1fmNzwxy7y5SGd7Sz4zFQ0dojPLBja/uNRCRuVQCFnltffQi3rMBDP/lZzxwEuC3r+WmjWRyfOZBXiwsZf7xXSzt4Ge7+fvvFmzftoUPfWBs5fJdR4699NyOFEzBGFeMMzD9h584cv01xx998mfwOdjppX/8ae68hc4WohGe3llwqJkk166jt5mXT3J0dHbjY8d45xpWd7CslX95iu1dLG9hSws/OP7aXBQRUQlY5LX27rXUR9h3hh8cBdjUwXu2Umfz+EG+VRga77yadT1MTvG1++ce5LLtNEaZPMsjjwzC9+AfWPIpbv5z/uiP7rrjb++64x9uufGnc3bZvZeoy9svn3uo7+7gwd3Ebf7oVjZ0nt/+2cdIpljTQVcbd+9nJs3WJXxCo5JEXjMKwCKvod/exopmDo3wKb/y+dJl9DUxMsk9O+cmvnkbdTbHB3n0pbkvbV5N3OXkcV7K943qAMBi7BwmzcTE3F2++j1SE/S3cc22uS89fYDJaZa2cvOm8xt/eow9AzRG+D+u5Z9f5plBmqNcv5SN9Ys5cRGpSAFY5DX08xtoirLnNDv8yt4PX0q9zZNHuO/A3MTrusim+fEzc7cDF62BND97JLBpPwAu2UnMDGv65+7y9B5OnqIpyvXb5770tefYdZQYvHNzwfa/+ikTM3TW8YdXct9hsg5b2nnb0gWcr4iEpwAs8lr5yytZ2sh4im/und3y+9ewuo3BcX42L/r+7V0sb+f0KJ+dW5fM1jX0tJKa5PScbskGLCJp6gyN8SIZOHYS22XbqiIv3b+HmTSrO7ik5/zGfRM8fIisw7vW8Z0T7DlD1OLOdaxNhD9pEQlLAVjkNfHWJdyykokUX93DjrHZjbesx7jsHeTueaNsL1mByfH4i0UOdc12Eob0FN/+SWCrV+Q1zCSJ5OhtK7LjC/s4c4b1vVy1Ye5L39rN0WESNr94acH2e14il6O7gY9t4/N7OTPF0gY+uDbsWYtIeArAIq+JD66nMcpUmn/wi7+/cyl9TUzO8Pmn5ia+82KaokxNn58fI6guQiTH1Jx1ilKzJeCJcawcfa1FdvznHzI9RRwuX1fk1bt34mR5+1ouXnJ+4yND7BwgAr+whUdHeXyQVJZ3LufKlrAnLiIhKQCLvCau7CFi+OGh81veuZY4vDjEs2fmJn7XRbTVcWCAx48UOVRzAjvHkTkvRWf/68xAhkyJJYzOjlJnc1mxIuzeAQbP0RTjzsJeWl/ZxdAEvQ38nxfz3UOcSdIe49IlRY4gIq+EArDIq+8PLqItxlSGv9szu+Wu9fQ1Mpbkm7vmJl7XxvoucHhsf/Gjre+FLG6ucKu3XoJhxoEs9dEiOwJ7DpNJs767yEvPDvPkUXC4blXB9ifO8MIgxuX2dTw5zt4REjbvXVn6bEVkURSARV59d6wClx8FyqzXLKUpyoERfnRsbuL3bqezjhNnebpY8RdY0oib5bk5Y5P6Z6eF3ncEckRN8X0f3oWbpSHKu+cNRgK+u4vxJF31/MFVBdt/7xGGJ+lr5JPb+clJRmZYVs+faG4skVeVArDIq+zPttIUZTjJPYdnt1zfzY0ryOb4XrE+Vlt7icHxMzx7qsirV6+lu5n0DCfPFr4wABGw6WomOU1bXfHMPHmY06M0RFlarJfW82d5cYgobO6a+9ILQ1guNyzjx2d59jRRi1uXsSpW7sRFZEEUgEVeZTf2EYMfHeXp8dktt6whbnFqnG8fmpt4bTNrO5hO8WCJ1YfW9FBnk07xwx2FL6yfXRUplSEBkRIlYODQKawcN20u/uo3d5HLsnUJVxW28v6v46Sy9NXzoVX83k7OztAR58rO4gcRkUVQABZ5NX18DX31jKT4sT+L8oYGrlsKhgePFkm/qoPWOJkMX91d/IDL24nB0Nl5LwwBYDF2tlwVNPDQXpws3Y3FX/3xcUamaYpyw6qC7fecZP8IbXHuWAPw4CAYfn19yXcRkYVSABZ5Nf1cP5bhudPsnJ7d8oub6KpjMsWP5hV/gQ9sJWHx0Lx5Oc4zWDlGxuZtnwQLbM4ayNFUohMW8K3nGZ+iu5H3bi2eYOcgxuWGeUsnPXaKmSwb23hvP/efZDpDXx2fXF06qyKyEArAIq+a/30F65pJZvlBoKfVtf1Y8ORJds+brhnY1kMmw4nRIi95tvaDi+XOe8HrhAU7nieTpnT8BThxhrjFTSXKrz89TDZLR4LbC2PwQyc5Okadzc3LeHSSfWMkbK7WeCSRV4kCsMir5oolxG2OTPCAH1Df2kFvPdNpnhookv6WFTRHGRznyXldo/OaY7gZXp7fP2tgtg14cy84JCwumteRKm94DBxaS8woec8xBifoiM8tBO+c5LlhXIdtHQDfPsZkmnVNvKtYfy4RWSgFYJFXzUWtpLL8eyCavmcNCYsj5/ja0SLp79hI3GIixY6Rksdc0UYuy/Pzd18LLsCRNE4OHPbOm98j76VBnCxLS89m9eJZjFswJZbnh8cYS9MR5xeW8f0RTkwTt9habNYtEVkoBWCRV8dvr6Q9xlSW3YH65EuXgOHhE8V36W8El5eGi78K/PxFOBkyGX4yv5E4NjsZVmOaqMFyubKv5HGOj5BO054o0tDr+fqLTKXpbeD9hQmeneDoOPU2b+sF2DeObbihdFFbRMJTABZ5dbyjl7jN7lFemJndcn0HbXEmUjw5VCT9W9roa2Q6zZ7SAfi7e4nbpIpOMzk12wb8BGTSuDmWli6Y/uAA40nq7JKF4GfOcWCUGFzZM/el7x8hk+PSTt7SyNePk8yyqoFfUEuwyCumACzyKviVXlY2MpPlO0fPb/zQOpojHB7nuWLdr65eSnOUsSTfPFjysL91FVFDsmgA9sYdWZBgOkWdTV3ZjljjSaIWl5de3PfgKMZlzbwI/c0BjkzQEuWyTl5IcWCSOos7tUiwyCumACzyKriikzqbA+M8MH5+Y089jsMTg8V3ubIXy/DS6XKHXdlGJsvJUn2kDViQJA4mx5nJEskAeHkYHPqbSia45zDpLGubuXReDH7mNJbhfSsAfjJE2mFFPduLrUAsIuEpAIu8CtY24ro8EegG9ZEVrG0ileOREgG4PYFx2F+6+xXQ1UC9zeliBWi8im4HmklnwKGvdHAFjo7g5ljWXDLBs+OcTZKwuWVeO/HDQ2QduhNc28xzY5xL0xjh59QSLPLKKACLvFK/s5Rl9YymeDFQ/F3TRNTi6ATPTxXZ5eYeeuvJOhwoPQIYcByMQy5b7LV2sAAYxbhEDKb0ZFjAZ19gOkVrnDtKz6QxMEXCZvO8UUYPj3F8isYoVyzhhTS7xonCuhJTa4lISArAIq/U2iYSFmdS3BeYr+q6bqKGJ4p1vwI2tNEa48wU95UoH3uWN+Pm2Fe0mtp7LxtamZ7BuDRVWikhncU2rCrdV2vHadJZljUUeWnPOSJwXRfAt0+RynFNOzcpBou8AgrAIq/UJa1kHf5XoDPzFU20xsjkeLxED+e39mEZ9pWtfwYaIqQynC7auDsCLlgww3Qa26Wx/GxYcPIcEcPGjpIJnhxiIk1rjBvb5710mlSOtY38/BIemmY4RcKir8TMHiIShgKwyCtybQMJi7EMjwfWS7iph8YIZ1I8UaJj1NIGHIeXygbgSzqIGCKGn5QYRowFLjhEwDKsnxc15zg7DS4xq2SCpyc5PE5zhF9aO/elu89yYJx6m6s7AV6cxIYPlB55LCIVKQCLvCJ3LaUjxqFJdqTPb7y0Awy7SrTvXtVKY4RMjudLz10FJGyaoqQy3FgqznmNvv2ks1gu55IVsjqaBJclJVYO9uwZJevQUyzNYBIMG1sAvj9IKse6ej6oWbFEFksBWOQV2d6CcXkqEGvf0cayerIOT5QYYrS+lcYIZ5M8ca7ckS/vI2KYTPFQqXZiCywYIm7hOmRzFbJ68BzZHG1lhw+dTuK6dBVL84kXyTr013FpHQ/PcCRJnc0VCsAii6UALLJ4t7XQGeVchmOB0mdfHe0xzs7wjRINwBtasQ1T6eKv5kUsbMNQ0TFIwFawwIDBuODSWmlg7ssjZHM0R3l36arjfz3Jy2M0RfjV/iKvHpui0eaqDoBjSaZzbFA/LJHFUgAWWbz1DcQt9k3yw0CYvKEL23BquvRezViG4Uo1xps7sE3pcu0e/48hTk1gu7RXCsBPjzKVoj7ChrKtxefSGMMN3UVeOjw1OwsH8NAIGFbW8YHSazyISBkKwCKLd0sXruF0YVm2rw7X5fnSA3zbYzhO8Qmig6KAy5lSgXwZADb0kc3hukRCZHh4iijUl016aIKERV+xZuAjU9RZbGoGuHuS0QwJi4tUCBZZFAVgkUW6NkFPnNMpdgfm33hPO311ZBz+2/Hie727k+46sg4TmQrHb41hXEZKFZSd2fWAGaTexnapNAwYYCyFcVhdej4s4PQM6RwdxcrT/98pzqVZXc/trQAHprENq+tDvLGIzKMALLJIH+jFgpEM/xboS7W8njqLE6Xrn9e1EIeJNN8tvQiSJ2FhHKZLxWnLnwkLRmYwLumiE2YVmkiBS2PZEvDnBziXIga3F6upnswSt1jVAPDMODnDunquCBP8RaSQArDIIrVEqbcYLlyq6LI2ojBVOhaubCRiMVSpARiIWRiXZKkA3DTbA4s1NEUxLnbpAb55Q0lcl6XF5roKmszSGJntbDXHoWlsuKEDYPcUYxlaI/QqAIssnAKwyCJd1MhUjhcLp3puj5JzOFhs/mfP8npwOVE6Qd6SBJksGafEy1m/EBxnIo1lsMvOBe15bphcjsYI28vG4KEZLJdtxYYYfX+AjMPKeoCnszw5RgyuUj8skYVTABZZjOsTRC2M4bHA/M+XJuhNcCbNl0vNXQVdcXIOz50tmcCzuR5cjEuyVGHa6/llA2QdTNkprvIeHGUmS8LmirJrGR2cxHVZUqxcu3+GkTQNNr/cDnBkBsflxnnrN4hIRQrAIotxZzdNNgemeSEQIK9rpzWKZThQYuzQO1rpiOG47Cs1utd3UTt1Nqks3xsokcLb7sIYlgGXTKWJODyD09imwp1/cJJciabiIw7TOWJwdRvAiRQpl3iI2C8icygAiyzGsgQY9hd2ttrWTMxwoHT18ooGojCTxa1UXRyzsAyDpTtzsR4M2FBHOodxaQx3NzsuEcPKsmOHDk4znSUBv1asoPz8GMalOw7w7zOkHeotfq1sz2oRmU8BWGTBLovSEyPj8FJhrO1PkHXZMVZiN0jmsAwnp3khVTKNpzFKxNBUprvyNFjgQJST0zg5llSaiMNzeBwbWsp2m9qV4XSKCDQXy8DzE2RcliW4NgGwe4qExds0J6XIAikAiyxYe4S2CNM5Dgbi6KUxuuOky/bAurYD21Qu/gLtcTDMlBlZZMACGyaJglWpVjkvAq5LV6WVBI9M4xouLRZWvzvBSJqWCKvqAA7NkHHpUUdokQVSABZZsGtbiMH+JM8FAuTyOuptxrP8uHS9cWccDC+Pl0yQl3PAZbzMfNFx8DpID+MabIPrhsr8gXGMQ3ulxYMPTpXshwWMZIhaXN4C8FKSlENvjHeGK4KLiEcBWGTBLm4i6TBR2OlpUyM2TJftCdUZw3E5WqZl19dbh2WYLDNb1jH/9l2Na3BcrHABeDKL69JaqcD6kxEyDr1xbm0q8upQGlz64wD3pJhxiFusrlSqFpEgBWCRBfMCzwOFsz2/pZmo4XTpxt3r6lkSI+vwSOlpovNWNIBbehosYA0ANjSScsANezP/2xDTWUylaP1yjukcdTaris00uXsSx9DnF3kH0sThUk0KLbIQCsAiC3NVlEaLZI7vzxRsb4uQdTk6U2I36IwThYxTcpBSkAWWN9FVKYcBcGGEzhiWCTsMyduryeZjSyukOpsmymxv5zmemWQyS5PN+xoBdk2TdVmqZmCRhVAAFlmYd7QShZHC7lFXxFjiDfAt3QOrN45Vtogc1BLBccrNKU2nH57TnEpiDJVadc87lSTtlJ5jK59sBuMWX+zo+RznMtRZs5UBf32OjEt/nLUaECwSmgKwyMJsrCMKZworhzfVE4dkjm9OltxxbQNRwobJuIXjMF6mCnrEX4yhnoRNMkvc4pJwrbCZHI125ZmzHhnFcekokeOJHDG42A/PpzLUWfySBiOJhKYALLIwHRGyLrsKS7q9cWIWk2UrgVcvV0RnAAAgAElEQVTWYRl2V5oDy5OwyLpzV3oosNH/I8O9Y7OhPWT5M+MQNVxSKVh+fYKcUzIAvzSNY1jiv3omQ9SwUh2hRUJTABZZgPU2XTEcw6nCsulVzUQMeystseC6DJZuJM67KOYvrlCmDXjQ/6OFVRFcFxO6H9ZMDteEWj847ZJ1+HCxtRYeHMNx6YvzlgjA3iRZl3UKwCKhKQCLLMD2OHUWaYfvFJZNu2JkXYbLDNuF9ggph5Mh2oD3ZmmK4DiMljlgz/nwfNTBNthmdpGiig5P4YaL1hmHeqt4tfmMIeWSsGiwAAay5AyNNleWXWxYRPIUgEUWYEOCqOFwYRC9JUGTTcbhQOlVft9RT2eMKJwIEYA3x8DFMjxbJgB7pW0DhwCyDpYhHW4ocNrFuMXHF81xNkMUNhRbu/CZHENp6uHKJoCvpphxSFhsVl9okXAUgEUWYEMdURgurH+edInDmQz7SsfLgQwxmMzyZJl+Vb5VddiGqTLzUAIt/mIMy9geI2JwXZxwAfjgFFmnbP2279gMwKoSfbvGskQD82+cyBCDFj1URMLRvSKyAEtjOIZdhSXdbfVY4Br2lB7Yc2MrEROq/hlY30jGYab8uN5JsGdno9yVxXWxDL3hekGfyeC4tIaoK96fxHXPT7gxx84pcg7r/DedcYnCpcWKyyIynwKwyAK02GRdThSWTbuiRM3cmSnnmHGwoDFcN+WcS71Nonxir7RtQx2bojguGMZCFK+Bp9O4LrZhc6VBUUdSpEuH6nEHx9DsP0X2p3AMzRoKLBKOArBIWB+I02CRdDhbWNLdXI89r156jlV1WGbu9B2lNNhYbqXZIqPgJThIysUyUHrthPnGM0RhXV2FZEmvc7Xh3cXK1jtmyLq02VwXBRjKkHOJKwCLhKMALBLW5jpiFiM5nisMjT1Rsi6HylYvL40TMZVnYPbELTAcLz9gaRi8cUo9NNqzg4AfC7HOksc1WIauSgF7zGUyh2WKt+w+6zCSJe3SHwV4KYtraLW5Qs8VkRB0o4iEtSFBxBSsAexpjeC4Faqgu6JkXY6GawNOWKScSuN01/mzag2zN4fjYlwyIfpVeYbTRGBFpRLwLofxLAlYUaJ1OedSb7EiBvCYy7RDm81GjUQSCUEBWCSsOnANU4X1zzdEiRmSDp8vu8hgFFyXl0IsRAi0R2mwKhWXbXDBhl6AmRyW4dpiSwcWNZUjlSvZuypo2sE2dJaIqaezxGC1f5zRHDFo03NFJATdKCJhdUVxDQOF7bjdESKQrLSwQVuE0SzJhYzTHShfXM43OecAHLfy4goFeztgCFNSPZTCMfSU6K51IIVj6PVfPZ4hClepI7RICArAIqFcZpOAjMuhwgC8KUGUucXiOe5IkHZosrg35DCkOjCVJnYe96d+jgC4hgZ7dm2iMI4mSXC+A3MZZ7OYwJzPczw/Q85lRYytNsChDDlDj6qgRUJQABYJ5Zo4jRaTDvcVtvW2RbBgsGwX6DGHzgjHy05UGTTjQNmlhQHOgQEHGgHGsxg37ExYQGMEy9AWIlI+npzt2/yWYk+LH2RJucRgjwuwN0vOxQ3dFC1yIVMAFgllWZTYvGWAgSURjMt42R5YvVEizM6ZHEa9hTG45aPpSjBgzc5J2WhjhVtfwXMoiWsqzfUBQLfNtENHhFKFfOMShQ/HAbpsHMMSPVdEQtCNIhJKq03EcGZexOqLMpLlTNkBvtvrsUMPAgYiYBlayhdPMxDxkgKkHEzFWuuACQfHpT3E6sQ/yBIxRE3JLtbnHCxmpw3Z75B1icMterSIVKK7RCSUegtgfF4xsDNCs81I2aJknY3FAoqnTTYzToVSNU3ggoFhgDhY4TpVeVrs2ZJrGGmXKLy1xEikM1lihivrAZ5ySRpiFtdqXUKRShSARUJZFSXrMlhYir0zRhSyLi+WLd0uj4HLydBtwEkH23C6/LySSbDOL6iQdcEQbipogBMZUg5WuDbjwSxZd/YnyHzPJskZ2v1nyeE0EUPoxmiRC5cCsEgoMy45l92FgbbdxhimHJ4sH3AMrmEy9DChJREsQ135CuVTfgBeB97MVi5jIdp0PUlDBGy4PURRNe1Qb5XsCN1i+3NyAZBywbAhZOFa5AKmACxS2cUW3VFyMFrYDjrjkqg0BgnIuRjDqXArJWyPYBtcl6HybcZL/OUIRwFiXqet0N2Px11wabQ4F6KsmnGJlJ6LY9Qh589GCezKkjHENCO0SCUKwCKVrbOJGpIOzxZGuM1xLFN5+YHuCK5LKly1bLsFJkSX5noADCwHcA22IRu65veQYdrBhOuxfDaH49Ja4jRfzpF2STqzU0BPGxzDUg0FFqlEAVikslURooaZeeHNmwbraKXG3YRFyg1bOj3rgovjUl/+7vSK3RaMA6SdBY++zbhYhqkQMXskB4bOEvk57TKaoytClwXw5RwZl4wagUUqUQAWqazPJgLH5tUhN1rkKtYVw6RDHPaHG4a0IUoEzmW5u3xc96Z9dmdnwmqPgMEsJAanHaKGK+orp5z0onuJg+80RCxiFkkDsAJOZc/3yRKRUnSXiFTWaGEMo/PaetttXFMhsr4zQl+E6YodtXxtEeJhZu0YBQus2aKwMVilm2mLMxjDRIiuYUcdzubojLCuRK6SDhHDlgjAcX988kdUCy1SlgKwSGWOwRjmd6Jqssi5Faa/WB6ZnUQ6pKyLMeQqpj/rjwNOA6RdLMNo6F7QwLSD5bI0RHflUUOLRc7lYIlC8FAODCv8iOsYbEO7+mGJlKUALFLZyigWnJ5XWKwD43KgbCGy2cIy2KErhydcjEusYvq1/jRYDkAyRyRcNM2bzBExoXZ5yiVnyq3ckHWJmPNPk2lDjNkmYREpRQFYpIJtFnGLrGGqMCh+MEIUsvBk2WCZ9UrPoUvAMbANuYoB+JBfAq4HGMnhmIW1vE66OKHn5xpziBg+ViL1sMO0e/7dB3M4Lkv1dBEpS7eISAW7DR1gXJ4pDKLeeJuKc190R7CKTSJdSsJbiaFiAO4HwJq9ibMulik5WVVRaQfLhB0clXLKjYyaNsShxX/3Yy4uWpJBpAJNVyNSwTUWaRcsb92/8zos6sCqFCnrYCb0IGBgwgVDqmLfqAH/jyOz/zVhys0BMUjmiIbbxTFYhmyJxGnDtEOfH3GPuVguIbpXi1zQ9BtVpILlFq02LjxVGH7ikHGLdI2eoztCs7WAKTJWRmbLphX0AGDBWgBjwBBZSAA+mKGOsOs3eHNNd5d4YAy71HN+NeI4ZIwCsEgFCsAiFfRbWC6n59Uhr7CxKblKbl7GhXDzXXjWxWig8thiov4MzF7Fr+s3CYeWNhhDV7hnwJhLlJLNujGYMrTabLcAjsJZh7jhGvXDEilNAVikglawYXr+NFg2jXC8UuOuN4AnzHBbT8RgGZoqhq5Tfrg9CF5jUugGXc9soTlc4lGHqdILIn3VBRfXr3Z+1tBiEZ8tnItIcQrAIhVEYMZlfF7hcsLFcUlUKnQ2WeAyETo0pt1wvab7/dWQOgEaLCLmfDeoMJosjBt2AsvjDvVlq7izELFY6mcgY7AtLckgUo4CsEgFDjTA/FX71kQABioFMG8ppPB3mjeLReVJJRvAPh+AZ1xYyIqHwC6HlEtduMRZcA2dpQPqtEO9mW2YBgZdKN1pS0RQABapaL0NzE50XMCEanZNWFgsYIF6b3RT5dbcKX85QgNgm9nlBcObMkTg3EImzyojaxhxyU9fHfMq0l+dY4u8OSkAi1QQA0OReSgT4MLuSqG1zuC6PBO6bNoEEROiYDrtB+kZ8O/k1EJKnA7kXJZGeHuIsG0bjItT+vhpQ7fFCv9QY4a04SI9YERK0zhgkQqawLjsmRdB4+C4Xvgr6RJo8hb+C/12k26FUOen86PuSYARF7PAYUgGbIMFR0Ps5fW27ih9GmkD5vxQ6TNQb6kELFKOfqCKVOCVgOf3Fm4GTIU24Be8/5gF/NRttsBgKtZZr/GHIQHQZWE5C5uIY6fBBsfwlhCPgcfc2QmfbysRg70ZN/OLOKUNlus1T4tIcQrAIhV4sz0m523PGSzYVrZse7E3VVbpSRzn67aIuAxXrLJO+g3LvcDsAKRE6HcBNoHrgglVAl5tYQwRGCyROGWIGhr9V5MG1xS5aCKSpwAsUs5FkDA4LmPzXnINx50ibcNBO/1xvUOhy6at4LgMVgzAFsF67QTgLmDJB2CVBS6WS1uIxMcMjsE17CiR4CWD4xXfARgHoH0B2RG54CgAi5SzBnKGjOFwYQS9CNoseiyOVYqsWcM5w/LQjcDe8Sr3Z57wu0APnd8WWUgv6PsMjkvUMBnix8Eyi6whaviFEm/RZrACo4qnIas+JiJlKQCLlHOY2cGvvYXbN1lEXXIuuysdwRjiLidCl4CnXTCMVUzf5bcBbwSv/7MbouW40JiLY7gixGPgh4aoIQLfKJGxl8AJdMLKMVtlLSKlKACLlLPOa2k1nCncfsAU75k1x8ctGsGBnaHfscUCQ0/FsuxRv3XaBWg04C6sFzSAIW7C9tBOGyIuHy5VAgYnsCLhmCFiwi61JHJhUgAWKccALiNmbkl3s4UJsRKDS7lV/IqyDZgiM1/OtdZ/+wiEyElRMbAhEy57WYNh7pqMeVNAoFW63u8ltmFRGRO5ECgAi5SThjhFVtZLG6KGYZeLyu5+DswCC4IZg3Gpr7jLjF/8PQcwanDdBU/9aIMVeojUJACp0ocy5vwqit+CONiweWE5ErmAqJOESDlbLSwYmRfYeiFiaIS9ZXdvAWuBxdOUi22FCKVRv79WDqAd7IXfz1kXY1Fx5UNPm8Gi5I8Jy8xO65E3A03QusAsiVw4VAIWKSdpzi+5G9RmQbHAPMeAwWIBg4CBBgtjKoxuAjgOgA1LAJLgmArTcs23oEWEU4ApmfgMuLM/BmY5YAq3iEiQArBIOd5SCvOjjlfWbKi0ews47vmK2TCmDRSbenoub+EhAy/PvlHEKtlAW0rKxXbD1pCfMsRgXYlOWGnIghVo9I2x4B8fIhcUBWCRckyJeSgTBRNBluSAHWLJwqAGwC3S6jzXsD8TVh/AlAEHe4FtwBFrtid1GN5vjlKrRDwELjTDfn+LMdjGa6EWkSIUgEXKSRksU2Qx4FGwYLrS7jbYpYNWUWkTrua2N9AM7MfiUj2kymTP5fwaguV59dtlRl55vzY+6P/fiMEsvFZc5MKhTljzvAdWwX74cbVzIm8AbRQv5y4xEGIE7Q0WURhdyDt6U0yGmsTRi7qD4Mf4hQ5GSoAFS4pWss/TbMCUPRd3doVij/eHFkQSboF18Cgk4UC1M/NGogAMwErYApfCMmiGen8qW7ngtYHF3Fk4gAkw3oJIZQ0ZHLOwNRKiYJXu63Rekx+rN8LLuOF+EMzhzRwZcp5qL7qXmTja69I1CMB2GIUeaFlgluRNqAfeDm+FaTgJL8IeeKnauXoDuLAD8E2wHjZANzRAAqIQBVu/22VWK0RgYt52L6ZWrLz1uiBV7lEVMG0whKi5TfqdsyfA/3NBkR5IG9xifbyL8h4WZWrdo96awQDsggjE/HgsF7oZ6IF26IWtcDskYQ88Qcn1PS4AF14A3gyrYDushg5ogBhEIQIRf26/Z/x5feSC55Vx5/d2rodIiD4UXhl1QZ2Tvfrgyr8ApyDnf2/9NuOFNrha5w9QmVdOLtOenYV6v17cy5JXPr9/gbmSN5smZidR8wa9OeBADpbDjXAaTsAOOA77qp3V19cFE4C3wqVwGXRDCyQg5j9BbT/uegUB75mkddQE8At880vAo+CECHiJ87M1h9UEdphm4yb/0CfBD8AVq8Tn8PIf8veBHfjfUglMYQA2foFYLnTj/sgB/G+J8RstumEjXA/TcBQehB9VN6+vnzd7AH4fXAxboRnqIO4HXavwnyf/hw0ObPG/NGlohGn/2zMAS6EJJsAFGwahDwahH4Dm2dkBaYEZOAVbYcavjhz1q+pOQV9g3MYaGIVGyMJpWA8HYCl0Qgr2w1L/++rlMwIZODPbBMgmSEMOZuAs9MMA9EDWb7WLQ87v89riN3JPQgc4kIVWmAQLBgDoBwMNfpwxftcgL/LEoQFsOAD9/hVrhX3QD2kYgaX+M7gBsnAKVvnvbmASBqEH2uBl6IIE1PkL33fBEHRTMAVzJ7RDElJ+hk8CsMz/YymMwSrYC2sgAtMw4F+NPsiCgTaYgg7IwBh0ehMiw2FYB2ehCbIQgQHcftwB/z7pAQOtMI09AF7d8ibIwTS4MAyb/catZVBH8iA5L9l6f3LkfJ3sGsjCCeiFHJyd3WvyJAbMMnBgClrB+Bf8FDT6jc8Z/5twCbxAehO5fST7oQEMHIJev1XlEGyAM9ANL/u16jFoYeQAUdjuXboULIF6mIEJSEMPOLAflkEC5xBAegOMQDdMguOfTidb2rAPYXvV4EvhFE7w9loJx2AtHII2qIMZsGAF7PI/oFVgQQJcGPGXe3L9bjtrIQL7YQOMQRa6YQjGoQNG/VvGO9RFsBcu8X9iTPlVoBlIgeP/6PCuUiMcgrfAUcjBJKyFIWiGFIzBJtjnf2+9K3/Wz9IhAFrBgl5IQg6ykIUef3C0d3e0QBKOQze0wVmIwmnY4CfzTjYLo9APOTjtf5G8G78PUnAOGiEOo9AMM9AAdbAbeiELI7ARJiAFrXAUlvglDa+uw7v4jt+LfgA6wEAvDPnPLu/JcwoS/r2QhB4YD9zObTAG3RD1c7sJYrAb+iDnP2mz/lM3H21M4I/gPxeWwVXwH/1i8eNv8mD8Jg3AvwE3wQqwoX5eoKVSq5d3773X73Zi+V8OO7Cv699XUT84mUB0tP00+d29xBE/jfGLLy5EAy/Z/nfRm8DQgqj/Dc5vpzCZ8d8l/7dbeL75vfJfdNv/v3Ygb1YgWT5vEf+wlv+j1Q4cNn8x8we5o/AgBE4quFc+z3Ygt/kzsgKn9i4/Tf4gwaO5xa6ABbf7W9zA55Wv53ADh7L8M7UDbxq4dN/LcPV9fOc2/2eTf0bDk6R+ykPvDNTJ5q/Me85veTjHzY/xxWshFvg+BD8Xq/BzNHwtR9MRnlxd+EsxeO6W/9F43g3v4CGXGzr4ztX+FzJ/asF3nPPFsPhalsvv5dffX/j1CF6EQPb+1aVjF//XpYHEgQ/uRYt7x7l5J9+5Yfbw35nh9hM8sgpigY/MDhzc+2RvKfxSBc/XKvysgci8+yj40Tv+ZbHgtvOnWfAdY941z2+xC29PAsne4x9hzvvm5b+WbuDiz3nO5DMZDTxAgs8KO5Bszncy/3b5Kz/ntn134aMpuKNdmG38748buLmCe+WvT/4BGLx6wczkL3jwa3Nr4NzzH7RVqcEm2BOwDfrhYrgZfhcGYQ/8Xdnda9ObdBxwAgYD0XGh+oB5c/RF5z0KCVS6mUC4Csaq4HMzErhFCdw8pvCBEgyTduB+tgpDSL7OnMCOwRhZ9L6NBO7YOSHHCtwzbuD3cn5fAgmC//K35ZxbzgmcSzCTBM6IwB8Erhjz4uKclMErnL8yBBLMiVh24LmW//jyz5dIIHGwytjCsQEeiPj58a/8YB3A8Zj/MdnzzsICQ8RlqJcoxX5nBC9s4CvRm6V3CItAE0n+AgYDaj6TLoDt0jDCj/MNK3PkdyyMInFIdvLOrP9h5S+UHTiIv681J6JYgTzYAE0pzvTy1uzseXmzgjwQD3w6c26QXOAbOOc054TGfJYc/z6KBO6C/Dcw+ClbgQ8l+IM1fynsed/D4HfJDRw2+OWPztvXBPJsAr+Vg1cveDpRP3v5l/Jf3fy3omjAC94IbiBZ/sgE0tjFEgTvUAJX0irMj+vnPPjMyeckeNEihZ/anLMI3ryL4GUgDWloghsq71Fz3qQl4H/0//glWAvroa+w3XdOCPEEn93T8PfQD12wE5b63UsOQx/UQRo6YQiWwxgcBvyxF+egGSLQAbthHZyGCVgGGTgNfZCEKKyESZiCQVgLURgHy6+X7oBuf26FOBgYn63/nK0lBi6GQYhBDhphxj/UpJ+yHYYhCkOwDc7AkF9N1+V//pZf2+aFk7MQgTj0wA5ohDaIQBQmYAkM+3WbGYj4lZDevBCdsBPWwzg0+U+rl6ET6mEK6sCFHpgGGw4Cfh3yapiBIVgKHZD2KznbYBo6IOXXW3rVenX+ZffquxJwEjbCpN/0amAMRqELItAEg36j7iYYhlao9+u6z8FpaPer4JYBkOKXz9IE/zDOf6yDBEzMXvy7oAluOsMDCchBK4z7w5W801kFR3kHbIWPwH3b4TQYiEIdTEKLv+NJ/0fAMMDP9bFkkBuO8tN888R6yMERWOq3OyyDJPwGAHfDNO8/STP8z//Bb/b4V7UXYnDC/6qsgqOwCpJQP3uJPrSPJfC+x7l/FYz4dYw29MCuwLm0wRgfhlXw1z/mT/v9xoIInIKLYYytx7gJ6qF+CZyFHu4YZgP8GfzVZj9g7Id1fpPHydmaaoAV/tTW3kczCU3QAGegwb8IqyHinxrQ5T+XvW9yh9+Qcca/Si2Q8nuCRfw36gs0GG2cncWTbmj2q6MtaIVT0AUtMAVj0AJnYcL/emwEYAIGYQm0QhqaA+NqOmEEgHZIQASawYIYjMApWA8ZZhfBGIC1kIJTsA6AKRjyzwVohBboh+f8hoxT/sHrIAUtcMyvJwfWwCSc8RYMgVZ4EVb5l2scDKRhAtb4NcMODEMc2iAK+2AZGDjln0sPuNDiPzEmoBUa4Tn/m9YNJ6ELXvbf1xvM2QmT8POBuE7g54gp/OeCAxn/ETcID8MDvIkt+sdJDfoQXAQb/BFHZTphPQEH4L9XOb/yRvBZeDs8CP+hcPsfw12wE/63srv/33AXPAW/Fvodvw7b4V74w/LpeuBXwYYvwhD/A34OHoCPh34j4IuwBu6Dvw6R+AtwLXwf/qDYqxvg89AI/wW+BcC34SL4EnxqIVmSN6Hfgi1wTaBMH+yE5f0gy8EonIJ9sBuernKWXx9v0hJwUV8FYA1sgMthJbRBfbFhSE6IOQblwuCNrJ0/MWQO6kLM1+K1fi5iMorKX8CIXwHYAkOMgLvAOS+BUbi82DQjRUX9ty1qP3iVzfmT9XKnAX1CDBxIzxuGlIUkTMBxeBx+Wu18vu4upADsOQyH4T4AboAroAv6/Y6F3kQcW/26Kbng5bv+zOFt6au0+xTYCwxCKbAtuiv2XfDaJszsNB9Nhf1MQ7ocDGHXAy4/XHi9355zMpA+ucD8yJvTNHTCNOT8kRpTMAovwWm4p9rZq54LLwAH/Qx+BsAmuAyWwlpohHq/CUcueN4wlvlr6oUMWvWFHZbDaLZwLSap1H+w1f91kAOos8DCmEp7FfI6yrSGS+wVykv12zzg9yLI/9roBLfYEGq54GSgHs74U1HuguN++/GF7cIOwHn7AjOwvA+WawI9mbUXrik2VUUHGC/glZUGy8JZYFC0wlRBn/N7LUyDP9vlgua8BFwbF6bcyimBiAUW9aVjvFfnvAoeg3WQ9Tv9yIVuDyTgHHy72jl5g1EAnufuamdA3khSXs/eedunLCwLQ4USZ8YCi3ilZEFp/30rGPNjdQ+cmy1kV17EsJA3WULI0nw0kL35rvR7yns90+v9Bun5s3jKBWfHBT3hcxlv0nHAIq+SLnAt4vOGC8QLh7CVYgIzuIQU8d+3giX+H+cgMC59QTKFY1bL8+rhS6338DS4NtisAWA3uBbGWnChXOTCoQAsUk7KwrKIzIu0M4EZFMqYNDjQv5Dhfq6FsUN0X8qvCzgMfnRcaC9oy8L8/+3dd5Rc933f/ffvd++U7b1g0XshQIIExSKSlmSRkmhLshRTeSwpUR5LsR87dopz5Cdyjp0iP09OnOTYjxw/iR1HsaXYkWXJomRZpAolUuwEiEICRCGIvou22F6m35s/Zi8wg51yZwlyPAef18GRFnPvzP4W51x+9te+P0tLuOblo7TC6oj8SPYpALaCa/Dzk9kiUooCWKSSfDc3saiTmB8irnr8X6K4slYYuYIiSJXMBu0bAsgacEjXuLE/EhSkCqMlqDRcThNcgS1AUN3B6jAGkfIUwCKV5IMwsijY8gPFuWpvjxh8S7aWXDQGLLmqb0kGj2/k6v/UXFjHFpccrSzfUZ4v86vBe2AK4sFk38FgX32tBzSJ3DwUwCKV5EtJLO7Cjvhg6aj2AOWHbRfndwXGYAzdVd/SCR7BWUvkDFj8GhPYu7qULIT8LxOZMt+iB7oNnuHA1Q8HDw7X1iKRm4gCWKSSV33iBfOtV12CTIjNQnMUn9MQQr5YULmO5jVRANyF0lP520OuZ74qafBtiBXXQNDdLzec3mpwzLVfNX4GmiAbYpBA5KalABapZApylrjl1uLXuww+9FVL1tZ8Ef5aAtgajMVUfctUwWFc+XMrDM019oDb8rUsw73LCcr4l5RvyNUsz6/RzsILtbVI5CaiABapZKDgTMhCcUL1a31CRGmpd1VfGOUG67XmAFoMjsGr8XvNG4wNWynTtVhbdk73aR+/YFVaxmBr7PqL3GwUwCKVHA7GhAeKXx8DDGnDxopvT0LCDzGeXMAJ2QO++kvBPECq4HD0GliMDTtKnAI/2GW02MMGF8aDn3ThbFkFsEh5CmCRSs7BmI9TcK5u3hAkwTVVthhlwLH4lneGjiLf4NsQT2ZrwVH2EDU4YWK7WIchV5CalbmGrGG2zM3tBuuQDBqQ71WfreU3D5GbjQJYpIq4YQ5WFb/4ZZ8oNMP6ipl33Cdj6XFCVLYKTIFjGajamU0G3cw+gCs+XvklyiVtAwzpcAVp74JOi2fKVu3Iz0dfPbXilmDOWETKUQCLVGPoKVWaatbg2oXjkspJBPHWFToa88Wzxqpm17kg3yYBBmzNFT9WG2IOORNqFfQys1C0slxpyZgBc23Hc2u+glgt7RG52SiARaqYBd+wbWLv1fAAACAASURBVFGCxgy5akWbvWA50nDozuC8jzGsqJpdm4LH9xLkDzu3ZZcolxQ1eIaow3dDtO1bPo7FWJ4qfxSSNVwpvBqmnIjITUwBLFKFY3EsbYuyJAnW0loxY07BrI9HibeX02TxDNGqj+bVc/7WAbRYrCVeywP9LQ9TbQ77qk9ajGGmfFTnO/pXr+dbd05j0CLlKYBFqjibw1IiaOcgYhiolqw5Q9Remxytat7HDzMs7Ad/ZhYa45maq15Yw7TPPSF+ObjkM+VXqqm5wmEezgSJu8Pig6cAFilPASxSRYvBK7XDJwO+qb7TNQ3GsCb0ozYLxhCren8UfHBgFPJrp8O8q8AnHDIQt7wYIiZXGXqdSgFvIFZQVKvDguFyDc0RuekogEWqOOhhLX2LnpVpH2yJo4KvE7fkbA3H4rZbTJj7M2Aht7BDudXWXIij29LiMBGuk2ryZwuWv3m5g1dQQdMFz3BJPWCR8hTAIlXMU3oSN2twDcuqrZZyLY6pYYdut4MNc6xCc1Ai5BLkFyebsAcL5qXBDzPZDORXWZffBAw4FizjwV/jBkPYk4ZFbk4KYJEqspQe3b3ikQ2xuupcFsfQE/pRS0DOEKu6CnoaDLjQD+AZsERqeaBTPoQete6xVYpWxi1pc63jHjH4hpPqAYuUpwAWqS5XKkGb7UJ55MrmfDD0h6l2AcCFHIkw5zcMgQGP/ERrzEIt/WzANZjQO4UiFmvLpvVPWVyDz8I2pO1gDEm4UkNzRG46CmCRKp7yAOKWncVZdSDLLNXPG0iDZ5gL3Rd8PUvc4lR9NDNFdTd8oMZV0ClIG86GG7bOT4FfKfMNHLAWz/C0D3AImi0Jn/3qAYuUpwAWqaIZ0oaoXTiB96q0jwvLq3VtUzAPA6H3IeV8fOiuOgR9NHh8l0N+CNrg1NID3hylKXT56A6XNJwvs2u4yxCBqeDqhw2uvVYXWkRKUgCLVDEHaZ+oYVvx4zLuh6p9cTFHzJIO3RdcFcGzmKqP5joA/IUzeDPghSsqeVWHC6G75vlTijvLtKrf4huOB53ptRZLDT+yyM1JASxSxSGfBPjm+locp30yMOfz/oqPUQKMoT90D3jaw5pKJS+C+wCwMASQA99woZZl0CkfY3DD9VONAcOTmdJXuyyOJRZ8VLcFqz1IIlUogEWqG/dLnFB00GfGp8elq+JjNJwl4dcwO9vqQIi1XUSD6ec5gLUxEpCsJfN6IlhLJNyBCe3Owsbj0h/lYg1jwRD0cge3xsaI3IQUwCLVnctiKRG0WcgaLlSsp9zl4Fo6Q6+Czvc1q/dL24J9wCfzbyPuVNqnu1irA4ZD4YatI5Y0NJdp1jKXrGHKL/rrvpp2JYvcfBTAItXljzpYvMRpDlosOyoOL7+cwTO0OOwK97RdyuIZ4lUD+xgY8KELoM2tYVNvnmfImlD91I+6C2uuHy/Tkc+PYx8OEjdmSPpVBgZERI+ISHUTHlnoXRSKF3Ng6as4iptfoBS+TNWMRxouZfnpaMX7BsCAWTgQMRWygnSBThdjaArxFtfQ4nC+/DB6kyVnFv5rstUQs1jDc+oBi1SkABap7lwO39C+6HHxwDUsr9gDvuQx7mEN8XDf6/UcGVgeo69yJ7gjv/IKJgGmPTAkahmCTvpgmA1xHqEDri0b1XdZHEsWDnsAR3zilnGfFfqvi0hFekREqpv0SfslEuh8Fs9Uqdt82KfZkoWhcMudhn0iFtdwKFXxvjFwwF94iNtcfENty54MWcN0iADeEceWP4lh3KfdYdyjGYBHHCKWJsvrIY8aFrlZKYBFqtuXwzO0LUrQ4Sw5WBer8vbLOeIureECmGBt14bKH+sGlbAMgGsxYepnFVgWwwu3c6k/gjUMl9mD9HCUiMWx7PUBNjlEYc7ngFZBi1SkABap7ohPDlzLzxUPC5/PkYFmh+0VO8EpcG21QC1gLRjmKvcgL4EFB9ogX/XC1Fb7IvwRwm0O2fId/S4XaxgJ4rnPxVgu1FQVU+SmpAAWCWUkh2sZKg7gSY+ET6RaDcixLDkYrLyoqkAWXMvyqvf74EMKIAW+IR66kw3kIAcTIZIyZkn6PD9f+mqXAwVHEbsGY5jQ+LNINQpgkVCmcjiwsjiAX/BwDL7hzophOZ7DsbSF3gp8OoV1GKzcYx4AwId5AMdibA1zwP+glWaXhM/zIYagcxCz9JZJ9+VRXMPF4HM2RvFhXAEsUo0CWCSUaZ85n95FC56PpojbsscE5V3K4hF2GxLQGsFaKq/Bogu8a09wfklz+NSbyhGx1Y9yyuuJ4JU/aqnVIW04FRwF3OGQhd3pMneLSEABLBLK8TRRW+JQo/4ojmVbxT1GCR/f0B66B+wbPJitPDg8DT7kYBTAGjC4oR/o/ii+DfUfgHsdulwy8GKZ3wg6XHIwGoS/Z0jBY5oDFqlGASwSynCWDLQ7PFQ8Evt6EmPY0lTpvSkf39AZ+jyGy2ls1SXNOXDBwlqAqLOwqCqkZoeUz6UyC5sLbYuRhakch8r0r5ssHnwtDfDJCK0Ol3Ns0VmEItUogEVCeTRL0mfOo6P4oZn08A39FeeAPbiUpclyf7gMdh3mvSqfmV/8jAcRAMdgLLnQk8CdEVpdMiHu74/S7FJuRPnTcSKWWZ+tBiBmcAye4aj2IIlUowAWCSsNvRFaix+affN4kK2YN38+R7sLNuSUK/MezS5DlUtnJfKFuBb+5jp4hunQA78DMXzDlRDz0kNRMIyX+eQNcRzLRI4jPsA7moiG+1gRUQCLhDXj4Vo6ioegv5lm1mcgxpaKD9NImianRCmPkuIOjq12TIIPBjx4nW0OEUsWLoRe+tTkYiynQxyF1B/FtUyWCeDuCK7leDA93B8hazhYZf2YiIACWCS8KY+IYXvxdO8DLtYQt6ytuMbKdXAtGypOFV9lLNgQK6q8YBmzocXF1FKIozeKF+7I3ksZ0j4vz5W+mt8Ane/xbjcsi5CBPVoCLRKCAlgkrP0JsrC6eHvuM1kuZ7GWe9oqvfdMsoYDkQ7OMuvRX3UIOr+PaDndDo6DZ2gJvdC6J4ZnmAsxZL2thYkc42WavjpGBi5nATa4tDokfJ7WEmiREBTAImGdzZA1tLvsKh5JPpcm7dNScXg5DdayvTXUN0r7xFzmK+/qnQAPPMjxbAoPMHSFXmh9JYMX4jDgn4jQG6UtwtkyAdwVIe3z+CyAa0j71bYvi0hAASwS1ldTjOdocdlSnHPG0OrQXTH8zqbwCHseQ9YHQ3flSlgtkP+0doCcwVhengn1+UB/nKRftrbGVSmYypH02VsqgD/RTNwhERxEuCZCk3OtKLSIVKYAFqnBRI6Iobt4pPflObKG1orDv+NZPENfuDOBL6VI+kQqP53LWag8OQ8QccCG2lYE3BPHM0QcJqqF5aY4fTEulLntthZcw6nUtZsjtnqoi0ieAlikBmM5rGVHS9GLWZ+Mz+qK4XpsHiytLmtDPHOjaVxLFjZUCPWR4IsOgBx4sKml/P0FlsdpiXAxzQ/KLK26ak0TxpAqMxi+Mo5jGQ06x50uGXih2meKSJ4CWKQGr86T8a8/p+hPZ8gaopZ/2FH2jS+mSfm0ufxE+Xuuilh8S1eMlgoPaGfQA77IJhdrsbbKoUzXGDJ+qLqVt7fh2IU1Vostj5KBvbMLf10TI0eoA4ZFBAWwSE2OpsiWOljwYpqIg1sx/yazRMLNAXsw7+FbuisUwxoJ9gHHiVmMwVguh1sBtbqZJrdarWkAWl2SPiNlPrbZJQcHkwCPxIhaEj6Pag+SSDgKYJEafD/NnE/M4aPFA86XskQdNjVXem8OsKGO7J3IgMEYkhUWQq8LanFkOZgGg2+q7VwK9MWhaqEPALqiZOF3Rktc+nQ7MYcpjwM5gLUxopYxzQCLhKYAFqnNlQzWcEfxhqITCbJ+leKRY1kilu3t1b/FwTRZaI5we1f5m04GAexwSxTfMu8xnAj1I3THMIbxEMuV833xkna1EbEMB53jDXFcy7C6vyKhKYBFanMmjbXXj0IfT5CD9S1sLt/BPTqHZ1gRrhgWtto+3e5rX76WJuaS9tkTorQksKIFbPWjk/7PLlrcssVDBmI4ZuE8pZ0OO1rIwDOh90GJiAJYpDYnU3iwtjhHD86TgoilwphxwsM3eOHWSSU8Uj5DFca089uOfYBfWg4GG/ppbo+SgTeqLVfe2IJjyx7w0BMlBYfmAVa6tDjMeBxUD1gkNAWwSG1GUsx4dEV4R8EeoSM5xjJELJ8YKPvG43NkYTBcD3g0RUuE9sonEhow4DObJeIQCV2HMubiG56/UuW21c04lhOlhrU/0ESbS8bnwBzAO9uIWkYyC/PBIhKGAlikNn85z7xH1HJb8abb4RTW0lm+Hta+GdI+zS4/W2FmN+AZsBUXYQXpywVchyxhK2DsiOFYcrB7vsqdQ3EycGS21KUYMYfRLC9nAPqiOLZsvQ4RKUkBLFKzyVyJo432zeAbbi2/xupMjgtpIk6otcpTGRIePRXuzM/NGhgk7i7MGYexsnVhH3BlP9lMS4QM15ZZFVrbjGsZDRJ3dZwsjGj8WaQWCmCRmh2ZB8u64gnaySwZ6IxwS/mh4IkMUafi2uZAwqMpQmeFctD5SwYcxlILhTjCyPpEXSaqlctY3kSTS8LjVKkh6E0tYDgXZHNPhLS/sCBLREJSAIvU7HyaHAzFubtgwPnPJ5jziLnXB3OhwzNkTagKGMPz+Ib2CgHcBj74MEKTi2+x4ap8PDiE6zBfrQ2bWolYJrIcKRXVfTE8w8vTAL/cgQfTOf6HlkCL1EIBLFKzl2aY84hZ7imuK3k5jbOoZ1wo6pD0WBGiYvOJGXKGXIVR5fMA+LCJqQyOZTL8CLDhbLUJ4E1tOA7PjJe49MFWOqPMewtnOWxoptUtW65SRMpRAIvUbE+Gk0laI/QXr1IeTuIbdnaWfePBSSIOq0KcCuxaPEN/hW1IrZCP53HWdWAc5iqfHxzY2k3CY6xaWq9qIQejpW67vYO4w3CSHyUAbm8jC4+OhfruInKVAlhkKY7P48P2tqIXzyTAsKz8RiPHkvaJu3y4p8rnuxbfVBxVbgqGoHvpiWNMpS3IhVqjxKOcrDhcfGuMpggZeGWqxNU7OsFycBbgoy10uqR99lfrUovIdRTAIksxnMI39BTP0b4+S8qnP877y6yF/too0zkiDv3VdgMPz5H0MJY7yi2rngUDFqbJ+hBuDvi9PbREyfi8VKq881V3deFaEmVKa3VHyfgcnAHoixB3mchxWDuARWqkABZZiqenuJCmM8KuglHob00xmsZaVpYfOh6ex1q2lh+mzruYwDdE3fInDHYFW4EBi29D7QMengPLvMcbFc9NWt5SdlL5kS7aIkxlGUkC3N4Ohsta/yxSOwWwyFIczXJiHsfyK6uKXj8xR8zh3vIjzBNZMGyodiTDmTQpH8dhU7k9S1OQAx/swj7gWIhKWJu7iEfIVtsEvLUdY3m+1Aqs2zuJuUxkeT4JcGsbPjw1Uf1bi8h1FMAiS3QigW9oLY69lyfxYHP5fN07Rs7QXL5g1lUTaRyHleU+6gJY8KET32AsYRZBb+zCcZioduvKVsbSnC+1A3hDK8ZwbBbgkTbaI8x7HKpWVlpEFlMAiyzRS5Okfda38vGCTuqxGRIeMYctZco4775CwqOv4snBeTmwFSo8DwZfnMd1Qp1uBGztwXE4VzEvP9RD1MFY/vj89ZceaGZdK1l4agzgPT24lmPzPKcaWCK1UwCLLNGPE5xL0uyyqyCAn50j6dMc4cEypzIcSTKTJeLwuc1VPj/p4zisKTdbPBfUgo4Sj4AlEWISuCtOjio94J9aQWuU06VWNa9uIeow7/HYFMD2Njw4qO6vyJIogEWW7tlxkj4bi/f17p3Acbiz/DTwlRSRCBuqFaScyeBV6Nf2L0wA42ItxjJecV0VsKON9jg5w4sXK9022AyGM6Vi9ZYOog7HZwF+Ir5wWvCPS00Vi0hVCmCRpTufAkN/jNsLBpxPzOFRaZnVkxfIwbJq9bCmM2AZLPc5J8AHDwxdzeBUH4Je2U53M55htuKi5Y4YWUoH8L195OBvLgA82Euzw2iGFzX+LLIkCmCRpfufV5jK0hYtqn61b5zT87RH+Uhf6XfNZclCZ7UzkY6OkTPMlAvL3qAH3IHr4BumqgXhXAY3wkyWpy6Xvedjy+htIuXzh2evv/T3BuiJk/L56hWA1c14hq9X7EyLSAUKYJE35dgsjuW+3muvPDvDsRmiLu8bKv2WL59hNktzlH+wttInH7zMfK58D/gKWDAwiW/xQxxH+K41uG75RAfg3gEiDidKnQH8jh5cy9EZgJ1RtrcznuWYJoBFlkoBLPKmPD/OXI5NbWwoWK58dg4Pbuku+67xNFGXZRWLQnc247pEI2wpmcFbgi+G8QxYllXbW7y8A8fldMUilMta8A2vlJrWXd8OlsPTAPd1YS2X0zypCpQiS6UAFnlTvnSZpE9vnI+vvPbiwQky0B7jgY7S7xpNYh22VKwIfW6K2SzNMbb1l7p8tfTzVlriYDlcsboksKGHLLxW/radLaxuI+NzdlG/9sO99MSZzXJgAuChAWIOx5W+Im+CAljkzdo9jjFsKcjaJya4lKTZ5SeXl37LgSv4pnyRDQCOTuMbIlHWl5xLzs/4ejCJb/AsqWqnMbQ3k4Nnz5W9YWcf7TEmMxxeVNnqvgE6Y5xP8p1xdsVY1UQOnrpS5TuKSAUKYJE366+GScOmdt5ZsLB5/zie4fYy67CeGWE6S38rD5dJ6LzhWTI+uZKVI08HhaD7aI3jGZIV9wH/3c3EIszneKV82ch7luG6nJ5l36I54F29+IZ94wDv6aMpwrkk36s4mi0ilSmARd6s5+Y4M49ji1ZdvTRKyqctyt8ptRRr7ySXE7gOq8uMUeeNzGAcNg+WurZzYQ8SCbAYp8rqqhUdRCKMlTrd6KpNXWTh9UVHEP50Lx0x5nM8PgLw3gE8dX9F3jQFsMgN8PwV4m5R8Y1vXOD0LIMt9JY5efDkNDjcu6LSxyayxGK0laxbeQAMeJCgKYax9FZc0rV9CCznSy1vzvv59XTEmcny+KINSA8sozXC8Rl2z/JzvQw1kfbZO1np24lIVQpgkRvg358k6bOyldsK4nY8TcThQ2X2Gl2YYyzFxl5uLd8J3jdCztJVrmSHDw60kvJJ+mQrzgGv6cEzHCm/AuuOAZojnJ5h3/T1l+7sxzO8PAbw4eU0uRyd4YcafxZ5cxTAIjfG6TkiDp/eeO2Vr59gLkd3nJ8otdr528fxDbNZ1pSvSXlukhy0NLNx8T2bFx7f9TNEo7hR9i86O6FQS5yM4WD5uhkdTfiWp0auf/03tzHYwkyWx85xfwsb20n5/NHJSt9LRMJQAIvcGE9dxDfc2s09bQuvPHaJs3N0N/O+1SXuPzzNeJLlnTywpuxnWkPW0NbM8cUrp/JLrjx2rKK5ibTPq+V7t/cMEY+TyDJWZuPQJ9dzaz/zHrsXJfTGTqIOe65wJMnDy2mNcCXNU+WHskUkJAWwyI3x+6e5kKAnzq6CqlivTWAst5c7GWkMx2Gw/BD0XIbpNJEYv3jfomuZ/IGFdMSxEXAqte3dm4lEGU2wv0wRyp0DxCIcHGN/8QqsLU1s7SHl8+xFgHcNMJfjr4YrfS8RCUkBLHLDnJjFdXi4oCLHSxdJ+qxs5+FlJe7/xhFSHut62FRmQ/DhUc5PY12aFheOPgMWfEYvgoNfMYBvWYF12Fs+ODf14lsOXLr+9YdW0RpleJ6vDvPZtXTFSHj8QfmdxCISngJY5Ib567MkPHqauC/o1H5rmN0XiUb4O6VO/919mfEUXS383O1lP9O3GJeexQkdJHpLHONyoeKSqKEeMoazZXYAPzDImi4uzHF20Yf89FqM5dA4wP0D2KAWtIi8eQpgkRvmO6McGKcrzscLlmI9O0IWdg7yYKlO8PAMxmGgrcSlvPF5jMO6xVuBLyz8f0cvToRs+VZ9eAdtLWTgdJmDez91G9EokykePV30+v+1iWWtTKb57hnub2dlKzNZni5/kpKI1EQBLHIj7R4lC3cPclfQCf7yG1ycx3G4tVRJ5yfewLdsK5XNea8O4zmsXDyLHHxaLIZxMe6iGwLdrbQ0k/J58o3SNwy0MZfl6UUDy+9dTVOUI5P8cJRf2ExblOMzHCiT4iJSKwWwyI30Ryd5fYqWCJ/YdO3FU9M0x3hXqQ3BX3qFsST9HfziO0p/4OgMaWhdXGQjCj5AJAsusfKnC9+9GRvhtUX7i/I+voU13aQ9vrC/6PWHBuluZjbLo2/wni62dTGV4RtnOVixlpaIhKcAFrnBXroMlrsLBo2/cYyEx4oOdpXaEHxiHBvh3vWlP+31C2QtToyV120FzheeNGQimAgVylD2deJbnjxS+uqHtxGPc3DRFqaPbWF5O0cn+PYIH1xFW5QLSf5q0SotEVkyBbDIDfY7h5lMg+U3diy88vg53pigtYlfLtXN3T9CFtYPcM/KElePnCfp4UT5wF3FF/JdXo+1gxDBiZVuzHu2sn4FyRwTiRJXH1zNhn5SOX5UXFhjezs7BsDyNyfZ1co7l3ElxRMVC32ISK0UwCI33r4rNLm8qyBQv/oaaY+tg2xYNJj8hWe4PE93B1vKzARfnsGNs/a6qtFpMGBpa4YIw2WmZm/fQCTKTJbHXilx9f1baIlzdpqvHi16/ePbaW/iyBj/6wQ/s5aOGBn4/ROVfmQRqZUCWOTG+6+vkoHVHXx+18IrX3+Dk5P0tvEv3l3i/pEpnAj3bipxCRidhQhD1y2EblmYA85YiDJfZhn0bZtworx0vNSlPt6xBs/w+LGi1z+ylnevAcveS2xt4oEhjOHAWLmfVUSWSAEscuO9Nsfey/jwUMHCqx+fJpVjfR/vWnQC0lPHSPvsXMfdpYpWHjiB79B/3SLqN8AHj4E+iDJbahL47s1sWUPGcLhUCY73bKa/g8tz/MHuotc/sJGeFoan+X/38Es7GGzhUpJf21f9pxaRmiiARd4Sv/QM5+dpj/NbQSf4d3dzZorlXfzcHdff/PJpzk3Q1srfXVxyEhIZMtDbz13biy9YMExn8CyJUj3gLatpbWcmy58+UeLq+2/BWF4tLv78gdVsHyQD3zvFvV3cs4wc/HmpDrSIvEkKYJG3yvfPgOWDm1gT7BH6i1dI51jby6bi+s+vXmDPKSYT7NrIxkXbhV88yOgcze3Fo9CD4IFh5QqIMZcu0YA7t5OB46UWT/3q/SzvYT7H1w4UvX7/anpaGZnld/fy97bSGeN8gj/W7K/IW0ABLPJW+Q8HODNNZzP/8v6FV/78EG+Ms7aPf/qe62/+t98El74efvK26y8dPc14AhvnjsJLF8GCR8rBNJEudRjwjq3MZHi61PKrezYScdk/wnNnr7147yAPbWY+y1de5cMreOdy0vDXp2r7qUUkJAWwyFvoG6+TynHHEI8E23wfO0I6x87VfHxR0L5wHOPwvjtLfM6pC5g4Q4X7lIYWVkG3tpExnF+0SOo3foHeAfwIBxYNIP/C/dyyioTHo8Xd38/cRXsz52c5Msqnd9AS49gkf3D0+reLyA1R8QgVEXlz9o3yofUsa8NY/uYEwFyCnatY04t1+Pahopvb4tyxif5ekileKd6Yu2oV995HJMbz/4td8BA82Mq7EjwwQjzH2EVO7+NQcZWMZ/exfiUHjvL1H1zfqs//LMs62HuG//D9ay/+o3fwoe1MzvPfd7OmiQ+sI53jP+3huE5fEHlrKIBF3lodUbb0s66HRJIDo4wlyWR4YAsDnVye5GjB2QZHhnnvLnp78A2PvVj0IftfxWR48f/mfbAL1sHgLINn6D+N/wKZH7P8Eu+GPjhY8K7vP8dzxQUmgc99kHs3k8ryhz/kWEFmf+69rOjk2CjfPcSvvYPOOM8O83uHrn+7iNwoCmCRt9aeS+zsZ10PO5ZxZYaj4xwd5YENrOpj63KePcJEQXXlnk7u2c5AHxdGeb1g79BdcMcBboF+aIdmaIZ48CcGTdAMQ7AdklChaNW/fITBTn50mP+voPv77x7m/vWMzvL7T/OJLWzv4+QUv7uH86XqZ4nIDaE5YJG33K98n2NjdLfySDDv+4k/YXyeoW4+dm/Rnf/tbzgyQlsHn/rItRfXwiOwCrqhDVqhBZqhpeCLVmiDbtgMH4VSpw8D/MMHWdHHVJqvF/Sw37ueB7eC5ZlTNHvcOkDW588OsUcHH4m8ldQDFnk7DLawsZ+BNrIZ9l8EWNXN9tVsGGJ8miMXrt0Zi7FrO4NDzCd45RjAZ2Eo6PjGIQpRiIAb/HGK/8RhNZyH6xZmbV7Gv/4Urc1852X+5Mlrr3/+Q6zr4/QYn/kqv/sQKzo5coXfeuGt/0cRubkpgEXeDi+NcGs/WwdZ18uX9gE8dYyP7GKgl9X9nL3M2aC7efANNm9g9SrWruHPvsk/gY1B+sYgBpGC9M0nrg3+mOBPFLrgugz90m+wqpcTl/jMF669+EsP8KGdzKf5zkE+voE7hhif598/y2mtvRJ5iymARd4m0/Pct57BDu5eyaOHAGYS7FrPYC99nXxrz7U7xyZ530/Sv4y1yxl8hs6Cvm/V6DXBh7gwBeeCv37+M9yzg/kMv/lFzgaHDz5yO7/6EPEoTxxmdowPbsFa/uwVvvL62/RvInIzUwCLvE3OzTA9z11rWNHNUDtPvsHRC1jDnVvo76YlzgvBjtsLo6xcyfqNJF7HfY5WaCru+5ZL3zwfPPAhCy8DcM8t/PLHaI7znRf4k+9ea9LvfJwV3Rw5zzee5R/dS1ucH5/iN59+O/9VRG5eCmCRt8/hUfpauG0Va3oZn+XwJfae4q5NrB5i4wrmUxw8vXDn5Cy33ML5z9IdDD5Hiqd7r0vfq/wggHMAPAH338Gv/zxrBjl1hed7dQAACgtJREFUif/yl9dKdnz5l9m+ksvT/Isv8W8/wEAbZyb44z2cmnob/0VEbmIKYJG31TOnePcGVveybSVvXODsJIdPc+cWVi1j/Uq+/L2F2y5e4vKj3Aod0ATRUiPPpkz65gM4Cxm4BP/419mxhUvjfO4LvBycPPif/j73bSWT5T8/zidvY0s/4wl+72keV+FJkbeLAljk7XZ6lLvWs7ybnWt54zyvjHD2Ig/sorebe3bw6I8XbrsFdkALxIsHn8uNPBOkbw5ykIEE3P5b3HMviSRf/BrfCcaW/+kH+dgD+D7feIlbO7h7DekcX9nHn5aqGi0ibxHtAxZ5u+0Z4QtPMDLN6gH+zccAnj/KN54ma7hjO//j3yzcllk47giKU9YUvHjdFxTcbMCBLdvxXB5/hi/+1cLVX3yYT70PN8IrZ2lNccdq0h5ffJHfe/4t+nFFpDT1gEXq4Nhlrkxz53r6u3jnJg6e5tsvMdTPhrWsWsmt2zh1jsg4twUTwG5B99cE/1vIC0ae84PPWUjBHCz7NPsP8c//1cJtn3wvv/oR2lrYfZQLp3h4BzMpfnCU33nq7f4XEBEFsEh9HL/MfIqtK9iwnG1rOHyGv3yKzg56eti+hS2b2PcYm6A1mAB2Fq26yv+vX7DsOb/2Kj/7m4JJeLGDv/wmFy8D/Nr/wS9+hLZm9h3j9SN89HaMZfcpPvvX9ftXELmJKYBF6ubQCC0xNq9kzSCbV3F+jL/4AVvWs3Ydy1ew8j7Gv70QwIvTN88vXnh1NX3TMA+j8Nt7FtL3cz/PJz9Aa5zdRzl7lI/cgbE88Rr/5Gt1+uFFbnoKYJF62nMS13DrelYMcN8OZhP8wVdxI9x2GytXkVpH8ocLAeyUGnn2C0aec0EApyEFs3AQ9gLw//8WH3wXjuHlw0yc4uGdGHj+OL/ylTr8yCKSpwAWqbO9J5maZ8tqhgZ4xzY6WvnCnzE7x8AQT+8hup9mcEste/YWTf1mCrq/E/Ad+NAn+e1/zs6tzCX57jOstty5nkyO7+znnyl9RerKVL9FRN56D+3kV3+W9avAZf9xPv/fOT4M8I9hJ3RCS0E5DlvQG74ugPNrryZhH6z7bd51F+2tjF3i0b/mfetY2c1cgj/5Af/5+xVbIyJvPfWARf5WOHmRv3iSDStZMciqFTx0P81N7DnERdgETcXjz9cNOxeOPCdgFo7Co/CxT7FsGadO8dpzvP8WBtq5NMkfPs5/+1E9f1IRyVMAi/wt8r3dtDSxZiWDfdy2nTdOsX+Yc7ABokEAe4smfTMF6TsDZ+ErcArOvE5vM+ujbFmJYzh8hi98k6+/VOefUUTyNAQt8rfO7Zv5mQe5ayc/9QvXXvwsrIf2oDCWDcro5MM4A0mYg7Pw/wDwmY/yoftojbGskyvjfPUJ/su36/HDiEgZCmCRhvEx2AYDLCzLyo9f5TvBCRiDp+AJ+Pw/476d9LbRFGVynB8+z2/8UX0bLiIluPVugIiEld+yezesgjW92NuIzTJ+jsR5jsJzwW3LB1nWz9wse/fx6X9dt9aKSGUKYJEG8xK8BETgTvDgP15/w4+eY+QC3/oBew/XoXkiEpKGoEUa1q+Dhe+CTjESaUA6DUmkMa0BwIPR+rZDRJZIQ9AijSkVnIskIo1Jj69IY2oBG9SiFJEGpAAWaUwz4IGBbL1bIiJLogAWaUwdweMbrXNDRGRpFMAijWlWg88ijU0BLNKYYkEAN9W5ISKyNApgkcaUXwXt6yEWaVTahiTSmExQRydX54aIyNLol2eRxhSHHPiQqndLRGRJFMAijSkJgIH5OjdERJZGASzSmPK7j3wYr3NDRGRpFMAijSkZzAFvrnNDRGRptAhLpDFdCAJYu4FFGpN6wCKNaQMAHlypc0NEZGkUwCKN6eroVV89WyEiS6YAFmlMR4MqHJfr3RIRWRIFsEjD8sGHjno3Q0SWRAEs0sgMJOrdBhFZEgWwSGNaDhYMROrdEhFZEgWwSGPqCzYgeXVuiIgsjQJYpDFdPQ0pWe+WiMiSKIBFGtORoO/bW+eGiMjSKIBFGtNGcMAP6mGJSKNRAIs0pnnwNAEs0sAUwCKNaQQAVz1gkUalABZpTMuCIejperdERJZEASzSmPLjzz401bslIrIkCmCRxhQDwEJznRsiIkujABZpTBlAxxGKNDAFsEhjynd8DVyoc0NEZGkUwCKN6URQilKFOEQakwJYpGHlS1H21LsZIrIkCmCRxnQr+GBhpt4tEZElUQCLNKZJMJALBqJFpNEogEUakxsUgtZDLNKY3Ho3QESW5GRQhHKkzg0RkaXRL88ijalPVaBFGpsCWKQxjQazvyvr3BARWRoNQYs0psEggHUioUhjUg9YpDF1BIuwEvVuiYgsiQJYpDEdC76I1bMVIrJkCmCRxjQU7ERSLWiRxqQAFmlMLmTBh6F6t0RElkSLsEQaUwwsGG1GEmlUCmCRxjQNPvgqxCHSqDQELdKYtPZKpMEpgEUalq/xZ5EGpgAWaUw2OA94Y71bIiJLogAWaUyZIIPn690SEVkSBbBIw8ovwpqudzNEZEkUwCKNKRXUgm6vc0NEZGkUwCKNqQts0AkWkQakfcAijWYnAA74YGEFROE0rIHT9WyXiNREmxhEGsTfhzUwHvzVgzgYSBY/x0n44tvfOBGpmXrAIg2iA3ZAD7jQBnPQBEmIgQOj4MIYfL/e7RSRcDQHLNIg9kIKchABHzogDp0QAx/aIAeJgmMKReRvNwWwSIN4AaYhC16pqx5kYRKeebvbJSJLowAWaRxHIQNe8eLn/NceZHQwg0gjUQCLNI5DkITcoq1HPuQgCUfr0y4RWQIFsEjjeBLmIRd0gvPy3d8cTMMr9WydiNREASzSUM5CtlQPOAvn4Ux9GiUiS6AAFmkohyFdMA18dQI4Ba/Vu20iUgsFsEhDOQezkCt+MQezcLY+LRKRpVEAizSU5+Fc8Tqs/AqsKXixnu0SkVopgEUazaVgM1KeB2l4vZ4tEpElUACLNJrHYLa4B5yCI/VskYgsgQJYpNEchUTBNHAOxuB79WyRiCyBAlikAR0LpoHzE8Bv1Ls9IlI7BbBIAzoGqeDrtApgiTQkBbBIA/oLSAQ94Cn4er3bIyK1UwCLNKazQQAP17slIrIkCmCRxnQgOJrwcL1bIiJLogAWaUx7IQUJeKLeLRGRJXHr3QARWZKX4TK4KsEh0qgUwCIN6xwk690GEVkqBbBIw9oNpt5tEBERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERkRvhfwNRz9lgXUxESAAAAABJRU5ErkJggg==" } }, "cell_type": "markdown", "id": "d6f7e47b", "metadata": {}, "source": [ "![dipole-gauss.png](attachment:dipole-gauss.png)" ] }, { "cell_type": "code", "execution_count": 16, "id": "1e9bfe06", "metadata": {}, "outputs": [], "source": [ "x,y,h, r, th = sp.symbols(\"x y h r theta\", real=True)" ] }, { "cell_type": "code", "execution_count": 3, "id": "b3238233", "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle - \\frac{1}{x^{2} + \\left(h + y\\right)^{2}} + \\frac{1}{x^{2} + \\left(- h + y\\right)^{2}}$" ], "text/plain": [ "-1/(x**2 + (h + y)**2) + 1/(x**2 + (-h + y)**2)" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "e = 1/(x**2 + (y-h)**2)-1/(x**2 + (y+h)**2)\n", "e" ] }, { "cell_type": "code", "execution_count": 11, "id": "3bf8bd25", "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\frac{2 x}{\\left(x^{2} + \\left(h + y\\right)^{2}\\right)^{2}} - \\frac{2 x}{\\left(x^{2} + \\left(- h + y\\right)^{2}\\right)^{2}} - \\frac{- 2 h - 2 y}{\\left(x^{2} + \\left(h + y\\right)^{2}\\right)^{2}} + \\frac{2 h - 2 y}{\\left(x^{2} + \\left(- h + y\\right)^{2}\\right)^{2}}$" ], "text/plain": [ "2*x/(x**2 + (h + y)**2)**2 - 2*x/(x**2 + (-h + y)**2)**2 - (-2*h - 2*y)/(x**2 + (h + y)**2)**2 + (2*h - 2*y)/(x**2 + (-h + y)**2)**2" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "e2 = (sp.diff(e,x)+sp.diff(e,y))\n", "e2" ] }, { "cell_type": "code", "execution_count": 22, "id": "c4de9567", "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\frac{2 r \\cos{\\left(\\theta \\right)}}{\\left(r^{2} \\cos^{2}{\\left(\\theta \\right)} + \\left(h + r \\sin{\\left(\\theta \\right)}\\right)^{2}\\right)^{2}} - \\frac{2 r \\cos{\\left(\\theta \\right)}}{\\left(r^{2} \\cos^{2}{\\left(\\theta \\right)} + \\left(- h + r \\sin{\\left(\\theta \\right)}\\right)^{2}\\right)^{2}} - \\frac{- 2 h - 2 r \\sin{\\left(\\theta \\right)}}{\\left(r^{2} \\cos^{2}{\\left(\\theta \\right)} + \\left(h + r \\sin{\\left(\\theta \\right)}\\right)^{2}\\right)^{2}} + \\frac{2 h - 2 r \\sin{\\left(\\theta \\right)}}{\\left(r^{2} \\cos^{2}{\\left(\\theta \\right)} + \\left(- h + r \\sin{\\left(\\theta \\right)}\\right)^{2}\\right)^{2}}$" ], "text/plain": [ "2*r*cos(theta)/(r**2*cos(theta)**2 + (h + r*sin(theta))**2)**2 - 2*r*cos(theta)/(r**2*cos(theta)**2 + (-h + r*sin(theta))**2)**2 - (-2*h - 2*r*sin(theta))/(r**2*cos(theta)**2 + (h + r*sin(theta))**2)**2 + (2*h - 2*r*sin(theta))/(r**2*cos(theta)**2 + (-h + r*sin(theta))**2)**2" ] }, "execution_count": 22, "metadata": {}, "output_type": "execute_result" } ], "source": [ "e3 = e2.subs([(x, r*sp.cos(th)),(y, r*sp.sin(th))])\n", "e3" ] }, { "cell_type": "code", "execution_count": null, "id": "fa5a4f76", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.10" } }, "nbformat": 4, "nbformat_minor": 5 }