Skip to content

Instantly share code, notes, and snippets.

@amcdawes
Created October 30, 2019 18:22
Show Gist options
  • Save amcdawes/437f8756fcb8ff8ca4c1b720600adf43 to your computer and use it in GitHub Desktop.
Save amcdawes/437f8756fcb8ff8ca4c1b720600adf43 to your computer and use it in GitHub Desktop.
Harmonic Oscillator states in the position basis
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"source": [
"from qutip import *\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt"
],
"outputs": [],
"execution_count": 1,
"metadata": {
"inputHidden": false,
"outputHidden": false
}
},
{
"cell_type": "code",
"source": [
"N = 50\n",
"psi2 = fock(N,0)\n",
"x = create(N) + destroy(N)\n",
"eigs = x.eigenstates()[1]\n",
"psi2_x = psi2.transform(eigs)\n",
"plt.plot(np.abs(psi2_x.full()))\n",
"plt.plot(np.real(psi2_x.full()),'.-')\n",
"plt.grid()"
],
"outputs": [
{
"output_type": "display_data",
"data": {
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
],
"image/png": [
"iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xl8lNW9+PHPdyYLSwgEAmELBCFsAoJBwIUKKu6KvVfv1dqW9tZLbfX+2tvlV21v9Xpb7+1mt19XW610Ua7Vqki1iJgoopV9DyFhDWsIGCCBLJP5/v54JjAJWSYzk8zyfN+vV16Z55nzPOeczGS+c855nnNEVTHGGOM+nlgXwBhjTGxYADDGGJeyAGCMMS5lAcAYY1zKAoAxxriUBQBjjHEpCwDGGONSFgCMMcalLAAYY4xLpcS6AO3Jzs7WvLy8sI6tqamhd+/e0S1QArB6u4vV211Cqfe6desqVXVgKOeL6wCQl5fH2rVrwzq2qKiIOXPmRLdACcDq7S5Wb3cJpd4isi/U81kXkDHGuJQFAGOMcSkLAMYY41IWAIwxxqUsABhjjEtZADDGGJeyAGBMKMpXw8onnN+h7DcmAcT1fQDGxIXy1fh/dzPi9+EXD8XZN1GdnkNG3VEmVL6ORxVJSYcFSyB3RqxLa0zILAAY04ETHzxLVmMDIuDxNzKxYikggOIRJ4021iF7V1oAMAklKl1AInKjiJSISJmIPNROustEpFFE7oxGvsZ0NfX7qdrxNgioeJHUnnjuW47nsSo89y1HvemogiqQNzvWxTWmUyIOACLiBX4O3ARMBO4RkYltpPsusCzSPI3pLutff5qLfLvZPurTyDXfaN7NkzsD+dRSygdcgQc/m0r3xrSsxnRWNFoAM4AyVd2tqvXAYmB+K+n+DXgRqIhCnsZ0uZrTVeSueZwy72jG3/sDmP3lC7t4cmcw+LMvsd8zjAErH6Gu9kxsCmtMGERVIzuB051zo6reF9j+BDBTVR8MSjMMeBa4BngKWKqqL7RxvoXAQoCcnJyCxYsXh1Wu6upqMjIywjo2kVm9o6du7TPcUP0Sr4z5b/oOv7jdtCf3rGX+vm/x174fo/e0f45qOdpjr7e7hFLvuXPnrlPV6SGdUFUj+gHuAn4btP0J4P+1SPNnYFbg8TPAnaGcu6CgQMNVWFgY9rGJzOodHftLN2vdI1m6+od3hXzM+u/epDWPDNQj5WVRLUt77PV2l1DqDazVED+/o9EFdADIDdoeDhxqkWY6sFhE9gJ3Ar8QkTuikLcxUad+Pyde+HfqSWPUPT8I+bicu36IBz8H/vcrXVg6Y6InGgFgDZAvIqNEJA24G1gSnEBVR6lqnqrmAS8An1fVl6OQtzFRt+mt/+WS2jVsHfs5sgePCPm4oaPGs2Hkpyk4/RbbVv21C0toTHREHABU1Qc8iHN1TzHwvKpuE5H7ReT+SM9vTHeqPVvDwFWPss+TS8FdbV7R3KZpdz/KIRlErxUP42uo74ISGhM9UbkPQFVfU9WxqjpaVR8P7PuVqv6qlbSf0jYGgI2JtdJFn2eYHsV36adITUvv9PE9emVw9PJHGeXfx6Ff/6NNEWHims0FZEzAqe1vMunwyygweuMPwv7wnjpuDH6E3Mp30EW3WRAwccsCgDEB5WuWIuJM8kBjPexdGdZ5ZP8qBHHO46sL+zzGdDULAMYEvFWVA4DiAW9a+FM75M2GlHQUUMSmiDBxywKAMUDF6Vr+djQTAJl8Z2Qze+bOQBYsobz3ZBpVONVndBRLakz0WAAwBnh9yxHSNXDVzpR/inxWz9wZnJn7GKnSSEnRc5EX0JguYAHAGGDp5kOM6Z/qbKT0iMo5x106l8MMJG2H3fJi4pMFAON6h6rOsmbvh1yV19vZEaUAIB4Pe4fcwMSz66mqPBKVcxoTTRYAjOu9tuUwADNHBCbZSo1OAAAYOOseUqWRnUXPRu2cxkSLBQDjeq9uOsSkYZnk9AzMjBulFgDA6MlXUC5D6VX6StTOaUy0WAAwrrb/+Bk2HTjJbVOGgq/W2RnFACAeDweG3ciE2k1UHimP2nmNiQYLAMbVXt3sTFx7y5Qh0BD9AAAw5Mp78Yqy6+0/RfW8xkTKAoBxtaWbD3PpiH4Mz+p1vgUQxTEAgLwJ09nrGUGfslejel5jImUBwLhWWUU1xYdPceuUoc6OLugCanI492YmNmzl6IFdUT+3MeGyAGBca+nmQ4gEun/ACQDiBW9q1PMaftXHANhj3UAmjlgAMK6kqry66RAz8vqTkxn4xt9Q2yXf/gFy8y+hzDua/nusG8jEDwsAxpV2HDnNrmM13HrJ0PM7fbVR7/8PVjnyFsb6dnJoz44uy8OYzrAAYFxp6eZDeD3CTZMGn9/p67oWAMCIj3wcgH0rrRvIxIeUWBfAmO7mdP8c5orRA8jOCFr1q4sDwNC8cZSkjCN37wuwsp8zTXSkk84ZE4GotABE5EYRKRGRMhG5YCFVEZkvIptFZKOIrBWRq6KRrzHh2HrwFPtPnHFu/grWhWMA57IYUsBw/yF0xbdh0e22WpiJqYgDgIh4gZ8DNwETgXtEZGKLZCuAS1R1KvAvwG8jzdeYcL1ZfBSPwLyJOc2f6OIxAICRg7IAEPwRrTpmTDREowUwAyhT1d2qWg8sBuYHJ1DValUNTLRCb0AxJkaKSiqYNiKLrN5pzZ/o4i4ggD5T78CPOP8Akaw6ZkwURCMADAOCJzk5ENjXjIh8VER2AH/FaQUY0+2Ona5j04GTzB038MInuyEAkDuDbQOup1GFM3f+0cYATExFYxBYWtl3wTd8VX0JeElEPgJ8C7iu1ZOJLAQWAuTk5FBUVBRWoaqrq8M+NpFZvdu36mADAH2q91NUdLDZc9OrKqntkcLWLv77nex7FZOPL+Ov735A3yOR/Qva6+0uUa+3qkb0A1wOLAvafhh4uINj9gDZHZ27oKBAw1VYWBj2sYnM6t2+B/60Ti/79nL1+/0XPvmTaap//nR0C9aK+rpaPfVIjn7wk3sjPpe93u4SSr2BtRri53c0uoDWAPkiMkpE0oC7gSXBCURkjIhI4PGlQBpwPAp5GxMyX6Ofd3YeY864gQTeji0SdEMXEJCalk5pn8vIO7EK9fu7PD9j2hJxAFBVH/AgsAwoBp5X1W0icr+I3B9I9o/AVhHZiHPF0D8HIpUx3WZDeRWnan3MHTeo9QTdFAAAfKPnMYgT7Nm+plvyM6Y1UbkRTFVfA15rse9XQY+/C3w3GnkZE67CHRWkeIQr87NbT9AN9wE0uWjmfNj0TY6uW8JFk2Z2S57GtGRTQRjXKCw5xvS8LDJ7tDHbZzfcB9Ake+hIyryj6XugqFvyM6Y1FgCMKxw5WUvx4VNtd/80NoA2dlsLAODYkKsZW7+dkyeOdVuexgSzAGBcoaikAoC549vp/4duDQBZl9xCivgp+/uSjhMb0wUsABhXKCypYFi/nuQPymg9QRetB9ye/GlzqCIDf8kb3ZanMcEsAJikV+/z825pZduXf0KXrQfcHm9KCmWZs7jo5Pv4Gxu7LV9jmlgAMElv7d4T1NQ3tt3/DzHpAgJgzDwGcJJdm1d1b77GYAHAuEBhSQVpXg9XjBnQdqIYBYDRl8/Hr0LlhqXdmq8xYAHAuEBhyTFmXtSfXmnt3PYSgzEAgKyBQyhNHUf/Q0Xdmq8xYAHAJLnyE2coq6huv/sHYjIG0OTEsDnkN+zkRMXBjhMbE0UWAExS6/DyzyaxGgMAsqfdikeU3XY5qOlmFgBMUissOUbegF6Myu7dfsIYBoDRk6+gkn5Qurzb8zbuZgHAJK3ahkbe21XJnI66fyBmYwAAHq+XPf0uZ/zpVfjf/r6tE2y6jQUAk7T+vvs4tQ3+jrt/IKZjAAB9howlg1qk8L9tsXjTbSwAmKT11o4KeqR6mDmqf8eJY9gFBJCXlYqqLRZvupcFAJOUVJUVxRXMzh9Ij1RvxwfEOAD0mHADKrZYvOleFgBMUio+fJqDVWeZNyEntANiOAYAQO4MtufegwDH5v3UFos33cICgElKbxYfRSSEyz+b+GpBPOBtY62AbpA190EAdu3aGbMyGHexAGCS0orio0zN7cfAPumhHdC0HGRbk8V1g2EXXcxeTy4Z++xyUNM9ohIARORGESkRkTIReaiV5+8Vkc2Bn/dE5JJo5GtMa46eqmXTgZNcF2r3D3TresDtOTx4LuNqt3Dyw8pYF8W4QMQBQES8OAu93wRMBO4RkYktku0BrlbVKcC3gCcjzdeYtqwodu7+7VQA6Mb1gNuTNfV2UqWR0vdeinVRjAtEowUwAyhT1d2qWg8sBuYHJ1DV91T1w8Dm34HhUcjXmFatKD5Kbv+ejM1pY/GX1nTjesDtyb90LifIhJK/xbooxgWiEQCGAeVB2wcC+9ryGeD1KORrzAXO1Pt4t6ySa8fntL34S2vipAvIm5LCrn5XMvbU+zTU18W6OCbJtTM/bsha+y/TVhOKzMUJAFe1eTKRhcBCgJycHIqKisIqVHV1ddjHJjK313v9UR91Pj+DGo5QVBT6YuuTKw6R2tDA+jj4253sPZnLql7nled+Sd+RU9tN6/bX222iXm9VjegHuBxYFrT9MPBwK+mmALuAsaGeu6CgQMNVWFgY9rGJzO31/r9/3qSTHv2b1vsaO3eCp29WferG6BcsDNWnPtTaRwbo+79Y2GFat7/ebhNKvYG1GuJnbDS6gNYA+SIySkTSgLuBZvPaisgI4C/AJ1TVLnI2XcLvV1bsOMrVYweS6u3kWztOxgAAevfpx46eUxleUYT6/bEujkliEQcAVfUBDwLLgGLgeVXdJiL3i8j9gWSPAAOAX4jIRhFZG2m+xrS08UAVldX1zJvYiat/msTJGECT2otuYLgeYX/JhlgXxSSxaIwBoKqvAa+12PeroMf3AfdFIy9j2rKi+ChejzBnbIh3/waLswCQd8U/wPZvc2jNy4ycUBDr4pgkZXcCm6Tx5vYKLsvLom+vMKZziJP7AJrkDB9NmXc0WeVvxrooJolZADBJ4dgZPyVHT3fu5q9gcTQG0OTYsGsZW1/M8aMHYl0Uk6QsAJiksLGiESC8/n+Iuy4ggIEFd+ARZZfdFWy6iAUAkxQ2HPMxZlAGIwd0sPZvW+IwAIyefDlHGUBqmd0VbLqGBQCT8E7VNlBywh9+90+jD/y+uAsA4vGwd8BsxlWvofZsTayLY5KQBQCT8Da/9waf9bzCHdlh9pXHeD3g9vScdBu9pI4Pn3/Q1gk2UWcBwCS28tXMeueTfCXlfxm37OPhfUjGeDnI9ozLHYgqDN7zF1ss3kSdBQCT0OrK3sarPjwCEu5i6nEcANIPrwVxJtxSWyzeRJkFAJPQ3m0Yj0Jki6nHej3g9uTNRj1pAKh4bLF4E1UWAExC+/WuATSKF78nDRYsCW8x9TgeAyB3BvrJJVTTk10p+bZYvIkqCwAmYe07XkPxvoOk0ojH74Phl4V3ojjuAgLw5l3OlqF3MapuB5VH9se6OCaJWAAwCevF9QcZ6DkJgOCH+jAvlYzzAAAw9OrPkCJ+yt58OtZFMUnEAoBJSH6/8uK6A1wzLGjtobpT4Z0snscAAkaOm0pJynhydv/Fpog2UWMBwCSkv+85zsGqs9yQF/QWrg0zAMTzGECQqnF3Mcq/j7LNq2JdFJMkLACYhPTCugP0SU9halbD+Z3htgASoAsIYPx1n6JOUznx7u9iXRSTJCwAmIRTXefj9S1HuPWSIaTWVp5/ItIWQJwHgL5Z2WzNnM24ymXU1Z6JdXFMErAAYBLOa1sOc7ahkTsLhkP10fNP1J0M74QNZ53fcR4AAFILPk4/qtlW9Hysi2KSgAUAk3BeWHeAUdm9uXREFlQfg14DnCfCbgHUOb/jfAwA4OKr5lNBfzybnot1UUwSiEoAEJEbRaRERMpE5KFWnh8vIu+LSJ2IfCUaeRp32ne8htV7TnBnwXBEBGoqYMAY58naMFsAvsRpAXhTUtg19FYmnVnNMbsnwEQo4gAgIl7g58BNwETgHhGZ2CLZCeD/AD+IND/jbi+uP4gIfHTaMGdH9THIykPxRDAIXAeIM5VEAmi6J2DXmzYYbCITjRbADKBMVXeraj2wGJgfnEBVK1R1DdDQ2gmMCUXTtf9XjclmaL+eoOq0ADIG4UvpGX4XUMNZ59u/SHQL3EVGjpvKjpTx5Ox+0e4JMBGJRgAYBpQHbR8I7DMmqpqu/b+zYLizo+6UcwVPRg6+lN6RtQASoP8/2Mmxzj0BOduftCmiTdhSonCO1r42aSv7QjuZyEJgIUBOTg5FRUVhnae6ujrsYxNZMtf7N5vr6OGFHsd3UlRUSs8zB5kJFO8/zlDpQc3B3WwNo+7jyvfQv1F4P4H+bj3woArjK1+n8Xcr2HTJtzjVd3ysi9Vtkvl93p5o1zsaAeAAkBu0PRw4FO7JVPVJ4EmA6dOn65w5c8I6T1FREeEem8iStd4Vp2pZt6KQOy7N5YZrpzg7970Hq2HCZR+h6vAbZGekhVf343+Eur6J9XdbuQ7d7jThVRu5tH8NzJ4T61J1m2R9n3ck2vWORhfQGiBfREaJSBpwN7AkCuc15pyfvlWKr1G5/+rR53dWVzi/ew/Cl9Ir8jGARJI3GzzO9zcfKbZOgAlLxAFAVX3Ag8AyoBh4XlW3icj9InI/gIgMFpEDwJeA/xCRAyKSGWnexh32VtaweHU598wYQV527/NPNAWAc2MA4V4GmnhjAOTOQK75DwAebfgURzKnxLhAJhFF5T4AVX1NVceq6mhVfTyw71eq+qvA4yOqOlxVM1W1X+BxmF/XjNs8sXwnqV4P/3bNmOZP1FSAeKBXfxq9EbQAfLWJ1wIAGHcLAPV4+cmK0hgXxiQiuxPYxLWtB0/y6qZD/MtVeQzKbPEhXV0BvbLB43W6gOpOO5eGdlaiBoD+o/CLl9uG1fD82nJ2H6uOdYlMgrEAYOLa95eV0LdnKgs/MvrCJ2uOQcYgAKcLSBvDWxSmIUEDgDeV2h45zMo8TnqKhyeW74x1iUyCsQBg4tb7u47z9s5jPDB3NH17pl6YoLoCeg8EcFoAEN69AL7axBsDCDjTazjpJ3dz31Wj+Ovmw2w5EOY4iHElCwAmLqkq31u2g8GZPfjk5XmtJ6qugIwcAGcMAMIbB0jULiDgTK9hcHwX9101kqxeqXxv2Y5YF8kkEAsAJi69sf0oG/ZX8cXr8umR6r0wwblpIJpaAIGrg8JtASRsABgOjXVk1h7mgbljWFlayXtllR0faAwWAEwcavQr319WwkUDe5+f9qGlutPOB3fvpjGAphZAGF0giToGQKAFAFBZysdnjWRo3x5892870HAGw43rWAAwcecv6w9QVlHNV68fR4q3jbdozTHnd/AgMIQXABJ6DCAQAI6X0iPVyxfnjWXTgZMs23YktgUzCcECgIkrZ+p9/PjNUqYM78uNkwa3nbBpJbBAADg3BtDZLiB/I/gbErYF4EvNhJ79odK5Augfpg1jzKAMvreshDpfY4xLZ+KdBQATNxr9yhcWb+TwybN8/eYJzoIvbQmaBgKCu4A6GQASZD3gdmWPhUrnRrAUr4dv3DyB3cdq+NoLm60ryLTLAoCJG4//tZjl24/y6G0XM+uiAe0nbtEF1Ojt6dwV3NkWQEMyBID8cwEAYO74QXz1hnG8vPEQP3rT7hA2bbMAYOLCovf28vSqPXz6yjwWXJHX8QHVTdNABAKFCKRnht8CSNAxAMAJADUVcLbq3K7PzxnNP00fzk9XlPLiugMxLJyJZxYATMy9teMoj726jesm5PAft7RcTbQNNRXOh78n6BLRHpmdbwEkSxcQwPGyc7tEhMc/Opkrxwzgob9s5v1dx2NUOBPPLACYmNp68CQPPruBi4f25af3TMXrCXFZxqCbwM5J7+vOMYAB+c7vyuZTQaR6Pfzi3gLyBvTms39YS1mFzRVkmrMAYGKjfDWnln+XH/3uj/TrmcpTC6bTK60T6xMFTQNxTjgtgGQYA8gaCZ7UZuMATfr2TOXpT11GWoqHTz+zmpM734WVT9gykgaIzopgxnRO+Wp00W309tXxM03l6K3PXzjTZ0dqKmBAiwni0jPhZCf7u5NhDMCbCv1HXdACaJLbvxe/XXAZ33nyGXo9+21UGhFvOixYArkzurmwJp5YC8B0K1Wl+O+vob5avCg9xEfe6fWdPQlUH2ujBdDJG8GSoQsInHGAoDGAlqbm9uP7o7eQSgOiftRXB3tXdmMBTTyyAGC6zYEPz/CZRWt5akMNTT394vF0fjnDutPgO9vKGEAEVwElegAYMAaO74JGX5tJcnsrCigg+ClZ9QrHDu3trhKaOGRdQKbLNfqVZ97byxNvlODBT2H2KqjPgpQ06NG3890QLe4BOKdH5vlFYdq7iSxYsgSA7LHOHc1V+y7sGmtyZBMybDqN+TewvWQHYw8tof7Xs/jg4i9x2czZePa/5wRj6xZyDQsApsvUNjSyqqySH79ZypaDJ5k7biA/HL2erLe2wT/8Fqr2wlvfDlzRM6jD851z7i7gFl1A6ZnnF4VJzwjtXA1JMAYAzr0A4AwEtxYAju9yuohu+h7emZ9l8hwoL9vCyT8/yMztj+PfLihASjqy4FULAi4RlS4gEblRREpEpExEHmrleRGRnwae3ywil0YjXxNHylfDyieo2vkuz68tZ+Hv1zLtv5bzmUVrOXyylp99bBpP35VH1nv/43zLnHwn5N/gHFu6vHN51TQtBt9KCwA6dyVQsrQABgTWSz7exp2/TX/j/HnnduWOmczFXyukfMj1CIqg0FDLwT/8K2te+QVVlUfOva521VByirgFICJe4OfAPOAAsEZElqjq9qBkNwH5gZ+ZwC8Dv7tG+WpG7HsByns1/yZTvtoZ+GrZzO3s/mieK8p5R6XeIRzjGzaTChnAhwd2InveZtyeZxD1k66pLK7/Okcyp3DX9OFcNyGHmRf1Jz3FCy9/3vl2fssTThfN4MnQZwiUvgHT7iVkLeYBOic9EABqT0Lm0NDOlSwBoFd/Z33kNq4EovQN536B/hc12y0eD7k3fxVd9A7qq8cv0Kv+OJdteJjG9Q/jF0EAv3jZO/0b9Lz4FgYNG0XKkY0x/Z/p1Ps8ynnHNI8oi0YX0AygTFV3A4jIYmA+EBwA5gO/V2dmqr+LSD8RGaKqh6OQf3Plq/H/7mby/D78v1tM6fjPU9NnFL1P7yF/xy8Qvw/1pIS9H4jauboi76Z67xz/ANWB/WN3/ByP34ffk0LJuPup6TWS3tW7GbfzN3jUh1+87Mj7BNU9BqONPjJq9jPh0It4tBE/Hkr6zOSsppFZf5TR9TsQFK/CUIGWH7M9pJ5fzjrBoNuvaT6Z2773YOOf4KovwcBxzj4R5xvptpehscG5nDEUTdNA9M5ukXlf53dnBoKTJQBAs0nhmqmvgb3vwmX3tX5c7gyn22fvSrx5s+k7tICdm94lvegxRpxahwBe9TF6zWOw5jF86sEvfgRQPOzsNZVT6UPoRS0TqorwqB+/eCke/k/UZIwETwp9ag8zbveiwPsthZ35/0pNnzx6n9nPuB2/RtR5P+8c9zlq+lzkvG9LOn6fx8P/a1flISk9uvxSXYl0tkARuRO4UVXvC2x/Apipqg8GpVkKfEdV3w1srwC+pqprWznfQmAhQE5OTsHixYs7VZ4R+15g1O4/hDwGaNqnCtX05Lhk0YuzDORD5xshsKP3TEoH3Uy2t4ZZu3+Mx98AKPWpmRRP/CpVWVMAEL+PgnVfIsV3ltUzfobfm37u/NnH3mfStu+wYerjnOw3KaQyjS35BdmVf+e9K39/bl91dTVDG8u5dMNDbJ78CCcGFIR0rrw9fyJv3/MUXf1y6APHcaS6upqMDGe8Y2zJz8iuXN3s7wIwoHI1k7c+zsZL/ouqrEtCPnfmyR1csumbzoeUeFk1+JNU+VIZf/Id8uu3BwIAVNGHBk2hL6dJw5eIf8a45MfD3lH3sn/knef2Bb/ebZk7d+46VZ0eSh7RaAG09nK3jCqhpHF2qj4JPAkwffp0nTNnTudKU94L//7n0cA3yiNXf5+67ImkV25n8NtfRfwNqCf8/UDUztUVeTd9kz569fcC+7eR8/bXkMYG1JtKxTVP0DBwMmmV2xm44ouIvwE8qVTd+ht02HS8KWmkVmyi5wv3QmMDkpJGnwVL6JM7w2maLrodGuvxeNOYePe3mdj07aR8ntNs7T2I9FU/ZuqmR+CKB+Gab8LqJ6FmH9z9HB8Zf0Pz16uuAIqfYFqvozDnQUJy+EloHE7we6OoqIhLJ86BDTBl7EiYHOL7pn4FHOzBnLlzQ0sfZ4qKis7/HdK2wBvLmTNjitMl1OTVlyEtg6m3fc658ipkc+DSS891R3zk3Gt9/n0g3jSyFvzF+ZZavhpddDvaWA/eVGo++nvqsi/G3+jDc2gdWa9/7tz79ti1P6Iu+2JSj20l560vndt/dM73qcu+2HnfFnX8Po+H/9euysPjTeeiaz7JRUEtgGavdzSoakQ/wOXAsqDth4GHW6T5NXBP0HYJMKSjcxcUFGhY9n+gu555QHX/Bxfs13d+EPn+aJ4rynlHpd7hHtOkrlr11S+qPpqp+uOpqo8NUP3t9W2nf+Y21Z/NbPv5ln5zreqi25vtKiwsVK064OS55qnQz/XXr6j+z4jQ08eZwsLC8xslf3PqH/za+P2qT0xUfe5j0c04xv8znXqfRznvWObR7PVuA7BWQ/38DjVhmydwWhG7gVFAGrAJuLhFmluA13FaArOA1aGcO+wAEOIfKhnFVb3f/YnzgfRopuq3BrYdNN77mZPmw32hnfdHk1Rf/NdmuwoLC1VrTzvnWfmj0Mv4yoOq3x8bevo40+z1rixz6r/+D+f3Hdnq7Fv7TLeXrSvF1fu8G0U7AER8Gaiq+oAHgWVAMfC8qm4TkftF5P5AstcCQaIM+A3w+UjzNQnA3+AM1oJzh2pbUw/kX+/8Ln2j43O2NQ0EQFpvEG/nLgNtSNz1gC/QbyR405oPBDf9TZv+xsYEicqNYKr6Gs6HfPC+XwU9VuCBaORlEkjebPCmQ2O988GlfchmAAAOdElEQVTU1pQPA8ZA1ijY+UbbV6o0qa8OTAPRyo1jIpDep/NXASXDFUAA3hTnMs/gALDzDedy28whsSuXiVt2J7DpOrkznMvYOrqmWcT5hrr+99BwFlJ7tn3Otu4BaNLZKaGTKQCAc0fwsRLn8dkPofwDuOrfY1smE7dsMjjTtXJnwOwvd3wt89jrnW/2e99tP925eYBa6QKCzi8Kk2wBYEA+nNjtXCWz6y1nagzr/jFtsABg4sPIqyC1V8fjANVHnd8tZwJt0tkWQDKNAUBgUjgffLjPmf6hZxYMD+mScONCFgBMfEjtAaOuhp3LnIHetnTUBdTZKaGTrQXQNCncsR1OABhzXfN1k40JYgHAxI/8ec50xq1NZ9Ck5hggzoLwrensojDJFgCaJoXb8jycqTw/4Z4xrbAAYOLHuctBl7WdprrC+fD3tnH9gttbAD37Oa2j4ledS3DHXBvrEpk4ZgHAxI9+uTBoYvvjANUVbff/w/kxgFDnuEq2MQBwxgHUD8Mvaz4lhDEtWAAw8SX/eti7Cgr/u/U56Gsq2r4CCAKLwvid+wVC4atLrhYAnF8XYdCE2JbDxD0LACa+9M11Ll18+3vOhGMtg0B1RdsDwHD+wy/UbiDf2eQKAOWrzy/+svE5W8jFtMsCgIkvtVWBB+rcQdxy+oiaY+0vH5neiVXB/H4nj2QKAHtXOi0gcC4HbWv6DWOwAGDizaiPONNGgDOIGTx9RF01NJxpfR6gJp1ZFMaXJOsBB8ub7fz9xNv+9BvGYFNBmHiTOwMWLIUX74P605ATtEhMRzeBQedaAMm0GliTUKffMAZrAZh4NGIm/ONvnLlsPvjl+f0dTQMBQS2AEO4FSMYAAKFPv2FczwKAiU8jZsHYm+Ddn8CZE86+ju4ChvODwG5tARjTCRYATPy69pvOB/mqHzvbNYEAEMogcChjAA1JOAZgTCdYADDxK+dimPLP8MGv4eTBQAtAoFd228d0ZlEYawEYl7MAYOLb3IfB3whvf7fjaSCgc4vCWAAwLmdXAZn4lpUHl30GVv8GBo5vv/unSahTQlsAMC4XUQtARPqLyHIRKQ38zmoj3dMiUiEiWyPJz7jU7K84H9IV25wVwzq6uzXURWFsDMC4XKRdQA8BK1Q1H1gR2G7NM8CNEeZl3CpjIFz8Uefxh3tanyIimLUAjAlJpAFgPrAo8HgRcEdriVT1HeBEhHkZN+s77Pzj1qaICBbqlNAWAIzLRRoAclT1MEDgdwgdtMaEYcx1zgd1KFMc9Mh0941gxoRItIN500XkTWBwK099A1ikqv2C0n6oqm2NA+QBS1V1UmvPB6VbCCwEyMnJKVi8eHG75WtLdXU1GRkZYR2byJK53pknd9CvaitV/SZxqu/4Zs8F1zt/568ZVPE2q656tt3zDS9fwphdT/HulX/Cl5qYf7Nkfr3bY/Vu29y5c9epamgLQatq2D9ACTAk8HgIUNJO2jxga2fOX1BQoOEqLCwM+9hEZvVW1TcfU/3Pfqp+f/sHvfOE6qOZqvVnurRsXcleb3cJpd7AWg3xMzbSLqAlwILA4wXAKxGez5jIhboojHUBGZeLNAB8B5gnIqXAvMA2IjJURF5rSiQizwHvA+NE5ICIfCbCfI1pW6iLwvhqwZvu3DxmjAtFdCOYqh4HLlh1WlUPATcHbd8TST7GdEqzKaGHtZ0uGdcDNqYTbCoIk3xCXRTGV2vdP8bVLACY5NMUADq6GcwCgHE5CwAm+ZybErqDewEsABiXswBgkk+oi8LYGIBxOQsAJvmEuiiMtQCMy1kAMMkn1EVhLAAYl7MAYJJPqIvCWAAwLmcBwCSnUKaEtjEA43IWAExyCmVRGGsBGJezAGCSUyhTQlsAMC5nAcAkp/QQuoAsABiXswBgklOPEFYFszEA43IWAExySs+Euna6gFShsc5aAMbVLACY5NQjE+pOOx/0rbG1AIyxAGCSVEeLwlgAMMYCgElSHS0K0xAIADYGYFzMAoBJTh1NCW0tAGMsAJgkld7BojAWAIyJLACISH8RWS4ipYHfWa2kyRWRQhEpFpFtIvKFSPI0JiQdTQltAcCYiFsADwErVDUfWBHYbskHfFlVJwCzgAdEZGKE+RrTvo4WhbExAGMiDgDzgUWBx4uAO1omUNXDqro+8Pg0UEy7K3UbEwXWAjCmQ5EGgBxVPQzOBz0wqL3EIpIHTAM+iDBfY9rX0aIwFgCMIaWjBCLyJjC4lae+0ZmMRCQDeBH4oqq2eY++iCwEFgLk5ORQVFTUmWzOqa6uDvvYRGb1DlDlajzs37mFPb6iC9IPrFjPxcCaDVuoKe1gyog4Zq+3u0S93qoa9g9QAgwJPB4ClLSRLhVYBnypM+cvKCjQcBUWFoZ9bCKzegf5nxGqS7/c+gEbnlV9NFP1+K4uLVdXs9fbXUKpN7BWQ/yMjbQLaAmwIPB4AfBKywQiIsBTQLGq/jDC/IwJXXtTQlsXkDERB4DvAPNEpBSYF9hGRIaKyGuBNFcCnwCuEZGNgZ+bI8zXmI6l97VBYGPa0eEYQHtU9ThwbSv7DwE3Bx6/C0gk+RgTlvamhLYAYIzdCWySWHuLwjRYADDGAoBJXo31UFUO5asvfM5XC9408Ni/gHEve/eb5FS+GnYXOYvCLLr9wiDgq4WUnjEpmjHxwgKASU57VzrrAYDTEti7svnzvlpISe/+chkTRywAmOSUN9vp4gHweJ3tYLYesDEWAEySyp0Bn3wFvOkweq6zHcxXawPAxvUsAJjkNfJyGH8zHNp44drAFgCMsQBgklz+9VB9FI5sbr7fAoAxFgBMkhtznfO79I3m+20MwBgLACbJZQyCodOgdHnz/dYCMMYCgHGB/OvhwBo4c+L8PgsAxlgAMC6Qf71zT8Cut87vswBgjAUA4wJDp0GvAc3HAWwMwBgLAMYFPF5nMLjsTfA3OvusBWCMBQDjEvnXw5njcGiDs20BwBgLAMYlRl8D4nG6gVQtABiDBQDjFr36w/DLnADgq3P22RiAcTkLAMY98uc5XUAny51tawEYl7MAYNwj/3rn946lzm8LAMblIgoAItJfRJaLSGngd1YraXqIyGoR2SQi20TksUjyNCZsg6dAxmAoftXZtgBgXC7SFsBDwApVzQdWBLZbqgOuUdVLgKnAjSIyK8J8jek8Eci/Dg6uc7ZtDMC4XKQBYD6wKPB4EXBHywTqqA5spgZ+tGU6Y7pFUzcQWAvAuJ5oy3nSO3OwSJWq9gva/lBVW+sG8gLrgDHAz1X1a+2ccyGwECAnJ6dg8eLFYZWturqajIyMsI5NZFbv9nl9NVy56hN4tJFNU/6TD/tP64bSdR17vd0llHrPnTt3napOD+V8KR0lEJE3gcGtPPWNUDIAUNVGYKqI9ANeEpFJqrq1jbRPAk8CTJ8+XefMmRNqNs0UFRUR7rGJzOodgrJJcGQTlwzvAzNCPCZO2evtLtGud4cBQFWva+s5ETkqIkNU9bCIDAEqOjhXlYgUATcCrQYAY7pU+Wqo2OY8XvZ1GDLlwuUijXGJSMcAlgALAo8XAK+0TCAiAwPf/BGRnsB1wI4I8zUmPHtXOjODAvh9zrYxLhVpAPgOME9ESoF5gW1EZKiIvBZIMwQoFJHNwBpguaoujTBfY8KTN9tZKF684E1zto1xqQ67gNqjqseBa1vZfwi4OfB4M5DYI20meeTOgAVLnG/+ebOt+8e4WkQBwJiElDvDPviNwaaCMMYY17IAYIwxLmUBwBhjXMoCgDHGuJQFAGOMcSkLAMYY41IRTQbX1UTkGLAvzMOzgcooFidRWL3dxertLqHUe6SqDgzlZHEdACIhImtDnREvmVi93cXq7S7Rrrd1ARljjEtZADDGGJdK5gDwZKwLECNWb3exertLVOudtGMAxhhj2pfMLQBjjDHtSLoAICI3ikiJiJSJyEOxLk9XEpGnRaRCRLYG7esvIstFpDTw+4I1mhOZiOSKSKGIFIvINhH5QmB/ste7h4isFpFNgXo/Ftif1PVuIiJeEdkgIksD226p914R2SIiG0VkbWBf1OqeVAEgsPj8z4GbgInAPSIyMbal6lLP4CyvGewhYIWq5gMrAtvJxAd8WVUnALOABwKvcbLXuw64RlUvAaYCN4rILJK/3k2+ABQHbbul3gBzVXVq0OWfUat7UgUAYAZQpqq7VbUeWAzMj3GZuoyqvgOcaLF7PrAo8HgRcEe3FqqLqephVV0feHwa50NhGMlfb1XV6sBmauBHSfJ6A4jIcOAW4LdBu5O+3u2IWt2TLQAMA8qDtg8E9rlJjqoeBufDEhgU4/J0GRHJw1lt7gNcUO9AN8hGoAJnaVVX1Bv4MfB/AX/QPjfUG5wg/4aIrBORhYF9Uat7sq0IJq3ss8uckpCIZAAvAl9U1VMirb30yUVVG4GpItIPeElEJsW6TF1NRG4FKlR1nYjMiXV5YuBKVT0kIoOA5SKyI5onT7YWwAEgN2h7OHAoRmWJlaMiMgQg8LsixuWJOhFJxfnw/5Oq/iWwO+nr3URVq4AinPGfZK/3lcDtIrIXp0v3GhH5I8lfb+Dc+uqoagXwEk43d9TqnmwBYA2QLyKjRCQNuBtYEuMydbclwILA4wXAKzEsS9SJ81X/KaBYVX8Y9FSy13tg4Js/ItITuA7YQZLXW1UfVtXhqpqH8//8lqp+nCSvN4CI9BaRPk2PgeuBrUSx7kl3I5iI3IzTZ+gFnlbVx2NcpC4jIs8Bc3BmCDwKPAq8DDwPjAD2A3epasuB4oQlIlcBK4EtnO8T/jrOOEAy13sKzoCfF+eL2/Oq+l8iMoAkrnewQBfQV1T1VjfUW0QuwvnWD053/bOq+ng06550AcAYY0xokq0LyBhjTIgsABhjjEtZADDGGJeyAGCMMS5lAcAYY1zKAoAxxriUBQBjjHEpCwDGGONS/x/CUcak9Ho3XgAAAABJRU5ErkJggg==\n"
]
},
"metadata": {
"needs_background": "light"
}
}
],
"execution_count": 2,
"metadata": {
"inputHidden": false,
"outputHidden": false
}
},
{
"cell_type": "code",
"source": [
"eigs2 = []\n",
"for eig in eigs:\n",
" if eig[0] < 0:\n",
" eig2 = eig*-1\n",
" else:\n",
" eig2 = eig\n",
" \n",
" eigs2.append(eig2)"
],
"outputs": [],
"execution_count": 3,
"metadata": {
"inputHidden": false,
"outputHidden": false
}
},
{
"cell_type": "code",
"source": [
"psi2_x = psi2.transform(eigs2)\n",
"plt.plot(np.abs(psi2_x.full()))\n",
"plt.plot(np.real(psi2_x.full()),'.-')\n",
"plt.grid()"
],
"outputs": [
{
"output_type": "display_data",
"data": {
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
],
"image/png": [
"iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xt4XPV95/H3d2Z0s+WrJMvWzZJs2fgCGDDGEJw4bJNwSSBJk23YNEu2aVl2kybd3WyXbrq73d6S9Elv+4Q2pUkaWtqQtARwCZQkBIEbLsLGFzCyLdmWdbMly7Zutq4zv/1jRkQIyxpJM3Nmznxez6NHM+ecOef705n5+PjMOb+fOecQERF/CXhdgIiIJJ7CXUTEhxTuIiI+pHAXEfEhhbuIiA8p3EVEfEjhLiLiQwp3EREfUriLiPhQyKsNFxcXu+rq6jm99sKFCyxcuDCxBWWIbG272p1d1O7p7d27t8c5VzLTujwL9+rqavbs2TOn19bX17Nz587EFpQhsrXtand2UbunZ2Yn41mXTsuIiPiQwl1ExIcU7iIiPqRwFxHxIYW7iIgPKdxFRHxI4S7S1gC7/zj6ezbzRNKYZ9e5i6SFtgbcQx+C8REiFqCx+DYG80oBKBzpYkPP0wRwWDAP7tkFlds8LlgkPgp3yW4tu2F8BMMRdGE2dj8JWGymIzDxMDwaXVbhLhlC4S7ZrXoHbw0RHyogMPnovK2ByLc+gBGBYA5WvcOjIkVmT+fcJaudyysn4qB18bXY1NMuldto2vR5DDi+5b/rqF0yisJdslrT898lZI7w+79yyfBefcd/56LL4+zx/R5UJzJ3cYW7md1qZkfMrNnM7r/MctebWdjMPpa4EkWSp7BpF62Bcmo2Xn/J+fkLCnlz8c2sO/ccY6MjKa5OZO5mDHczCwIPALcBG4G7zWzjNMt9FXgm0UWKJEPP6VauGDlIR/ltWGD6j0Lwyo+ylEEaf/bPKaxOZH7iOXLfBjQ7544750aBR4C7LrHcrwOPAt0JrE8kaY49/w8EzVF207+77HIb3/1RBlwBwwceTVFlIvMXz9Uy5UDbpOftwA2TFzCzcuAjwC3Apf9/G13uXuBegNLSUurr62dZbtTg4OCcX5vpsrXtyWh30ZHHOEYlbV0DnOi6/LoDeddz9fl6nv3JjwiGchNax+Vof2eXRLY7nnC3S0xzU57/GfA/nHNhs0stHnuRcw8CDwJs3brVzbUz/mztyB+yt+2Jbnd3xwmKnzvCK9X/Ma71HoicZskLL1Bk59iy8xMJq2Mm2t/ZJZHtjue0TDtQOel5BdA5ZZmtwCNm1gJ8DPgLM/twQioUSYLjzz9MwBwVN1/+lMyEDe/6MH0sZPzAPyW5MpHEiOfI/VWgzsxqgA7gE8DbPhHOuZqJx2b2HeBJ59zjCaxTJKGWHv8hx4K1rKm7Oq7lc/Py2b/0PWw8/xzDQxfIL8i+8T0ls8x45O6cGwc+R/QqmEbg+865Q2Z2n5ndl+wCRRLtdGsTV4w30l11+6xeV3DNxym0IRp3P5akykQSJ67uB5xzTwFPTZn2jWmW/fT8yxJJnpYXHmYlULXjl2f1ug03fZDzzy0i8vqj8P7ZvVYk1XSHqmSdopYf0hSqo7x2w6xeF8rJ5ejy97Kh/2cMXRhIUnUiiaFwl6zScbyRuvEmzlbfMafXL7z24yywERpf0Berkt4U7pJVWnc/DED1u+d2WmXD9tvpYSkc+kEiyxJJOIW7ZJUVrU9xOLSBlVV1c3p9MBTiWPEtbBh4mQsDvQmuTiRxFO6SNbpe/HvWhI8zvuqaea1n0XX/lgIb5fwj92n4PUlbCnfJDm0NFP/48wBs6nx0XqF8RdlSHFDe8TQ8dKcCXtKSwl2yQ8tuAm4cAIuMR4fMm6NA64vR9cDPh98TSTMKd8kOseH0HEAwF+YzZF71DiIWvUXEBXPmty6RJFG4S1boK6wlALQu2w5Th9ObrcptnNjymwCcuOZ+Db8naUnhLlmh/c1XAOi7+t6EhPGy7Z8EoLu3f97rEkkGhbtkhYGWPQCUb7hhhiXjU1RaQRdFhLoOJmR9IommcJesEOo6SBdFFJVWJGydnQvWUzJ4JGHrE0kkhbtkhZLBI3QuWJ/QdQ4XX0lluJ2Lg30JXa9IIijcxfcuDPRSGW5nuPjKhK63YPW1BMzReuiVhK5XJBEU7uJ7bW82EDBHweprE7re8g3bAeg9vieh6xVJBIW7+N5E+E6EcaIUr6yih6UE9aWqpCGFu/hesOsgPSyleGVVQtdrgQAd+XUU9TcmdL0iiaBwF98r6m+kI78OCyT+7X6xaDNV4VaGhy4kfN0i86FwF18bHrpAVbiVi0Wbk7L+/KprCFmE1kadd5f0onAXX2tt3EPIIuRXza+b3+mUro+exz9/7NWkrF9krhTu4msToTsRwom2qqqOPhbCqQNJWb/IXCncxd9OHaCPhaya48hLM7FAgLa8OpbrS1VJMwp38bXl/Y205SXny9QJg8s2sXrsBGOjI0nbhshsKdzFt8ZGR1g9doLBZZuSup1Q5TXk2jitR/YldTsis6FwF99qPbKPXBsnVJmcL1MnlK6LdiF8tlnD7Un6ULiLb02E7UT4Jkt57WYuuHxcx/6kbkdkNhTu4luuYz8XXD7ltcm5xn1CIBjkZO5alvTpS1VJHwp38a0lfY2czF1LIBhM+rb6l22kavQY4fHxpG9LJB4Kd/Gl8Pg4VaPH6F+2MSXbC5ZtYYGN0N6sTsQkPSjcxZfamw6wwEYIlm1JyfaKY+f1zxxV3+6SHhTu4ktnmqJfphYn+cvUCZV1VzPschjXl6qSJhTu4kvjHfsZdjlU1l2dku2FcnI5mVPLovOHUrI9kZko3MWXFp0/xMmcWkI5uSnbZu+SDVSNNBMJh1O2TZHpKNzFdyLhMFUjzfQu2ZDS7VrZFhbZEJ0ndEmkeE/hLr7TeaKRRTaElSX3ztSplq+Nnt/v0peqkgYU7uI7g3u/B0DZssKUbrfqiusYcwGKG/8O2tQVgXhL4S7+0tZAXeMDOAflL34ppSGb232QoDmqBvbBQ3cq4MVTCnfxl5bdBAhjBhYeg5bdKd224TCA8Ghqty0yRVzhbma3mtkRM2s2s/svMf8uMztoZvvNbI+Z3Zz4UkXiUL0DAAcQzH3reaq2HbFQdPvBnNRuW2SKGcPdzILAA8BtwEbgbjObek/3s8DVzrktwK8A30x0oSLxGFi8FgNal22He3ZBZWpuYgKgchsntvwmACeu/a3UbltkiniO3LcBzc654865UeAR4K7JCzjnBp1zLvZ0IbEDJ5FU62yK3iF6bsOnPAnXwmt/EYAzg6Mp37bIZPGEeznQNul5e2za25jZR8zsMPBDokfvIinX1/YGACW1qelTZqrS8loGXQF0H/Zk+yITQnEsY5eY9o4jc+fcY8BjZvZu4PeAX3jHiszuBe4FKC0tpb6+flbFThgcHJzzazNdtrY93naPNb3CsMvhaOtpmtt7kl/YJawIlJN37nBC9pP2d3ZJZLvjCfd2oHLS8wqgc7qFnXMvmNkaMyt2zvVMmfcg8CDA1q1b3c6dO2dfMVBfX89cX5vpsrXt8bb7wMu/T0eokltuecexRco07F9Hbe+LbEnAftL+zi6JbHc8p2VeBerMrMbMcoFPALsmL2Bma83MYo+vBXKBswmpUGQWVgy3cH5hrac1RIrXU0wvvT2nPa1DstuM4e6cGwc+BzwDNALfd84dMrP7zOy+2GK/CLxhZvuJXlnzS5O+YBVJiQsDvaziDGPL6zyto6B8EwCdzer+V7wTz2kZnHNPAU9NmfaNSY+/Cnw1saWJzE5H0wHWAfllyR0zdSala7bA8zDQ9gZsv9XTWiR76Q5V8Y3ek9Eh7oprrvS0jtKKNVxw+dCt3iHFOwp38Y1IVyOjLsSq6tR29TuVBQJ05FSxsL/Z0zokuyncxTcK+pppD1akdICO6fQVrqF05KTXZUgWU7iLb5QMn/D8SpkJ4aJ1lHCevnNnvC5FspTCXXzh4mAfZa6b0eXrvC4FgILy6Je6p5r2eVyJZCuFu/hCR3P0y9T8sql92nmjpDY6MHd/rDsEkVRTuIsv9LZEw3159dUeVxK1snItF10eEV0xIx5RuIsvjHcdZtQFKa9NjyP3QDAYvWKmT1fMiDcU7uILBX1NdKTJlTITehfWUjrS4nUZkqUU7uILJUPHObegxusy3ma8aD0rOEffeW96p5TspnCXjDd0YYBVkW5Gl6/3upS3mehj5pT6mBEPKNwl43U0HyRgjrxV3t6ZOlVJTeyKmVZdMSOpp3CXjDfRp0xRTXpcKTNh1ep1DLlcXTEjnlC4S8YbO93ImAtSVrvJ61LeJhAM0hGqZIGumBEPKNwl4+X3NtERLCMnN8/rUt6hd2EtpcMnvC5DspDCXTJeydAJzi1Ijz5lphorWk8pZxnoO+d1KZJlFO6S0YYvDrIqcpqRZenRp8xUBWWxUZmadMWMpJbCXTJaR/NBgubITbMrZSYU11wFQH/r6x5XItlG4S4Z7fzJaGim25UyE1ZVX8GwyyHcpStmJLUU7pLRfn6ljLfjpk4nGArREaqkQFfMSIop3CWj5Z8/SmewjNy8fK9Lmdb5hbWUDrd4XYZkGYW7ZLTioROcLUivPmWmGl++npWcYbD/vNelSBZRuEvGGhm+SFnkFKPL6rwu5bLyYgOIdMYGFBFJBYW7ZKyO5tcJmiNnVXr04T6d4tioTBPdJIikgsJdMta5kxOjL13pcSWXt2r1ekZcDhFdMSMppHCXjDV26k3GXYCyNekd7qGcXDqC5RT0NXldimQRhbtkrOU9exkMFJJ3Jv271B3JL6Zm6BC0NXhdimQJhbtkprYG6oYPssT1w0N3pndotjWwbmgfi7mAe+hD6V2r+IbCXTLSWNNPMecwgPAotOz2uqTptewm4CLRx+NpXqv4hsJdMlKXlWAGDoNgLlTv8Lqk6VXvwAWjA3dHLJDetYpvKNwlI50eGAWg98rPwD27oHKbxxVdRuU23KeeYMwFaVy6M71rFd8IeV2AyFyMnXqDURei8EN/CGk4SMdUweobOZ5Ti7t41utSJEvoyF0y0oLeo7QHK9Jy9KXp9BauZdWIRmWS1FC4S0YqHTrOucK1XpcxK+GSDRTTy/kzp7wuRbKAwl0yTn/vWVbSw1jRFV6XMisLK6I3W3Uefc3jSiQbKNwl43Qe2QvAgsr0vjN1qlV11wEw0HrA40okGyjcJeP0tUb7lFmx5hqPK5mdopWV9FKInVEfM5J8CnfJPF2HGHQFrKxM765+p7JAgI7cGpb0q48ZST6Fu2Scwv4m2nOqsUDmvX0HF9dRMdaCi0S8LkV8Lq5Ph5ndamZHzKzZzO6/xPxPmtnB2M+LZpaeoxVLxnORCOWjJ+hbnFlH7W8p3UShDXG6TUfvklwzhruZBYEHgNuAjcDdZjZ1dIQTwHucc1cBvwc8mOhCRQDOnm5jKYO4kg1elzInS6quAqD72D6PKxG/i+fIfRvQ7Jw77pwbBR4B7pq8gHPuRefcxACRLwMViS1TJOpUU/RKmUVVmfmfw7L10StmLra97nEl4nfxhHs50DbpeXts2nQ+Azw9n6JEpnOhPRqKq+q2eFzJ3CxeWsRpisk5e9jrUsTn4ulbxi4xzV1yQbP3Eg33m6eZfy9wL0BpaSn19fXxVTnF4ODgnF+b6bK17RPtttbX6HFLeOPNJngzM89bFwYrWNLfFNd+zPb9nW0S2e54wr0dqJz0vALonLqQmV0FfBO4zTl3yd6RnHMPEjsfv3XrVrdz587Z1gtAfX09c31tpsvWtk+0u+lfv8ip/NqM/hu8dGQzqzsfofqmG2fsGyfb93e2SWS74zkt8ypQZ2Y1ZpYLfALYNXkBM6sCfgB8yjl3NCGViUwRCYepGDvJhSXrvC5lXkIrN5Fr43QeP+R1KeJjM4a7c24c+BzwDNAIfN85d8jM7jOz+2KL/W+gCPgLM9tvZnuSVrFkrc4TjRTYKMGVm7wuZV6W10S/L+g5ritmJHni6s/dOfcU8NSUad+Y9PhXgV9NbGkib3fm+H4qgCXVmXmlzITyuqsZdwHGOnXkLsmTebf4SdYa7oheKVOxLrP6lJkqv2AhHcEy8s7pihlJHoW7ZIzcs4fpsFIWFC7xupR5O7uglpKh416XIT6mcJeMUXzxGN0Fa7wuIyFGlm+gLHKai4N9XpciPqVwl4wwPj5KebiD4WXrvS4lIfLLNxEwR0eT+naX5FC4S0YY6mklZBFyyzZ7XUpCFK+5FoDzLfs9rkT8SuEuGSHS2wpAUW1mdjswVVn1BoZcLpHTb3pdiviUwl0yQv5AC6MuSPmazBpabzrBUIj2UBULe494XYr4lMJdMsLy4TY6gpUz3q6fSXoL17Jy5ITXZYhPKdwlI1SEWzlbuNbrMhIqXLKBEs7T23Pa61LEhxTukvYG+s6xih7GivxxpcyEBRXRU0wdR1/zuBLxI4W7pL2OI9EBOgrK/XG+fcLKuugVM4OtuhxSEk/hLmmv7+RBAFasvdbjShKrZNVq+lgIZxq9LkV8SOEu6a/rEBdcPisr/XXO3QIBOnJqWNKfmYOOSHpTuEvaK+xv4mSgkkAw6HUpCTewZB3loydwkYjXpYjPKNwlrblIhMrRZvIDYWhr8LqcxFuxgUU2xOAPf9uf7RPPKNwlrfUdeJLFXKQ6fAIeutN3AVi6uACAwr1/4cv2iXcU7pLWLu77RwACOAiPQstujytKrLJQP86B+bR94h2Fu6S1sxfHcA4iBCCYC9U7vC4poXLXvw9nhgNftk+8o3CXtGYDp+gIltNS80m4ZxdUbvO6pMSq3MbhRTcxQg7uU4/7r33iGYW7pK3w+DjVw0foWHY9ras/5tvgG6x+H/mM0T5a4HUp4iMKd0lbbU37KbQhApXXe11KUpVccRMAp9/8mceViJ8o3CVtnWmMht2KDTd7XElyVa2/jgsun0jbq16XIj6icJe0FenYSz8LqFzrrz5lpgqGQrTkrWNZ7xtelyI+onCXtFXc+zon867w5Z2pU/UXbaF67BjDQxe8LkV8QuEuaWnowgCrx1sYLPHHsHozyaveRq6FaTn0steliE8o3CUtnXj9Z4QsQkG1P6+Qmapic/R7hd6jL3lcifiFwl3SUn9zNOQqr8yOm3pWlNfQzXJCpzRwhySGwl3SUu7p1+i0FRSVVnhdSsq0L9zEqkF9qSqJoXCXtFQ2+CanCjd5XUZKjZZeQ7nr4lx3h9eliA8o3CXt9HSeZCU9jK3018hLM1m09kYA2t5Q52Eyfwp3STsT4ba07kaPK0mt6itvIuyMiyd0M5PMn8Jd0s5wSwNjLkj15uwK94WLlnIyuJqFZ/Z7XYr4gMJd0s6iswdoCdWQv6DQ61JSrmfplawebtSwezJvCndJKxM9QZ5butnrUrxRfh1LuED7sde9rkQynMJd0kq29AQ5HfUQKYmicJe0ki09QU5HPURKoijcJa1E2vdkRU+Q05noIXJ5r07LyPwo3CWtFPdlT0+Q0+kv2sLqsePqIVLmReEuaePiYB+rx09mTU+Q01EPkZIIcYW7md1qZkfMrNnM7r/E/CvM7CUzGzGzLya+TMkGLW+8FO0JsuYGr0vx1ERnaeohUuZjxnA3syDwAHAbsBG428w2TlnsHPB54GsJr1CyxkRPkFVZ0hPkdErKqumiSD1EyrzEc+S+DWh2zh13zo0CjwB3TV7AOdftnHsVGEtCjZIloj1BlrJ8RbnXpXiuQz1EyjyF4limHGib9LwdmNP/m83sXuBegNLSUurr6+eyGgYHB+f82kzn57avHzzE8Zx1HL1E+/zc7ku5kFPFte4FXjnTmVXtnpBt+3tCItsdT7jbJaa5uWzMOfcg8CDA1q1b3c6dO+eyGurr65nrazOdX9veu28XSznL2dJqNl+ifX5t93QO5V6EHz3M+vbvsWn7l6AyO0akmpBt+3tCItsdz2mZdqBy0vMKoDMhWxcBaGtg8T//CgAbOx+FtgaPC/LemuJ8nIMNffXw0J36m8isxRPurwJ1ZlZjZrnAJ4BdyS1LskrLbiwS/brGXBha1J95ftc+sNgHNDyqv4nM2oynZZxz42b2OeAZIAh82zl3yMzui83/hpmtBPYAi4GImf0GsNE515/E2sUnwpU3ESB6rs+CuVCd3VfLAFC9g4iFCLpxIoEQAf1NZJbiOeeOc+4p4Kkp074x6fFpoqdrRGatqWuAK4COitsp/8B/ybrzy5dUuY2e93+d0mfu483yj7NZfxOZJd2hKp47v/8JRl2QRR/7uoJ9ktIb7+YE5biuQ16XIhlI4S6eK+uq50j+1SxeWuR1KWnn6MKtXDF8kL7zPV6XIhlG4S6eams6wOpIOxdr3u91KWlpeNUN5FiYphcf87oUyTAKd/FUxyvR0Kq68aMeV5KeFpWt5xyL4fDTXpciGUbhLp5a3PoTjgVrWLV6vdelpKVAIETz0ptZN/ASY6MjXpcjGUThLp7p7TnN+pE3OLPqFq9LSWs5G+9gMRc50vCM16VIBlG4i2eafvYDguYouu7DXpeS1tbf9CFGXA6DB//Z61IkgyjcxTPBpqfpZjlrrnqX16WktQWFSzi84FqqzjyPi0S8LkcyhMJdPDEyfJH1Aw2cKNqR1UPqxWu49gOUuS5aDu/1uhTJEAp38cSRl59moQ2Tv+mDXpeSEWrf9YsAdDXokkiJj8JdPDH0xpNcdHmsv/EOr0vJCCVl1RwNrWNZ+0+8LkUyhMJdUs5FIlT3PM+RwuvJL1jodTkZ42z5LawfP0LP6VavS5EMoHCXlDv2+kuUcpaxtbd6XUpGKb0+eqPX8Z896nElkgkU7pJyZ/Y+TtgZa9+lu1Jno2bj9XTaCnKP6Xp3mZnCXVKupONZjuZu1EDYs2SBAK3F7+GKC3sYujDgdTmS5hTuklLn9v6AteFjhFds8rqUjFR41YfItzHOf/+zGnpPLkvhLqnT1sCSJ38NgI2nn1A4zcEV5UU4YNXJJzS2qlyWwl1SJnysnoAbByAQGde4oHMQ6ngFAAOcxlaVy1C4S8oc7TwXDSUCoLFS56Z6By6Yh3MQcehvKNNSuEtKhMfHWdi8iw5biXvvl+CeXRpSby4qtxH49JN05Ncx7qAnZ6XXFUmaUrhLSuz7l29TFemga9v9BN7zRQX7fFRuw338bwgRofmxP/S6GklTCndJukg4TMlr/4+WQBVb3v/vvS7HFyrXXsm+Jb/AVad/wNmudq/LkTSkcJek2/fMd1gdaePsdV9QD5AJVHLHb5PPKEcf/6rXpUgaUrhLUkXCYZbv+XNOBirY8oFPe12Or6xev4V9i9/L1Z3f4/yZU16XI2lG4S5JdeAnD1MTOcmZaz5PMBTyuhzfKbr9S+QzyuHHv+J1KZJmFO6SNJFwmMUNf0qblXHNbZ/xuhxfqt6wlf2L3s2V7d+j72yX1+VIGlG4S9IcePa7rAmf4PTVn9NRexItvfVLFNoQbz7+R16XImlE4S5J4SIRFr3yJ7TbSq6549e8LsfXajffwL6FN7Op7R/oO9/jdTmSJhTukhRNT3yZteFj9K/5EKGcXK/L8b3FH/gSi7lI73fuVn8zAijcJQn6DjzJ2gN/hAM2tPydwiYF1iwLESFAVW8Dke/cob+5KNwlsS4M9DL2xOcxF+3cysJj6twqFVp2Y0b0JzxK32s/8Loi8ZjCXRJmdGSYYw98lKXhc7hgDlhQHYSlSvUOLJiHs+hHemTf93XnapbTJQySEJFwmIMPfJKtw3tp2PJ7bNu6PXrEXr1D/cikQuU2uGcX1rKblqF8Sn/2O7Q/eBd5X/gJhYuXeV2deEDhLvPmIhEaHvzPbO//CS/VfJYbP/L56AyFempVboPKbVQD+4NFbH7hP9H4wEdY/1//hdy8fK+rkxTTaRmZn7YG2r5+B9u7HuHlko+z/VO/73VFAmz5N5/gtS2/y5Uj+zj253cQef5r+pI1y+jIXeaurYHwt2+jyo0TIcC2D34GC+h4IV1s+8iv83r3Pq489SjuuT243QWY+tHPGvokypw0H/hXTj30aYKxYfPMjEDrix5XJVNt3rARR/TKJcaHOPXsA7hIxOOqJBUU7jIr7c1vsPePP8zax+6gcPwcYQvhLIjpqpi0ZDXvhlABDosOrN3yOIe/fDOHX/mR16VJksV1WsbMbgX+HAgC33TOfWXKfIvNvx24CHzaOfdagmsVD53f9wSDP/4Kqy4cZjk5vFT5K2z62G8THDimq2LSWeW26KmYlt2Ey2/g1T0vsebNr1P89MfZ//yNrNj+S5TZOe0/H5ox3M0sCDwAvA9oB141s13OuTcnLXYbUBf7uQH4y9hvySRtDW8F9ZngCtoOPMdoy8us6nmRqvFWlgIRCzL64W9y45Y7o69ZWqRQSHexq2hCwA21N3Nx8F5eevSPuOrEX7Pwpy/hgAgBjtT9Gotu+BQVtZuwjj36RzvDxXPkvg1ods4dBzCzR4C7gMnhfhfwt845B7xsZkvNbJVzLvEjCLQ1UHXyn6BtwTvfdJPC6W3zkj09hdu4ZNvjXJeLRBgdHWb8xIu44/X0FlRzPrickb4u8rv3s6H9HwkQxgElsZ9hl0NfcDlY9Lxt0GDpQNM794tkjAWFS7jxnj9g+Jlc3Et/iuEIEGFj019B01/RzwIKGcaIELEQR9b8B8bKtlKwdBXLw90Unj+M1ewgp3bHz0fWSvBnZj7v83i3kXaZkGDxhHs50DbpeTvvPCq/1DLlQGLDva2ByN/cTk1kjMi3HqbLShixPADy3Ail7gyGw2FvzUv29ERte9TyyJ1mOvDWvBockW89TLcVMUYuuW6UEnowwAFnWUaYELkMs4yB6N/NwThBcixMHpAX+3MWAhWT/rzORW9fB2gt3sHwTV+ketN2SrsPwkN3QnhUd5z6SP7G2+DVv3xrv556z1dp6zpLWfN3WTR8NPqPuRtnY/NfQ/Nfv/3FL34NgDEXZJwg+Tb61qxzLGGUXEKMU8T5t96bZyhmJBB99+W4UUpdz7SfgYn3+eU+G6n4XCYvE8BC+ZDEq5fiCXe7xDQ3h2Uws3u4f8u0AAAFtElEQVSBewFKS0upr6+PY/M/V3Xyn6iOjMc25hi1XE4FygAoi3Rgzr1jXrKnz3fbDseo5dERLKc8fOnpwDvmDduCn78mwlsfoMHgEtpC1VSMt7IsPIABEYMTOes4lr+Z2tGj1I0cIABEMA4vfx/Hy++iONzN9Ye/jEXGcYEQPZUfpL8fOl96GYDFV/4OS3vfoHfpZvqPXYRjs9t38zU4ODjr94sfJLvdb9uv4dVQvJrunAWsOvC/3novvFz7G5y15aztfpor+l4ggCOC0ZS3mZactdSONLJm7HDsPQV9weW0B6uoDJ+kOPzzcB8KRN+zxjvfzzN9BiZPT8XnMhWZEBkfoeWnf0vr6otv7Y9E7m+Lnkm5zAJmNwK/45z7QOz5bwE45748aZm/Auqdc9+NPT8C7LzcaZmtW7e6PXv2zK7atgZ46E4i4yMEQnlv/1cvNu+to8uJecmenoptX67tid5Gmp5nra+vZ+fOnV6XkXKetftS74UUfmbm/T5PZL2p2HZMPPvbzPY657bOuA+dc5f9IXp0fxyoAXKBA8CmKcvcATxN9AByO9Aw03qvu+46Nyetr7hj3/msc62vXHKee+Fr75yX7Okp3MYl257IbaSp5557zusSPJF27U7RZyYh7/NE1puKbbv49jewx82Qr865mcM9ui5uB44Cx4AvxabdB9wXe2xEr6g5BrwObJ1pnXMO9zj/AH6VrW1Xu7OL2j29eMM9ruvcnXNPAU9NmfaNSY8d8Nl41iUiIsmnO1RFRHxI4S4i4kMKdxERH1K4i4j4kMJdRMSHZryJKWkbNjsDnJzjy4uBngSWk0myte1qd3ZRu6e32jlXMtOKPAv3+TCzPS6eO7R8KFvbrnZnF7V7/nRaRkTEhxTuIiI+lKnh/qDXBXgoW9uudmcXtXueMvKcu4iIXF6mHrmLiMhlZFy4m9mtZnbEzJrN7H6v60kWM/u2mXWb2RuTpi03sx+bWVPs9zIva0wGM6s0s+fMrNHMDpnZF2LTfd12M8s3swYzOxBr9/+NTfd1uyeYWdDM9pnZk7Hnvm+3mbWY2etmtt/M9sSmJazdGRXukwbrvg3YCNxtZhu9rSppvgPcOmXa/cCzzrk64NnYc78ZB/6bc24D0bEBPhvbx35v+whwi3PuamALcKuZbcf/7Z7wBaBx0vNsafd7nXNbJl3+mLB2Z1S4M2mwbufcKDAxWLfvOOdeAM5NmXwX8FDs8UPAh1NaVAo45045516LPR4g+oEvx+dtj3XVPRh7mhP7cfi83QBmVkF0wJ9vTprs+3ZPI2HtzrRwn24g7mxR6mJDF8Z+r/C4nqQys2rgGuAVsqDtsVMT+4Fu4MfOuaxoN/BnwG8SHYJ1Qja02wE/MrO9sfGlIYHtjmuwjjQS10DckvnMrBB4FPgN51y/2aV2vb8458LAFjNbCjxmZpu9rinZzOyDQLdzbq+Z7fS6nhR7l3Ou08xWAD82s8OJXHmmHbm3A5WTnlcAnR7V4oUuM1sFEPvd7XE9SWFmOUSD/e+dcz+ITc6KtgM453qBeqLfufi93e8C7jSzFqKnWW8xs4fxf7txznXGfncDjxE97ZywdmdauL8K1JlZjZnlAp8AdnlcUyrtAu6JPb4HeMLDWpLCoofo3wIanXN/MmmWr9tuZiWxI3bMrAD4BeAwPm+3c+63nHMVzrlqop/nnzrnfhmft9vMFprZoonHwPuBN0hguzPuJiYzu53oObog8G3n3B94XFJSmNl3gZ1Ee4nrAv4P8DjwfaAKaAU+7pyb+qVrRjOzm4HdRAdanzgH+z+Jnnf3bdvN7CqiX6AFiR50fd8597tmVoSP2z1Z7LTMF51zH/R7u82slujROkRPj/+Dc+4PEtnujAt3ERGZWaadlhERkTgo3EVEfEjhLiLiQwp3EREfUriLiPiQwl1ExIcU7iIiPqRwFxHxof8PrBCEgw6uPOsAAAAASUVORK5CYII=\n"
]
},
"metadata": {
"needs_background": "light"
}
}
],
"execution_count": 4,
"metadata": {
"inputHidden": false,
"outputHidden": false
}
},
{
"cell_type": "code",
"source": [],
"outputs": [],
"execution_count": null,
"metadata": {
"inputHidden": false,
"outputHidden": false
}
}
],
"metadata": {
"kernel_info": {
"name": "python3"
},
"kernelspec": {
"name": "python3",
"language": "python",
"display_name": "Python 3"
},
"language_info": {
"name": "python",
"version": "3.6.7",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"nteract": {
"version": "0.15.0"
}
},
"nbformat": 4,
"nbformat_minor": 1
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment