{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib widget"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pysketcher import *"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "from ipywidgets import FloatSlider, AppLayout, Label, HBox"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "from IPython.display import SVG, display"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "import sys\n",
    "from ruamel.yaml import YAML\n",
    "from math import tan, radians, sin, cos"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "myfig={}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "yaml = YAML()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "gw = \"\"\"\\\n",
    "libraries: [\"from math import tan, radians, sin, cos\",\"from pysketcher import *\"]\n",
    "fontsize: 18\n",
    "g: 9.81              # constant gravity\n",
    "theta: 30.0          # inclined plane angle\n",
    "L: 10.0              # scene length\n",
    "a: 1.0               #\n",
    "xmin: 0.0\n",
    "ymin: -3.0\n",
    "r: 1.0\n",
    "setframe:\n",
    "    action: \"drawing_tool.set_coordinate_system(xmin=xmin, xmax=xmin+1.5*L,ymin=ymin, ymax=ymin+L+1,instruction_file='tmp_mpl.py')\"\n",
    "setblackline:\n",
    "    action: \"drawing_tool.set_linecolor('black')\"\n",
    "B: point(a+L,0)\n",
    "A: point(a,tan(radians(theta))*L)\n",
    "normal_vec: point(sin(radians(theta)),cos(radians(theta)))\n",
    "tangent_vec: point(cos(radians(theta)),-sin(radians(theta)))\n",
    "help_line: Line(A,B)\n",
    "x: a + 3*L/10.\n",
    "y: help_line(x=x)    \n",
    "contact: point(x, y)    \n",
    "c: contact + r*normal_vec\n",
    "outer_wheel: \n",
    "    formula: Circle(c, r)\n",
    "    style:\n",
    "        linecolor: blue\n",
    "        filled_curves: blue\n",
    "hole: \n",
    "    formula: Circle(c, r/2.)\n",
    "    style:\n",
    "        linecolor: blue\n",
    "        filled_curves: white\n",
    "N: \n",
    "    formula: Force(contact - 2*r*normal_vec, contact, r'$N$', text_pos='start')\n",
    "    style:\n",
    "        linecolor: black\n",
    "mg: \n",
    "    formula: Gravity(c, 3*r, text='$Mg$')\n",
    "    style:\n",
    "        linecolor: black\n",
    "wheel: \n",
    "    formula: \"Composition({'outer': outer_wheel, 'inner': hole})\"   \n",
    "    style:\n",
    "        shadow: 4\n",
    "body: \n",
    "    formula: \"Composition({'wheel': wheel, 'N': N, 'mg': mg})\"\n",
    "    style:\n",
    "        linecolor: black\n",
    "wall: \n",
    "    formula: Wall(x=[A[0], B[0]], y=[A[1], B[1]], thickness=-0.25,transparent=False)\n",
    "    style:\n",
    "        linecolor: black    \n",
    "angle: \n",
    "    formula: \"Arc_wText(r'$<bslash>theta$', center=B, radius=3, start_angle=180-theta, arc_angle=theta, fontsize=fontsize)\"\n",
    "    style:\n",
    "        linecolor: black\n",
    "        linewidth: 1\n",
    "ground: \n",
    "     formula: Line((B[0]-L/10., 0), (B[0]-L/2.,0))\n",
    "     stlye:\n",
    "         linecolor: black\n",
    "         linestyle: dashed\n",
    "         linewidth: 1\n",
    "x_const: \n",
    "    formula: Line(contact, contact + point(0,4))\n",
    "    style:\n",
    "        linestyle: dotted\n",
    "    transform: rotate(-theta, contact)\n",
    "x_axis: \n",
    "    formula: \"Axis(start=contact+ 3*r*normal_vec, length=4*r,label='$x$', rotation_angle=-theta)\"\n",
    "fixed: \n",
    "    formula: \"Composition({'angle': angle, 'inclined wall': wall,'ground': ground,'x start': x_const, 'x axis': x_axis})\"\n",
    "fig: \n",
    "    formula: \"Composition({'body': body, 'fixed elements': fixed})\"\n",
    "\"\"\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "gwd = yaml.load(gw)\n",
    "#print(gwd['g'])\n",
    "#yaml.dump(gwd, sys.stdout)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "for _k in list(gwd.keys()):\n",
    "    if _k == \"stop\":\n",
    "        break\n",
    "    _c = gwd[_k]\n",
    "    _t = str(type(_c))\n",
    "    if _k == \"libraries\":\n",
    "        for l in _c:\n",
    "            exec(l,myfig)\n",
    "    #print(_k, _c, _t)\n",
    "    if _t == \"<class 'ruamel.yaml.scalarfloat.ScalarFloat'>\" or \\\n",
    "       _t == \"<class 'str'>\" or _t == \"<class 'int'>\":\n",
    "        _formula = f\"{_k} = {_c}\".replace(\"<bslash>\",\"\\\\\") \n",
    "        #print(_formula)\n",
    "        exec(_formula,myfig)\n",
    "    elif _t == \"<class 'ruamel.yaml.comments.CommentedMap'>\":\n",
    "        #print(_c)\n",
    "        _keys = list(_c.keys())\n",
    "        #print(_keys)\n",
    "        if 'formula' in _keys:\n",
    "            _formula = f\"{_k} = {_c['formula']}\".replace(\"<bslash>\",\"\\\\\")\n",
    "            #print(_formula)\n",
    "            exec(_formula,myfig)\n",
    "        if 'style' in _keys:\n",
    "            for _style in _c[\"style\"]:\n",
    "                #  x_const.set_linestyle('dotted')\n",
    "                _param = _c[\"style\"][_style]\n",
    "                __t = str(type(_param))\n",
    "                #print(__t)\n",
    "                if __t == \"<class 'int'>\":\n",
    "                    _style = f\"{_k}.set_{_style}({_param})\"\n",
    "                else:\n",
    "                    _style = f\"{_k}.set_{_style}('{_param}')\"\n",
    "                #print(_style)\n",
    "                exec(_style,myfig)\n",
    "        if 'transform' in _keys:\n",
    "            #print(_c['transform'])\n",
    "            if str(type(_c['transform'])) == \"<class 'str'>\":\n",
    "                _t = f\"{_k}.{_c['transform']}\"\n",
    "                #print(_t)\n",
    "                exec(_t,myfig)\n",
    "            else:\n",
    "                for _transform in _c[\"transform\"]:\n",
    "                #  x_const.rotate(-theta, contact)\n",
    "                    _t = f\"{_k}.{_c['transform']}\"\n",
    "                    #print(_t)\n",
    "                    exec(_t,myfig)\n",
    "        if \"action\" in _keys:\n",
    "            _action = _c[\"action\"]\n",
    "            #print(_action)\n",
    "            exec(_action,myfig)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "d = myfig['fig'].draw()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "drawing_tool.display()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "def position(t):\n",
    "    \"\"\"Position of center point of wheel.\"\"\"\n",
    "    global tangent_vec,c,myfig\n",
    "    return myfig['c'] + 0.5*myfig['g']*t**2*myfig['tangent_vec']\n",
    "t = 0\n",
    "def move(change):\n",
    "    global fig,t,label,x,x0,dt\n",
    "    dt = change.new - t \n",
    "    t = change.new\n",
    "    drawing_tool.erase()\n",
    "    x = position(t)\n",
    "    x0 = position(t-dt)\n",
    "    displacement = x - x0    \n",
    "    myfig['fig']['body'].translate(displacement)\n",
    "    myfig['fig'].draw()\n",
    "    label.value = f\"p:{x}\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "slider = FloatSlider(\n",
    "    orientation='horizontal',\n",
    "    description='Time:',\n",
    "    value=0.0,\n",
    "    min=0.0,\n",
    "    max=1.0,\n",
    "    step = 1.0 / 30\n",
    ")\n",
    "slider.observe(move, 'value')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "label = Label(value=\"\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "0d24f9e65d9642e49d10453a9b1f0cfe",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "AppLayout(children=(HBox(children=(FloatSlider(value=0.0, description='Time:', max=1.0, step=0.033333333333333…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "AppLayout(\n",
    "    center=drawing_tool.mpl.gcf().canvas,\n",
    "    footer=HBox([slider,label]),\n",
    "    pane_heights=[0, 6, 1]\n",
    ")\n",
    "#drawing_tool.mpl.ion()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "drawing_tool.mpl.ion()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "myfig['fig'].graphviz_dot('fig')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "!dot -Tpng -o fig.png fig.dot"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "![](fig.png)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "drawing_tool.mpl.gcf().canvas.print_png(\"fig1.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "![](fig1.png)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" height=\"345.6pt\" version=\"1.1\" viewBox=\"0 0 460.8 345.6\" width=\"460.8pt\">\n",
       " <metadata>\n",
       "  <rdf:RDF xmlns:cc=\"http://creativecommons.org/ns#\" xmlns:dc=\"http://purl.org/dc/elements/1.1/\" xmlns:rdf=\"http://www.w3.org/1999/02/22-rdf-syntax-ns#\">\n",
       "   <cc:Work>\n",
       "    <dc:type rdf:resource=\"http://purl.org/dc/dcmitype/StillImage\"/>\n",
       "    <dc:date>2020-08-01T10:34:52.138957</dc:date>\n",
       "    <dc:format>image/svg+xml</dc:format>\n",
       "    <dc:creator>\n",
       "     <cc:Agent>\n",
       "      <dc:title>Matplotlib v3.3.0, https://matplotlib.org/</dc:title>\n",
       "     </cc:Agent>\n",
       "    </dc:creator>\n",
       "   </cc:Work>\n",
       "  </rdf:RDF>\n",
       " </metadata>\n",
       " <defs>\n",
       "  <style type=\"text/css\">*{stroke-linecap:butt;stroke-linejoin:round;}</style>\n",
       " </defs>\n",
       " <g id=\"figure_1\">\n",
       "  <g id=\"patch_1\">\n",
       "   <path d=\"M 0 345.6  L 460.8 345.6  L 460.8 0  L 0 0  z \" style=\"fill:#ffffff;\"/>\n",
       "  </g>\n",
       "  <g id=\"axes_1\">\n",
       "   <g id=\"line2d_1\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 221.310338 137.819034  L 221.179916 135.33042  L 220.790076 132.869072  L 220.145092 130.461957  L 219.252029 128.135448  L 218.120671 125.915034  L 216.763415 123.825043  L 215.19513 121.888372  L 213.433 120.126242  L 211.49633 118.557957  L 209.406338 117.200701  L 207.185924 116.069344  L 204.859415 115.17628  L 202.4523 114.531296  L 199.990952 114.141457  L 197.502338 114.011034  L 195.013725 114.141457  L 192.552377 114.531296  L 190.145262 115.17628  L 187.818752 116.069344  L 185.598338 117.200701  L 183.508347 118.557957  L 181.571677 120.126242  L 179.809546 121.888372  L 178.241262 123.825043  L 176.884006 125.915034  L 175.752648 128.135448  L 174.859585 130.461957  L 174.2146 132.869072  L 173.824761 135.33042  L 173.694338 137.819034  L 173.824761 140.307648  L 174.2146 142.768995  L 174.859585 145.176111  L 175.752648 147.50262  L 176.884006 149.723034  L 178.241262 151.813025  L 179.809546 153.749695  L 181.571677 155.511826  L 183.508347 157.080111  L 185.598338 158.437367  L 187.818752 159.568724  L 190.145262 160.461787  L 192.552377 161.106772  L 195.013725 161.496611  L 197.502338 161.627034  L 199.990952 161.496611  L 202.4523 161.106772  L 204.859415 160.461787  L 207.185924 159.568724  L 209.406338 158.437367  L 211.49633 157.080111  L 213.433 155.511826  L 215.19513 153.749695  L 216.763415 151.813025  L 218.120671 149.723034  L 219.252029 147.50262  L 220.145092 145.176111  L 220.790076 142.768995  L 221.179916 140.307648  L 221.310338 137.819034  L 221.310338 137.819034  \" style=\"fill:none;stroke:#808080;stroke-linecap:square;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_2\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 209.406338 137.819034  L 209.290489 136.162317  L 208.945198 134.537847  L 208.377183 132.977241  L 207.597503 131.510875  L 206.621331 130.16729  L 205.467669 128.972638  L 204.158971 127.950171  L 202.720708 127.11979  L 201.180877 126.497657  L 199.569446 126.095882  L 197.917782 125.922286  L 196.258032 125.980245  L 194.6225 126.268634  L 193.043021 126.781837  L 191.550338 127.509868  L 190.173504 128.438554  L 188.939317 129.549821  L 187.8718 130.822038  L 186.99173 132.230444  L 186.316237 133.747626  L 185.858469 135.344053  L 185.627336 136.988653  L 185.627336 138.649415  L 185.858469 140.294015  L 186.316237 141.890442  L 186.99173 143.407623  L 187.8718 144.81603  L 188.939317 146.088247  L 190.173504 147.199514  L 191.550338 148.1282  L 193.043021 148.856231  L 194.6225 149.369434  L 196.258032 149.657823  L 197.917782 149.715782  L 199.569446 149.542185  L 201.180877 149.140411  L 202.720708 148.518278  L 204.158971 147.687897  L 205.467669 146.66543  L 206.621331 145.470778  L 207.597503 144.127193  L 208.377183 142.660827  L 208.945198 141.100221  L 209.290489 139.475751  L 209.406338 137.819034  L 209.406338 137.819034  \" style=\"fill:none;stroke:#808080;stroke-linecap:square;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_2\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 217.310338 133.819034  L 217.295835 132.988147  L 217.252343 132.158272  L 217.179916 131.33042  L 217.078641 130.505601  L 216.948641 129.684818  L 216.790076 128.869072  L 216.603139 128.059357  L 216.388057 127.25666  L 216.145092 126.461957  L 215.87454 125.676218  L 215.576732 124.9004  L 215.252029 124.135448  L 214.900827 123.382294  L 214.523555 122.641855  L 214.120671 121.915034  L 213.692667 121.202716  L 213.240065 120.505769  L 212.763415 119.825043  L 212.263298 119.161366  L 211.740324 118.515547  L 211.19513 117.888372  L 210.62838 117.280607  L 210.040765 116.692992  L 209.433 116.126242  L 208.805826 115.581048  L 208.160007 115.058074  L 207.49633 114.557957  L 206.815603 114.081307  L 206.118656 113.628705  L 205.406338 113.200701  L 204.679517 112.797818  L 203.939079 112.420545  L 203.185924 112.069344  L 202.420972 111.744641  L 201.645154 111.446832  L 200.859415 111.17628  L 200.064713 110.933315  L 199.262015 110.718233  L 198.4523 110.531296  L 197.636554 110.372731  L 196.815772 110.242732  L 195.990952 110.141457  L 195.1631 110.069029  L 194.333226 110.025537  L 193.502338 110.011034  L 192.671451 110.025537  L 191.841576 110.069029  L 191.013725 110.141457  L 190.188905 110.242732  L 189.368123 110.372731  L 188.552377 110.531296  L 187.742662 110.718233  L 186.939964 110.933315  L 186.145262 111.17628  L 185.359523 111.446832  L 184.583705 111.744641  L 183.818752 112.069344  L 183.065598 112.420545  L 182.325159 112.797818  L 181.598338 113.200701  L 180.886021 113.628705  L 180.189074 114.081307  L 179.508347 114.557957  L 178.84467 115.058074  L 178.198851 115.581048  L 177.571677 116.126242  L 176.963912 116.692992  L 176.376296 117.280607  L 175.809546 117.888372  L 175.264352 118.515547  L 174.741378 119.161366  L 174.241262 119.825043  L 173.764612 120.505769  L 173.312009 121.202716  L 172.884006 121.915034  L 172.481122 122.641855  L 172.10385 123.382294  L 171.752648 124.135448  L 171.427945 124.9004  L 171.130136 125.676218  L 170.859585 126.461957  L 170.61662 127.25666  L 170.401538 128.059357  L 170.2146 128.869072  L 170.056035 129.684818  L 169.926036 130.505601  L 169.824761 131.33042  L 169.752333 132.158272  L 169.708842 132.988147  L 169.694338 133.819034  L 169.708842 134.649921  L 169.752333 135.479796  L 169.824761 136.307648  L 169.926036 137.132467  L 170.056035 137.95325  L 170.2146 138.768995  L 170.401538 139.57871  L 170.61662 140.381408  L 170.859585 141.176111  L 171.130136 141.961849  L 171.427945 142.737668  L 171.752648 143.50262  L 172.10385 144.255774  L 172.481122 144.996213  L 172.884006 145.723034  L 173.312009 146.435352  L 173.764612 147.132299  L 174.241262 147.813025  L 174.741378 148.476702  L 175.264352 149.122521  L 175.809546 149.749695  L 176.376296 150.35746  L 176.963912 150.945076  L 177.571677 151.511826  L 178.198851 152.05702  L 178.84467 152.579994  L 179.508347 153.080111  L 180.189074 153.55676  L 180.886021 154.009363  L 181.598338 154.437367  L 182.325159 154.84025  L 183.065598 155.217523  L 183.818752 155.568724  L 184.583705 155.893427  L 185.359523 156.191236  L 186.145262 156.461787  L 186.939964 156.704752  L 187.742662 156.919835  L 188.552377 157.106772  L 189.368123 157.265337  L 190.188905 157.395336  L 191.013725 157.496611  L 191.841576 157.569039  L 192.671451 157.612531  L 193.502338 157.627034  L 194.333226 157.612531  L 195.1631 157.569039  L 195.990952 157.496611  L 196.815772 157.395336  L 197.636554 157.265337  L 198.4523 157.106772  L 199.262015 156.919835  L 200.064713 156.704752  L 200.859415 156.461787  L 201.645154 156.191236  L 202.420972 155.893427  L 203.185924 155.568724  L 203.939079 155.217523  L 204.679517 154.84025  L 205.406338 154.437367  L 206.118656 154.009363  L 206.815603 153.55676  L 207.49633 153.080111  L 208.160007 152.579994  L 208.805826 152.05702  L 209.433 151.511826  L 210.040765 150.945076  L 210.62838 150.35746  L 211.19513 149.749695  L 211.740324 149.122521  L 212.263298 148.476702  L 212.763415 147.813025  L 213.240065 147.132299  L 213.692667 146.435352  L 214.120671 145.723034  L 214.523555 144.996213  L 214.900827 144.255774  L 215.252029 143.50262  L 215.576732 142.737668  L 215.87454 141.961849  L 216.145092 141.176111  L 216.388057 140.381408  L 216.603139 139.57871  L 216.790076 138.768995  L 216.948641 137.95325  L 217.078641 137.132467  L 217.179916 136.307648  L 217.252343 135.479796  L 217.295835 134.649921  z \" style=\"fill:#0000ff;stroke:#000000;stroke-linejoin:miter;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_3\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 205.406338 133.819034  L 205.399087 133.40359  L 205.377341 132.988653  L 205.341127 132.574727  L 205.290489 132.162317  L 205.22549 131.751926  L 205.146207 131.344053  L 205.052739 130.939196  L 204.945198 130.537847  L 204.823715 130.140496  L 204.688439 129.747626  L 204.539535 129.359717  L 204.377183 128.977241  L 204.201583 128.600664  L 204.012947 128.230444  L 203.811505 127.867034  L 203.597503 127.510875  L 203.371202 127.162402  L 203.132877 126.822038  L 202.882818 126.4902  L 202.621331 126.16729  L 202.348734 125.853703  L 202.065359 125.549821  L 201.771552 125.256013  L 201.467669 124.972638  L 201.154082 124.700041  L 200.831173 124.438554  L 200.499334 124.188496  L 200.158971 123.950171  L 199.810497 123.723869  L 199.454338 123.509868  L 199.090928 123.308426  L 198.720708 123.11979  L 198.344131 122.944189  L 197.961655 122.781837  L 197.573746 122.632933  L 197.180877 122.497657  L 196.783525 122.376175  L 196.382177 122.268634  L 195.977319 122.175165  L 195.569446 122.095882  L 195.159055 122.030883  L 194.746645 121.980245  L 194.332719 121.944031  L 193.917782 121.922286  L 193.502338 121.915034  L 193.086895 121.922286  L 192.671957 121.944031  L 192.258032 121.980245  L 191.845622 122.030883  L 191.43523 122.095882  L 191.027358 122.175165  L 190.6225 122.268634  L 190.221151 122.376175  L 189.8238 122.497657  L 189.430931 122.632933  L 189.043021 122.781837  L 188.660545 122.944189  L 188.283968 123.11979  L 187.913749 123.308426  L 187.550338 123.509868  L 187.194179 123.723869  L 186.845706 123.950171  L 186.505343 124.188496  L 186.173504 124.438554  L 185.850595 124.700041  L 185.537008 124.972638  L 185.233125 125.256013  L 184.939317 125.549821  L 184.655942 125.853703  L 184.383345 126.16729  L 184.121858 126.4902  L 183.8718 126.822038  L 183.633475 127.162402  L 183.407174 127.510875  L 183.193172 127.867034  L 182.99173 128.230444  L 182.803094 128.600664  L 182.627493 128.977241  L 182.465142 129.359717  L 182.316237 129.747626  L 182.180962 130.140496  L 182.059479 130.537847  L 181.951938 130.939196  L 181.858469 131.344053  L 181.779187 131.751926  L 181.714187 132.162317  L 181.66355 132.574727  L 181.627336 132.988653  L 181.60559 133.40359  L 181.598338 133.819034  L 181.60559 134.234478  L 181.627336 134.649415  L 181.66355 135.063341  L 181.714187 135.475751  L 181.779187 135.886142  L 181.858469 136.294015  L 181.951938 136.698872  L 182.059479 137.100221  L 182.180962 137.497572  L 182.316237 137.890442  L 182.465142 138.278351  L 182.627493 138.660827  L 182.803094 139.037404  L 182.99173 139.407623  L 183.193172 139.771034  L 183.407174 140.127193  L 183.633475 140.475666  L 183.8718 140.81603  L 184.121858 141.147868  L 184.383345 141.470778  L 184.655942 141.784365  L 184.939317 142.088247  L 185.233125 142.382055  L 185.537008 142.66543  L 185.850595 142.938027  L 186.173504 143.199514  L 186.505343 143.449572  L 186.845706 143.687897  L 187.194179 143.914198  L 187.550338 144.1282  L 187.913749 144.329642  L 188.283968 144.518278  L 188.660545 144.693879  L 189.043021 144.856231  L 189.430931 145.005135  L 189.8238 145.140411  L 190.221151 145.261893  L 190.6225 145.369434  L 191.027358 145.462903  L 191.43523 145.542185  L 191.845622 145.607185  L 192.258032 145.657823  L 192.671957 145.694036  L 193.086895 145.715782  L 193.502338 145.723034  L 193.917782 145.715782  L 194.332719 145.694036  L 194.746645 145.657823  L 195.159055 145.607185  L 195.569446 145.542185  L 195.977319 145.462903  L 196.382177 145.369434  L 196.783525 145.261893  L 197.180877 145.140411  L 197.573746 145.005135  L 197.961655 144.856231  L 198.344131 144.693879  L 198.720708 144.518278  L 199.090928 144.329642  L 199.454338 144.1282  L 199.810497 143.914198  L 200.158971 143.687897  L 200.499334 143.449572  L 200.831173 143.199514  L 201.154082 142.938027  L 201.467669 142.66543  L 201.771552 142.382055  L 202.065359 142.088247  L 202.348734 141.784365  L 202.621331 141.470778  L 202.882818 141.147868  L 203.132877 140.81603  L 203.371202 140.475666  L 203.597503 140.127193  L 203.811505 139.771034  L 204.012947 139.407623  L 204.201583 139.037404  L 204.377183 138.660827  L 204.539535 138.278351  L 204.688439 137.890442  L 204.823715 137.497572  L 204.945198 137.100221  L 205.052739 136.698872  L 205.146207 136.294015  L 205.22549 135.886142  L 205.290489 135.475751  L 205.341127 135.063341  L 205.377341 134.649415  L 205.399087 134.234478  z \" style=\"fill:#ffffff;stroke:#000000;stroke-linejoin:miter;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_4\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 181.598338 154.437367  L 180.183307 161.352273  L 178.260648 160.242225  L 157.800648 195.679984  L 157.780029 195.66808  L 178.240029 160.230321  L 176.31737 159.120273  z \" style=\"stroke:#000000;stroke-linejoin:miter;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_5\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 193.502338 205.243034  L 191.270338 198.547034  L 193.490434 198.547034  L 193.490434 133.819034  L 193.514242 133.819034  L 193.514242 198.547034  L 195.734338 198.547034  z \" style=\"stroke:#000000;stroke-linejoin:miter;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_6\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 81.408 96.592448  L 319.488 234.048  L 319.488 240  L 81.408 102.544448  z \" style=\"fill:url(#h02041ea892);stroke:#000000;stroke-linejoin:miter;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_3\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 157.790338 195.674032  L 181.598338 154.437367  \" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_4\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 193.502338 133.819034  L 193.502338 205.243034  \" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_5\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 257.633002 198.336  L 255.754697 201.807439  L 254.071024 205.377335  L 252.587126 209.034788  L 251.307533 212.768628  L 250.236154 216.567451  L 249.37626 220.419658  L 248.730478 224.313484  L 248.300779 228.237038  L 248.088475 232.178338  L 248.064 234.048  L 248.064 234.048  \" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_6\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 295.68 234.048  L 200.448 234.048  \" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_7\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 152.832 137.829114  L 200.448 55.355782  \" style=\"fill:none;stroke:#000000;stroke-dasharray:2,3.3;stroke-dashoffset:0;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_8\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 188.544 75.974115  L 271.017331 123.590115  \" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_9\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 267.446131 117.404615  L 271.017331 123.590115  \" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_10\">\n",
       "    <path clip-path=\"url(#p1e644576b7)\" d=\"M 263.874931 123.590115  L 271.017331 123.590115  \" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-width:2;\"/>\n",
       "   </g>\n",
       "   <g id=\"text_1\">\n",
       "    <!-- $N$ -->\n",
       "    <g transform=\"translate(147.26308 204.436824)scale(0.14 -0.14)\">\n",
       "     <defs>\n",
       "      <path d=\"M 74 57.53125  C 75.09375 61.8125 76.703125 64.796875 84.296875 65.09375  C 84.59375 65.09375 85.796875 65.1875 85.796875 66.890625  C 85.796875 68 84.90625 68 84.5 68  C 82.5 68 77.40625 67.796875 75.40625 67.796875  L 70.59375 67.796875  C 69.203125 67.796875 67.40625 68 66 68  C 65.40625 68 64.203125 68 64.203125 66.09375  C 64.203125 65.09375 65 65.09375 65.703125 65.09375  C 71.703125 64.890625 72.09375 62.609375 72.09375 60.8125  C 72.09375 59.921875 72 59.625 71.703125 58.21875  L 60.40625 13.34375  L 39 66.296875  C 38.296875 67.890625 38.203125 68 36 68  L 23.796875 68  C 21.796875 68 20.90625 68 20.90625 66.09375  C 20.90625 65.09375 21.59375 65.09375 23.5 65.09375  C 24 65.09375 29.90625 65.09375 29.90625 64.203125  C 29.90625 64 29.703125 63.203125 29.59375 62.90625  L 16.296875 10.15625  C 15.09375 5.28125 12.703125 3.1875 6.09375 2.890625  C 5.59375 2.890625 4.59375 2.796875 4.59375 1  C 4.59375 0 5.59375 0 5.90625 0  C 7.90625 0 13 0.1875 15 0.1875  L 19.796875 0.1875  C 21.203125 0.1875 22.90625 0 24.296875 0  C 25 0 26.09375 0 26.09375 1.890625  C 26.09375 2.796875 25.09375 2.890625 24.703125 2.890625  C 21.40625 2.984375 18.203125 3.59375 18.203125 7.171875  C 18.203125 7.96875 18.40625 8.859375 18.59375 9.65625  L 32.09375 62.90625  C 32.703125 61.90625 32.703125 61.703125 33.09375 60.8125  L 56.90625 1.796875  C 57.40625 0.59375 57.59375 0 58.5 0  C 59.5 0 59.59375 0.296875 60 2  z \" id=\"CMMI12-78\"/>\n",
       "     </defs>\n",
       "     <use transform=\"scale(0.996264)\" xlink:href=\"#CMMI12-78\"/>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"text_2\">\n",
       "    <!-- $Mg$ -->\n",
       "    <g transform=\"translate(182.897978 215.361434)scale(0.14 -0.14)\">\n",
       "     <defs>\n",
       "      <path d=\"M 90.796875 60.71875  C 91.703125 64.09375 91.90625 65.09375 99 65.09375  C 100.90625 65.09375 101.796875 65.09375 101.796875 66.984375  C 101.796875 68 101.09375 68 99.203125 68  L 87.203125 68  C 84.703125 68 84.59375 67.890625 83.5 66.296875  L 47 8.859375  L 39.5 65.796875  C 39.203125 68 39.09375 68 36.5 68  L 24.09375 68  C 22.203125 68 21.296875 68 21.296875 66.09375  C 21.296875 65.09375 22.203125 65.09375 23.703125 65.09375  C 29.796875 65.09375 29.796875 64.296875 29.796875 63.203125  C 29.796875 63 29.796875 62.40625 29.40625 60.90625  L 16.59375 10.15625  C 15.40625 5.375 13.09375 3.1875 6.40625 2.890625  C 6.09375 2.890625 4.90625 2.796875 4.90625 1.09375  C 4.90625 0 5.796875 0 6.203125 0  C 8.203125 0 13.296875 0.1875 15.296875 0.1875  L 20.09375 0.1875  C 21.5 0.1875 23.203125 0 24.59375 0  C 25.296875 0 26.40625 0 26.40625 1.890625  C 26.40625 2.796875 25.40625 2.890625 25 2.890625  C 21.703125 2.984375 18.5 3.59375 18.5 7.171875  C 18.5 8.171875 18.5 8.265625 18.90625 9.65625  L 32.703125 64.5  L 32.796875 64.5  L 41.09375 2.6875  C 41.40625 0.296875 41.5 0 42.40625 0  C 43.5 0 44 0.796875 44.5 1.6875  L 84.703125 64.984375  L 84.796875 64.984375  L 70.296875 7.375  C 69.40625 3.890625 69.203125 2.890625 62.203125 2.890625  C 60.296875 2.890625 59.296875 2.890625 59.296875 1.09375  C 59.296875 0 60.203125 0 60.796875 0  C 62.5 0 64.5 0.1875 66.203125 0.1875  L 78 0.1875  C 79.703125 0.1875 81.796875 0 83.5 0  C 84.296875 0 85.40625 0 85.40625 1.890625  C 85.40625 2.890625 84.5 2.890625 83 2.890625  C 76.90625 2.890625 76.90625 3.6875 76.90625 4.6875  C 76.90625 4.78125 76.90625 5.484375 77.09375 6.28125  z \" id=\"CMMI12-77\"/>\n",
       "      <path d=\"M 33.796875 12.671875  C 33.40625 11.0625 33.203125 10.671875 31.90625 9.171875  C 27.796875 3.890625 23.59375 2 20.5 2  C 17.203125 2 14.09375 4.578125 14.09375 11.46875  C 14.09375 16.75 17.09375 27.921875 19.296875 32.421875  C 22.203125 38 26.703125 41.984375 30.90625 41.984375  C 37.5 41.984375 38.796875 33.8125 38.796875 33.21875  L 38.5 31.8125  z M 40 37.40625  C 38.703125 40.296875 35.90625 44 30.90625 44  C 20 44 7.59375 30.328125 7.59375 15.453125  C 7.59375 5.078125 13.90625 0 20.296875 0  C 25.59375 0 30.296875 4.1875 32.09375 6.1875  L 29.90625 -2.75  C 28.5 -8.125 27.90625 -10.578125 24.296875 -14  C 20.203125 -18 16.40625 -18 14.203125 -18  C 11.203125 -18 8.703125 -17.8125 6.203125 -17.03125  C 9.40625 -16.15625 10.203125 -13.40625 10.203125 -12.328125  C 10.203125 -10.765625 9 -9.203125 6.796875 -9.203125  C 4.40625 -9.203125 1.796875 -11.15625 1.796875 -14.390625  C 1.796875 -18.40625 5.90625 -20 14.40625 -20  C 27.296875 -20 34 -11.84375 35.296875 -6.5625  L 46.40625 37.90625  C 46.703125 39.09375 46.703125 39.296875 46.703125 39.5  C 46.703125 40.890625 45.59375 42 44.09375 42  C 41.703125 42 40.296875 40 40 37.3125  z \" id=\"CMMI12-103\"/>\n",
       "     </defs>\n",
       "     <use transform=\"scale(0.996264)\" xlink:href=\"#CMMI12-77\"/>\n",
       "     <use transform=\"translate(104.780079 0)scale(0.996264)\" xlink:href=\"#CMMI12-103\"/>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"text_3\">\n",
       "    <!-- $\\theta$ -->\n",
       "    <g transform=\"translate(240.657146 214.021618)scale(0.18 -0.18)\">\n",
       "     <defs>\n",
       "      <path d=\"M 44.296875 50  C 44.296875 60.140625 41.09375 70 32.90625 70  C 18.90625 70 4 40.859375 4 18.984375  C 4 14.40625 5 -1 15.5 -1  C 29.09375 -1 44.296875 27.421875 44.296875 50  z M 14 36  C 15.5 41.859375 17.59375 50.203125 21.59375 57.265625  C 24.90625 63.21875 28.40625 67.984375 32.796875 67.984375  C 36.09375 67.984375 38.296875 65.203125 38.296875 55.671875  C 38.296875 52.09375 38 47.125 35.09375 36  z M 34.40625 33  C 31.90625 23.25 29.796875 17 26.203125 10.734375  C 23.296875 5.65625 19.796875 1 15.59375 1  C 12.5 1 10 3.375 10 13.21875  C 10 19.671875 11.59375 26.4375 13.203125 33  z \" id=\"CMMI12-18\"/>\n",
       "     </defs>\n",
       "     <use transform=\"scale(0.996264)\" xlink:href=\"#CMMI12-18\"/>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"text_4\">\n",
       "    <!-- $x$ -->\n",
       "    <g transform=\"translate(274.009725 127.558115)scale(0.14 -0.14)\">\n",
       "     <defs>\n",
       "      <path d=\"M 47.40625 41  C 44.203125 40.390625 43 37.984375 43 36.09375  C 43 33.6875 44.90625 32.890625 46.296875 32.890625  C 49.296875 32.890625 51.40625 35.484375 51.40625 38.1875  C 51.40625 42.09375 46.59375 44 42.40625 44  C 36.296875 44 32.90625 38.296875 32 36.390625  C 29.703125 43.59375 23.5 44 21.703125 44  C 11.5 44 6.09375 31.140625 6.09375 28.890625  C 6.09375 28.484375 6.5 27.96875 7.203125 27.96875  C 8 27.96875 8.203125 28.578125 8.40625 28.984375  C 11.796875 40.1875 18.5 41.984375 21.40625 41.984375  C 25.90625 41.984375 26.796875 38.09375 26.796875 35.6875  C 26.796875 33.484375 26.203125 31.140625 25 26.21875  L 21.59375 12.15625  C 20.09375 6.09375 17.203125 1 11.90625 1  C 11.40625 1 8.90625 1 6.796875 1.96875  C 10.40625 2.59375 11.203125 5.59375 11.203125 6.796875  C 11.203125 8.796875 9.703125 10 7.796875 10  C 5.40625 10 2.796875 7.890625 2.796875 4.6875  C 2.796875 0.890625 7.5 -1 11.796875 -1  C 16.59375 -1 20 2.390625 22.09375 6.484375  C 23.703125 1 28.703125 -1 32.40625 -1  C 42.59375 -1 48 11.734375 48 14  C 48 14.515625 47.59375 14.921875 47 14.921875  C 46.09375 14.921875 46 14.40625 45.703125 13.59375  C 43 4.6875 37.203125 1 32.703125 1  C 29.203125 1 27.296875 3.1875 27.296875 7.296875  C 27.296875 9.484375 27.703125 11.125 29.296875 17.90625  L 32.796875 31.859375  C 34.296875 37.890625 37.703125 41.984375 42.296875 41.984375  C 42.5 41.984375 45.296875 41.984375 47.40625 41  z \" id=\"CMMI12-120\"/>\n",
       "     </defs>\n",
       "     <use transform=\"scale(0.996264)\" xlink:href=\"#CMMI12-120\"/>\n",
       "    </g>\n",
       "   </g>\n",
       "  </g>\n",
       " </g>\n",
       " <defs>\n",
       "  <clipPath id=\"p1e644576b7\">\n",
       "   <rect height=\"261.888\" width=\"357.12\" x=\"57.6\" y=\"43.584\"/>\n",
       "  </clipPath>\n",
       " </defs>\n",
       " <defs>\n",
       "  <pattern height=\"72\" id=\"h02041ea892\" patternUnits=\"userSpaceOnUse\" width=\"72\" x=\"0\" y=\"0\">\n",
       "   <rect fill=\"#ffffff\" height=\"73\" width=\"73\" x=\"0\" y=\"0\"/>\n",
       "   <path d=\"M -36 36  L 36 -36  M -24 48  L 48 -24  M -12 60  L 60 -12  M 0 72  L 72 0  M 12 84  L 84 12  M 24 96  L 96 24  M 36 108  L 108 36  \" style=\"fill:#000000;stroke:#000000;stroke-linecap:butt;stroke-linejoin:miter;stroke-width:1.0;\"/>\n",
       "  </pattern>\n",
       " </defs>\n",
       "</svg>"
      ],
      "text/plain": [
       "<IPython.core.display.SVG object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from io import BytesIO\n",
    "f = BytesIO()\n",
    "drawing_tool.mpl.savefig(f, format=\"svg\")\n",
    "display(SVG(f.getvalue()))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}