{ "cells": [ { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "

The Quantum Eigenvalue Problem

\n", "\n", "The quantum eigenvalue problem is solved by matching the left and right wavefunctionsfrom a potential well, and using bisection algorithm. The plots are made with MatPlotLib. " ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Iteration, E = 0 -17.077272727272728\n", "Iteration, E = 1 -16.26590909090909\n", "Iteration, E = 2 -15.860227272727272\n", "Iteration, E = 3 -16.06306818181818\n", "Iteration, E = 4 -15.961647727272727\n", "nL 501\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEDCAYAAAA2k7/eAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3de3xU1d3v8c8vM7knEJBbINxFBBS5xCBqbVFLLVqoWqsea6PYR3xOta31VPHxaa32ctqjr3pqbY/QVq08itpaKk/FG2rBSpWLQAS5X5RAuMmdQEgy6/yxJzFArkxm9ly+79drXjOzZ81evzUMv6xZe+21zTmHiIgkvzS/AxARkdhQwhcRSRFK+CIiKUIJX0QkRSjhi4ikCCV8EZEUEfcJ38yeMLOdZraiFWUvMrMPzKzGzL52wmulZrYufCuNXsQiIvEp7hM+8BRwWSvLfgLcBDzbcKOZdQbuB8YAJcD9Ztap/UIUEYl/cZ/wnXPzgT0Nt5nZQDN71cyWmNk7ZnZmuOxm51wZEDphN18C3nDO7XHO7QXeoPV/REREkkLQ7wBO0XTgNufcOjMbA/wOuLiZ8r2ALQ2el4e3iYikjIRL+GaWB5wP/NnM6jZntvS2RrZpTQkRSSkJl/DxhqH2OedGtOE95cAXGjwvAv7RjjGJiMS9uB/DP5Fz7gCwycyuATDPOS287TVgvJl1Ch+sHR/eJiKSMuI+4ZvZTOBfwGAzKzezW4AbgFvMbDmwEpgULnuumZUD1wDTzGwlgHNuD/ATYFH49mB4m4hIyjAtjywikhrivocvIiLtI64P2nbp0sX169fP7zBERBLGkiVLdjvnujb6onMubm+jR492IsmotLTU4U0NPu42ZswYv0OLe9u2bXPXX3+9Gzx4sEtLS3OlpaUnlXnyyScb/XyPHDnS7L7nzp3rxo4d6/Ly8lyPHj3c3Xff7aqrq48rEwqF3COPPOIGDx7sMjIyXI8ePdw999zTpvicc27//v3ujjvucIWFhS4jI8MNHDjQPf/8823+PE4ELHZN5NS47uGLJLNLL72UGTNmHLctIyMjqnXW1NQQCARocA5LwqmqqqJLly5MnTqV6dOnN1kuJyeHDRs2HLctKyuryfJlZWVMmDCBe++9l6effpqtW7dy2223UVtby8MPP1xf7q677uLvf/87Dz30EGeffTb79++noqKiTfFVV1czfvx4OnXqxAsvvEBRURHl5eVkZrZ0SlGEmvpLEA839fAlWZWWlrrLL7+82TKAmzZtmvva177mcnJyXP/+/d2MGTOOK1NeXu6uvfZaV1BQ4AoKCtyECRPc2rVr61+///773bBhw9yTTz7pBgwY4NLS0tzBgwfdoUOH3I033uhyc3Ndt27d3M9//nN3+eWX1/dGH3jgATds2LCTYjr//PPdHXfcEfkH0E4axtzQk08+6XJzc9u0r3vvvdeNGDHiuG2zZ892WVlZ7sCBA84551avXu2CwaD76KOPIopv2rRprn///q6qqqpNMbYGzfTwddBWJI49+OCDTJo0ieXLl3PttdcyefJkPv74YwAqKysZN24cWVlZzJs3j3/9618UFhZy6aWXUllZWb+PTZs28eyzz/LnP/+Z5cuXk5WVxV133cW8efOYNWsWb731FsuXL+edd96pf8/kyZNZvXo1CxcurN+2Zs0aFixYwC233ALAM888Q15eXrO3Z555Jkaf1MmOHDlC3759KSoq4oorrmDp0qXNlq+qqjrpF0B2djZHjx5lyZIlALz00ksMGDCAV199lQEDBtCvXz9KS0vZuXNnm2L729/+xgUXXMAdd9xBjx49GDp0KD/+8Y+prq5uWyPbqqm/BPFwUw9fklVpaakLBAIuNzf3uNvdd99dXwZwU6dOrX9eXV3tsrOz63v5f/zjH93pp5/uQqFQfZmamhrXuXPn+rHg+++/3wWDQbd9+/b6MgcPHnTp6elu5syZ9dsOHTrkCgoKjuuNXn755W7KlCn1z++++27X8P/kgQMH3Lp165q91fWMo6WpHvSCBQvcU0895ZYuXermz5/vrr76apednX3cr58Tvfbaa87M3NNPP+2qq6tdeXm5+9znPucA9+yzzzrnnJsyZYrLzMx0JSUlbt68eW7+/PmupKTElZSUuNra2lbHN3jwYJeZmeluvvlmt3jxYveXv/zFde/e3d11112n/FnUQWP4IvHnoosuOmmMt6Cg4Ljnw4cPr38cDAbp2rVrfW9yyZIlbNq0ifz8/OPeU1lZedzYdVFREd27d69/vmHDBqqrqykpKanflpuby1lnnXXcfv7t3/6N0tJSHnnkETIyMpgxYwY//OEP61/Pz88/qe72lpeXV//4G9/4Bo8//nir3jd27FjGjh1b//z8889nxIgR/OY3v+HRRx9t9D3jx4/n4Ycf5vbbb+fmm28mMzOTH/7wh7zzzjsEAgEAQqEQVVVVzJgxgzPOOAOAGTNmMHjwYBYtWsSYMWNaFV8oFKJbt278/ve/JxAIMHr0aD799FPuvPNOHnrooagdY2mXhG9mlwG/BgLAH5xzvzjh9UzgaWA08ClwrXNuc3vULZKocnJyOP3005stk56eftxzMyMU8lb/DoVCjBgxgueee+6k93Xu3Ln+cW5u7nGveZ1AWkwql19+OTk5Obz44ot07NiRffv2cf3119e//swzzzBlypRm9zFt2jRuuOGGZss0Z9myZfWPO3TocMr7CQQCFBcXs27dumbLff/73+fOO++koqKCTp06sXnzZu6991769+8PQGFhIcFgsD7ZAwwaNIhgMMgnn3zS6oRfWFhIenp6/R8SgCFDhlBZWcnu3bvp2rXxWZWRijjhm1kA+C3wRbxFyhaZ2Wzn3EcNit0C7HXOnW5m1wG/BK6NtG6RVDZq1ChmzpxJly5dTvpl0JzTTz+d9PR0Fi5cWJ/IKisrWbFiBQMHDqwvFwwGuemmm3jiiSfo2LEjV1111XH1TJw4scUE1/CXxalo6Q9iaznnKCsr45xzWlp2y/tD2LNnTwBmzpxJ7969GTVqFAAXXHABNTU1bNiwof6z2rhxIzU1NfTt27fV8VxwwQU8++yzhEIh0tK8Q6lr164lJyeHLl26tLV5rdYePfwSYL1zbiOAmT2Ht7ZNw4Q/Cfhx+PFfgMfMzFxdV6OdPfrmOmpqG1wDpUFPpmGfpmEHx05YQfn415rY3kwPqal9n/iW1sTTXEfMWtW2Zt7TVLkTKm1s3wEzuuZn0rMgm0Hd8ggGNAegLaqqqti+fftx2wKBQKt7dzfccAMPP/wwkyZN4sEHH6RPnz5s2bKFl156idtuu41BgwY1+r68vDwmT57MPffcQ5cuXSgsLOSnP/0poVDos+/GwR2wey3furCQX/7yH6RZGq//5k5Y9If6/eQ7R5MDOnX/tT+t33Dya5z6a8vWlgOOA+WrSduXzbKn7iYjGGRo/x4APPDEK5w3rC+Dirpy4PBRHv3LfMqWL+P/TbkI/vkIAI+9+A6P/fWfrP6vqfX7fWjmW1xWciZpacZf55Xxi6ff5IUff4PAu48AjkszQow6oxeTr/4i//f2iQB877HZjBnSm+LDb8G8f4BzLFu/DYADW1aStjebZX/4DhnBAEP7dgPg38+p5bFfb+e7V1/A7V89j8079nH/Q3/lf14xCnv755CRCxd+r6lP95S1R8Jv7OIiJ/7Zry/jnKsxs/3AacDuE3dmZrcCtwL06dPnlAJ6fN4GjlTX4tV3SruQNsrPDDJxRE++e8kgunVoeq6zfGbu3LkUFhYet61Xr16Ul5e36v05OTnMnz+fqVOncs0117B//3569uzJuHHj6NSp+St4Pvzwwxw+fJiJEyeSl5fHnXfeyY4dO8g6tAV+dz7sXAnAAODzfdL4eH+IL2yfBi/Hx/z9kQ8cOO75f7+7gr4djc3f8/4E7Vt5lFtfrGb7IUfHTGNkYRrzS7Mo2fYEeLmY3UuPsuaTY/DmA/X7eeWVw/zsiZepqoVzugd46esZfLl6Frw1C/DWovn7FSG+8+pRLvr2r8kOGl8cGOBXE7JIm/e/m45vwarj4usNvH59gO+/togRf3+PHnnG5OHp/OcZZTD/Q8jrFpWEH/HiaeFlir/knPtW+PmNQIlz7o4GZVaGy5SHn28Il/m0sX3WKS4udosXL44ovqY0bPeJH4Frqtxx2098T9P7a+t7WlP/Sfs77v3t1zZ3Yq0NnlaHHDsPHOWTPZXMW7uLvy+voEN2kBm3jGFI4amPt0rsVVUepm/PrvygJMRdVxXDWVdDz5GQ35Ohn7uCG75+Nff9IJyATvrJedzP3sa3N/q+pl5r7n1J8FqUT3ozsyXOueLGXmuPHn453h+sOkXU/w09qUy5mQWBjpxwndpYa2pYo5GSUY8lkfUqyGZkn05MGtGLf//8QG7840JueWoRr3zvIjpmp7e8A/HF0qVLWbVqFSUlJRw8eJBffv+bHKw8wrVT/hOuehDM2LlzJzOfncnmT8qZ8p27ID96Y8sSG+2R8BcBg8ysP7AVuA74HyeUmQ2U4q1r/zXgrWiN34t/BnXP5/EbR3Pl797lt2+v5z8mDPE7JGnGr371K9asWUMwkMaITpXM/9mVFF39k/rXu3fvTpcuXZg2bVpUDyRK7ESc8MNj8rfjXUEqADzhnFtpZg/inQAwG/gjMMPM1uP17K+LtF6JTyN6F/CV4T2Z+f4nfPeSQeRm6lSPeDRy5Ejqh0tf/Q9YOA2m/Pa4MuqTJZ92+d/onJsDzDlh248aPD6KdxUqSQE3ju3L7OXbeP2j7Vw5ssjvcKQ5oVooex7OvAI6FLZcXhKa5tFJuxvdpxNd8jJ5c1Xb1hcRH5QvhsrdMHSi35FIDCjhS7tLSzPGDe7KO+t2EwppWCCubXgLMBh4id+RSAwo4UtUFPfrxP4j1Wz+9LDfoUhzyhdBt6GQ3fozdSVxKeFLVJzT20sgy8v3+RyJNCkUgq1LoKjRKduShJTwJSoGdcsnOz3Ah+UHWi4s/ti7CY7ug16j/Y5EYkQJX6IikGb075LLxt2H/A5FmrK9zLsvHN58OUkaSvgSNQO65rJxl8bw41ZFGaQFvTF8SQlK+BI1A7rkUr63kqqaWr9DkcZsL4OuZ0IwyhfOlrihhC9R079rLiEHW/ZUtlxYYq+iDHpoOCeVKOFL1PTsmA1Axf6jPkciJzm4HQ7v1Ph9ilHCl6gpDCf87Ur48acifMBWPfyUooQvUdOtgzc2rIQfh7Yv9+57nO1vHBJTSvgSNVnpATrnZrD9gBJ+3Kkog079IUsXqkklSvgSVd07ZKmHH4+2l2n8PgUp4UtUde+Qyc6DVX6HIQ0d3Q97N2v8PgUp4UtUdc7JYG/lMb/DkIa2f+jdF57jbxwSc0r4ElUFORnsq6z2OwxpSDN0UpYSvkRV59x0DlXV6GzbeLK9DPK6Q353vyORGFPCl6jqlJsBoF5+PNEZtilLCV+iqnOOl/D3HNY4flyoPgq7VmuGTopSwpeoKggn/L1K+PFh50fgatXDT1FK+BJVncNDOns0Uyc+aA38lKaEL1FVkJMOwIEjNT5HIoA3JTOzAxT08zsS8YESvkRVXmYQgENVOmgbFyrKvPVz0vRfPxXpX12iKicjgBkcPKoevu9CtbBjhcbvU1hECd/MOpvZG2a2LnzfqYlytWa2LHybHUmdkljMjLzMoBJ+PPh0A1RXavw+hUXaw58KvOmcGwS8GX7emCPOuRHh28QI65QEk58Z5FCVEr7vtusM21QXacKfBPwp/PhPwFcj3J8kofysdA6ph++/iuUQyISug/2ORHwSacLv7pyrAAjfd2uiXJaZLTaz98ys2T8KZnZruOziXbt2RRiexIO8LPXw48L2Mug2BALpfkciPgm2VMDM5gI9GnnpvjbU08c5t83MBgBvmdmHzrkNjRV0zk0HpgMUFxe7NtQhcSovM8i+I5ql4yvnvBk6Q67wOxLxUYsJ3zl3aVOvmdkOMyt0zlWYWSGws4l9bAvfbzSzfwAjgUYTviSfvKwgW/ZW+h1GajuwFY7s0fh9iot0SGc2UBp+XAq8dGIBM+tkZpnhx12AC4CPIqxXEkiHrKDG8P1WtySy1sBPaZEm/F8AXzSzdcAXw88xs2Iz+0O4zBBgsZktB94GfuGcU8JPIZqWGQe2lwEG3Yf5HYn4qMUhneY45z4FLmlk+2LgW+HHC4CzI6lHEltORpAj1bWEQo60NPM7nNRUUQZdBkFGrt+RiI90pq1EXXZGAICqmpDPkaSw7VoDX5TwJQZywgn/SLWueuWLyj2wf4vOsBUlfIm+rHQv4Vce0zi+L3SGrYQp4UvUZYcT/lH18P2hGToSpoQvUVeX8I8c0xi+L7aXQYciyOnsdyTiMyV8ibpsjeH7q6JM4/cCKOFLDNSN4Svh++BYJXy6TuP3AijhSwzUz9I5poQfcztWgguphy+AEr7EQP0YfrVm6cTc9uXevXr4ghK+xED9GL4O2sZeRRlkd4KORX5HInFACV+iTmP4Pqo7w9a0pIUo4UsMaB6+T2qrYcdHGr+Xekr4EnUZwTSCaaaDtrG2ey3UVkEPnXAlHiV8iYns9ICGdGJt2zLvvucIf+OQuKGELzGRmR6gUj382Nq2FDLyofNAvyOROKGELzGRGUzjmJZHjq2KZd74fZr+m4tH3wSJicxgGsdqlfBjprYGtq+AQg3nyGeU8CUmMoJpHKvRkE7M7F4DNUeg50i/I5E4ooQvMZERTNMVr2JJB2ylEUr4EhMaw48xHbCVRijhS0xkKOHHlg7YSiP0bZCYyAhoSCdmdMBWmqCELzGhHn4M1R+wVcKX4ynhS0xkBgOalhkr9QdsNUNHjqeELzGREUyjSksrxIYO2EoTIkr4ZnaNma00s5CZFTdT7jIzW2Nm681saiR1SmLK0IlXsaMDttKESL8RK4CrgPlNFTCzAPBb4MvAUOB6MxsaYb2SYDI1Dz82dMBWmhGM5M3OuVUA1vzFFUqA9c65jeGyzwGTgI8iqVsSiw7axogO2EozYvGbrxewpcHz8vC2RpnZrWa22MwW79q1K+rBSWxkhqdlOuf8DiW51R2wVQ9fGtFiwjezuWa2opHbpFbW0Vj3v8n/9c656c65YudccdeuXVtZhcS7jKD3VauuVcKPqm1LISMPTjvd70gkDrU4pOOcuzTCOsqB3g2eFwHbItynJJjMoHeZw2O1ofrkL1FQsQwKz9EBW2lULL4Vi4BBZtbfzDKA64DZMahX4khdktfUzCjSAVtpQaTTMq80s3JgLPCymb0W3t7TzOYAOOdqgNuB14BVwAvOuZWRhS2Jpi7ha2pmFOmArbQg0lk6s4BZjWzfBkxo8HwOMCeSuiSxZdYlfM3UiR4dsJUWaKBPYiIjmIYRouagZl5FjQ7YSguU8CUmMgJpPJn+EIVzbvI7lOSlA7bSAn0zJCYG98inc7/hZH/6EdQc8zuc5KMDttIKSvgSE31Py2X4mHFYbRXs1EnW7U4HbKUVlPAldnqN9u63LvE3jmS09QPvXksiSzOU8CV2CvpCzmmw7QO/I0k+W5dAVkctiSzNUsKX2DGDnqM+641K+9m6xOvd64CtNEPfDomtXqNh12qoOuR3JMmj+gjsWPnZkJlIE5TwJbZ6jQYXgorlfkeSPCrKwNUq4UuLlPAltnqN8u514Lb91H2WSvjSAiV8ia3cLlDQRwm/PW37ADr0gvwefkcicU4JX2Kv12gduG1PW5d89stJpBlK+BJ7PUfB/k/gkNbViVjlHtizUcM50ipK+BJ7dclJ8/EjV/cZKuFLKyjhS+wVngOWpnH89rD1A8C0ho60ihK+xF5mHnQ9U+P47WHrEug6GLI6+B2JJAAlfPFHr1FesnK6qPkpcy58wFbDOdI6Svjij16j4cge2LvZ70gS1/4tcHiXFkyTVlPCF3/0DE8j1IHbU6cTrqSNlPDFH92HQSBT4/iR2LoEAhnQ/Sy/I5EEoYQv/gike7N1NFPn1G39AHoMh2CG35FIglDCF//0GuUtolZb43ckiSdUC9uWaThH2kQJX/zTqxiqK3XJw1Oxaw1UH1bClzZRwhf/FBV79+UL/Y0jEemArZwCJXzxT6d+kNsVtizyO5LEU74ofEnDAX5HIgkkooRvZteY2UozC5lZcTPlNpvZh2a2zMwWR1KnJBEzKCpRD/9UbFkIRefqkobSJpF+W1YAVwHzW1F2nHNuhHOuyT8MkoJ6n+ut9nj4U78jSRxH93uXiew9xu9IJMFElPCdc6ucc2vaKxhJQUUl3n25hnVarXwx4LwevkgbxOr3oANeN7MlZnZrcwXN7FYzW2xmi3ft0nrpSa/nSEgLalinLcoXAaYDttJmwZYKmNlcoLFrp93nnHuplfVc4JzbZmbdgDfMbLVzrtFhIOfcdGA6QHFxsVbWSnYZOd6ZoluU8Ftty/vQbahWyJQ2azHhO+cujbQS59y28P1OM5sFlNC6cX9JBUXnwrJnvROwAi1+JVNbKATlS+Csq/yORBJQ1Id0zCzXzPLrHgPj8Q72inh6l3gnEe1a5Xck8W/3Gqja731mIm0U6bTMK82sHBgLvGxmr4W39zSzOeFi3YF/mtlyYCHwsnPu1UjqlSRTd/BRwzot2/K+d1+khC9tF9HvZ+fcLGBWI9u3ARPCjzcC50RSjyS5uhOwyhfBubf4HU1827IIsjvDaQP9jkQSkM7aEP/VnYClHn7LysMnXJn5HYkkICV8iQ+9z4U9G3QCVnMq98DutRq/l1OmhC/xoW5MeqtW3mhSefizUcKXU6SEL/Gh5wiwgIZ1mlO+ECzts8tDirSREr7Eh4xc6HH2Z7NQ5GRbFnqXhszM8zsSSVBK+BI/+pznrfNeW+13JPEnVOt9NlowTSKghC/xo89Y7wpYFcv9jiT+7FgJxw5p/r1ERAlf4kefsd79J//yN454VPeZ9DnP3zgkoSnhS/zI7+5dweljJfyTfLwAOhRBQR+/I5EEpoQv8aXP+V5v1mmh1HrOeZ9J37E64UoiooQv8aXPeXAkfIKRePZshEM7PhvyEjlFSvgSX/qe791/vMDfOOJJ3WdR99mInCIlfIkvnQd4C6l98p7fkcSPT/7lLZjWZbDfkUiCU8KX+GLmDV18oh5+vY8XeJ9Jmv67SmT0DZL402cs7PsE9m/1OxL/HaiAvZu8A7YiEVLCl/jTV/Px69X90umj8XuJnBK+xJ/uZ0NGnhI+eOckpOdC4XC/I5EkoIQv8ScQ9C7yoQO33h+93udCIN3vSCQJKOFLfOp7gbd+TOUevyPxz5F93meg4RxpJ0r4Ep/6XQg4+PhdvyPxz5b3AacDttJulPAlPvUaDek5sOkdvyPxz8fvQlo69Cr2OxJJEkr4Ep+CGd7a75v/6Xck/tn8T+g1CjJy/I5EkoQSvsSvfhfCzpVweLffkcTe0f2wbSn0/7zfkUgSUcKX+NX/Iu8+FXv5Hy8AF/rsMxBpB0r4Er96jvTmoG9OwXH8TfMhmOVNTxVpJxElfDN7yMxWm1mZmc0ys4Imyl1mZmvMbL2ZTY2kTkkhgXRvhkoqHrjdNN87hpGe5XckkkQi7eG/AZzlnBsOrAXuPbGAmQWA3wJfBoYC15vZ0AjrlVTR73Owew0c3OF3JLFzeDfsWKHhHGl3ESV859zrzrma8NP3gKJGipUA651zG51zx4DngEmR1CsppP/nvPtUGtapa6sO2Eo7a88x/MnAK41s7wVsafC8PLytUWZ2q5ktNrPFu3btasfwJCH1OAcyO6RWwt80HzLyvWMYIu0o2FIBM5sL9Gjkpfuccy+Fy9wH1ADPNLaLRrY1ecFS59x0YDpAcXGxLmya6gJB70pPqTRTZ9N8r82BFv97irRJi98o59ylzb1uZqXAFcAlzjV65elyoHeD50XAtrYEKSmu/0Ww9lXYXw4dGxs1TCL7t8Kn62H0zX5HIkko0lk6lwH3ABOdc5VNFFsEDDKz/maWAVwHzI6kXkkxA8Z59xve9jeOWKgfv9cBW2l/kY7hPwbkA2+Y2TIzexzAzHqa2RyA8EHd24HXgFXAC865lRHWK6mk2xDIL4QNb/kdSfRtmg/ZnaD7WX5HIkkookFC59zpTWzfBkxo8HwOMCeSuiSFmcHAi2HNHAjVQlrA74iiwzlY/yYM+IKuXytRoW+VJIaBF8ORvVCxzO9IomfHSji0HU5v9rCZyClTwpfEMOAL3n0yD+tseNO7H3ixv3FI0lLCl8SQ2wUKz0nuA7fr34Ruw6BDT78jkSSlhC+JY+DF3lWgqg76HUn7O3bYu37t6erdS/Qo4UviGHgxhGqS8ySszf+E2mMw8BK/I5EkpoQviaP3GO+yh8k4jr9+rte2Prp+rUSPEr4kjmCmdxWs9XP9jqT9rX/Ta5uWQ5YoUsKXxDJoPOzZCLvX+x1J+9mzCfZs0HRMiTolfEksZ3zJu1/7qr9xtKf66Zgav5foUsKXxFLQx5u6mEwJf+1r0Kk/nDbQ70gkySnhS+I540veFMYj+/yOJHJVh2DjPBg8wVtCQiSKlPAl8ZxxmTc9s24oJJFtfBtqq2Dwl/2ORFKAEr4knqJiyDnNGwpJdGtegawC6HOe35FIClDCl8STFvBm66x73Vs9M1GFar1jEYPGQyDd72gkBSjhS2I640ve6pnli/yO5NSVL4LKT2HwZX5HIilCCV8S08CLIS0dVr/sdySnbs0cSAtq/r3EjBK+JKasjjDg87BqtnfhkES05hXv7Nqsjn5HIilCCV8S15CJsHczbP/Q70jabvc62L0WztDsHIkdJXxJXGdeDpbm9fITzcq/efdDvuJvHJJSlPAlceV2gb4XwEcv+R1J2330N2/1z469/I5EUogSviS2oZO8oZGdq/2OpPV2r4MdK2DYlX5HIilGCV8S25lXePeJNKxTP5wz0d84JOUo4Uti61AIvc9LrGEdDeeIT5TwJfEN+6o3RLJzld+RtEzDOeKjiBK+mT1kZqvNrMzMZplZQRPlNpvZh2a2zMwWR1KnyEnOuhosAGUv+B1Jy1bO8u41nCM+iLSH/wZwlnNuOLAWuLeZsuOccyOcc8UR1ilyvLxuMHAcfPgXCOsL3gAAAAoNSURBVIX8jqZpzsHy57yZRRrOER9ElPCdc68752rCT98DiiIPSeQUnP112P8JbHnP70iaVr7Yu5ThOdf7HYmkqPYcw58MvNLEaw543cyWmNmt7ViniOfMyyE9J76HdZY/C8FsbyqpiA9aTPhmNtfMVjRym9SgzH1ADfBME7u5wDk3Cvgy8G0zu6iZ+m41s8VmtnjXrl1tbI6krMw8L+mvnAU1VX5Hc7Lqo7DiRRhyBWR18DsaSVEtJnzn3KXOubMaub0EYGalwBXADc41voqVc25b+H4nMAsoaaa+6c65YudccdeuXU+lTZKqhl8HR/d5i5LFm7WvwtH9Gs4RX0U6S+cy4B5gonOusokyuWaWX/cYGA+siKRekUYNHAcde8OSp/yO5GTLZ0J+IQz4gt+RSAqLdAz/MSAfeCM85fJxADPraWZzwmW6A/80s+XAQuBl59yrEdYrcrK0AIy80btO7N7NfkfzmQMVsO4NGP51L0YRnwQjebNz7vQmtm8DJoQfbwTOiaQekVYb+Q2Y9wv44Gm45Ed+R+P54E/gamH0TX5HIilOZ9pKcunYy7tG7NJnoLba72i8GJY85V3VqvMAv6ORFKeEL8lnVCkc2h4fB2/XzIGDFXDut/yOREQJX5LQoPHQsQ+8/7jfkcCiP3gHkgeN9zsSESV8SUKBIIyZAh+/C9uW+hfHrjWwaT4U36yDtRIXlPAlOY26ETLy4V+/8y+Gdx/1zqwdVepfDCINKOFLcsrq6CX9lX+FA9tiX//+rVD2PIz6pncpRpE4oIQvyWvMFG+FygWPxb7u934HLgRjvx37ukWaoIQvyatTPxh+LSx+Ag7uiF29lXtg8ZPeOv2d+sauXpEWKOFLcrvof0FtFSx4NHZ1vvtrqK6EC78XuzpFWkEJX5LbaQO9Xv6iP8aml39gmzcd9OxroPuw6Ncn0gZK+JL8LvoBhKrh7Z9Gv655v4RQLYz7j+jXJdJGSviS/E4bCCVT4IMZUFEWvXp2rfXqKL4ZOvePXj0ip0gJX1LD538A2Z3gtf/wZu60N+dgzl2Qkef9ohCJQ0r4khqyO8HF98Hmd7z58e3twz97Z9Ve+iPvouoicUgJX1LH6MnQ+zx45Z72PYB7+FPvl0OvYq8OkTilhC+pIy0NJj0G1Ufgv78DoVDk+3QOXvq2d/nCr/zaq0MkTunbKamlyyD44oPeNWYX/Dry/b0/Dda+Al/8CfQ4K/L9iUSREr6knjFTYOhX4c0HYf3cU9/P+rneUM4ZX/b2KRLnlPAl9ZjBxN9At2Hw/DehfHHb97FtKbxwE3QbClf/3tunSJxTwpfUlNUBvvEi5HWFp78KG//R+vd+vAD+NNGb+fM/nofM/KiFKdKelPAldeV3h5vmQEFv+K+r4Z1feWfJNiVUCwt+A3/6CuR1h8mvetfQFUkQSviS2jr2gptfgTOvgDcfgMcvhOXPw9EDn5WpOgRlf/Zee/0/4YzL4Ftzlewl4QT9DkDEd9kFcM1T8NHf4K2fwaxbAYMO4YR+YCvgoPMA+NqTMOxKjdlLQlLCFwEvgQ+7EoZMgi3ve2P6ezd7r3XuD/0uhD7na569JDQlfJGG0tKg71jvJpJkIu6umNlPzKzMzJaZ2etm1rOJcqVmti5801WdRURirD1+nz7knBvunBsB/B340YkFzKwzcD8wBigB7jezTu1Qt4iItFLECd8512A6A7lAY2vPfgl4wzm3xzm3F3gDuCzSukVEpPXaZQzfzH4GfBPYD4xrpEgvYEuD5+XhbSIiEiOt6uGb2VwzW9HIbRKAc+4+51xv4Bng9sZ20ci2Rq9CYWa3mtliM1u8a9eu1rZDRERa0KoevnPu0lbu71ngZbzx+obKgS80eF4E/KOJuqYD0wGKi4ujcGkiEZHU1B6zdAY1eDoRWN1IsdeA8WbWKXywdnx4m4iIxEh7jOH/wswGAyHgY+A2ADMrBm5zzn3LObfHzH4CLAq/50Hn3J52qFtERFrJXDQu6NxOzGwX3h+RU9EF2N2O4fglWdoBaku8Spa2JEs7ILK29HXOdW3shbhO+JEws8XOuWK/44hUsrQD1JZ4lSxtSZZ2QPTaooVBRERShBK+iEiKSOaEP93vANpJsrQD1JZ4lSxtSZZ2QJTakrRj+CIicrxk7uGLiEgDSvgiIikiqRJ+U2vzm+dRM1sffn2U37G2xMweMrPV4XhnmVlBg9fuDbdljZl9yc84W8PMrjGzlWYWCp+Q1/C1RGvLZeFY15vZVL/jaQsze8LMdprZigbbOpvZG+HrVLyRKMuWm1lvM3vbzFaFv1vfDW9PuPaYWZaZLTSz5eG2PBDe3t/M3g+35Xkzy4i4Mudc0tyADg0efwd4PPx4AvAK3iJu5wHv+x1rK9oyHgiGH/8S+GX48VBgOZAJ9Ac2AAG/422hLUOAwXjrJxU32J5QbQEC4RgHABnh2If6HVcb4r8IGAWsaLDt/wBTw4+n1n3P4v0GFAKjwo/zgbXh71PCtSecl/LCj9OB98N56gXguvD2x4F/j7SupOrhu6bX5p8EPO087wEFZlYY8wDbwDn3unOuJvz0PbwF58Bry3POuSrn3CZgPd5FZeKWc26Vc25NIy8lWltKgPXOuY3OuWPAc3htSAjOufnAiUuaTAL+FH78J+CrMQ3qFDnnKpxzH4QfHwRW4S25nnDtCeelQ+Gn6eGbAy4G/hLe3i5tSaqED97a/Ga2BbiBz66+lejr8U/G+4UCid+WhhKtLYkWb2t0d85VgJdEgW4+x9NmZtYPGInXM07I9phZwMyWATvxLhC1AdjXoNPXLt+1hEv4p7g2f6vX44+lltoSLnMfUIPXHkjgtjT2tka2+d6WZiRavEnPzPKAF4HvnfALP6E452qdd5nYIrxfkkMaKxZpPe1yxatYcqe2Nn850LvBa0XAtnYOrc1aakv4Yu9XAJe48EAeCdqWJsRlW5qRaPG2xg4zK3TOVYSHOXf6HVBrmVk6XrJ/xjn31/DmhG0PgHNun5n9A28Mv8DMguFefrt81xKuh9+cZtbmnw18Mzxb5zxgf93PvnhlZpcB9wATnXOVDV6aDVxnZplm1h8YBCz0I8Z2kGhtWQQMCs+eyACuw2tDIpsNlIYflwIv+RhLq5mZAX8EVjnnftXgpYRrj5l1rZuFZ2bZwKV4xyTeBr4WLtY+bfH7CHU7H+1+EVgBlAH/DfRqcBT8t3jjYh/SYKZIvN7wDmBuAZaFb483eO2+cFvWAF/2O9ZWtOVKvN5xFbADeC2B2zIBb0bIBuA+v+NpY+wzgQqgOvzvcQtwGvAmsC5839nvOFvZlgvxhjjKGvwfmZCI7QGGA0vDbVkB/Ci8fQBeB2g98GcgM9K6tLSCiEiKSKohHRERaZoSvohIilDCFxFJEUr4IiIpQglfRCRFKOGLiKQIJXwRkRTx/wG16XdQY48JKwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Iteration, E = 5 -15.9109375\n", "nL 501\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAEWCAYAAABFSLFOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deZgU1bn48e87G7OwCTMi27CjgiLLXATcE+KCRozGuMULQaPGmO3mXqM/rxqNZjHE3GvMjZpojIi4LyRBcQ2ooDLIIovIOjLsyDoss76/P07N2DNMz/TMVHf18n6ep56prjpV9VbT9NvnnKo6oqoYY4wxfkkLOgBjjDHJxRKLMcYYX1liMcYY4ytLLMYYY3xlicUYY4yvLLEYY4zxlSUWkxJE5BsislFEykRkRAyPe5WIvB6r4wVFRDaIyHif9vW4iNzjx75MMCyxJBnvP/gh7wt0q/eftH0z23zb226fiHwoIr2aKf9zEXkyglgeFpH/C3mdKSIHwiwbE8n5tcFU4CZVba+qi6JxABHpKyIqIhm1y1R1uqqeHY3jNRHHBhGpEJH8BssXe/H1jWAfR5xLIvAzwbUxjitFpMT7bL8sIl2CjimWLLEkp6+rantgODACuDVcQS/p/BW4DugM3AQc9imOucAZIa+LgM+B0xssA1jo0zHD6QMsj/Ix4sl64IraFyJyIpATXDipQ0SGAg8DVwPdgIPA/zW5UZKxxJLEVHUrMBuXYMIWA6qA9apao6oLVHWnTyHMAY4P+eV8GvA0kNdg2XxVrQQQkee8mtZeEZnr/SdFRMZ4y9Nrd+41by315tNE5BYRWSsiX4jIsyLSRUTaiUgZkA4sEZG1XnkVkYEh+6prfhGRM0WkVER+KiLbRWSLiHwnpGyOiPzO+0W6V0TeE5EcXCIF2OPVGMeKyGQReS9k23EissDbboGIjAtZ9y8R+YWIvC8i+0Xk9dr3SUSyReRJ79z2eNt2a+K9nwb8e8jrScAToQVE5HwRWeTVVDeKyM9DVh9xLt423xWRlV58K0RkZMg2w0VkqXduz4hIdsixLvBqTHtEZJ6IDAtZN0JEPvb2+QyQTRgiMkBE3vbeh50iMl1EOnvrpgGFwN+9mG9uZPuficgHtTUxEfmeiCwPjdUHVwF/V9W5qloG3A5cLCIdfDxGXLPEksTENWmdB6xpolgFsBh4VkSO8vP4qloKlOCSB7iayrvAvAbL5oZs9iowCDga+BiY7u3rA+AA8JWQslcCT3nzPwQuwtWQegC7gT+qarlXewM4SVUHRBj+MUAnoCdwDfDHkPdnKjAKGAd0AW4GaviyJtbZa3KbH7pDrznkn8ADQFfgfuCfItK1wTl9xzv/LOA/veWTvHh6e9veABxqIv4PgI4icryXjC8DGjZfHsAln87A+cD3ROQib90R5yIilwI/97bpCFwIfBGyv28B5wL9gGHAZO+8RwKPAdd7sT8MzPSSfhbwMi4RdgGeAy5p4rwE+BXu3/h47/34OYCqXo2rEX/di/m+Rrb/Le4z/98iMgj4JfBtVT2ili4ihV4iDDddGSbGocCS2hequtY75uAmziu5qKpN7nlpjwHbgWURlD0d96VXBXyzwbpJwGpvmhTAeWwAyoD9uNrIW7gvh3DlH/Kmm3HNUUd5y+8Ffhdmm58DT0YYz+PA73E/YrYDubgvxdplu4Ezwmzb2TuHTt7re4DHvPkOuC/GPt7rlcBXQ7btDlQCGd5rBQaGrG/4+nHgHm/+TNyXdkbI+u3AGC/mQ7gk1TDevt5+Q7ebDLznzV8NfNRgm/nAZG/+X8B/h6y7EXjNm5+CS8jDIvwMjAf+G/clfC7wBpDhxdc3zHb/A/y+iXOZDfyoiWN+O+T1fcBD3vyfgF80KL8K9yPgdGAzICHr5tX+W0RwrhcBixqeezPb9AV2eZ+ZW6Pwf/At4IYGyzYBZ/p9rHidrMbypcdx/wEj8TnuC+Op0IXeL9I7gZOB0cCdftcCInSRqnbAfUEeB+Q3VkhE8nC/xu9T9+vuDeBNL+ZxwJs+xDIX9+VxIrBOVQ8C74UsywE+9OJJF5Ffe81Z+3BfEoTE/xSuSaEdcDHwsaqWeOv6AC/V/prEfWlU49q4W+MLVa0KeX0QaO/Fkg2sbcU+e+BqcKFKcLWiWlsbOSa4X/SzgadFZLOI3Ccimc0cbxquBjSZBs1gACJysoi8IyI7RGQvLuE3+lnx9Kbp8w4Xex/gp6G/9r199fCmTep9+3oavkehMR8tIk+LyCbvM/JkMzEfQVU3AO/gEswfW7JthMpwNbpQHXE/9lKCJRaPqs7F/Yqp47XnviYiC0XkXRE5ziu7QVWX4po/Qp0DvKGqu1R1N+6LOtJk5TtVnYNLmFPDFEnD9T1UeeVvAYpxzSi5wGs+hDEXOAnX1PKut2w57ovlfGCBftkMcSUwEfdruxPuPz645g9UdQXuS+c86jeDAWwEzlPVziFTtqpuChPXQdw51jomwvPZibu4obEmteYeFb4Z9yUbqhD3a7ZJqlqpqnep6hBc0r+A+n0ojW1TguvEnwC82EiRp4CZQG9V7YSruUrt5o2U30jj592cjcC9Df5tclV1BrAF6CkiElK+sIl9/cqLbZiqdgS+HRJzuLjrEZEJwFhczeK3TZQr9Ppqwk1Xhdl0Oe4zX7uf/kA74LPmYksWllia9gjwA1UdhWvrbu7Kjp64/0S1Sqn/azQI/wN8TUSO6MBX1f245PF/ItLNa+9+G/flUQE09Ys4zetQrp3aNVZIVdcA24Af4SUW79fph96y0P6VDkA5rt0+F9f+3dBTuP6U03Ht8bUeAu4VkT4AIlIgIhObiH8xcKVXSzqX+levhaWqNbhm0/tFpIe3/Vjv/Hfgfmz0D7P5LGCwuEtRM0TkMmAI8I/mjisiZ4nIiV5/yT5cM191BCFfA3xFVQ80sq4DsEtVD4vIaFyyrtXYufwF+E8RGSXOwNr3uxl/Bm7wakgiInniLhzogGsKrAJ+6L0nF+Nq++F0wNUI9ohIT+C/GqzfRvj3H3EXQzwKXItrtv66l2iOoKqfq+urCTdND3OY6d5+T/NaBe4GXvT+v6UESyxhiLsMdxzwnIgsxnU4dm9us0aWBTrgjaruwDWD3B6myLdx/xmX4JLiVbiOacF9gYZzBa6voXZqqolkLlAAvB+y7F1cB3VoYnkCVyPZBKzA1ZwamoFr4ntb61+99r+4X9+vi8h+b9uTm4jpR8DXgT24c365ibIN/SfwCbAAV8v9DZDmNfPdC7zvNfnUuzdHVb/A1TR+ikueNwMXaGRX4R0DPI9LKitxV9w1ey+Rqq5V1eIwq28E7vberzuAZ0O2O+JcVPU5b9lTuGadl3Ed7s3FUAx8F3gQ16e2Bq9jX1UrcM2ak711l9F47arWXcBIYC/uQoiGZX+F65jfIyL/2XBj3I/FV1R1lvfvcQ3wlwYXULSJqi7HNStOx/XNdcC91ylD6jdtpjZxN479Q1VPEJGOwCpVDZtMRORxr/zz3usrcB1013uvHwb+5VX5jTEmJViNJQxV3Qes9y6xxKvCn9TMZrOBs0XkKK8D/GxvmTHGpAxLLB4RmYFr7z1W3M1x1+CaSK4RkSW4DrmJXtl/E5FS4FLgYRFZDqCqu4Bf4JpIFgB3e8uMMSZlWFOYMcYYX1mNxRhjjK8S6sml0ZKfn699+/YNOgxjjEkYCxcu3KmqBY2ts8QC9O3bl+LicFdkGmOMaUhEwj4hwZrCjDHG+MoSizHGGF9ZYjHGGOMrSyzGGGN8ZYnFGGOMryyxGGOM8VVCJRYReUzcGOTLwqw/U9x424u96Y5Yx2iMMakuoRILkY3y+K6qDvemu2MQk4kja3eUMXPGQ+yd+1DQoRiTshIqsTQ2yqMxodbtOAArXqb9u/dAecqMq2RMXEmoxBKhsSKyREReFZGh4QqJyHUiUiwixTt27IhlfCaKKqtreLTqPNIr98OicAP8GWOiKdkSy8dAH1U9CfgDTYwKqKqPqGqRqhYVFDT6uBuTgCqra1iiAzl8zCj48CGoiWT0XmOMn5IqsajqPlUt8+ZnAZneGNcmRVRU1QBQNvw62L0ePrNx1oyJtaRKLCJyjIiINz8ad35fBBuViaWqGje+UNWxF0D7Y2Dh48EGZEwKSqinG3ujPJ4J5HsjON4JZAKo6kPAN4HviUgVcAi4XG0ks5RSWe1qLJmZmTDiKnjv97B3E3TqGXBkxqSOhEosqnpFM+sfBB6MUTgmDtU2hWWkp8GIq+Hd38Hi6XDGzQFHZkzqSKqmMGNqm8Ky0tOgSz/odwZ8PA1qagKOzJjUYYnFJJVKr8aSmS5uwahJsPdzWPd2gFEZk1ossZikUlldgwikp3mJ5bgLILszLHkm2MCMSSGWWExSqaxRMtPS8C4OhIx2MPQi+PSfUHEg2OCMSRGWWExSqayq+bIZrNaJl0LlAVj1ajBBGZNiLLGYpFJZXUNmRoOPdeE46NADPnk+mKCMSTGWWExSqaxRMtMbfKzT0uDES2DNG3DQnmFqTLRZYjFJpbKqhsw0OXLFiZdCTRWseCX2QRmTYiyxmKTSaFMYwDHDIH8wLHsh9kEZk2IssZik0mhTGIAIHH8hlLwPB+zxccZEkyUWk1TcVWFhPtbHfx20Bj6zq8OMiSZLLCapVFY3crlxre4nQadCWPn32AZlTIqxxGKSSkV1De0a62MBrznsAlj7tg1bbEwUWWIxSaW8soascIkF3CNeqitg9RuxC8qYFGOJxSSViuoa92TjcArHQG6+NYcZE0WWWExSKa+soV1GevgCaelw3ARY/TpUHo5dYMakEEssJqlUVDfTFAauOayizF16bIzxnSUWk1TKK6vDd97X6nsaZGS7WosxxneWWExSiajGkpXrkoslFmOiwhKLSSrN9rHUGnwO7FoHX6yNflDGpBhLLCaplEdSYwEYON79tVqLMb5LqMQiIo+JyHYRWRZmvYjIAyKyRkSWisjIWMdogqOqVFQ1cYNkqC793EMpP5sd/cCMSTEJlViAx4Fzm1h/HjDIm64D/hSDmEycqKiuAYisxgIw6Gx3ZVh5WRSjMib1JFRiUdW5QFMjNU0EnlDnA6CziHSPTXQmaOVVLrFEVGMBl1iqK2D93ChGZUzqSajEEoGewMaQ16XeMpMCKlqaWArHQlYHWG3NYcb4KdkSS2OPtdVGC4pcJyLFIlK8Y8eOKIdlYqG2xhJxU1hGFgw40z03TBv9mBhjWiHZEksp0DvkdS9gc2MFVfURVS1S1aKCgoKYBGei68saSwSXG9caOB72bYKdn0UpKmNST7IllpnAv3tXh40B9qrqlqCDMrFRXlUNtKDGAtD/TPd33b/8DseYlJVQiUVEZgDzgWNFpFRErhGRG0TkBq/ILGAdsAb4M3BjQKGaALS4jwXgqL5ussRijG8ygg6gJVT1imbWK/D9GIVj4kyL+1hq9T8LPnkeqishPTMKkRmTWhKqxmJMU1rVxwKuOaxiP2z62PeYjElFllhM0mhVHwtAv9MBgXXv+B+UMSnIEotJGuWVrehjAcjtAj2GWz+LMT6xxGKSxqFKV2PJyWxhUxi4fpbSBVC+3+eojEk9llhM0jhY4RJLblZrEsuZUFMFG2xUSWPayhKLSRqHvRpLdmsSS++TISPH+lmM8YElFpM0amssrWoKy8yGPmOtn8UYH1hiMUnjUGU1melCZnorP9b9z4Qdn8K+Rp8CZIyJkCUWkzQOVVS3rrZSq/+Z7u+G9/wIx5iUZYnFJI1DFdXktKZ/pVa3E6BdJ9jwrn9BGZOCLLGYpHGwsprcrDY8pSgtHfqMsyvDjGkjSywmaRyqqCa7LU1hAH1PhV1rrZ/FmDawxGKSxqHKKnIy2/iR7nuq+2u1FmNazRKLSRqHKtrYFAZwzInWz2JMG1liMUnjoB9NYXX9LHZlmDGtZYnFJI3DldWte5xLQ31P8fpZbPBRY1rDEotJGocq23gfS63afpYS62cxpjUssZikcbC8jfex1DpmGLTraP0sxrSSJRaTFGpqlLKKKjpm+zDatvWzGNMmllhMUiirqEIVOmT7NGZ931PhizWwf6s/+zMmhVhiMUmh7HAVAO39qLFAyP0sVmsxpqUssZiksN9LLB38Six1/SyWWIxpqYRKLCJyroisEpE1InJLI+sni8gOEVnsTdcGEaeJvbLySgDat/MpsaSlQ+FYSyzGtELCJBYRSQf+CJwHDAGuEJEhjRR9RlWHe9NfYhqkCcy+uhqLT30s4PWzrIb92/zbpzEpIGESCzAaWKOq61S1AngamBhwTCZO1DaF+XJVWK0+49zfjR/4t09jUkAiJZaewMaQ16XesoYuEZGlIvK8iPQOtzMRuU5EikWkeMeOHX7HamLM9857cP0sGTnwuSUWY1oikRKLNLJMG7z+O9BXVYcBbwJ/C7czVX1EVYtUtaigoMDHME0Q9h92fSy+NoVlZEGvIvh8vn/7NCYFJFJiKQVCayC9gHqDZqjqF6pa7r38MzAqRrGZgJWVVyECuX480iVU4RjYshTKy/zdrzFJLJESywJgkIj0E5Es4HJgZmgBEeke8vJCYGUM4zMB2nOwkk45maSlNVaxbYPCMaDVsKnY3/0ak8QSJrGoahVwEzAblzCeVdXlInK3iFzoFfuhiCwXkSXAD4HJwURrYm3XgQq65Gb5v+Neo0HSrJ/FmBbwsacz+lR1FjCrwbI7QuZvBW6NdVwmeLsOVNAlLwqJJbsjdBtq/SzGtEDC1FiMacrugxUcFY3EAu5GyY0LoLoqOvs3JslYYjFJYdeBCrpGLbGMgcoDsO2T6OzfmCRjicUkPFWNbo2l9xj31/pZjImIJRaT8PaXV1FZrdHpvAfo1BM6F1o/izERssRiEt6usgqA6HTe1yoc62os2vCeXGNMQ5ZYTMLbuu8wAMd0yo7eQQrHQNk22L0+escwJklYYjEJb8veQ0C0E8tY99f6WYxpliUWk/C27HU1lu7RTCz5x0J2Z+tnMSYCllhMwtu69zCdcjLJzYri/b5padD7ZKuxGBMBSywm4W3Zezi6tZVahWNg52dwYGf0j2VMArPEYhLexl0H6dk5J/oHqu1n2fhh9I9lTAKzxGISWnWNsm7nAQYc3T76B+sxAtKzrJ/FmGZYYjEJbfOeQ1RU1TCgIC/6B8vMhh4jrZ/FmGZYYjEJbc0ONwDXgIIY1FjA9bNsXgwVB2NzPGMSkCUWk9BWbN4HwKCjO8TmgIVjoaYSNn8cm+MZk4BalVhEZLCIfENErheR67z5QX4HZ0xzFn2+h/4FeXTK9XGs+6b0Hu3+Jkg/y+TJkxGRI6YxY8YEHVrc27JlC1deeSXHHXcc6enpTJ48+Ygyjz/+eKPv7+HDh5vc91tvvcW4cePo0KED3bt352c/+xlVVV8Oy3D48GEmT57MsGHDyMzM5Mwzz2x0P3PmzGHUqFFkZ2fTv39/HnrooXrr9+/fz49//GP69OlDTk4O48aNY8GCBS1+L1oq4gv/ReR44AbgUqBb7WLvr3pltgHPAg+rqg0LbKJKVVm8cQ+nD86P3UFzu0DB8QnVzzJ+/HimTZtWb1lWVhSfqwZUVVWRnp6OiM9DRcdQeXk5+fn53HLLLTzyyCNhy+Xm5rJ27dp6y7Kzw1/+vnTpUiZMmMCtt97KE088waZNm7jhhhuorq5m6tSpAFRXV5Odnc1NN93ErFmz2LNnzxH7Wb9+PRMmTGDKlCk8+eSTvPfee9x4440UFBRwySWXAHDttdeydOlS/va3v9GrVy+efPJJxo8fz4oVK+jZs2dr3pbIqGqTEzAAeB6oBsqAV4E7gKuAc4HzvPk7vXVlXtnngP7N7T8eplGjRqlJPJ+U7tE+P/uHPrPg89geeOaPVH/ZS7W6KrbHbYVJkybp+eef32QZQB9++GH95je/qbm5udqvXz+dNm1avTKlpaV62WWXaefOnbVz5846YcIE/eyzz+rW33nnnTp06FD961//qv3799e0tDTdv3+/lpWV6dVXX615eXl69NFH6y9/+Us9//zzddKkSaqqetddd+nQoUOPiGncuHH6gx/8oO1vgE9CYw7117/+VfPy8lq0r1tvvVWHDx9eb9nMmTM1Oztb9+3bd0T573//+3rGGWccsfzmm2/WgQMH1lt2zTXX6JgxY1RV9eDBg5qenq4vv/xyvTIjR47U2267rUUxNwYo1jDfqZE0ha0A+uLGj++mquep6t2qOl1VX1PVV735u1T1PFxtZgrQ39vWmKh4Y8U2ROArxx0d2wMXjoXyfbA9eT7ed999NxMnTmTJkiVcdtllTJkyhZKSEgAOHjzIWWedRXZ2NnPmzGH+/Pl0796d8ePHc/DglxcxrF+/nqeeeornnnuOJUuWkJ2dzU9/+lPmzJnDSy+9xNtvv82SJUt4991367aZMmUKn376KR999FHdslWrVjFv3jyuueYaAKZPn0779u2bnKZPnx6jd+pIhw4dok+fPvTq1YsLLriARYsWNVm+vLz8iBpNTk4Ohw8fZuHChREfd/78+Zx99tn1lp1zzjkUFxdTWVlJVVVVXc2n4bHee++9iI/TKuEyTu0ETGyuTDS2jeVkNZbEU1lVrWN/+aZe+ef5sT/4rg2qd3ZU/fCR2B+7hSZNmqTp6emal5dXb7r55pvrygB6yy231L2urKzUnJyculrLo48+qgMHDtSampq6MlVVVdqlSxd95plnVNXVWDIyMnTr1q11Zfbv36+ZmZk6Y8aMumVlZWXauXPner/+zz//fL3++uvrXt98880a+n9y3759unr16ianxn7p+ylcjWXevHn6+OOP66JFi3Tu3Ll6ySWXaE5OTr3aXEOzZ89WEdEnnnhCKysrtbS0VE877TQF9KmnnjqifLgay6BBg/Suu+6qt2zOnDkK6ObNm1VVdezYsXrqqadqaWmpVlVV6bRp0zQtLU0HDx7csjegETRRY2m2j0VVX2lD0mr1tsY05fmFpWzee5ifXzg09gfvXAgdergO/NHfjf3xW+j0008/oo+gc+fO9V4PGzasbj4jI4OCggK2b98OwMKFC1m/fj0dOtS/8u7gwYP1+hZ69epFt27d6l6vXbuWyspKRo8eXbcsLy+PE044od5+vvvd7zJp0iR+//vfk5WVxbRp07j99tvr1nfo0OGIY/utffsvL1f/9re/fUQneDhjx45l7Nixda/HjRvH8OHD+cMf/sADDzzQ6DZnn302U6dO5aabbuI73/kO7dq14/bbb+fdd98lPT29RXE37MNy3/dfLp82bRpTpkyhV69epKenM3LkSK644go+/ji6VzW2+Kl9IvItVX02GsFEcOxzgf8F0oG/qOqvG6xvBzwBjAK+AC5T1Q2xjtNE17JNe/nFP1Ywum8XvjakW/Mb+E3E3c9SMt8N/BXnHdS5ubkMHDiwyTKZmfWvqhMRampqAKipqWH48OE8/fTTR2zXpUuXuvm8vPo3qTb8kgvn/PPPJzc3lxdeeIFOnTqxZ88errjiirr106dP5/rrr29yHw8//DBXXXVVk2Wasnjx4rr5jh07tno/6enpFBUVsXr16ibL/cd//Ac/+clP2LJlC0cddRQbNmzg1ltvpV+/fhEf65hjjmHr1q31lm3fvp2MjAy6du0KwIABA5gzZw4HDhxg3759dO/encsuu6xFx2mN1jwO9ikROUpVH/Y9miaISDrwR+BrQCmwQERmqmpoQ/c1wG5VHSgilwO/AS6LZZwmOmpqlHU7y/j7ki08MncdnXIyeeCKEcFddVQ4Fpa/CHs3uhpMEhs5ciQzZswgPz//iJpOUwYOHEhmZiYfffRR3RfZwYMHWbZsGQMGDKgrl5GRweTJk3nsscfo1KkTF198cb3jXHjhhZx88slNHiu0ptQazSXeSKkqS5cu5aSTTmq2rIjQo0cPAGbMmEHv3r0ZOXJkxMcaO3YsL7/8cr1lb7zxBkVFRUf8UMjLyyMvL4/du3cze/Zs7rvvvoiP0xqtSSyPA/8nIvmqem/DlSIyFrhPVU9ra3ANjAbWqOo67zhPAxOpf4HARODn3vzzwIMiIlr708ln332imPIq96uuqUOErlK00eUtKhfmhdZf02B/ocsblItgm4Yrw23TMI56+2viPJp6//YfrmL7/sNUVrsyZw/pxi8uOoFuHWPwRONwCr37QD7/IO4TS3l5+RG/bNPT0ykoKIho+6uuuoqpU6cyceJE7r77bgoLC9m4cSOvvPIKN9xwA4MGNX4LW/v27ZkyZQo/+9nPyM/Pp3v37txzzz3U1NS4HwRfrHXv396NXNtvC7/5zb9IE+H1n38dXv5+3X46oNRrCKv3WfHmP2liXVPbNbNu8fovAGXfumLStmax+LcXkpWRxpBeLvHd9exixgzKZ1D3juw7VMEDr37K0iXr+NPFR8MzVwPw4GurePD1VXx6/9frjvXbv6/g3JN6kCbw4oKN/PrF5Tz7o1NIf/bqumOvKN1LRVU1OxeupGzbfhb/8hwAhvc9CoAbuu/nwZK1/Pi8wVx/Vn/e/2wnjz/+MTO+PwamXwrA7KVbqalRjuvegTXbyvivpz/h2K5pfCfrVXhythtj6JuP4rcWJxZVvVZEdgC/8JLLTwBE5FjgV7gv90P+hglAT2BjyOtSoOHPmLoyqlolInuBrsARzzkXkeuA6wAKC1v3xbDvUGVdYnH7DNn/kcdrdF3DH9xC4zupt03odhK6jRy5v3r7kEaXNxQu1iP31/g2R64Lt6ap/X05n5OZzjGdcuiXn8u4Afn07pIbPvhY6TYUsjq4fpZh3wo6mia9+eabdO/evd6ynj17UlpaGtH2ubm5zJ07l1tuuYVLL72UvXv30qNHD8466yyOOuqoJredOnUqBw4c4MILL6R9+/b85Cc/Ydvna8le/yb84UWvlNA/tytn9M+hZHclZ+asgnWfUf8/wBGfxEZnj1zX1HbNrxtxc/2+iL8v3EifLllsuMv1Se3ZuYXr3lrJ1v2VdMpOZ0SvXOb+6DhGd94FO3cBws5tm1m1eZ9LpJ5Xi9dx70ufUF5Vw0k983jlhmM5b2A17Cmpi2PCr4sp2VVet82I214HQP90KpLhkCgAABbBSURBVAD9smHW94fwk+fW8ac319CjUxYPXDaAS4ZkQtl2EGHv7i+49aX1lO4pp0tuJpeMzOfei/qTWbkXKhu+b/6R1v6YF5EfA78DZuDuXZmCS/WPAner6tYmNm/N8S4FzlHVa73XVwOjVfUHIWWWe2VKvddrvTJfNLXvoqIiLS4u9jNckwqmXQz7t8CNiXEXfjwoX/k6fUafy3+dmc9Pb70DBn0NjuoHGVkMGTKEq666ittuuy3oME0ERGShqhY1tq4tQ+79GTgfuBKXUJ4Gbq9tqoqCUqB3yOtewOYwZUpFJAPoBOyKUjwm1RWOhXfugUO7IafpX+6patGiRaxcuZLRo0ezf1sJv7nxIvZXCJfd/xYMOhFwHc4zZsxgw4YNzXbSm8TQmqvCMoEbgVuBAuBjYDiQSf2mKr8tAAaJSD9gE3A5LqmFmglMAuYD3wTejlb/ijF1/SwbF8Dgs5sum8Luv/9+Vq1aRYZWMrygmrn/fI5eXlIB1/Gen5/Pww8/TH5+DB/PY6KmNTWW1bhawQrgGlX9p4h8C3eZ7ywR+YaqlvkZJNT1mdwEzMZdbvyYqi4XkbtxN+rMxDXDTRORNbiayuV+x2FMnZ6jIC3D9bNYYmnUiBEjKC4uhsN74XfHw9CLYPzF9crYb7/k05rEkg58F3hcVWsAVPVZEdkDvAC8IyLnqarvA4Or6ixgVoNld4TMH8Y9JNOY6MvKhe7DE+qBlIFZ9RpUHoBRk4OOxMRAax6bP0hVH6tNKrVU9XVgPNAPmOdHcMbEvcIxsGkhVJU3XzaVffYqtO8GPRvt6zVJpsWJxasVhFv3IXAaEOANBsbEUOEYqC53o0qaxtXUwNp3YODXIM3GFkwFvv8rqxuH5RS/92tMXOpde6OkXXIc1o6VcHgP9D016EhMjDSbWETkqy3dqapu9LYd35qgjEkY7Qug60DY+GHQkcSvEq9lvM+4YOMwMRNJjeU1EXlbRC7wntfVJBHJ9IYqnkODjnZjklLhGNeBX1PTfNlUVPI+dOwV94++Mf6J5KqwEcD9uHtEdorIG8BHwFrcJb0CdAEGAWOAr3ivZ+PubzEmuRWOhUVPwheroeDYoKOJL6quxtLvjLh/CrTxTyTjsSwDzvYeLnkj7llgtc+0Dr0AXYB9wIvAn1R1gc+xGhOfQvtZLLHUt2sdlG2zZrAUE/F9LKo6H5gvIl2BAcAQ3J33CuwAlgGLGl6GbEzS6zoAcvNdc5jdp1HfBm8I3D52PU8qac0NktuBDcAiYLE3LVLVTT7GZUziqB34y64MO1LJPJd08xt/tL5JTq253PhPuIc/fg24G3gF+FxENojIPSLSvsmtjUlGhWNh9wbY7+tDvRNfyTzXDGb9KymlNTdI3qSqp6lqJ2Ag7hEq9+JqMv8P+EhE7FGvJrUUeuOe2+NdvrTnc9j7ud2/koLadIOkqq5T1RdV9Q5VHQ2cC/QA7vQlOmMSRfdhkJFjiSVUidc0aB33KafFiUVEdorIWyLyOxG5WkROqL2/xXte2N3Ad0Uky+9gjYlb6ZnQq8j6WUKVvA/ZneDoIUFHYmKsNZ33c4CTgDNxlxgrUOGN3rgY2IN7VtgpwDv+hGlMAigcA+/eD+Vl0M66GimZ55oI05q9r9okmdaMeX8JgNdJP7zBdBXQziv6hois48urxn7lS8TGxKvCMaDVsKkY+p8ZdDTBKtvubhgdeXXQkZgAtLqPRVXLVPU9VX1QVa/1xj5uD1zmFXkDN6LkWcA9bQ/VmDjXazRImvWzQMjzwez+lVTUmqGJl+PuYVmCdx+Lqu4AUNVq76GVCtyoquu9bXr6F7IxcSq7I3Qbav0s4PpXMnOh+0lBR2IC0Jo+lkPAxbjx5hVARLbi7m3pCvQBXq9NKgB286RJGYVjYdF0qK6C9Nb890oSJfOg92h3UYNJOa25j6W2yesE4GpgKvAJLqlkAE8C3/YxRmMSR+EYNwTvtk+CjiQ4B3fBtuXWDJbCWvWTynse2ApvesrXiIxJZHUPpPwAeowINpagbPwQUEssKczGCTXGT516unFHajuvU1HJ+5CeBT1HBR2JCUhCJBYR6SIib4jIau9vo4+MEZFqEVnsTTNjHacxABSOcx34qs2XTUYl86BnEWRmBx2JCUhCJBbgFuAtVR0EvOW9bswhVR3uTRfGLjxjQvQ9BQ7sgJ2rg44k9srLYPNie4xLikuUxDIR+Js3/zfgogBjMaZptX0LJe8FG0cQSj9yN4laYklpiZJYuqnqFgDv79FhymWLSLGIfCAiTSYfEbnOK1u8Y8cOv+M1qaxLf2h/DGx4P+hIYq9kHki6u9TYpKy4udBeRN4Ejmlk1W0t2E2hqm4Wkf7A2yLyiaqubaygqj4CPAJQVFSUoo3hJipEXHNYyTzXz5JKY5FseN/dFNmuQ9CRmADFTWJR1fHh1onINhHprqpbRKQ7buyXxvax2fu7TkT+BYwAGk0sxkRVn3Gw7AXYvd7VYFJB5WH3nLTR1wUdiQlYojSFzQQmefOTcKNW1iMiR4lIO28+H/d05RUxi9CYUH28wa1SqTls00KorrD7V0zCJJZfA18TkdW4IZF/DSAiRSLyF6/M8UCxiCzBPa7/16pqicUEo+BYyO2aWvezlMwDxD19wKS0uGkKa4qqfgF8tZHlxcC13vw84MQYh2ZM40Rcc1gqXRlW8r57CGdul6AjMQFLlBqLMYmnzylu3Pc9G4OOJPqqK2HjR3aZsQEssRgTPXX3s6RAc9iWpe7hm5ZYDJZYjImebkPdmO+p0BxW4l2kUGiJxVhiMSZ60tLd+CypcGVYyTzoOhA6dAs6EhMHLLEYE019ToFda2H/1qAjiZ6aavh8njWDmTqWWIyJprp+liSutWxfAYf32v0rpo4lFmOiqftJkNUeNiRxP0vtxQlWYzEeSyzGRFN6hvslv25O0JFET8n70Km3G+DMGCyxGBN9/U53/Sx7S4OOxH+qrsZizWAmhCUWY6Kt/xnu7/q5wcYRDV+scYOaWTOYCWGJxZhoO3qoe25YMiaW2osSrMZiQlhiMSba0tKg72mun0WTbOifknmQdzR0HRB0JCaOWGIxJhb6nwH7N7umo2RS4t2/kkqDmZlmWWIxJhb61fazJNHVYbtLYO9GawYzR7DEYkwsdOkPHXsl12XHdf0rY4ONw8QdSyzGxIKIaw7b8C7U1AQdjT/W/Qty893FCcaEsMRiTKz0Ox0O7YZtnwQdSdvV1MDad2DAWe7iBGNC2CfCmFip7WdJhuaw7cvhwHYY8JWgIzFxyBKLMbHSsTvkD06ODvy1b7u//c8KNg4TlyyxGBNL/c9y47NUHg46krZZ+zYcPcQlS2MasMRiTCwNHA9Vh9z4JYmq4iCUzLdmMBNWQiQWEblURJaLSI2IFDVR7lwRWSUia0TklljGaExE+p4C6e1gzVtBR9J6n8+D6nLXcW9MIxIisQDLgIuBsA9bEpF04I/AecAQ4AoRGRKb8IyJUFaeu1N9zZtBR9J6a99xydHGtzdhJERiUdWVqrqqmWKjgTWquk5VK4CngYnRj86YFhr4VdjxaeI+Rn/t2+6myKzcoCMxcSohEkuEegIbQ16XessaJSLXiUixiBTv2LEj6sEZU2fgePc3EZvD9m91QxFb/4ppQtwkFhF5U0SWNTJFWuto7Cl4YR8lq6qPqGqRqhYVFBS0LmhjWqPgOOjYMzGbw2ovM7bEYpqQEXQAtVR1fBt3UQr0DnndC9jcxn0a4z8R1xy2/BWornLDFyeK1a9D+2PsMS6mSXFTY/HBAmCQiPQTkSzgcmBmwDEZ07iB46F8L2wqDjqSyFVXuua7wWfbY1xMkxLi0yEi3xCRUmAs8E8Rme0t7yEiswBUtQq4CZgNrASeVdXlQcVsTJP6nQGSDqvfCDqSyH0+H8r3weBzg47ExLmEqIOr6kvAS40s3wxMCHk9C5gVw9CMaZ2czlA4Bj6bDV+9PehoIvPZbHeZce0zz4wJIyFqLMYkpcHnuicd7/k86Egi89lr0O80aNc+6EhMnLPEYkxQjvUq26teCzaOSOxc44ZVHnRO0JGYBGCJxZig5A+EroNgVQK03q6e7f4OPjvYOExCsMRiTJCOPQ82vAeH9wYdSdM+ew0Kjoej+gYdiUkAlliMCdKxE6CmMr7vwj+8F0rmWW3FRMwSizFB6j0acrrAqleDjiS8Va9BTRUc9/WgIzEJwhKLMUFKS3dXh62e7W5AjEcrXoEOPaDnqKAjMQnCEosxQTvufNfctD7sqBDBKd/vnmk25EK7295EzD4pxgRt4HjIag8rXg46kiOtft0N6nX8hUFHYhKIJRZjgpaZ7a4OW/n3+GsOWzET8grcUwKMiZAlFmPiwZCL4NDu+GoOqzjonmV2/NddX5AxEbLEYkw8iMfmsLVvQeUBawYzLWaJxZh4EI/NYctecJdC9z016EhMgrHEYky8qG0OWzcn6EjcVWqrXoUTLoH0zKCjMQnGEosx8WLQ1yC7Myx9OuhIXKd91WE46fKgIzEJyBKLMfEio52rIaz8BxzeF2wsS5+BLgPspkjTKpZYjIknJ10BVYfc3e5B2VvqHow57DIQCS4Ok7AssRgTT3oVuZrCkgCbwz55DlAYdmlwMZiEZonFmHgiAsOvgJL3YHdJ7I+vCoufgl6joUv/2B/fJAVLLMbEm2GXAwKLpsX+2CXvw87PYNTk2B/bJA1LLMbEm869YfA5sPBvUFUR22MX/xWyO8HQb8T2uCapWGIxJh7927VwYDt8+o/YHbNsh7to4KQrISs3dsc1SSchEouIXCoiy0WkRkSKmii3QUQ+EZHFIlIcyxiN8dWAr0LnPrDg0dgdc/F0N5pl0Xdid0yTlBIisQDLgIuBSJ7Qd5aqDlfVsAnImLiXlgZFU1wn/vZPo3+86ioofhT6nAIFx0b/eCapJURiUdWVqroq6DiMiakRV0NGDsz/Q/SPteJl2PM5jP1+9I9lkl5CJJYWUOB1EVkoItc1VVBErhORYhEp3rFjR4zCM6YF8rrCyKthyTOwb3P0jqMK8x6AroNg8HnRO45JGXGTWETkTRFZ1sg0sQW7OUVVRwLnAd8XkdPDFVTVR1S1SFWLCgoK2hy/MVEx9ibQGpj/x+gdY/0c2LIExv3Ahh82vsgIOoBaqjreh31s9v5uF5GXgNFE1i9jTHw6qo97fljxX+G0n0JuF3/3rwpzfgt5R7tHuBjjg6T5eSIieSLSoXYeOBvX6W9MYjvtP6DyILx3v//7XvuWu0Dg9P90Y8IY44OESCwi8g0RKQXGAv8Ukdne8h4iMssr1g14T0SWAB8B/1TV14KJ2BgfHX28e3z9h4+4B0T6paYG3robOhXanfbGVwmRWFT1JVXtpartVLWbqp7jLd+sqhO8+XWqepI3DVXVe4ON2hgfnfX/AIV3funfPpe/6PpWzrrVPbLfGJ8kRGIxJuV1LoSTb3A3MZbMb/v+Du+D2bfBMSda34rxnSUWYxLFmbdAp97wj5+0/Rli79wLZdvggv+FtHR/4jPGY4nFmESRlQcTpsKOlTD3vtbvZ/278NEj8G/XQC8bIdL4zxKLMYnk2HNh+FUwdyqsm9Py7Q/shBeudWOtjL/L//iMwRKLMYlnwm8hfxA8/x34Ym3k21Uehucmw6HdcOnj0K59tCI0Kc4SizGJJisPLp/h5qd9I7JLkKsq4IVrYMO7MPFB12lvTJRYYjEmEeUPhCufc7WPv4yH0oXhyx74wiWgT/8B5/4Ghn0rdnGalGSJxZhE1WsUTHkNJB0e/RrMurl+09ihPfDBQ/BgEZR+BBf/GcbcEFy8JmXEzbPCjDGt0G0ofO99eOMON57KRw+7535ltPOayBQKx8H5U11ZY2LAEosxiS6nM1z4AJxxM3w6C7YuhepKd+XX4LOhx4igIzQpxhKLMcmiUy84uclhiIyJCetjMcYY4ytLLMYYY3xlicUYY4yvLLEYY4zxlSUWY4wxvrLEYowxxleWWIwxxvjKEosxxhhfiaoGHUPgRGQHUNLKzfOBnT6GE5RkOQ+wc4lXyXIuyXIe0LZz6aOqBY2tsMTSRiJSrKpFQcfRVslyHmDnEq+S5VyS5TwgeudiTWHGGGN8ZYnFGGOMryyxtN0jQQfgk2Q5D7BziVfJci7Jch4QpXOxPhZjjDG+shqLMcYYX1liMcYY4ytLLK0gIr8QkaUislhEXheRHt5yEZEHRGSNt35k0LE2R0R+KyKfevG+JCKdQ9bd6p3LKhE5J8g4IyEil4rIchGpEZGiBusS7VzO9WJdIyK3BB1PS4jIYyKyXUSWhSzrIiJviMhq7+9RQcYYKRHpLSLviMhK77P1I295wp2PiGSLyEcissQ7l7u85f1E5EPvXJ4Rkaw2H0xVbWrhBHQMmf8h8JA3PwF4FRBgDPBh0LFGcC5nAxne/G+A33jzQ4AlQDugH7AWSA863mbO5XjgWOBfQFHI8oQ6FyDdi7E/kOXFPiTouFoQ/+nASGBZyLL7gFu8+VtqP2fxPgHdgZHefAfgM+/zlHDn430vtffmM4EPve+pZ4HLveUPAd9r67GsxtIKqrov5GUeUHsFxETgCXU+ADqLSPeYB9gCqvq6qlZ5Lz8AennzE4GnVbVcVdcDa4DRQcQYKVVdqaqrGlmVaOcyGlijqutUtQJ4GncOCUFV5wK7GiyeCPzNm/8bcFFMg2olVd2iqh978/uBlUBPEvB8vO+lMu9lpjcp8BXgeW+5L+diiaWVROReEdkIXAXc4S3uCWwMKVbqLUsUU3A1Lkj8cwmVaOeSaPFGopuqbgH3ZQ0cHXA8LSYifYERuF/6CXk+IpIuIouB7cAbuJrxnpAfl7581iyxhCEib4rIskamiQCqepuq9gamAzfVbtbIrgK/nru5c/HK3AZU4c4HEvhcGtuskWWBn0sTEi3epCci7YEXgB83aLFIKKpararDcS0To3HNx0cUa+txMtq6g2SlquMjLPoU8E/gTly27x2yrhew2efQWqy5cxGRScAFwFfVa2glQc8ljLg8lyYkWryR2CYi3VV1i9c8vD3ogCIlIpm4pDJdVV/0Fifs+QCo6h4R+Reuj6WziGR4tRZfPmtWY2kFERkU8vJC4FNvfibw797VYWOAvbXV5XglIucCPwMuVNWDIatmApeLSDsR6QcMAj4KIkYfJNq5LAAGeVfrZAGX484hkc0EJnnzk4BXAowlYiIiwKPASlW9P2RVwp2PiBTUXvUpIjnAeFyf0TvAN71i/pxL0FcqJOKE+/WyDFgK/B3oqV9edfFHXLvlJ4RcmRSvE64jeyOw2JseCll3m3cuq4Dzgo41gnP5Bu7XfjmwDZidwOcyAXcF0lrgtqDjaWHsM4AtQKX373EN0BV4C1jt/e0SdJwRnsupuKahpSH/RyYk4vkAw4BF3rksA+7wlvfH/dBaAzwHtGvrseyRLsYYY3xlTWHGGGN8ZYnFGGOMryyxGGOM8ZUlFmOMMb6yxGKMMcZXlliMMcb4yhKLMcYYX1liMSbOiMi/eePjZItInjd2xglBx2VMpOwGSWPikIjcA2QDOUCpqv4q4JCMiZglFmPikPeMsAXAYWCcqlYHHJIxEbOmMGPiUxegPW7UwuyAYzGmRazGYkwcEpGZuJEj+wHdVfWmZjYxJm7YeCzGxBkR+XegSlWfEpF0YJ6IfEVV3w46NmMiYTUWY4wxvrI+FmOMMb6yxGKMMcZXlliMMcb4yhKLMcYYX1liMcYY4ytLLMYYY3xlicUYY4yv/j872OIvYych/gAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Iteration, E = 6 -15.936292613636363\n", "nL 501\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAEWCAYAAABFSLFOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3xV9f348dc7gwRCIEACZjAlKkNmRMCCVakDLFjH14EW626/Vn/WVqHWUatttVa/tdqqdSvgaB1YqbhF62CDILJXmEFmgOz374/PSbgJ2bnJueP9fDzu495zzuec+/5cwn3fzzjniKpijDHGBEuM3wEYY4yJLJZYjDHGBJUlFmOMMUFlicUYY0xQWWIxxhgTVJZYjDHGBJUlFhMVRORHIrJJRPJFZHALvu9EEXm3pd7PLyKyXkTGBOlYz4rIPcE4lvGHJZYI4/0HP+R9gW7z/pO2rWOfS7399onIVyKSVUf5u0TkxXrE8riI/C1gOV5EDtSwbnh96tcEDwDXq2pbVV3YHG8gIj1EREUkrnydqk5V1dOb4/1qiWO9iBSJSGqV9Yu8+HrU4xhH1CUcBDPBNTGOS0Rkg/e3/YaIdPQ7ppZkiSUy/VBV2wKDgMHAlJoKeknnGeAaIAW4HigIUhyzgZMDlnOAjcDoKusA5gfpPWvSHVjWzO8RStYBF5cviMjxQGv/wokeItIPeBy4DOgCHAT+VutOEcYSSwRT1W3ALFyCqbEYUAKsU9UyVZ2rqjuDFMInQJ+AX86jgJeApCrrvlDVYgARedVrae0Vkdnef1JEZLi3Prb84F731hLvdYyITBaRNSLynYi8IiIdRSRBRPKBWGCxiKzxyquI9A44VkX3i4h8X0RyReRmEdkhIltF5CcBZVuLyJ+9X6R7ReQzEWmNS6QAe7wW4wgRuVxEPgvYd6SIzPX2mysiIwO2fSwivxOR/4rIfhF5t/xzEpFEEXnRq9seb98utXz2LwA/DlieBDwfWEBExonIQq+luklE7grYfERdvH2uFpHlXnzfiMiQgH0GicgSr24vi0hiwHud7bWY9ojI5yIyIGDbYBFZ4B3zZSCRGojI0SLyofc57BSRqSKS4m17AegGvOXFfEs1+98qIl+Wt8RE5Kcisiww1iCYCLylqrNVNR+4HThXRJKD+B4hzRJLBBPXpXUWsLqWYkXAIuAVEekQzPdX1VxgAy55gGupfAp8XmXd7IDd/gNkA52BBcBU71hfAgeAUwPKXgJM817fAJyDayFlALuBR1W10Gu9AQxU1aPrGf5RQHsgE7gSeDTg83kAGAqMBDoCtwBlHG6JpXhdbl8EHtDrDnkbeBjoBDwIvC0inarU6Sde/VsBv/TWT/Li6ertex1wqJb4vwTaiUgfLxlfCFTtvjyASz4pwDjgpyJyjrftiLqIyAXAXd4+7YDxwHcBx/sf4EygJzAAuNyr9xDgaeBaL/bHgRle0m8FvIFLhB2BV4HzaqmXAH/A/Rv38T6PuwBU9TJci/iHXsz3V7P/n3B/878RkWzg98ClqnpEK11EunmJsKbHJTXE2A9YXL6gqmu89zymlnpFFlW1h7te2tPADmBpPcqOxn3plQDnV9k2CVjlPSb5UI/1QD6wH9ca+QD35VBT+ce8xy247qgO3vp7gT/XsM9dwIv1jOdZ4CHcj5gdQBvcl2L5ut3AyTXsm+LVob23fA/wtPc6GffF2N1bXg6cFrBvOlAMxHnLCvQO2F51+VngHu/193Ff2nEB23cAw72YD+GSVNV4e3jHDdzvcuAz7/VlwJwq+3wBXO69/hj4TcC2nwHveK+vwCXkAfX8GxgD/Ab3JXwm8B4Q58XXo4b9/g94qJa6zAJurOU9Lw1Yvh94zHv9d+B3VcqvwP0IGA1sASRg2+fl/xb1qOs5wMKqda9jnx7ALu9vZkoz/B/8ALiuyrrNwPeD/V6h+rAWy2HP4v4D1sdG3BfGtMCV3i/SO4ETgWHAncFuBdTTOaqajPuCPA5Ira6QiCThfo3fr+7X3XvA+17MI4H3gxDLbNyXx/HAWlU9CHwWsK418JUXT6yI/NHrztqH+5IgIP5puC6FBOBcYIGqbvC2dQdeL/81ifvSKMX1cTfGd6paErB8EGjrxZIIrGnEMTNwLbhAG3CtonLbqnlPcL/oZwEvicgWEblfROLreL8XcC2gy6nSDQYgIieKyEcikicie3EJv9q/FU9Xaq93TbF3B24O/LXvHSvDe2xW79vXU/UzCoy5s4i8JCKbvb+RF+uI+Qiquh74CJdgHm3IvvWUj2vRBWqH+7EXFSyxeFR1Nu5XTAWvP/cdEZkvIp+KyHFe2fWqugTX/RHoDOA9Vd2lqrtxX9T1TVZBp6qf4BLmAzUUicGNPZR45ScD83DdKG2Ad4IQxmxgIK6r5VNv3TLcF8s4YK4e7oa4BJiA+7XdHvcfH1z3B6r6De5L5ywqd4MBbALOUtWUgEeiqm6uIa6DuDqWO6qe9dmJm9xQXZdaXZcK34L7kg3UDfdrtlaqWqyqv1XVvrikfzaVx1Cq22cDbhB/LPBaNUWmATOArqraHtdylfLdqym/ierrXZdNwL1V/m3aqOp0YCuQKSISUL5bLcf6gxfbAFVtB1waEHNNcVciImOBEbiWxZ9qKdfNG6up6TGxhl2X4f7my4/TC0gAVtYVW6SwxFK7J4Cfq+pQXF93XTM7MnH/icrlUvnXqB/+D/iBiBwxgK+q+3HJ428i0sXr7/4Q9+VRBNT2izjGG1AufyRUV0hVVwPbgRvxEov36/Qrb13g+EoyUIjrt2+D6/+uahpuPGU0rj++3GPAvSLSHUBE0kRkQi3xLwIu8VpJZ1J59lqNVLUM1236oIhkePuP8Oqfh/ux0auG3WcCx4ibihonIhcCfYF/1/W+InKKiBzvjZfsw3XzldYj5CuBU1X1QDXbkoFdqlogIsNwybpcdXV5EviliAwVp3f5512HfwDXeS0kEZEkcRMHknFdgSXADd5nci6utV+TZFyLYI+IZAK/qrJ9OzV//oibDPEUcBWu2/qHXqI5gqpuVDdWU9Njag1vM9U77iivV+Bu4DXv/1tUsMRSA3HTcEcCr4rIItyAY3pdu1Wzztcb3qhqHq4b5PYailyK+8+4GJcUJ+IGpgX3BVqTi3FjDeWP2rpIZgNpwH8D1n2KG6AOTCzP41okm4FvcC2nqqbjuvg+1Mqz1/6C+/X9rojs9/Y9sZaYbgR+COzB1fmNWspW9Uvga2AurpV7HxDjdfPdC/zX6/KpdG6Oqn6Ha2ncjEuetwBna/1m4R0F/BOXVJbjZtzVeS6Rqq5R1Xk1bP4ZcLf3ed0BvBKw3xF1UdVXvXXTcN06b+AG3OuKYR5wNfAIbkxtNd7AvqoW4bo1L/e2XUj1ratyvwWGAHtxEyGqlv0DbmB+j4j8surOuB+Lb6rqTO/f40rgySoTKJpEVZfhuhWn4sbmknGfddSQyl2b0U3ciWP/VtX+ItIOWKGqNSYTEXnWK/9Pb/li3ADdtd7y48DHXpPfGGOigrVYaqCq+4B13hRLvCb8wDp2mwWcLiIdvAHw0711xhgTNSyxeERkOq6/91hxJ8ddiesiuVJEFuMG5CZ4ZU8QkVzgAuBxEVkGoKq7gN/hukjmAnd764wxJmpYV5gxxpigshaLMcaYoAqrK5c2l9TUVO3Ro4ffYRhjTNiYP3/+TlVNq26bJRagR48ezJtX04xMY4wxVYlIjVdIsK4wY4wxQWWJxRhjTFBZYjHGGBNUlliMMcYElSUWY4wxQWWJxRhjTFBZYjHGGBNUllhMRFmTl8/b0/7K3tl13TrHGNNcLLGYiLJ6Rz5ly/9N6zmP+B2KMVHLEouJKAXFpSwsy6ZV/mbYv63uHYwxQWeJxUSUQ0WlLCrzbsuea5fpMcYPllhMRCkoLmWZ9kBj4iF3rt/hGBOVLLGYiHKouIxCWqFdjofN8/0Ox5ioZInFRJRDxaUASNZQ2LwAykp9jsiY6GOJxUSUguJSWsfHIlknQPEB2LHc75CMiTqWWExEOVRUSutWsZCV41ZstgF8Y1qaJRYTUQ55LRY69oLWHW0A3xgfWGIxEeVQcSmJ8TEg4lotuTaAb0xLs8RiIkpBUSmJ8bFuITMH8r6Fgn3+BmVMlLHEYiJKQYnXFQaQNRRQ2LLA15iMiTaWWExEqRi8B8gc6p7tDHxjWpQlFhNRDhWXHe4Ka90BOmVbYjGmhVliMRGl/DyWClk5bsqxqn9BGRNlLLGYiHKoqJrEciAP9mz0LyhjoowlFhNRDhSVHB5jgcPjLHaipDEtxhKLiRiqSn5hCe0S4w6v7NwPYlu564YZY1qEJRYTMQ4WlaIKSQkBiSWuFRx1PGxZ5F9gxkQZSywmYhwoLAGgbWCLBSBjCGxdZFc6NqaFhFViEZEzRWSFiKwWkcnVbL9cRPJEZJH3uMqPOI0/9pcnloQqiSVzCBTlw85VPkRlTPSJq7tIaBCRWOBR4AdALjBXRGao6jdVir6sqte3eIDGd/kFLrEkV9diAXcGfufjWjgqY6JPOLVYhgGrVXWtqhYBLwETfI7JhJD8ihZLfOUNqdnQqq0N4BvTQsIpsWQCmwKWc711VZ0nIktE5J8i0rWmg4nINSIyT0Tm5eXlBTtW44P9XoslKSG28oaYWEgfZNcMM6aFhFNikWrWVT2d+i2gh6oOAN4HnqvpYKr6hKrmqGpOWlpaEMM0filvsSRXbbEAZAyCbV9DSVELR2VM9AmnxJILBLZAsoAtgQVU9TtVLfQW/wEMbaHYTAiocVYYuAH80iLYUXVIzhgTbOGUWOYC2SLSU0RaARcBMwILiEh6wOJ4wG54HkXKWyxHdIVB5QF8Y0yzCpvEoqolwPXALFzCeEVVl4nI3SIy3it2g4gsE5HFwA3A5f5Ea/yw71AxreJiSIirJrF06OFuVWwD+MY0u7CZbgygqjOBmVXW3RHwegowpaXjMqFh98EiOrSpZnwF3K2KMwbDloUtG5QxUShsWizG1GX3wWI6tGlVc4HMIbBjORQdbLmgjIlCllhMxNhzsKj2xJIxBLQUti1puaCMiUKWWEzE2H2wmA5JNXSFgesKAxtnMaaZWWIxEWP3gSJSamuxtEuH5HQbZzGmmVliMRFBVdlzqLjmwftyGUNsyrExzcwSi4kI+wpKKC3T2sdYADIHw3er4dCelgnMmChkicVEhF0H3KVaOibVkVjKT5Tcajf+Mqa5WGIxEWHb3gIAurRLrL2gDeAb0+wssZiIsH1fPRNLm47QoaeNsxjTjCyxmIiwrSKxJNRdOHMIbLaZYcY0F0ssJiJs21tAUqtYkhPrmBUGbpxlXy7k72j+wIyJQpZYTETYvq+ALu3r6AYrVz7OYuezGNMsLLGYiLB5zyEy2reuX+H0gSAxNoBvTDOxxGLCnqqybucBeqYm1W+HhLaQeqwN4BvTTCyxmLC360AR+wtK6FHfxALeAP4C0Kp3tzbGNJUlFhP21u08AECvhiSWjMFwcCfs3dRMURkTvSyxmLC3ekc+AL3SGthiARtnMaYZWGIxYW/J5r0kJ8bRrWOb+u/UpT/ExNvMMGOagSUWE/a+zt3LgKz2iEj9d4pLgKP6R/wA/uWXX46IHPEYPny436GFvK1bt3LJJZdw3HHHERsby+WXX35EmWeffbbaz7egoKDWY3/wwQeMHDmS5ORk0tPTufXWWykpKanY/s0333DKKafQpUsXEhMT6dWrF7/+9a8pKiqqdJyioiLuuOMOevbsSUJCAt26dePhhx+u2P7qq6+Sk5NDSkoKSUlJDBo0iOeee65pH0w9hNU9742pKr+whOVb93HN6F4N3zljMHz9Tygrg5jI/Y01ZswYXnjhhUrrWrWq42KdTVRSUkJsbGzDkn2IKSwsJDU1lcmTJ/PEE0/UWK5NmzasWbOm0rrExJrPqVqyZAljx45lypQpPP/882zevJnrrruO0tJSHnjgAcD9+0yaNInBgweTkpLC4sWLufrqqykpKeH++++vONbFF1/Mpk2beOKJJ8jOzmb79u0cOnSoYnunTp34zW9+w3HHHUd8fDz//ve/ufLKK0lLS2Ps2LGN/WjqpqpR/xg6dKia8PTO0q3a/dZ/6xdrdjZ85/nPq97ZTjVvZfADCxGTJk3ScePG1VoG0Mcff1zPP/98bdOmjfbs2VNfeOGFSmVyc3P1wgsv1JSUFE1JSdGxY8fqypWHP7c777xT+/Xrp88884z26tVLY2JidP/+/Zqfn6+XXXaZJiUlaefOnfX3v/+9jhs3TidNmqSqqr/97W+1X79+R8Q0cuRI/fnPf970DyBIAmMO9Mwzz2hSUlKDjjVlyhQdNGhQpXUzZszQxMRE3bdvX4373XTTTTp8+PCK5VmzZmm7du00Ly+vQe8/ePBgnTx5coP2qQ4wT2v4To3cn2kmKsxato3khDiGdOvQ8J1tAL/C3XffzYQJE1i8eDEXXnghV1xxBRs2bADg4MGDnHLKKSQmJvLJJ5/wxRdfkJ6ezpgxYzh48GDFMdatW8e0adN49dVXWbx4MYmJidx888188sknvP7663z44YcsXryYTz/9tGKfK664gm+//ZY5c+ZUrFuxYgWff/45V155JQBTp06lbdu2tT6mTp3aQp/UkQ4dOkT37t3Jysri7LPPZuHC2sftCgsLj2jRtG7dmoKCAubPn1/tPqtXr+add97h5JNPrlj3xhtvcMIJJ/Dggw+SlZVFdnY2N9xwA/n5+dUeQ1X54IMPWLFiBaNHj25gLRuopowTTQ9rsYSn3QcK9ZjbZuqvX1vSuAOUFKvec5TqzFuCG1gImTRpksbGxmpSUlKlxy23HK4zUOkXbHFxsbZu3bqi1fLUU09p7969taysrKJMSUmJduzYUV9++WVVdS2WuLg43bZtW0WZ/fv3a3x8vE6fPr1iXX5+vqakpFT69T9u3Di99tprK5ZvueUWDfw/uW/fPl21alWtj9p+6QdDTS2Wzz//XJ999llduHChzp49W8877zxt3bp1pdZcVbNmzVIR0eeff16Li4s1NzdXR40apYBOmzatUtkRI0ZoQkKCAnr11VdraWlpxbYzzjhDExISdOzYsfrll1/qO++8o9nZ2XreeedVOsaePXs0KSlJ4+LiNCEhQZ988smmfRgeammx2BiLCVt/+WAVRaVl/HhEj8YdIDbOXd4lwlsso0ePPmKMICUlpdLygAEDKl7HxcWRlpbGjh3uIp3z589n3bp1JCcnV9rn4MGDlcYWsrKy6NKlS8XymjVrKC4uZtiwYRXrkpKS6N+/f6XjXH311UyaNImHHnqIVq1a8cILL3D77bdXbE9OTj7ivYOtbdu2Fa8vvfRSHnvssXrtN2LECEaMGFGxPHLkSAYNGsRf//rXSoPogU4//XQeeOABrr/+en7yk5+QkJDA7bffzqeffkpsbGylsi+//DL79+9n8eLF/OpXv+K+++5jypQpAJSVlSEiTJs2jfbt2wPwyCOPcMYZZ7B9+/aKf4vk5GQWLVpEfn4+H3zwAb/4xS/o0aMHp512Wv0/oAYKq8QiImcCfwFigSdV9Y9VticAzwNDge+AC1V1fUvHaZrf6wtzeea/65l4YjeOPaoJXzoZQ2DeU1BaDLH1uDJyGGrTpg29e/eutUx8fOW6iwhlZWWA+wIbNGgQL7300hH7dezYseJ1UlLl84jcj1rqHMAfN24cbdq04V//+hft27dnz549XHzxxRXbp06dyrXXXlvrMR5//HEmTpxYa5naLFp0+I6i7dq1a/RxYmNjycnJYdWqVbWW+8UvfsFNN93E1q1b6dChA+vXr2fKlCn07NmzUrmuXbsC0LdvX0pLS7nqqqv41a9+RVxcHOnp6WRmZlYkFYA+ffoAsHHjxorEEhMTU/HvP2jQIJYvX87vf/97SywAIhILPAr8AMgF5orIDFX9JqDYlcBuVe0tIhcB9wEXtny0JtgKikvZvq+Ab7bs441Fm5m1bDvDenbk9rP7Nu3AmUPgy0dhx3JIH1B3+Sg0ZMgQpk+fTmpq6hEtndr07t2b+Ph45syZU/GFefDgQZYuXcrRRx9dUS4uLo7LL7+cp59+mvbt23PuuedWep/x48dz4okn1vpegS2lxqgr8daXqrJkyRIGDhxYZ1kRISMjA4Dp06fTtWtXhgwZUmP5srIySkpKKC0tJS4ujpNOOolXX32V/Pz8ihbXypUrAejevXutxyksLGxItRosbBILMAxYraprAUTkJWACEJhYJgB3ea//CTwiIqLlP52C7IyHZlNQUgoceckp5fCKI7bVEE3VMLXStgYcv9b3qu9+DYglYEUtb3fEttr200r7KQXFZRXL7RLjuPG0bP73lN60imvi/JPAS+hHaGIpLCxk27ZtldbFxsaSlpZWr/0nTpzIAw88wIQJE7j77rvp1q0bmzZt4s033+S6664jOzu72v3atm3LFVdcwa233kpqairp6encc889FV045M6DdZ/ArrVc1Wk19338ETExwru/OAGePqviOMkoldqklf4Aa/vP1fRtizblgyr7VqwmZlMci34zhFaxQt8M1zr77VsbGN4rmey01uwrKObhj7ayZFEef/9BCfz9ewA88vFWHvlkK9/eObjiuH96N5cz+6YQI8Jri7/jj//ZzCtXZhP7+PdAlRfm5JEYF8PxGa1pFSfM23iAKf/ayPkD25Pwj1EAXFJYwu8SSvjJqG7cdVY6ew6VcuPLGzh/UAqdX3af372ztnJi9zb06pRAYYky85u9vPDWFv56XhY8PBjadIKr3ifYwimxZAKBF3bKBar+jKkoo6olIrIX6ATsrHowEbkGuAagW7dujQpoQFZ7iksPf+FVbfJLjQsgASuq9hQELgZukyoHqbTtiN6Ghh+/YXHV3L3R+GNWPc7hFe0S4+jSLpFeaUkMyEohPjZIExo79oKE9i6xDJ0UnGOGmPfff5/09PRK6zIzM8nNza3X/m3atGH27NlMnjyZCy64gL1795KRkcEpp5xChw61z8Z74IEHOHDgAOPHj6dt27bcdNNNbN+0lsS1s+DJ11yh5HR6Jadz8rGd2PDdIb7fJ7W6P4YqR67tjz942wbf80mlEm8t2UX31Das/8vZAOwp3cY109axbW8B7dvEM7h7B2bfcRrDjkmtOObOkr2s2H4IOvSoOM5/Vq7h3neXU1hcxsDuKbz5y5M5a3Bmxfa49nH84a1vWLVtP6rQPS2J/z2rHzeN6wOtXLdlW+D9u7rz86e/4oQ/r6BDUgLnnNCDP04cCq3jQYT8+GJ++tp6cr87SOtWsRyXmcLzPz+Zi0d5LbSExnf71Uaa6cd80InIBcAZqnqVt3wZMExVfx5QZplXJtdbXuOV+a62Y+fk5Oi8efOaL3gT2p4bDwV74NrZfkcS8Qo3zKd7/2H86vsdufmuP0Gf8dDGjdP07duXiRMnctttt/kcpakPEZmvqjnVbQunFksu0DVgOQvYUkOZXBGJA9oDu1omPBO2MofA549AcQHE1/MulKZeFi5cyPLlyxk2bBj7d3/Hfdeeyf5C5cL7/wN93HfSjh07mD59OuvXr69zkN6Eh3BKLHOBbBHpCWwGLgIuqVJmBjAJ+AI4H/iwucZXTATJGAxlxbB9GWQN9TuaiPPggw+yYsUK4qSMQZ0Kmf3KI2T1OfxDt0uXLqSmpvL444+Tmppay5FMuAibxOKNmVwPzMJNN35aVZeJyN24E3VmAE8BL4jIalxL5SL/IjZhI8ObibNlgSWWIBs8eDDz5s1zg+N/Gw6tkuCcn1UqY7/9Ik/YJBYAVZ0JzKyy7o6A1wXABS0dlwlz7bMgKc0uod+cdnwDed/C2Q/5HYlpAXatMGNEXHdYhJ+B76uVs9zzMWfVXs5EBEssxoDrDtu5Agqrv4CfaaJ1s6FzP2iXXndZE/YssRgDrsWiZbBtid+RRJ6yMjd+1fUEvyMxLcQSizFQ+Qx8E1y71kDBXsi0iRHRwhKLMQDJXaBdpo2zNIfN3j1GMqs9l85EIEssxpTLGGwtluaweT7EJ0HasX5HYlqIJRZjymUMdt02h/b4HUlk2bwAMgZBTGzdZU1EsMRiTLnyWxVvXVR7OVN/JUWw7evDY1gmKlhiMaZc+iD3bOMswbNjGZQWHk7aJipYYjGmXJuO0KGnjbMEU3mSthlhUcUSizGBbAA/uLYsgNYdIaXmOxqayGOJxZhAmUNg7ybIz/M7ksiweaH7TOu4772JLJZYjAlkJ0oGT9EByFt++OrRJmpYYjEmUPpAQFwXjmmarYvdZXJsfCXqWGIxJlBCMqQeYy2WYKgYuLcWS7SxxGJMVZlD3Jei3YCqabYsgHZZ0Laz35GYFmaJxZiqMobAgR2wb4vfkYS3zQsg006MjEaWWIypqmIA38ZZGu3gLti9zsZXopQlFmOqOqo/xMTZOEtTlCdlmxEWlSyxGFNVfGvo3Mcu7dIUm72knDHI3ziMLyyxGFOdjCGuxWID+I2zZQF0yobE9n5HYnwQ15idROQYoB/QGVAgD1iqqquCGJsx/skYDAuec+MEHXv5HU342bwAep3sdxTGJ/VOLCLSB7gOuADoUr7ae1avzHbgFeBxVV0exDiNaVnl515sWWiJpaH2bYH8bTZwH8XqTCwicjRwH/Aj4BDwKfAFsAb4DpdcOgK9geHAVcDPReQ14FZVXds8oRvTjDr3hdgE98u7/3l+RxNeym9FbAP3Uas+LZZvgK+By4HXVPVAbYVFJAk4H7jB2zexiTEiIh2Bl4EewHrgf1R1dzXlSr1YATaq6vimvreJUrHxcNTxNjOsMTYvcLPqjjre70iMT+ozeP8/qpqjqi/UlVQAVPWAqj6nqkOBC5seIgCTgQ9UNRv4wFuuziFVHeQ9LKmYpskcAlsWQVmp35GEly0LXIsvvsm/KU2YqjOxqOqbjT14U/atYgLwnPf6OeCcIB3XmJplDoXiA5D3rd+RhA9V18qz8ZWo1uDpxiLyP80RSB26qOpWAO+5posPJYrIPBH5UkRqTT4ico1Xdl5ent17w1QjM8c9l48ZmLrtWgsFe+3Ck1GuMeexTBORa4MdiIi8LyJLq3lMaMBhumiRmX0AABaVSURBVKlqDnAJ8H/exINqqeoTXhdfTlpaWpPjNxGo09HuPIzceX5HEj5s4N7QuPNYngX+JiKpqnpv1Y0iMgK4X1VHNeSgqjqmpm0isl1E0lV1q4ikAztqOMYW73mtiHwMDMbNXjOm4URcl461WOpv8wKIaw1px/kdifFRg1ssqnoVcD/wOxF5qHy9iBzrTTH+DPeFHkwzgEne60nAEWM3ItJBRBK816nASbhZacY0XmYO7PjG3Q3R1G3LAncZl9hGnXttIkSjLumiqlOAXwA3iMiLIvIYbprvOOBx3DktwfRH4Acisgr4gbeMiOSIyJNemT7APBFZDHwE/FFVLbGYpsnKcXdB3LLI70hCX2kJbF1i3WCmcZd08fwDl0guwZ15/xJwe3OcEKmq3wGnVbN+Hu6ETFT1c8AmzpvgKp/dtHke9DjJ31hCXd5yKDlkA/emUbPC4kXkRtzYxanAAlxiiQc2BTc8Y3yWlAop3W2cpT4qBu7t5l7RrjFdYauAB3EXnhwfMAtrPDBTRNoGMT5j/JeVA7mWWOq0eQEkpti11UyjEksscDUwUFXfBlDVV3CJZTjwkTd4bkxkyBwK+3Jh/za/IwltWxa4bjCRusuaiNaYxJKtqk+ralngSlV9FxgD9AQ+D0ZwxoQEO1GybkUHYfs3NnBvgMZNNy6oZdtXwCiCcOFJY0JG+gB3UUU7UbJmWxeBlkLXYX5HYkJA0O8g6d2HxabPmMgR3xq69LMWS21y57rn8tadiWp1JhYROWKab11UdZO3b41n0xsTVjJz3MUVy8rqLhuNcue6QfukTn5HYkJAfVos74jIhyJytojE1lXYm478IxH5BJjZ9BCNCQGZQ6FwH+xc6XckoUcVNs2FrBP8jsSEiPqcIDkYN714BrBTRN4D5uDOY9nF4TtIZuNmhZ3qLc8CBjVDzMa0vKyAAfzOdh2sSvbmulsRW2IxnjoTi6ouBU73Li75M9y9US4u3xxQVIB9wGvA31V1bpBjNcY/nbIhoZ07A3/wRL+jCS3l4ytZNr5inHpf0kVVvwC+EJFOwNFAXyANl1zygKXAwqrTkI2JCDEx7oxymxl2pNx5EJcIXfr7HYkJEY2ZFbYDmA6cDSQA3wLvq+p8SyomomXlwPZlUHzI70hCS+5cl3Rj4/2OxISIxiSWvwNbcFcZvht3CfuNIrJeRO6xS7qYiJWZ487V2LrY70hCR0mh+zysG8wEaMwJkter6ihVbY+7PP4FwL24lsyvgTki0iG4YRoTAsq/PHNt+LDCtqVQWmgD96aSJp0gqaprVfU1Vb1DVYcBZwIZwJ1Bic6YUNK2M3ToAZu+8juS0FExcG9n3JvDGnPZ/J0i8oGI/FlELhOR/uXnt3jXC7sbuFpEWgU7WGN8lzUMNs1x524Yl1jaZUG7dL8jMSGkMTf6+gQYCHwfN8VYgSIRWQYsAvbgrhV2Eu5OjsZEjq7D4OtXYM9G6NDd72j8lzvXxlfMERqcWFT1PABvkH5QlcdE3EwxgPdEZC0u2SxU1T8EJWJj/NT1RPe8aY4llvwdsGcDDLvG70hMiGn0GIuq5qvqZ6r6iKpe5d3wqy1woVfkPdwdJU8B7ml6qMaEgM59IT4Jcuf4HYn/ys/psYF7U0WDWyxel9dCYDGuNbJIVfMAVLXUu2ilAj9T1XXePpnBC9kYH8XGQdZQG8AHl1xj4t1tBYwJ0JgxlkPAubjbESuAiGzDndvSCegOvFueVABUdXPTQzUmRHQ9ET59EIoOQKskv6PxT+48OOp4d1sBYwI05jyW8i6v/sBlwAPA17ikEge8CFwaxBiNCS1Zw9yJkpsX+B2Jf0pLXP2tG8xUozEtFrxLt3zjPaYFNSJjQl35LKhNX0HPUf7G4pcd30DxAUssplpBv4NkcxCRC0RkmYiUiUiNcxtF5EwRWSEiq0VkckvGaKJIm46Qekx0n4FfPsbU7UR/4zAhKSwSC+7KyecCs2sq4J2k+ShwFu7KyxeLSN+WCc9Ena7D3JdrtJ4oufFLSM6A9l39jsSEoLBILKq6XFVX1FFsGLDau8xMEfAS7t4xxgRf1xPh0G74brXfkfhj45fQbTiI+B2JCUFhkVjqKRN33ky5XG9dtUTkGhGZJyLz8vLymj04E2HKr421KQrPZ9mzCfblQrcRfkdiQlTIJBYReV9EllbzqG+ro7qfTjX2U6jqE6qao6o5aWlpjQvaRK/UYyCxfXSez7LxS/fcbbi/cZiQ1ahZYc1BVcc08RC5QGCHbxbu3Bpjgi8mxs2IisYB/I1fQKtk6NLP70hMiAqZFksQzAWyRaSnd2Xli4AZPsdkIlnXE2HHcjfWEk02fQVdT4CYWL8jMSEqLBKLiPxIRHKBEcDbIjLLW58hIjMBVLUEuB6YBSwHXlHVZX7FbKJA95GAwsYo6g47tMfdntnGV0wtQqYrrDaq+jrwejXrtwBjA5ZnAjNbMDQTzTKHumtlbfwcjj3T72haRu5cQG18xdQqLFosxoSk+NYuuWz43O9IWs7GLyAmztXbmBpYYjGmKbqPgC0Loeig35G0jI1fwlEDovvim6ZOlliMaYruJ0FZSXTMDispgs3zbXzF1MkSizFN0XUYIK6LKNJtXQwlBTa+YupkicWYpkhs7+5JsuG/fkfS/MqTpyUWUwdLLMY0VfeRsGmu6yqKZBu/gI69oG1nvyMxIc4SizFN1X0klBxyXUWRqqzUtcp6fM/vSEwYsMRiTFOVD2ZvjOBpx9uXQsFe6BGlNzYzDWKJxZimatsZOmVH9vks6z9zz9ZiMfVgicWYYOg+wo1BlJX5HUnzWP8ZdDwa2mX4HYkJA5ZYjAmG7ie5rqIdEXh5OhtfMQ1kicWYYOh+knsu7zKKJNu+tvEV0yCWWIwJhpSubiruutl+RxJ8FeMrJ/kbhwkblliMCZaeo92XcGmJ35EEl42vmAayxGJMsPQcDYX7Iut8lrJSN9utp3WDmfqzxGJMsPQY7Z7XfeJvHMG07WsotPEV0zCWWIwJlrZp0LlfZCWW8vGV7ja+YurPEosxwdRztLtnSUmh35EEx7rZ3vhKut+RmDBiicWYYOo52l1aPhLuz1JS5FosR5/idyQmzFhiMSaYuo8EiYmMace5c6H4APSyxGIaxhKLMcHUOgUyBsPaCBhnWfMhSKzNCDMNZonFmGDrORo2z4PC/X5H0jRrP4KsE9zNzIxpAEssxgRbr1OgrCS8L+9ycBdsXmDjK6ZRLLEYE2zdhkN8Eqx+3+9IGm/dJ4DC0af6HYkJQ2GRWETkAhFZJiJlIpJTS7n1IvK1iCwSkXktGaMxFeISXHfYqvdA1e9oGmfNR5DQHjKG+B2JCUNhkViApcC5QH2m2pyiqoNUtcYEZEyz630a7NkAu9b6HUnDqbrE0nMUxMb5HY0JQ2GRWFR1uaqu8DsOY+qt92nuORy7w3athb0bbXzFNFpYJJYGUOBdEZkvItfUVlBErhGReSIyLy8vr4XCM1GjYy/3WP2B35E03JoP3bONr5hGCpl2roi8DxxVzabbVPXNeh7mJFXdIiKdgfdE5FtVrbb7TFWfAJ4AyMnJCdOOcBPSeo+BhS9CcQHEJ/odTf2tehc69HSJ0ZhGCJkWi6qOUdX+1Tzqm1RQ1S3e8w7gdWBYc8VrTJ16j4Hig7DxC78jqb+iA+7kzmPP8jsSE8ZCJrE0lYgkiUhy+WvgdNygvzH+6PE9iG0VXuMsaz+B0kI45gy/IzFhLCwSi4j8SERygRHA2yIyy1ufISIzvWJdgM9EZDEwB3hbVd/xJ2JjgFZJ7tphK2f5HUn9rXwHEtpBt5F+R2LCWMiMsdRGVV/HdW1VXb8FGOu9XgsMbOHQjKndsWPhP7fAzlWQmu13NLVTdUnw6FMhrpXf0ZgwFhYtFmPCVvlYxYqZtZcLBVsXQf42OOZMvyMxYc4SizHNKaUbHHU8fBsGiWXlLEAg+wd+R2LCnCUWY5rbsWNh01eQH+LnS618B7oOg6RUvyMxYc4SizHN7dixgMKqEB7E37cVtiy02WAmKCyxGNPc0gdCu8zQ7g779t/u+dhx/sZhIoIlFmOam4gbxF/zIRQd9Dua6i17A9KOg87H+R2JiQCWWIxpCcedDSWHQvNkyfwdsOG/0HeC35GYCGGJxZiW0GMUtEmFZa/5HcmRls8AFPqe43ckJkJYYjGmJcTGQd/xbkpv0QG/o6ls2RvQKRs69/E7EhMhLLEY01L6nesuSrkyhK40lJ/nusH6nePGgowJAkssxrSU7iOhbRdYdsTVifzz7VugZTa+YoLKEosxLSUm1o1jrHoPCvf7HY3z9b+gU2/o0t/vSEwEscRiTEvqfy6UFMC3b/sdCezeABs+gwEXWTeYCSpLLMa0pKxh0KEHLJrmdySw5GX3PPBCf+MwEccSizEtKSYGBl4C62bDno3+xaEKi6e7adAp3fyLw0QkSyzGtLRBF7vnxS/5F8OmObBrLQy82L8YTMSyxGJMS0vpBj1Hw6KpUFbmTwyLp0N8G3dujTFBZonFGD8Mmgi718PGz1v+vQvzYem/oM94SEhu+fc3Ec8SizF+6PNDSGwP855u+ff++hUo3AcnXNny722igiUWY/zQqg0Mvgy+edPdC6WlqMKcJ+GoAZB1Qsu9r4kqlliM8UvOFVBWCvOfbbn33PgF7FgGJ1xl566YZmOJxRi/dDoaeo+B+c9ASVHLvOecf7guuOMvaJn3M1HJEosxfjrxWsjfDkv/2fzvtXu963obfJnrijOmmVhiMcZPvce463R99lDzTz3+78MgMTDif5v3fUzUC4vEIiJ/EpFvRWSJiLwuIik1lDtTRFaIyGoRmdzScRrTYCLwvZtg50p3peHmsn87LHzRnZzZLqP53scYwiSxAO8B/VV1ALASmFK1gIjEAo8CZwF9gYtFpG+LRmlMY/T7EXTsBZ/+2c3aag5fPAJlxXDS/2ue4xsTICwSi6q+q6ol3uKXQFY1xYYBq1V1raoWAS8BdpMJE/piYmHUzbB1MSxvhlbL3s0w5wnof76bMGBMMwuLxFLFFcB/qlmfCWwKWM711lVLRK4RkXkiMi8vLy/IIRrTQAMugrTj4P27oLQ4uMf++PfuZl6n/ia4xzWmBiGTWETkfRFZWs1jQkCZ24ASYGp1h6hmXY39Cqr6hKrmqGpOWlpa0ytgTFPExsEP7oZda2DeM8E77vZlsHAqDLsGOnQP3nGNqUWc3wGUU9UxtW0XkUnA2cBpqtV2ROcCXQOWs4AtwYvQmGaWfbq7OOWH97iLQyYf1bTjlZXB2790562Mujk4MRpTDyHTYqmNiJwJ3AqMV9WDNRSbC2SLSE8RaQVcBMxoqRiNaTIRGPeQu8Pkf25p+vHmP+MucnnGvdCmY9OPZ0w9hUViAR4BkoH3RGSRiDwGICIZIjITwBvcvx6YBSwHXlHVZX4FbEyjpPaGk29xJzIu/Vfjj7N7Pbx3J/T6vruSsjEtKGS6wmqjqr1rWL8FGBuwPBOY2VJxGdMsTroRVs6CGTfCUQNdsmmI4kPw8qXubpU/fNiuCWZaXLi0WIyJHrHxcMEz7vmli+HAd/Xft6wM3roRtn0N5/7DBuyNLyyxGBOK2mfBhS/Cno0w9Xw4tLvufcrKYObNsORlN7X4mDOaP05jqmGJxZhQ1eMkuOBZ1/p46nTIW1Fz2YK98PJEd+Owk/4fjPpli4VpTFWWWIwJZceeBT9+Aw7shMdGwXt3wO4Nh7cX7IW5T8Ejw9y4zJn3wZi7bFzF+CosBu+NiWo9vgf/+xXM+jV8/lf4718gqTPEJcDeXEAhcyhcNA2yhvodrTGWWIwJC207w3lPwim3wcp3YPtSd+mXDj3diZWZQ6yVYkKGJRZjwknHnjD8p35HYUytbIzFGGNMUFliMcYYE1SWWIwxxgSVJRZjjDFBZYnFGGNMUFliMcYYE1SWWIwxxgSVJRZjjDFBJdXf5Te6iEgesKHOgtVLBXYGMRy/REo9wOoSqiKlLpFSD2haXbqralp1GyyxNJGIzFPVHL/jaKpIqQdYXUJVpNQlUuoBzVcX6wozxhgTVJZYjDHGBJUllqZ7wu8AgiRS6gFWl1AVKXWJlHpAM9XFxliMMcYElbVYjDHGBJUlFmOMMUFliaURROR3IrJERBaJyLsikuGtFxF5WERWe9uH+B1rXUTkTyLyrRfv6yKSErBtileXFSJyhp9x1oeIXCAiy0SkTERyqmwLt7qc6cW6WkQm+x1PQ4jI0yKyQ0SWBqzrKCLvicgq77mDnzHWl4h0FZGPRGS597d1o7c+7OojIokiMkdEFnt1+a23vqeIfOXV5WURadXkN1NVezTwAbQLeH0D8Jj3eizwH0CA4cBXfsdaj7qcDsR5r+8D7vNe9wUWAwlAT2ANEOt3vHXUpQ9wLPAxkBOwPqzqAsR6MfYCWnmx9/U7rgbEPxoYAiwNWHc/MNl7Pbn87yzUH0A6MMR7nQys9P6ewq4+3vdSW+91PPCV9z31CnCRt/4x4KdNfS9rsTSCqu4LWEwCymdATACeV+dLIEVE0ls8wAZQ1XdVtcRb/BLI8l5PAF5S1UJVXQesBob5EWN9qepyVV1RzaZwq8swYLWqrlXVIuAlXB3CgqrOBnZVWT0BeM57/RxwTosG1UiqulVVF3iv9wPLgUzCsD7e91K+txjvPRQ4Ffintz4odbHE0kgicq+IbAImAnd4qzOBTQHFcr114eIKXIsLwr8ugcKtLuEWb310UdWt4L6sgc4+x9NgItIDGIz7pR+W9RGRWBFZBOwA3sO1jPcE/LgMyt+aJZYaiMj7IrK0mscEAFW9TVW7AlOB68t3q+ZQvs/nrqsuXpnbgBJcfSCM61LdbtWs870utQi3eCOeiLQF/gX8vyo9FmFFVUtVdRCuZ2IYrvv4iGJNfZ+4ph4gUqnqmHoWnQa8DdyJy/ZdA7ZlAVuCHFqD1VUXEZkEnA2cpl5HK2FalxqEZF1qEW7x1sd2EUlX1a1e9/AOvwOqLxGJxyWVqar6mrc6bOsDoKp7RORj3BhLiojEea2WoPytWYulEUQkO2BxPPCt93oG8GNvdthwYG95czlUiciZwK3AeFU9GLBpBnCRiCSISE8gG5jjR4xBEG51mQtke7N1WgEX4eoQzmYAk7zXk4A3fYyl3kREgKeA5ar6YMCmsKuPiKSVz/oUkdbAGNyY0UfA+V6x4NTF75kK4fjA/XpZCiwB3gIy9fCsi0dx/ZZfEzAzKVQfuIHsTcAi7/FYwLbbvLqsAM7yO9Z61OVHuF/7hcB2YFYY12UsbgbSGuA2v+NpYOzTga1AsffvcSXQCfgAWOU9d/Q7znrW5Xu4rqElAf9HxoZjfYABwEKvLkuBO7z1vXA/tFYDrwIJTX0vu6SLMcaYoLKuMGOMMUFlicUYY0xQWWIxxhgTVJZYjDHGBJUlFmOMMUFlicUYY0xQWWIxxhgTVJZYjAkxInKCd3+cRBFJ8u6d0d/vuIypLztB0pgQJCL3AIlAayBXVf/gc0jG1JslFmNCkHeNsLlAATBSVUt9DsmYerOuMGNCU0egLe6uhYk+x2JMg1iLxZgQJCIzcHeO7Amkq+r1dexiTMiw+7EYE2JE5MdAiapOE5FY4HMROVVVP/Q7NmPqw1osxhhjgsrGWIwxxgSVJRZjjDFBZYnFGGNMUFliMcYYE1SWWIwxxgSVJRZjjDFBZYnFGGNMUP1/6dnZQB84gu8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Iteration, E = 7 -15.948970170454544\n", "nL 501\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAEWCAYAAABFSLFOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXwV1fn48c+TBQIkECAkJOwCCogISFHwp61K3SsuVVS0UG21i/12V/haa7W1rdbq99tqFb+tdSmK2opSpSKuqCCb7JvsEhKSsO9LyPP748wNN8tNbpKbO3d53q/Xfd07M2dmnrmE+8w5Z2aOqCrGGGNMpKT4HYAxxpjEYonFGGNMRFliMcYYE1GWWIwxxkSUJRZjjDERZYnFGGNMRFliMUlBRK4SkS0isl9EhkRxv2NF5O1o7c8vIrJJREZFaFvPiMhvIrEt4w9LLAnG+w9+yPsB3eb9J82sZ52bvPX2ishcEelaT/lficg/wohlkoj8JWg6XUQOhJh3VjjH1wQPA3eoaqaqLmqOHYhITxFREUkLzFPVyap6YXPsr444NonIURHJqTZ/sRdfzzC2UeNY4kEkE1wT47hRRDZ7f9uviUgHv2OKJkssielrqpoJDAaGABNDFfSSzt+B24Bs4A7gcITimAV8OWh6GPAFcG61eQALI7TPUHoAK5p5H7FkI3BDYEJETgNa+RdO8hCRU4FJwM1AHnAQ+EudKyUYSywJTFW3ATNwCSZkMaAc2KiqFao6X1W3RyiED4H+QWfO5wBTgDbV5s1R1WMAIvKKV9PaIyKzvP+kiMhZ3vzUwMa95q2l3ucUEZkgIutFZIeIvCwiHUSkpYjsB1KBJSKy3iuvItInaFuVzS8i8hURKRSRn4pIqYgUi8g3g8q2EpE/emeke0TkYxFphUukALu9GuMIERkvIh8HrTtSROZ7680XkZFByz4QkV+LyCcisk9E3g58TyKSISL/8I5tt7duXh3f/fPAN4KmxwHPBRcQkctEZJFXU90iIr8KWlzjWLx1vi0iq7z4VorI0KB1BovIUu/YXhKRjKB9Xe7VmHaLyGwRGRS0bIiIfOZt8yUggxBEpLeIvOd9D9tFZLKIZHvLnge6A//2Yr6zlvXvEpFPAzUxEfmuiKwIjjUCxgL/VtVZqrofuAe4WkSyIriPmGaJJYGJa9K6BFhXR7GjwGLgZRFpH8n9q2ohsBmXPMDVVD4CZlebNytotf8AfYFc4DNgsretT4EDwPlBZW8EXvA+/xdwJa6GVADsAh5X1SNe7Q3gdFXtHWb4nYF2QBfgVuDxoO/nYeAMYCTQAbgTqOBETSzba3KbE7xBrznkTeBPQEfgEeBNEelY7Zi+6R1/C+Bn3vxxXjzdvHW/AxyqI/5PgbYi0t9LxmOA6s2XB3DJJxu4DPiuiFzpLatxLCJyLfArb522wBXAjqDtXQdcDPQCBgHjveMeCjwN3O7FPgmY5iX9FsBruETYAXgFuKaO4xLgd7h/4/7e9/ErAFW9GVcj/poX80O1rP8H3N/8L0SkL/Bb4CZVrVFLF5HuXiIM9boxRIynAksCE6q63tvnyXUcV2JRVXu556U9DZQCy8Moey7uR68c+Hq1ZeOAtd5rnA/HsQnYD+zD1Ubexf04hCr/pPe6E9cc1d6b/wDwxxDr/Ar4R5jxPAM8ijuJKQVa434UA/N2AV8OsW62dwztvOnfAE97n7NwP4w9vOlVwAVB6+YDx4A0b1qBPkHLq08/A/zG+/wV3I92WtDyUuAsL+ZDuCRVPd6e3naD1xsPfOx9vhmYV22dOcB47/MHwC+Cln0PeMv7fAsuIQ8K829gFPAL3I/wxcBMIM2Lr2eI9f4HeLSOY5kB/LCOfd4UNP0Q8KT3+Qng19XKr8GdBJwLFAEStGx24N8ijGO9ElhU/djrWacnsNP7m5nYDP8H3wW+U23eVuArkd5XrL6sxnLCM7j/gOH4AveD8ULwTO+M9F7gTGA4cG+kawFhulJVs3A/kP2AnNoKiUgb3Nn4Q+rO7mYC73gxjwTeiUAss3A/HqcBG1T1IPBx0LxWwFwvnlQR+b3XnLUX9yNBUPwv4JoUWgJXA5+p6mZvWQ9gauBsEvejcRzXxt0YO1S1PGj6IJDpxZIBrG/ENgtwNbhgm3G1ooBttewT3Bn9DGCKiBSJyEMikl7P/p7H1YDGU60ZDEBEzhSR90WkTET24BJ+rX8rnm7UfdyhYu8B/DT4bN/bVoH32qrer6+n+ncUHHOuiEwRka3e38g/6om5BlXdBLyPSzCPN2TdMO3H1eiCtcWd7CUFSyweVZ2FO4up5LXnviUiC0XkIxHp55XdpKpLcc0fwS4CZqrqTlXdhfuhDjdZRZyqfohLmA+HKJKC63so98pPABbgmlFaA29FIIxZwOm4ppaPvHkrcD8slwHz9UQzxI3AaNzZdjvcf3xwzR+o6krcj84lVG0GA9gCXKKq2UGvDFXdGiKug7hjDOgc5vFsx13cUFuTWn2PCi/C/cgG6447m62Tqh5T1ftUdQAu6V9O1T6U2tbZjOvEvxR4tZYiLwDTgG6q2g5Xc5XA6rWU30Ltx12fLcAD1f5tWqvqi0Ax0EVEJKh89zq29TsvtkGq2ha4KSjmUHFXISKXAiNwNYs/1FGuu9dXE+o1NsSqK3B/84HtnAS0BD6vL7ZEYYmlbk8BP1DVM3Bt3fVd2dEF958ooJCqZ6N++B/gqyJSowNfVffhksdfRCTPa+9+D/fjcRSo64w4xetQDrxa1lZIVdcBJcAP8RKLd3Y615sX3L+SBRzBtdu3xrV/V/cCrj/lXFx7fMCTwAMi0gNARDqJyOg64l8M3OjVki6m6tVrIalqBa7Z9BERKfDWH+EdfxnuZOOkEKtPB04WdylqmoiMAQYAb9S3XxE5T0RO8/pL9uKa+Y6HEfKtwPmqeqCWZVnATlU9LCLDcck6oLZj+SvwMxE5Q5w+ge+7Hv8HfMerIYmItBF34UAWrimwHPgv7zu5GlfbDyULVyPYLSJdgJ9XW15C6O8fcRdD/A34Fq7Z+mteoqlBVb9Q11cT6jU5xG4me9s9x2sVuB941fv/lhQssYQg7jLckcArIrIY1+GYX99qtczzdcAbVS3DNYPcE6LITbj/jEtwSXEsrmNacD+godyA62sIvOpqIpkFdAI+CZr3Ea6DOjixPIerkWwFVuJqTtW9iGvie0+rXr32v7iz77dFZJ+37pl1xPRD4GvAbtwxv1ZH2ep+BiwD5uNquQ8CKV4z3wPAJ16TT5V7c1R1B66m8VNc8rwTuFzDuwqvM/BPXFJZhbvirt57iVR1vaouCLH4e8D93vf1S+DloPVqHIuqvuLNewHXrPMarsO9vhgWAN8GHsP1qa3D69hX1aO4Zs3x3rIx1F67CrgPGArswV0IUb3s73Ad87tF5GfVV8adLL6uqtO9f49bgb9Wu4CiSVR1Ba5ZcTKuby4L910nDanatJncxN049oaqDhSRtsAaVQ2ZTETkGa/8P73pG3AddLd705OAD7wqvzHGJAWrsYSgqnuBjd4llnhV+NPrWW0GcKGItPc6wC/05hljTNKwxOIRkRdx7b2niLs57lZcE8mtIrIE1yE32iv7JREpBK4FJonICgBV3Qn8GtdEMh+435tnjDFJw5rCjDHGRJTVWIwxxkRUXD25tLnk5ORoz549/Q7DGGPixsKFC7eraqfallliAXr27MmCBaGuyDTGGFOdiIR8QoI1hRljjIkoSyzGGGMiyhKLMcaYiLLEYowxJqIssRhjjIkoSyzGJLDx48cjIjVeZ511Vv0rJ7ni4mJuvPFG+vXrR2pqKuPHj69R5plnnqn1+z18uMaAlFW8++67jBw5kqysLPLz87nrrrsoLy+vtezatWvJysoiMzOzxrLHH3+c/v3706pVK0455RSee67GsDv861//YsCAAbRs2ZIBAwYwderU8L6AJrDEYkyCGzVqFMXFxVVe06dPb9Z9lpeXE+9P9Thy5Ag5OTlMmDCBM88M/aDs1q1b1/h+MzIyQpZfunQpl156KRdeeCGLFi1iypQpTJs2jQkTJtQoe/ToUa6//nrOPffcGsueeOIJ7rrrLn75y1+yYsUK7rvvPr7//e/z73//u7LMnDlzGDNmDGPHjmXx4sWMHTuWa6+9lrlz5zbw22ggv4ewjIXXGWecoSYxbNtzSF97/n911weP+x1KTBg3bpxedtlldZYBdNKkSfr1r39dW7durb169dLnn3++SpnCwkIdM2aMZmdna3Z2tl566aX6+eefVy6/99579dRTT9W///3vetJJJ2lKSoru27dP9+/frzfffLO2adNGc3Nz9be//a1edtllOm7cOFVVve+++/TUU0+tEdPIkSP1Bz/4QdO/gAgJjjnY3//+d23Tpk2DtjVx4kQdPHhwlXnTpk3TjIwM3bt3b5X5P/rRj3T8+PG17mfEiBH6ox/9qMq8n/zkJ3r22WdXTl933XU6atSoKmUuuOACvf766xsUc22ABWpDE5tk8ObSYlLXvEmbhU/6HUpcuf/++xk9ejRLlixhzJgx3HLLLWze7O5/O3jwIOeddx4ZGRl8+OGHzJkzh/z8fEaNGsXBgwcrt7Fx40ZeeOEFXnnlFZYsWUJGRgY//elP+fDDD5k6dSrvvfceS5Ys4aOPPqpc55ZbbmH16tXMmzevct6aNWuYPXs2t956KwCTJ08mMzOzztfkyaHG3Gp+hw4dokePHnTt2pXLL7+cRYsW1Vn+yJEjNWo0rVq14vDhwyxcuLBy3ptvvskbb7zBn/70pwZtZ968eRw7dgxwNZYLL7ywSpmLLrqI2bNnh318jRIq4yTTy2osiaGiokJH/fED/cfvv6t6bzvVI/v9Dsl348aN09TUVG3Tpk2V15133llZBtAJEyZUTh87dkxbtWpVWWv529/+pn369NGKiorKMuXl5dqhQwd96aWXVNXVWNLS0nTbtm2VZfbt26fp6en64osvVs7bv3+/ZmdnVzn7v+yyy/T222+vnL7zzjs1+P/k3r17de3atXW+qp/pR1qoGsvs2bP1mWee0UWLFumsWbP0mmuu0VatWlWpzVU3Y8YMFRF97rnn9NixY1pYWKjnnHOOAvrCCy+oqmpRUZHm5+frnDlzVLX2mtHEiRM1NzdX582bpxUVFTp//nzNy8tTQIuKilRVNT09XZ999tkq6z377LPaokWLpnwdqlp3jcUe6WISxowV21hbup9eZw+DhZOhbA10Gep3WL4799xzeeqpp6rMy87OrjI9aNCgys9paWl06tSJ0tJSABYuXMjGjRvJysqqss7BgwdZv/7EwKFdu3YlLy+vcnr9+vUcO3aM4cNPjDTcpk0bBg4cWGU73/72txk3bhyPPvooLVq04Pnnn+eee04MeJqVlVVj35EW3DF+00038eST4dV4R4wYwYgRIyqnR44cyeDBg/nzn/8csqZx4YUX8vDDD3PHHXfwzW9+k5YtW3LPPffw0UcfkZqaWhnDd7/73TovsrjnnnvYtm0bI0eORFXJy8tj3LhxPPTQQ5XbARCpOrCtqtaYF2lxlVi8scn/F0gF/qqqv6+2vCVueNszcEO/jlHVTdGO00Tf5h0HuHvqcvp1zmL48N6wEChbbYkF17ncp0+fOsukp6dXmRYRKioqAKioqGDw4MFMmTKlxnodOpwYmbhNmzZVlrmT2po/bNVddtlltG7dmn/961+0a9eO3bt3c8MNN1Qunzx5Mrfffnud25g0aRJjx46ts0xdFi9eXPm5bdu2jd5Oamoqw4YNY+3atXWW+8lPfsKPf/xjiouLad++PZs2bWLixIn06tULgPfee48PP/yQ++67D3DfZUVFBWlpafzlL3/htttuo1WrVjz99NNMmjSJkpIS8vPzeeqpp8jKyiInJweAzp07s23btir7Li0trXIC0BziJrGISCrwOPBVoBCYLyLTVHVlULFbgV2q2kdErseNRT4m+tGaaNm25zDTlxXz5/fWUqHw+NihpHXIgNQWULrK7/ASwtChQ3nxxRfJycmpUdOpS58+fUhPT2fevHmVP5gHDx5k+fLl9O7du7JcWloa48eP5+mnn6Zdu3ZcffXVVfZzxRVX1HlVFtDkH8r6Em+4VJWlS5dy+un1DTbrEm5BQQEAL774It26dWPoUHcitGzZsiplX3/9dR544AHmzZtHly5dqixLT0+na9euAEyZMoXLL7+clBTXfT5ixAhmzpzJz3/+88ryM2fOZOTIkY0/yDDETWIBhgPrVHUDgIhMwY3oGJxYRgO/8j7/E3hMREQDp04R9sMpizha7s7qgvegVN1d1WWEXBa8tHrEwZPBh1PX9mouq/1rqLmvOuIIcZx1xVs9kFDr1RVv9WUHjxxn297D7DnkOimH9+rAg9cMoleOd9acc7KrsRiOHDlS46w1NTWVTp1qfeJ5DWPHjuXhhx9m9OjR3H///XTv3p0tW7bw+uuv853vfIe+ffvWul5mZia33HILd911Fzk5OeTn5/Ob3/yGiooKV4upOA7blkLxUr41sJwHH/yAFBHe/vW1MP3ED2EWkFXjbzf4D0dhdR3LQq1XY3nNZYs3lAGwd+NnpJS2ZPGjY2iRnsKArq6mdt9L8znr5Dz65rdj76Gj/OmNZSxd8jlPjOkBU78LwGPTl/HYf1aw+s9jKvf3h9eWcPGQbqSI8urcTfz+n4t4+ScXkPradwAYWC2eBcWfk1JxjIGrH/WOVfm8aA9z15Zx1smd2LX/CI+8sYLlS4t59rrL4ZXxAPzwtH2ce+87/O6GoVw1vDtT533B++8t5uP7LoKXvwEZ7eCKPxNp8ZRYugBbgqYLgeqnMZVlVLVcRPYAHYHt1TcmIrcBtwF07969UQFtKDvAkfLjJ7bJiSp/XbX/6k0DUmVZ7Z/r236VyaCF1cOosv1GxFQ9jhAfkeD1pNo6SMjjlGoRh/oeczJbMrxXB3p0bM3I3jkMKKjWfNGpH2yZV/vKSeadd94hPz+/yrwuXbpQWFgY1vqtW7dm1qxZTJgwgWuvvZY9e/ZQUFDAeeedR/v27etc9+GHH+bAgQNcccUVZGZm8uMf/5iSkhIy9m6ERwbAfpfwTgK+3LMFm3cf5yvMhqVSyz9+9T/E6n9xjVlWbXm1ZUN+tqbK9L/nb6JH+3Q2/aIPIOwu2sZtMxazbW857VqlMKRLBrO+34Ph6Wthk2sO2/5FCWu27obNn1Tu7z9zNvDAK/M5clw5vaAVr9/SnUsKSqCwpPZ4du0ArYCiE1ecHS85zCNTN7Gm9DDpqcJ5fbOY/cM+9KzYBNvcuiOzYcr4XvzijVXc+8pieue05KXxJ3Fm9k4o2wmtTjRlRlLcDE0sItcCF6nqt7zpm4HhqvqDoDIrvDKF3vR6r8yOurY9bNgwtfFYEsysP8B7v4GJhdCyeTt+TfiOHD5Mj/wO/PxLFfz0hgtg6Deg+5mQVcCAQYMZO3Ysd999t99hmjCIyEJVHVbbsniqsRQC3YKmuwJFIcoUikga0A7YGZ3wTEzp1N+9l62BrrX+7ZsoWLRoEatWrWL48OHs27ePB++8jX0HDjHm1p/BdQ+BCKWlpbz4+BNs2rSp3k56Ex/iKbHMB/qKSC9gK3A9cGO1MtOAccAc4OvAe83Vv2JiXK6XWEpXWWLx2SOPPMKaNWtIS0tjcIdDzLr7HLp6SQVcx3tOTg6TJk2qvJrJxLe4SSxen8kdwAzc5cZPq+oKEbkfd6PONOBvwPMisg5XU7nev4iNr9r3hLQM68D32ZAhQ04M+z3v/2D6z+DWR6r0Zdi5X+KJm8QCoKrTgenV5v0y6PNh4Npox2ViUEqquzLMLjmOHcv+CXmnWQ0yCdizwkziyu1vNZZYcWA7bJkL/S7zOxITBZZYTOLq1A/2boXDe/yOxGz6GFDo+1W/IzFRYInFJK7coCvDjL8K50NqS+g8qP6yJu5ZYjGJq1M/9166su5ypvkVLoCCwZDWwu9ITBRYYjGJK7sHpLeGUutn8ZUqlCyH/Pqfn2USgyUWk7hSUqDTKVBmV4b5as8WOLr/RNOkSXiWWExi69Tfaix+C3z/nSyxJAtLLCax5fZzDzo8tMvvSJJXoI8rt5+/cZioscRiElvgLNlqLf4pWw1Z+dCq7ichm8RhicUktspnhtmVYb4pXWn9K0nGEotJbO26QossuwPfLxXHoexz619JMpZYTGITcVeG2TPD/LFrE5QfshpLkrHEYhJfbj+rsfgl8L1bYkkqllhM4ssdAAfK3IMQTXQF+rY6neJvHCaqLLGYxBc86JeJrtLV0K67DQ+dZCyxmMSXO8C9W2KJvrLV1gyWhCyxmMSXmQcZ2XbJcbQdL4ftn9uNkUnIEotJfCKu1mI1lujauQGOHz1RYzRJwxKLSQ65/V1isfHVo6ey495qLMnGEotJDrn94cge2FvkdyTJo2w1IHZFWBKyxGKSg3XgR1/pSujQC9Jb+R2JiTJLLCY52DPDoq90tT3KJUlZYjHJoXUHyOxsNZZoKT8CO9bZpcZJyhKLSR65/a3GEi071oEet8SSpCyxmOSROwDK1rgn7prmFagZWmJJSpZYTPLI7e+etLtrk9+RJL7SVSCp0LGP35EYH1hiMcnDrgyLntKVkNMX0lr6HYnxgSUWkzwC91NYYml+NmpkUrPEYpJHy0zI7mEd+M3tyH7X3Jh7qt+RGJ9YYjHJxZ4Z1vwCg3vl2TPCklVcJBYR6SAiM0VkrffePkS54yKy2HtNi3acJg7k9ocda6H8qN+RJK5AjdCawpJWXCQWYALwrqr2Bd71pmtzSFUHe68roheeiRu5A6Ci3N1nYZpHyUpIbw3ZPf2OxPgkXhLLaOBZ7/OzwJU+xmLimT3apfmVrnRPNE6Jl58XE2nx8i+fp6rFAN57bohyGSKyQEQ+FZE6k4+I3OaVXVBWVhbpeE2syunr7q+wfpbmU7rS+leSXJrfAQSIyDtA51oW3d2AzXRX1SIROQl4T0SWqer62gqq6lPAUwDDhg2zQTqSRVpLd9OeJZbmsb8MDpTZFWFJLmYSi6qOCrVMREpEJF9Vi0UkHygNsY0i732DiHwADAFqTSwmieX2g+KlfkeRmKzj3hA/TWHTgHHe53HA69ULiEh7EWnpfc4BzgasId3UlDvA3Wdx9IDfkSSeQGLJsxpLMouXxPJ74Ksishb4qjeNiAwTkb96ZfoDC0RkCfA+8HtVtcRiasrtD6h7IKWJrNKV0LojtOnkdyTGRzHTFFYXVd0BXFDL/AXAt7zPs4HTohyaiUfBzwzrMtTfWBJNyUr3/Yr4HYnxUbzUWIyJnPa9ILUllFkHfkRVVLhkbc1gSc8Si0k+qWnQ6WS7MizS9nwBxw5Yx72xxGKSlD0zLPJKAleEWY0l2VliMckptz/s3QqHdvsdSeIoXeHec/v5G4fxnSUWk5wCHfiBJ/GapitdBdndoWWW35EYn1liMcnJnhkWeSUrrRnMAJZYTLJq1w1aZFo/S6SUH3XDEVjHvcESi0lWIu5HsMRqLBGxY60bjsAuNTZYYjHJLHeA63BWewZpk5XYM8LMCZZYTPLKGwiHdsG+Yr8jiX+lKyElDTr29TsSEwMssZjkFWi2KVnhbxyJoHQl5JwMaS38jsTEAEssJnkFBqMqWe5vHImgdKU1g5lKllhM8mrVHtp2tRpLUx3ZB7u/OHFvkEl6llhMcssbYFeGNVXgkm1LLMZjicUkt7xTYfsadx+GaZzKwb0ssRjHEotJbnkD3f0X2z/3O5L4VbLS3WzarrvfkZgYYYnFJDe7MqzpSldCp36QYj8nxrG/BJPcOvaB1BZ2ZVhjqbqkbM1gJoglFpPcUtMh5xR7GGVj7S+FQzut495UYYnFmLxTrSmssQI1PXtGmAliicWYvFPdY10O7PA7kvhTmVgG+huHiSmWWIwJnG2XWq2lwbYth7ZdoHUHvyMxMcQSizGBs21rDmu4kuXWDGZqsMRiTGYutM6xxNJQ5Ufc/T/WDGaqscRijIj3aBdLLA1SttrdXNrZEoupyhKLMeDOuktXQcVxvyOJH9sCHfen+RuHiTmWWIwB109Qfgh2bvQ7kvhRshzSWkHH3n5HYmKMJRZjIOjRLnYHfthKlrsxWFJS/Y7ExBhLLMaAe9aVpNgd+OFSdU1h1r9iahEXiUVErhWRFSJSISLD6ih3sYisEZF1IjIhmjGaOJfeyj03zDrww7Ov2D3KxfpXTC3iIrEAy4GrgVmhCohIKvA4cAkwALhBROwBRiZ8uQOsKSxcgY57q7GYWsRFYlHVVaq6pp5iw4F1qrpBVY8CU4DRzR+dSRh5A2HXJji81+9IYl/JMvduN0eaWqQ1ZiURORk4FcgFFCgDlqvq2gjG1lBdgC1B04XAmT7FYuJRZ69Zp2QF9BjhbyyxbttyyO4OGe38jsTEoLATi4j0B74DXAvkBWZ77+qVKQFeBiap6qqGBCIi7wCda1l0t6q+Hs4mapmndezvNuA2gO7dbeQ7A+QPcu/bllliqU/JcutfMSHVm1hEpDfwIHAVcAj4CJgDrAd24H7QOwB9gLOAbwE/EJFXgbtUdUM4gajqqMYcQJBCoFvQdFegqI79PQU8BTBs2LCQCcgkkax8aN0Rti31O5LYduwQ7FgHp17ldyQmRoVTY1kJLAPGA6+q6oG6CotIG+DrwH9562Y0McZwzQf6ikgvYCtwPXBjlPZtEoGIaw6zxFK30pWgFfaMMBNSOJ3316nqMFV9vr6kAqCqB1T1WVU9AxjT9BBBRK4SkUJgBPCmiMzw5heIyHRvv+XAHcAMYBXwsqrataOmYToPco92OX7M70hil10RZupRb40lzP6NiK9bbTtTgam1zC8CLg2ang5Mj8Q+TZLqPAiOH/We2mtXPNWqZDm0yITsnn5HYmJUgy83FpHrmiMQY2JC4MqwYmsOC2mbNwZLSlzcrWB80Ji/jBdE5PaIR2JMLMjp6x6suG2Z35HEJlV3Obb1r5g6NCaxPAP8RUTurm2hiIwQkY+aFJUxfklJdWOzWAd+7XZ/AUf2WP+KqVODE4uqfgt4CPi1iDwamC8ip3iXGH8MDIlciMZEWeDKMLWr0Fd9G5kAABViSURBVGsosTFYTP0adee9qk70bob8o4h0AvYDt+BuSJwE3B+5EI2Jss6DYOEzsGeLu7vcnLBtOSDucfnGhNCoxOL5P+Ay3L0iins21z3h3hBpTMzqHHQHviWWqoqXuH6olpl+R2JiWGOuCksXkR/i7rw/H/gMl1jSqfqsLmPiU94AQKwDvzbFSyD/dL+jMDGuMZ33a4FHcA+evEJVh+FqLVcA00XETmVMfGvRxo3NYpccV3VgO+wttMRi6tWYxJIKfBs4XVXfBFDVl3GJ5SzgfRHJiVyIxvggf5DVWKorXuLeLbGYejQmsfRV1adVtSJ4pqq+DYwCegGzIxGcMb7pfBrs+QIO7fI7ktgRSCyBPihjQmjM5caH61g2FziH6D140pjmEbgD32otJxQvgfY9oVW235GYGBfxZzJ447CcHentGhNVwVeGGcc67k2Y6k0sInJBQzeqqlu8dZs6xoox/sjMhczO1oEfcGg37NpoicWEJZway1si8p6IXC4iqfUV9i5HvkpEPsSeNGziWf7pJ/oVkl2g5maJxYQhnBskh+AuL54GbBeRmcA83H0sOzkxgmRf3FVh53vTM4DBzRCzMdFRMBjWzYSjB9wlyMmssuPeEoupXzjjsSwHLhSREcD3gNHADYHFQUUF2Au8CjyhqvMjHKsx0VUwxI2UuG0ZdD/L72j8VbwE2naBzE5+R2LiQNiPdFHVOcAcEekI9AYGAJ1wyaUMWA4sqn4ZsjFxK9+rcBcttsRiHfemARrzrLBSYBOwCFjsvRap6tYIxmWM/9rmQ2YeFC/2OxJ/HT3gRtQceLXfkZg40ZjE8gRwOvBV4Gq85jAR2QL8A/i9qu6PWITG+Cl/MBQt8jsKf21bDqjVWEzYGnOD5B2qeo6qtgP6ANcCD+BqMv8NzBOR9pEN0xifFAxxZ+tHD/gdiX/sjnvTQE26QVJVN6jqq6r6S1UdDlwMFAD3RiQ6Y/xWMPhEB36yKl4CrXOgbYHfkZg40ZjH5m8XkXdF5I8icrOIDAzc3+I9L+x+4Nsi0iLSwRoTdZUd+EncHBbouBfxOxITJxrTx/Ihro/lK7hLjBU4KiIrcB35u3HPCjsbeD8yYRrjk7b57g78oiTtwD92GMpWQd+v+h2JiSMNTiyqeg2AN+7K4GqvsUBLr+hMEdnAiavGfheRiI2JtoLByXtlWMlyqCh334ExYWr00MTelV8fey8AvCaxq4GXgJlAC+A84BrAEouJT/mDYe3bcGR/8g3Ju/Uz997lDH/jMHGlwYnFa/JaBCzBu49FVcsAVPW499BKBb6nqhu9dbpELmRjoiy4A7/HCL+jia6tC929PG3tv7AJX2NqLIdwtZIbOXEPyzagCOgI9ADeDiQVALt50sS1QAd+8eLkSyxFn0HBUOu4Nw3SmD6WYSKSAvTD9auc7r1O9rb3D+AnkQzSGF8lawf+4T3uHp7TrvM7EhNnGtXH4j0PbKX3eiGiERkTiwqS8A78QCLtMsTfOEzcifgIksYkpIKh7uz98F6/I4merQvde8FQf+MwcScuEouIXCsiK0SkQkSG1VFuk4gsE5HFIrIgmjGaBNf1DEBdn0Oy2LoQOpwErTv4HYmJM3GRWHCP5L8amBVG2fNUdbCqhkxAxjRY4HLbwiQ6XylaZLUV0yiNvo8lmlR1FYDYlSnGL63aQ8c+J5qHEt2+bbB3q92/YholXmos4VLgbRFZKCK31VVQRG4TkQUisqCsrCxK4Zm41mWYq7Go1l823lXeGGk1FtNwMZNYROQdEVley2t0AzZztqoOBS4Bvi8i54YqqKpPqeowVR3WqZMNt2rC0HUYHCiFPVv8jqT5bV0IkmqPyjeNEjNNYao6KgLbKPLeS0VkKjCc8PpljKlfcD9Ldnd/Y2luRZ9B3gBo0drvSEwcipkaS1OJSBsRyQp8Bi7EdfobExl5AyG1ZeL3s6i6pjDruDeNFBeJRUSuEpFCYATwpojM8OYXiMh0r1ge8LGILAHmAW+q6lv+RGwSUloLNy5Jol8ZtnMDHN5tHfem0WKmKawuqjoVmFrL/CLgUu/zBtyjZYxpPl2HwYKn4fgxSE33O5rmEaiRWce9aaS4qLEYEzO6nAHlh6Fkhd+RNJ8t86BFJuQO8DsSE6cssRjTEF29+263JnBz2Ja57jhTUv2OxMQpSyzGNER2D2idA4UJ2oF/ZL8bNbLrcL8jMXHMEosxDSHizuYL5/sdSfPYutANatbtTL8jMXHMEosxDdX1S7BjLRzc6Xckkbdlnnvvao/aM41nicWYhup+lnsP/Agnki1zoVN/aJXtdyQmjlliMaahCoZCShp8McfvSCKrogIK50E3618xTWOJxZiGatEa8ge7s/tEsv1zNxyx9a+YJrLEYkxjdD/LPfak/IjfkUROIFFaYjFNZInFmMbodiYcP3JiXPhEsGUetOoAHXv7HYmJc5ZYjGmMyg78T/2NI5IC/Ss2oJ5pIkssxjRGZq4bD/6LBEksB3e6PhbruDcRYInFmMbqPsL1SyTCiJKBBNntLH/jMAnBEosxjdXtTDi4A3as8zuSptv8iRtrxh6VbyLAEosxjRXoZ0mE+1k2feyeKJCe4XckJgFYYjGmsXJOdldRbY7zxHJ4D2xbCj3P9jsSkyAssRjTWCLQYyRs/tjvSJrmi7nuwZM9RvodiUkQlliMaYpe58LuL2DXZr8jabzNn0BKuj0q30SMJRZjmqLn/3Pvm+K41rL5EzcMcYvWfkdiEoQlFmOaolN/aN0xfhPL0QNQtAh6WP+KiRxLLMY0RUqK+1He9FF83s+yZS5UlFvHvYkoSyzGNFXPc2DPFtgdh/0smz4BSbUHT5qIssRiTFP1Ose9x2Nz2OZPIP90aJnldyQmgVhiMaapOvWD1jnxl1iO7IfCBe7KNmMiyBKLMU0l4q4O2xhn/SybP4GKY9D7PL8jMQnGEosxkdDrHNhbCDs3+B1J+Na/D2kZ9uBJE3GWWIyJhN7nu/d17/obR0Ns+MA9odmeD2YizBKLMZHQ4SRo3wvWv+d3JOHZWwxlq6wZzDQLSyzGRErv8939LOVH/Y6kfhs+cO8nWWIxkWeJxZhI6X0+HN3vhviNdRved1ey5Q30OxKTgOIisYjIH0RktYgsFZGpIpIdotzFIrJGRNaJyIRox2mSXK9z3c2Gsd4cpupqLCd9xT05wJgIi5e/qpnAQFUdBHwOTKxeQERSgceBS4ABwA0iMiCqUZrkltHWjRkf6x34patgf4lLLMY0g7hILKr6tqqWe5OfAl1rKTYcWKeqG1T1KDAFGB2tGI0BoPcFULwEDmz3O5LQ1r7t3vtc4G8cJmHFRWKp5hbgP7XM7wJsCZou9ObVSkRuE5EFIrKgrKwswiGapNX7fEBPdI7Hos9nQOdB0LbA70hMgoqZxCIi74jI8lpeo4PK3A2UA5Nr20Qt80LeBq2qT6nqMFUd1qlTp6YfgDEABYPdY/Q/n+F3JLU7uBO2fAonX+x3JCaBpfkdQICqjqpruYiMAy4HLlCt9bkZhUC3oOmuQFHkIjQmDCmp0PciWPMmHD8Gqel+R1TVunfdMMQnX+R3JCaBxUyNpS4icjFwF3CFqh4MUWw+0FdEeolIC+B6YFq0YjSm0imXwOE98MWnfkdS09oZ7jLjgqF+R2ISWFwkFuAxIAuYKSKLReRJABEpEJHpAF7n/h3ADGAV8LKqrvArYJPEep8PqS1gTW1dgT46Xg5rZ7rail1mbJpRzDSF1UVV+4SYXwRcGjQ9HZgerbiMqVXLTOj1ZVgzHS56wD39OBYUzoPDu60ZzDQ7O20xpjmccjHs2gjbP/c7khPWTIeUdHuMi2l2lliMaQ4nX+Le18RIBVoVVrzuHjqZ0dbvaEyCs8RiTHNo1wXyB8Oqf/sdiVO0CPZ8AQPsnmHT/CyxGNNcTr0Sti6EXZv9jgRWvg4paXDKpfWXNaaJLLEY01wGXOneV0z1Nw5VWPmau6CgdQd/YzFJwRKLMc2lQy93v4jfiWXbUti1yZrBTNRYYjGmOQ28GooXw471/sWw4jX3OP9+l/sXg0kqlliMaU6B5rCVr/mz/4oKWP5PN1ZMm47+xGCSjiUWY5pTdjfodiYsfcX1dUTbF7Nh9xdw+g3R37dJWpZYjGlup98AZaug6LPo73vxi9AiE/pbM5iJHkssxjS3gVdDWitYVNtoD83o6AHXBDfgSmjRJrr7NknNEosxzS2jHfT/muvrOHY4evtd9QYc3Q+DrRnMRJclFmOiYfCN7lH6q9+I3j4X/wOyu0P3kdHbpzFYYjEmOnp92f3IL/h7dPZXtgY2zoIzxtsj8k3U2V+cMdGQkgJf+hZs/hi2LW/+/c3/qxsTZui45t+XMdVYYjEmWobc7Drx501q3v0c3guLX4BTr4Y2Oc27L2NqYYnFmGhp3QEGXefuaTm4s/n2s2SK67Qfflvz7cOYOlhiMSaazrwdyg/B/L81z/aPH4M5f4auX4KuZzTPPoyphyUWY6Ip71ToexF8+jgc2Rf57S992d1pf+7PI79tY8JkicWYaPvynXBoV+RrLRXH4aM/QudB0PfCyG7bmAawxGJMtHUd5sadn/3nyNZalr0CO9e72opI5LZrTANZYjHGD+ffAwe3w0ePRGZ7Rw/Cu/dD/un2eHzjO0ssxvih6xkwaAzMeTwyQxfPeQz2boWLfmc3RBrf2V+gMX654F6QFHhrYtMeqb9zI3z8KPS/AnqeHbn4jGkkSyzG+KVdFzjvv2HNm7Dsn43bRkUFvPY9SEmDi34b2fiMaSRLLMb4acT3oetwmP4zNy59Q815zA3mdfHv3aBixsQASyzG+CklFa56ElCYcpMbQyVc696Fd+51nfWDb2y2EI1pKEssxvitY2+45mkoWQ4v3QTHDtW/zpb58Mp46NTfJSa7vNjEEEssxsSCvqNg9GOw/n14/mrYVxK67Jr/wHOjoXVHuPElaJkVvTiNCYMlFmNixZCb4Jq/QtEieGIEzJ3k7tAHd9VYyUp49TZ48XroeBLc8pb1q5iYJNqUyxyjRET+AHwNOAqsB76pqrtrKbcJ2AccB8pVdVg42x82bJguWLAgcgEb0xQlK2H6z93YLZLqrh47etDdUJnWyj3I8rz/hrSWfkdqkpiILAz1G5sW7WAaaSYwUVXLReRBYCJwV4iy56nq9uiFZkyE5Q2A8W9A0Wew5i13tVh6BhQMdR31mZ38jtCYOsVFYlHVt4MmPwW+7lcsxkSFCHQ5w72MiTPx2MdyC/CfEMsUeFtEFopInaMcichtIrJARBaUlZVFPEhjjElWMVNjEZF3gM61LLpbVV/3ytwNlAOTQ2zmbFUtEpFcYKaIrFbVWbUVVNWngKfA9bE0+QCMMcYAMZRYVHVUXctFZBxwOXCBhrjiQFWLvPdSEZkKDAdqTSzGGGOaR1w0hYnIxbjO+itU9WCIMm1EJCvwGbgQWB69KI0xxkCcJBbgMSAL17y1WESeBBCRAhGZ7pXJAz4WkSXAPOBNVX3Ln3CNMSZ5xUxTWF1UtU+I+UXApd7nDcDp0YzLGGNMTfFSYzHGGBMnLLEYY4yJqLh4pEtzE5EyoLHjw+YAiXCnf6IcB9ixxKpEOZZEOQ5o2rH0UNVaHwNhiaWJRGRBuM8ki2WJchxgxxKrEuVYEuU4oPmOxZrCjDHGRJQlFmOMMRFliaXpnvI7gAhJlOMAO5ZYlSjHkijHAc10LNbHYowxJqKsxmKMMSaiLLEYY4yJKEssjSAivxaRpd5zy94WkQJvvojIn0Rknbd8qN+x1kdE/iAiq714p4pIdtCyid6xrBGRi/yMMxwicq2IrBCRChEZVm1ZvB3LxV6s60Rkgt/xNISIPC0ipSKyPGheBxGZKSJrvff2fsYYLhHpJiLvi8gq72/rh978uDseEckQkXkissQ7lvu8+b1EZK53LC+JSIsm70xV7dXAF9A26PN/AU96ny/FDUImwFnAXL9jDeNYLgTSvM8PAg96nwcAS4CWQC9gPZDqd7z1HEt/4BTgA2BY0Py4OhYg1YvxJKCFF/sAv+NqQPznAkOB5UHzHgImeJ8nBP7OYv0F5ANDvc9ZwOfe31PcHY/3u5TpfU4H5nq/Uy8D13vznwS+29R9WY2lEVR1b9BkG9zIlQCjgefU+RTIFpH8qAfYAKr6tqqWe5OfAl29z6OBKap6RFU3Autw49vELFVdpapralkUb8cyHFinqhtU9SgwBXcMcUHd4Ho7q80eDTzrfX4WuDKqQTWSqhar6mfe533AKqALcXg83u/Sfm8y3XspcD7wT29+RI7FEksjicgDIrIFGAv80pvdBdgSVKzQmxcvgod9jvdjCRZvxxJv8YYjT1WLwf1YA7k+x9NgItITGII704/L4xGRVBFZDJQCM3E1491BJ5cR+VuzxBKCiLwjIstreY0GUNW7VbUbbpjkOwKr1bIp36/nru9YvDLVh32O22OpbbVa5vl+LHWIt3gTnohkAv8CflStxSKuqOpxVR2Ma5kYjms+rlGsqfuJi/FY/KD1DJUc5AXgTeBeXLbvFrSsK1AU4dAarL5jCTHsc1weSwgxeSx1iLd4w1EiIvmqWuw1D5f6HVC4RCQdl1Qmq+qr3uy4PR4AVd0tIh/g+liyRSTNq7VE5G/NaiyNICJ9gyavAFZ7n6cB3/CuDjsL2BOoLseqOoZ9ngZcLyItRaQX0Bc3Mmc8irdjmQ/09a7WaQFcjzuGeDYNGOd9Hge87mMsYRMRAf4GrFLVR4IWxd3xiEinwFWfItIKGIXrM3of+LpXLDLH4veVCvH4wp29LAeWAv8GuuiJqy4ex7VbLiPoyqRYfeE6srcAi73Xk0HL7vaOZQ1wid+xhnEsV+HO9o8AJcCMOD6WS3FXIK0H7vY7ngbG/iJQDBzz/j1uBToC7wJrvfcOfscZ5rH8P1zT0NKg/yOXxuPxAIOARd6xLAd+6c0/CXeitQ54BWjZ1H3ZI12MMcZElDWFGWOMiShLLMYYYyLKEosxxpiIssRijDEmoiyxGGOMiShLLMYYYyLKEosxxpiIssRiTIwRkS954+NkiEgbb+yMgX7HZUy47AZJY2KQiPwGyABaAYWq+jufQzImbJZYjIlB3jPC5gOHgZGqetznkIwJmzWFGRObOgCZuFELM3yOxZgGsRqLMTFIRKbhRo7sBeSr6h31rGJMzLDxWIyJMSLyDaBcVV8QkVRgtoicr6rv+R2bMeGwGosxxpiIsj4WY4wxEWWJxRhjTERZYjHGGBNRlliMMcZElCUWY4wxEWWJxRhjTERZYjHGGBNR/x9lfGvosX+YygAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Iteration, E = 8 -15.955308948863635\n", "nL 501\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAEWCAYAAABFSLFOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXxU9bn48c+TBMgGCZAAYUdAkUUWU3C3iFoBK6517cVq63bt7W3ttfizm9altmjv7aZYt6qAWiuKu6IIIiqgBGSVRZYQCEFkh0DI8/vjeyZMwsxkm+TM8rxfr3nNzFmfk+U857uc7xFVxRhjjImWFL8DMMYYk1gssRhjjIkqSyzGGGOiyhKLMcaYqLLEYowxJqossRhjjIkqSywmKYjIRSKyUUT2iMjQZtzv1SLyTnPtzy8isk5Ezo7Stp4SkXuisS3jD0ssCcb7B9/vnUC3eP+k2bWsc4233i4R+VREutay/G9F5Nk6xDJJRP4e9L2FiOwNM+2kuhxfI0wEblXVbFVd2BQ7EJGeIqIikhaYpqqTVfXcpthfhDjWichBEcmrMb3Ii69nHbZx1LHEg2gmuEbGcZWIrPf+tl8WkXZ+x9ScLLEkpu+qajYwBBgK3BFuQS/pPAncAOQCtwIHohTHbODMoO+FwAbgjBrTAD6L0j7D6QEsbeJ9xJKvgCsDX0RkEJDhXzjJQ0QGAJOA7wMdgX3A3yOulGAssSQwVd0CvI1LMGEXAyqAr1S1UlXnq+q2KIUwCzg+6Mr5dOA5IKvGtI9V9RCAiPzLK2ntFJHZ3j8pInKSNz01sHGvemux9zlFRCaIyBoR+VpEXhCRdiLSSkT2AKnAIhFZ4y2vItInaFtV1S8i8m0RKRaR20Rkq4hsFpEfBC2bISIPelekO0Vkjohk4BIpwA6vxHiyiFwrInOC1j1FROZ7680XkVOC5n0gIr8TkY9EZLeIvBP4OYlIuog86x3bDm/djhF+9s8A/xH0fTzwdPACIjJWRBZ6JdWNIvLboNlHHYu3zo9EZLkX3zIRGRa0zhARWewd2/Mikh60r/O9EtMOEZkrIicEzRsqIp9723weSCcMEektIu97P4dtIjJZRHK9ec8A3YFXvZhvD7H+L0Tkk0BJTERuFpGlwbFGwdXAq6o6W1X3AL8CLhaR1lHcR0yzxJLAxFVpjQZWR1jsIFAEvCAibaO5f1UtBtbjkge4ksqHwNwa02YHrfYm0BfoAHwOTPa29QmwFzgraNmrgCne5/8CLsSVkDoD3wB/U9Vyr/QGMFhVe9cx/E5ADtAFuB74W9DPZyJwInAK0A64HajkSEks16ty+zh4g151yOvAn4H2wEPA6yLSvsYx/cA7/pbAz73p4714unnr3gTsjxD/J0AbETneS8aXAzWrL/fikk8uMBa4WUQu9OYddSwichnwW2+dNsAFwNdB2/secB7QCzgBuNY77mHAE8CNXuyTgOle0m8JvIxLhO2AfwGXRDguAe7H/Y6P934evwVQ1e/jSsTf9WL+Q4j1/4j7m/+liPQF7gOuUdWjSuki0t1LhOFeV4WJcQCwKPBFVdd4+zw2wnElFlW1lxsv7QlgK7CkDsuegTvpVQCX1pg3Hljlvcb7cBzrgD3Ablxp5D3cySHc8o94r9tx1VFtven3Ag+GWee3wLN1jOcp4E+4i5itQCbupBiY9g1wZph1c71jyPG+3wM84X1ujTsx9vC+LwdGBa1bABwC0rzvCvQJml/z+1PAPd7nb+NO2mlB87cCJ3kx78clqZrx9vS2G7zetcAc7/P3gXk11vkYuNb7/AHwy6B5twBveZ+vwyXkE+r4N3A28EvcSfg84F0gzYuvZ5j1/hf4U4RjeRv4SYR9XhP0/Q/AI97nh4Hf1Vh+Je4i4AygBJCgeXMDv4s6HOuFwMKax17LOj2B7d7fzB1N8D/4HnBTjWmbgG9He1+x+rISyxFP4f4B62ID7oQxJXiid0X6G2AEMBz4TbRLAXV0oaq2xp0g+wF5oRYSkSzc1fgf1F3dvQvM8GI+BZgRhVhm404eg4C1qroPmBM0LQP41IsnVUR+71Vn7cKdJAiKfwquSqEVcDHwuaqu9+b1AKYFriZxJ43DuDruhvhaVSuCvu8Dsr1Y0oE1DdhmZ1wJLth6XKkoYEuIfYK7on8beE5ESkTkDyLSopb9PYMrAV1LjWowABEZISIzRaRMRHbiEn7IvxVPNyIfd7jYewC3BV/te9vq7L02qXf29dT8GQXH3EFEnhORTd7fyLO1xHwUVV0HzMQlmL/VZ9062oMr0QVrg7vYSwqWWDyqOht3FVPFq899S0Q+E5EPRaSft+w6VV2Mq/4I9h3gXVXdrqrf4E7UdU1WUaeqs3AJc2KYRVJwbQ8V3vITgAW4apRM4K0ohDEbGIyravnQm7YUd2IZC8zXI9UQVwHjcFfbObh/fHDVH6jqMtxJZzTVq8EANgKjVTU36JWuqpvCxLUPd4wBnep4PNtwnRtCVanVNlR4Ce4kG6w77mo2IlU9pKp3qWp/XNI/n+ptKKHWWY9rxB8DvBRikSnAdKCbqubgSq4SWD3E8hsJfdy12QjcW+N3k6mqU4HNQBcRkaDlu0fY1v1ebCeoahvgmqCYw8VdjYiMAU7GlSz+GGG57l5bTbjX1WFWXYr7mw9s5xigFfBlbbElCksskT0K/FhVT8TVddfWs6ML7p8ooJjqV6N++F/gHBE5qgFfVXfjksffRaSjV9/9Pu7kcRCIdEWc4jUoB16tQi2kqquBUuAneInFuzr91JsW3L7SGijH1dtn4uq/a5qCa085A1cfH/AIcK+I9AAQkXwRGRch/iLgKq+UdB7Ve6+FpaqVuGrTh0Sks7f+yd7xl+EuNo4Js/obwLHiuqKmicjlQH/gtdr2KyIjRWSQ116yC1fNd7gOIV8PnKWqe0PMaw1sV9UDIjIcl6wDQh3LY8DPReREcfoEft61+Adwk1dCEhHJEtdxoDWuKrAC+C/vZ3IxrrQfTmtciWCHiHQB/qfG/FLC//wR1xniceCHuGrr73qJ5iiqukFdW0241+Qwu5nsbfd0r1bgbuAl7/8tKVhiCUNcN9xTgH+JSBGuwbGgttVCTPP1gTeqWoarBvlVmEWuwf0zLsIlxatxDdOCO4GGcyWurSHwilRFMhvIBz4KmvYhroE6OLE8jSuRbAKW4UpONU3FVfG9r9V7r/0f7ur7HRHZ7a07IkJMPwG+C+zAHfPLEZat6efAF8B8XCn3ASDFq+a7F/jIq/Kpdm+Oqn6NK2nchkuetwPna9164XUCXsQlleW4Hne13kukqmtUdUGY2bcAd3s/r18DLwStd9SxqOq/vGlTcNU6L+Ma3GuLYQHwI+CvuDa11XgN+6p6EFetea0373JCl64C7gKGATtxHSFqLns/rmF+h4j8vObKuIvFV1T1De/3cT3wWI0OFI2iqktx1YqTcW1zrXE/66Qh1as2k5u4G8deU9WBItIGWKmqYZOJiDzlLf+i9/1KXAPdjd73ScAHXpHfGGOSgpVYwlDVXcBXXhdLvCL84FpWexs4V0Taeg3g53rTjDEmaVhi8YjIVFx973Hibo67HldFcr2ILMI1yI3zlv2WiBQDlwGTRGQpgKpuB36HqyKZD9ztTTPGmKRhVWHGGGOiykosxhhjoiquRi5tKnl5edqzZ0+/wzDGmLjx2WefbVPV/JAz/b71PxZeJ554ohqTiMaPH6+4Lu/VXiNGjPA7tJhXUlKiV155pR533HGakpKi48ePP2qZJ598MuTPd//+/RG3PWPGDD355JM1OztbO3XqpLfffrseOnSoav5XX30Vcrtvvvlm1TIzZ84Muczy5curlnnhhRf0xBNP1JycHM3MzNTBgwfrU0891fgfjqoCCzTMOdVKLMYkuLPPPptnnnmm2rSWLVs26T4rKipITU2l+g318aW8vJy8vDwmTJjAo48+Gna5zMxM1qypfhtXenr4wZIXL17MmDFjuOOOO3j66afZtGkTN910E4cPH2bixOqDZLz11lsMHnykM2q7dkffNrR06dJq0/PzjxQi2rdvzy9/+Uv69etHixYteO2117j++uvJz89nzJiQ94VGR7iMk0wvK7GYRDV+/HgdO3ZsxGUAnTRpkl566aWamZmpvXr10meeeabaMsXFxXr55Zdrbm6u5ubm6pgxY/TLL7+smv+b3/xGBwwYoE8++aQec8wxmpKSort379Y9e/bo97//fc3KytIOHTrofffdp2PHjq26+r/rrrt0wIABR8V0yimn6I9//OPG/wCiJDjmYE8++aRmZWXVa1t33HGHDhkypNq06dOna3p6uu7atUtVj5RY5s+fH3Y7gRJLWVlZvfY/dOhQnTBhQr3WCYUIJRZrvDfGcPfddzNu3DgWLVrE5ZdfznXXXcf69W4syH379jFy5EjS09OZNWsWH3/8MQUFBZx99tns27evahtfffUVU6ZM4V//+heLFi0iPT2d2267jVmzZjFt2jTef/99Fi1axIcffli1znXXXceKFSuYN29e1bSVK1cyd+5crr/+egAmT55MdnZ2xNfkyeFGV2l6+/fvp0ePHnTt2pXzzz+fhQsjP6C0vLz8qBJNRkYGBw4c4LPPqj/v7uKLL6ZDhw6ceuqpvPjiiyG3V1hYSEFBAaNGjWLmzJlh96uqvPfee6xcuZIzzjgj7HJRES7jJNPLSiwmUY0fP15TU1M1Kyur2uv222+vWgaodgV76NAhzcjIqCq1PP7449qnTx+trKysWqaiokLbtWunzz//vKq6EktaWppu2bKlapndu3drixYtdOrUqVXT9uzZo7m5udWu/seOHas33nhj1ffbb79dg/8nd+3apatWrYr4ClzpN5VwJZa5c+fqU089pQsXLtTZs2frJZdcohkZGdVKczW9/fbbKiL69NNP66FDh7S4uFhPP/10BXTKlCmqqlpWVqYTJ07Ujz/+WOfPn6+/+tWvNCUlpVpJcsWKFfrwww/rggULdO7cuXrzzTeriOisWbOq7W/Hjh2alZWlaWlp2qpVK33sscei8jMhQonF95N6LLwssZhENX78eB05cuRRJ+Lg6pPgE1pA9+7d9cEHH1RV1VtuuSVkchIRve+++1TVJZaePXtW20ZRUZECumbNmmrTTzvttGon6ZdffllzcnJ03759WlFRoQUFBfr3v/89mj+GWgUfV3CSCwiXWGqqqKjQgQMH1lqN9+CDD2qbNm00NTVVMzMz9f7771egKlGHcvPNN+ugQYMibnf06NH63e9+t9q0w4cP66pVq3ThwoU6ceJEbdOmjc6YMaPWY6lNpMQSV4333ii0/4cb6v0xVf19jfmtcAMZnogb5O9ydc9eMCZpZWZm0qdPn4jLtGhRfSBrEaGy0j0VorKykiFDhvDcc88dtV5wo3FWVla1ee7cQ60N+GPHjiUzM5N///vf5OTksGPHDq688sqq+ZMnT+bGG2+MuI1JkyZx9dXhRrGvXVFRUdXnNm1qPkql7lJTUyksLGTVqlURl/vZz37GT3/6UzZv3kzbtm1Zt24dd9xxB7169Qq7zogRI3jyyScjbnfEiBFH/Z5SUlKqfv9Dhgxh+fLl3HfffYwaNaqOR1V/cZNYvOHC/wacgxuOfr6ITFf3jI6A64FvVLWPiFyBG3X28uaP1pjEMWzYMKZOnUpeXh65ubl1Xq9Pnz60aNGCefPmVZ0w9+3bx5IlS+jd+8hjXdLS0rj22mt54oknyMnJ4eKLL662nwsuuIARIyINVA0dOzb0eW5HYo0GVWXx4sXVenKFIyJ07twZgKlTp9KtWzeGDRsWdvmioiIKCiIPsF6XZSorKykvL681vsaIm8SCe0bDalVdCyAiz+HG7gpOLOPwnn+NG2L8ryIiGrh0irI/v7eKisNBz/oKujKreY0WfNEmNeZWnxd6uvse/sov3PaP2kak7UdYL1wckY8z9DqRlqs5M3hem4wWdGqTTt8O2bTNatrusomkvLycLVu2VJuWmpparVtqJFdffTUTJ05k3Lhx3H333XTv3p2NGzfyyiuvcNNNN9G3b9+Q62VnZ3Pdddfxi1/8gry8PAoKCrjnnnuorKw88vewbztsW8UPT+/KAw98QIqk8M5f/hvm/aNqO629VzXV/qXV1U+Emxd2vRrzQ8wr+rIYgF3FK0jZkUHRU7fTMi2N/r1cIrvribc4aUAP+nbNY9fecv784mwWLyri4RtPhzl/AuCvL37IX6fNYcWzE6r298epMzlveD9SBF6a/QW/f/o9XvjtNaR+9BAA/3xrAS1SUxnatzMpIrw6dzl/e+wtHrhxNMz6A6jyvy/OoWentgzo2ZGDhyp4dsZCXn55Fv++62qYeT8A9z47kxHHd+WYgnaUH6rgjU+/5Jmn3+YvPz4f3r8XWmbBaf9d86fbaPGUWEI9RKvmZUzVMqpa4T1utT3uqX/ViMgNwA0A3btHemBdeI/MWsP+Q+5ZS02Tukw4Q7rl8p8j+3BO/8ZdqSaDGTNmHHUV26VLF4qLi+u0fmZmJrNnz2bChAlcdtll7Ny5k86dOzNy5Ejato385O2JEyeyd+9eLrjgArKzs/npT39KaWkp6QfK4PFzYeM8QDkGOLN7Cut3VvLtLY/CG7Fx/8vQu3ZV+/7qR0vokSOs+2+X6nYsPcAN/z7Elj1KTithaEEKs8enM7zkSfe8UGBb0QFWbjgI791VtZ0339zLvU+8TvlhGNwxlVe+15LRh6bB+9PcAssOcs9HB1m/s5JUgWPbp/DEd1tyTd5smOkeYXRwZTk/n3KQTbuVjDQY0CGV16/KYEzlqzDrVQD2rDzAzS9XULyrkow06JeXytPjWnFlziyYPQuyOzRJYombQSi94eu/o6o/9L5/Hxiuqj8OWmapt0yx932Nt8zXobYZUFhYqAsWhHsWUuMF/4xr/rg13HJHbSN4HQ07L9L04PWiEcdRF4Nhtn/0NhpwnAo79x+iZOcBvijewUufb2Lttr3ccMYx3DG6X1zfiJdMysvL6VHQnv8prOC2scfB4KugyzBo05n+Z17E1Zdfwp3/81Nv6aOK7DW2JhHm1XG9o+Y3dF5t8dRxXjTjaeL/CRH5TFULQ82LpxJLMe456QFdqbomOGqZYhFJwz033fdh66tVHUX8XdvJMZIObdLp27E1Zx6bz01n9uauV5fx6Oy19MrL4srhDSt1mqa1cOFCli9fzvDhw9m9ezcPTLiF3Xv2cvkVP4DxkyC1BVu3bmXq1Kms21DMjT++DbLz/A7bNFI8JZb5QF8R6YV7dO0VVH9GN7hH047HPVflUtzja+OjSGbqJS01hbsuGMDqrXv4/ZsrGDOogJyMFrWvaJrdQw89xMqVK0lLS2NIuwPMnnASXa/9B6SkAq7hPS8vj0mTJpGXZ0klEcRNYvHaTG7FPZExFXhCVZeKyN24/tTTgceBZ0RkNa6kcoV/EZumlpIi3Dn2eM7/yxxe/KyY608L31XT+GPo0KFUVTPP+we88XO4/v+qkgpUrxo1iSFuEguAqr4BvFFj2q+DPh/APdXRJImBXXIY0i3XEks8WPQcdDoBug33OxLTxGysMBP3Rg/sxPLNuyjZsd/vUEw4uzbDpgXQf5zfkZhmYInFxL2z+nUA4MNVZT5HYsJa5w082fdcf+MwzcISi4l7fTpkk5PRgqKNO/wOxYSz8VNomQ0dB/gdiWkGllhM3BMRBnfLZeEGSywxa+M86HJitUZ7k7gssZiEcEKXHL4s3U15xWG/QzE1HT4EW5dDQe3jZ5nEYInFJITeHbKoVNjw9b7aFzbN6+vVUHnIqsGSiCUWkxB652cDsKZsr8+RmKOULnXvHfr7G4dpNpZYTELoleeeBbJ22x6fIzFH2boMJBXyj/M7EtNMLLGYhNA6vQV52a1Yv82qwmJO6TJo3wfSWvkdiWkmllhMwuicm87mXQf8DsPUtHUpdLRqsGRiicUkjE5t0tmy0+6+jynlu2HHBmu4TzKWWEzCKMhJZ/NOK7HElK3L3XsHSyzJxBKLSRidcjLYfaCCPeUVfodiAgI9wqwqLKlYYjEJoyAnHcCqw2LJ1mVuKJccexBbMrHEYhJGXrbrdfT1noM+R2KqlC6DDsdDip1qkon9tk3CaJfVEoDtey2xxARV1yPMboxMOpZYTMKoSiz7LLHEhN1bYP831iMsCVliMQkjN9M98/4bK7HEhq02lEuyssRiEkZ6i1SyWqayfe8hv0Mx4NpXwEosScgSi0ko7bJb8o1VhcWGrcsguxNktvM7EtPMLLGYhNIus6U13seKUhvKJVlZYjEJJTfTSiwx4XAFlK209pUkZYnFJJQ2GS3Ytd/aWHy3fS0cLrf2lSRlicUklDbpaew+YEO6+M56hCU1SywmobROb8HuAxWoqt+hJLfSZSAp9nCvJGWJxSSU1ulpHDxcSXlFpd+hJLfSpdC+L7TI8DsS4wNLLCahtMlwN0nuOmDtLL4qXWLtK0nMEotJKG3S0zgnZQEsnOx3KMnrwC7Ysd4SSxKLi8QiIu1E5F0RWeW9tw2z3GERKfJe05s7TuO/1ulpXJg6h5wFf/Y7lOQVeLhXx4H+xmF8ExeJBZgAvKeqfYH3vO+h7FfVId7rguYLz8SK0/rkc+7IUbTctR7K9/gdTnIqXeLercSStOIlsYwD/ul9/idwoY+xmBjWMi2FFgUDERTKVvgdTnIqXQqtciCnq9+RGJ/ES2LpqKqbAbz3DmGWSxeRBSLyiYhETD4icoO37IKysrJox2v8FLhSDlw5m+YVaLgX8TsS45M0vwMIEJEZQKcQs+6sx2a6q2qJiBwDvC8iX6jqmlALquqjwKMAhYWFdtNDIsnt4R6HG3jeumk+lZXuHpbBV/gdifFRzCQWVT073DwRKRWRAlXdLCIFwNYw2yjx3teKyAfAUCBkYjEJLCXF3fEdGLbdNJ+dG+DgbuhkDffJLF6qwqYD473P44FXai4gIm1FpJX3OQ84FbAzS7LqOMBVydgd+M0rUEq0HmFJLV4Sy++Bc0RkFXCO9x0RKRSRx7xljgcWiMgiYCbwe1W1xJKsOg6AAztgV4nfkSSX0qWAQH4/vyMxPoqZqrBIVPVrYFSI6QuAH3qf5wKDmjk0E6sCV8ylSyGni7+xJJPSJdCuF7TK9jsS46N4KbEYUz+BB0xZz7DmVbrU7l8xllhMgkrPgZxu1jOsOR3cB1+vsfYVY4nFJLCOAyyxNKetywG1EouxxGISWMcBsO1LqCj3O5LkUDWUi5VYkp0lFpO4Og4APeyevW6aXulSd2Nqbg+/IzE+s8RiElfgynmr9TpvFqVL3Y2pKXZaSXb2F2ASV7vekNrKeoY1B1V7uJepYonFJK7UNOjQzxrwm8OuEndDqiUWgyUWk+g6DrTE0hxsKBcTxBKLSWwd+sOeUthjj0ZoUlU9wvr7G4eJCZZYTGKzZ7M0j9KlkNPd3Zhqkp4lFpPYOnnDx1liaVqlS2yofFPFEotJbFl50LoAtlhiaTKHDsC2VdZwb6pYYjGJr+NA2PKF31Ekrm0r3Y2olliMxxKLSXydBrmTnw3t0jSsR5ipwRKLSXydBkJlhQ3t0lRKl0JaOrQ7xu9ITIywxGISX0evAd+qw5pG6RLocDykpPodiYkRllhM4mvfG9IyrGdYU7GHe5kaLLGYxJeS6m7csxJL9O3ZCnvLjpQKjcESi0kWnQa5xKLqdySJpeqOeyuxmCMssZjk0HGgGyRxZ7HfkSSWLZZYzNEssZjk0OkE927tLNFVugRad4bMdn5HYmKIJRaTHAKDI1o7S3Rt+eLIsDnGeCyxmOTQqjW07WWJJZoO7Xf3BhWc4HckJsZYYjHJo9MgqwqLpq3L3VAuVmIxNVhiMcmj0yDYvhbKd/sdSWLYsti9d7ISi6nOEotJHlVD6C/zN45EseULaNUGcnv4HYmJMZZYTPIIDJIYuNI2jbN5sUvWKXYaMdXFxV+EiFwmIktFpFJECiMsd56IrBSR1SIyoTljNHEgpyuk51o7SzRUHnZDuVj7igkhLhILsAS4GJgdbgERSQX+BowG+gNXiog9gNscIXLkDnzTONvXwqG91r5iQoqLxKKqy1W1tjHPhwOrVXWtqh4EngPGNX10Jq50HOjaWCoP+x1JfKtquLcSizlaXCSWOuoCbAz6XuxNC0lEbhCRBSKyoKysrMmDMzGi0yCo2A9fr/E7kvi2eTGktID8fn5HYmJQzCQWEZkhIktCvOpa6pAQ08KOOKiqj6pqoaoW5ufnNyxoE386eQ34pVYd1ihbvoAO/SCtpd+RmBiU5ncAAap6diM3UQx0C/reFShp5DZNosnvBylp7sQ48BK/o4lPqq4qrO93/I7ExKiYKbFEwXygr4j0EpGWwBXAdJ9jMrEmrRXkHWcN+I2xp9Q9g8XaV0wYcZFYROQiESkGTgZeF5G3vemdReQNAFWtAG4F3gaWAy+o6lK/YjYxrGCwayMwDbPZGu5NZDFTFRaJqk4DpoWYXgKMCfr+BvBGM4Zm4lHBCbBoCuzeAq07+R1N/KnqETbQ3zhMzIqLEosxUVUw2L1vXuRvHPFqy2Jo2xPSc/yOxMQoSywm+QSqcCyxNMyWL+zGSBORJRaTfFq1hvZ9LLE0xIFd7q57Sywmgga1sYjIscAAoAPuXpEyYImqropibMY0nYLBsHGe31HEn1KvP4w93MtEUOfEIiLHAzcBlwEdA5O9d/WWKQVeACap6vIoxmlMdBUMhiX/hr1fQ1Z7v6OJHzaUi6mDWhOLiPQGHgAuAvYDHwIfA2uAr3HJpR3QBzgJ+CHwYxF5CfiFqq5tmtCNaYRAA/6WRdD7LH9jiSdbFkNmHrQu8DsSE8PqUmJZBnwBXAu8pKp7Iy0sIlnApcB/eeumNzJGY6Iv0EawebEllvoIPINFQo2gZIxTl8b773ljaj1TW1IBUNW9qvpPVT0RuLzxIRrTBDLbQW53a8Cvj4qDULbC2ldMrWpNLKr6SkM33ph1jWlyBYMtsdRH2XI4fNB6hJla1bu7sYh8rykCMabZFQyG7WtcF1pTu5Ii9955qL9xmJjXkPtYpojIjVGPxJjmVjDEvduAlHVTshBa5bubJNwAABX/SURBVEC7Y/yOxMS4hiSWp4C/i8idoWaKyMki8mGjojKmOdjQLvWzuQg6D7aGe1OreicWVf0h8AfgdyLyp8B0ETnO62I8B7Cysol92R1ct1lLLLWrOOhujgyU8oyJoEF33qvqHd7NkA+KSD6wB7gOd6PkJODu6IVoTBPqdMKRm/5MeFuXuYb7zpZYTO0aM2z+P4CxwFW4hPIc8Cu7IdLElYLBsPpdOLgPWmb6HU3s2mwN96buGtIrrIWI/AR35/1ZwOe4xNIC2Bjd8IxpYgWDQSvdFbkJr2ShGya/bS+/IzFxoCGN96uAh3ADT16gqoW4UssFwBsikh3F+IxpWlUN+EX+xhHrSopc+4o13Js6aEhiSQV+BAxW1dcBVPUFXGI5CZgpInnRC9GYJpTTFTLaHblHwxytotw13Fv7iqmjhrSx9FXVAzUnquo7InI28DowFzi2scEZ0+REXLtByUK/I4ldW5dB5SHrEWbqrCHdjY9KKkHzPgVOxwaeNPGkyzDYutw14Juj2R33pp6i/gRJ7zksp0Z7u8Y0mc7DQA9bt+NwShZCeq57zr0xdVBrYhGRUfXdqKpu9NY9uyFBGdOsugxz75s+9zeOWLW5yLWvWMO9qaO6lFjeEpH3ReR8EUmtbWGvO/JFIjILeKPxIRrTxFp3gtadocQSy1EqyqF0mbWvmHqpS+P9UFz34unANhF5F5iHu49lO0eeINkX1yvsLO/724D9NZr40GWYlVhCKV3qGu6tR5iph1oTi6ouAc4VkZOBW4BxwJWB2UGLCrALeAl4WFXnRzlWY5pO56Gw4jXYvwMycv2OJnbYHfemAerc3VhVPwY+FpH2QG+gP5CPSy5lwBJgoapWNkWgxjSpQDtLyULoPdLfWGJJoOE+t4ffkZg40pBeYVuBqcD5QCtgBTBDVT+zpGLiVuCK3NpZqispcj8ba7g39dCQxPIwUAKcgxvF+BVgg4isE5F7bEgXE5cy2rpxsKyd5YhD+93NkVYNZuqpITdI3qqqp6tqDtAHuAy4F1eS+X/APBFpG80gReQyEVkqIpUiUhhhuXUi8oWIFInIgmjGYJJAl2F2B36wzYugsgK6hv2XMyakRt0gqaprVfUlVf21qg4HzgM6A7+JSnRHLAEuBmbXYdmRqjrEGxzTmLrrPAx2bYLdpX5HEhuKvWuzLvavZOqnIcPmbxOR90TkQRH5vogMDNzfoqrv4KrHfiQiLaMVpKouV9WV0dqeMSFVNeBbdRgAmxZATjdo3dHvSEycacgglLOAwcC3cV2MFTgoIkuBImAHbqywU4GZ0QmzzhR4R0QUmKSqj4ZbUERuAG4A6N69ezOFZ2JawWCQFFcddtxov6PxX/Fn0OVEv6MwcajeiUVVLwHwGumH1HhdjespBvCuiKzFJZuFqnp/pO2KyAygU4hZd6rqK3UM71RVLRGRDt7+V6hqyOozL+k8ClBYWKihljFJpmUW5PezBnyAPVth5wYYcYPfkZg41OBHE6vqHmCO9wLAqxK7GHgeeBdoCYwELgEiJhZVbfS4Yqpa4r1vFZFpwHDq1i5jjNN5GHz5Jqgmdxdba18xjdCQNpalIvKsiPyPiJwjIvmBeap6GBiFq5K6RVVHqWo+0OR1TSKSJSKtA5+Bc3GN/sbUXZdhsO9r+Gad35H4a9MCkNQjT9g0ph4aUmLZjyuVXIU3pIuIbMHd29Ie6AG8o6pfBVZQ1U2NCVJELgL+grvT/3URKVLV74hIZ+AxVR0DdASmibvKTAOmqOpbjdmvSUJdv+Xei+dDuyR+vnvxAug4AFpm+h2JiUMNaWMpFJEUoB+uXWWw9zrW296zwM+iGaSqTgOmhZheAozxPq/14jCm4Tr0hxZZLrGc8D2/o/FHZaXrwDDoUr8jMXGqQW0s3tAty7zXlKhGZIyfUtNcddjGeX5H4p9tX0L5LmtfMQ0W9SdIGhP3un4LSpck76OKN3kN93bHvWkgSyzG1NT1W24ok8CQ8cmmeAG0yoH2ff2OxMQpSyzG1BTcgJ+MNi2ALkMhxU4PpmHsL8eYmrLz3UjHyZhYDu5zjyK29hXTCJZYjAml67dg43x3o2Qy2VwEetjaV0yjWGIxJpRuw2HPFthZ7HckzSvQG85KLKYRLLEYE0rgir04ybodb/gE2vdx1YHGNJAlFmNC6TgQ0jKOjJmVDFRh46fQ7SS/IzFxzhKLMaGktnCP5E2mGyW3rYL926H7CL8jMXHOEosx4XT7lns876EDfkfSPDZ+4t6txGIayRKLMeF0GwGVh9y4Wclgw6eQ0Q7y7MZI0ziWWIwJp/vJ7n3DXH/jaC4bPobuJyX3c2hMVFhiMSaczHbuiZLrP/Y7kqa3pwy2r3GlNGMayRKLMZF0P8n1lKo87HckTWvjp+69u7WvmMazxGJMJN1PcUPIly71O5KmtfETSG0JBUP8jsQkAEssxkTSI9DO8om/cTS1DZ9C52HQIt3vSEwCsMRiTCS53aFN18RuwD+03/V8s/tXTJRYYjGmNj1Odg34iTogZclC163a7l8xUWKJxZjadD/JDUj5zVd+R9I01nulMesRZqLEEosxtel+intP1G7H6+ZAhwGQ1d7vSEyCsMRiTG3y+0F6rruBMNFUHHRdjXue5nckJoFYYjGmNikp7i789QnYgF/yORzaB71O9zsSk0AssRhTFz1PdXem7yrxO5LoWvehe+9xqr9xmIRiicWYuuh1hnv/6kN/44i2dXPcs2cy2/kdiUkglliMqYuOg1w7y7rZfkcSPRUH3Y2R1r5ioswSizF1kZLiTsBfJVBi2fQZVOy3xGKizhKLMXXV60zYsQG+Wed3JNGxbg4g1r5ioi4uEouI/FFEVojIYhGZJiK5YZY7T0RWishqEZnQ3HGaBBfoOZUo7SzrPrT2FdMk4iKxAO8CA1X1BOBL4I6aC4hIKvA3YDTQH7hSRPo3a5QmseX3g6z8xKgOqyi3+1dMk4mLxKKq76hqhff1E6BriMWGA6tVda2qHgSeA8Y1V4wmCYi43mFfzY7/ccOKF0DFAdeN2pgoi4vEUsN1wJshpncBNgZ9L/amhSQiN4jIAhFZUFZWFuUQTcLqebobN+zr1X5H0jhrZ4KkuuMxJspiJrGIyAwRWRLiNS5omTuBCmByqE2EmBb2slJVH1XVQlUtzM/Pb/wBmOQQuJ9l7Qe+htFoq9+DroWQEbK50phGSfM7gABVPTvSfBEZD5wPjFINWQ9RDHQL+t4VSLDbpI3v2h0DOd1cYhn+I7+jaZh9291Q+d+2/i2macRMiSUSETkP+AVwgaruC7PYfKCviPQSkZbAFcD05orRJAkR6DMK1s6Cw4f8jqZh1n4AKPQ+y+9ITIKKi8QC/BVoDbwrIkUi8giAiHQWkTcAvMb9W4G3geXAC6qa4A8qN77oPQoO7obi+X5H0jBr3odWOe5RxMY0gZipCotEVfuEmV4CjAn6/gbwRnPFZZLUMWe6hu/VM6DHKX5HUz+qsGamO4bUuPj3N3EoXkosxsSO9Bz3tMXVM/yOpP62rYJdxVYNZpqUJRZjGqLPWbB5EeyJs67qa953771H+huHSWiWWIxpiD5eJ8bAiTperHkP2vWGtj39jsQkMEssxjREp8GQmedO1PHi0H43zlmfUX5HYhKcJRZjGiIlxbVTrH4PKiv9jqZu1s5yw+Qfe57fkZgEZ4nFmIY69juwb5t7rkk8+PJNaJltA0+aJmeJxZiG6nM2pKTBytf9jqR2qvDl266UldbK72hMgrPEYkxDZeS6h2StDDUmaozZXAS7N8Nxo/2OxCQBSyzGNMZxY6BsBXy9xu9IIlv5FiDQ91y/IzFJwBKLMY0RKAHEeqnlyzeh23DIyvM7EpMELLEY0xhte7jH+8ZyYtlV4m7mtN5gpplYYjGmsY4bDRvmuuHoY9EKr3PBcWMiL2dMlFhiMaax+o0FrTxyAo81S1+G/H7QoZ/fkZgkYYnFmMYqGOKGSFn6kt+RHG33Flj/EQy4yO9ITBKxxGJMY4nAgIvdne17v/Y7muqWTQcU+l/odyQmiVhiMSYaBlwEehiWx9hDS5e9DPnHWzWYaVaWWIyJhk6DoH2f2KoO270F1s+FAVZaMc3LEosx0SDiSi3r5sCerX5H4yx7BasGM36wxGJMtAy42PUOW/qy35E4i55z99hYNZhpZpZYjImWjv2h4yBYNMXvSGDrCij5HIZc5XckJglZYjEmmoZcBSULoXSZv3EsmuJGXh70PX/jMEnJEosx0XTC99wJvWiyfzEcroBFz7sBJ7Pz/YvDJC1LLMZEU1aeG5Nr8fNw+JA/MaydCXu2WDWY8Y0lFmOibcjVsLcMVr3rz/4XPgsZ7aDvd/zZv0l6lliMiba+50B2J1jwRPPve9dmWPEaDL4C0lo2//6NwRKLMdGX2gIKfwCr323+B4B99hRUHoZv/bB592tMEEssxjSFE691jfjzH2u+fVYchM+edCWm9r2bb7/G1GCJxZim0LqTu+N94WQo39M8+1w+HfaUwvAbm2d/xoQRF4lFRP4oIitEZLGITBOR3DDLrRORL0SkSEQWNHecxlQz4kYo3+ka05uaKnzyMLTrDb3Pavr9GRNBXCQW4F1goKqeAHwJ3BFh2ZGqOkRVC5snNGPC6DYcup8Cc//sqqma0lezYNMCOPk/ISVe/q1NooqLv0BVfUdVK7yvnwBd/YzHmDo74zbYtQkWTW3a/cz6I7QugKHXNO1+jKmDuEgsNVwHvBlmngLviMhnInJDpI2IyA0iskBEFpSVlUU9SGMA6D3KPWFyzp/cHfFNYf1cWD8HTv0JpLVqmn0YUw8xk1hEZIaILAnxGhe0zJ1ABRBuvIxTVXUYMBr4TxE5I9z+VPVRVS1U1cL8fBv2wjQRETjzdvjmKyhqgrYWVXj/XsjMg2Hjo799Yxogze8AAlT17EjzRWQ8cD4wSlU1zDZKvPetIjINGA7MjnasxtTLcWOg20kw8z4YeCm0yo7etle+4UorYx+ElpnR264xjRAzJZZIROQ84BfABaq6L8wyWSLSOvAZOBdY0nxRGhOGCJz7O9cVeO5forfdioPwzi8h7zgYdm30tmtMI8VFYgH+CrQG3vW6Ej8CICKdReQNb5mOwBwRWQTMA15X1bf8CdeYGroNdw8Cm/Mn2LYqOtv85O+wfS18515IjZnKB2NipyosElXtE2Z6CTDG+7wWGNyccRlTL+f9Hta8B6/+BMa/1rhuwWVfuqq1fudDn4i1yMY0u3gpsRgT/1p3hHPvhfUfwaePNHw7hw/BK7e4NpWxD7mqNmNiiCUWY5rT0Gvg2NHw7q9g4/yGbeOdX0LxfBgz0SUrY2KMJRZjmpMIXPQwtOkMz18D36yr3/qfP+1KOyfdAoMubZIQjWksSyzGNLeMtnDVC3C4HJ6+EHZsrNt6i5537TPHjIRzfte0MRrTCJZYjPFDh+Ph6hdh33Z4bBSs/zj8spWHXUP9tBuhx6lwxRTrBWZimiUWY/zStRCufwdaZMCTo+GVW2HLEnc3PUBFOSx/Df4xEmY9AIOvdCUduxHSxDi77DHGTx36wU0fwcx7Yf7jsPAZV1XWqg3s3gyHD0JOd7jkcRh4ifUAM3HBEosxfmuVDefdD6f9zA3RUrIQDu5xoxX3PB16j3SPOzYmTlhiMSZWZOfDiePdy5g4Zm0sxhhjosoSizHGmKiyxGKMMSaqLLEYY4yJKkssxhhjosoSizHGmKiyxGKMMSaqLLEYY4yJKtHAuERJTETKgPUNXD0P2BbFcPySKMcBdiyxKlGOJVGOAxp3LD1UNT/UDEssjSQiC1S10O84GitRjgPsWGJVohxLohwHNN2xWFWYMcaYqLLEYowxJqossTTeo34HECWJchxgxxKrEuVYEuU4oImOxdpYjDHGRJWVWIwxxkSVJRZjjDFRZYmlAUTkdyKyWESKROQdEensTRcR+bOIrPbmD/M71tqIyB9FZIUX7zQRyQ2ad4d3LCtF5Dt+xlkXInKZiCwVkUoRKawxL96O5Twv1tUiMsHveOpDRJ4Qka0isiRoWjsReVdEVnnvbf2Msa5EpJuIzBSR5d7f1k+86XF3PCKSLiLzRGSRdyx3edN7icin3rE8LyItG70zVbVXPV9Am6DP/wU84n0eA7wJCHAS8KnfsdbhWM4F0rzPDwAPeJ/7A4uAVkAvYA2Q6ne8tRzL8cBxwAdAYdD0uDoWINWL8RigpRd7f7/jqkf8ZwDDgCVB0/4ATPA+Twj8ncX6CygAhnmfWwNfen9PcXc83nkp2/vcAvjUO0+9AFzhTX8EuLmx+7ISSwOo6q6gr1lAoAfEOOBpdT4BckWkoNkDrAdVfUdVK7yvnwBdvc/jgOdUtVxVvwJWA8P9iLGuVHW5qq4MMSvejmU4sFpV16rqQeA53DHEBVWdDWyvMXkc8E/v8z+BC5s1qAZS1c2q+rn3eTewHOhCHB6Pd17a431t4b0UOAt40ZselWOxxNJAInKviGwErgZ+7U3uAmwMWqzYmxYvrsOVuCD+jyVYvB1LvMVbFx1VdTO4kzXQwed46k1EegJDcVf6cXk8IpIqIkXAVuBdXMl4R9DFZVT+1iyxhCEiM0RkSYjXOABVvVNVuwGTgVsDq4XYlO/9uWs7Fm+ZO4EK3PFAHB9LqNVCTPP9WCKIt3gTnohkA/8G/rtGjUVcUdXDqjoEVzMxHFd9fNRijd1PWmM3kKhU9ew6LjoFeB34DS7bdwua1xUoiXJo9VbbsYjIeOB8YJR6Fa3E6bGEEZPHEkG8xVsXpSJSoKqbverhrX4HVFci0gKXVCar6kve5Lg9HgBV3SEiH+DaWHJFJM0rtUTlb81KLA0gIn2Dvl4ArPA+Twf+w+sddhKwM1BcjlUich7wC+ACVd0XNGs6cIWItBKRXkBfYJ4fMUZBvB3LfKCv11unJXAF7hji2XRgvPd5PPCKj7HUmYgI8DiwXFUfCpoVd8cjIvmBXp8ikgGcjWszmglc6i0WnWPxu6dCPL5wVy9LgMXAq0AXPdLr4m+4essvCOqZFKsvXEP2RqDIez0SNO9O71hWAqP9jrUOx3IR7mq/HCgF3o7jYxmD64G0BrjT73jqGftUYDNwyPt9XA+0B94DVnnv7fyOs47Hchquamhx0P/ImHg8HuAEYKF3LEuAX3vTj8FdaK0G/gW0auy+bEgXY4wxUWVVYcYYY6LKEosxxpiossRijDEmqiyxGGOMiSpLLMYYY6LKEosxxpiossRijDEmqiyxGBNjRORb3vNx0kUky3t2xkC/4zKmruwGSWNikIjcA6QDGUCxqt7vc0jG1JklFmNikDdG2HzgAHCKqh72OSRj6syqwoyJTe2AbNxTC9N9jsWYerESizExSESm454c2QsoUNVba1nFmJhhz2MxJsaIyH8AFao6RURSgbkicpaqvu93bMbUhZVYjDHGRJW1sRhjjIkqSyzGGGOiyhKLMcaYqLLEYowxJqossRhjjIkqSyzGGGOiyhKLMcaYqPr/sodLIB5Ub4YAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Iteration, E = 9 -15.95213955965909\n", "Final eigenvalue E = -15.95213955965909\n", "Iterations = 9 , max = 100\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAEWCAYAAABFSLFOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAar0lEQVR4nO3dfbzt9Zz38de7eyo1dNx1jxOO3JRzpXggxKNCuRsqIZMaXOExwuQyY5qMaWSYuVxXKMa4TeLhykGmDFGornNI6dTVOCU6hU5RJEp8rj9+v+2s9tk3q9337L32Pq/n47EeZ6/f+v5+67O+Z+/1Xt/v72alqpAkqZWN5roASdLCYrBIkpoyWCRJTRkskqSmDBZJUlMGiySpKYNFC0KS5ye5NsmtSfaYxed9aZJzZuv55kqSa5Ls12hbH03yDy22pdFksMwz/R/4b/s30J/1f6RbTbPO4f16v0pyUZIdpml/fJJPDlHLKUneP3B/0yS/mWTZ3sO8vnvgn4Fjqmqrqrp4fTxBkl2SVJJNxpZV1aeq6lnr4/mmqOOaJHck2W7c8u/39e0yxDbWeS3zQcuAu4d1HJbkx/3v9plJ7jvXNY0Sg2V+em5VbQU8DtgDeOtkDfvQ+XfgaGBb4Bjgd43qOA946sD9pcBPgKeMWwbw3UbPOZmdgZXr+TlGyY+AQ8fuJHk0cK+5K2fDkeRRwCnAy4AHALcB759ypQ2MwTKPVdXPgLPpAmbSZsCdwI+q6o9VtbyqbmxUwjeBRw58cn4ycDqw5bhlF1TV7wGSfLYfad2S5Lz+j5Qke/fLNx7beD+9dWn/80ZJjktyVZKbkpyR5L5JNk9yK7AxcEmSq/r2leRhA9v60/RLkn2TrE5ybJIbkvw0ySsH2t4ryXv6T6S3JPlWknvRBSnAzf2IcZ8kRyT51sC6T0yyvF9veZInDjz2jSTvSPLtJL9Ocs5YPyXZIskn+9d2c7/uA6bo+08ALx+4/wrg44MNkjw7ycX9SPXaJMcPPLzOa+nXOSrJFX19lyfZc2CdxyW5tH9tn0myxcBzPacfMd2c5DtJHjPw2B5Jvtdv8zPAFkwiyUOTfL3vhxuTfCrJtv1jnwB2Ar7Y1/yWCdb/6yQXjo3EkrwmycrBWht4KfDFqjqvqm4F/hZ4QZKtGz7HvGawzGPpprQOAFZN0ewO4PvAGUn+rOXzV9Vq4Md04QHdSOV84Dvjlp03sNpXgMXA/YHvAZ/qt3Uh8Bvg6QNtDwNO639+PfA8uhHSg4FfAidX1e396A3gsVX10CHLfyCwDbA9cCRw8kD//DPweOCJwH2BtwB/ZO1IbNt+yu2CwQ320yFfBt4H3A94L/DlJPcb95pe2b/+zYA39ctf0dezY7/uq4HfTlH/hcB9kjyyD+OXAOOnL39DFz7bAs8GXpPkef1j67yWJH8OHN+vcx/gIOCmge29GNgf2BV4DHBE/7r3BD4C/GVf+ynAsj70NwPOpAvC+wKfBV44xesKcCLd//Ej+/44HqCqXkY3In5uX/NJE6z/brrf+b9Jshj4R+DwqlpnlJ5kpz4IJ7sdNkmNjwIuGbtTVVf1z7nbFK9rg2KwzE9nJvk1cC1wA/B3U7T9X3R/BJ8G/nPszTPJO5O8p0Et3wSekmQjYC+6N7zzB5Y9qW8DQFV9pKp+XVW3071hPDbJNv3Dn6af3uk//R3YL4PuTettVbV6YN0XZeb7CH4PnFBVv6+qs4BbgYf3Nf8F8Iaquq6q/lBV3+mfczrPBn5YVZ+oqjur6tPA/wOeO9Dm36vqv6rqt8AZrB1t/p7uTflh/XN+t6p+Nc3zjY1antk/z3WDD1bVN6rqB/1I9VK6vnzqupv5k1cBJ/Wj2qqqVVX144HH31dV11fVL4AvDtR+FHBKVV3U1/4x4HZg7/62KfCvfV9/Dlg+WQH9c361/8Cwhi6cp6p5/Pp/7Pvk9cCy/vVMuM+tqn5SVdtOcTttovWArYBbxi27BXDE0jNY5qfnVdXWwL7AI4DtJmqUZEu6T+Mn9Z/uvsracHki8J8NajmP7tPvo4Grq+o24FsDy+4FXNTXs3GSf+qns34FXNNvY6z+0+imFDYHXgB8b+CNbWfg/4x9mgSuAP5AN8c9EzdV1Z0D92+je8PYjm6q5qoZbPPBdCO4QT+mGxWN+dkEzwldSJwNnJ7k+iQnJdl0muf7BN0I6AjGTYMBJHlCknOTrElyC90oaMLfld6OTP26J6t9Z+DYwU/7/bYe3N+uq7te7XZ8Hw3WfP8kpye5rv8d+eQ0Na+jqq4BzgV2AU6+O+sO6Va6Ed2g+wC/Xg/PNS8ZLPNYVX0T+Cjd1M1ENqLb93Bn3/44YAXdqOLewH80KOM84LF0n9bP75etpHtjeTawfGAa4jDgYGA/ummfXfrl6eu7nO5N5wDuOg0G3ejsgHGfKLeoqrt8Sh9wG91rHPPAIV/PjXQHN0w0pTbdpcCvp3uTHbQT40YSE+k/zf99VS2hC/3ncNd9KBOt82O6nfgHAp+foMlpdJ/ad6yqbYAP0vc1E7+Wa5n4dU/nWuCd4/5v7t2P2H4KbJ8kA+13mmJbJ/a1Paaq7gMcPlDzZHXfRZIDgX2Ar9FNjU3Wbqd+X81kt5dOsupKut/5se08BNgc+K/pattQGCzz378Cz0yyzg78qvo1XXi8P8kD+vnur9O9edxBN0UxmY36Hcpjt80nalRVq4CfA2+gD5b+0+lF/bLB/Stb002R3ET3pv+PE2zyNLppjKfQzceP+SDwziQ7AyRZlOTgKer/PnBYP0ranyGnU/qplI8A703y4H79ffrXv4ZuX8tDJln9LGC3dIeibpLkJcAS4EvTPW+SpyV5dL+/5Fd0U2N/GKLkI4GnV9VvJnhsa+AXVfW7JHvRhfWYiV7Lh4E3JXl8Og8b6+9pfAh4dT9CSpIt0x04sDVwAd0Hm9f3ffICuinTyWxNNyK4Ocn2wJvHPf5zJu9/0h0M8W9003qvAJ7bB806+qmwraa4fWqSp/lUv90n97MCJwCf7//ehMEy7/Xz0B+nOzJlIofT/TFeQvfJ8qV0O6ZD9wY6mUPpdh6P3aaaIjkPWAR8e2DZ+XQ7qAeD5eN0I5LrgMvpRk7jfZpuiu/r445e+590n77P6fcvXQg8YYqa3kC3b+Nmutd85hRtx3sT8AO6fQG/AN4FbNRP870T+HY/5XOXc3Oq6ia6kcaxdOH5FuA5Qx6F90Dgc3ShcgXdfqlpzyWqqquqasUkD78WOKHvr7fT7dMZW2+d11JVn+2XnUY3rXMm3Q736WpYQbef5X/THVSxin7HflXdQTeteUT/2EuYeHQ15u+BPen2WXx5grYn0u2YvznJm8avDJwKfKGqzur/P44EPjzuAIp7pKpW0k0rfopuH+fWdH2tXsov+pIkNeSIRZLU1MgES5KPpDtZ7bJJHk+S9yVZle4krT0naidJmlsjEyx0RzftP8XjB9CdWLeY7vIkH5iFmiRJd9PIBEtVnUe3o3QyBwMf70/cuhDYNsmDZqc6SdKw5tOVTbenO6ppzOp+2U/HN0xyNN2ohi233PLxj3jEI2alQElaKL773e/eWFWLZrLufAqWTLBswkPaqupUusMOWbp0aa1YMdnRmJKkiSSZ9AoJ0xmZqbAhrKY7m3vMDnRnOkuSRsh8CpZlwMv7o8P2Bm6pqnWmwSRJc2tkpsKSjJ1xvV2S1XRX7N0UoKo+SHe5jAPpzuq9je7S45KkETMywVJVh07zeAH/fZbKkSTN0HyaCpMkzQMGiySpKYNFktSUwSJJaspgkSQ1ZbBIkpoyWCRJTRkskqSmDBZJUlMGiySpKYNFktSUwSJJaspgkSQ1ZbBIkpoyWCRJTRkskqSmDBZJUlMGiySpKYNFktSUwSJJaspgkSQ1ZbBIkpoyWCRJTRkskqSmDBZJUlMGiySpKYNFktSUwSJJaspgkSQ1ZbBIkpoyWCRJTRkskqSmDBZJUlMGiySpKYNFktTUSAVLkv2TXJlkVZLjJnh8pyTnJrk4yaVJDpyLOiVJkxuZYEmyMXAycACwBDg0yZJxzf4GOKOq9gAOAd4/u1VKkqYzMsEC7AWsqqqrq+oO4HTg4HFtCrhP//M2wPWzWJ8kaQijFCzbA9cO3F/dLxt0PHB4ktXAWcDrJtpQkqOTrEiyYs2aNeujVknSJEYpWDLBshp3/1Dgo1W1A3Ag8Ikk67yGqjq1qpZW1dJFixath1IlSZMZpWBZDew4cH8H1p3qOhI4A6CqLgC2ALableokSUMZpWBZDixOsmuSzeh2zi8b1+YnwDMAkjySLlic65KkETIywVJVdwLHAGcDV9Ad/bUyyQlJDuqbHQscleQS4NPAEVU1frpMkjSHNpnrAgZV1Vl0O+UHl7194OfLgSfNdl2SpOGNzIhFkrQwGCySpKYMFklSUwaLJKkpg0WS1JTBIklqymCRJDVlsEiSmjJYJElNGSySpKYMFklSUwaLJKkpg0WS1JTBIklqymCRJDVlsEiSmjJYJElNGSySpKYMFklSUwaLJKkpg0WS1JTBIklqymCRJDVlsEiSmjJYJElNGSySpKYMFklSUwaLJKkpg0WS1JTBIklqymCRJDVlsEiSmjJYJElNGSySpKZGKliS7J/kyiSrkhw3SZsXJ7k8ycokp812jZKkqW0y1wWMSbIxcDLwTGA1sDzJsqq6fKDNYuCtwJOq6pdJ7j831UqSJjNKI5a9gFVVdXVV3QGcDhw8rs1RwMlV9UuAqrphlmuUJE1jlIJle+Dagfur+2WDdgN2S/LtJBcm2X+iDSU5OsmKJCvWrFmznsqVJE1klIIlEyyrcfc3ARYD+wKHAh9Osu06K1WdWlVLq2rpokWLmhcqSZrcKAXLamDHgfs7ANdP0OYLVfX7qvoRcCVd0EiSRsQoBctyYHGSXZNsBhwCLBvX5kzgaQBJtqObGrt6VquUJE1pZIKlqu4EjgHOBq4AzqiqlUlOSHJQ3+xs4KYklwPnAm+uqpvmpmJJ0kRSNX43xsKydOnSWrFixVyXIUnzSpLvVtXSmaw7o/NYkuwGPAq4P90O9jXAZVX1w5lsT5K0cAwdLEkeCbwa+HPgAWOL+3+rb/Nz4AzglKq6omGdkqR5YtpgSfJQ4F3A84HfAucDFwBXATfRhct9gYcBewOvAl6X5PPAX1eVO9claQMyzIjlcuAHwBHA56vqN1M1TrIl8CLg9f26W9zDGiVJ88gwwfLiqvrCsBvsg+djwMeSjL8kiyRpgZv2cOO7Eyot15UkzU93+zyWJC9eH4VIkhaGmZwgeVqSv2xeiSRpQZhJsHwUeH+St030YJJ9kpx/j6qSJM1bd/sEyap6VZI1wDuSbFdVfwWQ5OHAiXTfofLbtmVKkuaLGZ15X1Vv7U+GfE+SRcCtwF/QnSh5CnBCuxIlSfPJPflq4g8BzwYOowuU04G/9YRISdqwzeSosE2TvIHuzPunA9+jC5ZNues3QEqSNkAz2Xn/Q+C9dBeePKi/+uVhwEHAWUm2alifJGmemUmwbAwcBTy2qr4MUFVn0AXL3sC5/ZdwSZI2QDPZx7K4qn43fmFVnZNkP+DLwHfovt1RkrSBudsjlolCZeCxi4An44UnJWmD1fyrifvvYXlS6+1KkuaHaYMlyTPu7kar6tp+3f1mUpQkaf4aZsTyH0m+nuQ5STaernF/OPLzk3wTOOuelyhJmk+G2Xm/B93hxcuAG5N8Ffi/dOex/IK13yC5mO6osKf3988GHrceapYkjbBpg6WqLgOelWQf4LV01wI7dOzhgaYBfgV8HvhAVS1vXKskaR4Y+nDjqroAuCDJ/YCHAkuARXThsga4DLi4qv64PgqVJM0PMzmP5QbgGuBi4Pv97eKquq5hXZKkeWomwfIB4LHAM4EX0E+HJbkW+CTwT1V1a7MKJUnzyky+j+WYsZ+TPIRuB/3jgP2B/wG8IMmTquqXzaqUJM0b9+gEyaq6uqo+X1Vvr6q96MLlwcDfNalOkjTvzOSy+Tcm+VqS9yR5WZLdx85vqapz6L7k66gkm7UuVpI0+mayj+WbdPtY9qU7xLiAO5KspNuRfzPdtcKeBJzbpkxJ0nwxk30sLwTov3flceNuLwU275t+NcnVrD1q7MQmFUuSRtqMv5q4P/LrW/0NgH5K7AXAZ4CvApsBTwNeCBgskrQBuNvB0k95XQxcQn8eS1WtAaiqP/QXrSzgtVX1o36d7duVLEkaZTMZsfyWblRyGGvPYfkZcD1wP2Bn4JyxUAHw5ElJ2nDMZB/L0iQbAY+g26/y2P62W7+9TwJvbFmkJGn+mNE+lv56YJf3t9OaViRJmteaf4PkPZFk/yRXJlmV5Lgp2r0oSSVZOpv1SZKmNzLB0h9RdjJwAN2Vkw9NsmSCdlsDrwcumt0KJUnDGJlgAfYCVvWXibkDOJ3uu1/GewdwEvC72SxOkjScUQqW7YFrB+6v7pf9SZI9gB2r6ktTbSjJ0UlWJFmxZs2a9pVKkiY1SsGSCZb96Rsq+yPR/gU4droNVdWpVbW0qpYuWrSoYYmSpOmMUrCsBnYcuL8D3bkxY7YGdge+keQaYG9gmTvwJWm0jFKwLAcWJ9m1vzLyIcCysQer6paq2q6qdqmqXYALgYOqasXclCtJmsjIBEtV3QkcA5wNXAGcUVUrk5yQ5KC5rU6SNKwZX4Ryfaiqs4Czxi17+yRt952NmiRJd8/IjFgkSQuDwSJJaspgkSQ1ZbBIkpoyWCRJTRkskqSmDBZJUlMGiySpKYNFktSUwSJJaspgkSQ1ZbBIkpoyWCRJTRkskqSmDBZJUlMGiySpKYNFktSUwSJJaspgkSQ1ZbBIkpoyWCRJTRkskqSmDBZJUlMGiySpKYNFktSUwSJJaspgkSQ1ZbBIkpoyWCRJTRkskqSmDBZJUlMGiySpKYNFktSUwSJJamqkgiXJ/kmuTLIqyXETPP7GJJcnuTTJ15LsPBd1SpImNzLBkmRj4GTgAGAJcGiSJeOaXQwsrarHAJ8DTprdKiVJ0xmZYAH2AlZV1dVVdQdwOnDwYIOqOreqbuvvXgjsMMs1SpKmMUrBsj1w7cD91f2yyRwJfGWiB5IcnWRFkhVr1qxpWKIkaTqjFCyZYFlN2DA5HFgKvHuix6vq1KpaWlVLFy1a1LBESdJ0NpnrAgasBnYcuL8DcP34Rkn2A94GPLWqbp+l2iRJQxqlEctyYHGSXZNsBhwCLBtskGQP4BTgoKq6YQ5qlCRNY2SCparuBI4BzgauAM6oqpVJTkhyUN/s3cBWwGeTfD/Jskk2J0maI6M0FUZVnQWcNW7Z2wd+3m/Wi5Ik3S0jM2KRJC0MBoskqSmDRZLUlMEiSWrKYJEkNWWwSJKaMlgkSU0ZLJKkpgwWSVJTBoskqSmDRZLUlMEiSWrKYJEkNWWwSJKaMlgkSU0ZLJKkpgwWSVJTBoskqSmDRZLUlMEiSWrKYJEkNWWwSJKaMlgkSU0ZLJKkpgwWSVJTBoskqSmDRZLUlMEiSWrKYJEkNWWwSJKaMlgkSU0ZLJKkpgwWSVJTBoskqSmDRZLU1EgFS5L9k1yZZFWS4yZ4fPMkn+kfvyjJLrNfpSRpKiMTLEk2Bk4GDgCWAIcmWTKu2ZHAL6vqYcC/AO+a3SolSdMZmWAB9gJWVdXVVXUHcDpw8Lg2BwMf63/+HPCMJJnFGiVJ09hkrgsYsD1w7cD91cATJmtTVXcmuQW4H3DjYKMkRwNH93dvT3LZeql4/tmOcX21AbMv1rIv1rIv1nr4TFccpWCZaORRM2hDVZ0KnAqQZEVVLb3n5c1/9sVa9sVa9sVa9sVaSVbMdN1RmgpbDew4cH8H4PrJ2iTZBNgG+MWsVCdJGsooBctyYHGSXZNsBhwCLBvXZhnwiv7nFwFfr6p1RiySpLkzMlNh/T6TY4CzgY2Bj1TVyiQnACuqahnwb8AnkqyiG6kcMsSmT11vRc8/9sVa9sVa9sVa9sVaM+6L+IFfktTSKE2FSZIWAINFktTUggkWLwez1hB98cYklye5NMnXkuw8F3XOhun6YqDdi5JUkgV7qOkwfZHkxf3vxsokp812jbNliL+RnZKcm+Ti/u/kwLmoc31L8pEkN0x2rl867+v76dIkew614aqa9ze6nf1XAQ8BNgMuAZaMa/Na4IP9z4cAn5nruuewL54G3Lv/+TUbcl/07bYGzgMuBJbOdd1z+HuxGLgY+LP+/v3nuu457ItTgdf0Py8BrpnrutdTXzwF2BO4bJLHDwS+QncO4d7ARcNsd6GMWLwczFrT9kVVnVtVt/V3L6Q7Z2ghGub3AuAdwEnA72azuFk2TF8cBZxcVb8EqKobZrnG2TJMXxRwn/7nbVj3nLoFoarOY+pzAQ8GPl6dC4Ftkzxouu0ulGCZ6HIw20/WpqruBMYuB7PQDNMXg46k+0SyEE3bF0n2AHasqi/NZmFzYJjfi92A3ZJ8O8mFSfaftepm1zB9cTxweJLVwFnA62antJFzd99PgBE6j+UeanY5mAVg6NeZ5HBgKfDU9VrR3JmyL5JsRHeV7CNmq6A5NMzvxSZ002H70o1iz0+ye1XdvJ5rm23D9MWhwEer6j1J9qE7f273qvrj+i9vpMzofXOhjFi8HMxaw/QFSfYD3gYcVFW3z1Jts226vtga2B34RpJr6OaQly3QHfjD/o18oap+X1U/Aq6kC5qFZpi+OBI4A6CqLgC2oLtA5YZmqPeT8RZKsHg5mLWm7Yt++ucUulBZqPPoME1fVNUtVbVdVe1SVbvQ7W86qKpmfPG9ETbM38iZdAd2kGQ7uqmxq2e1ytkxTF/8BHgGQJJH0gXLmlmtcjQsA17eHx22N3BLVf10upUWxFRYrb/Lwcw7Q/bFu4GtgM/2xy/8pKoOmrOi15Mh+2KDMGRfnA08K8nlwB+AN1fVTXNX9foxZF8cC3woyV/RTf0csRA/iCb5NN3U53b9/qS/AzYFqKoP0u1fOhBYBdwGvHKo7S7AvpIkzaGFMhUmSRoRBoskqSmDRZLUlMEiSWrKYJEkNWWwSJKaMlgkSU0ZLNIsS/Lf+u+22CLJlv13n+w+13VJrXiCpDQHkvwD3WVC7gWsrqoT57gkqRmDRZoD/TWqltN9B8wTq+oPc1yS1IxTYdLcuC/d9dq2phu5SAuGIxZpDiRZRvfNhbsCD6qqY+a4JKmZBXF1Y2k+SfJy4M6qOi3JxsB3kjy9qr4+17VJLThikSQ15T4WSVJTBoskqSmDRZLUlMEiSWrKYJEkNWWwSJKaMlgkSU39f2u2dfbXp7OqAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "\"\"\" From \"COMPUTATIONAL PHYSICS\" & \"COMPUTER PROBLEMS in PHYSICS\"\n", " by RH Landau, MJ Paez, and CC Bordeianu (deceased).\n", " Copyright R Landau, Oregon State Unv, MJ Paez, Univ Antioquia, \n", " C Bordeianu, Univ Bucharest, 2021. \n", " Please respect copyright & acknowledge our work.\"\"\"\n", "\n", "#QuantumEigen.ipynb: solution of eigenvalue problem for a potential well\n", "\n", "%matplotlib inline\n", "\n", " \n", " \n", "from numpy import *\n", "import numpy as np, matplotlib.pyplot as plt\n", "# from rk4Algor import rk4Algor\n", "\n", "eps = 1e-1; Nsteps = 501; h=0.04; Nmax = 100 \n", "E = -17.; Emax = 1.1*E; Emin = E/1.1 # Init E & limits\n", "# m/(hbar*c)**2 = 940MeV/(197.33MeV-fm)**2 = 0.4829 \n", " \n", "def f(x, y): # RHS for ODE\n", " global E\n", " F = zeros((2), float)\n", " F[0] = y[1]\n", " F[1] = -(0.4829)*(E-V(x))*y[0]\n", " return F\n", "\n", "def V(x): # Potential \n", " if (abs(x) < 10.): return (-16.0) \n", " else: return (0.)\n", "\t\t\t\t\n", "def rk4Algor(t, h, N, y, f):\n", " k1=np.zeros(N); k2=np.zeros(N); k3=np.zeros(N); k4=np.zeros(N);\n", " k1 = h*f(t,y) \n", " k2 = h*f(t+h/2.,y+k1/2.)\n", " k3= h*f(t+h/2.,y+k2/2.)\n", " k4= h*f(t+h,y+k3)\n", " y=y+(k1+2*(k2+k3)+k4)/6.\n", " return y \n", " \n", "def diff(h,E): # Change in log deriv\n", " y = zeros((2),float)\n", " i_match = Nsteps//3 # Matching radius\n", " nL = i_match + 1 \n", " y[0] = 1.E-15; # Initial left wf\n", " y[1] = y[0]*sqrt(-E*0.4829) \n", " for ix in range(0,nL + 1):\n", " x = h * (ix -Nsteps/2)\n", " y = rk4Algor(x, h, 2, y, f)\n", " left = y[1]/y[0] # Log derivative\n", " y[0] = 1.E-15; # Slope for even; reverse if odd\n", " y[1] = -y[0]*sqrt(-E*0.4829) # Initialize R wf\n", " for ix in range( Nsteps,nL+1,-1):\n", " x = h*(ix+1-Nsteps/2)\n", " y = rk4Algor(x, -h, 2, y, f)\n", " right = y[1]/y[0] # Log derivative\n", " return( (left - right)/(left + right) )\n", "\t\t\n", "def plot(h): # Repeat integrations for plot\n", " global xL,xR,Rwf,Lwf\n", " x = 0. # matching point \n", " Lwf = [] # left wave function\n", " Rwf = [] # Right wave function\n", " xR = [] # x for right wf\n", " xL = [] \n", " Nsteps = 1501 # Integration steps\n", " y = zeros((2),float)\n", " yL = zeros((2,505),float) \n", " i_match = 500 # Matching radius\n", " nL = i_match + 1; \n", " print('nL',nL)\n", " y[0] = 1.E-40 # Initial left wf\n", " y[1] = -sqrt(-E*0.4829) *y[0]\n", " for ix in range(0,nL+1): \n", " yL[0][ix] = y[0]\n", " yL[1][ix] = y[1]\n", " x = h * (ix -Nsteps/2)\n", " y = rk4Algor(x, h, 2, y, f)\n", " y[0] = -1.E-15 # - slope: even; reverse if odd\n", " y[1] = -sqrt(-E*0.4829)*y[0]\n", " for ix in range(Nsteps -1,nL + 2,-1): # Right WF\n", " x = h * (ix + 1 -Nsteps/2) # Integrate in\n", " y = rk4Algor(x, -h, 2, y, f) \n", " xR.append(x)\n", " Rwf.append (y[0] )\n", " x = x-h \n", " normL = y[0]/yL[0][nL]\n", " for ix in range(0,nL+1): # Normalize L wf & derivative\n", " x = h * (ix-Nsteps/2 + 1) \n", " y[0] = yL[0][ix]*normL \n", " y[1] = yL[1][ix]*normL\n", " xL.append(x)\n", " Lwf.append(y[0]) # Factor for scale \n", "fig = plt.figure() \n", "ax = fig.add_subplot(111)\n", "ax.grid() #j +=1\n", "\n", "for count in range(0,Nmax): # Main program\n", " E = (Emax + Emin)/2. # Bisec E range\n", " Diff = diff(h,E)\n", " Etemp = E\n", " E = Emax\n", " diffMax = diff(h,E)\n", " E = Etemp\n", " if (diffMax*Diff > 0): Emax = E # Bisection algor\n", " else: Emin = E\n", " print(\"Iteration, E =\", count, E)\n", " if ( abs(Diff) < eps ): break\n", " if count >3: \n", " fig.clear() #erases previous figure \n", " plot(h)\n", " plt.plot(xL,Lwf)\n", " plt.plot(xR,Rwf)\n", " plt.text(3,-200,'Energy= %10.4f'%(E),fontsize=14)\n", " plt.pause(0.8) # Pause to delay figures\n", " plt.xlabel('x')\n", " plt.ylabel('$\\psi(x) $',fontsize=18)\n", " plt.title('R & L Wavefunctions Matched at x = 0')\n", " \n", "print(\"Final eigenvalue E =\", E)\n", "print(\"Iterations = \",count,\", max = \", Nmax)\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.13" } }, "nbformat": 4, "nbformat_minor": 2 }