{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Visualizing event horizon and ergosphere of Kerr black hole\n", "\n", "### Importing required modules" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import astropy.units as u\n", "import matplotlib.pyplot as plt\n", "from einsteinpy.utils import kerr_utils, schwarzschild_radius" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Defining the black hole charecteristics" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "M = 4e30\n", "scr = schwarzschild_radius(M * u.kg).value\n", "# for nearly maximally rotating black hole\n", "a1 = 0.499999*scr\n", "# for ordinary black hole\n", "a2 = 0.3*scr" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Calculating the ergosphere and event horizon for spherical coordinates" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "ergo1, ergo2, hori1, hori2 = list(), list(), list(), list()\n", "thetas = np.linspace(0, np.pi, 720)\n", "for t in thetas:\n", " ergo1.append(kerr_utils.radius_ergosphere(M, a1, t, \"Spherical\"))\n", " ergo2.append(kerr_utils.radius_ergosphere(M, a2, t, \"Spherical\"))\n", " hori1.append(kerr_utils.event_horizon(M, a1, t, \"Spherical\"))\n", " hori2.append(kerr_utils.event_horizon(M, a2, t, \"Spherical\"))\n", "ergo1, ergo2, hori1, hori2 = np.array(ergo1), np.array(ergo2), np.array(hori1), np.array(hori2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Calculating the X, Y coordinates for plotting" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "Xe1, Ye1 = ergo1[:,0] * np.sin(ergo1[:,1]), ergo1[:,0] * np.cos(ergo1[:,1])\n", "Xh1, Yh1 = hori1[:,0] * np.sin(hori1[:,1]), hori1[:,0] * np.cos(hori1[:,1])\n", "Xe2, Ye2 = ergo2[:,0] * np.sin(ergo2[:,1]), ergo2[:,0] * np.cos(ergo2[:,1])\n", "Xh2, Yh2 = hori2[:,0] * np.sin(hori2[:,1]), hori2[:,0] * np.cos(hori2[:,1])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plot for maximally rotating black hole" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[,\n", " ]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAD8CAYAAACCRVh7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXlwXNd1p7/TjcZCgCBAEAQXENy1ULIlS7AsxfbEsR2JVlyWMuNUyTU1VhxNqeLIGacyqUSKUuUZ21MVx1XxUvESVeyJnHFGVpxkpPLI0ciO7MSxRQmUbIkiRREUNxDcQRACQWLrM3+c10IT6gZA9PJev3e+qlvovu81+vTbfveec+69oqo4juM4TiFSYRvgOI7jRBcXCcdxHKcoLhKO4zhOUVwkHMdxnKK4SDiO4zhFcZFwHMdxiuIi4TiO4xTFRcJxHMcpiouE4ziOU5S6sA0olRUrVuiGDRvCNsNxHKem2Llz52lV7Zxvv5oXiQ0bNtDX1xe2GY7jODWFiBxayH7ubnIcx3GKUjaREJG0iLwgIt8L3m8UkR0isk9EviMi9UF9Q/C+P9i+Ie9/PBDU7xWR28plm+M4jrM4ytmT+CSwJ+/954AvqOpW4CxwT1B/D3BWVbcAXwj2Q0S2AXcB1wDbga+KSLqM9jmO4ziXSVlEQkS6gV8D/ip4L8B7ge8GuzwM3Bm8viN4T7D9fcH+dwCPqOq4qh4A+oGbymGf4ziOszjK1ZP4IvCHQDZ43wEMq+pU8H4AWBu8XgscAQi2nwv2f6O+wGccx3GcEChZJETkg8BJVd2ZX11gV51n21yfmf2d94pIn4j0nTp16rLsdRzHcRZOOXoS7wQ+JCIHgUcwN9MXgTYRyaXYdgODwesBYB1AsH0ZMJRfX+Azl6CqD6lqr6r2dnbOm+brOI7jLJKSRUJVH1DVblXdgAWe/1lV/yPwNPDhYLe7gceC148H7wm2/7PaGqqPA3cF2U8bga3As6Xa5ziO4yyeSg6m+yPgERH5LPAC8I2g/hvA34hIP9aDuAtAVV8WkUeB3cAUcJ+qTlfQPqcUslm4cMHKxYswPg4TEzA5CVNTMD1tJZud+YwIpFKQTkNdHWQyVhoarDQ2QlOTvXaiw+QkjI3NnOdC5zqbBVUrIlbS6UvPdX29ldx5bmqybU6kEWvE1y69vb3qI64rSDYLr78O587Z39dfh9FRe2BU6tpJp6G52crSpdDaaqW5uTLf5xgXL9p5Hhmx83z+vJXJycp9Z0ODndeWlpnzvGyZi0cVEJGdqto7335+JpxLmZiAM2dgaMjKyMilvYFqMD1t3zsyAseOzdRnMtDWZqWjA9rb/WGyWLJZGB62c3z2rJXx8erbkeuZDA1dWt/SAsuXW+nogCVLqm+bA7hIOKp2g548aWVkJGyLijM5CadOWdm3z1wabW3Q2QkrV9prKZQk5wDWOzh50o7f0JCJcVQZHbVy+LC9X7Jk5jx3dlpv06kKLhJJJJu1B8WxY3D8eGXdCZVEdaYV/Oqr5u/u6oJVq+xhkkr41GS5BkDuPF+4ELZFi2dsDA4dspJKmVCsWWPnO5MJ27pY4yKRJM6ehSNHYHCwdoVhLiYm7PcdOWJuqFWrYO1ae6AkqYdx9iwcPWrnOQwXUqXJZuHECSuplAnFunXWMEjSea4SLhJxZ3LSHpqHDln3PSlMTcHAgJWGBujuhp4e83XHkfFxO8+HD1uwOSlks9ZTOnbMznNPD6xfb5lTTllwkYgro6Pw2mv2kIyy77kajI/D/v1WOjpg40brZcSh1XnmDBw4YO6kGs9ULJnxcYtV9fdb72LzZgt8OyXhIhE3hoftRjl+PGxLosmZM1aamkwsenpqz6edzZo76bXXop1oEBaqdv0fP24ZcFu2WKPAWRQuEnFheBj27rXsFWd+LlyA3bst4L1hA2zaFP1BfNPT5jbcv9/GNDjzc/YsPPecjb+48koXi0XgIlHrjI7Cnj3ec1gsU1PmnjhwwHzZW7ZETyympuDgQROHiYmwralNRkZMLNra4OqrYcWKsC2qGVwkapWJCes5HDrkvuhyMD1t7ptDh8wNtWVL+G6obNbs2bcvnllKYTA8DD/7mcUstm2LbyJDGXGRqDVU7cHxyivxTGMNm+lp61kcOgRbt5pghDHeYnDQeohjY9X/7iRw4oS5Zjdtgiuu8JH7c+BHppY4dw5efNFaQ05lmZy0mMWhQ9birJYve3gYdu0yX7pTWVTNhTc4CNde6/GKIrhI1ALZrAVY+/vdtVRtzp83X3ZnJ7zlLZWbZHBiwnoOuWkonOpx4YKd4zVr7BzX14dtUaRwkYg6587BCy/YvDtOeJw6BT/6kbmgtmwprwvqyBHrtXhQOlwGB+H0abjuOu9V5OEiEWX277fYQ7VnYXUKk81assDRo3D99ZaDXwpjY/CLX9iDyYkGExPWq+jpMReUTyToIhFJxset9+Drd0eT0VH4t3+zoPbVV19+ryKXfLB7t4+GjyqHD9vkiDfeaGMsEoyLRNQYGoKdO32wVNRRtZTZkyfhhhtsoZyFcPGiNQC89xB9RkfhJz+Bt77V5v5KKCU7VkWkUUSeFZFfiMjLIvLfg/qNIrJDRPaJyHdEpD6obwje9wfbN+T9rweC+r0iclupttUcBw9aDrcLRO2Qe5AsJKng2DGLa7hA1A7T0ybqu3Yl1u1bjujbOPBeVb0OuB7YLiI3A58DvqCqW4GzwD3B/vcAZ1V1C/CFYD9EZBu23vU1wHbgqyKSDIegqqW2vvRSYi/EmiabtcykZ54pPOhtetpiD319PralVjlwAHbsSGRyQckioUZuDupMUBR4L/DdoP5h4M7g9R3Be4Lt7xMRCeofUdVxVT0A9AM3lWpf5JmctIvv0KGwLXFK5fRp+PGPL40ljY7Cv/6rp7bGgdOnrdeYpKnYKVNMImjx7wS2AF8B9gPDqjoV7DIArA1erwWOAKjqlIicAzqC+mfy/m3+Z2Z/373AvQA9PT3l+AnhcPGiCYTP5FnTZLMzRcfH0R/vYHrjFliyhPSeXUh2GhGbmTyV8oSZmub8eROKm24qPbutRiiLSKjqNHC9iLQB/whcXWi34G+hSfx1jvpC3/cQ8BBAb29vbY4uGx0190QtLykZY7JZ8yyMj8/8nZyEySmYmrTX09NW3hyKUPjFPqaljrQWyF4SSKdsJoh02qaIypX6eqhvgIZ6m2fQBSWCTExY7LC311bDizllzW5S1WER+RFwM9AmInVBb6IbGAx2GwDWAQMiUgcsA4by6nPkfyZejIzYRZZA/2bUmJ42nb5wAcYuwIUxE4TxCYo0UcqAzgjMfGQy0NhopWkJLGmCJUt8qqHQmZ628RQ33ACrV4dtTUUp+VITkU5gMhCIJuD9WDD6aeDDwCPA3cBjwUceD97/LNj+z6qqIvI48Lci8ufAGmAr8Gyp9kWO4WHrQXgAs+qo2vi10VErY2Nw4SKVE4MyMBn0WmYPuK+vtxlCmpth6VITDu91VJls1tLV3/Y2W0s9ppSjPbIaeDiIS6SAR1X1eyKyG3hERD4LvAB8I9j/G8DfiEg/1oO4C0BVXxaRR4HdwBRwX+DGig8uEFVF1VzIIyNWzp+Pz9i1iQkruXkARUwwWltNNJYuDWfy2sShCs8/b69jKhSiNT5hXG9vr/b19YVtxvyMjMBPf+oCUWEmJ226q+Fh+xumKFhMYmr+HSuAiAlFWzu0t0VvHaXYIWKupzVrwrZkwYjITlXtnW8/92xWg/PnvQdRQcbHrUU9NASj54m0+6haqM70oA4fsiW929th+XJzTTllJtejqKuLXTDbRaLSXLxoQWpfWaysTE6aKJw5Y/EFZ25ywfnBQQuCd3TYCp7ewygjqjZg8pZbYpUe6yJRSaambByEp7mWBVUYPmdj1c4N+9Iai+XiRZvI9uhRW72zs9N6GB74LgPT03bPv+tdsVka1UWiUuRaFT5QrmQmJ20evZMn3WNXbnKZXocPW+9iZZel2TolkJtF4d3vjsUCRi4SlWLXLp/qu0TOn4fjx82t5L2GyjI9PSPEra3QtQrallk81lkEY2M2juKWW2o+zcxFohIcOmQzujqL4tw5mzDVO2HhkAt4NzXZAm0dHTX/nAuHoSGbtPO668K2pCRcJMrN0JD1IpzL5uwwDB5N3PxpkeXCBZv89OhRy+xcscLF4rI5fBja2mD9+rAtWTQuEuVkfNxGYPp035fFuXMwMODiEFUmJqxjPDho48VWrHA31GWxa5ctStXWFrYli8LbBeVC1RYn8QWDFszYmC3hvXevC0QtMDFhPYtdu6zX5yyQbLam1xJxkSgX/f0eqF4gk5PBw+ZljzvUIhcuwL5XTeDHPLt7YVy4YAtP1SAuEuVgeNiaw86cqMLxE7YI36lT+MjoGmdkxHoVhw7HZ06sinLsWE0uPuUxiVKZnrbh+J6jOSfnz1vvYWwsbEucsqJw4jgMnYENG2I10Lgy7Npl6WLNzWFbsmC8J1Eqe/a4Q30Oslk4cgRe3u0CEWcmJ2HfPis16nqvDtPT8POf11Sj0kWiFIaGrHnsFGR01BpOx47hrqWEcPasDQ0YGgrbkghTY88NF4nFks3WbCCq0qhauuTuPZ7slUSmpiyP47XXPFZRlFdeqZmutYvEYtm3z6cfLcDEhMXwBwbw3kPCOX3aepLujS3A9LR1uWqAkkVCRNaJyNMiskdEXhaRTwb1y0XkKRHZF/xtD+pFRL4sIv0i8qKI3JD3v+4O9t8nIneXalvFOH/emkrOJYyMwMue1urkMT4Ou3fbnFDOLE6etMnJIk45ehJTwH9V1auBm4H7RGQbcD/wQ1XdCvwweA/wAWz96q3AvcDXwEQF+BTwDuAm4FM5YYkcu3b5qOpZHD8Or+z1oKXzZlRtxPaBA37bvIlduyLvkytZJFT1mKo+H7x+HdgDrAXuAB4OdnsYuDN4fQfwLTWeAdpEZDVwG/CUqg6p6lngKWB7qfaVnRMnvFmURzYL+/cH6d/uXnLm4NQpc8V7QyKPCxfsBoowZY1JiMgG4G3ADqBLVY+BCQmQW9NvLXAk72MDQV2x+uigan1nB7Cb/ZVXbHU4x1kIo6PmkvSR2nn090c6w6NsIiEiLcDfA7+nqnN5pQtNDaZz1Bf6rntFpE9E+k5VcyqMw4c9WB2Q8zX74XAul4kJ2LPbY1dvMD0d6RkbyiISIpLBBOLbqvoPQfWJwI1E8DfnoxkA1uV9vBsYnKP+TajqQ6raq6q9nZ2d5fgJ8xPxE1lNxsZMIHzZbmex5G4nH08RcORIZFtc5chuEuAbwB5V/fO8TY8DuQylu4HH8uo/GmQ53QycC9xRTwK3ikh7ELC+NaiLBgcO+FMRu4737HG/slM6qtC/31JlE4+q+W4jSDnmbnon8J+Al0Tk50HdHwN/CjwqIvcAh4HfCLY9AdwO9ANjwMcAVHVIRD4DPBfs92lVjUY7Y2oq8sGlavD66/Dqq5FPxnBqCbVBd6pQLadAZMktx9jaGrYll1CySKjqTygcTwB4X4H9FbivyP/6JvDNUm0qOwcPmiM1wYyOukA4lSM3S0XiheLVV6G3N2wrLsFHXM/H9LQ1dRLM2Jj5j10gnEpy4IBnynHsWORiEy4S83HkSKJjERcvukA41WP/a7acbaKJmGvbRWIuVBPdi5icNIHwILVTNdSmRUv0fE8DA5EaN+EiMRcnTiT2as1m7WZNcCfKCYls1lzzib32slmLg0YEF4m5qKE538vNgQORc406CSK3iFFi3ZyHD0dmoisXiWKMjiY2gXtw0AOITviMjSW4nTY+HqzWFT4uEsU4dChsC0JhZAQGjoZtheMYQ0M1MZt2ZYjIM8hFohDZbLBqTrKYnAwSK3w2VydCHI7ujBWV5cyZSPxwF4lCnDiRyMFz+/d7JpMTQdSuzampsA0JgSNH5t+nwrhIFCICJ6baHD/hs3I60WV8PDLel+oyMGCp+CHiIjGbiQlbHSVBXLwIA8nTRafGOHMGzp4N24oqc/Fi6FkkLhKzOXYsMqln1eK11xL3k50a5eDBBLqdjoabSeIiMZuIpJ1Vi1OnIhEbc5wFMTmZQG/w8eOhupxcJPKZnEzU2IjJyWBtasepIRLXsJmYCNXl5CKRz4kToQeJqsnRowke0erUNIcOJepWDXWwiItEPgkatTN2AU4mKz7vxIjz5xO29OmJE6F9tYtEDtVEZTUNHMEHzTk1zZEjCUq4GBsLzcdWFpEQkW+KyEkR2ZVXt1xEnhKRfcHf9qBeROTLItIvIi+KyA15n7k72H+fiNxd6LsqxtmziUmbGB2F4eGwrXCc0khctvrJk6F8bbl6En8NbJ9Vdz/wQ1XdCvwweA/wAWBrUO4FvgYmKsCngHcANwGfyglLVUjQ1RZyRp3jlI3BwQT1JkJ6RpVFJFT1X4DZHsI7gIeD1w8Dd+bVf0uNZ4A2EVkN3AY8papDqnoWeIo3C0/lSEhW0/nzvvKXEx8SlZB45kwoiljJmESXqh4DCP6uDOrXAvmZzgNBXbH6NyEi94pIn4j0nSqHuk5PJ8b/krBhIE4COHYsIZlO09OhtPDCCFxLgTqdo/7NlaoPqWqvqvZ2dnaWbtHZs4nos46Pw1DSpjVwYs/4OAwnpXccwniJSorEicCNRPA3F3UZANbl7dcNDM5RX3kSMiHMqVN4RpMTS04kJXs9hGdVJUXicSCXoXQ38Fhe/UeDLKebgXOBO+pJ4FYRaQ8C1rcGdZUnASKRsAxfJ2GMjCRkTewQnlV15fgnIvK/gfcAK0RkAMtS+lPgURG5BzgM/Eaw+xPA7UA/MAZ8DEBVh0TkM8BzwX6fVtXqDJdJQDxi+JyvFeHEm1OnoLs7bCsqzPg4XLgATU1V+8qyiISqfqTIpvcV2FeB+4r8n28C3yyHTQvm4sVENEHOJCUDxEksZ84kQCTAgtdVFAkfcZ2AfNAEJW85CWZ8PCET/1X5meUikYDl2IaHE5G85TjJmM+pys8sF4nXXw/bgorjvQgnKSQgB6XqzywXiZiLhKqLhJMcxsctzBhrxsaq6hpwkTh/PmwLKsr5875mhJMsYh9mVDWhqBLJFomLF2P/BE1AyMVxLiER13wVG7fJFokqqnFYJCLbw3HySMQ17z2JKnHhQtgWVJxE3DCOk8fkZAKGPrlIVImYi8T4eGLWUXKcS4h5qLGq0flki0TMmxsJ8KY5TkFif+27SFSJmOfKxbyj5DhFGYv7tV/FBm6yRWJiImwLKoqLhJNULsb92neRqBIxF4mYe9Mcpyixv/anpqo2oC7ZIhHzubNjroGOUxTV2N/eVctKSbZIxDz1ZyLuN4njzEHsRaJKP9BFIqZMTeFLlTqJJvYiUaXZIiInEiKyXUT2iki/iNxfsS/KZq1PGlNirH+OsyBifw8k0d0kImngK8AHgG3AR0RkW0W+LOYLLMT85znOvMT+Hkho4PomoF9VX1PVCeAR4I6KfFOMexGQgBvEceZhOu73QJWeYVETibXAkbz3A0Fd+Ym5SMT85znO/MT9HkioSEiBujcdCRG5V0T6RKTv1KlTVTCr9nCRcBynHERNJAaAdXnvu4HB2Tup6kOq2quqvZ2dnYv7JimkR/Eh5j/PceYn7vdAlW7yqInEc8BWEdkoIvXAXcDjFfmmmD9FU1E7s45TZVLxvsWr9gyrq8q3LBBVnRKRTwBPAmngm6r6ckW+LOZP0XQ6bAscJ1xifw9U6RkWKZEAUNUngCcq/kWplClxTJ33sb9BHGceYn8PVOkHxrs5PR8xvorq6oi/T9Zx5qAuck3gMlOlH5hskYjxVSQCmfj+PMeZl/r6sC2oMC4SVSCTCduCihL7m8RxiiGxv72r9gOTLRIxf4o2NoZtgeOEQ0N9zBMY0+mqBa5dJGJMQ0PYFjhOOMS+gVTFZ1eyRSLmT9GmprAtcJxwiP21X0UVTLZIxLy5sWRJ2BY4TjjE/tqvYgM32SIR8+ZGY2Oss3wdpyjNzWFbUGGq+OxykYgxIgm4WRxnFul07J0ELhJVI/Z9UmhpCdsCx6kuLS0xz2yCqj67ki0SjY2xn8OptTVsCxynuiTimq+iiyDeT8j5SIA/pqUl9jroOJewbFnYFlQBF4kqEnN/TCqVkJaV42CDkGPvRW5qqmpGiovE0qVhW1Bx2tvDtsBxqkNbW9gWVIEqP7NcJBLQzG5rw2eEdRLB8uVhW1AFqvzMcpFIgEhkMtAa/w6Tk3AymUTczi4SVae5OdZThudYsSJsCxynsixfnoDUV6i6T60kkRCR3xCRl0UkKyK9s7Y9ICL9IrJXRG7Lq98e1PWLyP159RtFZIeI7BOR7wRrXFeHBDgy29t99LUTbzo7w7agCoQQmS+1J7EL+PfAv+RXisg24C7gGmA78FURSYtIGvgK8AFgG/CRYF+AzwFfUNWtwFngnhJtWzgJiOym09DREbYVjlMZmpsTkNUE1qCtcnepJJFQ1T2qurfApjuAR1R1XFUPAP3ATUHpV9XXVHUCeAS4Q0QEeC/w3eDzDwN3lmLbZZEAkQBY2RW2BY5TGVauDNuCKhFCZL5SMYm1wJG89wNBXbH6DmBYVadm1VeHhDgzlzQlJLDnJIpMJkG95BBEYt6IrYj8AFhVYNODqvpYsY8VqFMKi5LOsX8xm+4F7gXo6ekpttvCyWQs93hkpPT/FXFWr07Ez3QSRFdXQmYVSKVC8XrMKxKq+v5F/N8BYF3e+25gMHhdqP400CYidUFvIn//QjY9BDwE0NvbW1RMLosVKxLx9Fy2zPy358+HbYnjlE46nSBXU0jZJ5XS38eBu0SkQUQ2AluBZ4HngK1BJlM9Ftx+XFUVeBr4cPD5u4FivZTKkIjUCGNN9Rx5jlNRuroSkcFuhPSMKjUF9tdFZAC4Bfi/IvIkgKq+DDwK7Ab+CbhPVaeDXsIngCeBPcCjwb4AfwT8voj0YzGKb5Ri22XT0ZGQPiu0t8V+yionAaTTsKqQIzyuhCQSYo342qW3t1f7+vrK88927ICTJ8vzvyLOyAi88krYVsSbaakj/UYuhlNuenoSJBINDfCrv1rWBBsR2amqvfPtl4ym80LpSk6OaGsrtCUj89eJIQ0NCYpFgP3YkDIwXSTySUyzxFjfk4jMXyeG9KxPjHfYCPHZlKTDPD+NjYkZWAfWGluzJmwrHOfyaGuzuFpiSKdDTaxxkZjN6tVhW1BVVq9OwKLxTmxIp2H9+rCtqDJdXaFOvOYiMZuENa1TKdi0CV9vwqkJurutB5wo1oabs+4iMZumpgSN8TdaWmB1ssIxTg3S2pqwYDXYbBAh/2gXiUJ0d4dtQdVZuzYhs2g6NUk6DRs3JjDRYs2a0CP0LhKFWLMmcYsvpFKweXPo16PjFGTTpgS6mcAGg4SMPxIKUVeXuNgEmKdtw4awrXCcS+nqSlTS4QytrZFYEM1FohiJS6EwVqwI3QXqOG/Q0gLr1s2/XyyJyDPIRaIY7e02ZWoC6enxuZ2c8MlkYMuWhLpA6+oiExtN4uFfOBs3hm1BKKRSsHVrQn3ATiRIpeCKK6C+eivdR4vu7shMb+siMRdr1yb2SZnJ2E2asPi9EwXEkiiam8M2JEQ2bQrbgjdwkZiLVCrRkdymJhOKRHb3ndDYsD6hgeocq1ZFSiH99p+PjRsj0+0Lg6VLYfOWBOanO6HQ3e2JE2zZErYFl+AiMR+ZTGSyDMKivc2n7nAqz+rVicw8v5QVKyLXjXKRWAibNyfeOd/RAZtdKJwKsWpVglNd89m6NWwL3kSpy5d+XkReEZEXReQfRaQtb9sDItIvIntF5La8+u1BXb+I3J9Xv1FEdojIPhH5TrAGdjRoaEh0bCJHTijc9eSUkzVrIjGwOHw6OqwnETFK7Uk8BVyrqm8FXgUeABCRbcBdwDXAduCrIpIWkTTwFeADwDbgI8G+AJ8DvqCqW4GzwD0l2lZetmxJdGwiR0cHbNnqwWynPHR3R2Y4QPhcdVXYFhSkpFtdVf+f6huL+D4D5E73HcAjqjquqgeAfuCmoPSr6muqOgE8AtwhIgK8F/hu8PmHgTtLsa3s1Neb28mhvQ2uvDLxHjinFMRyQhIfg8jR1QXLl4dtRUHK2R78LeD7weu1wJG8bQNBXbH6DmA4T3By9QURkXtFpE9E+k6dOlUm8xfA5s2+Qk/A0qWwbVtih5E4JZBOwxVbQ11sLVqIwNVXh21FUeYVCRH5gYjsKlDuyNvnQWAK+HauqsC/0kXUF0RVH1LVXlXt7azmlZZOWxPaAWwcxbZtPoWHs3Dq6+15GIF566LD+vXW6ooo8zrZVfX9c20XkbuBDwLvU9Xcg30AyM9V6AYGg9eF6k8DbSJSF/Qm8vePFuvWwcGDcO5c2JZEgkzGXKmHDkE1O3VO7bF0qYX2MpmwLYkQmUzkG56lZjdtB/4I+JCqjuVtehy4S0QaRGQjsBV4FngO2BpkMtVjwe3HA3F5Gvhw8Pm7gcdKsa1iiMBb3hK2FZEilTL/8oYNnvnkFKZrlT0LXSBmcdVVkZ+gqtSYxF8AS4GnROTnIvJ1AFV9GXgU2A38E3Cfqk4HvYRPAE8Ce4BHg33BxOb3RaQfi1F8o0TbKkd7u+fsFWDlSo9TOJeSTlvvYX2PZ8S9iWXLamKgrsx4iGqT3t5e7evrq/4XT07C00/D+Hj1vzviTE+bR+7MmbAtCZdpqSP9Ri5G8mhpsVwPbzQUQATe/e5QlyMQkZ2q2jvffq7tiyWTgWuvDduKSJJO28PBB6onFLGxD1df7QJRlE2bama9Gh8dVgpr1sDgIBw7FrYlkaSjw4KVBw/C8HDY1jjVYMkSe/4tWRK2JRGmpSXywep8vCdRKm99qzeX5qC+3qYb37zZg5ZxRoLew7ZtLhBzIgLXX19TXWwXiVKprzehcOako8OSwjo78UkCY8ayZXZu16zx4PS8bN0auVle58NPaTlYtconAFwAdXWWKnuND8CLBQ0N9sy78kqfiGBBtLdHcpbX+fCYRLm45hoYGoKRkbAtiTzNzeaWOHMGBgY8QazWSKet19DV5T2HBZPJwA031OQBqz2Lo0oqBTfe6DPFXgY5F1TPeo9X1AKplC0MdN119rdGeh/hAAAQLUlEQVQGn3fhcf31NRus8SdaOWlpsTto586wLakZUilY1QWdK+DkSTh+3IagONFBBFZ2wZrVLuaLYssWc0nXKC4S5WbNGpvXqb8/bEtqinTaWqcrV9ocUMePw8RE2FYlm3TazseqVS4Oi6azM7LrRCwUF4lKcNVV8PrrcOJE2JbUHOm0PZRWrrQQz/HjMDY2/+ec8lFfb/GGzk73npZEc7O5oGt8QjO/BCqBiAWpfvITEwvnskmlbCXHFSvsEJ48aaJR47PIRJrWVhPn9vaaf66FTyYD73hHLLpgLhKVoq7OLpKf/AQuXgzbmppm6VIrPT1w+rS5o/yQlodMxhIIVq70NNaykUrB299uPYkY4CJRSZqaTCh++lOPxpaBTMbiFqtXw+iopdCeOQNTyZ1Db1GkUtZb6OiwgXDeaygzb3ubHdyY4CJRaVpbrVXxzDOQzYZtTWxoabHS02NDU4aGbH4o1+LCpFImCMuX26pwNTQrRG2RG3oeI1wkqkFHhwnFc8+5UJQZEXv4LVtm8YrRUROLc+c84F1fb8elrR2Wtfq4hopz9dWxnHnBRaJarFxpmQ59fR59rRAiM/GLdesshXZkZKbEPaU2nYalrdZ5XdZq3k6nSlxxhY2HiCEliYSIfAa4A8gCJ4HfVNVBERHgS8DtwFhQ/3zwmbuBPwn+xWdV9eGg/kbgr4Em4Angk1rrKyLNZtUqE4rnn/ceRRWor5/JkAKb/mN0dKaMjdW2Xjc1WWy0pQValkJTo8cXQuGKK2pq6u/LpaSV6USkVVVHgtf/Bdimqr8tIrcDv4uJxDuAL6nqO0RkOdAH9AIK7ARuVNWzIvIs8EngGUwkvqyq35/PhtBWpiuFEyesR+FCESrZrGVJnT8PYxfgwpi9L1ePo1wr09XVmSA0NtrMDrnicYUIcNVVNTlpHyx8ZbqSehI5gQhoxh78YL2LbwU9gWdEpE1EVgPvAZ5S1aHAyKeA7SLyI6BVVX8W1H8LuBOYVyRqkq4uy3p67jlPzQmRVGrmgZvP9DSMT8D4Ret9jE/A5ISJx+SknbLp6fLZkMlYqctAfcZ6QA0NVhobfUBbZHnLW2IZg5hNyZefiPwP4KPAOeBXguq1wJG83QaCurnqBwrUx5cVK+CXfgl27PBpUCNGOg1LmqwUQ9XEIps1wZietrps1lpK052r0aYl1B3aDwIpMUFIpez/54oHk2uQVMrSXGOWxVSMeUVCRH4AFJqd6kFVfUxVHwQeFJEHgE8An6LwsjK6iPpiNt0L3AvQ09Mz9w+IMsuWwbveZUIxOhq2Nc5lIFJkMG0qZdPG51qYpzotBhX3qHlSyGQsUzFG4yDmY952jKq+X1WvLVAem7Xr3wL/IXg9AKzL29YNDM5T312gvphND6lqr6r2dnZ2zvcTos2SJSYUueiqU7s0N9u5zHdBdHbCL/9yoh4qsSV3fhN2Lkvq7IpIfsTmQ8ArwevHgY+KcTNwTlWPAU8Ct4pIu4i0A7cCTwbbXheRm4PMqI8Cs0UovmQycPPNtmybU5t0d8O/+3fWO5xNYyPccotlwHj6UW3S2Qnvfncil1QsNSbxpyJyJZYCewj47aD+CSyzqR9Lgf0YgKoOBWmzzwX7fToXxAY+zkwK7PeJa9C6GCJw7bU2HPbFF8sXGXUqSyZja5zP558WsVTJzk544QVLqXJqgy1bLIspoQJfUgpsFKjJFNj5GBmxhYs8ThFtVq60RaYud2a86WnYvRsOHqyIWU6ZyGQsQN3VFbYlFaEqKbBOhWhtNdfFSy/BkSPz7+9Ul0zGgtPr1s2/byHSaUufXL3aeo3eq4geHR0mED5s3UUisqTTti5uV5c9SDw7JhqsWWNuwYaG0v/XihUW1H71Vdi/v7aHf8eFVMpiR5s3J9a9NBsXiaizerVN3fnii7ZMmxMOzc3W+i93Nl06bRPDrV1rPcehofk/41SGZcusYdbaGrYlkcJFohZoaLDc7MFB2LXLB99Vk7o6m3Zh06bKjnxrbYV3vhOOHrV4ha+qVD3SaUsq8N5DQVwkaok1a6wlu3evBT3dPVE5RCzmcNVV5XEtLZS1a20iyP5+c0F5lltl6eoy9+HsuVmcN3CRqDUyGbuoe3rg5ZdtPU+nvKxaZeKwdGk4359Om198/XqLVxw+7A2CctPSYskHK1eGbUnkcZGoVVpbbYDWiROwZw+8/nrYFtU+nZ32cG5vD9sSo7HRxmBs3gz79sHAgItFqTQ0mGupp8cnzlogLhK1TleXtYaOHrVWp6dTXj5dXRZ3iIo4zKa52QKqW7eaG2pgwKeZv1wyGRPbjRt9Wt3LxI9WHBCxaSHWrjWx2LfPB+LNRyplx2vz5vDcSpdLc7MN3rvySjhwAA4d8kW956OhwZIONmxwcVgkftTiRE4surstXfa11+DMmbCtihaNjebrX7++ugHpctLYaGmzV1xhvYoDB9zdOJuWFhOH7m5fnalEXCTiyqpVVkZGLBPq6NHkLnAkYi65nh5zLcUlzTGdnhG8oSELcA8OJjcjSsTO74YN5R/PkmBcJOJOa6sFP7dtg2PHbJqPpPQuli2bccPVaq9hoSxfbuXaa+08Hz1qmW9JCHQvXWrned26+J/nEHCRSAp1dXYTrVsHFy5Yi3NwEIaHw7asvLS12Sj11avNh5808s/z+Li5HY8ds4ZBnILdzc12jteu9RHSFcZFIok0NVnAdvNmG9l74oSVM2dqzyWVyZhrobPTXA3ekpyhoWHGHTU5aT2Lkyfh1ClrKNQSItZTWrnSznOtJBvEABeJpJMfyM1mrWdx5oyVs2ejJxqNjZaquny5zdTZ2hqfGEMlyWRmelhgqdL55zlqqdOplLkLly+3iRCXL/fspJDwo+7MkErN+La3BosOjo6acIyMWBkdrU4rNJUyl8LSpVaWLbNyuWs3OIVpbraSWyN+chLOnbMyMmLZUufPV6eRUF9v2UitrVaWLbO/PtgtErhIOHPT0vLmJRunp2FszMqFC+ayGh+36cwnJ61MT1vPJJu14KmIlVTKsnLq6qzU11tpaDA3WFOTzaPT2Og9hGqSyViLffZa6xcvXnqec+c6d56npmbOteqbz3XuPGcyM+e6sXHmXDc32zYnspRFJETkD4DPA52qejpYp/pL2BKmY8Bvqurzwb53A38SfPSzqvpwUH8jM8uXPgF8Umt92by4kk7PtPCdeNPY6L23hFNyf05E1gG/ChzOq/4AsDUo9wJfC/ZdDnwKeAdwE/ApEcnNhfC1YN/c57aXapvjOI5TGuVw+n0B+EMgv9V/B/AtNZ4B2kRkNXAb8JSqDqnqWeApYHuwrVVVfxb0Hr4F3FkG2xzHcZwSKEkkRORDwFFV/cWsTWuB/MWZB4K6ueoHCtQ7juM4ITJvTEJEfgCsKrDpQeCPgVsLfaxAnS6ivphN92KuKXpy2RmO4zhO2ZlXJFT1/YXqReQtwEbgFxanpht4XkRuwnoC6/J27wYGg/r3zKr/UVDfXWD/YjY9BDwE0Nvb68Ftx3GcCrFod5OqvqSqK1V1g6puwB70N6jqceBx4KNi3AycU9VjwJPArSLSHgSsbwWeDLa9LiI3B5lRHwUeK/G3OY7jOCVSqXEST2Dpr/1YCuzHAFR1SEQ+AzwX7PdpVR0KXn+cmRTY7wfFcRzHCRGp9aEIvb292tfXF7YZjuM4NYWI7FTV3vn283HvjuM4TlFqvichIqeAQyF89QrgdAjfWwq1aDO43dXG7a4uYdm9XlXnXZ2p5kUiLESkbyFdtShRizaD211t3O7qEnW73d3kOI7jFMVFwnEcxymKi8TieShsAxZBLdoMbne1cburS6Tt9piE4ziOUxTvSTiO4zhFcZEogoj8rojsFZGXReTP8uofEJH+YNttefXbg7p+Ebk/r36jiOwQkX0i8h0Rqa+C7X8gIioiK4L3IiJfDmx7UURuyNv37sC2fcGCULn6G0XkpeAzXw6mS6mUvZ8XkVcC2/5RRNrytkX+eBeimH0h2bJORJ4WkT3B9fzJoH65iDwVHKuncmu7LOZ6qbD9aRF5QUS+F7wveI5FpCF43x9s35D3PwpeRxW0uU1Evhtc13tE5JZaOd5vQlW9zCrArwA/ABqC9yuDv9uAXwAN2OSG+4F0UPYDm4D6YJ9twWceBe4KXn8d+HiFbV+HzZF1CFgR1N2OTXMiwM3AjqB+OfBa8Lc9eN0ebHsWuCX4zPeBD1TQ5luBuuD154DP1crxLvJ7itoX0vW8GptXDWAp8GpwbP8MuD+ovz/vuF/29VJh+38f+Fvge3OdY+B3gK8Hr+8CvjPXdVRhmx8G/nPwuh5oq5Xj/abfEtaFG+USXITvL1D/APBA3vsngwfpLdhEhZfsF5z003kPwEv2q5Dt3wWuAw4yIxJ/CXwkb5+9wYPjI8Bf5tX/ZVC3Gnglr/6S/Sps/68D366V413kNxS0r9p2zGHfY9hqknuB1UHdamDvYq6XCtvaDfwQeC/wvbnOce76CF7XBftJseuogja3AgcIYr6zj2OUj3eh4u6mwlwBvDvosv5YRN4e1F/uYkodwLCqTs2qrwgSj0WgfouZyR0jfbznoJh9oRO4YN4G7AC61GZgJvi7Mtjtco97JfkitvJlNng/1zl+w75g+7lg/2rbvQk4BfzPwE32VyLSTG0c7zdRqVlgI4/MvZhSHda9uxl4O/CoiGyi+OJIhcT2shdTWgjz2B3KIlALYS67VfWxYJ8HgSng27mPFbGjasd7kUTFjksQkRbg74HfU9WROcJMVbsu5kJEPgicVNWdIvKeeWyba1u1z0cdcAPwu6q6Q0S+hLmXihEVuwuSWJHQIospAYjIx4F/UOvjPSsiWWx+lWKLKVGk/jS2vndd0LKZczGlUuyWEBeBKsXuPPvvBj4IvC847sxhN0Xqy368F8lcdoeCiGQwgfi2qv5DUH1CRFar6jGxdeZPBvWXe71UincCHxKR24FGzI3zRYqf45zdAyJSBywDhqj++RgABlR1R/D+u5hIRP14F6ba/q1aKMBvY2tdgLmejmCqfg2XBsBew4KUdcHrjcwEKq8JPv93XBpk+50q/YaDzMQkfo1LA2PPBvXLMd9pe1AOAMuDbc8F++YC17dX0NbtwG6gc1Z9zRzvWXYXtS+k61mAbwFfnFX/eS4NpP7ZYq+XKvyG9zATuC54joH7uDRw/ehc11GF7f1X4Mrg9X8LjnXNHO9LfktYF26US3Bj/y9gF/A88N68bQ9i2RF7ycv4wTIUXg22PZhXvwnLFOoPLu6GKv2Gg8yIhABfCWx7CejN2++3Atv6gY/l1fcGv38/8BfMCsKV2dZ+TIh/HpSv19rxLvCbCtoXki3vwtwUL+Yd49sxf/0PgX3B31wD4bKvlyr8hvcwIxIFzzHW2/i7oP5ZYNN811EF7b0e6AuO+f8JHvI1c7zzi4+4dhzHcYri2U2O4zhOUVwkHMdxnKK4SDiO4zhFcZFwHMdxiuIi4TiO4xTFRcJxHMcpiouE4ziOUxQXCcdxHKco/x87gqqtCASKiQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "fig, ax = plt.subplots()\n", "# for maximally rotating black hole\n", "ax.fill(Xh1, Yh1, 'b', Xe1, Ye1, 'r', alpha=0.3)\n", "ax.fill(-1*Xh1, Yh1, 'b', -1*Xe1, Ye1, 'r', alpha=0.3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plot for rotating(normally) black hole" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[,\n", " ]" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAD8CAYAAACCRVh7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3XtwZGd55/Hvo27dNbqONKO5z3jGBNvYYAYDIReHqyEUZqvyh6mtjSvJFrUJSSWVTSU4/MFWsrW1JFULy24S4g3ZhV22gJBkoShY1hhIQhJsxsaM7x55PJ7R3CWN7tfufveP57TVM1ZLfY5afVH/PlWnpD59jvQe9dF5zvs+7/seCyEgIiKylqZqF0BERGqXgoSIiBSlICEiIkUpSIiISFEKEiIiUpSChIiIFKUgISIiRSlIiIhIUQoSIiJSVLraBdisnTt3hkOHDlW7GCIideWxxx4bCyEMbrRd3QeJQ4cOceLEiWoXQ0SkrpjZy6Vsp+YmEREpSkFCRESKUpAQEZGiFCRERKQoBQkRESlKQUJERIpSkBARkaIUJEREpKi6H0wnUjYhwMoKLC+TXVxhaXaFlYUMmcUMKwsZsksZMktZsstZcivRkskRsjly2UDI+UL03PgcRhP+vTWZL6kmmlKGpVOkmptoak6Ralld0m1pmjuaaW5P09LZTHNHMzQ3Q0sLNOmeTipPQUK2v2wWFhdhYQEWF8nMLrI4ucji1BLL04sszy6TmVtieW6FlRWPE7lcGX6tpUmFzKZ+htlqjEi3pWnubKFlRystO1pp7W6lrbeNlu42aGuD9nb/2ty8+cKLRBQkZHtYXIS5OV/m5wmzcyxOzLMwscDi1BKLi77J0pIHgXoRAiwv+8JsBsYywPx126RS0NrqMaK1Fdq7m2nra6d9oIOmrg7o6ICuLujs9I3MqnIsUp8UJKS+zM/DzMzqMjtLdmqW+ekMc/OwMO+bLCyUpzZQD7JZP+b5fOy4sAKsgE3T3uYxor3dY0RHVxPNvZ0eNHbsWF06O9WcJWtSkJDalMt5EJia8mV6GqanCSsZFhZgdtaXuTlYWISo6V8KBQ+WCwuFK3O0ts7Q2TlDV9dFuro8iDSlmzxw9PRAdzf09vrXtC4RjU5ngNSGuTmYnIRr13yZnoZcjlzO3yqoOJDNVruw9W1pyZeJCX9tBp2dObq7p9mxY5quLm/CAryW0dsLfX2+7Nih5qoGoyAhlReCB4HxcV+uXfOrVmR+frXyMDPTOM1G1RLCas0MAIOuTq9I9PTM0Dk1Q9O5c/5eOu3Bor/fl76+gogi25GChFTG9DSMjfkyPg6Z1V4/2exqq9LkZH0llrelgqBx4YLHgHwLVE9PhpbMVbh61bdtavJAMTAAg4O+kXIb24qChGyN5WW4csUvJlevXldTyL+db1mamXllaIHUoGx29bMCz3H39UFvH3S051ZrhC+84DWNgQEYGvKlo6O6hZdNU5CQ8pmehsuXfZmcfNWVf3nZ28EnJgqaNqTu5Hsaj476sIx8y1NHB15DzJ8D4MnwoSHYvds3Uj6j7ihISHIh+O3lxYtw6VJBH8xVmYwHhbExBYbtaHHRm6QuXPButgMDvrS2Rhvk261On/ZBfrt2wfCwBw41S9UFBQmJb2ICzp/34HBDMxJ47JichKtjMPXqCoVsUwsLXrsYHfUKxOCgVx5eyWuvrKxukE57wNizRwGjxilISGmmp/2f+8KFGzvev2Jx0dMPY2NKPje6fAXi5Zc9UAwOeu/ZV2QyfqNx/rzXMIaHYe9er4aoSaqmKEhIccvLHhjOnfMgsYZ8i9OVK0U3kQaWy612auvo8GCxc+cNvWZXVuDsWV/a22HfPti/3zPkUnUKEnK9ELw6cPasJx+LDFJYWfHNLl9WrUFKMz/vNYvRUQ8Uu3Z54vs6Cwtw6pQv/f1w8KDXMjQWo2oUJMQtLa3eza2RgM5bXPRUxPi4BrlJMtnsageo3j4Y3n1DU1RevivcU095zeLQIdUuqkBBotFNTXnPkwsX1r3qz856cLg2ieZJkrKZvOZLV5dXGHp710hJrKz4OXr6tLdXHTniX5W7qAgFiUYUgt/Gvfji6gQ+RczMeG5R+QbZSrOz3sLU0eHBouiQivzgzM5OuOkmz1+oKWpLKUg0klzOG4RHRnw01DoUHKQa5uf93uXCBe8dWzRYzM3ByZPw3HNw+LAvetjSllCQaATZrGcMX3zRkwrryI+knZqqUNlE1rCw4KfrxYuwdx/09RbZcHkZnn/eNz540GsXr4zkk3JQkNjOslk4c8b/gdYY9FZoacmDw/h4ZYomUor5eTj1gucsDhzwr2vKZPw8f+klT3AfPapgUSZlG+ZoZikz+5GZfT16fdjMHjGzU2b2JTNrida3Rq9HovcPFfyMB6L1z5vZe8pVtoaTy/k/y8MPwzPPrBsgMhnv0HTypAKE1K7ZWT+VR0Y2uN/J5TzB/fDD8Oyz6p9dBuUcC/+bwLMFrz8BfDKEcAy4BvxKtP5XgGshhKPAJ6PtMLNbgPuAW4F7gD81M2Wk4gjBqwPf/a53G1znvyk/HOLkSZ92SVNnSD2YmPBzdnR0g4dPZbMeUR5+2DPielJVYmUJEma2D/h54C+i1wa8HfhKtMnngA9G398bvSZ6/x3R9vcCXwwhLIUQXgJGgLvKUb6GMDYGf//38KMfrTvOAVbvyl566brHOojUhRA8sf3kkxt2zvOaxHPPwXe+4zMH6G4otnLlJD4F/C6QHxIzAEyGEPKXoFFgb/T9XuAcQAghY2ZT0fZ7gR8U/MzCfaSYuTm/4l+6tOGmmYzfgV25isY6SN1bXvbKQk+PpyHWTUEsLsITT/id0W23ebcpKcmmaxJm9n7gSgjhscLVa2waNnhvvX1u/J0fNrMTZnbiav4JWY0mm/U7pO99r6QAMTHhd15XrqAAIdvK1JSf2xcvllBRmJqCf/xHePzxDXv6iStHTeJtwAfM7H1AG9CN1yx6zSwd1Sb2ARei7UeB/cComaWBHmCiYH1e4T7XCSE8CDwIcPz48ca75F265DmHIrOxFlpZ8Q5O+aeKiWxHuZy3Jk1MwOEj0NG+wQ7nz/uA0te8xsdYaPR2UZuuSYQQHggh7AshHMITz98JIfxL4LvAL0Sb3Q98Nfr+a9Frove/E0II0fr7ot5Ph4FjwKObLd+2srgIP/yhLyUEiPFxv8NSgJBGMTcHTz/lOYsNaxWZDDz9NPzDP2hg0Dq2cpzE7wFfNLN/D/wI+Gy0/rPA/zSzEbwGcR9ACOFpM/sy8AyQAT4SQlCXhLyXX/bcQwmZ5kzGN1eXVmlE+U5+k5Mljq2bmvJAcdNNcPPNmubjBhbqPNt//PjxcOLEiWoXY+ssLHjCbWyspM1nZnxM0fLyFpdLNpS1NKmg7mPVlEr5QOydO0vcoasLXv966Ovb0nLVAjN7LIRwfKPtNOK6lp0757mHEmoPIXjibvQ8SkyLRLJZH1s3NeU9oDasJMzOemL76FGvVeixqgoSNWllBX78Y7/ql7j5iy9qMj6RYsbHPV9x9FgJSe0QfADe1atw550N/wwLhclaMzEBf/d3JQeI2VnPvSlAiKxvcRGeeTpGrm5y0v8Xz53b0nLVOtUkasnIiI99KDFPdPWqd2+t87SSSMXkcl7rnp2DA/tL6PmazXpOcHwcXve6hkxqK0jUgpUVn07j8uWSNg9h9RHUIhLf5UuwMO+ph3QpV8Fz57xm8aY3NVzzk5qbqm1mxrvflXjFz2bhhRcUIEQ2a3rae5WXPPB6ZsbnR2uwfz4FiWq6dAm+//0NnxKXt7zssx9r3I9IeSwueqCYnS1xh0wGHn3UE9sNQkGiWkZGfOR0idOwzi94gnqDCV5FJKZMxlOBsWYmeO45n/8pl9uyctUKBYlKC8G7tz777MbbRmZm4Nln9PwUka2Sy8GpEe8MUrLz5+Gf/3nbj1xVkKikTAYeecSzziWamvJH+OqZKSJbLPhM4iVMqrxqYsIH323jKr6CRKUsL/tdR4xblWuTnqRugBqtSM04e9YnCCxZfpT2Nh2spCBRCYuLfhJNTpa8y7VrMHJKYyBEqmF01FuTSra4CP/0T9tyymUFia02P+8BouTuE16DGBlRgBCppvPnY9YoVla8tWCbTb+sILGV5uf97iJGe+XUlGoQIrVidLTkGXJcNut5xxJnba4HChJbJR8gSng4UN7MjHe/VoAQqR3nzsXs9ZTN+liKbVKjUJDYCouLXu2MESDmF5SkFqlVL52JmW7I1yi2QY5CQaLc8r2YYjQxLS/DC+rmKlK7gucJZ2Zi7JMPFHXe60lBopzy1cwYSer8XEzbfDyOSN3LP2ai5LmewJPZP/hBXY+jUJAolxDgscdiVS9DgBdP1/X5I9JQMhm/qStxNh23tOQ1ijq9E1SQKJennoo9O+T58zBZ/02WIg1lcTFBF/XZWZ+rrQ6TjgoS5XDmjC8xTEzE7IMtIjVjetq7x8YyMQEnT25JebaSgsRmjY15LSKGhQWfI0ZE6tfFiwk6L507B6dPb0l5toqCxGYsLHgeIka9M5v1qqp6MonUv9OnPeUQyzPP1NVgOwWJpHI5OHEidjLq7NlYwydEpIblb/pipRpC8GdRxOomVT0KEkk9+2ysCfvAq6axRm6KSM2bm0uQX1xa8kBRB9MrKEgkcfly7HbFlRXlIUS2qwsXYw60A5+2ow4eg6ogEdfSEjzxROzdXjoTs2+1iNSP4PeNsXONL7xQ81N3KEjEdfJk7DzE+LjGQ4hsd0tLMBrnGRTgzU0/+lFN92TZdJAws/1m9l0ze9bMnjaz34zW95vZQ2Z2KvraF603M/u0mY2Y2Ukzu7PgZ90fbX/KzO7fbNnKbnQ05rMNvZnp5Ze3qDwiUlMuX441K4+bm4v1zPtKK0dNIgP82xDCa4G3AB8xs1uAjwIPhxCOAQ9HrwHeCxyLlg8DfwYeVICPA28G7gI+ng8sNWFpKfZ4CPC4omYmkQYRfFxt7Hz0Sy/5YLsatOkgEUK4GEJ4PPp+BngW2AvcC3wu2uxzwAej7+8FPh/cD4BeMxsG3gM8FEKYCCFcAx4C7tls+crm6ae9WhDD7Kx6M4k0mvl5uHIlwY4nT9bktB1lzUmY2SHgDcAjwK4QwkXwQAIMRZvtBc4V7DYarSu2vvrGxmI+8NapmUmkMY2Oxr6n9O5RNTgau2xBwsy6gL8GfiuEsN4E6rbGurDO+rV+14fN7ISZnbi61bfqISRqZhob86ZGEWk82ax3i43thRdqbpBdWYKEmTXjAeILIYS/iVZfjpqRiL7mK2CjwP6C3fcBF9ZZ/yohhAdDCMdDCMcHBwfLcQjFnTkTuwN0Lpdg8i8R2VYuX05wvc9m4bnntqQ8SZWjd5MBnwWeDSH8p4K3vgbkeyjdD3y1YP0vRr2c3gJMRc1R3wLebWZ9UcL63dG66slPHh/TlSt1O3W8iJRLSNRK7ZMA1tDT7MpRk3gb8K+At5vZE9HyPuA/Au8ys1PAu6LXAN8ATgMjwH8Dfg0ghDAB/CHww2j5g2hd9YyMxL7aZ7OaAlxE3Pi4P78+thrqEpve7A8IIXyftfMJAO9YY/sAfKTIz/pL4C83W6ayWF5ONI/GlSvq8ioiqy6ch6NHY+505Yp3ie3v35IyxaER18W8+GLsq30uF3usnYhscxMTCXPRzz9f9rIkoSCxlpWV2E+aA+/RFLvbm4hsexeT9HQaG6uJeZ0UJNZy5kyiNiPVIkRkLYlvIEdGyl6WuBQkbpTLJcpFTE3VXPdmEakRISScfeHSpaoPuFKQuNGFCwmeRwiXkwzDF5GGceVKwmcMVflBNAoSN0qQi1hejv2QOhFpMMvL3uIQ27lzVZ1KXEGi0PR0okTR2BhFJhAREVmVqMkpk0k4Kq88FCQKnTu38TZrGBsrczlEZFuanEw4jirhtakcFCTyQrIx9PPzSliLSGlCSNirdWLCLzZVoCCRNzaWKGFdo88JEZEaNT6ecMcqNTkpSOQlnHBJQUJE4pieSdjkVKVJ4RQkwOuACUbCLS6qqUlEYgoJezlNT1elyUlBAryRMMHc3ur2KiJJJL52VGFaBwUJSPhA2oR3AyLS8KamEg6sS3it2gwFCUjUeTmXi/3AOhERwHMSC0meMzE+7hefClKQyGQSVQnm5yv+WYnINjKd5CYzl6v4zLAKEhMTiep9qkWIyGbMJr2GVLhLpYJEwqg8O1vmcohIQ0k8uauCRIUl7GZQpcGPIrJNLC0lfMZEhbtVKkgkyEdkMokGZ4uIXCdR8np5OeGOyTR2kFhaSnS1r+DnIyLbWOIWienpspZjPY0dJBJmnzXKWkTKIfENZwWToo0dJBJmjhQkRKQcEjdbK0hUSMK6nvIRIlIOia8lFew5oyCRQIJpnkREXiXxtUSJ6wpJ2G6kICEi5RBCwm6wFWzzbuwgkbCul2gueBGRNawkuZ5ksxW7EDV2kEhQJcjlNGeTiJRPNum1vkJNGjUXJMzsHjN73sxGzOyjW/rLEkTibHYLyiEiDSuT9JrSiDUJM0sBfwK8F7gF+JCZ3bIlvyyXSzSxn4KEiJRTTkEilruAkRDC6RDCMvBF4N4t+U0Jr/a5JA8KEREpItHDh6Bid6y1FiT2AucKXo9G68ov6SejICEiZZQ4SCTeMZ5aCxK2xrpX/SXM7MNmdsLMTlxN8FQ5EREpTa0FiVFgf8HrfcCFGzcKITwYQjgeQjg+ODiY7DfZWvGolP2S7SYiUlZJr2Ex1VqQ+CFwzMwOm1kLcB/wtS35TU3JDr1JQUJEyijhpQhSqbKWo5h0RX5LiUIIGTP7deBbQAr4yxDC01vyyxL+gRN/oCIia2hKeq1vxCABEEL4BvCNivyy5ubYY+Ir9LmISINIJb3xTFfm8t3Y98XNzbF3SaUq1hQoIg0g8bW+paWs5SimsYNEa2ui3SoUwEWkASS4V/V270Q7JvhVFfkttaqtLdFuFQrgIrLdWcKbzoQ3uEk0dpBob0+0W7OChIiUQUtzwubrhNeuJBo7SHR0JNqtgkFcRLaxxNeShNeuJBo7SHR2JtqtTUFCRMogYYs3dHWVtRzraewgsWNHot0qWNMTkW1MQaLWtbcn6iGgICEi5ZD4WtLdXdZyrKexgwQk+mM3N1es95mIbGOJUgvptHISFdXbm2i3hOkMERHAbzQTdafv7a3oiF4Fib6+RLspSIjIZnQmTSskvGYlpSDR359ot4Q5bxERALqS3mgmvGYlpSDR2proit/ZqTmcRCS5RLlnMwWJqkjw4KJUSk1OIpJMKpUw99zXV/HJ4xQkAIaGEu3W01PmcohIQ9jRnfDZNAmvVZuhIAEwMJAoOifsGCUiDa436Q3mrl1lLUcpFCTAQ3qCCN3ZqfESIhJfohvMjo6KDqLLU5DIGx5OtFuFe6OJSJ3r6ko4PmLPnrKXpRQKEnm7diVqcqpwRwMRqXOJrxl795a1HKVSkMhLpRLVJnbsUJOTiJQuUZDYsaMqTU2gIHG9/ftj72IGO3duQVlEZNvp7k7Y1HTgQNnLUioFiUIDA4mm4FWQEJFSJBiS5R1r9u0re1lK/vVV+8216tCh2Lu0t2uaDhFZXyqVsKPLnj0Jqx/loSBxo/37EyWwqzDGRUTqyOBgwgF0hw+XvSxxKEjcKJ1O1P7X16cEtogUYQnHwfX3V33UroLEWo4ciR3ym5qqMhhSROpAX6/PJRrb0aNlL0tcChJraW9PlCgaGvJ2RxGRQonG6nZ318Sdp4JEMceOxZ4LPJ1WbkJErtfdnajTJLzmNWUvSxKbChJm9sdm9pyZnTSzvzWz3oL3HjCzETN73szeU7D+nmjdiJl9tGD9YTN7xMxOmdmXzKx66XzweVIS5CZ2706YnBKRbSnRQOmeHr+Y1IDNXs4eAm4LIdwOvAA8AGBmtwD3AbcC9wB/amYpM0sBfwK8F7gF+FC0LcAngE+GEI4B14Bf2WTZNu/mm2O3HzU318xnKyJV1t2dsHv8a19b9rIktakgEUL4fyGETPTyB0C+If9e4IshhKUQwkvACHBXtIyEEE6HEJaBLwL3mpkBbwe+Eu3/OeCDmylbWbS1eRI7pt27lZsQkYRj4IaGEo662xrlbBj5ZeCb0fd7gXMF741G64qtHwAmCwJOfn31HT3qwSKGdLpqEzaKSI3o70+QizCDW27ZeLsK2jBImNm3zeypNZZ7C7b5GJABvpBftcaPCgnWFyvTh83shJmduHr16kaHsDnpdKIPbdeuhF3eRKTumSWaCs4HztXY9A0bDi0OIbxzvffN7H7g/cA7Qgj5C/soUPgn2gdciL5fa/0Y0Gtm6ag2Ubj9WmV6EHgQ4Pjx40WDSdns3Qtnz8LYWMm7NDV53vvUqS0sl4jUpOHhBDeJbW0106Op0GZ7N90D/B7wgRDCfMFbXwPuM7NWMzsMHAMeBX4IHIt6MrXgye2vRcHlu8AvRPvfD3x1M2Uru9tvj91tqa9Pz8EWaTStrQnHRdx2W6IpgbbaZnMS/xXYATxkZk+Y2WcAQghPA18GngH+L/CREEI2qiX8OvAt4Fngy9G24MHmt81sBM9RfHaTZSuvzs5EUf7QIXWJFWkkBw4m6LgyPJz46ZhbzVZbiOrT8ePHw4kTJyrzy0KA738fJidj7XbxIpw7t/F2sr1kLU3qlb4Y0gj6+xPMpNHSAnffXfEkppk9FkI4vtF2useNwwze8IbYtwm7d3tFRES2r3QaDh5MsOPtt9d0LxcFibi6umL3djLz4RYxZ/kQkTpy6FCCmaAPHKjZZqY8BYkkDh2KPay6vR32V+8JhCKyhQYGEjy7uqvLk9U1TkEiqde/3ud3imH3LvV2EtluWlsTPNAylYI3vrEupmZQkEiquRmOH4/ddenIET2cSGS7MIObbkpwrb/9dp/YqQ4oSGxGTw/ccUesXZqb/aRac4y5iNSV/QcSTL1x+HDCSZ2qQ0Fis/btiz0JYHc37K+fc0RE1jAw4E3IsezcWXNzM21EQaIcbrkl9hOkhof9JBOR+tPRkSAP0dmZqIm62uqrtLXKDO68M3ZW+tCh2LlvEamydNofXBkrD9HSAm9+c10mJBUkyiWd9pMgxlU/lfLnGtXheSPSkMw8QMQa+5ZKwV131e2IWgWJcmpthbe8JdYZ1NLiU0LVQU84kYZ35EjMmbzNvImpr2/LyrTVFCTKrbPTA0WM6kFHh8/3ohHZIrXrwIGYecR8M/TQ0JaVqRIUJLZCd7c3PcWY9renJ+okpUAhUnOGhxM8u/6OO7bFIyoVJLZKX1/sQDEwAIeSTBAmIltmcDDBU+buuCPho+lqj4LEVurvj930NDTk1VoRqb6BgZhdXc18yp5t9E+sILHV+vrgrW/1DHWJdu/eNjchInVrYCDm7M35HMQ2++dVkKiEnh5429v8GbYlGh7eVjcjInUldoDId3PdBjmIGylIVEpXF/zUT8XqP7d7d4JRnSKyKYODMQNEc7M3K9d5L6ZiFCQqqb3daxQx+tENDWlCQJFKGR72+fdKDhAdHX7zF/thEvVDQaLS8ncdMdotBwbg5mN1N+WLSF3Zvz9mOqGvzwNE7Glg64suO9XQ1OQ9IF772pJvWXp74Sd+QlN4iJRb/pkQsZ4ium8f/ORP1vSzqctFQaKajh6FN72p5Ct//vHa7e1bXC6RBpFO+81XyS3AZv5P+IY3NEzVvjGOspbt2gU//dMlJ7RbW70CosegimxOW5tf70vuS9LS4k3FN920peWqNQoStaCz0wNFiU+rSqd99tjY0wSICODNt7fcEqNXen8//OzP+kODGkzpc0bI1kqlvAq7cyc8+SRks+tububjKDo74aWXIJerUDlF6tyePbB3b4weTEeP+lTNDdK8dCMFiVqzf7/3mnj8cZia2nDzgQFo74CRU7C4WIHyidSpVMpbinp7S9yhrc07mAwObmm5al1jhsZalx94V+L84R3tcOutehyqSDFdXXDbbTECxPCwNy81eIAA1SRqV1OTZ6h37YInnoC5uXU3z98ldXfD2bMbtlaJNIw9e3wpqbWoudmjSYn5wUZQlpqEmf2OmQUz2xm9NjP7tJmNmNlJM7uzYNv7zexUtNxfsP6NZvZktM+nzfQIHmA1YVbiPAGDg16r2Obje0Q21NLi3Vv37SsxQOzeDXffrQBxg00HCTPbD7wLOFuw+r3AsWj5MPBn0bb9wMeBNwN3AR83s/xz/f4s2ja/3z2bLdu2kUr5lf9tbyupv15bm1dCYiXnRLaRgQGvEHR3l7Bxayu88Y0+ZinGJJyNohw1iU8CvwuEgnX3Ap8P7gdAr5kNA+8BHgohTIQQrgEPAfdE73WHEP45hBCAzwMfLEPZtpe+PviZn/EIsMFDsc08SNx6q08vI9IImpvh2DFvei3peV8HD8LP/dy2nL21XDaVkzCzDwDnQwg/vqF1aC9wruD1aLRuvfWja6yXGzU1eUJ771546im4dGndzTs6vD/4pUtw/jyEsO7mInUr/wS5koJDTw+87nV+4yXr2vDPaWbfBtYatvUx4PeBd6+12xrrQoL1xcr0YbxpigON+tCF9navHl+96sFidrbopk1NfqPU3w9nzsD0dOWKKbLVWlt95taSmpZaWrwmvn+/2mJLtGGQCCG8c631ZvY64DCQr0XsAx43s7vwmkDhfIr7gAvR+rtvWP+9aP2+NbYvVqYHgQcBjh8/3tj3xoODnmx7+WV4/nlYXi66aVubJ/LGxuDcOVhZqVwxRcotf/Oze3cJiemmJu/8cfSoZsmMKXFzUwjhSeCVp2yY2RngeAhhzMy+Bvy6mX0RT1JPhRAumtm3gP9QkKx+N/BACGHCzGbM7C3AI8AvAv8ladkajpk/nWjfPhgZgdOn1+0Du3On17JHz8OVy2qCkvrT1+czDpQ0Ceu+fX53pJkxE9mqcRLfAN4HjADzwC8BRMHgD4EfRtv9QQhhIvr+V4H/AbQD34wWiSM/peWhQ3DqlA+YKDJfRyoFBw/A0KDXKiYnK1tUkSQ6Ojw4lNS0tGuX/z+UtLEUY6HObyOPHz8eTpw4Ue1i1Kb5eQ8W585tWF2YnvaYMj9fobI1gKylSYVMtYsQQyKeAAAK+klEQVSxLTQ3e4Vg584SUglDQz7XUsnDqxuTmT0WQji+0XYacb2ddXTAHXd4n8CREQ8WRWoW3d3er3x8HEZHYWmpwmUVWUMq5XmHoaENe317zeHYMfVYKjMFiUbQ0QG33+7zi7/4olcZMmvf4Q4MeC+osTHvMrtOHlxky6RSnpDevXuD4GDmUeToUTUrbREFiUbS1uaj626+2fvCnjmz5tSxZt5pamDAaxYXLqhmIZWRDw67dm0w3iGV8uTEkSMaLbrFFCQaUeGw1AsX/IEUa2Sum5o8WOzc6cHi4kVYWKhCeWXba272wLBr1wY1h44O75hx4IC6slaIgkQja2rybOC+fXDtmtcsLlx4Vd7CzAPFzp1wbRIuXYSZmeoUWbaXtjaflXtgYIOxDoODPmJuaEiD4CpMQUJcX58vt93mCe6zZ9eMBH29vszN+VQfExMaZyHxdXd7s1JPzzrX/NZWrzEcOKAmpSpSkJDrNTd7O++RI167OHvW25luGJ7d2emtVQcO+MwgV64oyS3rS6W8Njo0tM64tqYm3+DAAdUaaoSChBRXWLu4fNn7xl69el1zVHOzdy4ZHobJKbh6xb8Wn3lLGk1Xl7cW9fevk2/o6/Nmzz17fH4lqRkKErKxfGf1PXu8RnHhgi/j46+0NZmtNkUtL3sX2rExPXe7UTU3e55h56A/XndN3d0+m/GePWpOqmEKEhJPc7PPwX/woPeLvXx5NWBENYyWltWYMjvrb01MaELB7a6pySsEAwPr5Bp6e73aOTzsbZZS8xQkJLnCxGIm44mJS5f8axQRurp8OXDA8+Dj457qKDKWT+pMU5MHhP5+v/6/qjmpqckTEfnBD3ryW91RkJDySKdXqw8heCS4csWXqSnMvHWhu9u7uc/M+CbXrinhXW9SKQ8IfX0eIF4VGDo7PQkxNOQBYsP5NKSWKUhI+Zn5rWV/v8/CmU9SXL0K4+PY3NwrAePgQZ9UcHLSl9k5lPSuQW1tHhh6e71meN2YhtZWb2PKj7xUfmFbUZCQrVeYpAAftj0+/kqyooNZOjr87UwGpqZ8VtrpaU0HUi3ptAfxnh7/et1zG9raVif5GhiAHTuqVk7ZegoSUnnt7asjvcFrGlHbU/raNQbapxgY8JzG0tJqwJidVdDYKuk0dO2A7h0eFNrbo8RzPunQ1+fViP5+PbynwShISPW1tKxO3JM3NweTk7ROTTEYLayssLzs+YzZOZib9c004jsm826pnZ2rHQva2/HcQb4dsLfXg8OOHSU8G1S2MwUJqU2dnb7s3bu6bnGRlulpBmZmGIiqFrnpWRZmMszPe8CYn/fWrHWe3tpQzDwAdHRARyd0dkBHVxOp7k4PAPmlu9s30ghnuYGChNSPtjZfhl55tDpNQOfiIp2zswzOzXmkmJtjcWKehfF5FmYyLC7yyrJdu942Na3+edraoL2zifb+dtoGOmnq6lgNul1dCgYSi4KE1L/8lXHnztVV0dK3srJavVhcZGVmkcVrCyxOLbE8vcjyzBIr8yssLwWWl2s3iKRSPo6xpQVaOptp7myhtaeN1p422npaae1tj6JDuy+trQoEUhYKErK9NTd723pPj7+Mluv644TgyfPlZXJLKyzPLrM8t8Ly3AqZhRVWFjJklzJkl7OrX1dyhEyWXCZHyEZLLhByniAJuYClmrBsE9ZkYIY1GZZqoikVfW1OkWr2r02tadItKdJtadJtaZo7mkl3NNPS0UxrVzOp9pYoQrQoRyAVpSAhYuZ33q2tNO2Atp1eCxERb9IVERFZk4KEiIgUpSAhIiJFKUiIiEhRChIiIlKUgoSIiBSlICEiIkUpSIiISFEKEiIiUpSFOp9n2cyuAi+X8UfuBMbK+POqQcdQG3QMtUHHsLaDIYTBjTaq+yBRbmZ2IoRwvNrl2AwdQ23QMdQGHcPmqLlJRESKUpAQEZGiFCRe7cFqF6AMdAy1QcdQG3QMm6CchIiIFKWahIiIFNVwQcLMfsPMnjezp83sjwrWP2BmI9F77ylYf0+0bsTMPlqw/rCZPWJmp8zsS2bWUuHj+B0zC2a2M3ptZvbpqJwnzezOgm3vj8p5yszuL1j/RjN7Mtrn02aVed6lmf2xmT0XlfNvzay34L26+hyKKVbeajOz/Wb2XTN7Nvof+M1ofb+ZPRT9HR8ys75ofezzqoLHkjKzH5nZ16PXa54LZtYavR6J3j9U8DPWPN8qeAy9ZvaV6P/hWTN7a819FiGEhlmAnwO+DbRGr4eir7cAPwZagcPAi0AqWl4EjgAt0Ta3RPt8Gbgv+v4zwK9W8Dj2A9/Cx4fsjNa9D/gmYMBbgEei9f3A6ehrX/R9X/Teo8Bbo32+Cby3QuV/N5COvv8E8Il6/BzWOb6i5a32AgwDd0bf7wBeiP7ufwR8NFr/0YLPJPZ5VcFj+W3gfwNfX+9cAH4N+Ez0/X3Al9Y73yp8DJ8D/nX0fQvQW2ufRdVP2gp/IF8G3rnG+geABwpefyu6eL4V+NaN20Uf0ljBhe667SpwHF8B7gDOsBok/hz4UME2z0cXhA8Bf16w/s+jdcPAcwXrr9uugsfyL4Av1OPnsM4xrVneaperSFm/Crwrf75E64aB55OcVxUs9z7gYeDtwNfXOxfy51H0fTrazoqdbxU8hm7gJaLc8I1/41r5LBqtuelm4KejKuffmdmbovV7gXMF241G64qtHwAmQwiZG9ZvOTP7AHA+hPDjG96Kewx7o+9vXF9pv4zfHUEdfQ4bKFbemhI1u7wBeATYFUK4CBB9HYo2i/uZVMqngN8FctHr9c6FV8oavT8VbV/tYzgCXAX+e9Rs9hdm1kmNfRbpcv2gWmFm3wZ2r/HWx/Dj7cOram8CvmxmR/C7ihsF1s7ZhHW2L4sNjuH38eaaV+1WpExx15fFescQQvhqtM3HgAzwhfxuRcpUlc9hE2q1XK8wsy7gr4HfCiFMr5OOqsr5sx4zez9wJYTwmJndnV+9Tnlq7hgiaeBO4DdCCI+Y2X/Gm5eKqcpxbLsgEUJ4Z7H3zOxXgb8JXid71Mxy+Jwoo3g7f94+4EL0/Vrrx4BeM0tHdyaF22/ZMZjZ6/C20x9H/9T7gMfN7K51jmEUuPuG9d+L1u9bY/uyWO9zAE+0Ae8H3hF9HlBjn8MmrHccVWdmzXiA+EII4W+i1ZfNbDiEcNHMhoEr0fq451UlvA34gJm9D2jDm20+RfFzIX8Mo2aWBnqACar/OY0CoyGER6LXX8GDRG19FpVqf6uFBfg3wB9E39+MV9EMuJXrE1in8eRjOvr+MKsJyFuj/f+K65Nkv1aF4znDak7i57k+qfVotL4fb/fsi5aXgP7ovR9G2+YT1++rULnvAZ4BBm9YX5efwxrHV7S81V6iz/rzwKduWP/HXJ8s/aOk51WFj+duVhPXa54LwEe4PnH95fXOtwqX/x+A10Tf/7voc6ipz6LqJ22FP5AW4H8BTwGPA28veO9jeO+G5yno5YP3KHgheu9jBeuP4L2DRqKTs7UKx3OG1SBhwJ9E5XwSOF6w3S9H5RwBfqlg/fHob/Ei8F+5IYG2heUewQP0E9HymXr+HIoc45rlrfYC/BTeFHGy4O//PryN/mHgVPQ1fyMR+7yq8PHczWqQWPNcwGsbfxWtfxQ4stH5VsHyvx44EX0e/ye6yNfUZ6ER1yIiUlSj9W4SEZEYFCRERKQoBQkRESlKQUJERIpSkBARkaIUJEREpCgFCRERKUpBQkREivr/IsqNEra/8i0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "fig, ax = plt.subplots()\n", "ax.fill(Xh2, Yh2, 'b', Xe2, Ye2, 'r', alpha=0.3)\n", "ax.fill(-1*Xh2, Yh2, 'b', -1*Xe2, Ye2, 'r', alpha=0.3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - The inner body represents event horizon and outer one represents ergosphere. It can be concluded that with decrease in angular momentum, radius of event horizon increases, and that of ergosphere decreases." ] } ], "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.2" } }, "nbformat": 4, "nbformat_minor": 2 }