ex_17_semi_circles.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true,
    "execution": {
     "iopub.execute_input": "2024-12-20T01:45:41.289414Z",
     "iopub.status.busy": "2024-12-20T01:45:41.288644Z",
     "iopub.status.idle": "2024-12-20T01:45:41.528232Z",
     "shell.execute_reply": "2024-12-20T01:45:41.527859Z"
    },
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib as mpl\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Function to make a single semicircle pair"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2024-12-20T01:45:41.530434Z",
     "iopub.status.busy": "2024-12-20T01:45:41.530272Z",
     "iopub.status.idle": "2024-12-20T01:45:41.534190Z",
     "shell.execute_reply": "2024-12-20T01:45:41.533871Z"
    },
    "tags": []
   },
   "outputs": [],
   "source": [
    "def plot_semicircles(ax, x, y, i0, i1, scale=1):\n",
    "    HalfA = mpl.patches.Wedge((x, y), i0*scale, theta1=-90, theta2=90, color='C0', linewidth=0)\n",
    "    HalfB = mpl.patches.Wedge((x, y), i1*scale, theta1=90,theta2=-90, color='C1', linewidth=0)\n",
    "    ax.add_artist(HalfA)\n",
    "    ax.add_artist(HalfB)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Make a plot"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true,
    "execution": {
     "iopub.execute_input": "2024-12-20T01:45:41.535544Z",
     "iopub.status.busy": "2024-12-20T01:45:41.535439Z",
     "iopub.status.idle": "2024-12-20T01:45:41.690154Z",
     "shell.execute_reply": "2024-12-20T01:45:41.689826Z"
    },
    "jupyter": {
     "outputs_hidden": false
    },
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUMAAAE8CAYAAABaR9MTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy89olMNAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAgyUlEQVR4nO3deXwU9f3H8dfmWgiwS4gcIYSEU4lcEcNVFGpRoFAEFZVaRX8qQkGp2t+jxfprtR5RW7FaIi2oIFq1Uo1WgYLUE+WQQBAEDUIISxKOIGQTkmyu/f1BiaZcOWZ3Zjbv5+OxD9jJ7Pf7yWT3vd+Z+e6sw+/3+xERaebCzC5ARMQKFIYiIigMRUQAhaGICKAwFBEBFIYiIoDCUEQEUBiKiAAKQxERACLMLuBs7r77bkpKSnC5XGRlZTF79mwmT55sdlkiEoIsHYYOh4NFixYB8P7773PttdcqDEUkICy9mzxv3rza/2dnZ9O/f38TqxGRUGbpkSHAli1bePjhh/F4PGRkZJzy85qaGvbu3UtkZCQOh6N2udPpxOl0BrNUETGBz+fD5/PV3vf7/VRWVpKUlERYWAPGe36bWL16tb9bt27+kpKSOst3797tB3TTTTfd6tx2797doIyx7MiwurqasrIyWrduDcDll19OcXExmzZtYuTIkbXrRUZGArBx40bi4uJqlwdzZJiamsrnn38elL5Ox+v1kpCQgMfjweVymVaHtoO2wUnB3A7/PTIsKChg8ODBtdlQX5YNQ4/Hw3333ccrr7wCQH5+PsXFxSQlJdVZ7+SucVxcHF26dAl2mQCEh4eb+sQ7yeVymVqHtoO2wUlW2A7fP2xWH5YNw3bt2lFdXc0tt9xCTEwMO3bsYPHixSQmJppd2ilmzZpldgmWoO2gbXCSHbeDw++395Wu9+/fX7tbYNbI0Gxerxe3201RUZHp78Zm0nbQNoDGZ4Klp9bUx8njgs35zLHT6eR3v/tds94GoO0A2gbQ+Eyw/chQ74Qi8n2NzQTbjwxFRIygMBQRQWEoIgIoDEVEAIWhiAigMBQRARSGIiKAwlBEBFAYiogACkMREUBhKCICKAxFRACFoYgIoDAUEQEUhiIigMJQRARQGIqIAApDERFAYSgiAigMRUQAhaGICKAwFBEBFIYiIoDCUEQEUBiKiAAhFIapqakkJyeTnp5udikiYoL09HSSk5NJTU1t1OMdfr/fb3BNQeX1enG73RQVFeFyucwuR0RM1thMCJmRoYhIUygMRURQGIqIAApDERFAYSgiAigMRUQAhaGICKAwFBEBFIYiIoDCUEQEUBiKiAAKQxERQGEoIgJAhNkFiDG2eo5RVeMnMTaa81o7zS5HxHYUhiHiV298wVcHigGIc7egb7ybS3u3Z3JKPK2d+jOLnIteJSGooKicgqJy3ttxkMdW7GRSSjzTL+1OYmwrs0sTsSwdMwxxxyuq+duGfYz508c898keampsfS3fs9pzuITlXxSQ5TlGeWW12eWIzWhk2EyUV9bw8PKdrNhWQPoNFxHnbml2SYb7OPswD7yzA4CIMAfnd2rD5JR4pgxKwB0daXJ1YnUaGTYzm/cd45oF68g9ctzsUgKqqsbPl/leHl6+kyFpa3h0xU6NFuWsFIbNUN6xMn66aAP5x8rMLiUoyitrWPjxHsY9/QmZud+aXY5YlMKwmco7VsZdr26x9zHEI7vh3w/Biz+BJ3pw8ddPnnX1nMLjTF24gX9tPxCkAsVOFIbN2Kbcozy/NsfsMhquaD+8+lP48yD45I+Q8zGUFhLurzznQyuqa7jz1c188PWhIBQqdqIwbOb+uPprcgptdPxwx9vw7DD4ejnQuFFtZbWfe/6exeFin7G1ia0pDO1sfyb88054YRydq/Mb1YSvqobn1+4xuLAA2fYPWHYz+LxNbupoaSW/ydjW9JokZCgM7SgvE/56KTx3GWxeCvs+w0njRzkZm/Mo8VUZWGAAHM6Gt2eBv8awJlfvOMjGHJ1QkRMUhnaz5yNYMgEKthrW5PGKajK25BnWXkC8PQuqyg1vdum6vYa3KfakMLSTYx74+8+gstTwpj/7ptDwNg2zdy3s3xiQpld9eYBvj1cEpO2gqqmB7FWw+n6zK7Ety34CpaSkhDlz5hAVFUVUVBR79uzhySefpHfv3maXZp6NCw05XnY6X+wvCki7hshcErCmK6v9ZHmOctkFHQPWR8AdL4SXr4aCLIhoAVc8bHZFtmTZkWFhYSEVFRUsWLCAp59+mjFjxjB9+nSzyzJPZTlseTlgzecdK7PuCCl3XUCbt/QbQX28NfNEEHLi/PofV32tM+WNYNkwTEpKYunSpbX3u3fvTl6exY9rBZJnA5QF9mC/JT+RUvotePcHtItvDpUEtP2AOrIbdr1XZ9H8D77h2r+uo6xCHz9sCMvuJgM4HI7a/69cuZKZM2eecV2vt+7uo9PpxOkMoYuclh0NeBe+KuPO1BqmPPCjtvJKC/7e9XXwS0433zKn8DjZB4sZkNA26CUFm8/nw+f7biT831lQX5YdGX7fqlWrOHr0KL/4xS/OuE5CQgJut7v2lpaWFrwCgyEy8FeZiQx3nHulYAuPCngXlvy96yu63WkXR4WH0cndIsjFmCMtLa3Oaz8hIaFR7Vh6ZAiwevVqXn/9dZYsWUJY2Jmz2+Px4HK5au+H1KgQoGNfcISDP3C7Pu6WFrzMVZs4iGoDFcUB66Kjy8ah0XU4tE2EY7m1i7rEtOT+8X3s/Xs1wNy5c7nnnntq73u93kYFoqVHhu+++y4ZGRksWrSIiIgI5syZc8Z1XS5XnVvIhaE7Hs4fF7Dm2zgj6NouOmDtN1pYGHQeGNAu+ndxB7T9gAoLg5/8CcJPPN8dwNpfXcbYvnGmlhVMTqfzlNd/Y1g2DHNycrjqqqt488036dy5M506dWLhwoVml2WuwbcHrOkL4111jtFaSr9rAtq8rcMQoMdlcNsa6DcFomPNrsa2LLub3K1bNyoqLDrVwyzdR8Gw2bBuvuFNp3SNMbxNw/S7FtY8GJCz6Rd2dtGzQxvD2w26uP5w9XNmV2Frlh0Zyhlc8TBccq/hzV4zqIvhbRomKhrGPhaQpm8alhiQdsV+FIZ243DAj34Lsz6HITOhRdN38X7QM5Ye7VsbUFwADbjuxAjRQF1iWnLlwHhD2xT7UhjaVfveMO4x+OU3cPcOPGGNH9ndPLybgYUF0KRn4cLJhjTlcMAT1/SnRWS4Ie2J/SkM7S4iCtzxVDkaNy3miuSOXJ5sk8/lhkfCNYth/JMQ1bSR7C3DuzG8x3kGFSahQGHYjMVER/LI5H5ml9EwDgek3gZ3boYf3g8xSQ1u4sqBnbl/fB/jaxNbs+zZZAmsMAc8fnV/2rex6XzMNh1h5P+euJUchoPbyd5VAbvOfJFahwNu/UE37vtxH8LCLDqNSEyjMGyGwhzw2FX9ueLCTmaXYozW7aH1Dzl6IAfYcdpVuraL5olr+jO0u+bhyekpDJuZqIgw/nBN/2ZzFrV/Fzc/G5LIxIGddbJEzkph2IwMSGjLH6/pT6+OITDJ+DS6xERzRXJH4mNa0i/ezYCEttafMiSWoTBsBmJbRXHHyO7cOqI74SF8rGx0ckdG2+XMuFiOwjBEhTlgUGIMPxuayLi+cURFaOKAyNkoDEPEuL5xDO0eS9d20fTr4iY5zkUrp/68IvWlV0uImDO6l9kliNia9p1ERFAYiogACkMREUBhKCICKAxFRACFoYgIoDAUEQEUhiIigMJQRARQGIqIAApDERFAYSgiAigMRUQAhaGICBBCYZiamkpycjLp6elmlyIiJkhPTyc5OZnU1NRGPd7h9/v9BtcUVF6vF7fbTVFRES6Xy+xyRMRkjc2EkBkZiog0hcJQRASFoYgIoDAUEQEUhiIigMJQRARQGIqIAApDERFAYSgiAigMRUKOt7wSX1W12WXYToTZBYiIccorq+n/wGoAYltF0TfeTb94N33j3QzvGYurRaTJFVqXwlAkRB05XsFH2Yf5KPswAC0jw5k4oDM3Dkukb7zb5OqsR7vJIs1EWWU1f9/kYcKf1zIp/VOyPMfMLslSFIYizVCW5xhXL/iMtBU7Ka/U8UXQbnLIuPvvWRQUlZEU24q+/zlG1C/eTXiYw+zSxKKqa/z89eM9rNl5kD9PvYjkzs37EngKwxCxs8DLVweKWb/nW/jcA0CcuwXXp3Zl6uAEOrhamFyhWNXuw8e5buE6np+WyuBu7cwuxzTaTQ5hBUXlPLUmmx88/r52h+SsisurmPbCRjbt/dbsUkyjMGwGKqtP7A79+JlPyMw9anY5YlFlldXcsuRzvjlUYnYpplAYNiN7Dh/n+oXr+OfWfLNLEYsqLq/i3mVbqa6x9beBNIqOGdpZyWEoyIKje3HVdKI+722V1X5+8doWACYO6BzY+oLso+zDvLM1n54dWtMv3k1qUjuiIkL8/b7KB18th33roWArUeVFwP81qcmtnmMs/HgPM0f1MKZGm1AY2k2VD754HTY9D/lbahfHtnkGOK9eTdT44d7Xs+jaLpqBCW0DU2cw+Epgx1vg2QAHv6R1y2H8Y/vw2h+f1zqKay9O4I5Le+CODrFPXlSWw9p58PlzUHqkdrEjwpgTZU+tyWZ0nw706tjGkPbsQGFoJ57P4e1ZUPh1k5uqrPbzy2VbWX7XCJwR4QYUF0QVpfBhGmxaDBXFtYujE3rVWa2wpIJnP9zNPzL388jkflye3DHYlQZG/hZ4czoUZgesi4qqGv78/jc8MzUlYH1YTYjvQ4SQzS/BC2MMCcKTvjlUwrzVgXtBBUTBF/CXH8Bnz9QJwrM5VOzj9qWbeO6TPQEuLghyPobF4wMahCf9a/sBDhf7At6PVSgM7WD7G/DPO8Fv/NSYxZ/t5UiJTZ7w+VtgyQT4tnGh9vDynby+yWNwUUF0OBtenQqVx4PSXUV1Da9t3BeUvqxAYWh1RfvhnV8AgTm7V1FVw2uf2yAgyovg1Z+Cr6hJzTz4zy/xfFtqUFFBVFMNb82AiuBOe3lVYSiWsfJX4PMGtItXNuyjxupTKVb9BoqbPiXoeEU1v393hwEFBdmWlyEvM+jd5heV2/PNoxEUhlZ2dC98vSLg3eQdK+ObwxaeaOvNh62vGtbcv3cetN8LfOMi07relte00bhdNCgMn3rqqUDVcYrKykqeeOIJWrVqxfbt24PWr6VkvQr+mqB0tW2/hZ/wW/4GNVWGNVfjhzc35xnWXsAd/BIObjOte4Xhabz66qtUV5/+IH56erohBZ20cOFCLrnkEkpLbfYObqT9G4PWlaWf8LmfGt7k1v3HDG8zYPZvMrX7HfmBPUxjFQ0Kw/j4eH7/+9+fsnzHjh3MmzfPsKIAZs2axbBhwwxt03YOBG9EnHesLGh9NdgB40dFOwts9AI/ZO4xTm95pan9B0uDwjAjI4Pu3buzYsV3x7Heeecdhg0bxpgxYwwvriG8Xm+dm89nk+kiZ+Or3zw6Q7qqCs7ueKME4ARSic+43e6A85l7PNdXaeHnBuDz+U55/TdGg8LwX//6F9OmTWP16tXk5eXxyCOPMGXKFB5//HGeffbZRhVglISEBNxud+0tLS3N1HoMEREVtK4irXwR2HCn4U1Ghdvo3GG4uR8ljAy38HMDSEtLq/PaT0hIaFQ7Dfo43rx583C5XPz617/myiuvJDc3lzVr1jBixAgqKiqIigrei/e/eTweXK7vrtTrdBr/Agq62F6QF5zjRTGtzPvbnVP73oZPK+nevpWh7QVUbE9Tu28ZZe2Pa86dO5d77rmn9r7X621UIDYoDLOzs7n55pvZu3cv0dHRDBgwgA0bNlBVVcXSpUt54YUXGlyAUVwuV50wDAmdU4IWhn2tfMn3+EGGh6Gtvh2u80BTu+/ZobWp/Z+L0+k0ZPDToH2FBx98kOzsbAoLC3nrrbeYOHEimzdvZsaMGSxdurTJxXzf2rVrmT17NgCPPvooy5YtM7R9W7hwctC66telbdD6arD+1xve5Li+cYa3GTBdBkOr9qZ1389ObxxN0KCR4bRp04ATo7BRo0YxatSo2p/deeedhhY2YsQIRowYwfz58w1t11aSfgAdLoRDXwa0m1ZR4Vxo5ZFhl0GQMOTEpboMcEGnNvb6ro+IKLjoJvjkSVO67xff1pR+g82wo8gPPfSQUU3J9417DAjsAewrU+JpEWnt40JM+BOEG3Nc874f9zGknaAaOsuU0WF0VDi9O1p7N9kohoVh27ZtjWpKvq/bpTDkjoB2cdOwxIC2b4iOyTC26TMEpg7uyqW9zdvlbLRWsTDe2Lm89TFxQGci7HTmvQmax29pd2Mehb5XB6Tpkb3bc0EnC+8if1/qbTDuD+Bo3Ch21PnteXDihQYXFUTJE2Hkr4Pa5Y12eKM0iMLQDsLC4arnYNR9EGbcnLM2zggevaqfYe0FxZDpcNuaE8dSG+CGIV1ZeOPF9v9OlB/OhSseNvR5cCaDEmO4sHPzOHkCCkP7CAuDUb+C6R9Cn4kQ1vRvbPjN+D7Et23Z9NqCLf4imPkp3PAGJF8Jbc58Zjg5zsXfbhvCI5P72T8ITxp+J0z/ADpfFNBupl/aPaDtW42+A8VuOvWF614CbwHsfv8/346XS+mBhoXa1MEJXD+4a2BqDAaHA3qNPnEDKP2WLzPzOK/gAD3an/h2vHH94hiUGGNunYHSqd+JQNy3HjJfhH3r4GiOYc2P7xfHmAs7GdaeHTj8fr/Fr+p5dl6vF7fbTVFRUehNum6AsX/6mK8O1O+zzFMGdeHxq/sTZuWP4EnDlR3DV1zI+fOa9j0557WOYvXdI2ln5U8lnUVjM0Ejw2YkMtzBnZf14s7LeuJwKAhDTsu2+CPaAE0Lw4cn9bNtEDaFwrCZSI5z8ccpA0i28uRqMd3ccRcwtm/z2j0+SWEY4lK6tuVnQxK5cmDzmS8mjXPv5b25Y2QPs8swjcIwxHRytaBvvIu+8W5G9+lorwsSiCkiwx38dkIyNw5LMrsUUykMQ8TrM4YRGRZm+cstibVc2NnFH67R4RNQGIYMVwtzLwAq9hIVHsbsy3ry81E9dPjkPxSGIs1IZ3cLpg7uynWDE+jQpoXZ5ViKwlAkhEWGO+jdsQ394t388IIOjO7TkXDNLz0thaFICIkMD+OZqSk4I8Lo7G7J+Z3ahM7HEANMYSgSQsLDHEwc0NnsMmxJbxkiIigMRUQAhaGICKAwFBEBFIYiIoDCUEQEUBiKiAAKQxERQGEoIgIoDEVEgBAKw9TUVJKTk0lPTze7FBExQXp6OsnJyaSmpjbq8fp2PBEJKY3NhJAZGYqINIXCUEQEhaGICKAwFBEBFIYiIoDCUEQEUBiKiAAKQxERQGEoIgIoDEVEAIWhiAigMBQRARSGIiKAwlBEBFAYiogACkMREUBhKCICQITZBYgY5dNvClmxrYCu7aLpG++mXxc3rhaRZpclNqEwlJCx62Axf9uwr/Z+VEQYP+7biRuHJTEoMcbEysQOtJssIauiqoa3svK5esFnzH5lM98erzC7JLEwhaE0C+9+UcAVT33ExpxvzS5FLEphKM1GYUkFN72wgc92F5pdiliQwlDsL/czWPsUXQ5/fM5VyytrmL40k31HSoNQmNiJwlDsqaoCNi6C9KGweByseYD4I5/V66Elvir+9x9bsflXhovBFIZiP2XH4KVJsOKXcHhno5rYkPMtyzL3G1qW6Q7thB1vm12FbWlqjdhLZTm8ch141je5qcWf7uXaixMMKMpku9bAx0+AZwP+iBb0f70FAxPa8rOhiVzepyNhYQ6zK7QFjQxDxO/e3s4z/97FIW+52aUE1rr5hgQhwM4CL5m5Rw1pyzSZL8IrU8CzoXZRcXkVn+wq5I6XMrlt6SbKK6tNLNA+LB2G+/btY9KkScyYMYMJEyawfft2s0uypoNf8sWe/cx7L5vLnvwodM+W1lTDpsWGNvnZNzbeVgd3wLt3g7/mjKu8/9Uhfvu2Xjf1YekwnDlzJlOnTuUvf/kL999/PzfccIPZJVlLRSm8dBUsGE6nmoPAiZMDP//bZsoqQnA08PVK8Bp7nO+LvCJD2wuqjQvBf+6/85ub8zhQFOJ7DAawbBgeOXKElStXMn78eACGDh1Kfn4+WVlZ5hZmJRsWwO5/n7L4WGkl73yRb0JBAbZvneFN7j5UYnibQbNtWb1Wq6rx81ZWXoCLsT/LnkDJzc0lOjqa1q1b1y7r0KEDOTk5DBw48JT1vV5vnftOpxOn0xnoMs2Ve+apJLZ+kZ9JufGjONseT6sohYr6/41DeWTo8/nw+Xy19/87C+rLsiPDhs4BS0hIwO12197S0tICVJmFnNf7jD9q3yYE3wgiow1vMiLcsi+Bs4toAY7weq/epoVlxz1NlpaWVue1n5DQuBkCln0mJCUlUVpaSknJd+9+hw4dIikp6bTrezweioqKam9z584NUqUmGjYbXPGnLI4Md3DlwFOX217HZMObTGjX0vA2gyIsDHqPrffq4/rGBbAYc82dO7fOa9/j8TSqHcuGYWxsLGPHjmX58uUArF+/nri4OFJSUk67vsvlqnML+V1kAHc8/HwdTHiKo462tYvv+3Gf0BwZ9r0GnC5Dm+wX39bQ9oJq8G31Wm1It3YkdzZ2u1mJ0+k85fXfGJYeOy9YsIC77rqLDz74AI/Hw8svv2x2SdbTwg0X/w/Vmz5jZO8Ibh3RjUt7tze7qsBwtoYB1584i2qQi+18ncMel8GgmyFzyRlX6doumqevP/0AQuqydBgmJiby9tv6eFF9LJsx3OwSgmPEPbDzXShu+tnyTq4WjDrf5m8cE/4ErTvBunSoKK5dHBMdybUXJ3DHyB60axVlXn02YukwFDmFKw5ufBNemgzFBU1q6qdDutr3BMpJDgf8cC4MvxO2/wMKv+G9gZfSNTYaZ0T9T7CIwlDsqEMfuO3fsOZ3Jy5MUN3wK1h3bRfNbZd0C0BxJnG2hkE34wB6mV2LTSkMxZ7c8XD1czD2MdjyEng2UkRivR7qcMAfrulPdJSe/vIdPRvE3lqdByPuBuCrT3Ng645zPuT+8ckM6R4b6MrEZhSG0myEOU4E4f+MCKHdYzGMwlCahe7nteIPU/ozKLGd2aWIRSkMJaQlxkZzw5Cu3DQsiRaROrsqZ6YwlJDRNjqKCzq1oWu7aPrGu7k4MYZhPWJxOHSlZzk3haGEjEkp8UxKCcHPZEtQ2HzGqYiIMRSGIiIoDEVEAIWhiAigMBQRARSGIiKAwlBEBFAYiogACkMREUBhKCICKAxFRACFoYgIoDAUEQEUhiIigMJQRARQGIqIAApDERFAYSgiAoRQGKamppKcnEx6errZpYiICdLT00lOTiY1NbVRj3f4/X6/wTUFldfrxe12U1RUhMvlMrscETFZYzMhZEaGIiJNoTAUEUFhKCICKAxFRACFoYgIoDAUEQEUhiIigMJQRARQGIqIAApDERFAYSgiAigMRUQAhaGICKAwFBEBFIYiIoDCUEQEUBiKiAAKQxERQGEoIgIoDEVEAIWhiAigMBQRARSGIiKAwlBEBFAYiogACkMRESAEwtDn89X5tzny+Xw88MADzXobgLYDaBtA4zPBsmH41ltvkZSUxPz588+6nsLwxO/+4IMPNuttANoOoG0AIRaGa9asoaioiKSkJLNLqZf09HSzS7AEbQdtg5PsuB0sGYajR49m2rRpZpdRb3b8wweCtoO2wUl23A4RZhfQVH6/H4CCgoI6y51OJ06nMyg1VFdX4/V6g9LX6Zzs28waQNsBtA1OCuZ28Pl8dXaJT2bByWyoL1PCcMyYMXz99den/dnatWvp0qVLvduqrKwEYPDgwYbU1lhut9vU/gESEhLMLkHbAW2Dk8zeDiezob5MCcNVq1YZ1lZSUhK7d+8mMjISh8NRuzyYI0MRMc9/jwz9fj+VlZUNPudg+93ksLAwunfvbnYZImJzljyBcvDgQWbPns2uXbtYtmwZjz76qNkliUiIc/gbepRRRCQEWXJk2Fj1nagdSvbt28ekSZOYMWMGEyZMYPv27WaXFHSVlZU88cQTtGrVqln+/gAlJSXceuutzJw5kzlz5vCTn/yE7Oxss8syxd13383tt9/Ovffey49+9CMyMjLq9TjbHzM8yW4TtY0yc+ZMbrrpJq677jrWr1/PDTfcwNatW80uK6gWLlzIJZdcQmlpqdmlmKawsJCKigqef/55AObPn8/06dP58MMPzS3MBA6Hg0WLFgHw/vvvc+211zJ58uRzPi5kRoZ2m6hthCNHjrBy5UrGjx8PwNChQ8nPzycrK8vcwoJs1qxZDBs2zOwyTJWUlMTSpUtr73fv3p28vDwTKzLPvHnzav+fnZ1N//796/W4kAnD5ig3N5fo6Ghat25du6xDhw7k5OSYWJWY5ftTy1auXMnMmTNNrMZcW7Zs4eqrr+aFF17gpZdeqtdjbLObbORE7VChc19yOqtWreLo0aM8/fTTZpdimpSUFN544w3ee+89LrnkErZt20arVq3O+hjbhKGRE7VDRVJSEqWlpZSUlNSODg8dOtTsjpvKd1avXs3rr7/OkiVLCAtrfjt+1dXVlJWV1b4eLr/8coqLi9m0aRMjR44862Ob39YKIbGxsYwdO5bly5cDsH79euLi4khJSTG5MjHDu+++S0ZGBosWLSIiIoI5c+aYXVLQeTwepk+fXns/Pz+f4uLieg0QQmae4cGDB3nooYfIyMigZ8+ejBkzhvvuu8/ssgIuNzeXu+66i7i4ODweD2lpafU+YBwq1q5dy2uvvUZ6ejpTp05l8uTJTJkyxeyygionJ4fzzz+fmJiY2mOHRUVFlJWVmVxZcHm9Xm6//Xaio6OJiYlhx44dTJs2jalTp57zsSEThiIiTaHdZBERFIYiIoDCUEQEUBiKiAAKQxERQGEoIgIoDEVEAIWhiAigMBSbW7JkCYmJiYwePZoVK1YA1H4OtWfPnjz++OMmVyh2oU+giO316tWLF198keHDh9cuu+uuu2jZsqXCUOpNI0OxtaKiIvbt23fKxSk2btzIRRddZFJVYke2uYSXyOls3ryZ6OhobrnlltplNTU1ZGZmMmjQIBMrE7tRGIqtZWZmMmnSJBYvXly7bNu2baxatYoePXqYWJnYjXaTxdY2b97MwIED6yzbunUrKSkpdS6DL3IuCkOxtczMTAYMGFBnWVZWlo4XSoMpDMWWjh49yuzZs9m1axdLliyhvLwcgMcee4yMjAw2b97Mpk2bTK5S7ERTa0RE0MhQRARQGIqIAApDERFAYSgiAigMRUQAhaGICKAwFBEBFIYiIoDCUEQEUBiKiAAKQxERAP4fjGRoHImOpPkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 337.5x337.5 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(figsize=(3 + 3/8, 3 + 3/8))\n",
    "\n",
    "X, Y = np.meshgrid(np.arange(3), np.arange(3))\n",
    "I0 = np.random.rand(9) + .1\n",
    "I1 = np.random.rand(9) + .1\n",
    "\n",
    "for x, y, i0, i1 in zip(X.ravel(), Y.ravel(), I0, I1):\n",
    "    plot_semicircles(ax, x, y, i0, i1, scale=.3)\n",
    "\n",
    "\n",
    "ax.axis([-1, 3, -1, 3])\n",
    "\n",
    "ax.xaxis.set_major_locator(mpl.ticker.MultipleLocator(1))\n",
    "ax.yaxis.set_major_locator(mpl.ticker.MultipleLocator(1))\n",
    "\n",
    "ax.xaxis.set_minor_locator(mpl.ticker.AutoMinorLocator(2))\n",
    "ax.yaxis.set_minor_locator(mpl.ticker.AutoMinorLocator(2))\n",
    "\n",
    "ax.set_xlabel('$H$')\n",
    "ax.set_ylabel('$K$')\n",
    "\n",
    "fig.savefig('ex_17_semi_circles.pdf', dpi=1200)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}