{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# Symbolic calc with Sympy\n", "# PHY 341, Wang\n", "from sympy import *\n", "init_printing()" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Example 2.2\n", "# declare variables\n", "m, n = symbols(\"m n\", integer=True,positive=True)\n", "x, a = symbols(\"x a\", real=True, positive=True)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# stationary state\n", "def psi(n,x):\n", " return sqrt(2/a)*sin(n*pi*x/a)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAALsAAAA/BAMAAABKn6SXAAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAIpnNuzKriUQQ7912ZlSugimNAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAGQ0lEQVRYCa1XXYxTVRCe29vb9rbbHwzhgQd7BSOJolt/IoZEtmpMlCy2MSGSVbP1j0RDpMS4mohpY3zSh+0bCg9dlECICA0aCGCk4A+aoBSNiUZx1ycSEtNlBVz+nTNzzu0927tuy+556Mx8M/Od6TnnnjsXQB+BOtnGwzrctWUc8E35mdHRx329XYBf+8VaWwmNT/g5u8KCRZ/wzQ0C0xkfX3eQvdsnfjFjfXkfX5fQmUpbgnWeIOufNk/3QMxpyzELBIX/bfN0DyTOtuU0UwSFrpCw7zzGtoqLf6e0TuR9FHTT4kfd4M9ZC10luSU/qK/fisGyGzqz0sxjjFGAkkqyrnFSiOUiGNBI4rUF+nSat82IOgiFKhAdka7QOCuSfhK+lQ4JFzRzJiMg2KIXIaQ2oepwCtPbZ2HVh16OkHR7sf/RI2KJQzUI8XEEGBWrJTBenK3WZIMB/u2SHj7mtJhcE8hJMkn/zp/HvPTWK9sL9tKj5cjBH19LAUSG5m1cZ218zFsA6gQz1otROAbrJADOSSnppeWKUAGeTkE/JO9K47ThYD9Uf4Uzwh04RKOAKsHM28QoHHvpFydWT9P09PuxmHysZmPGmmgW0kUYlclKEJwnK10QIuCIXxyBSZZy7ZPXcYxLSIhQQZzcZNEsCKtSLUOzAi8K3TMELKfkE/mTcpoXpTZt9QksIJk1HYprAuANuEtlK4mwnDIo+BJFmMeuZE2GTEtP1TtmkeJeANgPkbNioTxrDwhfFRhAjziATwF8TxaksyzVwdw09LIEWODWfgQwnJf0O8CehFDmZi0GYAckxsOEJfClZOxddJusuunISFm9A3doqcEsDORhGfBVa5yH+DiYRf4rbiDC4cxJMi08iSbuX4adJ+oySNJPQjovESGslefqxpKhciS3bQzNSA16HEgsTQlfayBsrBoj277cglHrZdR9an/T6bXYjozrWlRO1RrkSwHgLc3ftbGn4k3pU//TlHecdcDr7l53CSnVnSw5wlTztdm7p8+VvTnuUg07BBsFS/0fb1jn+voxb6xL38dVz1/95OzK96cP86sW9lx3J/RW0bne2/DE2oqsWfegs1GRXlyLfC/al5iqx699u6FZSnVPmqLnZ97juGHVl96+d3bnpVWNRg9q7c2drYhZadrWuvT2XHSwoi5/eiiNzapoN1k/93sUnq4pbXYyN+bNP6yeUezc5mS4VzCx5dSJCVLLEN54i/8k1nLHz7GwMRU9rAjJ0VuWfn5brYAtmls6iwCxgtQ1IS9CD+auNmElNT3Rx2sQznpipWo4AMGCNDRha5YwxDVgLDlSFzp2sFgYDaJHjoTPHuD7dRp6mdwShnjXvpeidge14YJ0EX3SgYh6KUpciEGnY/rIBQD7UwjUOL2aYcl9ThXp5c0sYdEeJ9bvexWCzx1fh3976VF4f9fpu3+//HffA6FLJ5cX4Jk/viKYE+L4iRa4CDG5KG4jTtVje2rxJ5xk5/ZYXHjBVTCYgoFUTwYm/rpqX4FqGR7CDbezEGOYUgL43RDNCieb6ivFpVctM/u5PXaQfgSSDdFYX4NDUMG2L+pAEXckshPCDFOCOQ5QbfXQcbUWRI+OiFY9t8dUfRaSdQujl8H9SNRMnR6PlMWGH+5PMUz0ometjsEjiRSZBvaENNTWJrSt5fbYpU9MnDr1C9yDCdFicXeYNnxT3wTDRDOIi27W43vRR+MTKYkeO8q4djC5PY45YUAXVk9zL8OUwIH6rURvVOzR1a2STuTx5Kx7fcE3krZUYYXo4yMQyEgHC26PFT1sxz4TlwebzwupwSOi+kQdQg2CKWGllg2QrjNA9PAlbMbpPYPa43AmDKaoHgbKsIE/nD6DGBYSzCZ2QjBFsMiy8dhrw5TVMv2m4y9pXny+h8oAazdYuW3P5vaB8fwbqXcvPYExP0BPAxvobR+8ffxNgikvrr5nFUtcnkymV+iNSrM2JdOW53xu6MXXqT5uZ3Nu6Es6N1rVOkEB9Xy1BXQB4F0xdfSMEBKfuudT4zqx1UXpjf2CDdXxeF3d6ulye8ZwirBSvd3VLbLWJyHMJ9+cemR9QmeA4r6tXj9l2X3lGbJndC/M+4VsYVR9QfiFdIY96BtmlBle4+vtHLSZ5z+WT6ltZJ3U/wAAAABJRU5ErkJggg==\n", "text/latex": [ "$$\\frac{2}{a} \\left(\\begin{cases} \\frac{a}{2} & \\text{for}\\: m = n \\\\0 & \\text{otherwise} \\end{cases}\\right)$$" ], "text/plain": [ " ⎛⎧a ⎞\n", " ⎜⎪─ for m = n⎟\n", "2⋅⎜⎨2 ⎟\n", " ⎜⎪ ⎟\n", " ⎝⎩0 otherwise⎠\n", "─────────────────\n", " a " ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#Integrate wrt x from 0 to a, orthogonality check\n", "integrate(psi(n,x)*psi(m,x),(x,0,a))" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIUAAAA/BAMAAAA8pRYcAAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAIruZZhDNMomrdlTvRN05z604AAAACXBIWXMAAA7EAAAOxAGVKw4bAAADuUlEQVRYCc2WTWgTQRTH/9lms5tuPhZPCkLi10WFRk+iIBEEQcH26EFtRJBKK8lBUKQ2K4ioIEYLih/gakWwIip66K2lpYfag4snT21BFIui1dbvzzczzWY3Joo7Cr7DZt6b9/9lZnZn3gBkoUXsKWfNS+T0pA5PSyOQyskz8gVphvZOGgH9gzwj+qkWI3yzVrReLDrDe7pafAntRcfn/9qJfqZ+tTvvY4Qn5vn8XyPAGcAtnyZa+o3K312bkfEn/cb7Vwzt4khJXd3jxF/cOGLSGB6dmnvd3n/FqRoOD7NY7XGU0GViKZJrUlnKsRba6oxjDFATa18ym2RNFn7NGvUYz4BiITGhUormDCL0AeEpll8xHn7P/DoMjV550oqUWIqqfoXxFtE0cyrGwoJbh2F8JUY6kuESYwoNOcRs7lQeFBbcOgw+jkzE4opoDgkL42YLeZ71oHAyy2J1GBgCxguzjJiFVgdnQrc50X1QuB9WPYaSRmcBbUiIuaQKaAbadFctGhQ+EjZdxukHk3YlQ1v+xg6tOubEmwb5X48BPcDdy5UM3qLwnE2sNTuXqu4/cxW2byUtwr80OUhyQE7P1ONi6aVA/gMsEEoXx2kgbVnUb5dbgX8bRwNLXeHsx+z6QRrqLva9S1pkWBJAcvUvlGyc5DtTbjCpCTk9U8feyjMUOnwB/fiKWiito+Zump/1J4szqB3bf37JFpAo+bOFV71POSM8AT1dnR2iMSg1GaxweY0zKJWKUJU11mVUJYrzNJlB/KczsfhnjFZieA8BVvSNS5MXoJzqPU+38dU9GLvft2Pzt3v5PdEv6zpKOLDxGg/zIfG5pCxo3rtdFyv6bD8qi1E00Wk25jB9Z0b9xGrVAlppNY2ECDOIy/DeMZ/xos8YA1QO8Rj4jJdowRnEMlTZlFJ8GLoIu4xWC3HPOETR5+NII2lrNM027KbsfrNvKu6wt/VqqSnCLoPW1PCsqSj6LsOYXr9+A3ZSdsyyRqliEmNLflqEXYaSRtjzbkXRT2R0UAeNg+PbKLvhib2SM0ItavM+91/FNzaAhhzllG2IF/0yAyNAnGZDN6KPZpGqrlIybESzPMwknIHn2FooA+iXF309pyPCxoFOB4fYCgJPkcgRI20MQzF5mEUFY0vvWeaUjRV94NwhrWlwb9MkQt1HzYdfllHvVTRmoS0f3Ha49wQPc4VglMXBfv8XRoN3nwSbCsIfAwq9su9eJ2D7pB1Q6JFF+JXbEwjQVPP0Ocna37jG4KDUKH4ApNEdVJdZ6E8AAAAASUVORK5CYII=\n", "text/latex": [ "$$\\begin{cases} 1 & \\text{for}\\: m = n \\\\0 & \\text{otherwise} \\end{cases}$$" ], "text/plain": [ "⎧1 for m = n\n", "⎨ \n", "⎩0 otherwise" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "simplify(_) # simplify last result (_)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "def wf(x):\n", " return x*(a-x)\n", "\n", "I=integrate((x*(a-x))**2, (x,0,a))\n", "A=sqrt(1/I) # normalization constant" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAACcAAAAyBAMAAAAtlt8uAAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAMom7q0QizZkQVN1mdu8ilFydAAAACXBIWXMAAA7EAAAOxAGVKw4bAAABiklEQVQ4EcWSO0gDQRCGf4157OUMsbHTxFgJBhS1joUWFkJ6wVylNmIqOyUgomCRQxQUIV4ldhELwS594oNgYSGKKJY+ghKNSpwJyd6FC6R04Xa/+XZ2mbsboDoCZTlqChg0UVLLgEQT2qImS5qVZIFVuEIrlB261qQVcexALcFtiBMpPQbuDFFEFuiVcg+YjqKIYSCt12w/g/KDTyBmMHIxR/RgKiNYZghv40D7JUHXmN9VAALErs03wEuOzpQqMkgo3oEzdsDLLmXG+BSSGk5p6QAiBt2ZNlguJJQgHSjrJIeAC52l48vN6yMd18+BNXZw/h7ykuM6vYa45wBYXuTZuU4fRISv/BzQ1aOVpX7yavXxf0eyWUx4bVaT0GbsKerzRANpV4C6rdm1YtzYpcC+XfqQsksl12OXDQ2/Z8ONpnIuPE/9lR+3JjqfENHhq68zlcEGNSW6rZlJHcccj1jlA8QHxZ5qx1S2RAEq/4K+VkumKMERP4C7/G2RWEI2EbQK5q1856T1vj/E7XsJNEE6vgAAAABJRU5ErkJggg==\n", "text/latex": [ "$$\\frac{\\sqrt{30}}{a^{\\frac{5}{2}}}$$" ], "text/plain": [ "√30 \n", "────\n", " 5/2\n", "a " ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "A" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "cn = A*integrate(wf(x)*psi(n,x),(x,0,a))" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATUAAAA/BAMAAABz3ULlAAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAMom7q0QizZkQVN1mdu8ilFydAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAIKklEQVRoBd1Zb6wcVRX/zc7u7Mzu7B/qnw9G6dZg0MSmzxhjYoQuqBhDXrrxgzGiZYJYIYhvE4wGI3YlKi9Q8zaI8oHk7dpEY4z4HkZM8ZGw8YPxX+lSG1CkdQGVaEy7oNa+Unyee+femTM7Mzv7wL5GT97O/Z1zfufMmTv33pm5D5gmX5vmPH8+66Hs3G47m5PFcLx9WRTmV+zfMFMKvNhLcWzC7J78wObZZiMrxvlbFmMGvzsDJ6QotvPX0ERo20YgY+WwX+YtNWQmlcX9hhc5T0xJZP+gy3lv44qPt9fjtkSLe+ANv2UOk8eV+k8xVwJMZJc7jGnMMUXBz8RNZLGiVo9Ud++gKlolr9JAtA56XI3jRHbpNCPmW0zxYcQfep8MoUAfo9+Nv0etG5rfF0KgimWuxnEy+xZG/CRw3ZvfAeR2flVby5GutBZvqQtPqaH9gEHrUKVP+jz2hlarjRvuCIZI6dc7QpdGOhnpKewVeTKf/xU4HRzp4/O4z1MZtvcVkM31yP1TgHtD4z13rJNyK/2exu1GYLdbVhO7WoGeAHQycqWwq50gzGmjcha9truKnO6uPYFXAFqlTor2gDgoqYjafk7XPsaD39JGqj7fRXWkdDFnA5cGOhnpKWx7rLmw+yj8CbW22URJnJDEOeO36niqiz0e4ZHSRSNrq3ZRaOLoILDfieo68nwwBy4NdDLSU9iFc5qL+yRaGdQ6KJz1re7zgVeAq4E9NAYKbWaVtRUHzELQWUWe/qLRUYpONo39XBDyVoEKD6NHtamKzbBXFe0Rau2GUkQja8s1BQwlNxC4PIZ1++NDmjtHIxMx4IlkJAEbN+y8mbOXPPKWWnRYpZ9x7AosNGC8RJikNue3wVE6Ir0ka+MvBN8mclkGbB9ir2fS8Kgmr28iWYRd+DN20zAN2CsDynOiTR0wlAlX+qK2FyTGtQ2/DY5yrS73A131W0lcly8VcSOOS3wI9DS2yWXhYuWNNJQswsbyAHdx9kKTtLtpZPjXSmO410BB9dtSK5IMWBR6laz3rgl5QNVmjAPejbvq/i2A3YG9jnJDuGjVjAsli7Cx1MUPBU2z5UR3/gEcJqNRh3m21kFJzYX93WhGtyN0UVsg8p4WRoGOXhPuUKi0sFTn0BNkux76AySTcTaehSPWz4AtxgNAo+4hanrrVJs5B1etIS8KH5ODsDyg2GcmWRu7p7BPo9wlQqmBfXQLjgjuW3IsQkOZjLOdF+COOZuWW5KFDqWi62zSPXVHsNvCCOvvsgkORgMVOq05DCzqnrpzocU6g4NC2wb8rtfHO0seihuTF0l+QybjbOcszPZ3GLskl1dzvUgnhTsAzY2fYLlOCl16dOnFTYuLV5M516SDFtlvNq92d/dWclqHFr+8Why6h5K6TASrZIyN2/CrTkP4lBiybwovfVfqH16kcXj38St8Z0VNCc3dv7HxL8KFVW0A8pe/OE/zyAstuPbJDmlFek61nfkd1z/GXByqZIyNe46+dt7vFZ/o+N39xc/yOIVtNSUmXTQ9o3ITV4vPeVzNwFPZ/lN4pZmQw6SlJUloBkblAFdLNOlnl6nsU12RKHJXdOZiysPa9DTDb625iP7jiDahfHRCxzT2fm+SHej0LpEoTrQWmPLyEqkx4/0xyxTD7laqszZOcX0qav9ZVJ2qbaq2Pf3UXKm1pUZkO/5valsapF5sL+2epkZkOzbVb1QbLZOTIqvqjeS59DtHbZK0KV2kumtt7eTa2o8EzEomONg1lE3SQdWW5HrZtk312/9Obd8UL0FaLvg8jc6FR//Y1ZXRmJh5LlhPHauruJsPqzeFMA9HCfeUBTMogqK1DXie2WvLNyojFXklTvAckzihNhbMoAiMrr2/+BBLVp2534x6rq0CL8KlLEUMJtTGghkUkbv7PP57lUGollOepyEjROUw7uHQGkexZ72gsGAGqTY9UGSersP6ykx5D4mfEZXbAuMTwwDOCFgwgxT8iMczjJxRqNriq2dG+XpQkXPZjCEhLQwGg8CpkELo+z9l3Zg7F/FNVYzwPXTmfdggIQtmEPh3wIgBV33exxzMUDlFT67Oct06hw9edd2xofM0FgbMPx2eoOAzMrhwyb6d8/QRRXkCsSKfZ+IpFw64wgxv12/8yI7DF9EXqfsXNN4+dE7Td+6j3SB7Bsi96eDrH6vL4Jx5GXqehEHQtPM74qtquhj9ikdfUsbOY57RehbW89h2/KrpIcz7OjxeqctgvIb2ABbqIk/od/V3ntxfwv2XtEIf1CYUs8ShKbYzhTi0YSC/y311tuMc7W9K6dLmhNwDYHG2GsP+/hJGfA0RWxGZ8mq5OSZopTHyc5n8CMFolrW+ArxHY9UWxz5YlvtL8IwRIxwZMCUF3gmVAfk2agPVDSnkSbPdr2nT+4Fzjlb8trbqt0tqf8lsMf9ChynJ0FmHXgWrDaygkUxLsfa6NX2+Z6jfJ/Yntqtkan8JX+JZqqpwbpvAlTaeUaaFOi51vQn/dPUJ0AVJoVmUa098sP2hLl16fynXVVxpzWU/UM0+rpRc4Brg479UeMbmXcg3fCptsljv7kfDLvdVtb+E/Ru8qyw9iaMxW6SJHUspk/tLvvVzyntBmuA/CJP7S341R7wLUpV/0iK/hfE6eo24bcssGc9le7RllcRPtCtu4hbnAk6GzHOf4JVuLbanDzd6kx9ubUHsbAstpiTBwuyfDEnhr8T23szgqR90mdGvgEAvq1kS+UdHFvm/6f90PTOb80Am5fwQvjBD2k/MwDkPFMefCf8BSVGufFwAA/wAAAAASUVORK5CYII=\n", "text/latex": [ "$$\\frac{\\sqrt{30}}{a^{\\frac{5}{2}}} \\left(- \\frac{2 \\left(-1\\right)^{n} \\sqrt{2} a^{\\frac{5}{2}}}{\\pi^{3} n^{3}} + \\frac{2 \\sqrt{2} a^{\\frac{5}{2}}}{\\pi^{3} n^{3}}\\right)$$" ], "text/plain": [ " ⎛ n 5/2 5/2⎞\n", " ⎜ 2⋅(-1) ⋅√2⋅a 2⋅√2⋅a ⎟\n", "√30⋅⎜- ─────────────── + ─────────⎟\n", " ⎜ 3 3 3 3 ⎟\n", " ⎝ π ⋅n π ⋅n ⎠\n", "───────────────────────────────────\n", " 5/2 \n", " a " ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "cn" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAALEAAAAvBAMAAABeazbFAAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAMpndq3ZmzYkiRO+7EFRAb64nAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADoElEQVRYCa2WW2gTQRSGT67N5tZVoRQfbMRLRRADKYIPSpRWK/QhQqtIQGOhiBcw4A0RNQ+tgiAEpYogJGDVB0EqiHhBGh8iiCgBFQQLDYrQt6ZWC6KlnpnZ2cwJiU3ZnIed85//zJfN7GYyAGq0zptRUuvW81fWEdUJtlz1uvWqM2udUZ1w0yjrEEhX71hkdansX4/JsgLA2JkIXhsQHQYjkITA9ovIjMxcbQAWwCXJ9iLyDiA52hAuQMtJA/SBjY0kpyT5UoPJrjSSg1mkhiV5R99KllqNG4DkeBLXO8ZQbDXeQSZrFYvzU0jWPk8DODiMkQH8YS4sXdxpds+BnwADnCPI9t+WoHzyMDAybEzDFkn2J8GO38FqjHR2/ukGyISCKY7Ce/aFwVmqk9t/sC9Uu/UxWv65Jr7A7Al6AHy52u3EOfIEnpECEX9ReX99FTUkayGIi4+BNGmkSkNP+xSGOdIjhOhsn8+jPHeW13rGOmJwbfUG0eAqirGaAm8My86Yd1TtMfK7Sm0qqgiZTsqEj6piW+J+rPp0e0onXWDDV0G9JUea+kzh5qeEosSWyJ5MBhwPKfljO1ue88rEKqk9oRaJYq+9OwQwAi6+3WDjPaPZxcj9hqgx+EidKEaGl8T/bihO9jG/dtwnFlGcfJr4hNyUIF6lWEcKRHHyIPEJ2RMlXqWg90QUJx8iEwjZnSRepdhKCkRx8jfiE3IwDOZxqyIp4azn+OKdiLDQKxUnt5XJqyKRF5FIByvwJ2hj82sHkpUgipMfKC4Aueeqv8xy+4VyihlRdawGmVwh8kQTxcnHiE/u2Z0zveGBbUau7elly4pB3yqiOHmNaDOuhGyPmV43xI3cmXKNivSO6bJEVXxLhM3EN8jOtzNd+KeaNr0lcNnIbbonKVKHbtqYUIUFLaza5hPkxeuqtckUjoRI3SExVlNYI1tUeUfi3WuVmXtjUrhOyeyNTPhIFcBt4hKh5RQZyJtiXH7IuFliCVXswFMz/AXVatOlsskzgZaQJTZSBXhSqRmHuRPH3/WPwCPIJLyDt4a64IuuzdacshjDc2WyZbeO55ijBY8/DxNpR8FNt/TF0NTeFbDPhevQ2rcLluMxI6Pbhnr/8z3VqQvmOWdB9BQmsnB8wfb6G2xRcQrFGVMAPfVPXLDTXmyWPTsBZgNSWB8nCs1Zg/IUgiWPdaIk7AVfSuTaNHiS4lWUpqXxPThTAuANg/a6KHLL13+8efsvq/xnQAAAAABJRU5ErkJggg==\n", "text/latex": [ "$$\\frac{4 \\sqrt{15}}{\\pi^{3} n^{3}} \\left(- \\left(-1\\right)^{n} + 1\\right)$$" ], "text/plain": [ " ⎛ n ⎞\n", "4⋅√15⋅⎝- (-1) + 1⎠\n", "───────────────────\n", " 3 3 \n", " π ⋅n " ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "simplify(_)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# other integrals\n", "$\\int_0^\\infty \\frac{\\sin^2 x}{x^2} dx$" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAA0AAAAlBAMAAABrOn4UAAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAdu+ZiURUMmYiu6sQ3c2sJJ3AAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAk0lEQVQYGWNgEPr///9nBgYW15QK8QIGhjIGSUYgBQQLuA+AKPYNvCCKgTmBH0yfP8AfAGLIMTA5ACmGZQzcYBrEpgAAXQEEH0gwYbLdcpBqzg0M/iAHAV3LpACkmf4wcP8E8S8wcH8H0kDACzW2XgDCvwyhmA0gdDOEYnNgmAlilTAwyAApzsvG1heANB/QWQ8YABXUJZ6Hm5PQAAAAAElFTkSuQmCC\n", "text/latex": [ "$$\\frac{\\pi}{2}$$" ], "text/plain": [ "π\n", "─\n", "2" ] }, "execution_count": 25, "metadata": {}, "output_type": "execute_result" } ], "source": [ "integrate((sin(x)/x)**2, (x,0,oo))" ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAKQAAAA1BAMAAAA9oV2UAAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAIpmJdkTvuxBm3VTNqzIaxHoRAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAEH0lEQVRYCbWXS4wURRjH//PYne7p7WEEE4kJ2cl4Q5E+eFNgjAe9KJuYmMBlR00YTTDsRQaIjz4ZL2QnJFzksJM12ZgJhDVwZ4ImcNtR8RESYI0iCcHsKBCX8Fjr9XV3ddXsDgTq0P19v/9X/+nqrq6uAQa3zIHB2iMqV795xI7Jbt7zx4JE/lgsX0ddWFbE0bT8WfBBB8d2py6XpkS9PBqWfmuQm+Q3QlNv548KOMDyJ7OHRry2lvLE6XkvCygs/3xvXq9w23qezK6HPNuSRHosr1JnwPhimlDu/LgQ8rjYIoJrTdbKWOEN2Waz0TyIdSzuRxU8+FrL9GRnyHPvA50mMutVuvVERTqUljib5ji9bUIwYUkJVY0GFFnOynKyzLTMc6/WopK60xexsKSEVNnr7R82vBl3IA3KcqzCyDvlwnIklEt1EQtLSkg9IYKgWXP+IRSflaXHvJwPkW3HSmFCxBVxVIlSSw94UJpoICMihdVJWTq3gOwyikGsvhiHgJbAr3PNcW7D7SerZKwscQgYncK0vDIu5cKEvZawyp7o7PXhT0mb5JEsZ8qYDjATSwsr7QEJsNQSkt/FulYYV6mILKdbmJ7HJ17ZqLCAHbKKTaUzCAydLNlYxmrubM4osIH9ErK38hXXvAayLLThXNi+6S2bg8EakuwGnjlmiL/uO1ETMNc3NANcm2tc2v8+e9T/GpINuLdtNMUme2LRKrEZN0QbqmycWRbYInPXYrjLYEMNhltmAdfyyuC0YYn7JjIIt/SAvOXFtln+ZxiYgFuylosXl7jGcpV7Y3VgpCzZYmA2i2XHrDKIsvT7hgLbwOfk50Z8c+wH9hGTA1eW/IvEGrff3Ons6XTEKiqY/M05yy+n0ZMbeG7IJ/4Qjyc/5Lz8KD1KS74k7+WwU32IpeDvk/cu8x/ybOuBOYmcm8mrylTPLSZzPS7dUfn6Fz6LlMgygvqysQk525VQfzWiTA87JojRshFDfXX5DdhDtZbzccn8EKP1tBxDfSnYF+JKOV0c57/IcHQZvjGfYqi/ZBeBK4uxRTqaCQXx2/CN+RRD9hHX20lg/UtP64yypXmKin2KEmcFx2sJxkL2uEbexeehTlWmtgYsW6pZChRkWwOtFSv4ooXNGouSbJ3CWQqSZwXZ3dNaFdgaDtrsjtAMy1a0TjIheEvXXFbbANtH2Rtd3BGbrKCbWqhPIbPxpnVfJkzUI/cCPGWYEiz2NKkUIP/lXRS632o4StTG+g/gjYhRQHBSfzrPVqsX8T2+qgRUqJ/loDKz1S1tXWC7c4Kf6srCyso9HH7t9wuLOo+y73g0xj4K3QipgOBI6lam64z8L31Uhs62tfMWuBrKd1dTuXZprQJDX/NvadfoshYotFavuBGurtvUj20wYs55Hv4PyWcxs/XZsTEAAAAASUVORK5CYII=\n", "text/latex": [ "$$\\frac{a^{- \\frac{n}{2} + \\frac{1}{2}}}{2 a} \\Gamma{\\left(\\frac{n}{2} + \\frac{1}{2} \\right)}$$" ], "text/plain": [ " n 1 \n", " - ─ + ─ \n", " 2 2 ⎛n 1⎞\n", "a ⋅Γ⎜─ + ─⎟\n", " ⎝2 2⎠\n", "─────────────────\n", " 2⋅a " ] }, "execution_count": 28, "metadata": {}, "output_type": "execute_result" } ], "source": [ "integrate(x**n*exp(-a*x*x),(x,0,oo))" ] }, { "cell_type": "code", "execution_count": 31, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAAgAAAAOBAMAAADgeEClAAAAIVBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADdcGRXAAAACnRSTlMAMt0Qq5nNdrvvxbMB0AAAAAlwSFlzAAAOxAAADsQBlSsOGwAAABtJREFUCB1jYGBUZmAwCVvMwMDAThbB0rnIHQBAqQoWQ8kakgAAAABJRU5ErkJggg==\n", "text/latex": [ "$$1$$" ], "text/plain": [ "1" ] }, "execution_count": 31, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Chris' contribution\n", "integrate(DiracDelta((a-x)),(x,-2*a,2*a))" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAABYAAAAlBAMAAAC0dAQLAAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAdu+ZiURUMmYiu6sQ3c2sJJ3AAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAu0lEQVQYGWNgYGAQ+v///2cgDQQsrikV4gUQdhmDJCOUCRRYwH0AIszAwL6BF8ZkYE7gh7PPH+APgHHkGJgcYOxlDNxwNkyMJjTQHzBAqfnFSAbcRrAZkdgVrxHiDgg24wQgm9NEUgAoO4MByA6awPUHyHYAsnk0GJgvMDBwTACymf8w8DowMExkALKZFjCcBwZd9O7dn7acd2DohxiszXA+gWEd2wQQ7xcDnwDHZRYQ0+7/Qp6N7rPFGQDJYzcFylU0uQAAAABJRU5ErkJggg==\n", "text/latex": [ "$$\\frac{\\pi}{4 a}$$" ], "text/plain": [ " π \n", "───\n", "4⋅a" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "integrate(x**2/(a**2+x**2)**2,(x,0,oo))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## momentum w.f. in ground state of infinite pot. well\n", "$$ \\phi(k)= \\int_0^a \\sqrt{\\frac2a} \\sin\\left(\\frac{\\pi}{a} x\\right) e^{-ikx} \\, dx$$ " ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "k = symbols(\"k\", real=True)\n", "a = symbols(\"a\", real=True, positive=True)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1MAAADACAMAAADWfMneAAAARVBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAeYafpAAAAFnRSTlMAEFRmq5l2u0Tdic0yha8i72xg31xAOe2oOAAAAAlwSFlzAAAOxAAADsQBlSsOGwAAIABJREFUeAHtXenCs6COdqm2tc5nZ86M93+pE5aEfVPb17b0R0UICUQiCISnaTJ/K/4y6StZ1UDVQFQDfR9NrolVA1UDrgbazo2jmCuFaqBqoGogTwP9egkTjl3TTtPl0oZIulsopca/VgPtspX/tDVjzZengSFmUg10UwPwGe4BZh1Lrb+/0MA06lL7u3GrJznhuQ4+HJ0cGTGuMQVDN9Ws0BMt6+wV2j4wOv3STFMgL+sazBhMsBhk3wYZYgJeszl6CYNcMAGv3tyByGfBgKHPt7+AtBod1sBtXYPDuqZpmb31QBCyqSvOYDDji//SFIH84Yy3g5tGWtIhEl8iZlkD6vNG3yMP3ZvhSyL7seDNk1nn5S5mxulFeF+9g+tF2MqEJZj8Y7/lKeV2Npf2+jSzOBR2iZcn2qeZ4mTUWIdbp+IxZo9OcyTlSFSyvaFjxNxG/hXcDv04safZX5slWdUZbUlk9hbvmyMf3cxGXcf+hhasRWv/gW7qNvHei3dTrAD0LKzSDJKVGgIqgofRln0UipaFbquRARN9GRXrR/p9e9eqi0y91zxJGRK93CnyIDFTyx5Re5+bmT+ry7h07ZPelSROBdrL8MBBCc+skn4k1N+bJTzbtlEJ8BSmpdXa7mWlDyKT5YO1xEF+RC3XQNNFox88XcxqjAZ9FKbExi/Dl1Gx7rW6cHY3T0ksOaHbPEm7JW4UY5V6aTs2Dhig/jMfLDwnqLn/o1dlHdGmmsd2PSl2nxa6RCa499Tl0cxaO7z7O4em6UD72E0tUJTZ9wYUb0goDQ4BtYLd6PHxSA+FRhwJejLqrJ+WJXbbP7EyJe2VuFEMqGi+4o89jCurKXup8XHcsg74BowoU9nUEJuZinD46KSL1vIPrMjt2tw0a13XUBuEIZN8SMt1nufBaru8SKP8ZLp5OrvJaDs+iqxa+TLqrO337XabypW0U+JWMba2YPQ2N+0K1nXt25Z9Tt1gZGNTWffKpkbj8Vhk33r7IpuCkYcYKwi9hW2qe87S9p58WsOn54u0JfysapZhGhdhpQ/I3Y1oikThspnHfhwfzfJ4su7w/hyWcbyQpVPGAOtJe0Ew3pZNteOArNphmiY2Hm2HcRzYl6qUzLLBL1dSQiKwj8rZIKadoMD2O617sO9i+ICCr24Yx4HG5rv+pSwqZf0rm9L7eovoW2/ny/PJRn/4/G/359hND/Ue0ppYkQ7YxA9O/kDGdQ2Oq+/RD14u9CoHEPjqvrEsVz7oa4Fx19FcHlK4ZRUDTNbh3bl5PJk9zez9y3+YMcS6t7pIy6am5iZfyMsT7On2hE/6O+MNIycluUhSQmJKTm6FlJgbzEM0sz3T0sKkBFRo6LuRfUmxpnHp87+n5sCCo9D6l/7z7yl6/tD6h25RBqA1MVb/4a79rDd3VD8Rm8rYIcH6Iva7i+mIls+Giy0X3brcumbGt4CkEOTG/41/n7AWc+XcrtwGqKYp1vxLXWNo2tTSNYM0uge7dmD0Yh4Q2pSSzPPnSopLbFJyisWIp37MuorqpxacXtJ09/VBblP0/KHdwquc+n+99e7ShGFT6NSRvKLIu/zmk0v4A++hRJmHp/HSlBTtQ7N9+VK9r1f2xiWb4n0ffealWM+01CneK88nlyCNHRQmObW4DjfLSUO4J8m8PrmSEhJTcorFXNdxnGihEFUfv/r0zHIom2p/1aa0599gn8CVqbfesHZdzWrWIrIZNhXmFEjBMsl28mTmIB/W/TFM2khEUvj4tNNj5bOPsp/y2lSQNQ7tkLPZT4n5S5Y2o5V2sjHBliySzDOnKoGS8OqXmJRTLCa6HRMLkXlVNhXaGJPJ6DPJWD+lPf9GrXNCffQmtqt6x9jUg4/9RFfQ846DWdZNG7ILCl9RF+hL2p6Rem0qxToxEqNXUYuLBvSeGpRkXq5cSXGJSTnFYrDj8ymvNE7ZlPpiLeXxwfTMptTzh49Q7SNVb728iq/4npK6u6i+zdHmRc5R8KUS6KGYaV0fTbc0zLLYPU66CQqHAUSIjuUCn9lem0qxhjlk42f1U6wMQnH8OwfENU8+YAWDV5I5h1xJCYkpOcVixFgaCn7AT9nUD877wSQO+ySg52/blGq9hZq2XALMfsrjKWW4+VqZG1zkkNPLD7CgHr6xoBEza2uh88EXgT0BrUrd8TmKK3xRiVnEhzn2S7GmqXLJ0bKp29rIrYtiChK2xs0w98fVqSTzvLmSEhJTcorF8AnAduCfnEpt20LKpvpfXJ/iLYyeP8xR6GvAWustVa75ZbOu/6UxYCIsTymzGzAzw+hOzJrAaiP7wfTuNC+XkY3jmNHD2hDOqkgKTmb+3SZYnoJP8Pn6fF7b+bo+r3N/XcXEOkwWJ1jziXGNo2VTMGeOiwWw3W3iu01hGWK6wPoUSpa5cyWlJCbklIsZYXkK+3utpuXBy/UJehWsBjmLU87lY3NAC4PpMPCxwOfP2prq/rXWW1hFyyXAtCn2+W5+vJrbsK3M0I+ypVP4sRWU+C9NEcgfz7gIk1N5bZtSKclQnqTdEjeKSRa/jACXycpyVWqfBiyXAHPs53pK8W4q7E+AmwKcjaWO6DSFk0VExDM6G1J37KHNk7Rb4kYxAfVsjP7NfekblZXKZrkEmDbFMhueUrybivgT0DYn9j0U/6UpAvljGRdzaBrgkBudI+kAiW8SE611/3tDv6g+diVaLgGOTRmeUmL3Tsyf4CLHo7iFPVy2NEUgbyxjyAclwCoRnSPpAIlvEhOtrP1RGCWuiVEN2C4Btk1JTynDzZd9ZYX8CcjBTuaICE9TBDKHM/bJ3jHAMhCdlnSIxDeJCVSSRdfzKCLKKU2yXQIsm5KeUoabb9yfIGNXYGkZK/2rNXDAAPbVRfwg/rZLwL9/+oyd8pTS3Xzj/gT1fL8PevyyqBMucXxe0U9YYtslwOynlKeU7uZb6E9wwlrXIlUNvE8Dpk1pcmFHTIYvtpahBqsGqgaYBoI2pdx8q6KqBqoGCjQQtKkmw823QE4lrRr4FQ2EbapO6P1KG6j1PFYD5n4/L2/y9fCm1siqgaoBQwPhfgrJDGcPjKzXqoGqgYAG0jZ16I66QClqdNXA92ggaVOwi9bjqagpoK75asp4WbDit71MtYczTtoUdFMeT0VVjjqVoXTxwlDFb3uhcg9mnbIp5uzheiqqQtAe2uSxwuDtqLKFQ2GicEqYWywlyI8SKBBjk0oLM8EUvKY4aenoC6pFBYN1v2xQNa9JSMz7cbcE11NRleVoTDfTr1jJgdAhcGqKY1gSCaKAylUcCouhCpWLcfyto8XSjh2O0v1C4hGYbg6Gm6W4UD8l3RHoqEbDU1Hx2IrpRmc+K1Y8FIM7ayLNE9kkQc2QsIlJIkEUoGylgZgYqlBSTMVvK9V7gP4QTDcbw82WFbApw9kD8hieihoP9J1XQ0CVaBxI2JgUm/DbmjScmjy3SxUiFDKLI6hUgUkQBUJ8EvFxMVShlJiK35bQc2byIZhusKffxHCzhQdsSh59hrtoj8d08zsbOGc9QHHlUcys4PGTHOyqRe/jkkgQBaLMwolxMVShhJgs/DbnIA517lg9z0U+IX6+X/hpZaeYGG52tqBN6c4eb8N085w5h8ed8YLbcGp2bfLvE5JIEAXyWeuUCTENsseryluM3+YcGKVs6ifx25QqVegg/CkLw03xF6GgTTH4i3djuqXgzo4DzExJojc7BWzFZd2nxFCFEmJy8NsiNoVHm2aV+ZuJDrIpGH3Q4UYedYXn/ZSzh/JUtBkUY7rtw29rwqfZspJp+G1HAaslJO4VQ+xVxTbjt0Vsyujr7Wf4O/eHYbpZGG62BsP9VI6zRymmmzgxiB1zuQm/rVE4Z3Y92L3Cb2uOAlZLSNwrhthTxbbjt0Vs6ifx23wthH9PfRmmm0JR24Tf1tiQGobeNPy2w4DV4hJ3iyH2WLEd+G0Rm/pJ/DajbcgbblNibpi/Z74C041Q1LwYHnJ3gA6uZULDEZyai6wFQz/CbwP0A4l9QMAoFoBbShIJooBX4l4xANwinjZet+C3SZBME8wO2Ko5ip/EbwvZlNYkCEiJ0+rNzpcZ48j3aSU4QUwS18jYzyT03hVjuhGKWsymdHAtExrOhlOzCsXmKvnvMGA1QgO2JB0khthjxXbgt0X7Kf1wLG9VfiOS9VPfhummUNS8NpWCO6OhUqAFoIlDPyVBUOilZAG4pSThWCwucbcYYo/ysOMLVDAWHbGpn8Rv8+mK2ZRqEjb+FBvbCERCnvdDMN3Ec9+M39bYcGqm3hR+m/jQiQC45QKrJSTuxW8j9lgxuQ9EAbiYNYzdRWyqzvtJxfHvqf2YbjYMm/VYzDMzPZ5ShpuvzQwXPlI4ZTAnx8+8Vyhqm/DbGhtOzayMwm+DeT+GOhUGcEsVmARRwBQl7/aKoQqhmB34bRGb+kn8Nt8DOwrTLe4WYH5PeTylTDdfixm+AFM4ZQjNhihqG/HbmsTJ+Qq/DeDlogBuqQKTIAr4HhFDsdslhipEYrbjt1k29eP4bb7HdRimW8ItwLQp11PK3DLtMEMbi+OUsQqmKdJEDpyaT3N5cfHikCAK5HF1qOJiGmSPVyd/QYRlU3rOxC4NnbSGMzRgYbjZOUybcj2leDf1Nky3xF5S345Uu0KZ93FJJIgCmWxtsriYBtnj1c5ecu/soaXMFb+NVHFMwMJws5maNsVSDU8p3k29EdMtBnfWHIpGEZNEgihgqy37PiaGKrRfTLQ8Fb8tqp7yRAvDzWbg2JThKSX2Er0R0y0Gd9YcAKemqh+TRIIooPIVhmJiqEL7xURLRR9rUaqamKsBG8PNzmfvof1rTLcw3FlzCJyaqn9YEgmigMpVHAqLoQodISZSrnoeRUQ5W5JsDDebh9VPVUw3W0Eff//igeXH66e8AjaGm83BtKmK6Wbr5/PvK37b0c/QxnCz+Zs2pTyldDffiulma63eVw2ENWDalEan3Hy1yBqsGqgaSGkgaFPKzTfFoqZXDVQNaBqw5/1UUsV0U7qooaqBfA0E+6nmaw5CHx756qiUVQO7NRC2qTLW7Rb3hDIRBdSGFX3Ny6FAAZX0DzWw2aYeN73U/Q4HVZ3PK8JHbKZ7Rbkqz6M0sJzrRPitNjVwfyhUynAqk+rGkRWs78Xp6Y9lqVsJ8El951Ud5XaG+m20KeFciBUYV9PLCuP/6Nrx0V4PrtD8GNh16aK77v6olFXsgRrgLqgH8tvFKjzvF2PbGnW4rXi2SizPi9OWQfzYkPTKjjSB1e6Zmfq8XvyHs7+4QJX9GzVwWzeAeL2qfNv6KTx/VpTqvuNkklfUq703YFkwxh6Zq8o4wXDb0rh16ylEmsKTiUUFMwYTAoyS0UGGmIDXJKcoQZALJuA1yuXViZcTjZQ22dTM/IHpt62b0oCxyOWRWJYE2GHI7PAJ7XcZoWfqpm5k0xPD0ky9OStpOi9rGSmYpiBSMxDOWI7cZnK279KSDpG4X8wRMGp25bXWI5IWDf7FJn73/SabkscFybJeVmPqOrMGChhLuTxmZjXIWO7GmDAxkn03UZw1nsGhsNlswaQLt07F/VhMuhyJSrY35GiivT7Z2dzsh+zxKmLt/0Ng1GymqvVgyhVLhRF/d91iU3jqoCz1HQ/TK6kFdCQIjKVcHksYIC0cL7VMVj+Faf5rAmcNMvkoTF6bMOlSyG0gQp5FZgrz3fmK6GLSZUj0MVdx+8UcAqOmCiRCWuvBJDomFSP+7rrFph7mK2Fd+dR1cR0QGIsNJMfcnsY5deE59F3ZHIRvvcocOfgorOr5RfoyKtbO4RROZSwhkds8SbslbhSjFZwfqafdHxXE1kP8rh6sL0p8a2DDvJ/9RthoUwiM1bIpm2vf+luprQzndKC1+Nhij+7xTDUpzUNhl8N/78mos7aR25zK+Ln6YjMl7ZW4UYxW4oMgnzSOPIitR8V36uWlIv8ktKGfulgz5xtYsKrCK1AAYz1h0nvNfWU7zfDKph9KOsokzlrjo8h6OL6Mk9Ys7TPBnMpkSWFEuZJ2StwqRqvHi2yKWo8StemzXmU/LlRuEK29wFvOghefgLHI5TGnVk4zhEnDnk3zZf8Q1rtZhmlchDWyg9ZhllCyIQqX5z5MOnu536qMhknXwlGcE3tbtMM4DmyWVUqWRaIiBiqBkvCKNbEk7sWkQ/Z4RTF0PQxGjTjKALUelXApH7KozEeGyg2it/vYchY7KmA3inJW+Oq+Qf/YXHmf265903U0l4cULm+xH2MzJh0ht0nOVmUUJl0KK44wS0OVQEl4xapYElNyUBOlYlAcXPn31Mth1JjA21mWScsNQjRDTWvlLLTMpUGrUZRmB/q7WKsSW0EGPt3Srcuta2Z2zgD7SQpxY/zvxKRD9A7kaVZGw6TjixUdGD3hjynJPHOqEigJr36JuzHpkD1eUYx25TZF1QDlgsZpWKE/BC3PxuCqDbM3sjgkW7FBtKs56wdAavCWx5+DfKYBYG0MImsBXWIjl+UzRT4pODfYJShWtJ26AId9mHSI3AYfk3f2MysDLU3OoCYh6bCIOgyZgX6HkvAakIhvdsKP2Qh9p8SglunKbEpjvwVGzX0Q2lMnQWzQUTxfpeU+LqgbRBbXHp8DURezoJwbAuarfQMDbI5h4Dik8DHfh0mHyG3I2a4MYtLRzCquBMInLEnmmaXZByuBkvDql5iUs1UMioMrsymtGo1aRINEvfxano3B0WmaGxntzFZsEBdjXxKTnsmiG6eLvqWJzQrYMUZlvMl2MzRyZN2k4Nxg9625k0njuhOTzh4i2ZVBTLokVpwsoujPBAwZwxe90YsaJeEV62BKTMpJ6QrZ4xXFaFejnwLcmIdEhWUkevl5ls0wajz3fJKZv0yD4GXmf08HxDSTBQx2W5EXUaxUjGKvhbzJZqPQqHmQ8LNQhE0A9yk4N6TwZIW3AIvdjEmHyG3I2qqMwqRLYcUlK4GS8BqQmJKT0hWyxyuK0a78e2o/jJrGkQe9DzgEsGtnfvF9pkFQKRbnc6r5z3/+l5IjAZgAEOtR9MGixfjyeZOtZmjlU/hZctRipbNbOe/7APPo70MDr022p7mFVzy+QYMzw81OTDpEbsNSWZVRmHQprLhkJVASXgMSU3JSukL2eEUx2vUoGDWNJQ/6HvDDGUPZud5x3672zHhCamevTmWP/YDxVYz9NBQrGRMS6iZbzdDKyHS6sAGQJsIiaVJwbog/Z+dj9zsx6WwwALsyCpMugRWXrARKwitWxpaYkJPSFbLHK4qh62EwasRRBrwPeNi4S85mvu9+XmkInsdocL8DE12d8uQY5LQMjBQWsQEbYzTZ+i5+T7LdKLScECT8LCXCJGB3CZy1LAqXK4+Js3aQ2+KVCcgQ0XmSdkvcKCZa9EMSvQ/4HKei9KVD0Ifbr8VtSnlywMHdYnRFkFgUo6lZ28XvS07vYeL4WSRC44xBZ2MpJtA1TUGkZiCe0dmQmq6MyV67y5O0W+JGMVpBXxT0PuD5FCtUl9LDWZ5uvxa3KfLkeMCigljaQkgsFaMUD4s06APiS1aEoZDAz0IRXqoozhrPkabwMobRrVw29iUfC7CRI+kAiW8S49NXNM77gOFLJprpLYlQiMJlMs/RE3GbYp84pieHA4k1X/HHW6Szi79EFQI/yxFhsEif+JKmMBiqm1jGY5HbciQdIPFNYpQG80KBB7ye4FiKqdSdcPG8CaLuIh5PDhgKw7dv8IXu7uLPUzOnkvhZCRERnDUpK00RKFQ449HIbWlJh0h8k5iAOkPRgQd8//uJP8/EeKgSMr5zp9ITa76uJ0cCEsuziz9RKJWM+FkJESpDDX2mBgIP+KFtk/ubirXlMAOjZ6k6PvZzPTkSkFieXfzZ6kH8rISIbH6V8KQaCDzgizsr/d4KtPfy8fbkLk8l+qn31qlK+2kN/Pl5yF3p/AQ8Lt88Ybyf+vtnXHE9/v4ZvKkEvlf+m0RLMVP56PPqmdQ4oU0Zh6NVXI/3Nqs/lDZ5Pk3eXJyh2KgeXpv6rzeXu0ycs/BZlr1Sf44GRs8U2rtLfy8d/j08H4Gn66cqrse729FZ5PUnsKlbaRl8R6Ofz6YqrsdZGvmby9GfYXPSs3D0d/fs/D2HTVVcjze33zOK6zxbEt5ezqHwo85ng+ewKU11aVwPjbgGv0gDp7CpUkwOn/344v72OaVwPcIbo6jcwa1TRAGBMFE4Rc+fHw7yowQK5DN1KcNMMAWvbt4zxEBz/vtigP+UeUZEokjr+t8ORXS/n0ONEZqflPKwwsSSqwcrJ5U9DkTBcqcp4kSHwNSoaoSLQ4IooHIVh8JiGmSP1wLeR2DlaK0lJvn/PFu8Y/QvSYON6eogsQwJ6/o/DlUhC5lf+UkpDyuHdUbE67Fy6CxNqzQxGJkcozwWFydHolUB+zarQhG7s/mJ+0OwclRr8QuRsaewqdKNRT5X+002pflJkYdVVF2hRDhD5MVYOZtwcZo0TM3BuDgZEkM6FPFxXBxin66YIeYQrByttRjM7RsYdgHpX/8KDeIwm4J6o5+U62EVUYrjGvsGrBz/Y/ItJWtO0HEP2UgN3aS4JBJEAZdDVkxcTIPs8ZrFU5zvl0kaJcPWEiX6dZtCPymPh1VEb84RDqWr1hKM3hShA9pAStYZwR4ig40NU2NKLLlLSCJBFChhrmgTYhpkj1eVMxY6CNcDW0tMFDvw9rf7KXgtfjJWTgpGhmAC4q0gIzUlCYECdkpMiSH2JC+j7NBPsfPh9v+otURZfZFNbdrvR35SrodVRG9OP/U6rJx9uDh4BF+gLhouzlGANQmJe8UQ+/Dxh25dD8PKodbiytBivsimiqYONRVsCDo2VcwDX7Mh/Bf40OP1EScxsFNp7hwY9Qm+pTC4wAWaJBsbpsYsqMLFaY4CrGniEveKIfZxMWY15ffUW7ByQHK1KVv9Off7bSoFMwM2xM9LV+g0V25TAi+Wlh6SbCJniMNqcdeQW1cKGCcliQRRwKvIvWIaZI9XrxQnkp/t/B6snC+1KQfWREM12RjEp+SDl2GH4Ob+kI88FDgIM0NHT+/DxYF3ppTo6ATiYSoRd04mgXFSBSZBFPBKRI8CDbiGashKmhJDFSIxsn7xi4FBAOc7Is4Cz6U/hBgb7RnHyE7fT12oIlo1jpxL19jmB/f3U7Lt6DAt98cwaWeySQpCp5H9FJyqriwBm2CYjQ1TY9URcXGgHQj8Ux1RhkSzTKkCkyAKWKL47W4xDbLHq0+KG8ds6k1YOWfpp/790xqTrpHe+5G0w6Zc6Bs3Ri+AH0pnv02l8F/gVcrHfjtxcWioZNRJ3dD7+jDAmrjE3WKIffnYj7rG12LlnMWmgmu+/K2smoAM7bApWP74KKwcYb2bcXEaOIIu8lO4OAIAFF7mDSHKKNGMQy5gTULiXlwcYh+vmF1n/j1FNYMXFqKnACFTAtTv0XQ482PnpnsvNg6lYuDccxRsWxehOWGRAcHyn1v7oFmqbCwEGT8KK2cnLk4TwY8BbShcHJj3Y0OCaQH1sDEDNDolGu5xDjsI7kOCKMByOb+9YqhCcTG23KOwcuQI2GZv3rf//pkRf3IXMgj2klVoTlS0Hf0U8PgorJyduDhNED9GKFPh4jQHAdakJO4UQ+wTFaO2wgKHYeV4sXEMUVzcmfdR8N3HbCMeR3NSZd9gU8qTA6FvYPDwBVg5cRiZxoGpUUosDcUlkSAKlPKX9HExVKG9YraVTrWZWP4zj/3EeiehOalqeGwqAQunPDkI+oZwTihGCWi0Xf2+ZGcPrZZVBN+FlZPYS+rbkeqUNS8iLokEUSCPq0MVF9Mge7w6+V8aQW0mKuUsNmVsLJInzk/kqMhbqKqHx6agHoGpQ56NPDkU9A3inKgYJUDb1e9LVoSh0PuwcmIwMs0BMDWqhjFJJIgCKl9hKCaGKrRfTGGpBDm2mXjmUpuaH2L2JM5Vpob86JzM5vfUbeL7ehVAimihKpfHphKwcK4nh4Nz8qFYOUpLSj8UKj82m7K6gZgkEkQBN39mTExMg+zxmsnzIDKnzfj5FtiUnH4U+8783KzYgB+dRQW3pk3JqU788IGXk+U75LEp33HPSo7HkwOGxt+BlROGkWkOgalRWgxLIkEUULmKQ2ExVKEjxBSXC0Bl420GOebb1CK3y4v1fMwfv1q2ECS2bYot79P7SqI5aZldm0rBwlWsHE19NbhRAwFsHJtbvk3hwQUlNmVLC93bNtXAdkfsphDNScvr2lQKFs715AjgnKCUvF39SG1eK1aOqY+vuUu0Gaxntk3NeEb5W2yqe858GzYUE1solhiujk2Vw8Jp3GqwauBQDaRsqh3GcYD5t+6y3i8Xvonjsowj8+KB38JSYTn+/hy76bH09/Xe9Y/1cQM4m2e/PJ7MMKRXHVHznOafe5BYFOjNtqkNsHCm/HpXNXCcBhI2pblx3fF76gqbWIRX3AyfKc3IIGHWoVvAqadlu5tbflgfnxJl8xniuwi86hS1U3xn7NdEoWUsm9oCC+cUoUZUDRykgYRN6W5caFP8FA7uFSfMjG1BvIPNsAmJCzMwPtvOZwnZOJG86hS1U3bXphwSPcI6M3MLLJzOroarBo7UQPx8P9oeDxai+ikmn3VIs3A3ZXt80WegA6/ubgI/uBtfr+XfXtLnTKN2KmDYFHlMmQEtk9VPNRtg4TRuNVg1cKgG4v2U7sbl2lS3Dj37wQaGh+zEmufUTM2zgz/24zYlXds0aqcGhk05qU6EbVNNOSycw7NGVA0cpIG4TVE/BSbDbIpZCt82L/spNmfBf+SPMtzbsRmumk2hVx3ywiz6NWxTqq/S6B2baoph4TRuNVg1cKgG4jalHNREV2TaVMOn9WCcJzc+sILN6wVccNZRODixfopc2xS1U4OgTeX6+RbDwskiaKfKq53rTvFqRNUA04DWWmKP9jEvAAAKq0lEQVQKSdgUOahBHwXTD8ymHupEhJm7aIEfG31sQYjNUvA/uLI+jVzbFLVToKDvvN9L1e2nYBXL62XvSLIi1P5ztXPdIqm3X6GBI3A9VGuJqiRhU7ACNU0XsT/8cmFeoNf1eZ37K6xWAd/lAslLc2ORchw4srbN/5gr2PPaolcdUXvKE+qnAvANHpsqhYUThdD2n9POdU/xatTHayCK69Fen+wExdRPay1R0pRNRTMflhiyKdZNecZkHpsqhYWjouOp8u7OdSKpgY/XQArXg+bYEjXF1hIlO7VN8eMowPFzts5099gU1EP0ptHauol4qrxn57pLXGM+VQP8jJdI4TUolAgVDMiaGxufxX9ntim+A8M3JvPYVCksHGoF/EU3YBBg7nr9EA0kMAgMKJRIlai1RGjYPJ3VB0SpX5Vojf10P1/fmMxjU44LVmZRaf+5u3M9k0Ml+wANJGxq8iD0+GpFrcWXSHFnsSndd1738/WOyQ60KVJEDXyvBpK4HmwbUDcOue5+CU2dxabMmXC+hCwcqFxvQo+vB1TS6uoS1a7Jv6UB/j2lbQiXe76lElpYiOk6ddKDOBT/Ln7p7ydblee0Kc3P1zcm23Fmpl3/ev8TGuA2pW8Il3u+ReVhl+qta2bmcnjE76RnZmp+vp5apsZ+DqKE2uG0NeQpRY36HA0wm8LNcQxaBPd8ixoMT7ZnNfnT2k6U9pz9FGy+ID9fT/FTNuXJUqN+WgPMpvQN4eZ6lIWosldTJ7WppsjPlykh83vKRfFwYwydJpIN2npzVg0Y/RRsBqc936zAzIn2ph8P+SXfU/q8H6tniZ8vo8+1qWNwPbjE+vc5GuDfUwFcj6YHfz8G7iERuHbX6qz9VLRiO8Z+8Bl6AK5HtHQ18XwaiOJ6XGALXAv9FNshfsTva2wquLXd0dIBuB4OzxpxZg0kcD349u9xGA9anjrnPorU89nQT6mduHhuoMJowBhNrO4p40nWKGvwlzSgmlGk1l/TT8XnKJR3FMF0EEYDxWhq0jxlfMkaZQ3+kAZUM4pV+kdsinbiKpgOxGhQMUpPmqeML1kR1tBPaYCaUbTWZ7Epe94vWujysZ+7E9fBaDgS1yNa+pr4sRpwm5GvKmexKXO/n6+kWlyxTXl24iYwGvI8ZbQy1eD3a8DTjHyV/iKbinV17k7cBEZDnqeMT6M17ns14DYjX10/1KZcn974HIW7EzeB0ZDnKePTaI37Xg24zchX1w+1KTpYkOoUtykiq4GqgVdrAGzq1SIy+BcaxNNzfEAhCyzU9u39Azt0DX7bOYj89f/bNADnDZ2gSi5WTrRQd89hfttsKnDYWUC6tCKRijsSb0dtEwsIrdGfpoGz2FTRvN9hNtXt2OIFsxjiV2aXn9ZAanmLNfCRNvXw9lOxeT+/Xlqj4/HTqNhu5D1S3wsU1sey9KKLehy1UUzJqqEP1kC3Zh4Z89I6Fg7crgypx/oVsuC5qauxePlvOz7a66FrEwhcSydhvHsENfHnq7E/poH+M23KHbNtsSn5Omn7cH+1DOLHJu+vzMEa5uBndjwugC1Q9/Sk0I+1nlpdnwbGNeekaF/OI+MKDeKCeN1aGQpZsJw3YUrzOCodRLqc9s4QTGDhamRfUOO03Fs55/fATyutPDX4sxqY1vA7+n1KKTSIYXUPiCpkweo2YG/XKZvyA4kITVyYh003wUFwYEPD0ky9XCab3OK8T3dV0tk0MKyxVvSu0uY7FPISjZ5Cx21K94ciHxjqXnw21U6jQCfJ0EFk7JiRu5J8mQZ8w6j3VzFuEE55Ok/nGl/i0vyhlA/MHXdjeGyKUTW53Y+W3ylrjfg5Dfim0N6vhEKbmj0TK1EWMGabllZ8KykfGI5Ixyqr2QT22nAoyDLlfibNZ1g3f/9TqxL9Gnh4tvn4KV8ZGzUIVzCAeDiRKRaEHKR8YCybGq7we7I/sKXn0HfZs3m3M6xHOBqpEX+kgad+rtkflSH7IDEq37o6u+wSNkX+UJoPzCMy9itSi9bPURFr4Fc1AG/87LfxC3WUMAhH8p11NuYvwUL5QykfmCuuHGs2gWM/cZaOKSJ4B/6N9Vc1IDXg+zL5A+XEJxjcAl3czUkJm1L+UMoHRk6BL9NjHdC80DpgcrDPPpyKcrslrTGn1MD8EAdoZhVOAX5kkfeeGbSsjMcSJQzCETa6i76lLIDnTS1LkQS0KYrICdydkWhOrkrzJxqQi5Ic4z2vADe3tcUyTgzk4O9/pQYxu6+CwiUuXmc2X279IvsoLEp1u3hsU6XW0Kk0sMgnfM1dKIHSl30eXU8x7Zd72Ll6OO7EX6lZMl4HbX4t24qrKlFDf6AB4VQAmzcLbKqsmPdTTFGU29TdmfjbYlPN9YhB27JpvFj2nCr1QRqYcRz3MptqPYunBxW+iE2xQUzOJEUxC1ZA6axRVFaH+Fo2NHDy14jXa6AdxnGAueLust4vF7aGcr0s4wiHZ7HfwlLZ9/Vz7KbH0t/Xe9c/1setuazPfnk8WZ82j/04ss2xkppl9P06z2ZUH92r44oN4uYUvJgFr5NEuN9Tv/6Ivm5PAWrepAY0JN87fk+xlf2ZD3dmWF5puLlImN+W+ecxVCowPfZ02XyGeP3Cl7Oi9ou9nONzqnzs1zztrQvbbMqvlhr7ZRrQkXzRpoRjKdt+JsyMmcIdbIaNOi6sP+Kz7XwKj40Tb9xJDm4VtV9JT3ePj5/wxbHlBjHYGx10Fi/A89WwXP3BF2uost+hAR3JV/VTjCHrkGbxHcGgExHmF0CzwakHjOPG9xbwb6/7eh2h09KovSUCX1Vv/NsjdYPIE36zFwHKWeQJqlSfrwEdyde1qW6F9X34wcIKwfw+J0B4e3YS5o3bVAtbA2CCQ6P2KmY8ydAPxn64j8FbTl+kPfjbalMOXq8TYUpPJJvE9e4UGqB+CsZ9zKbYgI4jJ8p+Cvd9Kpjf4d6OzXDVbGqBMWHbrzPyClXsYX+UhAhfHX/HGc58QaO15W+rTYEKWj4C7u/CrinCX5ZEsj9Tjf1bDehIvo5NNXxaD8Z5GnQ2jODg3JGVDfbgx/qpjreOS69R++o02+MnH9Fb4i7uvoikXGvcuqGr4yLYEFm8WqTrh4rwFiGR7M1TI/9YA/OTH8/D+qcJph/Y9cFXFfn4aOZPfoIniwNDCLFZCv4HV9andXyOAqYBFTWkOL/LORZ8oVxwymDxIs9kZtnQ1aFCBMDvouZrJOIvptvXRLJNXu//XgPLME0XPt/QXC5gPfN1fV7n/gqrVVC45QLJS3NjkXIcyM9NEIcnABrw89reJliemhgLSe2tVGu96b1Eb4q8l39QtWYvW9rV0akUjcTrrfi+b3rW3yxmdPb3/F1tb2v5GXmj0VEVdnXqVArE6634vn/3+L9GcoEHyevrPG3oNGHCU/uVdXV0KgXh9VZ8X02ZNbhJA+ZrfhOLIzNtOMCJrcypX1lXp06lkBwqvq9SZQ1t00Dr7ELdxuewXGN5T3UxdoSXdHXaqRSyAhXf97An+bOMzPZ4BjUstPSWXZq7cVpYSVenTqWQwiq+b7bWK6FfA2JLoD/tc2Jnc2ajoKtTp1LI2lZ838957Cct6V1M1p+0dNnFuulfVLB0UN7VZYuqhFUDUQ20ZzGp/wdk1hJvd4vE3QAAAABJRU5ErkJggg==\n", "text/latex": [ "$$\\begin{cases} \\frac{\\sqrt{2}}{2 \\sqrt{a} k} \\left(\\left(a k \\left(- \\sin{\\left (a k \\right )} + i \\cos{\\left (a k \\right )}\\right) + \\cos{\\left (a k \\right )}\\right) e^{- i a k} - 1\\right) & \\text{for}\\: a = - \\frac{\\pi}{k} \\\\\\frac{\\sqrt{2}}{2 \\sqrt{a} k} \\left(\\left(a k \\left(\\sin{\\left (a k \\right )} - i \\cos{\\left (a k \\right )}\\right) - \\cos{\\left (a k \\right )}\\right) e^{- i a k} + 1\\right) & \\text{for}\\: a = \\frac{\\pi}{k} \\\\\\frac{\\sqrt{2}}{\\sqrt{a}} \\left(\\begin{cases} - \\frac{\\pi a}{a^{2} k^{2} - \\pi^{2}} + \\frac{e^{- i a k}}{2 k} \\left(a k \\left(- \\sin{\\left (a k \\right )} + i \\cos{\\left (a k \\right )}\\right) + \\cos{\\left (a k \\right )}\\right) & \\text{for}\\: a = - \\frac{\\pi}{k} \\\\- \\frac{\\pi a}{a^{2} k^{2} - \\pi^{2}} + \\frac{e^{- i a k}}{2 k} \\left(a k \\left(\\sin{\\left (a k \\right )} - i \\cos{\\left (a k \\right )}\\right) - \\cos{\\left (a k \\right )}\\right) & \\text{for}\\: a = \\frac{\\pi}{k} \\\\- \\frac{\\pi a}{a^{2} k^{2} - \\pi^{2}} & \\text{for}\\: k = \\frac{\\pi}{a} \\\\- \\frac{\\pi a \\left(1 + e^{- i a k}\\right)}{a^{2} k^{2} - \\pi^{2}} & \\text{otherwise} \\end{cases}\\right) & \\text{otherwise} \\end{cases}$$" ], "text/plain": [ "⎧ ⎛ -ⅈ⋅a⋅k ⎞ \n", "⎪ √2⋅⎝(a⋅k⋅(-sin(a⋅k) + ⅈ⋅cos(a⋅k)) + cos(a⋅k))⋅ℯ - 1⎠ \n", "⎪ ────────────────────────────────────────────────────────── \n", "⎪ 2⋅√a⋅k \n", "⎪ \n", "⎪ ⎛ -ⅈ⋅a⋅k ⎞ \n", "⎪ √2⋅⎝(a⋅k⋅(sin(a⋅k) - ⅈ⋅cos(a⋅k)) - cos(a⋅k))⋅ℯ + 1⎠ \n", "⎪ ───────────────────────────────────────────────────────── \n", "⎪ 2⋅√a⋅k \n", "⎪ \n", "⎪ ⎛⎧ -ⅈ⋅a⋅k \n", "⎪ ⎜⎪ π⋅a (a⋅k⋅(-sin(a⋅k) + ⅈ⋅cos(a⋅k)) + cos(a⋅k))⋅ℯ \n", "⎪ ⎜⎪- ────────── + ───────────────────────────────────────────────── for a \n", "⎪ ⎜⎪ 2 2 2 2⋅k \n", "⎪ ⎜⎪ a ⋅k - π \n", "⎪ ⎜⎪ \n", "⎪ ⎜⎪ -ⅈ⋅a⋅k \n", "⎨ ⎜⎪ π⋅a (a⋅k⋅(sin(a⋅k) - ⅈ⋅cos(a⋅k)) - cos(a⋅k))⋅ℯ \n", "⎪ ⎜⎪- ────────── + ──────────────────────────────────────────────── for a\n", "⎪ ⎜⎪ 2 2 2 2⋅k \n", "⎪ ⎜⎪ a ⋅k - π \n", "⎪√2⋅⎜⎨ \n", "⎪ ⎜⎪ -π⋅a \n", "⎪ ⎜⎪ ────────── for k\n", "⎪ ⎜⎪ 2 2 2 \n", "⎪ ⎜⎪ a ⋅k - π \n", "⎪ ⎜⎪ \n", "⎪ ⎜⎪ ⎛ -ⅈ⋅a⋅k⎞ \n", "⎪ ⎜⎪ -π⋅a⋅⎝1 + ℯ ⎠ \n", "⎪ ⎜⎪ ─────────────────── other\n", "⎪ ⎜⎪ 2 2 2 \n", "⎪ ⎜⎪ a ⋅k - π \n", "⎪ ⎝⎩ \n", "⎪─────────────────────────────────────────────────────────────────────────────\n", "⎩ √a \n", "\n", " \n", " -π \n", " for a = ───\n", " k \n", " \n", " \n", " π \n", " for a = ─ \n", " k \n", " \n", " ⎞ \n", " -π ⎟ \n", "= ───⎟ \n", " k ⎟ \n", " ⎟ \n", " ⎟ \n", " ⎟ \n", " π ⎟ \n", " = ─ ⎟ \n", " k ⎟ \n", " ⎟ \n", " ⎟ \n", " π ⎟ \n", " = ─ ⎟ \n", " a ⎟ \n", " ⎟ \n", " ⎟ \n", " ⎟ \n", " ⎟ \n", "wise ⎟ \n", " ⎟ \n", " ⎟ \n", " ⎠ \n", "────── otherwise \n", " " ] }, "execution_count": 18, "metadata": {}, "output_type": "execute_result" } ], "source": [ "phik = integrate(psi(1,x)*exp(-I*k*x),(x,0,a))\n", "simplify(phik)" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhwAAAB+CAMAAABYgzg2AAAASFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACCgUnDAAAAF3RSTlMAEFRmq5l2u0Tdic0y31yFIu+vucMwIPkEVOsAAAAJcEhZcwAADsQAAA7EAZUrDhsAABIZSURBVHgB7T3pgoQsjt737OrsN+P7v+kGFE0gKFJoqV3+6LY4QggxQEhCFC1PnCyv/i9BgPg3/6t5DgXSMfMB3Fa4licQDOKZ730dPxNxJ6xzP97ICUdRIHHv1PI7CjWEEO/ok+pFMZbq9cj/ZMAfjAakKY6AenrZIX16D2z4V+OIR9lWTE+PCUU8gehAH/q7Gt8qKOux8RmTrMW1PIFgEE9+z7xk7/177PnNdyNejWpAquI+s3BaYESDjEdfj/JZxUU/htjuBUEuKJBsJCLAFXZLamlAmthrpnJt+0i5NqFsfKSurWwe92lExW1Z2wo/Or0ecw/8EyI4Igqkj5O70Cqtoz7zWVJt0ATANX1MqdaNr1yBj17dauno60DKu5AqO2d+a6OOMkdUDhvs9NgsfVydOqJ/KBoQmFU6JzinF8q0MQzTYFVGlcZ1yStXHePosUnPtDWFBiRptSk5zJh4QDmHOfI06nRJ4bd08+jSlVW0cXVqOtbVZhqQOCvWlbwTxHMKddkwiIklzosih01LVQ9F0rQLbn3eFL3HDCj05brOPBtvIixDklIbVyfQqS5DfYA4tfRpIbnm6GrBDnIdNOZJv4jKaoD0chWCeY0ebdrYRaTyUxftwv1qAZ9xRQSdcPcBckmvJXPUctvZiU+7hoW02rxMOt6cLq39sRr1icYf1G1qeoxrPOoEJUAmDRH8vUEfBXN0s5wTn3aLBEIuZcbEORuoxi2SJ2LqWPpHe1i+cF4h47pBI5SVGhLUAwiCd+KrYA6lkxELpRbtXgah846pwuYDVAqDKh8Au0lVj3HNDIJ6ALmm+0RyAGO0q2IzloOZIvn20ZoDBBRRGl/Tv5Nb8RjXARF0wo4CiZsmC62W9KSCXHMMUlxUQuwT5hDnIWUbJcvuxbmRtGb2OHSWcYZ154L/+tf/HESvN5YcEQUixiLYOu8gclrxUkwd3SCWCpIvajSttDDAKSSswkSrvPFzYE7zWkOgbgB4Rhb96F1wTnQth1ijYU2aIFJ/h+VZVw6j2MGCPqPJhJ6jHIdyOT8FbUzT9VkhOOfg0xuyEwDkXgcRB1u+tjgdV5e2c3PlRYGksFm8BXO49MarTFpGfYGEkATyQiNaOq4upGp1FZguOQSQRt/tukB+TBnQACexUKHhp3ufpuM4cwzmjGEA6XTtMqbi89+HBqZRfTqK37ciNcZ1d+gYm1MdSF8qNeQutCcW6Mc811kD+jG+zpR0HP/32Pj0zAeiAemFXlITuscauXdpWHJUYEWkI1m/bruif/R6j43fibmT1XYrfdl1Xf5i2ZEVUVebW+CW7NkMwj0w4TBzFIwmkAIZ5OnDA2nhirI49M9SY2LJzG2cK8SblqPj6oBkY6o5NMnhAOSVRajT3xu6uMEc6/Ej7mjGOE9uAMFVX/7OfTbP7rJ9XFOs9Vw7WTLm6nYga733vzXMhPvsXmsbDdQZixdX68Mcf8KxumCW6oieD3y1fvTFcghBe9Uyyy4rEFX1TzhWp3+HOaQpTJ4Wjbaf59xid5lDMcmr/6ev05/bxlUIjrjuwIZFU1jUzOGjDcirecHoXMKoB41Cj0qwjGssBEcOB4+Gg8ZgqnrgVMEyB02kuJNjtcfgmI7TPJB95rjEq5tHzilV76nOHOAjLJ5GWLMIdbDhMK9XgHLdaJ7FCSDq+aZjdYDxYBynVc/If4g1QH4bPy7y6jbadU7Qe6rtVqpGTiNScMTiIKlMYzqvMFIi3abKNx2rA4wHdN9wnGbJDd8IJZVW6iKv7vXQszto5Gb01BAE0pxuOnQEiwUIYKCpO/7vn/9onY44vRgp8zXH6kDjYTpOk+7NP/77zz9c8pImDVqXX8FedOTSxb61NDT8u41SYAZzJMD+UnCAWV2eJoUwXcCPKTlgybGNxfccq8OMB+M4jSmi3vckxzmOuyZyys6qP24OrwEzmCMCNx/OWkGRwGSOhlGLqdLy//ccq8OMB+c4TTo4/fgOc5jINfPuIN/cJTAdgC0IdRE3mSMZOuQXZsAwmIOxRtcqfc+xOgxzcI7TWhfFz23muM6rW0kMJUEANVeXca2nJnNEtW4dSQihM0e8XZzUvfjHWeNh6cY2c8Apv/jmrvDqzuQsXyxrj8jXZVzbrYh+J7pdNSGGxhxxva6OSblb/LhsPERvnZjjCq/uSfKvZrzeLuNEcqyH9PQND7TGHMnOYhRXvf5dMscV4yG65sIc13h1C2sk9Ik7uoyDWc7yTENFmMNl9DTmiBqqMTWc0tcmXaAHLiOY48PxWPHf1ua4Mcc1Xt2JCCiwnol5u4x/zBzg6aVtdgMP8CfgBHNcMx6A5SHJAVM3cdwV+qpwXt111K0WF/4u43bmyBYZQz4ZXXJARJT7zixEcpw7Hm7MEV3j1V2kwv1zfmJ58uXmMk5dxK3MYTEEAymrb5+rHTuGL3rdyzXH5+NBSaaIrv/fkxwXenXHNVaAHXEZJy7iVuZYxRIlgskcEXNQi+uIzc+XvO5DjQchGe4bft9mjmu9upUiTOJ3wGWcuojbmMNmCMZIDlh1YD7F9JLvX/O6DzYelGRGB+eEbeaw1fo4nfXq9oVKgdmYQwgO/Xhftjiafl1aXHwdsed73VOS6f1Tv7/EHKxXt8Lp6H8KzMIcUnDox/uyJWZaAaqsyx8enWd73VOS8T2UC1Jb1pnprFe3b4MEmMWIa7LnYA0ZBmNBClHX9nazD/e6JySz0n3X2Mda85MMi1e3H0gKDD56ctyODcEAPj3elw3WDCPsMcfDve4pyaxk/w5zWLy6rVhuZlBgmhEXNgQDKNrxvoTrwRxP97qnJLNS9zvMYfHqtmK5mUGBGUZcyBDMON6XcDlf8m3J8Xive0oyK3WTr7gmWLy6rVhuZhBgphEXMgSbor+vh3sT2PKwa8Ljve4JyezEfZvfCmPEtWMIBpHkTbPVbclhJ+e7cl7mDskZce0YgoE1sWnucT1z5Jtqt+8w3bscqXkjrh3LLi4KxSXMQdgBWSt8hxOYVnMmcglT7BlJFiOuHbprd0/Lrl7CHISoYAl7u4eTqbdD0hUhPyOuhDlIuYA5kskkMk2n4LBt36/+Ga4dPrkct1Q/uckTwWtGXG4tcbFY//1vokmTgOKUzANuwO2lEinQUlgMy0tvRogSazs3tgM5N4eJ3ntug+dC9zHigt2vgRQjObqiQKbxRg23hD6fHnFyI923YFMpTZy68SbXMuB+MNF7cfbj3n2MuJhYrAxzgHne58yB6BnXEbAIKF0KcSYIUUPqu0ULgq9m01UWdeYZr3tGXFwvmFisYZgDWQjETVFoS86sANInTVKItWjeR02qW6RxyF6Yxs23FzZ/QlM7Rlxci5l58haGOVYLARE1Jtpyu+Pwcknr2slm0KVs1A8af27XSpmV+naNu+fuGHFx6BemFiwIc6AAAGD62esO3BwqR9JmvW7tznOV2dGtBhtGdbxV/v55O0ZcXAeYy8wYtzmvNYeyEBjAvT/w/N3Pet3SnTmWS0Y5KphppWnoYhZ6VIqDEZfRH3O7EkRyIAsBamNiIOCToO7OOcIcx9qpX7YeFVdnmiuIPZrUxtURLJDju5UlAIC8a2v1A95DyCG/U1PEacwRf+XA3qHrHxRhx3Ubnqk7Y4D0DUSzPTi+SwCAOAdNaIh5Rd1in2RjnWVie1NmfVGAgYZ4enHHPQis+Yr7tB7rJG3HtoLTxSHt20FMQV2RFoXQ582lRUXuSZiw31y5J6Ux47qHfmWQwQPIXiNB8lG8A3UvZFnCFqSTsq8ToSbkuM9X3MfCVGW6hrgUnqZi8TrpYUFls5bmUctet+Tgrmfj+45TB/3a9rsyB/avVwvSSfsuNqkTv4gxVVfcZ0JCyM2u3NuIKagahASDn2tpTIr13bxtd8176pvPuOb6gpEAQS626nUmjumBLzNUKf//FuJj//pVcojCQkR002pL2EWqK+4TkChJA8cDlXS1kOuTGm4fBTGCSrOtgUKfTX90IhlXx55U+o7eB4hjW58Uw/71JnMkI6xs4AF923LF/dCAkBiS+RYmyRwxLJ5gNYtKsygVL5xVwLvxkBpwIow+r3gBYWgMsQuRdzgUMBJopZ3sRXLAlCKYQ8wV0nt2lhyL/n0JhpDXcRHlJWKOHmaVOB07BYsisP5q9al2zXrum9e4FpqHGwXi71UPBI7lma9ya18SeALvZJN4BwZzRHIzAlMIipQBk4OIvSrmEXiE5EjkxiZLUWkOl04Xplyhx6XRcXVFX5tgR2KRLpZ+fl71Ym6fvsDZrX1NYDHbyca32Dew1hSSo5VWIBLfTjYiboVQcw68iSWpil8gpEwiF6SweVlLc6hkr9OAiV76MccUBHuhUq1UTDLlI6/6Ui4FkVv7lLA0pb9sZ6tb7KFWlgEbdOI2+y4tQesBSX0Gl9z35Ip7eRY8HQiDj/5QxlUDag4ZCX4urSMgf8fa18IWel4iZ8S134spktBSLiMHksqrHh3BLyX5F1B6zfe6zOFxV7d2Jl4uhstk802cm1oYSuNz27sIup/kiAoiRs3jO+FVvx7B7/Rlvdeliae14OLWrhIwBASXy8ZFr3o/YAhwFUoh2vFkjgg2fejRb+gRXvWwzHe7bWC91wU2jVMIKeXWviSgthBcLhuVvOyVfiqXNXt6Q77MIfRF61ONct2mEpRXvTqCV+mW/2KRQu51Mdzau1I9slVHuJbmgifHPuqA4FicANCXOaKMGH83eEWmvOo5J32mC+a9Ljtu7Y5wmabOSaKkOKeNr0D1Zo6oJtoz5NGzeNUvR/A7PTPuddlxa3eFu9NsqOzp+CUUtDvBYY24nBDsyEwCS1R1uLB41S9H8DvwjHtddtzaXeHuNBsqu5a771DQ7gTHX3JEFV51gNZg0UbfqX/n4xK/ljc8lWDnk/zXwg0o8IHksGBP5YmlEJs8x1XwB8BC/SV6UyA4c1ij27IoiqOM5Zn9+6vJiG9J/r18iwKhmUPc5uD7CF828Rzjr6nO7+8JFHDfrVhuUQBj3OWBE3ciCvbwtcVVaIVp3u/5OgWcJYftFgWaPn/8jvEXbHEVUnHq/3u+TgFn5iAKUYQ2TZ+sMUj8BX2gXeIqUBUKau33eiUFXJnDtg6g6ZWaVZBHE+Ue0jdrXIWWaF9Jnd+P6yjgyhxiiLEhhcJQDv1ij5Gr5SjHHAfiKjRK26qa+f3/BgUcmUMKCGRIoVCV6as9xvLFM8xxJK5C2JBRCtnf/4MUcNut2G5RmPzBclhWTLaftVKhM8xxJK4Cqn6wP7/iASmwIzlstyiQ9NUeY7nuahrdXFhhDOIPrCGOxFXozLhjAfv8A+VIgW3msN2iQNKRPYbGHBIHtSDV3eS2EKymXc9WkV/e+RTYZo7ZpUOZ8SJDCny7wmqP0W5MK0fiKvymlfNH3qGFPeaw3aKA01d7jHI6FiHxF5TkOBJXAUzB3v+cGaAsDPX2mCOy3aLApnNbUJ9xPhrZIwwtroMy7/jPC1AWpiu7zGG7RYFNr5jYo7qG1AXv+t3H9ucHKHMh8n6Zfacm2y0KbLqIEfn50zM89jnU+0A4P0BZmL7uSo7IdosCmx7myEyd3Yfp4u2gnB+gLFCX95njWEMyXtKxKkbp3meZYkC5XcKFAcoC9T00c4S4y6B8pTnHlQHKbsoc0aw7/QC99J2r0SsDlH1AfVwVSQ5LyC5c+vfuTQEVGkhEeVHxQHCYIWmkEChAmTeSWkXEHFoO/Fzs/77yYuLz5JQrA5SFotMmc7g2YsTmMhI0SHv5WvE3/FwkB+h9zg5QFopeQZjDiM1lJGjo7uVrxV/xc4rhUYnzRxG9UChJ12kFhRxbotd5BygLRq4gzAEryMmeQwV6WxIseO7lW6o9OrkbhIJQDv3JAcqCkYnGevMHO8fmWs7so+1gXfDVvNfF1ELFqwKUWZr3SKax3jwATFXmQ/010Js65bdB3Mu31fulX0gBGuvtYMOLYbGKzbUEelMJGCC2T+bycdnf+x0oYMZ6c8dqMSxeYnOpQG9LAgaG7JPZfFz2934LCuix3g4gtRoWq0oq0Jv6jf+jQG84+fd+Ywposd6OYLoaFs+1Hhbo7Uhf/2hZEuvtCA2QYfFc7WGB3o509q+WRbHejpFgNSye6z0r0Nuxzv7V0must2MUWA2L53rPCvR2rLN/tvRfjfX2Zwd8t+P/D5vv144xoWQVAAAAAElFTkSuQmCC\n", "text/latex": [ "$$\\begin{cases} \\frac{\\sqrt{a}}{\\sqrt{\\pi}} \\left(\\begin{cases} \\frac{i}{2} & \\text{for}\\: a = - \\frac{\\pi}{k} \\\\- \\frac{i}{2} & \\text{for}\\: a = \\frac{\\pi}{k} \\\\\\frac{1}{2 \\pi} - \\frac{\\pi e^{- i a k}}{a^{2} k^{2} - \\pi^{2}} & \\text{otherwise} \\end{cases}\\right) & \\text{for}\\: a = - \\frac{\\pi}{k} \\vee a = \\frac{\\pi}{k} \\\\- \\frac{\\sqrt{\\pi} \\sqrt{a} \\left(1 + e^{- i a k}\\right)}{a^{2} k^{2} - \\pi^{2}} & \\text{otherwise} \\end{cases}$$" ], "text/plain": [ "⎧ ⎛⎧ ⅈ -π ⎞ \n", "⎪ ⎜⎪ ─ for a = ───⎟ \n", "⎪ ⎜⎪ 2 k ⎟ \n", "⎪ ⎜⎪ ⎟ \n", "⎪ ⎜⎪ -ⅈ π ⎟ \n", "⎪ ⎜⎪ ─── for a = ─ ⎟ \n", "⎪√a⋅⎜⎨ 2 k ⎟ \n", "⎪ ⎜⎪ ⎟ \n", "⎪ ⎜⎪ -ⅈ⋅a⋅k ⎟ \n", "⎪ ⎜⎪ 1 π⋅ℯ ⎟ \n", "⎨ ⎜⎪─── - ────────── otherwise ⎟ \n", "⎪ ⎜⎪2⋅π 2 2 2 ⎟ \n", "⎪ ⎝⎩ a ⋅k - π ⎠ -π π\n", "⎪─────────────────────────────────── for a = ─── ∨ a = ─\n", "⎪ √π k k\n", "⎪ \n", "⎪ ⎛ -ⅈ⋅a⋅k⎞ \n", "⎪ -√π⋅√a⋅⎝1 + ℯ ⎠ \n", "⎪ ───────────────────── otherwise \n", "⎪ 2 2 2 \n", "⎩ a ⋅k - π " ] }, "execution_count": 19, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# give a little help, x ==> a*x\n", "phik = 1/sqrt(2*pi)*a*integrate(psi(1,x*a)*exp(-I*k*a*x), (x,0,1))\n", "simplify(phik)" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdwAAAA/BAMAAAC4KmReAAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAIruZZhDNMomrdlTvRN05z604AAAACXBIWXMAAA7EAAAOxAGVKw4bAAALoUlEQVRoBeVafawcVRU/s7uzsx9v941STArYrgqSAKar74m1KG4/tLHGslIlapC3iJCimLciEdHKm8YKRUO6RPCjCmzzCpSiyRI1qdXY12IDCMIL+odCQ5c0pWBIWWhrv2if55z7MftxZ3b78qj9OMnce+75/e7d+5s7M7t3zgKgWe+j8rSxy847baSi0Nie00ktzCydVnKHy6eTXOe/AWqtgLgIh6OqqyudUDaBoQQ12lTUiQP+KEsWLfcbj/iuwQtHVYfnpRPKzlQB6Gi22G+bW2G+9frExBthBIHNmJiYoJOfPKS59i5nn26kPe2Cs9P3hReOpvKSn6qxQ+zOMUYlCehktJ2Qb4yMK7RLHf1cwetCIfjD0xLMSu7V5GjeGdONK7QHX7rYPwsyGo5mPdX5WnaQbRjjFUWKFwHo8C1Wn97S9pEOr5g4M89B/jr9bgcuAjYsX8Je8rBm9I+nKroxqD2AVIdcM/rPpj7Cnc8VsfUYmYKA/DKG86XDt6Q/Dz8Y4MUd0XUR4ecHkABK32aoSe6QG/lkUfHHlIO1nqqOmdGNGldOH49HbD1GuqpAXW9Fjw5tSaFAt8Oc5zIVhrNYWqUgZib3b4aa5Pa5I3p9Ms0d9VTVYAFop9x4AbswW4+RrqpRdP019OjQdixyB+UdmKgBRPEwW7SwjIEmuZnBjSsUOeopD2s9VRULQKXcL88B+A0+CfHnWoLOPLP1GFKuPVqF1LPIygPcgSQ6lDm/frJizxkdz7z2yG0uBq3Za6oK66h3wk6O0cpm2Tvjx2d1sLLFbJGCTXKbOfFCU0tPVcUCUCnX2waJOzfc9FAZb8kSdmG2HkPKTWTGYOULty+6GQmfkgdWwvDevdWFC6D/kpkFjGxzI/z0uHQ32auS1F79HGA9xTKbYRZJizJ7d6WFFyA3W2ti6amqWAAq5Fq1i+H7sCaFaiFdx4LZegwp9+xoyaml3IdoyL/Ig3w2lIvP7ZFytm5j2/4rROsCCCtvAbHOVxbg0kCekNuPV5VvSO4bB7iCT88ufTHrthHFU/80n3obDuEAuWQRS6chx5I3BJ7x157ZTWe82F8FiMA9iMOIPMhnS1Yc/MLo9+JIxGXaB/5XmyD4c+NZ06dA3E3nqf5Ekc5Vqwl9GAtYXZKrTa+MigSg8mJO45XnVPg+witWnDr//perC5uKALdAnYbcJg/y2ZKV9BGUm4vz/PtyMNQ8HcVqq9Ne3KXQFrCxc4C1y73p02VixmtN/A65AaiUS2sSrfXTAHwxM1uPoeRej/Al0CCW4WLm1c3HPYKHPLiMagi/d6F+I5HsAxDjUcPv3RJx8Ys6XaIqUqVSmp6qCgSgUu5QzYOhYv84smM5LJitx5By7cMJsPcBr8M8JNGhDe/dxwC2l6XcGrw/7WqwzVlcVYG195JnH4ZI6VEVU3XsP8JTq5viM4L31JIaAfwNorh6qioQgEq5260CPAh9HrKjVSyYrceQcjONaZAqwROIwxx5kM8WycGyMiyFbJ6a8WpsR0IAhtK/sa4qMfx1WJz3Ooh/FxEl9wUp15pfJiBTFzCVyYuOrvZb5AWgUm5kLcCFkPSQmHUl2x9DXcwXFCBSg18gDi/Jg3wy56I3q9bsO8Yzs7bUsGmvvn36QwyYipk6mCAy/jxf963VrEED6FgN0ZJy7QcawJskF68wsl0CDijNqJTb1GcG+y1sJddnWXWcDh6Ts3nWaA8dY7mz8kSTclN2gzdJdsMScueHjmFG9U9Q3XeQvRZ25xYhWsCrHo/J2eqzr3S790yu+cNmYkXEjuhRuyE2SQseHufOkdAxwlH96Vauh7HgnUiiY3L2RMHpoWPfCigTLc6vAqwCrq7aJHFvO8dVQBGO6k6RIrtd2OchiY5JmbXnI730mznAaqF/jNjOyh8eqKlNkug+TVQBZTiqOs2VTig75uIXFh6Ts1juRauHnvMydf6I7XlBxltWbZJ66H3iUJLVW4O/pPxprnbqHrWGiyI272BNbZJE4CQp425yYQ9TvRruLiMtsbcH7qlD2VQ9dbR0V5J6qjvnFGKIH6WnkKBwKfZ1/HgOJ51CaHzrKSSmuxQ7KCPWvetJyVhVM0+bt0f24LqAiz3kl0wLZO5ujponMrXRmXXzeJxDuqo8UjTDj5jDFG2BnjfyzFEjdYqDfWK/1z6q2B4NgHj33o5C2gM44+LzO+IYIIjt8mewkjlAGcLKGvhquTPq4y2ec0O+pS0b5xRM0Z5ikaC3dpRDOgIPmwfBnJ5VgVXjBtRPDl5D6LVtlOmQwE9sj7aRRNPDtyEVE6J+55uwLjH18qadRtsj+w049wvtALcH8b1AEfrGDCBC0ngnrbb2ald/D8DTACqqqKbawpWNGOXaJnpvsSC5vD3a7BwpGIcZw7fH+yBpSpsjJCzDt4nIAeI1XhXRZ4twvwsqKqnGKhUo10jvLRgkl7dH//jMuoJpmEwJV7cOyf2dIEHCEmNUxwtU+nLvA7i/rKMCM5cjx1OueQoyGvXYyTbwNeHslTVYsv7j4KxdX8IoQ5yDi1fSo0X9Bjdd5S5UvN72XlcD2qH0X/r+V38FkbvWr8bHxJxRePFvG65efPT3w99MvjX3hgp8b9GDHNZdenSCVje0u1yykSrAObXM485L+Or8coA6diJI5OCGCn/CR2ishEE0X65zSGYGBWAqOf1Hv+cj58KIC8tcfCO959G99iHKnrwHH2B2DrIibOoeEpuU3GyNR9yB5XWQ3pz1wNm/rQgetgm6knNwm95d+44r0iYY9+WSDE6mYDTAXuH0H8kdg/4C4Occht1QhJ9BXx4/JVLJbIWECAeMEBQWclszgH4u0OThSOLNfRTngxkZnBpKPJB9cw19BkEiB/exhdR2Glj4iT9sDKgo1mYT6T9e3Rz0Vx18BbEUrkfuJndDIzNOz+vXL3BF2DxCYPTYVnfjHahGyaVvE7oyYRMG37Rvm8hLSOTgXp5dxnZmDAs0vbr8pxMVFVh7KdJ/Wm56z4IFC+EaZPV53lP4ZgrlfnF4jwi39+3SPja5Y/R+Fm/QGhb48+kdYOM/Enh195wJn6cHNUIiB2cfyVBbXbZa7gaw/EscCQYT6b9sPgERXl3+Al+KxOjO6odYrlW0L7tZvCE39A8JKblpl0jpdZxFa3JkV1rWKoDrjFEggi7gf5geoHsX0nhDZfZvBLhbQCIHF2vYh1MyB4hxJdfxIFUUmUGMBthjnP5TcuFJvEjwcsZr6aA7grdMpIKjJQscDhghKKzkfpYJfTAEsLiKlw05TUbLGsV2ZJyClNOzdgx8tI5P5gIsdx6HuPfnIj4zBcQ5uGTO3o8rwTlAjCu5MwYG7vOjCJiM03+JUgLitLqwbByW04MJ4GXIlnASufRWiLgcNnUPiUm5NgpAs2hPQE8bdpq60bJmPADxyj5Tx6sWH2MlfFb9dKULZ6y8F/61aB2dCoI4B4fpvVFcfc4BYlzJHZ6YwC22iiJgMkr/AfxyuTNry02zXgULP+O5ty5E5lpIFTBBuOUrK9bfyWFT77CYlBulqZJ9AECkENFpNlrWayFR9Di4qxlq9dugGRJVckVTRVu7HoeWlJvFdaK9fLSMiXRKIZIjTPxZlJb1JzA8UefgfAGZyjZI7RjUFkF0UVHTAG9rTMp9FwDv5X+EtxunEMkRxht9XtbfqRDeONptd1ohK9eOU9vKYfF/MSn3Ofr6X4Y35FEASiGywxMSG31eViRJs3PK66hboUixg4ABc9TEnOpYdC+P+Ed/L69TiM570fBpK/8sCnCj/+HTfLfda4HmtqPcNkeN1CkOxg7ygDf6e/n2FKKfB/Uv5imexHEcboI/Cx9Vai/fnkL086Drj+O03q6PWlWlkfEBo/by7SlEPw+KF/ZJb3H6Zcs/IbpJaX0IdWOfoLg9PE4z+2D36SWMT9nu/U4shsho96BF/Kw+sSY/idn8gPt0X7rujEl8+PHs8j/QDFlrpDbEYwAAAABJRU5ErkJggg==\n", "text/latex": [ "$$\\begin{cases} \\frac{a}{4 \\pi^{3}} \\left(-1 - \\left(-1 + i \\pi\\right) \\left(1 + i \\pi\\right)\\right) & \\text{for}\\: a = - \\frac{\\pi}{k} \\vee a = \\frac{\\pi}{k} \\\\\\frac{\\pi a}{\\left(a^{2} k^{2} - \\pi^{2}\\right)^{2}} \\left(2 \\cos{\\left (a k \\right )} + 2\\right) & \\text{otherwise} \\end{cases}$$" ], "text/plain": [ "⎧a⋅(-1 - (-1 + ⅈ⋅π)⋅(1 + ⅈ⋅π)) -π π\n", "⎪───────────────────────────── for a = ─── ∨ a = ─\n", "⎪ 3 k k\n", "⎪ 4⋅π \n", "⎪ \n", "⎨ π⋅a⋅(2⋅cos(a⋅k) + 2) \n", "⎪ ──────────────────── otherwise \n", "⎪ 2 \n", "⎪ ⎛ 2 2 2⎞ \n", "⎪ ⎝a ⋅k - π ⎠ \n", "⎩ " ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" } ], "source": [ "simplify(abs(phik)**2)" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAasAAAA/BAMAAACsv/q/AAAAMFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAIruZZhDNMomrdlTvRN05z604AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAI2klEQVRoBd1afWwcRxV/97G7d747e5WklRqgvhpClVDUq2rAKlAW9wPJqMnRQNWgqr5SQAX+8OFKhUJbX0QICRXEiFZgSxEXJS2u+89GCNGGijhNIgpN6SlUiI+QXhWlKRQlTp00aVLneO/NzN7ueS+3rnLnj5G88+Y3b968n9/u7Ny8BcAS+iBdF11Zs2LRUUJC0anFyAo6c4uS1kB+MdIy3l6MrCB2VtFa27deiS2r+0pBp7qyUqmYDZVDJyqVk6QVPyd1tWPG6YbDLrGC/rvpoBY/tiwWQDXyBavAavFTUjuSMSak2LJqdejxgHNpsH6tUOXX0XfrDMvGLstwV/y81OgotQ3XUW4aPDILy7n7hXIfVVfXG6gbNbT6zfAt2XraTcJ/HtxuKv0PodyOVShXb+DBlIiNE62kOfRKPeVm4bOgFbEeEl7EbIAI/vmXbvkkObRS3bs2+Ks2DX3wwAisHf80HP7zM/fQJKGebcV6k7Vn28W9RJGiiAEs/fFyrt2Xo3CUmw4td2er5BEwXoNwEaZ28sJ1yAzzYnzjcSr/qePFLwDGqSu1B64lphHWPj7sUZ9jWu0FMM7AcSD/tH0QKXuc8208CCIid1lwo68CgYJWB77sRKmr2IyOERiyAc7CN9h45DQgS2+5g0NxDISD3KebiQwJn8nCG15lAOYwiegcR2t3CeAtuJfdS6ahH5uNSqKgm6TzPGj1X+ZzTIuiNQX3MZX+Aqxh4eLPFpQHSUs7C9FJEi7+bOVIpcVlBJIZSJ1RtGz4cIID4efG6qJCd2wlSTsP4dxOhak6+l8hqWi1MXPV26K6C4w/gl6AIzyfXoweqb/xS5aUU3fnWPomrM4UFObULwlJ0frXHNAafHkVLN20FQ6++xFyRhv90RW/dvyrFTodgF7IWNaNfXs074BSCEkakpb2xCQE2/vXGmpVuze0PcBU0fTyDKlJWm3aJATb+wew3RSV0ffdZTY2HN/22z2kFRY7+J3aZHXv33j0HGi8YBkBZk1uAL4xdf4xGbIwWqD2/gGGt1wlNPWJIHN2Xs+soGOCtI1NPzhrO3v/IONbrRNNHw4FmLM3VeZb9dWMUA6dBmfvH2B4y1XixUfqL/5Vb0aNcoFaA1mB9b5jO3t/gcyvq27GPxfAo3vgMboLY+ooI8CQBaSyu7iAnA3satuLgVUXkmK7XDAWks8BfNW+bvpq3bk56AGe7/A5B/X9vi5kLvJ72nfAPAM1/8zCNHTSOrlwyxbbz/dbkBZv6LXuMf/b1Af17AI8DWcKn1FO3yUVOsv+5n4qNvR354eyfgpt9kz0aTeUKrpbSv6nEppdJ/3TJMZRsaG/HuRhao0fz9S0sZkoSOz2v5DgIUmAgQeTfv8N7HpgLynUltC+WiR4O+x/erMEY0Qb+ml4ys+WlpuJ3uFAfJCk10R53cfpH/g1R8kRaKvne9SZmnJ0Zi2oH/3egaFhw6QNvXYSPvQlbxe3ItZMsNuB+CdcNCPbr8i6jWjd5Cg5wlKUVjott1Dzj3F3NZL9aS35zv1Z3tDvMaYtHxNJH2xCYSlxX6tXxy6JM63kTFfpt7w/LWXwPdT+tF6uVIA39H/9/JjlY/VvM7FUTmGxCZa+KtseWrolUacy6FZrES1n0jpCD+JrNz0K2g3bS/C9viexGSnghZMY+nBiOwZlM7apeGjFhgUorpQuefbCiA0rlz9VxNHjmdT/nt647MSbt75d+tW+y/8Exo7xHMHuQYFk/2g1GoqRCL0Ih/KPmLBKS3NihuIgkhj91u/fj41bpREPrWhOolSJdMkRlLos4yQYe+BJ6Phkp9VfxGPosIkLye0AZYZdowKJ743WmwB6Dp4z3wAYeng/0GLWbgOIJMbuD9gPmAB/kPN7aCXKLq9EuoRoXQNwCnDf3WG3lzWI52Av2VwJh7JQYNg1KpAoaFUzJipzUr8ms3h2iifmYOCa11E4sQpJAJ26iiTGp8TP2CFSxOP0AzJVxUuGMUmoLCJdwtECOAdDFp7w6sP4BjyZ2AL/Jlrtb20TsBoStJ5dtBJjW9kw0hoq4vz40utI3zlATz3REkmM13vypHSILlg80UpNMCYuIl3i0NrdffPNeT2DfS/EknaBaGkbKxmGXaMCibOjlYR+tvotjBbOz9G6LqutQR66LZMY2jQmC7A0vglFuuQIFKCLo2XjKDyRxxh/OZLGeiVcBl88Q1qzLoqWTFXs2lxiEyosyh7hRVwn5F5qFT9bYO8FePX7RYhbeI5alEmM6KR2vg3H9cqxnmhF0xKlSqRL7lW0ksMAJaY19DPjnEm0cOxjDLtGBRIVrduE9gQeg1IyRYVF2SA8Qo2PMoKLN+7YYuZDebgvsZ9WLaClm5MY8bR2htaQG1iz5iaMFCVKlUiXbERaX6FoGS+hRVwgcPnJwetYdcFzWbiKYVKfTZG0NHSKimlM8IKgwsKgxFMFfDvlGTmM13U96zE3v7mU2jD+E2ymyjKJEbZh+xMIvIZ/VGS04tdcGEWHTQGKK6VLIP5Lc/Ddq3svDMMVPetT1z5v447Ygh0AgwdW/L1vrESwe1AgWdKKlKR2uAQimSLC4thAnPapP6Q4oHd5p8MRjjkSC6GybEtaonWlBJtdSVr0pQN/ooYn3ZxMkWFBWOSH6AT8UdArF9ghZx/LLXG5ySWjGLFkW211uVndD8veJlWS1uWYNaNP1GLZAnAyRYYFZ+X8EOHwm6oPV1VFJdED5ipLXLIjhtKO2FxB0jqI/176RG2gUgZKpqiwIFv+NoxwQCVVcFdTW7S0B1nhaclGOOuHNgGLnGKjz+J7VX2ippIpRhcWDIuTHxqszh+yqrKSlimB6qg3drLrs26VZsrRd9g6eux8oqaSKXLaan7IdRM206VLYrvCVnDJcD5RG5XJFGm+mh8avyQTtsbIliLNg8+784maTKao6av5IZ91QinNu1rnDRy9ahuVmjWhkfrc9msDJXLgusZexFq1jDV2JYCGyNwF8Pm2AMbmkcrD7EvjUDTWmCek/g/oNIWJvqklUgAAAABJRU5ErkJggg==\n", "text/latex": [ "$$\\begin{cases} \\frac{a}{4 \\pi} & \\text{for}\\: a = - \\frac{\\pi}{k} \\vee a = \\frac{\\pi}{k} \\\\\\frac{2 \\pi a}{\\left(a^{2} k^{2} - \\pi^{2}\\right)^{2}} \\left(\\cos{\\left (a k \\right )} + 1\\right) & \\text{otherwise} \\end{cases}$$" ], "text/plain": [ "⎧ a -π π\n", "⎪ ─── for a = ─── ∨ a = ─\n", "⎪ 4⋅π k k\n", "⎪ \n", "⎨2⋅π⋅a⋅(cos(a⋅k) + 1) \n", "⎪──────────────────── otherwise \n", "⎪ 2 \n", "⎪ ⎛ 2 2 2⎞ \n", "⎩ ⎝a ⋅k - π ⎠ " ] }, "execution_count": 21, "metadata": {}, "output_type": "execute_result" } ], "source": [ "simplify(_)" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEFCAYAAAAYKqc0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3XlclXXe//HXBYdNQRYFZZPFo6C4kB7U1CxNRa1wbCytyb2smbpr6r5Lp7nHmcb5jc3UzD0z2WZjZZtYVtKmZmZppiGaKy4goGwhIIsiy1m+vz9QxwUUFc51zuHzfDzIDud7Lt7n4vDm4jrf67o0pRRCCCFci5veAYQQQrQ+KXchhHBBUu5CCOGCpNyFEMIFSbkLIYQLknIXQggXJOUuhBAuSMpdCCFckEHvAMI5aZrmBgwAwoBaYL9SqkTfVEKIszQ5QlVcDU3TegDzgTFAFlAKeAO9gNPAq8BypZRNt5BCCCl3cXU0TVsBvAxsVhe9eDRNCwHuBSqUUsv1yCeEaCTlLlqNpmleSql6vXMIIeQNVXGNNE17/aLbvsAXOsURQlxEyl1cq0JN014G0DQtEPgSeEffSEKIs9pyt4zs73Fx8+fPp6qqih07drBgwQJ+/vOf6x1JCFegtcpCpNzF1fjoo4/O/b9SikWLFjF48GDGjx8PwJ133qlXNCFchZS7sL/Zs2c3e5+mabz++uvN3i+EaBEpdyGEcEGtUu7yhqq4Kn/605+oqKho9v6vv/6azz77zI6JhBBNkdMPiKvSr18/br/9dry9vRk4cCDBwcHU1dWRlZXFrl27GDNmDE8//bTeMYVo92S3jLgq06dP5+233+avf/0rISEhFBcX4+PjQ+/evRk5ciQ+Pj56RxTC2bXKbhnZchdXZceOHRw9epR3332XjRs3XnBfbW2tlLsQDkLKXVyVhx56iPHjx5OTk4PJZDr3eaUUmqaRk5OjYzohxFmyW0Zck1/+8pe8/PLLesdoMatNkVtWw4Hiag4UV9O1kzeRQT4kRQfh5+2hdzwhzidTIYVojlKK3QWVpKbnc6C4mkMlJ6kzN56F2F2D4x/8ns5TnsHdTWNAhD/DenRhmLEzA7sH4u3hrnN60c5JuQvRlMMlJ1n0WSbfZZcRGdiB8AAf+oR1ondoJ3qH+mEM8WXYkCG8sHItW46U8f2RcvYUVGG1KTwNbswZHs2dAyPo1dVP76ci2icpd+F81q5dy2OPPYbVauX+++9nwYIFF9xfX1/PjBkz2LFjB507d2blypVER0ezfv16FixYQENDA56enjz33HOMHj0aaHyTd9asWZy2aHSf+BDHvKLp6OnOE2N78YuhUXi4X3o4h8lkIiMj49zt6joz23NP8F1WGd8eLqWgspYnx8UxZ0QM7m6t8rMmREtJuQvnYrVa6dWrF+vXryciIoKkpCRWrFhBnz59zo156aWX2LNnD6+88gqpqal8/PHHrFy5kh9//JGuXbsSFhbGvn37SE5OprCwEICkIUO447Fn+TjbTEVNPSPD3fjH3LEEdfRsNsvF5X6+0pP1PP3xXtZnljA4Joi/3TWAyKAOrbsyhGieHKEqnEt6ejpGo5HY2Fg8PT2ZNm0aaWlpF4xJS0tj5syZAEyZMoUNGzaglOKGG24gLCwMgISEBOrq6qivr+fI0QLKe0xg1eEG4rt14tG40/js/+SyxX4lwX5eLJ0+iOem9OdAUTXj/7GJ1PRjyIVthDORchd2U1hYSGRk5LnbERER57a+mxpjMBjw9/envLz8gjEffvghN9xwA/XKjV+9n4ktchALxsfz3gNDSOoZdskyz1q6dCkmkwmTyURpaells2qaxl2mSNY+PpIBkQEs+Ggvz645SJ3Zei1PXQi7k3IXdtPUlq+maVc1Zv/+/cyfP5+//vMlfvHaD2SfMNPtyOfcOzTq3LiLl3nWvHnzyMjIICMjg+Dg4BZlDg/w4Z25Q/jNhHhe35LLI+/txGyVa38LxyflLuwmIiKC/Pz8c7cLCgrO7WppaozFYqGqqoqgoKBz4ydPnswLr73J0+t/4lDJSf56Rw8q931z2WVeLzc3jQdv7sHCOxL46sBxHl+5C6tNdtEIxyblLuwmKSmJrKwscnNzaWhoIDU1lZSUlAvGpKSksHz5cgBWrVrF6NGj0TSNyspKbrvtNhY8s5jnf7RytPw0b8xK4ufD4vHz82Pbtm0opXjrrbeYNGlSm+SfPjSKpyfG89meYuZ/uAebFLxwYHL6AWE3BoOBJUuWkJycjNVqZc6cOSQkJLBw4UJMJhMpKSnMnTuX6dOnYzQaCQoKIjU1FYAlS5aQU1LFn7bWYPME34y36PlfrwHw8ssvM2vWLGpra5kwYQITJkxos+cwb2QPTjdY+cdXWXTwdOeZlIRmdwMJoSeZCimcQkVNA7/9eC9bc8p5c/ZgBkQGXNfyLjcV8kqUUixec5Clm3J48OZYFoyPl4IXrUnOCinaB5tN8fj7u/g+u5z3Hhhy3cV+vTRN4zcT4jndYOHVb3Po5ufN7BExumYS4mJS7sLhvfRNNt8cKmXRz/piig7SOw7QWPB/TOmLm6bx9/WH6RHiy8heLZuBI4Q9yBuqwqFtyS7j7+sPMykxjPuGdNc7zgXc3DR+M6E33fy9eeL93ZSdqtc7khDnSLkLh/VTVR2Ppf5IbLAvf57czyH3a/t4uvPCvTdQXWfmyQ92y1GswmFIuQuHZLba+K8VOzndYOWV+wbS0ctx9yDGd+vEbyf2ZuOhUt7Ykqd3HCEAKXfhoJ5bd4jteRUsvrMfxhDHP/XujBujGNM7hGfXHGR/UZXecYSQcheOZ3NWKUs35TB9aBSTEsP1jtMimqbx1ykDCOzowaMrfuR0g0XvSKKdk3IXDqXqtJlfp+5iYt9Q/vf23nrHuSpBHT35v7sTySmrYdFnmXrHEe2clLtwKM+uPUhlrZlfjeqBl8H5Lnc3zNiFh27uwYr0fL7YW6x3HNGOSbkLh7E97wQr0o8xZ3g0fcP99Y5zzZ4Y24uhsZ1ZkX6McpkeKXQi5S4cQr3Fym8+2kt4gA+Pj+2ld5zr4uHuxh9T+rD1SDl/WXtQ7ziinZJyFw7h1W9zyD5+ij/9rC8dPB132mNL9erWibkjYng/o4AdRyv0jiPaISl3obuc0lMs2ZjN7f1DGRUfonecVvPorT3p1smb363eh0Uu8CHsTMpd6EopxdMf78Xb4MbCO/pc+QFOpKOXgd/d3ofM4mre/eGY3nFEOyPlLnT1wY4CtuWcYMGE3oT4eesdp9VN7NeNEcYuPP/lIUpPypurwn6k3IVuquvMfLKrkNHxIUxLirzyA5yQpmk8MymBOrOVxWsO6B1HtCNS7kI3r3xzhC1HynlibC/c3BzvpGCtpUewLw/cFMtHOwtJzz2hdxzRTki5C10UVday7LtcJieGO/Wc9pZ6ZLSRMH9v3t6WJxfXFnYh5S508ff1h1HAE+Oce057S3XwNPCbifF8uruYT3YX6h1HtANS7sLuDhRX8+HOAmYPiyYisIPecezmtn5h9AntxN/XH6bBIlMjRduSchd2t3jNQTp5e/CrW4x6R7ErNzeNJ8fHkX+ilpXbZWqkaFtS7sKuNmeVsulwKf812oh/Bw+949jdLb2CGRwdxL++zpbTAos2JeUu7MZmUyz+4iARgT5MvzFK7zi60DSNp8bHUXqynje/z9M7jnBhUu7CblbvKiSzuJonk+Oc8nS+rcUUHcTo+BBe+eYIVafNescRLkrKXdhFndnK8+sO0S/cnzv6h+kdR3f/My6O6joLr246oncU4aKk3IVdvL89n8igDiyYEO/SByy1VJ+wTqQMCOONLXkcr67TO45wQVLuos3Vma28sDEbdzeN4cYuesdxGE+M7YXZauOFr7P1jiJckJS7aHOp6ccoPVnPo7f21DuKQ4nu0pGpSZGsSD/GsfLTescRLkbKXbSpOrOVl789wuCYIIbGdtY7jsN59NaeGNw1/u+rw3pHES5Gyl20qQ92FFBSXc9jstXepK6dvJk5LJrVPxZyuOSk3nGEC5FyF22mwWLj5Y3ZDIoKZFgP2WpvzkMjY0kI78S/NmTpHUW4ECl30WY+3FlAUVUdj97aE02TGTLNCezoxQhjMJ/vLSan9JTecYSLkHIXbcJstfHixmwGRAYwsqfMkLmSuSNi8HR345VvZd67aB1S7qJNfPxjIQUVtTx2q1G22lsg2M+LaUmRfLSzkMLKWr3jCBcg5S5aneXMVnvf8E6MigvRO47TmHdzDwBe25SjcxLhCqTcRav7ZHcRR8tP8+ho2dd+NcIDfJh8Qzip249Rdkoupi2uj5S7aFVWm2LJ19n0Du3E2D5d9Y7jdB66pQf1Fhuvf5erdxTh5KTcRav6cv9PeHm48V+je8hW+zXoEezLxH6hvL31KFW1csZIce2k3EWrUUqxZGM29RYbyX266R3Haf3qlh6crLfw9tY8vaMIJyblLlrN90fK2V9UzbybYnF3b/qltXbtWuLi4jAajTz77LOX3F9fX8/UqVMxGo0MGTKEvLw8AMrLyxk1ahS+vr488sgjFzzmlltuIS4ujsTERBITEzl+/HirPzd7SgjzZ3R8CK9vyZOrNYlrJuUuWs2rm3Lo4uvFz24Ib/J+q9XKww8/zJo1a8jMzGTFihVkZmZeMGbZsmUEBgaSnZ3N448/zvz58wHw9vZm0aJFPP/8800u+91332XXrl3s2rWLkBDnn6Hz8KgenKhpYEV6vt5RhJOSchet4kBxNZsOlzJ7eDTeHk1fZSk9PR2j0UhsbCyenp5MmzaNtLS0C8akpaUxc+ZMAKZMmcKGDRtQStGxY0dGjBiBt7d3mz8XRzAoKoifDwxn7b5izFab3nGEE5JyF61i6aYcOni6c9+Q5q+NWlhYSGRk5LnbERERFBYWNjvGYDDg7+9PeXn5Fb/+7NmzSUxMZNGiRSilms64dCkmkwmTyURpaWlLnpaubusfyva8Cj7fU6x3FOGEpNzFdSuqrOXT3UVMS+qOfwePZsc1VboXz6hpyZiLvfvuu+zdu5fNmzezefNm3n777SbHzZs3j4yMDDIyMggODr7sMh3BLb1CMIb4snRTTrO/sIRojpS7uG6vf5eLAuaMiL7suIiICPLz/7MPuaCggLCwsGbHWCwWqqqqCAoKuuxyw8Mb9/H7+flx7733kp6efvVPwgG5uWk8cFMMmcXVfH/kyn+9CHE+KXdxXapqzaxIP8bt/UOJCOxw2bFJSUlkZWWRm5tLQ0MDqamppKSkXDAmJSWF5cuXA7Bq1SpGjx592S13i8VCWVkZAGazmc8++4y+ffte57NyHJMSw+ni68lrm+WUBOLqGPQOIJzbez8co6bByryRsVccazAYWLJkCcnJyVitVubMmUNCQgILFy7EZDKRkpLC3LlzmT59OkajkaCgIFJTU889Pjo6murqahoaGli9ejVffvklUVFRJCcnYzabsVqtjBkzhgceeKAtn7JdeXu4M/PGaP62/jCHfjpJXDc/vSMJJ6G14b482Uno4uotVm76y0Z6dfXjnfuH6B3nqphMJjIyMvSO0SIVNQ3c+OwG7ugfxnN3DdA7jmh7rXJot+yWEdcsbVcRx0/Wt2irXVy7wI6e3DUosnF9V9fpHUc4CSl3cU1stsaTW/UO7cRNcjGONjd3RAxmm43lW/P0jiKchJS7uCabssqoPN3Aw7fICcLsIbpLR5L7dOOdbcfklASiRaTcxTV5fUseNgXjEuQEYfbywMgYqmrNpO0q0juKcAJS7uKqZZWcZNPhUmbcGIWnQV5C9jIoKogxvUN49dsj2GwyX0Fcnvxkiqv2+pY8vAxu3HuZUw2ItjEpMZy88tNsPOTcZ74UbU/KXVyVipoGPtpZwJ0Dwwnq6Kl3nHZnfN9uhPp78/oWuVKTuDwpd3FV3ks/Rr3FxuzhMXpHaZc83N2YcWM0W7LLOfhTtd5xhAOTchctZrbaeGtrHjf17EKvrnKkpF7uGRyJt4cbb27J0zuKcGBS7qLFvthbTEl1PXNGyFa7ngI6eHLnwAg+/rGQEzUNescRDkrKXbSIUopl3+USG9yRm3s6/ulyXd3sYdHUW2ysSD+mdxThoKTcRYvsPFbBnoIq5gyPwc1NDlrSW8+uftzUswtvbc2jwSJXahKXknIXLbLsu1z8fTy4c2DT10cV9jdnRAwl1fWs2SdXahKXknIXV1RQcZqyUw3MGR5NB085S7SjuLlnMLFdOjZeLEWu1CQuIuUurujtbUfZcbSCKabIKw8WduPmpjF7eDS7C6rYeaxS7zjCwUi5i8uqbbCSmp5PckJXwgN89I4jLnLnwAj8vA1yUJO4hJS7uKzVuwqpqjUza5hMf3REHb0M3DO4O2v3/URRZa3ecYQDkXIXzVJK8caWXPqEdiIpOlDvOKIZ04d2Z3B0ECu35195sGg3pNxFs7YeKedwySlmDY+Wc7Y7sMigjgR08GD51jzqzFa94wgHIeUumvXG93kEdfQkZUCY3lHEFcwcFk3laTNpuwr1jiIchJS7aFL+idN8daCEewd3x9vDXe844gqGxAQR382PN7bkybRIAUi5i2a8tTUPN03jvqFyznZnoGmN0yIP/nSSH3JP6B1HOAApd3GJmnoLqdvzmdC3G938vfWOI1poUmI4AR085GyRApByF0346MdCTtZZmD08Wu8o4ip4e7hzz+DufJn5EwUVp/WOI3Qm5S4uoJTizS259Av3Z2B3mf7obO4bGoWmaby97ajeUYTOpNzFBb7LLuNIaQ2zhsn0R2cUHuBDckJXUtPzqW2QaZHtmZS7uMCbW/Lo4uvJ7QNC9Y4irtGsYTFU1ZpZLdMi2zUpd3FOXlkNxVV1zBsZi5dBpj86q6ToQPqEduJNmRbZrkm5i3OWb83jcMlJfpYo52x3ZpqmMWt4NIdKTrI1p1zvOEInUu4CgFP1Fj7IKOC2/qGEdJLpj84uZUAYQR09Wf59nt5RhE6k3AUAH+4o4FS9hdnD5eyPrqBxWmQk6zNLyD8h0yLbIyl3gc2mWP59HomRASRGBugdR7SS+4ZGkRDWSd5Ybaek3AWbskrJKauRg5ZcTKi/D907d+S1TTmcbrDoHUfYmZS74I0teYT4eTGhr0x/dDVzhkdTXWfhw52y9d7eSLm3c0dKT/Ht4VLuGxqFp0FeDq5mYPdA+kf48+aWXGw2mRbZnshPczv31vd5eLq7cc/g7npHEW3g7Nkij5TW8F12md5xhB1Jubdj1XVmVu0o4PYBoQT7eekdR7SR2/qFEeznxRtyEe12Rcq9HVuVUUBNg5XZcvFrl+ZpcOMXQ7qz8VApOaWn9I4j7ETKvZ2y2hTLt+ZhigqkX4S/3nFEG/vFkCg83d3koKZ2RMq9nfr28HGOlp9mlp2nP65du5a4uDiMRiPPPvvsJffX19czdepUjEYjQ4YMIS8vD4Dy8nJGjRqFr68vjzzyyAWP2bFjB/369cNoNPLoo4/K+VSaEOznxe0DQlm1o4DqOrPecYQdSLm3U69tymGEsQvJCd3s9jWtVisPP/wwa9asITMzkxUrVpCZmXnBmGXLlhEYGEh2djaPP/448+fPB8Db25tFixbx/PPPX7LcX/7ylyxdupSsrCyysrJYu3atXZ6Ps5kzPIaaBivvb8/XO4qwAyn3duhAcTVbc04w3NgFD3f7vQTS09MxGo3Exsbi6enJtGnTSEtLu2BMWloaM2fOBGDKlCls2LABpRQdO3ZkxIgReHtfeN6b4uJiqqurufHGG9E0jRkzZrB69Wq7PSdn0jfcn6ToQJZvzcMq0yJdnpR7O/T6d7n4nDn3iD0VFhYSGfmfrxkREUFhYWGzYwwGA/7+/pSXN39mw8LCQiIiIi67TPEfs4fHkH+ilg0HSvSOItqYlHs7U3aqnrRdRfx8UDgBHTzt+rWb2hd+8dWeWjLmWscvXboUk8mEyWSitLT0SnFd0rg+XQkP8OENuYi2y5Nyb2fe2XaUBqtNl7M/RkREkJ//n/29BQUFhIWFNTvGYrFQVVVFUFDQZZdZUFBw2WWeNW/ePDIyMsjIyCA4OPh6norTMri7Mf3GKHYeO8HBn6r1jiPakJR7O1JvsfLOtqOMigumR7Cv3b9+UlISWVlZ5Obm0tDQQGpqKikpKReMSUlJYfny5QCsWrWK0aNHX3bLPTQ0FD8/P7Zt24ZSirfeeotJkya16fNwdlOTIgns4MWyzXJQkysz6B1A2M8nu4ooO9XA3BGxunx9g8HAkiVLSE5Oxmq1MmfOHBISEli4cCEmk4mUlBTmzp3L9OnTMRqNBAUFkZqaeu7x0dHRVFdX09DQwOrVq/nyyy/p06cPL7/8MrNmzaK2tpYJEyYwYcIEXZ6fswjs4Mm4MxfRfnJ8HCF+cnEWV6S14ZxgeTvegSilmPiv77DZFGt/fdNlt4bbA5PJREZGht4xdJNXVsOov33DI6OM/Pe4OL3jiAu1yg+n7JZpJ7bmlHOguJo5I6LbfbELiO7SkbG9u/LOtqPUNlj1jiPagJR7O/H6d3kEdfRkklz8WpzxwMhYKk6bWbWz4MqDhdORcm8H8spq2HCwhPuGdMfbw13vOMJBmKICGRDhz+vfybneXZGUezuwfGseBjeN+26M0juKcCCapnH/TbHkltWw4eBxveOIVibl7uIqTzew+sdC5oyIkVkR4hIT+nYjPMCHf2/O0TuKaGVS7i7u3R+OUXHazOQbZF+7uJTB3Y3Zw6P5IfcEewuq9I4jWpGUuwurt1h5Y0seI3sFE9+tk95xhIOamhSJr5eB12Tr3aVIubuw1T8WUnaqngdH6nPQknAOft4eTEuK5PO9xRRV1uodR7QSKXcXZbMplm7KoU9oJ4b16Kx3HOHgzl605U25UpPLkHJ3URsPHedIaQ0P3hwrBy2JK4oI7MCEvt1Y8cMxTsqVmlyClLuLenVTDmH+3kzsF6p3FOEk7r8plpP1Ft7PkIOaXIGUuwvalV9Jeu4J5oyIseuVloRzS4wM4LZ+oWzOKsVstekdR1wn+cl3Qa9tysHP28C0wd31jiKczJ0Dw/nmUCmf7i7SO4q4TlLuLuZoWQ0bDx3nF0Oi8PWSMzqLqzM6PoT4bn68/M0ROSWBk5NydzFLN+fg52Vg1nA51YC4epqm8ctbepB1/BRfyXVWnZqUuwspqa7jg4wCbu3TlW6dfPSOI5zUbf1CiQzy4aVvjjR5jVrhHKTcXchrm3KwKsVDI3voHUU4MYO7G/NG9mBXfiXbck7oHUdcIyl3F1FR08C7PxwjZUAY3Tt30DuOcHJ3DYqgi68XL32TrXcUcY2k3F3EG1tyqTVb+dUtstUurp+3hztzR8SwOatMTijmpKTcXcDJOjNvfp9HckJXenb10zuOcBG/GNodPy8DL38rW+/OSMrdBbyz7RjVdRYeHmXUO4pwIZ28PZh+YxRr9v3EkdJTescRV0nK3cnVma0s+y6Hm3p2oX9EgN5xhIuZMyIGT3c3Xv32iN5RxFWScndyK7fnU3aqgUdkq120gS6+XkxNimR73gkKK07rHUdcBSl3J9ZgsfHqt0cwRQUyOCZI7zjCRT04MpbCyjpelq13pyLl7sRW7yqkqKqOh0cb5bS+os2EB3bg5wMjWLk9Xy7m4USk3J2UxWrjy/0/0TesE7f0CtY7jnBxD49qnGL74kaZOeMspNyd1OpdRXx14DiPjukpW+2izUUEduAuUyTvZ+RTKFvvTkHK3QmZrTb+tSGLhLBOjO3dVe84op04O9VWtt6dg5S7E1q1o4BjJ07zxNhestUu7CY8wIe7TZF8kJFPgcyccXhS7k6m3mLlhQ1ZJEYGMDo+RO84op35z9a7zJxxdFLuTmbl9nyKqur473Gy1S7sLyzAh6lJsvXuDKTcnUid2cqSr7MZHB3ECGMXveOIdupXtxhx0zTZ9+7gpNydyDvbjnL8ZD1PyFa70FFYgA8P3hzLzmOV5JbV6B1HNEPK3UnU1Ft4+ZsjDDd2ZmhsZ73jiHZu+o1R5J84zfPrDukdRTRDyt1JvL89n8paM0+MjdM7ihCE+Hlz/02xfL63mF35lXrHEU2QcncC5afq+fv6w6QMCGVQVKDecYQAYN7IWLr4erL4iwNyrVUHJOXuBF74OpvTZqucr104FF8vA4/e2pMfck+w8dBxveOIi0i5O7ic0lO8s+0o05IiMYbIVZaEY7lncHeiO3fgL2sOYbXJ1rsjkXJ3cH9ZexAvgxu/HtNL7yhCXMLD3Y0nk+M5VHKSD3cW6B1HnEfK3YGl555g3f4SHrq5B8F+XnrHaRVr164lLi4Oo9HIs88+e8n99fX1TJ06FaPRyJAhQ8jLyzt33+LFizEajcTFxbFu3bpzn4+OjqZfv34kJiZiMpns8TTEeSb268aAyAD+b/1h6sxWveOIM6TcHZRSiv/3xQG6dvLi/pti9Y7TKqxWKw8//DBr1qwhMzOTFStWkJmZecGYZcuWERgYSHZ2No8//jjz588HIDMzk9TUVPbv38/atWv51a9+hdX6nyLZuHEju3btIiMjw67PSYCmafxmQjzFVXW8sSVP7zjiDCl3B/XZnmJ251fyP+Pi8PF01ztOq0hPT8doNBIbG4unpyfTpk0jLS3tgjFpaWnMnDkTgClTprBhwwaUUqSlpTFt2jS8vLyIiYnBaDSSnp6ux9MQTRga25nR8SG89E02FTUNescRSLk7pHqLlb+sPUh8Nz/uHBihd5xWU1hYSGRk5LnbERERFBYWNjvGYDDg7+9PeXn5ZR+raRrjxo1j0KBBLF261A7PRDRl/vh4+ob7888NWXpHEYBB7wDiUu9vz6egopa35w7G3c11TjPQ1Fzoi0+j0NyYyz12y5YthIWFcfz4ccaOHUt8fDwjR468ZPzSpUvPlX9paek1PQfRvLhufvQK8eWtrXlMGRRB33B/vSO1a7Ll7mBKquv4y9pD3G2K4KaernX5vIiICPLz88/dLigoICwsrNkxFouFqqoqgoKCLvvYs/+GhIQwefLkZnfXzJs3j4yMDDIyMggOdq116yieGBdHYAdPFqbtwyZTI3Ul5e5g/vT5ARqsNpc8YCkpKYmsrCzIIrmZAAAUrklEQVRyc3NpaGggNTWVlJSUC8akpKSwfPlyAFatWsXo0aPRNI2UlBRSU1Opr68nNzeXrKwsBg8eTE1NDSdPngSgpqaGL7/8kr59+9r9uYlG/j4eLJgQz85jlaySqZG6kt0yDmRLdhmf7i7isVt7EtW5o95xWp3BYGDJkiUkJydjtVqZM2cOCQkJLFy4EJPJREpKCnPnzmX69OkYjUaCgoJITU0FICEhgbvvvps+ffpgMBh48cUXcXd3p6SkhMmTJwONW/r33nsv48eP1/Nptns/HxhB6vZ8/rLmIMl9uuHfwUPvSO2S1obnhJC/ya5Cg8XG+H9uwmpTrPv1SLw9XGOGjKMymUwybbIN7S+q4o4XvuO+oVH8cZL8JXWVWuWNNtkt4yBe25xDTmkNf0hJkGIXTi8hzJ/pQ6N4Z9tR9hVW6R2nXZJydwAFFad54esskhO6MipOrosqXMPZN1d/J2+u6kLK3QH88dNMNDQW3pGgdxQhWo2/jwe/mdibH49V8sGO/Cs/QLQqKXedrdlbTHFVHU8mxxEe4KN3HCFa1Z03hHO3KYK/fXmY49V1esdpV6TcdVR+qp7/Xb0PhWL6jVF6xxGi1bm5aTx4cw+qas08/fE+uaiHHUm562jhJ/uprjPzt7sS8XCXb4VwTT2CfXkyOY6vDpSwelfhlR8gWoU0ik4+31PM53uK+fWYXsR1k4twCNc2e3gMpqhAfp+2nxLZPWMXUu46KDtVz+/S9tE/wp8HR7rG6XyFuBx3N43n7hpAg9XG0x/tld0zdiDlbmdKKX63eh+n6iw8f9cADLI7RrQTMV068mRyPBsOHufDnbJ7pq1Js9jZF3t/Ys2+n/j12J706iq7Y0T7MntYNIOjg3jmk30UVJzWO45Lk3K3o7yyGuZ/uIdJiWHMc5GrKwlxNdzcNP46pT/xof48sXI3FqtN70guS8rdTurMVh5+byfubhpPJsfJ7hjRbkV36ch9Q7uTnneC5788rHcclyUNYyd//uIA+4uq+dtdA4gI7KB3HCF0NSkxnF8M6c4r3x5hw4ESveO4JCl3O/h8TzFvbT3KAzfFMKZPV73jCOEQfnd7H/qGd+KJ93eTf0L2v7c2Kfc2dnY/+w3dA3hqfLzecYRwGN4e7rx07yBsSvHIezupt1j1juRSpNzbUG2Dhf9a8SPubhpL7h0oR6EKcZHunTvw3JQB7C6oYvEXB/WO41KkbdqIzaZ46sO9WKw2/n73ADkpmBDNGN+3G3NHxPDm93l8trtI7zguQy6z10b+76vDfLq7iKfGx3Frb9nPLsTlLJgQT1Wtmd+u3kd4oA83dA/UO5LTky33NvBBRj4vfJ3NVFMkv7y5h95xhHB4Hu5uLJgQj7+PB/cvz+BoeY3ekZyelHsr23qknKc/3stwY2f+NLkvmtYql0MUwuV18fXizdlJWJVi1hvbOVHToHckpybl3oqOlJ7ioXd2ENW5Iy/9YpC8gSrEVYoN9uXfM0wUVtbywFsZ1JllBs21kvZpJUWVtfz58wN4uGu8MSsJfx8PvSMJ4ZRM0UH8Y2oiO49V8PjKXXKKgmsk5d4KiiprmbZ0G+m5J1g+ezCRQXIEqhDXY2K/UP53Ym9Kqut46sM9WOUC21dNZstcp6LKWu55bRsVNQ28ff8QEsL99Y4khEuYe1MsNQ1W/r7+MFab4m9yiuyrIuV+HYqrGov9xKkG3po7mMTIAL0jCeFSHr21Z+OFPtYdwmpT/GNqohR8C0m5X6PiqsZdMWeLXeblCtE2Hh5lxOCmsXjNQWxK8c9pN8hkhRaQcr8GB3+qZv6qPaBguRS7EG3uwZt74O6m8afPDxDYYT+/va03HTylvi5H1s5V2njwOI+8txNfbwOvz0oiIUz2sQthD/ffFEtQR0+eXLWH3QWVvDbDRKi/nNajOfK3TQsppXhjSy5zl28nuktH0h4eIcUuhJ3dOTCC12YMIq/sNJOWbGF3fqXekRyWlHsLNFis/G71Pp75NJMxvbvywUM30s3fW+9YQrRLo+O78uEvh+FpcOPuV7fy6S452VhTNKXabP6oS0xMzSk9xeMrd+Hj6c6AyADmJ8fj5ianFHB2JpOJjIwMvWOI61B2qp6H3s7AaoOE8E78dmIffDzd9Y7VGlqlYGTLvRlKKd774Ri3/es78spPM2NoFL+Z0FuKXQgH0cXXixUP3IgpOpB3th3j9hc2s6+wSu9YDkO23JtQerKepz/ey/rMEoYbO/O3uxJlN4yLkS1317Ilu4z/fn835TX1PDG2Fw/cFOvM8+FbZQtSyv08FquNd384xr8351BnsfHgyFjmDI+RrXUXJOXueipPN/D0x3spP9VAVa2ZhXf0YViPLnrHuhZS7q1pS3YZz3y6n8MlpxjWozOLJiXQI8RP71iijUi5uyalFF/u/4lFnx+goKKWCX278fTE3s52vicp99awK7+S97fn8176MSKDfPjtxD4kJ3SV87C7OCl311ZntrJ0Uw4vfZONzaaYN7IH02+Momsnp9i9KuV+rWw2xcZDx3l1Uw7puSfoHtSBewZHMnt4DN4eLvFuu7gCKff2oaiylje/z+Xfm3Nxd9NIGRDOAyNjiO/WSe9olyPlfrWq68ys2/cTr27KIfv4KcIDfJgzIoapSZH4esnBuu2JlHv7cqz8NK9vyWXl9nxqzVZG9gpm7ogYRhi74O5476lJubdETb2Frw6U8OnuYjYdLiU+1A+LVfHgzbFM7BcqJyBqp6Tc26fK0w28+8Mx3tiSR3iAN4WVtYxL6MZt/UIZEhPkKDNspNybYrMpDpWc5IeccrKOn+LDnQXUmW107eTFbf3CuGNAKImRAbJPXSdr167lsccew2q1cv/997NgwYIL7q+vr2fGjBns2LGDzp07s3LlSqKjowFYvHgxy5Ytw93dnX/9618kJye3aJlNkXJv3+otVr45WMone4r4+sBxas1WAjt4MKFvKMOMnbmheyBh/t569YSUu9WmKKg4zeGSUxw5fpKMo5VszztBVa0ZgIHdA0gI8+eOAWGYogJlSqPOrFYrvXr1Yv369URERJCUlMSKFSvo06fPuTEvvfQSe/bs4ZVXXiE1NZWPP/6YlStXkpmZyT333EN6ejpFRUWMGTOGw4cPA1xxmU2Rchdn1TZY+fbwcb7Y+xPlp+rZcqQcgBA/L27oHsAN3QO5ITKAXl19CezoZY9IrVJUbbajuaKmgVqzFU+DW+OHuxue7hpubs3/2WO1KcxWGw0WK7VmGyfrzFTVWqiuM3Oy1kJhZS0l1XUUV9VSb7axLbecOnPj9RU1oHdoJ5ITujIkpjNDYoOICHSq6U8uLz09HaPRSGxsLADTpk0jLS3tgiJOS0vjD3/4AwBTpkzhkUceQSlFWloa06ZNw8vLi5iYGIxGI+np6QBXXKYQl+Pj6c74vqGM7xuK2WrjYPFJfsyv4Mdjlfx4rIJ1+0vwcNOwKIWvp4HwQB/CAnwID/AhPNCHyAAfAjp64udtwNfLgK+3gQ6e7ni6u2Nw05rdqFRKndd5inqrlQaLrdV6q83K/cWN2fz7u9wLPtcntBOZxdW4aeDupqFpGh083DhttmGx2jh7mcTuQT4cO1F7wWMTI/3ZlV+Fn5eBbv7eJMUE8YuQKHp19aVnVz+MIb508paLUjuywsJCIiMjz92OiIjghx9+aHaMwWDA39+f8vJyCgsLGTp06AWPLSwsBLjiMoVoKQ93N/pF+NMvwp8ZNzZ+rvxUPZlF1Rz86SSFlbUUVNRSWFlLRt4Jquss9O7mx4GfTl6wnMHRQaTnnQBA08DgpuHuphHbxZfs0lNYbercdWF9PNypNVvPPTbv2dta5bm02W4ZY+9+yoaGovE3lE01PknbmSekzv3nTBANNDQ0Ddy0xn/d3TTctMaV4u6m4eHeeLs1lZaWEhwc3KrLbAuukLOiooLq6mqioqIAKC8vp6amhu7du58bs3//fnr27ImnpycAe/fupXfv3hQVFeHr60tQUBAAeXl5+Ps3nnL5Sss8P1tZWRkAZrOZ/v37t9Kzbjuu8H13FG2R0aYUZovCYrNhUwqrDaxKoQEWmwIFCnWu7zQN1Jl/ATRNQ+NM/2kabkDOwb37lVJ9rzucUqpNPgYNGqScgeRsXZfL+f3336tx48adu/3nP/9Z/fnPf75gzLhx49T333+vlFLKbDarzp07K5vNdsnYs+NassyrzelIJGfrcYaMSikFZKhW6GCHmPcj2oekpCSysrLIzc2loaGB1NRUUlJSLhiTkpLC8uXLAVi1ahWjR49G0zRSUlJITU2lvr6e3NxcsrKyGDx4cIuWKUR7JEfuCLsxGAwsWbKE5ORkrFYrc+bMISEhgYULF2IymUhJSWHu3LlMnz4do9FIUFAQqampACQkJHD33XfTp08fDAYDL774Iu7ujUcTN7VMIdq91tj8b+rj1VdfbcM/XFqP5GxdkrN1Sc7W4wwZlVIKmKdaoYOdep67EEK4ILkSkxBCiKZdV7lrmnaXpmn7NU2zaZpmOv++xYsXYzQaiYuLY926dU0+Pjc3lyFDhtCzZ0+mTp1KQ0PD9cRpkalTp5KYmEhiYiLR0dEkJiY2OS46Opp+/fqRmJiIyWRqckxb+sMf/kB4ePi5rF988UWT49auXUtcXBxGo5Fnn33WzinhySefJD4+nv79+zN58mQqK5u+Gr1e6/NK66e+vp6pU6diNBoZMmQIeXl5dssGkJ+fz6hRo+jduzcJCQn885//vGTMN998g7+//7nXwh//+Ee7ZjzrSt9DpRSPPvooRqOR/v37s3PnTrtnPHTo0Ln1lJiYSKdOnfjHP/5xwRi91uecOXMICQmhb9//zHI8ceIEY8eOpWfPnowdO5aKioomH6tp2kxN07LOfMxs0Re8nn06QG8gDvgGMJ33+T79+/dXdXV1KicnR8XGxiqLxXLJvqW77rpLrVixQiml1IMPPqheeuml1tx1dUVPPPGEeuaZZ5q8LyoqSpWWlto1z/l+//vfq+eee+6yYywWi4qNjVVHjhxR9fX1qn///mr//v12Stho3bp1ymw2K6WUeuqpp9RTTz3V5Dg91mdL1s+LL76oHnzwQaWUUitWrFB33323XTMWFRWpHTt2KKWUqq6uVj179rwk48aNG9Vtt91m11xNudL38PPPP1fjx49XNptNbd26VQ0ePNiO6S5lsVhU165dVV5e3gWf12t9fvvtt2rHjh0qISHh3OeefPJJtXjxYqWUUosXLz7783NxzwYBOWf+DTzz/4EXj7v447q23JVSB5RSh5q4a1Jzh4qf/0vl66+/ZsqUKQDMnDmT1atXX0+cq6KU4v333+eee+6x29dsbecfzu/p6Xnu0Ht7GjduHAZD46SroUOHUlBQYNevfzktWT9paWnMnNm4ITRlyhQ2bNhw9gfKLkJDQxk4cCAAfn5+9O7d+9yRt84mLS2NGTNmoGkaQ4cOpbKykuLiYt3ybNiwgR49epw7wE1vI0eOPHcQ3lnnv/4u04HJwHql1AmlVAWwHhh/xS94pfZvyQeXbrkvAe477/YyYMpFj+kCZJ93OxLY1xp5Wph5JJc5WADIBXYCO2ild6+vMt8fgDxgD/A6TfymBqYA/z7v9nRgib2znvf1Pz3/+673+mzJ+gH2ARHn3T4CdNFp/UUDx4BOF33+FqAc2A2sARJ0ynfZ7yHwGTDivNsbzu8FHfK+DjzSxOd1W59nvsf7zrtdedH9FU085n+A/z3v9u+A/7nS17riPHdN074CujVx12+VUs1tJjb1bu/Fm0MtGXNNWpj5HmDFZRYzXClVpGlaCLBe07SDSqlNrZGvJTmBl4FFNK6TRcDfgDkXL6KJx7b6ZmdL1qemab8FLMC7zSymzddnE3R9HV4NTdN8gQ+BXyulqi+6eycQpZQ6pWnaRGA10NPeGbny99Ah1iWApmmeQArwmybudpT12VLXtF6vWO5KqTHXEKaAxi3xsyKAoovGlAEBmqYZlFKWZsZckytl1jTNANwJDLrMMorO/Htc07SPgcFAq5ZRS9etpmmv0bhVdLGWrOfr1oL1ORO4HbhVndm0aGIZbb4+m9CS9XN2TMGZ14U/cKKNc11A0zQPGov9XaXURxfff37ZK6W+0DTtJU3TuiilyuyZswXfQ7u8HltoArBTKVVy8R2Osj7PKNE0LVQpVaxpWihwvIkxBTT+tXFWBI17Sy6rraZCfgJM0zTNS9O0GBp/K16w0/1MCWyk8U9ngJmAvXYYjwEOKqWa3EGsaVpHTdP8zv4/MI7GP9/t5sw3+qzJzXz97UBPTdNizmypTKNx3duNpmnjgflAilLqdDNj9FqfLVk/n9D42oPG1+LXzf2Cagta49UglgEHlFJ/b2ZMtzPj0DRtMI0/t+X2ynjm67bke/gJMENrNBSoUkrptdO92b/MHWF9nuf8119zHbgOGKdpWqCmaYE0rvumpyCe7zr3H02m8bdKPVACrDvvvt/SuP/yEDDhvM9/AYSd+f9YGks/G/gA8LLTfq83gYcu+lwY8MV5uXaf+dhP4+4He+8vfBvYS+M+90+A0Itznrk9ETh8Zl3rkTMbyAd2nfl4xZHWZ1PrB/gjjb+MALzPvPayz7wWY+28/kbQ+Cf2nvPW4UTgobOvUeCRM+ttN7ANGKbD97nJ7+FFOTXgxTPrei867W8HOtBY1v7nfU739UnjL5tiwHymN+cCnWl8byLrzL9BZ8aauPD9ojlnXqPZwOyWfL22PEJVCCGETuQIVSGEcEFS7kII4YKk3IUQwgVJuQshhAuSchdCCBck5S6EEA5A07RoTdNa7fgPKXchhHBBUu5CCOFgNE2L1TTtR03Tkq51GVLuQgjhQDRNi6PxXEOzlVLbr3U5VzxxmBBCCLsJpvH8Mj9XSu2/ngXJlrsQQjiOKhrP1TT8ehckW+5CCOE4GoCfAes0TTullHrvWhck5S6EEA5EKVWjadrtNF4QpUY1f1Gky5KzQgohhAuSfe5CCOGCpNyFEMIFSbkLIYQLknIXQggXJOUuhBAuSMpdCCFckJS7EEK4ICl3IYRwQf8fN1gEwGSN9vsAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 23, "metadata": {}, "output_type": "execute_result" } ], "source": [ "def phi(k):\n", " a=1\n", " x = a*k\n", " return (cos(x)+1)/((x-pi)*(x+pi))**2\n", "\n", "plot(phi(k), (k,-10,10))" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.0" } }, "nbformat": 4, "nbformat_minor": 1 }