Skip to content

Instantly share code, notes, and snippets.

@ssugiyama
Last active May 31, 2017 15:37
Show Gist options
  • Save ssugiyama/ee185ee6de94b91e088698791145c762 to your computer and use it in GitHub Desktop.
Save ssugiyama/ee185ee6de94b91e088698791145c762 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Visualization of $(x+y\\mathrm{i})^\\mathrm{i}$"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"from matplotlib import pyplot as plt\n",
"from matplotlib import ticker\n",
"import numpy as np\n",
"import seaborn as sns\n",
"%matplotlib inline \n",
"\n",
"x = np.arange(-3, 3, 0.01)\n",
"y = np.arange(-3, 3, 0.01)\n",
"X, Y = np.meshgrid(x, y)\n",
"Z = (X+ Y*1j)**1j"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## absolute values"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAATYAAAD7CAYAAADgvbh3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFGtJREFUeJzt3X+sV/V9x/HnFRGCgd4bG0DpnayavY2Z60hIoEURF4w/\nsq7V6LowFZXaYGoiuhZEpQmJzri0qxiLNiKLWp2NWBrXpYXMaRt/1HaNW3Sm71kMc2MZbDeCWn8w\n9O6P77n45XLv98f5+Tmf83okJvd74XvOx7Y8+z7nc8+XgdHRUUREYnJM1QsQEcmbwiYi0VHYRCQ6\nCpuIREdhE5HoKGwiEp1j077RzKYA9wMGjAKr3f2VvBYmIpJWlont8wDuvgS4Fbg9lxWJSBTMbL2Z\nvWBmvzKzVW3fX2FmLxR57tRhc/cfAl9JXp4M7M9lRSJSe2a2DPgcsAQ4GxhOvr8AWAUMFHn+1Jei\nAO5+yMweBC4CLun0e3fv3l3rRxzefPPNqpdQayMjI1UvIXrLly/PHIsZM2b0/Of03Xff7XS+84CX\nge3ALODrZnYC8JfAGlq3sQqTefPA3VcCvwfcb2bHZ1+SiETgk8BC4FJgNfC3wFbgRuDtok+eOmxm\ndrmZrU9evgt8lPwTpaGhoaqXIFInI8AOdz/o7g58itYAdC/wGHC6md1V1MmzTGw/ABaY2c+AHcAa\nd38vn2WJSM09C5xvZgNmdhKwB/h9d18G/BnwqruvKerkqe+xuftvgT/NcS3BGxoa0r02kR64+4/M\nbCnwC1oD1Ffd/cOyzj9Q1scW1X3zYIzC1j9tHJQjsM2DSunJgz7pXptI+BQ2EYmOwpaCpjaRsCls\nIhIdhS0lTW290caBVEFhy0BxEwmTwiYi0VHYMtLUJhIehU1EoqOw5UBT28S0cSBVUdhEJDoKW040\ntYmEQ2ETicCePXuqXkJQMn00uBxJH2skVcgzaieddFJux6qSJrac6ZJUyqRJbWIKmxRCO6LFU9Qm\np7AVQFObFE1R60z32ERqREHrjSa2gmhqk7wpar1T2ERqQFHrj8JWoKZObdo4yJei1j+FrWBNjZvk\nQ1FLR5sHIgFS0LLRxFYCTW3SD0UtO4VNJCCKWj4UtpI0ZWrTxkF6ilp+FDaRAChq+dLmQYn06R8y\nnoJWjFRhM7OpwFZgPjANuM3dn8xxXSLRa0LUzGwRcKe7LzOz2cD9wBAwBbjC3XcVcd60l6KXASPu\nfhZwPnBPfkuKW1PutUlnDYnaWmALMD351l8Bj7j7UuBW4LSizp02bI8DG5KvB4BD+SynGRS3ZmtC\n1BK7gIvbXi8BPmVm/wD8OfBMUSdOdSnq7u8AmNlMYBut+krDaUe0u9CjNjw8nNux3P0JM5vf9q35\nwJvuvtzMvgGsA76R2wnbpN4VNbNh4GngYXd/NL8lNYOmtmbZs2dP8FErwQgwdi/+74CFRZ0oVdjM\nbA6wE1jn7lvzXZJIXBS0w54FLky+Xgr8a1EnSvvjHjfT2tnYYGZj99oucPf38llWM+jHP+KnqB3h\nL4AtZnYtcABYUdSJ0t5jux64Pue1iERFUQN33w0sTr7+d+DcMs6rJw8qFsu9Nm0cHElRq5aePBDJ\nkYIWBk1sAYhlams6RS0cClsgFLd6U9TCorCJZKSohUdhC4imtvpR1MKkzQPJrIk7ogpa2DSxBUZT\nW/gUtfApbCJ9UNTqQWELkKa2MClq9aF7bCJdKGj1o4ktUHWZ2mLfOFDU6klhC1hd4hYrRa2+dCkq\nMoGmRm3evHlVLyEXmtgCp6mtfE2NWkw0sYkkFLR4aGKrgVCntpg2DhS1uChs0niKWnwUtpoIdWqr\nO0UtTgqbNJaiFi9tHtSI/larfCho8dPEVjO6JM1GUWsGhU1SqeOOqKLWHApbDWlq65+i1iwKm0RP\nUWsebR7UlDYSulPQmqu0iW1wcJDBwcGyTicNp6g1W+mXomOBU+Syq+peW+gbB4qaZAqbmS0ys2fS\nvl+Rk7wpamExs9lm9h9mdlqZ500dNjNbC2wBpuexEAUuHe2QfkxRC4uZTQW+C7xX9rmzTGy7gIvz\nWsgYTXH9a3rc9uzZo6iF6ZvAfcB/lX3i1Lui7v6Emc3PcS1HGYvb/v37izyN1JiClq/Zs2fnchwz\nuxL4H3ffYWbrczloH2rxc2ya4rora2oLaeNAUQva1cC5yT34PwQeMrO5ZZ28dj/H1h43TXLNpaiF\nzd2Xjn2dxG21u/93WeevxcQ2GU1xR2rKvTZFTbrJNLG5+25gcT5LSU9TXDMoaPXk7svKPmetJ7aJ\nNH2Ki3VqU9SkH7W7x9YrTXHxUNSkX9FNbBNp2q5qUVNbFTuiipqk0YiwtWtK4GK4JFXUJK1oL0W7\n0aVquBQ0yapxE9tEYp3i6ji1KWqSh8ZObBPRFFctRU3yooltErFMcXlNbUVvHChqkidNbF1oiiue\noiZ5U9j6UNfIhfr3IyhoUhRdiqZUt0vV0DYSFDUpkia2jOo6xVVJUZOiaWLLUehTXNqpLc+NA0VN\nyqCJrQCa4iamqIVvzpw5VS8hFwpbwUKLXBUbCQqalE2XoiUK/VK1CIqaVEFhq0DVnzZS1g6poiZV\nUdgqFusUp6hJlRS2QJQ9xfU6taXZEVXUpGraPAhQWX+fat4bCQqahEJhC1hoO6qdKGoSEl2K1kRR\nl6p5bCQoahIaha2GQtpwUNQkRLoUrbG8LlUnu9fWaeNAQZOQaWKLRJlTnKImoVPYIpP2Xlyv99oU\nNakDXYpGLO9dVUVN6kITW0P0MsV1mtoUNamT1BObmR0DbAY+A3wAfNndf5PXwqQY3aa4sY2EsY0D\nBU2yMLNFwJ3uvqzM82aZ2L4ITHf3zwI3Ad/KZ0lSlm5TnKImWZjZWmALML3sc2e5x3Ym8BMAd/+5\nmS3s9JsXLuz4yxKga665puolSL3tAi4GHi77xFkmtlnAgbbXH5qZNiMisnfv3qqXIDXm7k8A/1fF\nubOE6C1gZtvrY9z9UMb1SMWGh4cBmDdvHvBx3GL5yGjpbO7cuVUvIRdZJrbngAsBzGwx8HIuK5JK\nDA8PH47amH379h3+WtOb1EmWiW07cK6ZPQ8MAFflsyQpy/iQjRmb1sbbu3evJjephdRhc/ePgNU5\nrkVKMFnMeqW4ST/cfTewuOzz6mZ/Q/QatMmmtXaKm4ROYYtY1umsE20qSMj0SFWEJtoI6MVE01r7\nBsJEtKkgIdLEFokip7NudGkqoVHYaizPmPVyb60TxU1CorDVUN7TWdaojVHcJBQKW01UeanZD20q\nSAi0eRC4tBsBveplWuu2gTARbSpIlTSxBagu01k3ujSVqihsgagiZnndW+tEcZMqKGwVi2U660Rx\nk7IpbBUIIWZlTGvttKkgZVLYShRC0KD8qLXT9CZlUNgKFkrMstq3bx+zZ8/O5ViKmxRNYStIqEGr\nclprp7iFKZT/fWSlsOUo1JiFSnGToihsGdUpZiH+v7E2FaQICltKdQpaHWh6kzwpbH2oc8zymNby\n3ECYiOImeVHYelDnoNWN4iZ5UNgmEVPMQry31oniJlkpbG1iitmYukVtjDYVJAuFjTiDFgtNb5JG\nY8PWhJjVdVobT3GTfjXugyaL/uDG2KX50Mk86IMrpR+NmNiaGLJYprV2mtykV1GHrYlBi502FaQX\n0YVNMYtzWhtP05t0kilsZnYRcKm7r8hpPakoZs2kuMlkUofNzDYB5wH/nN9y+qOgHa2Maa3oR6v6\nobiFycymAluB+cA04DZ3f7KH9y0GNgGHgJ3uvrHt1xYBd7r7sm7HybIr+jxwbYb3pzK2q6moHa0J\nl6AT0Y5pkC4DRtz9LOB84J4e33cfsAI4E1hkZgsAzGwtsAWY3stBuk5sZrYKuGHct69y9++b2bIe\nF5uZQiadaFMhOI8D25KvB4BDZnYGcHfyegS42t0PjL3BzGYB09x9V/J6B7AceAnYBVwMPNzLybuG\nzd0fAB7o9d8mT4pZ75o6rY2nS9NsTjjhhFyO4+7vAJjZTFqBuxW4n1bMXk0GprXALW1vmwW81fb6\nbeDTyfGeMLP5vZ4/uF1RxUyyUtzCYGbDwHZgs7s/amb3ApvNDGAq8JqZXQdckrxlJTCz7RAzgf1p\nzh1M2BS09KqY1kLaQJiI4lYtM5sD7ASuc/enkm87cIW7v2FmS4AT3X0bbfffzOygmZ0CvE5rc3Ij\nKWQKm7s/AzyT9v2KmRRJcavUzcAQsMHMNiTfuwF4yMyOBUaBVRO8bzXwCDCF1q7oi2lOPjA6Oprm\nfX079dRTD59IQctPlffWQp7YxmtC4NatWzeQ9RgvvfRSz0FYsGBB5vMVpbRLUcVMqqTprVka9+ke\nMdFOaH/0827NobDVlKKWjuLWDAqbpFbVZ7NltXfvXgUucgpbDWlay4fiFi+FTRpNcYuTwlYzmtby\np7jFR2ETQXGLjcJWIyFOa3XdQJiINhXiobCJjKO41Z/CVhMhTmsxU9zqTWGrAUWtGopbfSlsIh0o\nbvUUzOexycTqMK2F/tlsWTXpY8eHhoaqXkIuNLGJ9EjTW30obAGrw7TWNIpbPShsIn1S3MKnsAVK\n01rYFLewKWwiKelJhXApbAGq47QW06NV/VLcwqOwBaaOURPFLTQKm0hOFLdwKGwB0bRWf4pbGBQ2\nkZxpU6F6ClsgYpjWmryBMBHFrToKm0iBFLdqKGwBiGFak8kpbuVT2ERKoLiVK9XHFpnZJ4DvAbOA\n44Ab3f2FPBfWFJrWmqNJH38EYGbrgT+h1YjN7v5AD+9ZDGwCDgE73X2jmU0FtgLzgWnAbe7+ZKfj\npJ3YbgSecvezgSuB76Q8TqPFGDVtIHTXhOnNzJYBnwOWAGcDwz2+9T5gBXAmsMjMFgCXASPufhZw\nPnBPt4Ok/aDJbwMftB3j/ZTHEWmkvXv3xj65nQe8DGyndWX3dTM7A7gbGABGgKvd/cDYG8xsFjDN\n3Xclr3cAy4F7gW3JbxugNc111DVsZrYKuGHct69y91+a2Vxal6Rruh1HjhTjtCb9CTFug4ODeR3q\nk8DJwB8Dvws8Ceyn1Y5Xk66sBW5pe88s4K22128Dn3b3dwDMbCatwN3a7eRdw5ZcFx91bZzU9zHg\na+7+027HEZGjhRi3nIwAv3b3g4Cb2fvAHwCbzQxgKvCamV0HXJK8ZyUws+0YM2nFEDMbpjX9bXb3\nR7udPO3mwenA48CX3P1f0hyjyTStSbtINxWeBa43s78GTgSOB14BrnD3N8xsCXCiu2+j7Z6ZmR00\ns1OA12ldzm40sznATuA6d3+ql5Onvcd2BzAd2JTU94C7fyHlsUSEuKY3d/+RmS0FfkFrk/KrwP8C\nD5nZscAosGqCt64GHgGm0NoVfdHMNgFDwAYz25D8vgvc/b3Jzj8wOjqa379NB+ecc045JwpcU6a1\nmP/WqqKljdu6desGsp57//79Pf85HRwczHy+ougHdEvUlKhJNk34cZCiKWwiAVLcslHYSqJpTfql\njz9KT2ETCZzi1j+FrQRNnNb0aFW+FLf+KGwiNaG49U5hK1gTpzUpjuLWG4VNpGa0qdCdwlYgTWtS\nJMVtcgpbQRQ1bSCUQXGbmMImUnOK29EUtgJoWpOyKW5HUthEJDppP7ZIJqFpTepsxowZVS8hF5rY\nRCQ6CluONK0dTTujUgWFTUSio7DlRNOaSDgUthwoaiJhUdhEJDoKW0aa1rrTBoKUTWETkegobBlo\nWhMJk8ImItFR2FLStCYSLoVNRKKjsKWgaa1/2hmVMilsfVLURMKX6mOLzOx44FFgCDgIrHT3PXku\nTETqy8ymAluB+cA04DZ3f7KH9y0GNgGHgJ3uvjHNsdJObNcAv3L3pcD3gLUpj1MrmtZEenYZMOLu\nZwHnA/f0+L77gBXAmcAiM1uQ5lipJjZ3v8vMpiQvfwfYn+Y4IhKtx4FtydcDwCEzOwO4O3k9Alzt\n7gfG3mBms4Bp7r4reb0DWA7cO/5Y3U7eNWxmtgq4Ydy3r3L3X5rZPwJnAOd2O87TTz890O33iEi1\njjvuuFz+nLr7OwBmNpNWlG4F7qcVs1eTrqwFbml72yzgrbbXbwOfnuRYHQ2Mjo5m+hcws9OAv3f3\nUzIdSESiYmbDwHZgs7tvNbMDwEvJL08FXgP+Cbgk+d5K4Mfufnry/uuBqe7+zfHH6nbutJsH64H/\ndPeHgXeAD9McR0TiZGZzgJ3Ade7+VPJtB65w9zfMbAlwortvo+2emZkdNLNTgNeB84CNkxyro1QT\nW3KiB4HpwBTgJnd/ru8DiUiUzGwT8CXg123fvgW4ndZANQqscvd/G/e+xcBdtLqy091vmeRYF7j7\ne5OdP/OlqIhIaPQDuiISHYVNRKKjsIlIdEr5m+BDfATLzD5B66mJWcBxwI3u/kKVawIws4uAS919\nRUXnPwbYDHwG+AD4srv/poq1tDOzRcCd7r4sgLWkelyowPVMofUzYkbrpvxqd3+lqvWEoKyJLcRH\nsG4EnnL3s4Erge9Uu5zDO0l3UO0k/UVgurt/FrgJ+FaFawHAzNYCW2jtwocg7eNCRfk8gLsvofXD\nq7dXu5zqlfIHyN3v4uP/sEN5BOvbwHeTr48F3q9wLWOeB66teA1nAj8BcPefAwurXQ4Au4CLq15E\nm8eBDcnXPT3iUyR3/yHwleTlyYTx56tSuV+K5vUIVolrmktrilwTwHq+b2bLylrHJGYBB9pef2hm\nx7p7ZX943f0JM5tf1fnHS/OIT9Hc/ZCZPQhcxMc/yd9YuYfN3R8AHpjk1/5o7BEsoLRHsCZbU/JQ\n7mPA19z9p1WvJxBvATPbXh9TZdRCNe4Rn0erXg+Au680s3XAi2Z2urv/tuo1VaWUS1EzW29mlycv\ng3gEy8xOp3VJscLdf1z1egLyHHAhHP4p8JerXU542h7xWdfLc4slrOfy5DFHgHeBj5J/GquUXVFa\nO0gPJpdgU4CrSjpvJ3fQuhm9ycwADrj7F6pdUhC2A+ea2fO07h+F8N9VaG6mtcO/wczG7rV1fMSn\nYD8A/sbMfkbr4fI1Fa4lCHqkSkSiox/QFZHoKGwiEh2FTUSio7CJSHQUNhGJjsImItFR2EQkOv8P\nWja6Umw+sgwAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x105cee908>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.contourf(X,Y,np.abs(Z), locator=ticker.LogLocator(base=2))\n",
"formatter = ticker.LogFormatter(2, labelOnlyBase=False) \n",
"plt.colorbar(format=formatter)\n",
"plt.gca().set_aspect('equal')\n",
"plt.show()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## angles"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAATAAAAD7CAYAAADto8gwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAE69JREFUeJzt3X+wXGV9x/H3TYIJjSHK3KlkGIVeG7+BEWk71kCAip2i\nwKiAhf5IQWAiv6wtEBkEC6WdwaZUkB9VWo0XpdWoAxJkQFEHhh8mYNEOnVDL18Q7dJC5UAMS+ZVg\nSPrH7g2bZX+cPb+f53xeM3dmd8/e53l2z3k++5yzz54ztnPnTkREQjSr6gaIiKSlABORYCnARCRY\nCjARCZYCTESCpQATkWDNSfuPZjYbWA0YsBM4290fyathIiLDpA4w4AMA7n6YmR0JfAo4Lo9GiUjY\nhg1wzOzPgfOA7cAG4KPuvmPUelLvQrr7rcCZ7bv7Ac+mLUtEorNrgANcQmuAA4CZ7QlcDrynvXwh\n8P40lWQZgeHu283sRuAE4MRBz12xYoWm/IsUaHJycixrGXvvvXfifvrMM8/0rc/dbzWz29t3uwc4\n24Bl7v5i+/4cYOuobYUcDuK7+6nA24DVZjY/a3kiEoeOAc4/A1/teHyHuz8FYGZ/Bbwe+H6aOlIH\nmJmdYmYXt+++COxo/4mIAP0HOGY2y8yuBI4C/tjdU+2hZRmB3QL8rpndB3wXOM/dX8pQnohEIsEA\n5/PAPOD4jl3JkY2VdTYKHQMTKVadjoG1R1tfAvYB9gD+EZhPa3fxR+2/+2l9QwlwrbuvHbW9mQ7i\ni4j04u4vAH8y4Cm5TKLXTHwRCZYCTESCpQATkWApwEQkWAowEQmWAkxEgqUAE5FgKcBEJFgKMBEJ\nlgJMRIKlABORYCnARCRYCjARCZYCTESCpQATkWApwEQkWAowEQmWAkxEgqUAE5Fg6Zz4IpI7M9sD\nuAHYH5gLXO7ut/V43heAZ9z9ojT1aAQmIkU4GXja3Y8AjgY+2/0EMzsLOChLJRqBiUgRbgJubt8e\nA7Z3LjSzZcBSWteHXJK2EgWYiOxywAEH5FKOuz8PYGYLaAXZJTPLzGwRcBlwAoMvvTaUAqxBJiYm\nCit7amqqsLIlTGb2ZmAtcL27r+lYdBIwDnyb1oVvf8PMHnX3L49ahwIsUkWG1Sj1KdiayczeBHwP\n+Ji739W5zN2vA65rP+80YEma8AIFWDTKDqykutulQGuMTwJvBC41s0vbj60G5rv7F/KqRAEWqLoG\n1jCd7VaYxcvdzwXOTfC8L2epRwEWkFBDqx+FmWSVKsCSTlKTfMQWXL0ozCSNtCOwmUlqp5jZ3sDD\ngAIsR00IrX5mXruCTIZJG2ADJ6lJek0Orm4KMhkmVYANmqQm6Si4+lOQST+pD+IPmKQmI1BwJacg\nk25pD+L3naQmySi40lOQyYy0Z6PonKR2T/tvzxzbFTWFVz70PkraY2CJJqnJ7tTh8qfRWLPpfGAl\nUXgVS+9vM2kmfsHUscqj0VjzaARWIIVXNfS+N4cCrCDqRNWamJjQOmgABVjO1HHqResibgqwHKmz\n1JPWS7wUYDlRJ6k3rZ84KcByoM4RBq2n+CjAMlKnCIvWV1wUYBmoM4RJX7TEQxNZU1IHCN/ExIQm\nvRbMzJYCV7j7kV2P/z7wGVrnE3wSONndt45avkZgKSi84qF1WRwzuxD4IjCv6/ExWlcoOt3dDwfu\nBPZLU4cCbETa4EUS+xnwoR6Pvw14GjjfzO4F9nZ3T1OBdiFHoPCKk3YlX7VkyZLcynL3b5rZ/j0W\njQPLgI8Bm4DbzexH7n73qHVoBJaQwituWr+lehrY5O7/4+6/prUL+c40BSnAEtDG3Qxaz6WZAl5v\nZr/dvn8E8N9pClKADaGNulm0votjZsvN7Ex3fxlYAawxs4eAx939jjRl6hjYANqYm0nHxPLj7o8B\nh7Rvr+l4/G7gXVnL1whMRIKlAOtDo69m0/oPgwKsB228AtoOQqAAExlAIVZvCrAu2mBFwqEA66Dw\nkl60XdSXAkwkAYVYPSnA2rSBioRHAYbCS5LRdlI/CjARCVamADOzpWZ2T05tqYQ+VWUU2l7qJXWA\n9TvbokjsFGL1kWUE1u9si8HQhigSttRnoxhwtkUJzOLFixM/d+PGjQW2JBw6Y0U9NPZ0Ok0efY0S\nWMP+V4EmVWpsgDVNltBKWq7CTMrWyABr0uirqOAaVlcTwky7kdXLFGCdZ1uUeikzuAbV34Qgk+o0\nbiJr7KOvxYsXVx5enerUliLEvj3VXSN3IWNU56DQaKx5zGwWcD1wMLAN+Ii7b+pY/nFgObAD+Ad3\nX5umnsaNwGJU5/DqFEo7R6VRWE/HA/Pc/VDgIuCqmQVm9gbgXOBQ4L3ANWkraVSAxbihhRYKobVX\nUjuc1gVrcfcH2f3CtS8A/wvMb//tSFuJdiEDFmoYaJeyvnL8kN8L2NJx/xUzm+Pu29v3Hwd+AswG\nVqWtpFEjsJiEGl6dYngNM2Ic3Wf0K2BBx/1ZHeF1DLAI+C3gLcDxZpbqGpGNCbCYNrCYOn5Mr0V2\nsw44FsDMDgE2dCz7JfASsM3dtwLPAm9IU0ljAiwWMXb4GF+TsBbYambrgauB881spZl90N3vBx4C\nHjSzB4CfAt9PU4mOgQVEHb3eNDP/Ve6+Azi76+FHO5ZfBlyWtZ5GjMBi2H2MPbxif31SjEYEWOia\n0rmb8jolPwowqRWFmIxCAVZz6tBhieFwRUiiD7CQN6imhldTX7eMLvoAE5F4aRpFTVUxChk0Wi17\nesDixYv1UyMZSgHWcEl3sTufp7lOg2k+WHmi3oUM+fhX0SYmJlK/P1n+dxQ6FibDRB1goSq64+YV\nPvqAkKopwBom79BRiEmVFGA1U+Toq6iwKTLEtBspgyjAGqLokZJGYlIFBVgDlBUuCjEpmwJMai/E\n3UiFeTk0D6xGiuioo3akRYsWveax6enpkerTHCgpi0Zgskuv8Jp5vN8ykSpFG2Aawo/2HiQJqKQh\npvdeyhJtgElyGl1JqFIfAxt26XAJQ6/wGh8f33V78+bNr3n+KMfEpJmG5YOZnQGcBWwHLnf329PU\nk2UE1vfS4RKu8fHxXce8Fi1atFuYiYygbz6Y2T7AXwOHAe8DVpnZ3DSVZAmwQZcOl0BdeeWVu93X\n7qWkNCgf3gWsc/dt7r4F2AS8I00lWQKs56XDM5QnIvEYlA/dy54DFqapJEvgDLp0uATqggsu4Kqr\nXj0aMD09/ZrjYBKvHOciDsqH7mULaF2de2RZRmCDLh0ugdq8eTMbNmxgenqaDRs2KLwkrUH58B/A\nEWY2z8wWAgcAj6SpJMsIbC1wVPvS4WPA6RnKkopMT0/3PM7VL7j0DaQk9Jp8MLOVwCZ3v83MrgPu\npzWI+ht335qmktQB1ufS4VIjU1NTiSaV9guxXs8TSaJPPjzasXw1sDprPdFOZNXv8UYzLJxGCS+9\n91IWfWsou8yEVOdoTKMuqbNoR2AhKuIyYmlGQ9PT07v+yqhPJC0FmNReiNeHVJCXQwHWAGV1JnVa\nKZsCrCGKDheFl1RBAVYzRe4uFRUyRYZXiLuPUh4FWMPkHTYaeUmVog6wUDtX0aOOqampzO9NHmWI\nZKV5YA02E0CjnAK6zNDS7qMMowCrqY0bN5Z2ObHOUOoVZhppjUbvV3kUYLKbunQ+jb4kiaiPgYVO\nnVhksOgDrC4jCklOwS1JRR9goWtaZw799eoDs1wKMBEJlgIsAKGPSpJqyuuU/DQiwGIY1sfeuWN/\nfVKMRgRYLGLt5LG8rhg+KEOjeWCBKXOCaxliCS8Zzsz2BL4C/Cata0Ge6u6/6HrOZ2hdFHcH8HF3\nXzeozMaMwGL6dIyl08fyOiSxc4AN7n4E8G/AJZ0LzexgYBmwFDgFuG5YgY0JsNiE3vlDb3+3mD4g\nC3Q4cGf79neAP+pa/gTwIjCX1tW7fz2sQO1CBizU3cnYwktey8xWAOd3PfwUsKV9+zlgYdfy7bR2\nHR9tLztjWD2NCrCk10kMSWghpvCqt7z6h7tPApOdj5nZLcCC9t0FwLNd//Zh4Engfe3lPzCzB939\n5/3q0S5kBDZu3Fj7YAihjWlp9zGxdcCx7dvH0Loyd6dfAs+7+yu0RmjbgPmDCmzUCCx2dR2NxRpc\nMrJ/AW40sx8ALwPLAczsn4CbgTXAYWa2HpgNfNXdfVCBjQuwGHcjO82ERR2CrAnBpdFXcu7+InBS\nj8cv7Lh79ihlNi7AmqLKIGtCcEk9ZAowMzsBOMndl+fUnlLEPgrr1BkmRYdZ04JLo6/qpQ4wM7uW\n1rcFD+fXHClSEWHWtNCSeskyAlsP3AqclVNbStWkUVgvvYJnUKgpqHan0Vc9DA2wPhPSTnf3b5jZ\nkYW0SiqhkJLQDA2wXhPSYtH0UZiko9FXfWgiq4gEq/EBpk9TGYW2l3rJNI3C3e8B7smlJSIiI2r8\nCAz0qSrJaDupHwVYmzZOGUTbRz0pwEQkWAqwDvqUlV60XdSXAqyLNlaRcCjARAbQB1q9KcB60EYr\noO0gBAqwPrTxNpvWfxgUYANoIxapNwWYSBd9cIVDATaENuZm0foOi86Jn4BOu9MMCq9yDDoVvZkd\nA1wGjAE/Bv7S3Xf2K0sjsIS0ccdN67cc7VPRr6JH9pjZAuDTwPvdfSnwGDA+qDwF2Ai0kcdJ67VU\n64Fz+ixbBmwArjKz+4Gn3P0XgwrTLuSItDsZF4XX7hYtWpRLOSlPRT8OvAf4HeB54H4ze8Ddf9qv\nHgVYCgqxOCi8ipPyVPRPAw+5+5MAZnYfrTDrG2DahUxJG3/YtP5q6T+Bt5vZuJnNAQ4BfjLoHxRg\nGagThEnrrV7MbKWZfdDd/w+4GPgu8EPgFnd/ZND/ju3c2fcbylytWLGinIoqoN3JcMQcXpOTk2NZ\ny3jiiScS99N99903c31ZaQSWg5g7RUy0nuKjAMuJOke9af3ESQGWo6mpKXWUGtI6iZcCrADqMPWg\nD5T4KcAKoo5TLb3/zaAAK5BGANXQe94cCrASqEOVQx8YzaMAK4k6V7H03jaTAqxk6mj50gdDs+nH\n3BWY6XCawZ+eQksgZYCZ2ULgK8BewOuAle7+QJ4NawIFWToKL5mRdhdyJXCXu78bOA34XG4taiDt\nBiWj90m6pd2FvBrY1lHG1nya02wakfWm0JJ+hgbYgDMrPmRm+9DalTyviMY1lYJMoSXJDA2wfmdW\nNLODgK8DF7j7vQW0rfGaGGQKLhlF2oP4BwI3AX/q7v+Vb5OkW2enjjHMFFqSVtpjYKuAecC1Zgaw\nxd2Py61V0lcsYabQkjykCjCFVT2EFmYKLcmbJrJGojsc6hBoCiwpmgIsUv3Co4hgU1BJEkknwJvZ\nLOAO4Fvu/q+DylSANYzCRio0MwH+GmsdPP8a8Hs9nnc58MYkBSrARKQsQyfAm9mJwA7gziQFKsBE\nZJfx8fFcykkzAd7M3g4sB04E/jZJPQowEcldygnwHwb2Be4G9gdeNrPH3L3vaEwBJiKlGDYB3t0v\n7Hju3wFPDgovUICJSHl6ToA3s5XAJne/bdQCx3buTHwl8UxWrFhRTkUiDTU5OTmWtYxt27Yl7qdz\n587NXF9WOqW0iARLASYiwVKAiUiwFGAiEiwFmIgESwEmIsFSgIlIsBRgIhIsBZiIBEsBJiLBUoCJ\nSLAUYCISLAWYiARLASYiwVKAiUiwFGAiEiwFmIgESwEmIsFSgIlIsFJd1MPM5gNraF0992XgVHd/\nIs+GiUhckuSGmZ0P/Fn77rfd/e8HlZl2BHYG8GN3/wNaF6i8cMjzRUQG5oaZTQB/ASwDDgHea2bv\nGFRgqhGYu19jZrPbd98CPJumHBFpjgS58ThwtLu/AmBmewBbB5U59LJqQy4RfjdwEHCUuz+c7GWI\nSOyy5IaZjQGfBha4+1mD6sl8XUgzWwLc4e5vzVSQiDRGv9wws3nADcBzwEdnRmP9pD2IfzHwc3f/\nd+B5YGAlIiLDcqM98voWcLe7X5GkzFQjMDN7E3AjrcuEzwYucvd1IxckIo3RLzfMbCWwqf3Y14AH\nO/7tYnd/oF+ZmXchRUSqoomsIhIsBZiIBEsBJiLBSvUt5Kjq+NMjM1tIazbwXsDrgJWDDhaWxcxO\nAE5y9+UV1T8LuB44GNgGfMTdN1XRlk5mthS4wt2PrEFb9qD1Vf/+wFzgcne/rcL2zAZWAwbsBM52\n90eqak+ZyhqB1fGnRyuBu9z93cBpwOeqbQ6Y2bXAKqodGR8PzHP3Q4GLgKsqbAsAZnYh8EVa317V\nwcnA0+5+BHA08NmK2/MBAHc/DLgE+FS1zSlPKR3F3a/h1Te1Lj89uhr4fPv2HIb8ZKEk64FzKm7D\n4cCdAO7+IPDOapsDwM+AD1XdiA43AZe2b48B2ytsC+5+K3Bm++5+1KN/lSL3XcikPyHIu94MbdqH\n1qjwvBq05xtmdmRZ7ehjL2BLx/1XzGyOu1fWSd39m2a2f1X1d3P35wHMbAFwM61RT6XcfbuZ3Qic\nAJxYdXvKknuAufskMNln2R/O/IQAKO2nR/3aZGYHAV8HLnD3e6tuT038CljQcX9WleFVV2b2ZmAt\ncL27r6m6PQDufqqZfQL4oZkd6O4vVN2mopWyC2lmF5vZKe27tfjpkZkdSGtXYLm7f6fq9tTIOuBY\nADM7BNhQbXPqpz2j/HvAJ9z9hhq055T2z3QAXgR2tP+iV8q3kLS+sbmxves0Gzi9pHoHWUXroPC1\nZgawxd2Pq7ZJtbAWOMrM1tM6vlOHdVU3n6T1jfqlZjZzLOwYd3+povbcAnzJzO4D9gDOq7AtpdJP\niUQkWJrIKiLBUoCJSLAUYCISLAWYiARLASYiwVKAiUiwFGAiEqz/B4U3+akoGhx5AAAAAElFTkSu\nQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x108a07f60>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.contourf(X,Y,np.angle(Z))\n",
"plt.colorbar()\n",
"plt.gca().set_aspect('equal')\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"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.6.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment