manual_teachbooks/book/instructie.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "````{margin}\n",
    "```{attributiongrey} Attribution\n",
    ":class: attribution\n",
    "\n",
    "Deze pagina is aangepast van https://github.com/TUDelft-books/CT1000/blob/2025/book/week_2/session_1/intro.ipynb\n",
    "\n",
    "% Figures from https://github.com/TUDelft-books/CT1000/blob/2025/book/week_2/session_1/intro_data/Tekening1.vsdx\n",
    "\n",
    "```\n",
    "```` "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Instructie_oud\n",
    "\n",
    "In constructiemechanica 2 heb je geleerd hoe je met differentiaalvergelijkingen voor statisch bepaalde constructies de krachtsverdeling en vervormingen kan vinden. Precies dezelfde aanpak werkt ook voor statisch onbepaalde constructie, met als enige verschil dat we de inwendige krachten niet kunnen vinden zonder ook de vergelijkingen en randvoorwaarden voor verplaatsingen op te stellen. Dat doen we volgens de volgende stappen:\n",
    "\n",
    "::::::{prf:algorithm} Differentiaalvergelijkingen toepassen\n",
    ":nonumber: true\n",
    ":label: differentiaalvergelijkingen_algoritme\n",
    "\n",
    "1. Eerst definiëren we de segmenten zonder discontinuïteiten. Daarna kan elk segment afzonderlijk worden getekend, waarbij de belastingen, coördinatenstelsels, afmetingen en snedekrachten aan de randen van elk segment worden aangegeven.\n",
    "2. Nu kunnen we voor elk element de belastingsvergelijking opstellen en volgens de relaties van de differentiaalvergelijkingen alle krachtsgrootheden en verplaatsingen opstellen. Daarbij kunnen we direct ook de randvoorwaarden opstellen.\n",
    "3. Er zijn nog meer vergelijkingen nodig om alle onbekende integratieconstantes op te lossen, deze kunnen we vinden door de overgangen tussen de verschillende elementen te analyseren waarbij zowel verplaatsingen als krachten in overeenstemming moeten zijn.\n",
    "4. Er zijn nu evenveel vergelijkingen als onbekenden opgesteld. Deze kunnen nu worden opgelost en ingevuld in de opgestelde vergelijkingen voor het gewenste resultaat. Omdat het aantal onbekenden bij statisch onbepaalde constructies als snel uit de hand loopt wordt wordt aangeraden om je rekenmachine of een symbolische wiskundetool te gebruiken om dit op te lossen. Om het stelsel vergelijkingen op te lossen met behulp van numerieke hulpmiddelen, kun je de vergelijkingen opschrijven als een matrixformulering $Ax=b$ die eenvoudig kan worden opgelost door bijvoorbeeld je grafische rekenmachine of in python. Als alternatief kan een symbolische programmeertaal worden gebruikt om de vergelijkingen exact op te lossen. Als voorbeeld wordt het [`sympy`](https://docs.sympy.org/latest/index.html) pakket in `Python` gebruikt:\n",
    "\n",
    "::::::\n",
    "\n",
    "Als voorbeeld bepalen we de een simpele statisch onbepaalde constructie:\n",
    "\n",
    "::::::{prf:example}\n",
    ":nonumber: true\n",
    ":label: df_example_0\n",
    "\n",
    "Gegeven een constructie zoals hieronder weergegeven:\n",
    "\n",
    "```{figure} intro_data/structure.svg\n",
    "---\n",
    "align: center\n",
    "---\n",
    "Constructie, $EI = 5000 \\ \\rm{kNm^2}$\n",
    "```\n",
    "::::::\n",
    "\n",
    "1. Eerst definiëren we de segmenten zonder discontinuïteiten. Daarna kan elk segment afzonderlijk worden getekend, waarbij de belastingen, coördinatenstelsels, afmetingen en snedekrachten aan de randen van elk segment worden aangegeven.\n",
    "\n",
    "::::::{prf:example}\n",
    ":nonumber: true\n",
    ":label: df_example_1\n",
    "\n",
    "In dit voorbeeld moeten we element AB in twee delen opsplitsen omdat de verdeelde belasting een discontinuïteit veroorzaakt.\n",
    "\n",
    "Voor AC leidt dit tot:\n",
    "\n",
    "```{figure} intro_data/FBD_AC.svg\n",
    ":align: center\n",
    "\n",
    "Vrijlichaamdiagram element AC\n",
    "```\n",
    "\n",
    "Voor BC leidt dit tot:\n",
    "\n",
    "```{figure} intro_data/FBD_BC.svg\n",
    ":align: center\n",
    "\n",
    "Vrijlichaamdiagram element BC\n",
    "```\n",
    "::::::"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "2. Nu kunnen we voor elk element de belastingsvergelijking opstellen en volgens de relaties van de differentiaalvergelijkingen alle krachtsgrootheden en verplaatsingen opstellen. Daarbij kunnen we direct ook de randvoorwaarden opstellen.\n",
    "\n",
    "::::::{prf:example}\n",
    ":nonumber: true\n",
    ":label: df_example_3\n",
    "\n",
    "In dit voorbeeld is er enkel invloed door buiging, dus hebben we aan deze differentiaalvergelijking genoeg: $ -EI \\cfrac{{d}^4w}{{d}x^4} = -q_z$. De differentiaalvergelijkingen voor rek $ EA \\cfrac{{d}^2u}{{d}x^2} = -q_x$ en torsie $ GI \\cfrac{{d}^2\\varphi_x}{{d}x^2} = -q_T$ hebben we in dit geval niet nodig.\n",
    "\n",
    "**Element AC**\n",
    "\n",
    "Omdat dit element onbelast is, geeft de belastingsvergelijking voor dit element:\n",
    "\n",
    "$$q_{z,{\\rm{AC}}}=0$$\n",
    "\n",
    "wat leidt tot de vergelijkingen:\n",
    "\n",
    "$$ V_{\\rm{AC}} \\left(x\\right) = C_1$$\n",
    "\n",
    "$$ M_{\\rm{AC}} \\left(x\\right) = C_1 \\cdot x + C_2 $$\n",
    "\n",
    "$$ \\kappa_{\\rm{AC}} \\left(x\\right) = \\cfrac{C_1 \\cdot x + C_2}{5000}$$\n",
    "\n",
    "$$ \\varphi_{\\rm{AC}} \\left(x\\right) = \\cfrac{\\cfrac{1}{2}\\cdot C_1 \\cdot x^2 + C_2 \\cdot x}{5000} + C_3$$\n",
    "\n",
    "$$ w_{\\rm{AC}} \\left(x\\right) = \\cfrac{-\\cfrac{1}{6}\\cdot C_1 \\cdot x^3 - \\cfrac{1}{2} C_2 \\cdot x^2}{5000} - C_3 \\cdot x + C_4$$\n",
    "\n",
    "Vanwege de oplegging bij A kan een randvoorwaarde worden geformuleerd als:\n",
    "\n",
    "```{math}\n",
    ":label: BC_1\n",
    "w_{\\rm{AC}} \\left(0\\right) = 0\n",
    "```\n",
    "\n",
    "Omdat A ook een onbelast uiteinde is, kan een andere randvoorwaarde worden geformuleerd als:\n",
    "\n",
    "```{math}\n",
    ":label: BC_2\n",
    "M_{\\rm{AC}} \\left(0\\right) = 0\n",
    "```\n",
    "\n",
    "Als dit uiteinde belast werd kon met een krachten- of momentenevenwicht ook een relatie gevonden worden.\n",
    "\n",
    "**Element BC**\n",
    "\n",
    "Omdat dit element belast is, geeft de belastingsvergelijking voor dit element:\n",
    "\n",
    "$$q_{z,\\rm{BC}}=10$$\n",
    "\n",
    "wat leidt tot de vergelijkingen:\n",
    "\n",
    "$$ V_{\\rm{BC}} \\left(x\\right) = -10 \\cdot x + C_5$$\n",
    "\n",
    "$$ M_{\\rm{BC}} \\left(x\\right) = -5 \\cdot x^2 + C_5 \\cdot x + C_6 $$\n",
    "\n",
    "$$ \\kappa_{\\rm{BC}} \\left(x\\right) = \\cfrac{-5 \\cdot x^2 + C_5 \\cdot x + C_6}{5000}$$\n",
    "\n",
    "$$ \\varphi_{\\rm{BC}} \\left(x\\right) = \\cfrac{-\\cfrac{5}{3} \\cdot x^3 + \\cfrac{1}{2}\\cdot C_5 \\cdot x^2 + C_6 \\cdot x}{5000} + C_7$$\n",
    "\n",
    "$$ w_{\\rm{BC}} \\left(x\\right) = \\cfrac{\\cfrac{5}{12} \\cdot x^4 -\\cfrac{1}{6}\\cdot C_5 \\cdot x^3 - \\cfrac{1}{2} C_6 \\cdot x^2}{5000} - C_7 \\cdot x + C_8$$\n",
    "\n",
    "Vanwege de inklemming bij B kunnen twee randvoorwaarden direct worden geformuleerd:\n",
    "\n",
    "```{math}\n",
    ":label: BC_3\n",
    "\\varphi_{\\rm{BC}} \\left(4\\right) = 0\n",
    "```\n",
    "\n",
    "```{math}\n",
    ":label: BC_4\n",
    "w_{\\rm{BC}} \\left(4\\right) = 0\n",
    "```\n",
    "\n",
    "::::::"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "3. Er zijn nog meer vergelijkingen nodig om alle onbekende integratieconstantes op te lossen, deze kunnen we vinden door de overgangen tussen de verschillende elementen te analyseren waarbij zowel verplaatsingen als krachten in overeenstemming moeten zijn.\n",
    "\n",
    "::::::{prf:example}\n",
    ":nonumber: true\n",
    ":label: df_example_4\n",
    "\n",
    "Bij C kan een randvoorwaarde worden geïdentificeerd uit de consistentie van verplaatsingen: zowel de neerwaartse als rotatieverplaatsing moeten hetzelfde zijn als de neerwaartse en rotatieverplaatsing aan de linkerkant van element BE:\n",
    "\n",
    "```{math}\n",
    ":label: BC_5\n",
    "w_{\\rm{AC}} \\left(4\\right) = w_{\\rm{BC}} \\left(0\\right)\n",
    "```\n",
    "\n",
    "```{math}\n",
    ":label: BC_6\n",
    "\\varphi_{\\rm{AC}} \\left(4\\right) = \\varphi_{\\rm{BC}} \\left(0 \\right)\n",
    "```\n",
    "\n",
    "Ten slotte moet er ook consistentie zijn met de snedekrachten, zoals getoond in het vrijlichaamdiagram:\n",
    "\n",
    "```{figure} intro_data/FBD_C.svg\n",
    ":align: center\n",
    "\n",
    "Vrijlichaamdiagram knoop C\n",
    "```\n",
    "\n",
    "Verticaal evenwicht geeft:\n",
    "\n",
    "```{math}\n",
    ":label: BC_7\n",
    "\\sum {{{\\left. {{F_{\\rm{v}}}} \\right|}^{\\rm{C}}} = 0}  \\to -{V_{\\rm{AC}}}\\left(4\\right) + V_{\\rm{BC}}\\left(0\\right) = 0\n",
    "```\n",
    "\n",
    "Momentevenwicht geeft:\n",
    "\n",
    "```{math}\n",
    ":label: BC_8\n",
    "\\sum {{{\\left. {{T}} \\right|}_{\\rm{C}}^{\\rm{C}}} = 0}  \\to -{M_{{\\rm{AC}}}}\\left(4\\right) + M_{\\rm{BC}}\\left(0\\right) = 0\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "4. Er zijn nu evenveel vergelijkingen als onbekenden opgesteld. Deze kunnen nu worden opgelost en ingevuld in de opgestelde vergelijkingen voor het gewenste resultaat. Omdat het aantal onbekenden bij statisch onbepaalde constructies als snel uit de hand loopt wordt wordt aangeraden om je rekenmachine of een symbolische wiskundetool te gebruiken om dit op te lossen. Om het stelsel vergelijkingen op te lossen met behulp van numerieke hulpmiddelen, kun je de vergelijkingen opschrijven als een matrixformulering $Ax=b$ die eenvoudig kan worden opgelost door bijvoorbeeld je grafische rekenmachine of in python. Als alternatief kan een symbolische programmeertaal worden gebruikt om de vergelijkingen exact op te lossen. Als voorbeeld wordt het [`sympy`](https://docs.sympy.org/latest/index.html) pakket in `Python` gebruikt:\n",
    "\n",
    "\n",
    "::::::{codex-start}\n",
    ":nonumber: true\n",
    ":label: df_example_5\n",
    "\n",
    "::::::\n",
    "\n",
    "Er zijn 8 randvoorwaarden geformuleerd met 8 integratieconstanten. Deze kunnen nu worden opgelost. \n",
    "\n",
    "Het invullen van de randvoorwaarden geeft:\n",
    "\n",
    "$$\n",
    "\\begin{align}\n",
    "C_{4} = 0 \\\\\n",
    "C_{2} = 0 \\\\\n",
    "- \\cfrac{4 C_{5}}{1875} - \\cfrac{C_{6}}{625} - 4 C_{7} + C_{8} + \\cfrac{8}{375} = 0 \\\\\n",
    "\\cfrac{C_{5}}{625} + \\cfrac{C_{6}}{1250} + C_{7} - \\cfrac{8}{375} = 0 \\\\\n",
    "- \\cfrac{4 C_{1}}{1875} - \\cfrac{C_{2}}{625} - 4 C_{3} + C_{4} - C_{8} = 0 \\\\\n",
    "\\cfrac{C_{1}}{625} + \\cfrac{C_{2}}{1250} + C_{3} - C_{7} = 0 \\\\\n",
    "4 C_{1} + C_{2} - C_{6} = 0 \\\\\n",
    "C_{1} - C_{5} = 0 \\\\\n",
    "\\end{align}\n",
    "$$\n",
    "\n",
    "Het oplossen hiervan geeft:\n",
    "\n",
    "$$\n",
    "\\begin{align}\n",
    "C_{1} = \\cfrac{35}{8} = 4.375 \\\\\n",
    "C_{2} = 0 \\\\\n",
    "C_{3} = - \\cfrac{1}{150} \\approx -0.006667 \\\\\n",
    "C_{4} = 0 \\\\\n",
    "C_{5} = \\cfrac{35}{8} = 4.375\\\\\n",
    "C_{6} = \\cfrac{35}{2} = 17.5 \\\\\n",
    "C_{7} = \\cfrac{1}{3000} \\approx 0.0003333\\\\\n",
    "C_{8} = \\cfrac{13}{750} \\approx 0.01733\n",
    "\\end{align}\n",
    "$$\n",
    "\n",
    "Nu is de volledige krachtenverdeling en verplaatsingen van de constructie opgelost. Bijvoorbeeld de momentenverdeling in BC die kan worden uitgezet als:\n",
    "\n",
    "```{figure} intro_data/Momentenlijn.svg\n",
    ":align: center\n",
    "\n",
    "Momentenlijn BC\n",
    "```\n",
    "\n",
    "De krachtsverdeling en verplaatsingen kunnen no ook \n",
    "\n",
    "**Oplossen met digitale hulpmiddelen**\n",
    "Om het stelsel vergelijkingen op te lossen met behulp van numerieke hulpmiddelen, kun je de vergelijkingen opschrijven als een matrixformulering $Ax=b$ die eenvoudig kan worden opgelost door bijvoorbeeld je grafische rekenmachine of in python.\n",
    "\n",
    "\n",
    "$$\\left[\\begin{array}{cccccccc}0 & 0 & 0 & 1 & 0 & 0 & 0 & 0\\\\0 & 1 & 0 & 0 & 0 & 0 & 0 & 0\\\\0 & 0 & 0 & 0 & - \\frac{4}{1875} & - \\frac{1}{625} & -4 & 1\\\\0 & 0 & 0 & 0 & \\frac{1}{625} & \\frac{1}{1250} & 1 & 0\\\\- \\frac{4}{1875} & - \\frac{1}{625} & -4 & 1 & 0 & 0 & 0 & -1\\\\\\frac{1}{625} & \\frac{1}{1250} & 1 & 0 & 0 & 0 & -1 & 0\\\\4 & 1 & 0 & 0 & 0 & -1 & 0 & 0\\\\1 & 0 & 0 & 0 & -1 & 0 & 0 & 0\\end{array}\\right] \\left[ \\begin{array}{cccccccc} C_1\\\\C_2\\\\C_3\\\\C_4\\\\C_5\\\\C_6\\\\C_7\\\\C_8 \\end{array} \\right] = \\left[\\begin{matrix}0\\\\0\\\\- \\frac{8}{375}\\\\\\frac{8}{375}\\\\0\\\\0\\\\0\\\\0\\end{matrix}\\right]$$\n",
    "\n",
    "[Hier](rref_ti_84) is uitgelegd hoe je dit op een grafische rekenmachine kan oplossen.\n",
    "\n",
    "Als alternatief kan een symbolische programmeertaal worden gebruikt om de vergelijkingen exact op te lossen. Als voorbeeld wordt het [`sympy`](https://docs.sympy.org/latest/index.html) pakket in `Python` gebruikt.\n",
    "\n",
    "Beide opties zijn hieronder getoond"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "tags": [
     "remove-cell"
    ]
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAD0AAAASCAYAAADomNYKAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAASdAAAEnQB3mYfeAAAAs9JREFUeJzl102o1FUYx/HPaEh2q3sxCgKtqNRdS7tgQi8g0SZXLgLBCiEsYnITFPTwBLqMUBRq0S2iVhEGtRF6WVTULktLbwTtgsDK6nbpzXFxztTcaeY2M3eGS/SF4YHznJfn9z/nPOeZRqvV0ovM3IwHsBPXYQY/4CRex4sR8UvPwatEZm7E07gbV+EbHEdGxPftfo1u0ZnZQOJxrMOH+BTncb3yETbg/YjYMWkhg5KZNymxXoM3cAbbcAfOYntEnINLeox/AXvxOe6LiJNdk6/HAdw8ofhH5Zgi+NGIONJuzMxn8BgO4iFY0zkqMw8ogr/AbLdgiIjFiDiI/ZOKflgy80blBH6No13uwAL2ZOYUHaIz81rla/yB3RHx03ILRcTi+MJeMXdWeyIiLnQ6qo4PcBlmWXq8m7gUcxFxatJRZmZTSY6D8klEHO/j21rtfB//l8pJ2IK3O0XvqvblIQJZCU0lMQ7KS0om7sV0tef7+NvtM9SdzszLla/QwsdDBPIXmfmEcj2ORsQj/9Y/Im4YZZ0RaVTb4u87fXW1P47y9mbmLPYpT9tq0N7J6T7+Kzv7tY/3QrXrM3NtRPw56GqZOY1X8CCeGmJc0/ju9Nlqt/Txb652nrrTEfGtku7X4fblVs7MNV1Nz+O1iHhnuXE9aCrPyaC/XcvM9W61O7vjy8wrsB2L+Iil2fvZ+juWmfdGxJmuwQ2lvLsfu2vbPqVI2TOw1Mo473REfJWZJ5QM/TCOdLgTU3guIhZYKvowblHq7c/qJPP4HZtwGzbiVcjMrTiEHRHx27gErID9Shl6ODPvUgqsW5UydB5Ptjv2qr3vUZLSNiXB/awU7qfxFt6MiHOZuRdz6Lz/a5UMeQFTEfHrBMT1JTM36f+H47t2v3+IHmKBGWXnO5lTCoFDOB0Ro00+YUYW3YvMfA+nBnmnV5PuTPy/YKw7/V/hImUh8qiSnFTZAAAAAElFTkSuQmCC",
      "text/latex": [
       "$\\displaystyle C_{4} = 0$"
      ],
      "text/plain": [
       "C₄ = 0"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAD0AAAASCAYAAADomNYKAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAASdAAAEnQB3mYfeAAAAt9JREFUeJzV10+IVXUUB/DPTBGZmWIYBFpSaruWJZmRBSZtkhYGkZBBUAoxSRDS4nACV1HISANt+kPQMoxqIw1tKmqn/TGdCNwFQX+smYbIei1+99ad17vje773kL5wObx7fud3zvf3O+e8cyc6nY5eyMzNeAw7cQPW4GecxNt4PSJ+62l8iZCZ6/E8duFafIdjyIj4qV430U06MyeQeBZX4BN8jnO4UTmEtfgoIraPm0i/yMyblVivwzs4jduwA2ewLSJ+gMt72L+KR3EKD0fEya7NV+AgNo0p/ovFjEL4qYg4Wr/MzJfwNA7jCZhsWmXmQYXw19jaTRgiYjEiDmP/uKIfFJl5k5KBZ/FylzqwgL2ZuZIG6cy8XjmN89gTEb8u5ygiFkcX9tC4p5LHI+KvpqLi8TGuwlaWpvcUrsRrEfHluKPMzCmlOfaLExFxrEV3SyXnWvTfKJmwBbNN0rsr+eYAgQyDKaUx9os3lE7cC6srea5FX79fQ3XTmXm1cgodfNZvFJl5CA8qJ/07PsWhfjIlIjb262cEmKhkh39rel0lfxnwv/dupWveodTVeXyQmWuHj3Mg1De5ukV/TXNdnd4LlVyRmZdFxJ/9eIqI+5q/M3NvtfE2vLuc7Yhr+kwlt7ToN1dyjop0RHyfmWexUbm92TbPmTnZ3SEbWKVkz49t9g1MGV1Nf1jJnd3xZeYq5RIWlfJb0r2PVM9MZj4QEaebu1aT2i7sw54W50dwot58OYyypiPi28w8rnToAzjaUCdW4pWIWGAp6WncqszbX1SbzOEPbMCdWI+3ejnOzBdwF7b3Wx4jxn5lDJ3OzHuVAet2ZQydw3P1wl6z9/14XJlb12FeGdy/wvt4r55hGzYv4hHsiIhT4+F0YWTmBu0fHP+U3H9IX4SjaTzkEhMeBEORzswZ5YZ3Kx8oNeYjYn640MaHyQsvWRZPKh17Vkml+nlmyH3HiqHT+/+IvwHtQ/pzJbnwjAAAAABJRU5ErkJggg==",
      "text/latex": [
       "$\\displaystyle C_{2} = 0$"
      ],
      "text/plain": [
       "C₂ = 0"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAa8AAAAUCAYAAAA6EbIDAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAASdAAAEnQB3mYfeAAADP1JREFUeJztnXewVdUVxn+IscRGorEkGhUjUWwwo4AoiILEGt+zTWIsWFBHDXmWoViy+HRQnKggKmosoEaNiQoGxTqixJ4RsaFi4VkRFMVgj0L+WPvA4XLOLefe++59yf1m3mzO2eWsu9faa++91tqbDkuWLKFakHQEcGN4HGxm1xYovyVwDDAQ+DnQCVgIPA/cCUw0sy8z0LEB8AFwhZkNCe8OBnYDugE7AGsBN5vZ4QXa2hf4A9AVWBeYCzwLXGJmT+aUHQRMKEDeYjPrGKvTCmyaUnaemW1YoL2qoy35FMvrA7QAvYEfA58ALwJjzWxqKLMu0AzsC2wH/Az4NpSbAEwws8UJ322l/vu8LsZSW6G90l+Kfsipl6SjLgR2BLoA6wFfAW8Dk4HLzWxBrP4gStQ1oV4rdSj7kjYGzgX2Ylk/TgZkZp8CrFzFj28CXAZ8DqxZoGwHQMAwYBXgCeB24DO8YwcC/YHfAn0ykHMAsBIwKfbubHzS+hx4D9iqUCNBmIYCC/CO/Bj4RWj/IElHmtlfYlVmht+VhD7AHsC9CXmfAWMT3n9eiMZqokZ8QtLZwHl4f9+NC/J6QHegHzA1FD0EuDLkTwPeATYADgSuBfaWdIiZJa3Y6rLPoe7GUlXRnunPoB/iSJL9U4EZwIPAfGANoBcwEjheUi8zezeUnUk2XQN1JvuStsD5vj5wF/Aq0ANfFOwlaRczW1CVySsI4ASciXcCZxSocj0wCJgFHGZmz+e0tzpwGi4IWdAcaJkee3cqPmm9ge/ApuVrQNKG+O+YB2xvZvNjebsDD+MrhaXCaWYzcaFKai9ahf05IXuhmY3MR0+N0OZ8knQIPnE9BBxoZotyvvmD2ONs4NfAPfEdlqQzgWeAg/CJ7I6Eb9dln9fhWMqL2A5gdzN7JEMTNaO/HNqz6IccJOmotc3s64RvjQLOBEYAJ0FZugbqT/bH4xPXEDO7LHop6RJcb48CTqzWzmsIPtP3C2kqJJ2GC+srQK9c5QRgZl8Bo4LglgRJ6wQabjGz72NtTouVKaapTfGV0dNxwYzakrQI+EmRNG2Lr6DeB+4ppk6tUQs+SVoJuBD4EldkSd/8T+zfDye1bWYfSroKF/p+JE9e9Yq6GUvVRjunP7N+yKOjVpi4Av6GT15bFiKqvekaSZ3x3XUrcEVOtgHHA0dIOr3ik5ekrYHRwKVmNl1S6oCTtBGuUL4DDk0S1jiC4JaKfXHzw50Z6sbxOu4/6SFpPTP7OMqQ1Bf3mU0usq0TQnpdXFhjWFXS4bit/wvgBWB6Stmqo4Z86g1sjpuNPg3+hG2Br4Fn8vkQEhBNct+l5NdVn0NdjqWqob3TT3n6oVQdtX9IXyiibCFdA/Ul+5GMP5DrnzazRZIexye3XhWdvCStDNyE+xrOLKJKC7Aa7kh/qZK0xNCMM+TBchoxs08kDQMuAWZJmoxv87fATVUPskxQUhFWjIcDi3E/TBI2xPsxjjmSjjazR7P9grLQQm34tFNI5+G2/+3iFSRNBw42s4/yNRzk8sjweF9Ksbrq8zodS9VEC+2Y/jL1Q14dJekM3Ne5Dh7AsSs+wYzOR1ORugbKlH1JLXhATbGYaWaTU/J+GdLZKfmv45NXl5VK+GAx+CPuRB9U5MqoKaS5HVcRSFoNj1a5N88WvGiY2VjcZ7IyMBgYjgcJvItHP81Pr70Uh+KMvjfmbI1jAu6Q3hB30G4HXA1sBtwraYeyfkQ2NIW0rfm0fkhPBFYHBuAr2G2B+4G+wN+L+MToUGeqmd2fkF+PfV5XY6kN0BTS9kp/Jv1QpI46AzeZteAT133AwEKLNgrrGqiM7LcE+or9a8rT1joh/SwlP3rfabmdV4GwySQsDS2X1ANfIV5cjDlH0pp4COgS4OkSvomkkXgnxJEU1jkQX7FMogKQNBQ4HxgHXA58iEcpXgDcLKmbmQ0t0MzxIb06KdPMch1wLwEnSvocOB2PNGouQGcrGfmY0FZmPoX6G+GTxz44L94ATslZ0aXxKQrr7YDvsCLn/cuSmvHV2W6Sdk6TOUlD8H57FTgiqUwd9nlbjqVWkukeb2YnZ6wLMC3Bl3yDmQ1KaKcc+jvi/Dkc2AiPNL0ZGGlmaSbiitGe02YW/VBQR0V6LYTT98bH03OS9jOzGXlIyqtrQttly76ZbZYvv8LoENIluWbDN3F/QrH4AJYzccwGzimybuS8/HfG8xqv4U7sCEn22WbcDl22o1JSPzx4YJKZnRbLmhFTpKdLusrM3kppoysufO+xLLy7WFyFC1PfIspm4mMKMvNJUifgceAx3K7/EdAZD/uNI41Pn4b0rdyoMzP7StL9wLF4GO0KSl7SycCleORafzP7pBT6qUGf12As7cSyRQL4DvVBitvRjmVFc1E3POz7BtzpHsfMlHbKoX8YcDJwFH6eb3tgIvANHqWahrFUhnagLP1QtI4ys3nAJEkzQns34vxKoqccXQOlyX4lEe2s1knJXzsqt9zkZWb9M34wWjkBfJ0SvXeNpGtw53MLbuMFWF1SxwzOwe/M7MO0zLAi2x942MzStqClYL+QrhBSb2ZfSnoGF8TuQOLkRXHO0zRECn+NQgXL4GMSyuHTUGCumR0ZezcnXqAAn14L6cKU9qPJbYXIs2CHH4OvJPsXadLNRS36vE3HUq7pSdJwfCIu6OsIZrLlEMLND8DNZI8USUY5MtYb+IeZTQnPrZKmAD3zVaog7RFK1g9ZdZSZvS1pFtAtNzgkhnJ0DZQg+xX2eUVjvktKfhRhObtSPq9vgOtS/p4LZR4Lz08CBGXSikfZ9MvXeAiZzkVnSe9LmiPpVkmb5eT3xU9mV8RkCKwa0rRw+Oj9t0mZwbZ9BO48vS7D93cOadrEWBWUyacm4ClJt0iaL2mmpFPkZ5ci5OPTdDz6bEtJqyTkR6vO1hw6huET10z8zE6WiQtq0+e1GEtR3iq4+e16Sz7IXRWUSf8TwB6Stgr5XfGItSy7jXKQRT+Uo6N+GtIVJqYK6BooTfZbqJzPK5r8B+byWtJawC74TSNPVSTaMDiUj0vKC/6p7rjNODfiZWz4Gy/pADN7NaduB9yZeTTufIzwNM6c2fgNCiOAJyVtEzMNHYgz767MP2x5/BM4BT/ZfrWZvR+jc2+8U7/GB1MSDgF+BNyd5jyVtA2+U/kk5/2muA0d0g85VhNjycanzrhJZwzwK9wscxnu24jOcKTyycw+lnQb8Ds8gOHs2Df3DG1+RiyCUNI5+GHQZ3Gndl5TYb31eQ3GUhxN+Ap6Yhbay8RYstE/Gg/imSXpezxYYpSZjW8LomPIoh9SZT9MxgtzrUtBoZ+HBzM9YeGqpBwU1DWhrYrIfiV9Xmb2pqQHcF/gybi+WEoavhO82sy+qNr1UEViHG6jPgZ4MRA9Gz+TswkeWbMxcEu8kpnFrzl5UX6C/C38gOMl4X0Tztx5SR+W1MSyFUAU6LGzpInh3x+bWfw2g9vxWx4GAK9ImoQ7ZLfGTQYdgOEWu28sB5HzNO2UO7jQDZc0DTevLcJDbffFw4inAhflqV8tZOITfmjzWTMbEZ6fk9QFF8po8moiD5/w2xR6AmfJz8s8gzvam/FV52AzWwgg6Sh84voeVyZDEsxurWY2MfZcr31eKrLyKI5j8ci0fD7QaiEr/Qfji5vDgJfxBdKlkuaYWdZdRxZk0Q9NpMv+XsCf5MdB3sTD7jfAbwPqHNoenEJLMboG6lf2T8In+XGS+uOH1nsCu+MycRZU8W7DYhBME8dKugNnRA9gT/xOrbm4WeQe/D67fO18EWzAWwJI2gkX9DF5qnXDnbxxdA5/4BdgLp28zGyxpH1wxfsbXHn+EL8gdiowzsweSPqQ/LDprhR2nk7Dzzl0x7fta+D+nsdwJ/5NbWnOiVAGn+biCiWOWfitEUXxyczmS+qJ77qa8dsCFoXvXWBmT8WKbx7SjrgpIwmPsvzOoi77vFSUO5bCansAvhtoc5RB/8XARWb21/D8YvgtZ5LdZFYyStUPRcj+Q/jkswt+B2sn3Dc4G5fLcUlWhRJ0DdSp7Ifd144su5h3H1wGxgGKfneHat4q31YINt45wJVmdq6k83FTYmczm5O/dgPVgqRbgE3MrE/s3XnAQWbWtcGn+kEwSZ6A8ys1xLzeIGkBYGZ2eezdCOA4M9uidpTlR0P2y0etzYaZIOkiYAp++8D6eEjxGniIK/iq5/mGUNQcY4AnJJ0F3Iav8Iaw7MaIBp/qAMGPcjTuS2s3E1fAFNz0NQff5XfHTc035q1VezRkv0y0y8kL327fiv+3GB8BT+EXeb4NYGZb15C2BgLM7F/Bt3g+vsB4J6TjQ36DT/WBAfi9dtfXmpAM+D0ewBDdRD4XuAY3OdUtGrJfPv4nzIYNNNBAAw38f6HSdxs20EADDTTQQNXRmLwaaKCBBhpod/gv8eISpCWMfDQAAAAASUVORK5CYII=",
      "text/latex": [
       "$\\displaystyle - \\frac{4 C_{5}}{1875} - \\frac{C_{6}}{625} - 4 C_{7} + C_{8} + \\frac{8}{375} = 0$"
      ],
      "text/plain": [
       "  4⋅C₅   C₆                 8     \n",
       "- ──── - ─── - 4⋅C₇ + C₈ + ─── = 0\n",
       "  1875   625               375    "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVUAAAAUCAYAAAA9f65TAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAASdAAAEnQB3mYfeAAAC45JREFUeJztnHmwFcUVxn8sUYkLJBLARKNixF2xygU1IIrigonvqVhZRDGWS4lFXpBC3HL4tNwqFr7gSioqxmhiooKlgoqRSFypEnFDRQVcEUTEIGIZhfxxemAYZu6be+/c66vkflWv+t3u6e4zp0+f7j7n9HRYs2YNaZC0I/ArYDDwQ6AbsBx4AbgXmGRmn6dWLgFJPYEPgOvNbGSirD/QAhwIfBdYBrwEtJrZ1PDMlkAzMATYA/gB8GV47lbgVjNbndLvQmDbDLIWm1mvct+lSNSK36HtDXgu6QTgYKAvsBewOXCHmZ2UUr+uPJe0NXAJcCSwJbAImALIzD7J8861Qi3HqZaQNAT4NbAr63j6HDDezJ4uUW892ZF0FbAP0AfoDqwC3sbH5zoz+zhRfzguI6Ww2sw6JeotpB3O1zyy2TmlUgdAwHnARsBTwN3Ap/hLDgYGAT8H+ldA17FAR2Byot+LgEuBpcADgdjuwN7AQGBqeHQocGMonwG8A/QEjgP+CBwlaaiZpa0WnwKtKfmfVfAehaAO/IZ0nl+EK9PPgPeAnUvUrxvPJe2A86AHcB/wGrAfrhCOlHRQcuLWA3Uap5ogKMIxwMe4AlgK/AiXi+MlnWxmf86onpSd3wCzgenAEmBToB8wDjhDUj8zezdWfw7OtzT0Bw4FpmWUt6v5mlc2N1CqwC3AcGAu8AszeyHRcBdgFD4olaAZH9yZsTaH4gr1UeA4M1uR6PNbsZ/zgJ8CD8Z3R5IuAGYBx+OT/Z6Uvpeb2bgK6U5FbCU+xMz+WUETteY3pPAcnxzvAW/iO9YZJerXk+c34EI70syujfU1PtB8GXBWGe0VMUZQn3EqHJJ6AaOBxcCeZrYkVnYI8Bi+88pSqknZ2cLMvkjp5zLgAuB84Owo38zm4Io1jbZoh/yHjL4Ln69VIpdsdozXkDQKF5xXgX5JwQEws1VmdhkxxuWFpK74ynS/mX0d8joCVwGf48K6IlnPzP4T+/8xM7s/edw0sw+Bm8LPgeXS9k2g1vwOfWzA89DuDDN7I2N3maShLjyX1Bvf8S0Erk+SAawEhknatNq+yqSr5uNUQ2yL7zSfjStUcBkAVgDfS6uYJjtpCjXgbyHdMQ9RknbHd7jvAw/mqfNNohzZ7ByrtBWuab8CTkxTbuu1YraqAtqG4Eene2N5BwLb40epT4LtZ3fgC2BWKXtPCiLl+1VG+caSTsJtYSuBF4GZcWVTL9SJ35DO8yJRJM8PDekjKQp8haQnccHuB/yjaspzoI7jVCu8gdu/95PU3cyWRgWSBuC29CkZdcuRnZ+E9MWcdJ0Z0ptLzL92M18pQzbjx/8WYBPc6fByjQhrxpkzPZa3b0gX47aaPeIVJM0ETjCzj0o1LKkzcHL4+VDGY72A2xN5CySdamaPt01+oWih9vyGdJ4XghrwfKeQzsto6w1ccPtQJ6VK/capJjCzZZLOA8YDcyVNwY/zO+AmnemsU3BJZMqOpNHAZkBX3HH1Y1zpXdkWTcFUchKwGrfJZ6Gq+SqpBXck5sUcM5uSUZZbNuPH/6aQJl+iEEjaBPeYTUscIXqE9CygC3AYvnruDjwMDAD+nqOLK0OdqWb2cEr5rbgjoRduXN8DmAhsB0yTtFeZr1QtmkJaE35DSZ4XhaJ53jWkn2b0F+V3q5zkstEU0pqNU61hZq24zbszcDowFnc+votHKyxJ1skhO6PxY28LrlAfAga3tfkJOBEfw2kJp1YcRczXlkBj3r+mEm3lls3OAJI2w1f/NcCzOYgl1BsXiIkjK9xhML6yTU7kR6EUHfAdaWSvekVSM74yHCzpgCxTgKSRwLm4N25Y2jNmlvRAvgycJemzUHccvjJnoo0wjxnSBk7O28xseEo7FfE7Vn8rXKEdjfP0TeCclNU7i+dVo148T6BDSDPtwEWNUWir4nEqQccNZjaignpZSA2DS7Q5BrgcmABcB3yIR3tcAdwhqa+ZjUlUKyk70RwPIVcH4vL4vKRjzGx2GzSfEdKJWQ8UITtmtl0bdBSJtbIZHf8jQ/W/K4ixe531nRRZ9o5m3LaTNEpHcYfzkw4AM1sl6WHgNDx0YQOlKmkE8HvcKzvIzJaVSf9N+CANyPFsKxvukvriYSe34UbsOOZktFMxvyV1A54EnsBtXh8BvfHwliSyeF4VasjzaLXvSjq2SDyXhlaKGSOobl7sy7oNA/iOfjptn7rewv0JefFBqUJJA3FH8GQzGxUrmh3btJwr6SYzmx8rzyU7ZrYYmCxpdmjrT/i7ZtGzK66E32NdmGQ5KGe+Foncshkp1ZUh7SKpU5mG4K+CFzgTkjrhhuzHzCw5IV4P6fKM6pHS7ZLSbgtwDb6KDUo7xuRAVKdNj3I4RiVpGI5P2EllhOtUw+8xwCIzOzmWtyCFrlI8rxg15nkkC30y6kWe5Sy7VpFjBFWMU/IYLGksrjBL2gLNbFAZ9OXBMSHdIGTOzD6XNAtXoHsD86Ey2TGztyXNBfomHWIJ5HFQlULu+VqwTTW3bHYECBNjIe7pG1iq1xACFUdvSe9LWiDpL5K2S6k2AL99kHaUmIl7VneUtFFKebTqLUzQcR4+uefg8YeVTG6AA0I6v+RTBaJKfjcBz0i6U9ISSXMknROC0+MoxfOKUAeeRxN/cPK9JW0OHITf4Hmmwn7LQpXjFC/bCHfM3JInhK1gbBzS1LCpWP6XsbxKZef7IU1VlsFOOwx3UN1cZtsRypmvLRRnU80tm3Hvf2v4u0HSsWb2WqJiB9xwfSpuaAa3Mw3Ddw498cDfpyXtljgSHocz8r4kpWa2VNJdwC+B3+I3faI+DweOwLfeD8XyL8YDlp/DjeMlj5+SdsN3d8sS+dviNibIDn6uFVopn9/gR/0RuHI7Aj/aXovb/eLxc5k8rwT14LmZvSXpEdyeNwJ/r7VV8d3JRDNbSf3QSmXjFEcTvmOaVCMaS+FfwDn4baeJZvZ+VCDpKFwZfIHfFIqQKjuSdsYD8j9M5HfEL+/0AJ6y7KvEQ4HvAA+UcFAVNl+LtKmWI5txpToB2BO/1/xSaGAeHoe4De7h2xq4M9ZR/HrZS/IbEvPxQOnxsbImnNmLM2geBewPXBhi52bhxvpmfNU73cyWA0g6BZ/cX+MCMzLF+bDQzCbFfg8FxkqagR+VV+AhJUPwcJmpwNUZtNUKZfM7oCPwnJmdH34/L6kPPtBxpdpECZ5LamLdyhw5Fg+QNCn8v9TMRodn68nzs/EJPkHSIDzgfn/gEJw/F6a9Tw1R6TjFcRru6S5p/6wR7sZvKh4GvCppMu6o2gU3DXQAxtr6V3+bSJedI4HfhTDHt/DQrJ74jbzeod3TS9ASOaiyblBFaI/zFXLK5lqlGo4lp0m6B2fMfsDh+D3bRbiT6EH8Xn4qzGxlsKusvVUhaV9c6K4pUW+JpP3xXWozHty9IvR3hZnFj3vbh7QTvr1Pw+OsvyuYgceZ7Y0fHzbFbbhP4KEyt9f7WFYFvxcBryTy5gJrP06Th+f4DveURF7v8Af+kYzR4f+68TzsCPZh3UcrjsbfeQL+0YpynWJVodp5EXZXh+G7v7rDzFZLOhpfdH+Gz69v4x8rmgpMMLNHoufbkJ1HcYV4EP7diG643XkePqYTssZH0i74ApTHQdXu5ivkl80OWV+pqgTBZrIAuNHMLgl5l+Nmgd5mtoFDpYHyIOlOYBsz6x/LuxQ43sx2Db8bPG8nCGGHZ+JjlnXrrN2gITvVI+2DKrkh6WrgfvyrRT2Ai/FV5bbYY83AC40BKgzXAE9JuhC4C1/NR+Ifs4jQ4Hk7QLA1norHwrZ7hRrQkJ0qUdVOVdJfcU9hdzxe8hngYjObWwx5DaRB/n2Ey/Ej0ju48f7ab+JI1EA2JA3GbwXuZGaZYWAN/G+h0ON/Aw000MD/OzJj6xpooIEGGigfDaXaQAMNNFAg/gtSyfhJ0+EDjAAAAABJRU5ErkJggg==",
      "text/latex": [
       "$\\displaystyle \\frac{C_{5}}{625} + \\frac{C_{6}}{1250} + C_{7} - \\frac{8}{375} = 0$"
      ],
      "text/plain": [
       "C₅     C₆          8     \n",
       "─── + ──── + C₇ - ─── = 0\n",
       "625   1250        375    "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAAUCAYAAACajtrtAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAASdAAAEnQB3mYfeAAAC9JJREFUeJztnHu01UUVxz9XTDE0KE1xpYGSDwgVlspDA0EUH1ReVCwNFJ9hGt20RFDbfDUVl68rKomPAE1XrVAwFHyCkuGjNBRFxQc3X4ghQoCiCfTHngM/Dr/feZ97Dq3zXeuuub+Z+c3sM/s3s2c/ZurWrVtHuSBpCHBneDzTzG7PUn8P4DSgP/BtoA2wDHgRuA+YaGafFkDHTsAHwM1mNjzkHQ8cAnQB9gO2A+42s8FZ2hoA/ALoBGwPLAKeB64zs6fT6g4FJmQhb62ZtYi80wS0S6i72MzaZmmv7GhOPkXKegENwEHAN4ClwDyg0cymhzrbAwOBAcA+wLeAL0K9CcAEM1sb028T1T/mVTGXmgs1+isDSbsAlwJHsmF9mwrIzD5J1duyjATsCtwIrAS2zVK3DhAwAtgKmANMBpbjE7o/0A84EehVADnHAFsAUyJ5F+MCYyXwHrB3tkYkXQVcAHyMD+YS4Duh/eMknWxmf4i8Mjf8rjj0Ag4FZsSULQcaY/JXZqOxnKgQn5B0MXAZPt4P4B/zDkBXoA8wPVQdBPwulM8C3gF2Ao4FbgeOkjTIzOJ2S1U55lB1c6msqNFfOUjqgNO7I3A/8BrQDd8oHynpYDP7GMokOMLgTcAX2PuAX2V55ffAUGA+cJKZvZjW3jbAefgiXQgGBlpmR/J+iQuMN3HNY1amBiS1xX/HYmBfM/soUtYXmIlL6vWCw8zm4sIjrr2UdnJrTPEyMxudiZ4Kodn5JGkQLjQeA441sxVpfX4l8rgA+CHwYFSzkDQKeA44Dhci98b0XZVjXoVzKSMiWnZfM3uigCYqSn8JsDnTPw4XGsPN7MZUpqTr8PXycmAYlE/jGI7vpvuENBGSzsMH+lWgR/rCAGBmnwGXh0HPC5JaBxruMbM1kTZnRerk0lQ7fDf8bFRopNqStAL4Zo40dQZ6AO8DD+byTqVRCT5J2gK4CvgUn4Rxff438v/MuLbN7ENJt+Affh/iBUe1omrmUrlRo79ykLQ7rg01ATenFRtwFjBE0vlmtqrkgkNSR2AMcIOZzZaU+LFL2hmfzF8CJ8QNdBRh0PPFAFxlvK+Ad6N4A7eXd5O0g5ktSRVI6o37SKbm2NZPQ3pHVJhFsLWkwbhtdBXwEjA7oW7ZUUE+HQTshqv6nwT/UmdgNfBcuk8pC1IC5suE8qoac6jKuVQ21OivOFLf1iPpfkAzWyHpb7hg6QE8XlLBIWlL4C7ctjwqh1cagJa40/LlUtISwUB8IXi0mEbMbKmkEcB1wHxJU3HzQQfcPPIoGwRCIsJOYzCwFre7x6EtPo5RLJR0qpk9WdgvKAoNVIZPB4Z0MfAC7vBeD0mzgePN7N+ZGg7f5cnh8aGEalU15lU6l8qJBmr05wxJDbjDPVfMNbOpGcr3CumChPI3cMGxJ/D4Fnl0nAt+gzssh+YoUetDmj5hSwJJLfHogBlmtrrY9sysEbeRbwmcCVyIO2TfxaMkPkp+ez1OwBk+w8zejSmfgDvP2gKt8MVyPNAemCFpv6J+RGGoD2lz82nHkA4DtgEOwzW7zsDDQG/gzzl0MSa8M93MHo4pr8Yxr6q51AyoD2mN/tzQgJuQcv2rj2skgtYhXZ5QnspvA2k+jixhiXFYH74qqRu+M7o2FxOCpG1x6bUOeDaPPlOmoV8D+wM7A4PMbHJM1f54FMqUmLK8IekC4ApgLHAT8CEejXUlcLekLmZ2QZZmzgrp+LhCM0t3uLwMDJO0EjgfGI3vzjPR2USBfIxpqxg+jcQF7V7A58AzwMiYHVkSn1JhynW4ZpFyNL4iaSC+OzpEUs+kb07ScHzcXgOGxNWpwjFvzrl0Dq4ptw9ZrwC/NbOsvrcsv3lWjO9wkpkNjWmnYPrT2hmFm4tuNrNzs9RtojrmSAv8+xqMr2WLgLuB0WaWZFbFzNrn008JUBfSdbCpc/wt3H6cKz6AjdTqBcAlOb6bciT/p4B45lZ4PPQEMjs6B+J+iaId0JL64I7aKWZ2XqTohcgidr6kW8zs7YQ2OuF2+/fYEEKaK27BF7HeOdQtiI8JKIZPffBIjb/jH96lwGOSOpnZ0ki9JD6l4sbfTo9OMbPPJD0MnI6HDG6ywIZF8QY8wqVfWp+5oNnHvAJz6T1cc34D59EpwFRJ+5vZS1nebWRTc0kXPKx6Eu5ojWJuQjvF0A+ApB64FSAbzSlUyxwZAZyDj/s8YF9gIr7RuizPtopBSqNonVD+tWi9jQSHmfUrsNOUxAVYnRCldJuk23BHXwNuzwbYRlKLfJyQZjaDcP4hKSIqSPIfADPNLEn9ygffD+kmYbtm9qmk5/AFsCsQKzjI7hTPhJQZrFW2ikXwMQ7F8OmI6HM4xLYcOBiYFvIy8en1kC5L6CIlWDaJUAk24Otx7aFfjmbEdFRizJt7Lt2flnWRpLOBnmRZhIPpdiOEcNxjcNPtEzmSUTD9oc/W+C79dNzElxXVMkfwjeRfzGxaeG6SNA3onumlMvg4UnNtz4TyPUK6ADy8tBT4HLgj4e+foc5T4flpgDCRm/BImj6ZGg9hmfmiN37ysSRmKmDrkCaF3Kbyv4grDHb8IbhT/I4C+u8Z0iShVBaUmE/b4d9cdOefiU+z8SiVPSRtFVPeOaRNaXSMwIXGXPw8QSFCAyoz5hWbS5JaSDoRF15zivgNeaEE9N8KTE4Kxy43iqR/DnCopL1DeSc8wimbRaKB0vo4Uhvi/uk0StoO3+x9hpubS3OOIzjvzogrkzQa34VPsk2vSWgMf+MkHWNmr6W9W4c7TU/Fncr54Fh8kU7fURWKvwLnAmdJGm9m70foPAof2NUkT7hBwNeBBxKc4kj6LrAo3aQiqR3uU4HIAcNmRCOl4VMjvpg/E8lL5JOZLZH0J+An+E7y4kifhwNH4BrMQ5H8S3CT2PNA/2zmqWob80rMJUn74EKoJbACqDezecX+ljzRSAH0SzoTP0wX679qRjRS2PiPwTdU8yWtwdfky81sXKbOSu3jMLO3JD2C+xvPwW8qSEG41j3ezFZBGa8cyRFjcZveacC8QPgCPOZ+V+B7wC7APQW0XQ/MMbPFcYWS6tkghVN3EfWUNDH8v8TMoqd0J+Onlw8DXpU0BXeOd8TNWHXAhRaO5Mcg5RSPOymewiDgQkmzgIX4JO6An3Foie9CrsnwfrlQNJ8kXY1rF73SVPl6MvAJP2XbHTeh9MZPgLfDzYJr8HubloU+TsGFxhpc0A+PMfU0mdnEyHO1jnm+KIZHr+O+iTb46fo7JfVp5rDYvOmXtBcerNLLzGI1/WZEoeN/PL4xOgkPTOgC3CBpoZkVYpkoBj/DN75jJfXDDzJ2B/riv+WiVMWKCg7zO4NOl3Qv7tjqBhyO38mzCN8FPYjfT5QzJB2IM+n6DNW64A6pKHYPfwD/InK9g5mtlXQ0Lo1/jC9cX8XNLtOBsWb2SAI9HfEPJ5tTfBYegdQVN5O0wu37T+EO07ss/p6lsqJYPkm6Fo8a6Wtmb0bys/LJzD6S1B3XNgbiB5BWhP6uNLOo9rJbSFvgqnwcnsSdjylU5Zjni2J4FBbdFF/+EfjSQILmUw4USH9P/M6ylyMbhBZAb0nDgFZm9nkV0w9wLXCNmf0xPM8L2u4oCjNpF4ygdRzAhksOj8ZpH4tfcrheK68r5+24zQFJ60gLx5V0BTAS2N3MFlaMuBqQNBb4ES405qeV1fhUhZA0E/ggKfy0WiCpDb7xiGICHiF2BfBKtQt9SR8DZmY3RfJGAmeYWYfKUZYZlTZVFYQQNx29JKy9pC7AUjN7B9+ZvlhbjCoLSeNwTaMeWCq/KBJgpZmtpManikPSGHwn/C5uaz8Jd/AOqCBZOSGYKJdF8yStwteBzeX0+TTcVLoQN1V1xc2zd2Z8q8LYLAUHcAAbh8VeHdJJ+Enbjs1PUg0xODukj6flCz/gVONT5dEWd/63xQMNXgKOsvgT9jWUHj/Hz2ukbqZdBNyGm4uqFpu9qaqGGmqooYbmRanvqqqhhhpqqOH/HDXBUUMNNdRQQ174Hyg4Ej1GSXI4AAAAAElFTkSuQmCC",
      "text/latex": [
       "$\\displaystyle - \\frac{4 C_{1}}{1875} - \\frac{C_{2}}{625} - 4 C_{3} + C_{4} - C_{8} = 0$"
      ],
      "text/plain": [
       "  4⋅C₁   C₂                      \n",
       "- ──── - ─── - 4⋅C₃ + C₄ - C₈ = 0\n",
       "  1875   625                     "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAATQAAAAUCAYAAAADWsEdAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAASdAAAEnQB3mYfeAAACUVJREFUeJztnH+wVVUVxz9PiHCQYLIUJw3E0H5o4SghGgiiqDHleyQ2GgREFkZDL7QR0lp8LdPGzBcm2UzDjyxnmtRwCsgf8ZKKyn4MP8SUVF7mQBgiBIqZSH+sfeG+wzn33XPvOYdH3e/MnX3vOWfvvc53r7PP2mutfZv27dtHHCQNAT4BjAPeDvQHdgBrgfuAxWb2cmzlCpB0LLAZuMPMZkXOjQRagbOBNwPbgfVAm5ktD9ccDbQA44HTgLcBr4brFgGLzOz1mH47gIEJYm01swFp7yVL5MV3aPsgziVdCpwLDAXeB/QFfmhmk2LqF8q5pOOBG4CLgKOBLcBSQGb2YjX3nBfyHKcicLjKX61O9Iyp2AQIuBboBawG7gF24so5DhgLXA6MrEG2S4AjgJ9E+r0e+AqwDfhZEPgtwOnAaGB5uHQi8J1wvh14FjgWmAB8D7hY0kQzi5updwJtMcd313AfmaAAviGe8+vxiWw38Bzwzgr1C+Nc0kk4B8cA9wNPAO8HPgdcJOkcM3uhgqy5oKBxyg2Hs/xpdOKgCQ1YCEwFHgeuMLO1kcaPBGYD76hRvhbgBWBVWZsT8cnsYWCCme2K9PmGsp8bgQ8Dy8qtAklfBB4FPoI/aPfG9L3DzObVKHcsJE3FrZQxZvbLGprIm2+I4Rz4PD6RPYVbau0V6hfJ+QJccWeZ2e1lfX0zyHwjMCNFe1mMERQzTnnicJa/ap3oNKFJmo3f9F+As6ITC4CZ7QFuDASkgqR+wHnA3Wa2Nxw7Avg68DJOdFyf/yn7vjKubTP7h6Q7w82NJv7h6lbIm+/Qx0Gch3bby66p2EZRnEsajFsKHcAd0e6ATwGTJV1tZi/V01dKuXIfpzxxOMufVid6llU8DlfM14DL4m66U0tOQFqMx83d+8qOnQ2ciJu/L0oaD5wKvAI8ama/TdF+aeJ7LeH8GyVNwn0HLwHrgFXlD3pRKIhviOc8S2TJ+XmhfDDqkzOzXZJ+gyv3WcAv6pa8ChQ4TrngcJeflDpRbqG1Ar1xB+9jOQnXgiv1Q2XHhoVyK/Bn3Om8H5JWAZea2T8rNSypJ/Dx8PPnCZcNAO6KHNskaZqZPdK1+Jmilfz5hnjOM0EOnJ8Syo0Jbf0VV96TKWhCo7hxygutFCi/pFY80FAt1pjZ0grnU+lE+YTWHMqo8mUCSb3xCMUKM3ul7NQxoZwBbALOB36POypvBS4EfowvaSrhZtyyW25mD8ScXwT8CtgA7AIGA5/FTdYVkkZE/Qo5ozmUufANFTnPCllz3i+UOxP6Kx3vX6fcadAcytzGKWc0h7Io+VtJjmzHYQkerUxCKp3oCSDpKHyG24dPJlVB0ijgC8AZwHHARDO7J+HyccBRRKKbQI9QNuGWWEnBN0hqwWfmc4Pyxy4/Jc0CrsajH5PjrjGzqKPoMWCGpN2h7jzcmklEF2kI7TG+qCVmNjWmnZr4DnXn4g74U4B/A78D5ia8fZM4rxtFcR5BUyjjc43IboxCW/WM00zg08CgcGgD8FUzW9ZFvQ7STQixqTahrXrkT5JjgZnNTKpnZoPS9JMBOulEyUJ7ayj/lTIHpQ+ev7KIrh3CLXjuUnRASzkkz0QtJDPbI+kBYDoepj1oQguK8y08ejPWzLankB/gTvzhGlXFtW0cbB0MxdMiluCOy3KsSWinVr7BLdUFwB/wwbwBeFjSu2PuPYnzupAj56W3bT/i8abIdXFoI5sxgvrG6TlgDr4kagKmAEslnWFm6yrUexr3H1eLzRXO1SP/MA4YG+CW+EP4aqlIpNKJ0oRWihgdKalHtU5yM1sBrIDKkTJJPYAPASvNLKqMT4ZyR0L10oR3UPQlrNdvw9/8Y83s+WrkjqBUp09XF5pZW4wMU/GHZXGKlICa+A4yXBjpfzI+mOcAPy07XonzmpEz5yVdODmh3pBQJvlTshwjqG+c7o8cuk7SVcAIPDCSVG9sCvm6Qj3yd/JZS5qDT7YVfc05+NBS6cQRAEEpO/Bo2OhKvYc0i7QYhWf3xi19VuERmCGSesWcPzWUHRE5rsUfrDV4flEtDxa4ggE8U2P91MiY7774OEatpEqc14QCOC+lkoyL3rekvvikvQdfZueOrMZJUg9Jl+PL/9UZilgRGcrfC5gELExIni5HK55OUe2nuYv2UulEeVCgLXwWSLrEzJ6IVG7CHczTgMu6ECKKCcDreJZvJ5jZNkk/Aj4GfBnPYC/1eQEeFNhJWRRN0pfwpdafgHFdLXkkvQfYEr1O0kDg2+HnD1LeU71oIxu+2/AJJvqQJ3JeC4rg3MyelvQg7vubCdxeXhW36L5bZA4adYyTpNNwN0lvPCjSbGbrC5C5HG3Ur2fNuNW1uKvOsvahpdWJ8gltPvBefJ/X+tDIRjzP6ATgA8DxwN01yNUMrDazrQnnZwPDcbN8FJ59PhD3Ae0FrjSzHQCSpuAP1l48gjYrZrnbYWaLy35PBOZIascjqbuAk/Acrd74tqpv1HBf9aBuviXdgltiI2OWE81U4FxSMwfejqU9lSMkLQ7ft5nZNeHaIjn/DG7FzJc0Fk8GHQ6Mwfm5Lu5+ckQ94/Qk7r/rj++m+L6k0QWnf2TxXE/HI+WV/HV5omqd2D+hBVNyuqR7gStxJ/wF+J67LfibZhm+z7JqSBqGE3Zb0jVm9ryk4bh11oInye0K/d1kZuXWx4mh7IGbt3F4hM5vk3Y8Kng6vtzpg/vsfo2Hs++qwpTOFPXyLelWfBkwxsyeipzrknP8QZsSOTY4fAD+BlwTvhfGeXgjn8mBjcgfxPmYj29EThuAqAv1jJOZvYpvLQP4YxiXVuCT+Uu+X4Z69Wwgnko1oRCBY5BGJ5qS/m0jLSTtIyZtQ9LXgLnAYDPblEln/+eQNB/4KD6ZPR5zvsF5N4SklcDmpDSL7ghJ8/D0kxPMLGk3SLdB3Ob0qhHyXMo3sw6SNBTYbmbPhmMtwNrGg5UNJC3ALbNmYLuk0nJxt5mV/sGiwfkhhqSbccvn73jg5grcMT/+EIqVCsEJPw3P1ev2kxnUOaEBZ9L5XxpuCeUSfDMsZvauOvtooDOuCmV064/wRNUG590DA/CgxwA8qLUOuDhhR0V3xfn4HtyFh1qQapHZkrOBBhpo4FCjlpyyBhpooIFuicaE1kADDfzP4L+nwvob2OmQoQAAAABJRU5ErkJggg==",
      "text/latex": [
       "$\\displaystyle \\frac{C_{1}}{625} + \\frac{C_{2}}{1250} + C_{3} - C_{7} = 0$"
      ],
      "text/plain": [
       "C₁     C₂               \n",
       "─── + ──── + C₃ - C₇ = 0\n",
       "625   1250              "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAKoAAAAUCAYAAAAdtqaHAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAASdAAAEnQB3mYfeAAABSNJREFUeJztmmuoVUUYhp9jIpmaQhQIWqFp5A8pMjO8kB5QqR8eioxCRbtAqdjJbmjY62cYgVhHJSG01CQhsDAKBFG0sDQkVDKzk5L4I8HQsjSVLrsf36zarc5ae+3bWUfYDyxm77m+s5j5Zuab1VQoFEjCzKYD74S/j0tam5J3CPAIMBG4HugH/AwcBD4A1kv6LbGxnGnozwczGwAsASYD1wAngS2ASfopytc9pYKBwCrgHNA7JV8TYMALQA/gc2AzcBa4AX9xzcBDwNgq+lQXGvrzw8wG43qvAz4EjgAjgaeAyWY2WtJpSBioofPrgNP4bHw2pb23gZnAYeBhSQdjdfUE5gM3Vd6lZMxsZtA6XtKuCqrIVX8NuJz1r8YH6TxJq6JIM3sNeBpYCjwB0C2hgnnABGAWcD6pFTObj7+kb4BR8ZcEIOmCpKXA7Ep6Uk8a+vPDzAbh1v448EYsWfi4m25mvaCDgWpmtwCvAiskfZrSUH98xP8BTJX0a5owSReyd6P+NPTnzoQQbpP0V3FC6MtnwFXAKIgt/WbWHdgInAAWlmioFbgSWCfpUNWyO59WGvozY2at+AEtKwckbUlJvzmE7Qnp3+EWdyiwI75HfQm4DRiTYQa2hHBjiXxdlZYQNvRnoxU/nGVlA356T6JvCM8mpEfx/aDIoprZSNyKLpe0J02BmfXGR3oB+KKU4ljZccBzwO1Af+ABSZszlj1O8svaaWbxuA2SZnZQTzX6FwD34RbhErAXWJDFqpXQ3xHvSprWQT0V6w/l++Pbu3twj85RYK6kT5LKSLqx3HaqpCmEBQgDtWjJbwcWZajk2hD+UoFvrhfu21sHvF9m2Tb+v/zcCkzBZ/DxWNqBhHqq0X83flrdh7/MJcB2Mxsm6UyJsseAi2W09UNCfMX6zawfvv/bDdwL/AgMAk6VU08NiCxm34T0q4vzRRY1mqEAFzuwTABrzGwNsAJ4JcT1NLMrJP2ZVZ2krcBWgIR20sq2xeOCe2oK7tDelbGqyJNRif5Jsfan4y9zNPBRibLNWdspQcX6geeBk5JmFMV9X6pQHfao34ZwaEL6kBC2w7+n/kvAWwnP/pBnd/i/R9Ip3Hr1wC1MImaW5ALLjRrr74O/x1LWtGZUqb8F2Gtmm8zslJkdMLO5wXeeRivuNsr6tJSob2cIJ8Y1mlkffOJfwLdWblHDwemxhI4uxg9YG2JXqG3hWW1mUyQdiZVrwq/FZgFTS4jOgzZqo78N32LsrYfIEu22Ub7+QcAc4HVgEr51WoXvBeP+zH+o9R5V0jEz24af7OcEDRGGbxHflHQeUq5QM7ASGI7fL38VGm0HfgcGAmOAAcCmKtqoJ1XrN7NlwDhgbJnLby2oVH834EtJC8L//WY2FB8siQO1TszGr1BXmlkzfnFxJzAe78uLUcaKl2VJBUmP4hvyj/GZOQe3zMOBPfitybxK26gn1eo3s+XADKBZ0tFOkPwfqtB/Evg6FneY8rwRNUHSMWAEsB4foM8Ag/FJeFd0zw/QlPb1VL0xswJluKe6Cma2EngQ/77gcN56ysHMNgEDJY0tinsZuF/SsPyUpdPpAzX4AKMPJPbjPtXtwBlJJzpVTAWY2WpgGn5YKB6k5ySdy0VUGZjZHfhyuxh4Dz9/rAUWSurspT8zeZzIR+ADNPImLAu/l+SgpRKexE/6O/BlNHrSvjDrMkjah0+yqcAh/HuBRbhvuMuS69LfoEFWupyPs0GDjmgM1AaXBX8Db+ZLKnj2sa4AAAAASUVORK5CYII=",
      "text/latex": [
       "$\\displaystyle 4 C_{1} + C_{2} - C_{6} = 0$"
      ],
      "text/plain": [
       "4⋅C₁ + C₂ - C₆ = 0"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAG4AAAAUCAYAAABs4ghcAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAASdAAAEnQB3mYfeAAAA6lJREFUeJztmV2IVGUYx3+zimhWCn0KWy2adhGEkNSGFdWCCF20NwlFQqUXYmCTEhEh//6CV4aIktFNukR5Y2EUCEJ4VZt3iWF+sPkBfZOlohukTRfve7az48x4zsyZ2R3YHxwezvv5f+Y57+eUKpUKtbC9EHgZWAbcDcwF/gIOA58CuyVdrll5gulW7bZ7gU3AcuAW4GdgH2BJf6bLTq9RuQQYeAOYAXwN7AXOA/cQfowB4DngsXY50Qxdrn0BQe/twGfAMeAh4FVgue2lkv5Iyl8TOOAD4EXgKPC8pMNVHcwC1gP3tsOBFulm7TsJQVsnaUeSaHsr8BqwGViTpPeka9peT3D8e6C/2nEASaOSNgNr26G+Wbpc+3zCbHAaeLcqW8AlYKXt2UliT6ryPEJUrwArJF1s1Jmk0WJkt043a488Fe0BSf+mM6IvXwE3AP1JenqqLAMzgV2SvmuvzsIp02HttsuETU9WvpW0r07efdGeqJN/kjAiFwFfwvjADUb7YQ4xk4XBaDupvUzY8GRliLBDrMWcaM/XyU/S5yYJ0wFs30iIZgU4lFWJ7ceB14EHgXnAs5L2Zqx7mnyOfyTphRrtNKU91n2bsIak+VXSnderK6kvT18tUop27OyWjLjbor2Q83wzm3A22gV8klPMCPB3jvI/1UlvVnvCceCJ1PvVJtpolWREzamTf3NVubHAXYp2lu1pkjKJl7Qf2A9gO5dSSQO5KtSnKe0prkj6JW+nBa9xx6NdVCd/YbRja2APgKTfCFvRGYz/+q7Bdk+j/E5TgPb5tn+0fcr2Htt9GbsuE6bZrM9gg7YORrusWqPtm4ClwCjwTZKe3pxsi89O289IOlbVQIlwFfMSsCKjc51iG81pPwSsJHzJdwBvAsO275d0rlGHRa5xkkZsHyDsHF8BdqSyTViS3peUzC7jArcdeIBwx3ckNnQC+Ae4C3gU6AU+LkpwgTSlPU71CUdsDwM/EA7yW9svexxrCVde220PEC4SHgaeJPjyVrrw2LCUVJG0Cnga+AJYTIj+asKPMkxwaF27PchLUdrjF32U/9eUjiFpBFgC7CYEbAOwgPBRPpK+pwQo1ft3IC+2K+Q4DkxGbM8ETgHvSdo00XoaUeuSOTPxDJW+sO2zvRg4J+lsK213AtvvAJ8DZwkXvBsJ68nQROrKQkuBIwztg6n3LdEOEaamyU4vsAe4FfidsGvrl3RmQlVloLCpcorOMqnOZFNkZypwXcp/TvVqwtp6oi4AAAAASUVORK5CYII=",
      "text/latex": [
       "$\\displaystyle C_{1} - C_{5} = 0$"
      ],
      "text/plain": [
       "C₁ - C₅ = 0"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "C_{4} = 0\n",
      "C_{2} = 0\n",
      "- \\frac{4 C_{5}}{1875} - \\frac{C_{6}}{625} - 4 C_{7} + C_{8} + \\frac{8}{375} = 0\n",
      "\\frac{C_{5}}{625} + \\frac{C_{6}}{1250} + C_{7} - \\frac{8}{375} = 0\n",
      "- \\frac{4 C_{1}}{1875} - \\frac{C_{2}}{625} - 4 C_{3} + C_{4} - C_{8} = 0\n",
      "\\frac{C_{1}}{625} + \\frac{C_{2}}{1250} + C_{3} - C_{7} = 0\n",
      "4 C_{1} + C_{2} - C_{6} = 0\n",
      "C_{1} - C_{5} = 0\n"
     ]
    }
   ],
   "source": [
    "import sympy as sym\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "%config InlineBackend.figure_formats = ['svg']\n",
    "sym.init_printing()\n",
    "\n",
    "x = sym.symbols('x')\n",
    "\n",
    "q = 10 \n",
    "EI = 5000 \n",
    "EA = 20000 \n",
    "F = 150\n",
    "\n",
    "q_AC = 0\n",
    "q_BC = q\n",
    "\n",
    "C_1, C_2, C_3, C_4 = sym.symbols('C_1 C_2 C_3 C_4')\n",
    "C_5, C_6, C_7, C_8 = sym.symbols('C_5 C_6 C_7 C_8')\n",
    "\n",
    "V_AC = -sym.integrate(q_AC, x) + C_1\n",
    "M_AC = sym.integrate(V_AC, x) + C_2\n",
    "kappa_AC = M_AC/EI\n",
    "phi_AC = sym.integrate(kappa_AC, x) + C_3\n",
    "w_AC = -sym.integrate(phi_AC, x) + C_4\n",
    "\n",
    "V_BC = -sym.integrate(q_BC, x) + C_5\n",
    "M_BC = sym.integrate(V_BC, x) + C_6\n",
    "kappa_BC = M_BC/EI\n",
    "phi_BC = sym.integrate(kappa_BC, x) + C_7\n",
    "w_BC = -sym.integrate(phi_BC, x) + C_8\n",
    "\n",
    "eq1 = sym.Eq(w_AC.subs(x,0),0)\n",
    "eq2 = sym.Eq(M_AC.subs(x,0),0)\n",
    "eq3 = sym.Eq(w_BC.subs(x,4),0)\n",
    "eq4 = sym.Eq(phi_BC.subs(x,4),0)\n",
    "eq5 = sym.Eq(w_AC.subs(x,4)-w_BC.subs(x,0),0)\n",
    "eq6 = sym.Eq(phi_AC.subs(x,4)-phi_BC.subs(x,0),0)\n",
    "eq7 = sym.Eq(M_AC.subs(x,4) - M_BC.subs(x,0),0)\n",
    "eq8 = sym.Eq(V_AC.subs(x,4) - V_BC.subs(x,0),0)\n",
    "\n",
    "display(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8)\n",
    "for eq in [eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8]:\n",
    "    print(sym.latex(eq))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "tags": [
     "remove-cell"
    ]
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAzAAAAAVCAYAAACZgQsDAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAASdAAAEnQB3mYfeAAAEgxJREFUeJztnXu0H1V1xz8hLB4FBJXyUCg0FAhUIBREHoYmoDEYrTcgtlJAEEEKrHhFyrO480WBuHjFgCi0SgShtIKEUoI8GhArryUKBcIbAggRSpBHwkMo6R/7/JJf5s5vfvP6zfzuZb5r3XXWnZkzZ//23mefOWfvs8+opUuX0qBBgwYNGjRo0KBBgwbDASvVTUCDBg0aNGjQoEGDBg0apMXKSTcl7QdMA7YC1gR2MLN7KqCrQYMGDRo0aNCgQYMG7yFImg6cBCwC7gCONbNHos91nMBI2hr4CbAE+CnwLPD7hOc3B74MTAL+DFgHeBm4F/gZMNvMXs/zY7JC0neAHYEtgHWBN4CngDnAeWa2KKbOAmCTDq983sw2SGhvfeA54HtmNk3SFOBrwNbAB4GFwN3A2WZ2e86flRr9JItAz0bAKcBklvNjDiAz+0OP2+4rXhSFpM8Dfw2MA7YD1gIuNbP9U9SN6mnmd+XtJ40OrEBPbbwI7fcNP7La6jLtdLhWm60eznIIdRZQgiwAAVOBKcA2wIeBPwL3ARcBF5nZu3l+V1r0mSxKs/HhWh7ZZrZRVdUpC8Nd5nnk2la3pSfdorDeNbPRbfUWUM34fwuwGs6PzwEfAnaKPpTkgZkIjAJOMrNZnR6SNAo3QMcBqwC3AVcAr+A/dBKwJ/BFYHxCe2Xi68BvgBuBF4A1gJ2B6cBhknY2s2di6r0CzIy5vrhLe5/DFeGqoFTH4jPHOcCLwF+EZ/aRdKCZ/STj70mFfpSFpM0CHesBVwMP4Yr4NWCypN2SOlqBdvuOFyXhn3ADtxj4HTA2Q91lelrwXZn6SaMDK9BUCy9C233HD/LZ6sJ2GpZ9AFRuq0eQHKAcWewLfB//qLkZeBpYH9gb+BdgL0n7mlnpG3b7VBZl2njIKNs8NqqqOmVgBMk8b5+F5XpyK97nohgP7AFcF3Ov5+O/md2CT2KQdDewo6Q1zWyFNpImMB8I5YMJzwD8CDgImA/sZ2b3RohfHTgaHxiqwvvM7M3oRUmnAicCJwBHxNR72cym52hvKj4IPgLcBDwPbGtmL7S1PRGYh89CezKBoT9lcT6uuNPM7Nw2Ws7GO+CpwOE9aLcfeVEGvo4buMfwFZs449MJLT29teC7svaTRgeWoy5eQH/yI4+tLmqnb5W0AXAM9djqkSIHKEEW+ELp3wDXtntaJJ0I3AXsg09mrszRTjf0oyzKtPGQXbZ5bFRVdcrASJF53j4Ly/VkDzP7v5h3tLzPF8bUrXr8fwj4K+D9RCZJSe6jltvo7U4PSDoaV4QHgZ2jigBgZm+Y2al0ZmTpiBNqwL+HcvOy2pK0Nj5TvQbYCOfpne0DYqDpZuA14E/LajtCR9/JQtIYfDVjAR4qsAI5eHjiAZLWKLndvuNFWTCzm83s0ayrke162jJYed+Vsd1GB5bTVAsvQtt9x4/QZiW2Okb/N6EGW/1elwMMlYWZzTOza6JhYmb2e+AH4d8JZbXfRke/yqI0Gx/el1q2eWxUVXXKwEiSed4+20lP2u5/BPfkPAtcm5aeDm2VIefWHGR09EbiJv4uhG2Iz5zeAb5gZq8lPW9mb+Rtq63N2cCXgIPNbHaOV3w2lP/T4f6qkvbH4yGXhOdujRNyG6bgbsifAY/isbs7SVrXzF5so313PK5xTg66E1G1LDLIYY9Q3hAzOL0m6Ve4cu8M/FcRmtpoG456WQXa9bQosvSTRgeWo3JeBNr61T4kIclWF7XTUIOt7mO9TEKvx8xuaH28vJPi2dQYprLohqw2Pk62eWxUVXUKYYTKPA7d+mw3PflqKH/YoR/3zfifNIFZJ5SdhDiIb7K5yMzuz9pwFZB0DJ49bW18s9PHcWbP6FBlA+CSyLUnJR1sZr/oUGcqLsQbzexNSccBZwPzJc3B3XSb4S7yG1muHGVikP6UxZahHJI9IuBRXHm3oLwPtkH6kxd1Y5melvCuLP2k0YHlqIMX0L/8WIaMtrqQnQYws5dqsNWDjCw5QAmySKBlZeDA8O/Pu1OfCYP0uSxyIJGvKWWbx0ZVVacoBhl5Ms/TZzvqSQif2x94F99/Foeqx/+Wp2nt6I3YELKwyWk8sBTPahCHgVBGf0gvcQKe0vmqbg8GHIO7qQZxof4cmGRm/xvz7EX4xq0N8M1Q2wAXAJsC10naLlpB0mp4VoXrWu48M5uJx+uuDBwKHI9vUnwGz2rxQvQ9JWAglFXJIq0cWgr3Sof7revrlEBTCwOh7Ge9rBRxeloAWftJowPLUQcvoH/tQzvS2upS7DTUYqsHQtlvetmOysfMBMwAPgLMNbPrU9KfFgOh7GdZpEZKvqaRbR4bVVWdohgI5YiQeRtS99kUevIFnOfXWXwCgDrG/wWhnBC9sYIHRtJk4JPA7vimmQtCHCqR59bEZ0xLgTsTGh6C4J7/R2AHYENgXzO7Ik1dM1uIZypJBQsp3eQp43bFDeJvJX3GzH4TeVaR6vcDh0taDHwDz+wwNfLMJHzmu0w5JR0LnAbMAs7DU0+PBU4HLpU0zsyOTfsbuqGgLE7AB/AtgbfwfNsndFudyCqHBIwKZSl7MAry4kh8xXXTcOkB4Ntm1jUGNC0/lJyCMA6pUmemwBA9zYuc/SQJfaMDkfeciIcbfM/Mjur2fEl9olReQOE+MR0fGNuRmB4X8vEira0uy06Htiqz1UX1MoS/zAA+jf+Wx4CjErwcwPAYM+MgaVp4x0PAAWnpT4OCfWIB8Tb8fDM7MqluieNmHLryNYtsE5DHRlVVpyMKynw0rsv749+sC4FLgelmlhja2GOZt9rIItduenJYKC/o0FYd4/8l+Cb/MyXthntt/tnMFkQ9MJPx7As7BsLOJR6tzY2vWvY82Wvgeba7fhSUBTN73syuwoX3QeDiDNVbmwh3j7k3FY+jvhZA0gTgO8B/mNnRZvaEmb0elGgqvinqG2FjU1koIosJeIaIXfFYxXeAmyR9IKlSBrRm10NcfwHvizxXFEV48Tt8BXYHXP/nAXMkbVsSbQCPAw9n+HuupHZX0NMeoVM/GU46AICknfEV+U4xxHlRNS+gOD8exgft1t82ZREWhwK2OrWdhlpsdW45SFoH+BU+2E/Bz6wZxFOn9gRVjZlxCItJ38WzRE00s5cytJ0GRfrER1mxP3wyXP9pSbTlRWob30W2eWxUVXWKoIjMjwOOxA91H4un/z0C9670DVL22Y56Ij/7cVf8W2huxuZ7Nv6b2bP4otLbuIf8RMJC88qRBwclfRM/A+Zy4HpJm8RszlkSytUlje6yYS9KzHWE3NJSdDLXW5jZU5LmA+OiGzcT0BokVsiSEGblnwXmmVmL+Z8J5ZAUeGb2uqS7cAXaHngiz2+IQRFZfKr9f0kH4Iq0G55VrSgeDuUWHe63smR0io/MiiK8uDpy6SRJ/wDsQkkfsma2ZxnvyYIOetoLxPYThpEOwLIMLZcChwDfLImmFqrmBRTkB/BOnBe+18hhq7PYaajeVheRw7HAQjM7sO3akyXQ1BUVjJlEnhkEzsEXUPfsUch1kXFihbAcScfjC1OJnrBeIq+N7yDbPDaqqjpFUKT/7YovdLS+iRZIugb4WEm0lYpOfTaFnnTbvJ+Eno3/8sOGf4Av4nwVeMTM3oaYPTBm9mr4mLsSPxF3q5hnXsDj0lahS3pDSd1O+qwaHwplWgHtEsroILY7Psttd8WtGspO6Tdb1/+Ysu2uKFkWa+E6UdaKV+vjYFK0bUlr4ROlN/DQtcIoixeSRkv6Iu5qva0M2mpEnJ72Ap36yXDTgQuBK8xsXhn0RFApL6AUfoyR9KykJyX9q6RNy6ItBbLY6ix2Giq21QXlMADcIekySS9IukfSUWGvahXo5Zi5DPKkCucA9+Cel554mEocJ1bBw4p+ZD1MRZ8CRWx8VLZ5bFRVdXKjoMxvA/aQNDbc3xqPWMnqpagScX22o56EvTEH4Jv3f5ijvV6O/61MZtPN7IHW5AWSz4Fpbd7vFE40M5TntwQbIW6UpL1wT04pkLShpLFhlbTTM2Plh5RFr68kP+BnPeA2M/tD272/jAubkrQJHhsNQw802xsXdvvK/S9DeZikD0fetRcurDeJfBRLmi1pqaSDOv2uLpgZyqKymIkPHolGI40cAMzsceAG3N0XjQ8WPlu/2MyWrHCjGD9mhjIzLyRtE2I538JD6wbM7L5uDablR02I09NcyNNPhpMOSDoUP7zs5KwNptGBmngB+fvEnfigNhkPqVsPuD1OByLvS9UfstrqEu001GOrZ4YyqxzG4PryFPApPLxqBinOqRgGY2ar3snhN92Ne166enrqGifaMIBvPp6dpsEejhNJfM0k2zw2qqo6ge46ZD4D98rPl/Q2vj/2x2Z2frcGeyXzPH2W5G+BffGDIuda/Ob9Ssf/CN4fygXRG0lplFuznE6TnFnAtsCXgfsk3YC7gd4GNsazIWwEXJbQRlacTsipTWejMRk4Q9KtuGt3EbA+frrpGHyj5qGROvsCx0u6GXfNv4an05yCp92bC5wZqTOAK8jzbdeuAG4CPgE8KOmq0N5WeMjCKOB4M1sUeVeLx3lz3ReWhaQz8Bn6+BTuwzRyaOEI/CNglqQ98QOkPoaHKT4CnBRTpwg/ivDiYWAcPijtA1wsaYJ1T7mYhR+FIGmA5dlUWgZsF3m+eYAXzeyYtioDDNXTvO/K009gGOiApC3xDd3jzSzPqntaHaiaF5CzT4Rw3xbuk5/O/AR+ENzZCe2l5UVWW12WnYZ6bHVe27QScLeZtWLufytpC/yDIHo4XBT9PmYi6UvAKfhq8S+BaRoaXr7Ahp6jUdc40cIheLamtHsVU/WLMm08+WSbx0ZVVacOmX8e+HtgP3zyMg74rqQnzaybt6JXMs8j1wE660lr8/6FCb+lyvG/HR1lnjSBSXSJBpfpIZKuxBm1E76hbTGedeF2fKPQf3YhrmzchAthN2A7/GN0Cc6oS4BZMZsCb8YzcW2Pu8LWAF4G/jvUuaTdRSzpo7iin9P+EjN7V9Kn8YHl7/AY6j/BQ7LmhrZviKF5G1wZcm2yLioLSWfhrvCJZvZYHhoSaHtc0o74ADUZz6KzEDcm6rBBMzc/ivAifLi2fv+vg5wHga9kpaOHGIcbxHaMCX/gq7THQGc9zfOugEz9pIVhogO7AOsC97d9PI0Gdpd0OLCGmb2VlZYY2irlRWizFFttZkvksdWbJz2XAVltdSl2OvyWym11ATksxD+c2jEf31RcBmobMwP+PJSjcXsbh18w9AOwlnEClq06fwJf1S4b4yjPxmeWbR4bVVUd6pH5WcCZZtbyzNwX5H8i+cKt4jCObGNxJrkm6YmkrfDJW7fN+1WO/3EY8u5RS5fGz1Pkm9NOB/Y2z2xQOiQtJUMa5X6BpNPwDBRjzKzQZkp5hplFwFlWYnrlDO3PAv4Wn7zMr7r9GHrWoUZ+RGiZBzxn5aQyrhxl6mmVqEMHQpsbRS5fhKdsPA14oI449z7rD6vhK2/fN7NT6qQlDcrW/7pkIekyYGMzG9927VvAPma2dVV0FMFIkUVb+9PxDcUbW5dUuj2mY1ja+Dyosf8tAszMzmu7dgLwFTPbrCo6imA464mkq/HDhTeIeo+SPDCPhvIQSb8GnjWzd0sgZk08zryFTSWNA14ys6eLvr8iTAXuLUkRxuMuzKSQjJ5A0vm452UAeKktpnKxmS2ump6AWvghaQa++vIMnsxgP3yz35Qq6SgZZepplahcB8zsZXw1aRkkLcHtUp2nNtdpH87EsxE+jcdUn4yvuv24alpyomz9r0sW5wC3SToJ+Dd8BXQavgI8XDBSZNHa5H0wvg+itslLwHC18XlQl8yvwUOnnsQ9odvjx41kSS1eN4adnkhaFfcuTQBeBYYczJnkgVkduAs/CbeF7c3snoJETSAmdSVuDA4q8u4G2RA8YLG3zGx6lbTUjRBrOhGPP30FT518hpV/AnSDYQJJtwD3W4qDLEciJF2O74tbFx887gBO7gdP7XsN8lSip+EhHE/jm2bPrcMr+F6HpEnA9cCWZlZmyvMGfQh5tqxv4ZOA9fDwp8uBUyz+NPsGBaGhhyhPt6GHaHaewISXrIqnMNsKTyl7odVwJkCDBg0aNGjQoEGDBg1GNoKj4+N4yODtnRwn/w8FyZdc/s0UKwAAAABJRU5ErkJggg==",
      "text/latex": [
       "$\\displaystyle \\left\\{ C_{1} : \\frac{35}{8}, \\  C_{2} : 0, \\  C_{3} : - \\frac{1}{150}, \\  C_{4} : 0, \\  C_{5} : \\frac{35}{8}, \\  C_{6} : \\frac{35}{2}, \\  C_{7} : \\frac{1}{3000}, \\  C_{8} : \\frac{13}{750}\\right\\}$"
      ],
      "text/plain": [
       "⎧                                                                        13 ⎫\n",
       "⎨C₁: 35/8, C₂: 0, C₃: -1/150, C₄: 0, C₅: 35/8, C₆: 35/2, C₇: 1/3000, C₈: ───⎬\n",
       "⎩                                                                        750⎭"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "'\\\\left\\\\{ C_{1} : \\\\frac{35}{8}, \\\\  C_{2} : 0, \\\\  C_{3} : - \\\\frac{1}{150}, \\\\  C_{4} : 0, \\\\  C_{5} : \\\\frac{35}{8}, \\\\  C_{6} : \\\\frac{35}{2}, \\\\  C_{7} : \\\\frac{1}{3000}, \\\\  C_{8} : \\\\frac{13}{750}\\\\right\\\\}'"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sol = sym.solve([eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8],(C_1,C_2,C_3,C_4,C_5,C_6,C_7,C_8))\n",
    "display(sol)\n",
    "sym.latex(sol)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "tags": [
     "remove-cell"
    ]
   },
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?>\n",
       "<!DOCTYPE svg PUBLIC \"-//W3C//DTD SVG 1.1//EN\"\n",
       "  \"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd\">\n",
       "<svg xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"397.070085pt\" height=\"280.512pt\" viewBox=\"0 0 397.070085 280.512\" xmlns=\"http://www.w3.org/2000/svg\" version=\"1.1\">\n",
       " <metadata>\n",
       "  <rdf:RDF xmlns:dc=\"http://purl.org/dc/elements/1.1/\" xmlns:cc=\"http://creativecommons.org/ns#\" 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>2025-07-02T16:58:07.214407</dc:date>\n",
       "    <dc:format>image/svg+xml</dc:format>\n",
       "    <dc:creator>\n",
       "     <cc:Agent>\n",
       "      <dc:title>Matplotlib v3.10.3, 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-linejoin: round; stroke-linecap: butt}</style>\n",
       " </defs>\n",
       " <g id=\"figure_1\">\n",
       "  <g id=\"patch_1\">\n",
       "   <path d=\"M 0 280.512 \n",
       "L 397.070085 280.512 \n",
       "L 397.070085 0 \n",
       "L 0 0 \n",
       "z\n",
       "\" style=\"fill: #ffffff\"/>\n",
       "  </g>\n",
       "  <g id=\"axes_1\">\n",
       "   <g id=\"patch_2\">\n",
       "    <path d=\"M 32.750085 273.312 \n",
       "L 389.870085 273.312 \n",
       "L 389.870085 7.2 \n",
       "L 32.750085 7.2 \n",
       "z\n",
       "\" style=\"fill: #ffffff\"/>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_1\">\n",
       "    <g id=\"xtick_1\">\n",
       "     <g id=\"line2d_1\">\n",
       "      <defs>\n",
       "       <path id=\"m49336a3f35\" d=\"M 0 0 \n",
       "L 0 3.5 \n",
       "\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use xlink:href=\"#m49336a3f35\" x=\"48.982812\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_1\">\n",
       "      <!-- 0.0 -->\n",
       "      <g transform=\"translate(41.03125 205.450546) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-30\" d=\"M 2034 4250 \n",
       "Q 1547 4250 1301 3770 \n",
       "Q 1056 3291 1056 2328 \n",
       "Q 1056 1369 1301 889 \n",
       "Q 1547 409 2034 409 \n",
       "Q 2525 409 2770 889 \n",
       "Q 3016 1369 3016 2328 \n",
       "Q 3016 3291 2770 3770 \n",
       "Q 2525 4250 2034 4250 \n",
       "z\n",
       "M 2034 4750 \n",
       "Q 2819 4750 3233 4129 \n",
       "Q 3647 3509 3647 2328 \n",
       "Q 3647 1150 3233 529 \n",
       "Q 2819 -91 2034 -91 \n",
       "Q 1250 -91 836 529 \n",
       "Q 422 1150 422 2328 \n",
       "Q 422 3509 836 4129 \n",
       "Q 1250 4750 2034 4750 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "        <path id=\"DejaVuSans-2e\" d=\"M 684 794 \n",
       "L 1344 794 \n",
       "L 1344 0 \n",
       "L 684 0 \n",
       "L 684 794 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-30\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-2e\" transform=\"translate(63.623047 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(95.410156 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_2\">\n",
       "     <g id=\"line2d_2\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m49336a3f35\" x=\"89.564631\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_2\">\n",
       "      <!-- 0.5 -->\n",
       "      <g transform=\"translate(81.613068 205.450546) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-35\" d=\"M 691 4666 \n",
       "L 3169 4666 \n",
       "L 3169 4134 \n",
       "L 1269 4134 \n",
       "L 1269 2991 \n",
       "Q 1406 3038 1543 3061 \n",
       "Q 1681 3084 1819 3084 \n",
       "Q 2600 3084 3056 2656 \n",
       "Q 3513 2228 3513 1497 \n",
       "Q 3513 744 3044 326 \n",
       "Q 2575 -91 1722 -91 \n",
       "Q 1428 -91 1123 -41 \n",
       "Q 819 9 494 109 \n",
       "L 494 744 \n",
       "Q 775 591 1075 516 \n",
       "Q 1375 441 1709 441 \n",
       "Q 2250 441 2565 725 \n",
       "Q 2881 1009 2881 1497 \n",
       "Q 2881 1984 2565 2268 \n",
       "Q 2250 2553 1709 2553 \n",
       "Q 1456 2553 1204 2497 \n",
       "Q 953 2441 691 2322 \n",
       "L 691 4666 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-30\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-2e\" transform=\"translate(63.623047 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-35\" transform=\"translate(95.410156 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_3\">\n",
       "     <g id=\"line2d_3\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m49336a3f35\" x=\"130.146449\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_3\">\n",
       "      <!-- 1.0 -->\n",
       "      <g transform=\"translate(122.194886 205.450546) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-31\" d=\"M 794 531 \n",
       "L 1825 531 \n",
       "L 1825 4091 \n",
       "L 703 3866 \n",
       "L 703 4441 \n",
       "L 1819 4666 \n",
       "L 2450 4666 \n",
       "L 2450 531 \n",
       "L 3481 531 \n",
       "L 3481 0 \n",
       "L 794 0 \n",
       "L 794 531 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-31\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-2e\" transform=\"translate(63.623047 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(95.410156 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_4\">\n",
       "     <g id=\"line2d_4\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m49336a3f35\" x=\"170.728267\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_4\">\n",
       "      <!-- 1.5 -->\n",
       "      <g transform=\"translate(162.776705 205.450546) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-31\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-2e\" transform=\"translate(63.623047 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-35\" transform=\"translate(95.410156 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_5\">\n",
       "     <g id=\"line2d_5\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m49336a3f35\" x=\"211.310085\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_5\">\n",
       "      <!-- 2.0 -->\n",
       "      <g transform=\"translate(203.358523 205.450546) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-32\" d=\"M 1228 531 \n",
       "L 3431 531 \n",
       "L 3431 0 \n",
       "L 469 0 \n",
       "L 469 531 \n",
       "Q 828 903 1448 1529 \n",
       "Q 2069 2156 2228 2338 \n",
       "Q 2531 2678 2651 2914 \n",
       "Q 2772 3150 2772 3378 \n",
       "Q 2772 3750 2511 3984 \n",
       "Q 2250 4219 1831 4219 \n",
       "Q 1534 4219 1204 4116 \n",
       "Q 875 4013 500 3803 \n",
       "L 500 4441 \n",
       "Q 881 4594 1212 4672 \n",
       "Q 1544 4750 1819 4750 \n",
       "Q 2544 4750 2975 4387 \n",
       "Q 3406 4025 3406 3419 \n",
       "Q 3406 3131 3298 2873 \n",
       "Q 3191 2616 2906 2266 \n",
       "Q 2828 2175 2409 1742 \n",
       "Q 1991 1309 1228 531 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-32\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-2e\" transform=\"translate(63.623047 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(95.410156 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_6\">\n",
       "     <g id=\"line2d_6\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m49336a3f35\" x=\"251.891903\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_6\">\n",
       "      <!-- 2.5 -->\n",
       "      <g transform=\"translate(243.940341 205.450546) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-32\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-2e\" transform=\"translate(63.623047 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-35\" transform=\"translate(95.410156 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_7\">\n",
       "     <g id=\"line2d_7\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m49336a3f35\" x=\"292.473722\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_7\">\n",
       "      <!-- 3.0 -->\n",
       "      <g transform=\"translate(284.522159 205.450546) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-33\" d=\"M 2597 2516 \n",
       "Q 3050 2419 3304 2112 \n",
       "Q 3559 1806 3559 1356 \n",
       "Q 3559 666 3084 287 \n",
       "Q 2609 -91 1734 -91 \n",
       "Q 1441 -91 1130 -33 \n",
       "Q 819 25 488 141 \n",
       "L 488 750 \n",
       "Q 750 597 1062 519 \n",
       "Q 1375 441 1716 441 \n",
       "Q 2309 441 2620 675 \n",
       "Q 2931 909 2931 1356 \n",
       "Q 2931 1769 2642 2001 \n",
       "Q 2353 2234 1838 2234 \n",
       "L 1294 2234 \n",
       "L 1294 2753 \n",
       "L 1863 2753 \n",
       "Q 2328 2753 2575 2939 \n",
       "Q 2822 3125 2822 3475 \n",
       "Q 2822 3834 2567 4026 \n",
       "Q 2313 4219 1838 4219 \n",
       "Q 1578 4219 1281 4162 \n",
       "Q 984 4106 628 3988 \n",
       "L 628 4550 \n",
       "Q 988 4650 1302 4700 \n",
       "Q 1616 4750 1894 4750 \n",
       "Q 2613 4750 3031 4423 \n",
       "Q 3450 4097 3450 3541 \n",
       "Q 3450 3153 3228 2886 \n",
       "Q 3006 2619 2597 2516 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-33\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-2e\" transform=\"translate(63.623047 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(95.410156 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_8\">\n",
       "     <g id=\"line2d_8\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m49336a3f35\" x=\"333.05554\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_8\">\n",
       "      <!-- 3.5 -->\n",
       "      <g transform=\"translate(325.103977 205.450546) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-33\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-2e\" transform=\"translate(63.623047 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-35\" transform=\"translate(95.410156 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_9\">\n",
       "     <g id=\"line2d_9\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m49336a3f35\" x=\"373.637358\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_9\">\n",
       "      <!-- 4.0 -->\n",
       "      <g transform=\"translate(365.685795 205.450546) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-34\" d=\"M 2419 4116 \n",
       "L 825 1625 \n",
       "L 2419 1625 \n",
       "L 2419 4116 \n",
       "z\n",
       "M 2253 4666 \n",
       "L 3047 4666 \n",
       "L 3047 1625 \n",
       "L 3713 1625 \n",
       "L 3713 1100 \n",
       "L 3047 1100 \n",
       "L 3047 0 \n",
       "L 2419 0 \n",
       "L 2419 1100 \n",
       "L 313 1100 \n",
       "L 313 1709 \n",
       "L 2253 4666 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-34\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-2e\" transform=\"translate(63.623047 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(95.410156 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"text_10\">\n",
       "     <!-- $x$ -->\n",
       "     <g transform=\"translate(208.310085 219.128671) scale(0.1 -0.1)\">\n",
       "      <defs>\n",
       "       <path id=\"DejaVuSans-Oblique-78\" d=\"M 3841 3500 \n",
       "L 2234 1784 \n",
       "L 3219 0 \n",
       "L 2559 0 \n",
       "L 1819 1388 \n",
       "L 531 0 \n",
       "L -166 0 \n",
       "L 1556 1844 \n",
       "L 641 3500 \n",
       "L 1300 3500 \n",
       "L 1972 2234 \n",
       "L 3144 3500 \n",
       "L 3841 3500 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "      </defs>\n",
       "      <use xlink:href=\"#DejaVuSans-Oblique-78\" transform=\"translate(0 0.3125)\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_2\">\n",
       "    <g id=\"ytick_1\">\n",
       "     <g id=\"line2d_10\">\n",
       "      <defs>\n",
       "       <path id=\"m81db698713\" d=\"M 0 0 \n",
       "L -3.5 0 \n",
       "\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use xlink:href=\"#m81db698713\" x=\"48.982812\" y=\"38.35779\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_11\">\n",
       "      <!-- −40 -->\n",
       "      <g transform=\"translate(20.878125 42.157009) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-2212\" d=\"M 678 2272 \n",
       "L 4684 2272 \n",
       "L 4684 1741 \n",
       "L 678 1741 \n",
       "L 678 2272 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-2212\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-34\" transform=\"translate(83.789062 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(147.412109 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_2\">\n",
       "     <g id=\"line2d_11\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m81db698713\" x=\"48.982812\" y=\"76.48137\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_12\">\n",
       "      <!-- −30 -->\n",
       "      <g transform=\"translate(20.878125 80.280588) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-2212\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-33\" transform=\"translate(83.789062 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(147.412109 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_3\">\n",
       "     <g id=\"line2d_12\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m81db698713\" x=\"48.982812\" y=\"114.604949\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_13\">\n",
       "      <!-- −20 -->\n",
       "      <g transform=\"translate(20.878125 118.404168) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-2212\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-32\" transform=\"translate(83.789062 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(147.412109 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_4\">\n",
       "     <g id=\"line2d_13\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m81db698713\" x=\"48.982812\" y=\"152.728529\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_14\">\n",
       "      <!-- −10 -->\n",
       "      <g transform=\"translate(20.878125 156.527748) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-2212\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-31\" transform=\"translate(83.789062 0)\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(147.412109 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_5\">\n",
       "     <g id=\"line2d_14\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m81db698713\" x=\"48.982812\" y=\"190.852109\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_15\">\n",
       "      <!-- 0 -->\n",
       "      <g transform=\"translate(35.620312 194.651328) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_6\">\n",
       "     <g id=\"line2d_15\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m81db698713\" x=\"48.982812\" y=\"228.975689\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_16\">\n",
       "      <!-- 10 -->\n",
       "      <g transform=\"translate(29.257812 232.774907) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-31\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(63.623047 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_7\">\n",
       "     <g id=\"line2d_16\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m81db698713\" x=\"48.982812\" y=\"267.099268\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_17\">\n",
       "      <!-- 20 -->\n",
       "      <g transform=\"translate(29.257812 270.898487) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-32\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" transform=\"translate(63.623047 0)\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"text_18\">\n",
       "     <!-- $M_{BC}$ -->\n",
       "     <g transform=\"translate(14.798437 149.556) rotate(-90) scale(0.1 -0.1)\">\n",
       "      <defs>\n",
       "       <path id=\"DejaVuSans-Oblique-4d\" d=\"M 1081 4666 \n",
       "L 2028 4666 \n",
       "L 2572 1522 \n",
       "L 4378 4666 \n",
       "L 5350 4666 \n",
       "L 4441 0 \n",
       "L 3828 0 \n",
       "L 4622 4091 \n",
       "L 2791 897 \n",
       "L 2175 897 \n",
       "L 1581 4103 \n",
       "L 788 0 \n",
       "L 172 0 \n",
       "L 1081 4666 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-Oblique-42\" d=\"M 1081 4666 \n",
       "L 2694 4666 \n",
       "Q 3350 4666 3675 4422 \n",
       "Q 4000 4178 4000 3688 \n",
       "Q 4000 3238 3720 2911 \n",
       "Q 3441 2584 2988 2516 \n",
       "Q 3375 2428 3569 2181 \n",
       "Q 3763 1934 3763 1522 \n",
       "Q 3763 819 3242 409 \n",
       "Q 2722 0 1819 0 \n",
       "L 172 0 \n",
       "L 1081 4666 \n",
       "z\n",
       "M 1234 2228 \n",
       "L 903 519 \n",
       "L 1919 519 \n",
       "Q 2491 519 2800 781 \n",
       "Q 3109 1044 3109 1522 \n",
       "Q 3109 1891 2904 2059 \n",
       "Q 2700 2228 2247 2228 \n",
       "L 1234 2228 \n",
       "z\n",
       "M 1606 4147 \n",
       "L 1331 2741 \n",
       "L 2272 2741 \n",
       "Q 2775 2741 3058 2959 \n",
       "Q 3341 3178 3341 3566 \n",
       "Q 3341 3869 3150 4008 \n",
       "Q 2959 4147 2541 4147 \n",
       "L 1606 4147 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-Oblique-43\" d=\"M 4447 4306 \n",
       "L 4319 3641 \n",
       "Q 4019 3944 3683 4091 \n",
       "Q 3347 4238 2956 4238 \n",
       "Q 2422 4238 2017 3981 \n",
       "Q 1613 3725 1319 3200 \n",
       "Q 1131 2863 1032 2486 \n",
       "Q 934 2109 934 1728 \n",
       "Q 934 1091 1264 756 \n",
       "Q 1594 422 2222 422 \n",
       "Q 2656 422 3056 561 \n",
       "Q 3456 700 3834 978 \n",
       "L 3688 231 \n",
       "Q 3316 72 2936 -9 \n",
       "Q 2556 -91 2175 -91 \n",
       "Q 1278 -91 773 396 \n",
       "Q 269 884 269 1753 \n",
       "Q 269 2309 461 2846 \n",
       "Q 653 3384 1013 3828 \n",
       "Q 1394 4300 1883 4525 \n",
       "Q 2372 4750 3022 4750 \n",
       "Q 3422 4750 3780 4639 \n",
       "Q 4138 4528 4447 4306 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "      </defs>\n",
       "      <use xlink:href=\"#DejaVuSans-Oblique-4d\" transform=\"translate(0 0.09375)\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-Oblique-42\" transform=\"translate(86.279297 -16.3125) scale(0.7)\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-Oblique-43\" transform=\"translate(134.301758 -16.3125) scale(0.7)\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"line2d_17\">\n",
       "    <path d=\"M 48.982812 257.568374 \n",
       "L 52.262151 258.211157 \n",
       "L 55.54149 258.791704 \n",
       "L 58.820829 259.310015 \n",
       "L 62.100168 259.76609 \n",
       "L 65.379507 260.159929 \n",
       "L 68.658846 260.491531 \n",
       "L 71.938184 260.760898 \n",
       "L 75.217523 260.968028 \n",
       "L 78.496862 261.112921 \n",
       "L 81.776201 261.195579 \n",
       "L 85.05554 261.216 \n",
       "L 88.334879 261.174185 \n",
       "L 91.614217 261.070134 \n",
       "L 94.893556 260.903846 \n",
       "L 98.172895 260.675323 \n",
       "L 101.452234 260.384563 \n",
       "L 104.731573 260.031567 \n",
       "L 108.010912 259.616334 \n",
       "L 111.290251 259.138866 \n",
       "L 114.569589 258.599161 \n",
       "L 117.848928 257.99722 \n",
       "L 121.128267 257.333043 \n",
       "L 124.407606 256.606629 \n",
       "L 127.686945 255.81798 \n",
       "L 130.966284 254.967094 \n",
       "L 134.245622 254.053972 \n",
       "L 137.524961 253.078613 \n",
       "L 140.8043 252.041018 \n",
       "L 144.083639 250.941188 \n",
       "L 147.362978 249.779121 \n",
       "L 150.642317 248.554817 \n",
       "L 153.921655 247.268278 \n",
       "L 157.200994 245.919502 \n",
       "L 160.480333 244.50849 \n",
       "L 163.759672 243.035242 \n",
       "L 167.039011 241.499757 \n",
       "L 170.31835 239.902037 \n",
       "L 173.597689 238.24208 \n",
       "L 176.877027 236.519887 \n",
       "L 180.156366 234.735457 \n",
       "L 183.435705 232.888792 \n",
       "L 186.715044 230.97989 \n",
       "L 189.994383 229.008752 \n",
       "L 193.273722 226.975377 \n",
       "L 196.55306 224.879767 \n",
       "L 199.832399 222.72192 \n",
       "L 203.111738 220.501837 \n",
       "L 206.391077 218.219518 \n",
       "L 209.670416 215.874963 \n",
       "L 212.949755 213.468171 \n",
       "L 216.229093 210.999143 \n",
       "L 219.508432 208.467879 \n",
       "L 222.787771 205.874379 \n",
       "L 226.06711 203.218642 \n",
       "L 229.346449 200.500669 \n",
       "L 232.625788 197.72046 \n",
       "L 235.905127 194.878015 \n",
       "L 239.184465 191.973334 \n",
       "L 242.463804 189.006416 \n",
       "L 245.743143 185.977262 \n",
       "L 249.022482 182.885872 \n",
       "L 252.301821 179.732245 \n",
       "L 255.58116 176.516383 \n",
       "L 258.860498 173.238284 \n",
       "L 262.139837 169.897949 \n",
       "L 265.419176 166.495377 \n",
       "L 268.698515 163.03057 \n",
       "L 271.977854 159.503526 \n",
       "L 275.257193 155.914246 \n",
       "L 278.536532 152.26273 \n",
       "L 281.81587 148.548977 \n",
       "L 285.095209 144.772989 \n",
       "L 288.374548 140.934764 \n",
       "L 291.653887 137.034303 \n",
       "L 294.933226 133.071605 \n",
       "L 298.212565 129.046672 \n",
       "L 301.491903 124.959502 \n",
       "L 304.771242 120.810096 \n",
       "L 308.050581 116.598454 \n",
       "L 311.32992 112.324575 \n",
       "L 314.609259 107.98846 \n",
       "L 317.888598 103.590109 \n",
       "L 321.167936 99.129522 \n",
       "L 324.447275 94.606699 \n",
       "L 327.726614 90.021639 \n",
       "L 331.005953 85.374343 \n",
       "L 334.285292 80.664811 \n",
       "L 337.564631 75.893043 \n",
       "L 340.84397 71.059038 \n",
       "L 344.123308 66.162797 \n",
       "L 347.402647 61.20432 \n",
       "L 350.681986 56.183607 \n",
       "L 353.961325 51.100658 \n",
       "L 357.240664 45.955472 \n",
       "L 360.520003 40.74805 \n",
       "L 363.799341 35.478392 \n",
       "L 367.07868 30.146497 \n",
       "L 370.358019 24.752367 \n",
       "L 373.637358 19.296 \n",
       "\" clip-path=\"url(#pcc06866bbc)\" style=\"fill: none; stroke: #1f77b4; stroke-width: 1.5; stroke-linecap: square\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_3\">\n",
       "    <path d=\"M 48.982812 273.312 \n",
       "L 48.982812 7.2 \n",
       "\" style=\"fill: none; stroke: #000000; stroke-width: 0.8; stroke-linejoin: miter; stroke-linecap: square\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_4\">\n",
       "    <path d=\"M 389.870085 273.312 \n",
       "L 389.870085 7.2 \n",
       "\" style=\"fill: none\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_5\">\n",
       "    <path d=\"M 32.750085 190.852109 \n",
       "L 389.870085 190.852109 \n",
       "\" style=\"fill: none; stroke: #000000; stroke-width: 0.8; stroke-linejoin: miter; stroke-linecap: square\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_6\">\n",
       "    <path d=\"M 32.750085 7.2 \n",
       "L 389.870085 7.2 \n",
       "\" style=\"fill: none\"/>\n",
       "   </g>\n",
       "  </g>\n",
       " </g>\n",
       " <defs>\n",
       "  <clipPath id=\"pcc06866bbc\">\n",
       "   <rect x=\"32.750085\" y=\"7.2\" width=\"357.12\" height=\"266.112\"/>\n",
       "  </clipPath>\n",
       " </defs>\n",
       "</svg>\n"
      ],
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "M_numpy = sym.lambdify(x, M_BC.subs(sol))\n",
    "x_numpy = np.linspace(0,4,100)\n",
    "plt.figure()\n",
    "plt.plot(x_numpy,M_numpy(x_numpy))\n",
    "plt.xlabel('$x$')\n",
    "plt.ylabel('$M_{BC}$');\n",
    "ax = plt.gca()\n",
    "ax.spines['right'].set_color('none')\n",
    "ax.spines['top'].set_color('none')\n",
    "ax.spines['bottom'].set_position('zero')\n",
    "ax.spines['left'].set_position('zero')\n",
    "ax.invert_yaxis()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "tags": [
     "remove-cell"
    ]
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\\left[\\begin{matrix}0 & 0 & 0 & 1 & 0 & 0 & 0 & 0\\\\0 & 1 & 0 & 0 & 0 & 0 & 0 & 0\\\\0 & 0 & 0 & 0 & - \\frac{4}{1875} & - \\frac{1}{625} & -4 & 1\\\\0 & 0 & 0 & 0 & \\frac{1}{625} & \\frac{1}{1250} & 1 & 0\\\\- \\frac{4}{1875} & - \\frac{1}{625} & -4 & 1 & 0 & 0 & 0 & -1\\\\\\frac{1}{625} & \\frac{1}{1250} & 1 & 0 & 0 & 0 & -1 & 0\\\\4 & 1 & 0 & 0 & 0 & -1 & 0 & 0\\\\1 & 0 & 0 & 0 & -1 & 0 & 0 & 0\\end{matrix}\\right]\n",
      "\\left[\\begin{matrix}0\\\\0\\\\- \\frac{8}{375}\\\\\\frac{8}{375}\\\\0\\\\0\\\\0\\\\0\\end{matrix}\\right]\n"
     ]
    }
   ],
   "source": [
    "A,b = sym.linear_eq_to_matrix((eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8), (C_1,C_2,C_3,C_4,C_5,C_6,C_7,C_8))\n",
    "print(sym.latex(A))\n",
    "print(sym.latex(b))\n",
    "A = np.array(A).astype(np.float64)\n",
    "b = np.array(b).astype(np.float64)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "\n",
    "A = np.array([\n",
    "    [0, 0, 0, 1, 0, 0, 0, 0],\n",
    "    [0, 1, 0, 0, 0, 0, 0, 0],\n",
    "    [0, 0, 0, 0, -4/1875, -1/625, -4, 1],\n",
    "    [0, 0, 0, 0, 1/625, 1/1250, 1, 0],\n",
    "    [-4/1875, -1/625, -4, 1, 0, 0, 0, -1],\n",
    "    [1/625, 1/1250, 1, 0, 0, 0, -1, 0],\n",
    "    [4, 1, 0, 0, 0, -1, 0, 0],\n",
    "    [1, 0, 0, 0, -1, 0, 0, 0]\n",
    "])\n",
    "\n",
    "b = np.array([\n",
    "    0,\n",
    "    0,\n",
    "    -8/375,\n",
    "    8/375,\n",
    "    0,\n",
    "    0,\n",
    "    0,\n",
    "    0\n",
    "])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ 4.37500000e+00],\n",
       "       [ 0.00000000e+00],\n",
       "       [-6.66666667e-03],\n",
       "       [ 0.00000000e+00],\n",
       "       [ 4.37500000e+00],\n",
       "       [ 1.75000000e+01],\n",
       "       [ 3.33333333e-04],\n",
       "       [ 1.73333333e-02]])"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.linalg.solve(A,b)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "tags": [
     "remove-cell"
    ]
   },
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$\\displaystyle \\left[\\begin{matrix}4.375\\\\0\\\\-0.00666666666666667\\\\0\\\\4.375\\\\17.5\\\\0.000333333333333334\\\\0.0173333333333333\\end{matrix}\\right]$"
      ],
      "text/plain": [
       "Matrix([\n",
       "[               4.375],\n",
       "[                   0],\n",
       "[-0.00666666666666667],\n",
       "[                   0],\n",
       "[               4.375],\n",
       "[                17.5],\n",
       "[0.000333333333333334],\n",
       "[  0.0173333333333333]])"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\\left[\\begin{matrix}4.375\\\\0\\\\-0.00666666666666667\\\\0\\\\4.375\\\\17.5\\\\0.000333333333333334\\\\0.0173333333333333\\end{matrix}\\right]\n",
      "\\left[\\begin{matrix}1 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 4.375\\\\0 & 1 & 0 & 0 & 0 & 0 & 0 & 0 & 0\\\\0 & 0 & 1 & 0 & 0 & 0 & 0 & 0 & -0.00666666666666667\\\\0 & 0 & 0 & 1 & 0 & 0 & 0 & 0 & 0\\\\0 & 0 & 0 & 0 & 1 & 0 & 0 & 0 & 4.375\\\\0 & 0 & 0 & 0 & 0 & 1 & 0 & 0 & 17.5\\\\0 & 0 & 0 & 0 & 0 & 0 & 1 & 0 & 0.000333333333333334\\\\0 & 0 & 0 & 0 & 0 & 0 & 0 & 1 & 0.0173333333333333\\end{matrix}\\right]\n"
     ]
    }
   ],
   "source": [
    "import sympy as sym\n",
    "\n",
    "A_sym = sym.Matrix(A)\n",
    "b_sym = sym.Matrix(b)\n",
    "\n",
    "sol_sym = A_sym.LUsolve(b_sym)\n",
    "\n",
    "display(sol_sym)\n",
    "\n",
    "print(sym.latex(sol_sym))\n",
    "\n",
    "identity = sym.eye(sol_sym.rows)\n",
    "rr_augmented = identity.row_join(sol_sym)\n",
    "print(sym.latex(rr_augmented))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "C_1 = 35/8 ≈ 4.37500000000000\n",
      "C_2 = 0 ≈ 0\n",
      "C_3 = -1/150 ≈ -0.00666666666666667\n",
      "C_4 = 0 ≈ 0\n",
      "C_5 = 35/8 ≈ 4.37500000000000\n",
      "C_6 = 35/2 ≈ 17.5000000000000\n",
      "C_7 = 1/3000 ≈ 0.000333333333333333\n",
      "C_8 = 13/750 ≈ 0.0173333333333333\n"
     ]
    }
   ],
   "source": [
    "import sympy as sym\n",
    "sym.init_printing()\n",
    "\n",
    "x = sym.symbols('x')\n",
    "\n",
    "q = 10 \n",
    "EI = 5000 \n",
    "EA = 20000 \n",
    "F = 150\n",
    "\n",
    "q_AC = 0\n",
    "q_BC = q\n",
    "\n",
    "C_1, C_2, C_3, C_4 = sym.symbols('C_1 C_2 C_3 C_4')\n",
    "C_5, C_6, C_7, C_8 = sym.symbols('C_5 C_6 C_7 C_8')\n",
    "\n",
    "V_AC = -sym.integrate(q_AC, x) + C_1\n",
    "M_AC = sym.integrate(V_AC, x) + C_2\n",
    "kappa_AC = M_AC/EI\n",
    "phi_AC = sym.integrate(kappa_AC, x) + C_3\n",
    "w_AC = -sym.integrate(phi_AC, x) + C_4\n",
    "\n",
    "V_BC = -sym.integrate(q_BC, x) + C_5\n",
    "M_BC = sym.integrate(V_BC, x) + C_6\n",
    "kappa_BC = M_BC/EI\n",
    "phi_BC = sym.integrate(kappa_BC, x) + C_7\n",
    "w_BC = -sym.integrate(phi_BC, x) + C_8\n",
    "\n",
    "eq1 = sym.Eq(w_AC.subs(x,0),0)\n",
    "eq2 = sym.Eq(M_AC.subs(x,0),0)\n",
    "eq3 = sym.Eq(w_BC.subs(x,4),0)\n",
    "eq4 = sym.Eq(phi_BC.subs(x,4),0)\n",
    "eq5 = sym.Eq(w_AC.subs(x,4)-w_BC.subs(x,0),0)\n",
    "eq6 = sym.Eq(phi_AC.subs(x,4)-phi_BC.subs(x,0),0)\n",
    "eq7 = sym.Eq(M_AC.subs(x,4) - M_BC.subs(x,0),0)\n",
    "eq8 = sym.Eq(V_AC.subs(x,4) - V_BC.subs(x,0),0)\n",
    "\n",
    "sol = sym.solve([eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8],(C_1,C_2,C_3,C_4,C_5,C_6,C_7,C_8))\n",
    "for k, v in sol.items():\n",
    "    print(f\"{k} = {v} ≈ {v.evalf()}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "::::::{codex-end}\n",
    "::::::"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Meer voorbeelden\n",
    "Het voorbeeld in hoofdstuk 2.7.3 van het boek Mechanica, spanningen, vervormingen en verplaatsingen {cite:p}`Hartsuijker2013` richt zich op extensie. Voorbeeld 3b in hoofdstuk 8.2 is een buigingsprobleem.\n",
    "\n",
    "## Instructies in collegevorm\n",
    "\n",
    "Dit onderwerp is [in les 3](https://collegerama.tudelft.nl/Mediasite/Channel/public-ceg-ctb2210/watch/757766716fcb45718fca678dadaa396d1d?sortBy=most-recent) gepresenteerd in collegevorm van 0:40:00 tot 0:48:30\n",
    "\n",
    "## Opgaves\n",
    "\n",
    "- Opgaves 2.48 en 2.50 in hoofdstuk 2.9 van het boek  Mechanica: spanningen, vervormingen en verplaatsingen {cite:p}`Hartsuijker2013`, antwoorden zijn beschikbaar op [deze website](https://icozct.tudelft.nl/TUD_CT/boekantwoorden/vol2/Chapter2/).\n",
    "- Opgaves 8.7, 8.10, 8.11 in hoofdstuk 8.6 van het boek  Mechanica: spanningen, vervormingen en verplaatsingen {cite:p}`Hartsuijker2013`, antwoorden zijn beschikbaar op [deze website](https://icozct.tudelft.nl/TUD_CT/boekantwoorden/vol2/Chapter8/)."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "base",
   "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.12.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}