{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Mechanics Diagrams"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%load_ext tikz_magic"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Quarter Car Model\n",
    "\n",
    "via https://tex.stackexchange.com/a/53885"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfoAAAOECAQAAADjqztaAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAAAmJLR0QA/4ePzL8AAAAJcEhZcwAAASwAAAEsAHOI6VIAAAAHdElNRQfhDBYMAAZBN9uUAAAmMElEQVR42u3dP3IqXfTe+2e53sSBA3wTVzlDnsDPrfwmzRA4Q0BDgCG0hoCqnLvEEGAI4t5y5sCiyoETV/kQOHC4HYA4jYR0kIC9dvf6fqh63wNSN5umH+3Vf3a3JQGI5F95NwBAXoQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEPigbercAXgh9VGOrvZsAH4Q+qkqVdxPg4x/vBsAJkQ+L0IdklYYaeLcCPijvY6olDYzePiRCH9Mu7uzKC8m42UVE+699ne69W4L86OkDsvH+H5WxXR8QoY/oT1k/vmAu6ChCH9GfqLNVHxChD8fGrYN1Ywr8eAh9PMe9OwV+OOy9D8YGej06LYc9+OHQ00czfncmXsV4u2gIfTQfy/mJd5OQF+V9KDbU64cXN+nOu13IiZ4+llO9+tDo60Mh9IHY4JNSnj34oRD6SMafDKetGW8XCaGPZPrpTyjwAyH0YdhEnx+cm3DgLg5CH8fXW+7TM+eCzuOQXRBWa/mXX/m3aevdSuRATx/F33ty+vog6OlDsLGez/i1u7Txbiluj54+huas36KvD4FLYAdgAz0cvfCo/7j/13/Wf/JuHXIj9AGkrVbt5/b78M//nlbfnh06jvIeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgjGkncLimFDTbzbkMn/q3+3/9d/1X/xbkweaebdgnIQ+gOrtfRuA24lmXcLyvGPdwMKs9HCuwm4uql3A8pC6I9tKAP7xwj9EXbkAcEQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEKPL9jAGht4twLXRejxKav1oqkq73bguv7xbgBKZJUqjVV7twO3QE+PI1ZbsqQXTbX2bgtug54ex9YaSdqkjWRT78bgFgg9jqStVt5twG1R3gPBEHogGEIPBMM2fU/YUENt09Eed6s0kLRO26NXB6qk97+LOOjpe8EmWmqspS1358/ZwBr7rUa1Gv22+dtZdVbbUi+aaKylvdrYu91wkXjsH6qVtPRuxQ9b/qpBkqa7T6BKL2o02P+0UdLL/l8vqvevDvVbSeMv57v7X+39+S5ePknJuw0lPejp+2Cqx7SVtJZUW61nzdLsUNQvJFU2tkYDjdL+gFzaaC1p4t105Mc2fefZUFUaSZKGkqRnzVLrWHtamyRNJY2Otu7XqjnRNiJC333jw+k0u9A/paf2j203YKbS3fEOPQbSREV5333V4Sz5Xb/9+O7nb38KNide3wjh0NN330oL6XAobvGuP3/r0RfHL9pQQ4lBNRER+s47FPO7fv5j311LUlqdepXQR0R53x+7Hv1duPf9/8dBNLvQL/42U/QPoe+P7/XotaR1Yps+IELfE29b9B9+cLr/n2gg+vmgCH1ffNWjf9b/70NvHLwLhdD3xTe26G2gsQ7FvQ304t145ETo+6LWqZFzp/v/3cm3b3v9x5T5sRD6Xtj36B+L+5P9/7sj92M9CYEQ+n743hZ9a5S9jU/8HL3GyTn9cPqsu8+O0S9Uaxf5So1+vZuq3s/xbfdeY4v9n5P1h7P90EGEvh+GkjYftuh3Z91/2GJPT1ZpYo2ksX59mGrZ+vdKVSv+I66U2weEvh9m2nzcMk9re2idptv+yYM9aayN7j/23cm8Pwxuy5J3C4phtZZa7Uemo0cs8aesjR15QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4I5h/vBhRmYLV3E4DbIvTHKi29mwDcFqH/Y6uVw7v+B0nSb/32/vgZ/Hv9a0n/R//DuyGxWfJuQXD7L+AxzbxbkuGzLlVLWqWRd0tiY0ceEAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0LuysXcLEA+hd2NDW+rZuxWIh9C7sIE1elV9eGHt3SLEQegd2Fgvmh6ernSXFt5tQhz/eDcgGhtq3urhN5oReORFT5/Rh6L+UfdEHrnR02djYzUaHp6u9JA23m1CRIQ+E3tW+/DcLD16twhRUd7n8qCn1rOpNd4NQlSEPpO0TQ8aHQ7NDTS1V6svmiXwI4Q+o7RK95ppu3861NKebXjRLIFvI/SZpUfd6c/2/FgvFPrIi9Bnl7Zppnut9k8p9JEZoXeR1mmkBwp9eCD0btIThT48EHpHFPrwQOidUegjN0JfAAp95EToi0Chj3wIfTEo9JEHoS/LRu2Rd5UG3g1C/xD6YtjA5lqq2j/d6jHdJS6jhatjaG0hbKppq19faMZoe9wGoS+A1WoOPby00UNaXTA74EuE3pkN1GhyeLrVU5p5twn9RuhdUdQjP0LvhqIePgi9C4p6+CH0Dijq4YnQZ0ZRD2+EPiOKepSA0GdDUY8yEPpMuNkFSsG597k8qB3zxhpjMA1cEPpM0jbNdKc/O+2merXJBTMEfojQZ5Q2aaRfh8GzA81tadVFswS+jdBnlha6bxX6tV4o9JEXoc+OQh++CL0LCn34IfRuKPThg9A7otCHB0LvjEIfuRH6AlDoIydCXwQKfeRD6ItBoY88CD0QDKPsimFDzfXn/nVbzdKTd5vQR/T0RbCBNXptRf5Rd0Qet0FPXwAbq9GfW1WuNON2VrgdQu+Moh65Ud47oqiHB3p6NxT18EHoXVDUww/lfXYU9fBFT58ZRT28EfqMKOpRAsr7TCjqUQp6+lzm3OwCZaCnzyT90kzbw9PGnm14weyAHyP02aRH3elPQT/WizXebUJEhD6jtE0PGultb/1AU3u1+qJZAt9G6DNLq3TfKvSHWlLoIy9C74BCH54IvQsKffgh9G66VuhbY0tb2tK1DUtb2pJLhl6G0LvqVKFfqVYt33pk14KC/zR2QY9OzrFK1X51qLpzRnva6sEWarS77u1AUxvrIa0unC2+xYYa7v+cDbXq+5mSPQq9mlYv9O1bRdhAUz2m7Xenu4a00r1NNd23eqilLTRLmwtni/NNND38+8vl7rmeXEuPyvs0SqZf+tHXYbVeNJXjVeY7Vej3TpolO1r+n/BfT66hR6GXpLTQNwtjq2xiSy39txNP7tHv+OrVJWmjxec/LWc9uVzPQq+/FGdtVluypBdNVc72f3W0Uq1/VrfgmopcTy7Sp23671prJGmTNpJNL57bxaw+7M6TpE2XdudZ+sYv12f99iqNvD/VXmHryeUChz5tv7spcDs2UKM/R5+3ekoz7zZhp6T15DoCh74crT33ksSee9wUoXfW5aK+5fqXBOHP3s0Qekf9Keq72u6YCL0bivoS7IY5dbK6+rHeh94qHZ/kUsQJuj0p6jvNak001lprDWyulcJUK70PvbbS4fTchTb+R777U9R3lw0011hr3b1VVzbVy9/PyOuH3oc+bexRtaQnPZZQPlPU+7OBlqqOzwRIj7bSi3fL8ujfGXkf1dpqlB78w2W1vag5RH6jUfrl36qAnlVpq1/HL6b1DY5BFKn3obdGtUb+W8w2sLmWh+34rR7TnX+rIrLJrvI7MVIuyPfR8/Le5qo0KmIg5PHNLh652YWb3am0iwvn0mE97ultYC+qC4m80i/9ap1w0tiy5Etj9ZfthzSVcAzHS29Db5WWqjQsZyhkWui+tdVY69Ua+/bFPnChXb1VREfgpa+hr7Tcb6HNywlW2qaZ7lpbjlO92PiCGeKnAvfz/d2mn+ohLUyaqtK0pNMu0kaj1j3qh3q2lQo4snA+LtXddX0N/SitpTSzSrWmtiprP3la2ErTw1XZar3aY4euu+Z6EewrKWajz0NPy/vDbpoHbSU9l1Pi79tHoe9lt2YQ+v5KGz1IGmju3ZJTbUuj1h79oZ7Zo397abHbiRd5Sfe1vD9IC1torLFNSzwy3qlCf/b9C4u/86j/KOn/06UXnbpsD8iTppLqKGfan5B69lCjpHT0ykC/lfRb1RdT7f5XO7V5qGXr6avG3kvxRp9z9ymXzq3YrQ+vJ34yV1JSU+p6cq1H/8r7D2Vb2upJ0qCkg3fvWkihn1HaaqSthu/vK2CT/TH83i/5HpX3VmmgajeM1p71JGmdtpK93X+t0os9ab17VTocfKoOZ8Q3ttjv6FnnLrE7Veh3XFrbSEtNrXXFH2s01qMaSWNrtNL2zzl7Ja0n11kAvXkclcitMuzw7PV9cXb0u0v9fj+lw2cY6qX19LXrheTJb8i5vN+3ZaC5futVjRo1etVcA9Xt9aH1u8WtJ5c9etTTf3ad9GSfTvHpT9xMjorLRewzx24pbfWgh30VuN7VVLbR7FQPXuB6cpFv3aUAt2S15q3Ir8q4sNdVP+FStUq6jUVQ/duR10k2tOfWfdK2ekijvkUepehRed9d1mjSOgbODjzcFKF31v+iHqWhvHdEUQ8P9PRuKOrhg9C7oKiHH8r77Cjq4YuePjOKengj9BlR1KMElPeZUNSjFIQ+l8nRzS6eUtxLOMAZoc8kzTRqDZ+Z2gtXlYUPQp9NWqV7zQ63Wai0tGIv64E+I/RZpUfdta7NNtGrXXq9OOCbCH1maZseWoX+QA2FPvIi9A4o9OGJ0Duh0IcXQu+GQh8+CL0rCn3kR+jdUegjL0JfAAp95EToC0Ghj1wIfUHS49Htqyd65fbVuD5CXxAbaKo/Zf1WT2nh3Sb0D+Ppi2ETNa3Layw0S5fdkhk4idAXwSo1rT5+o4e0umB2wBcIvTsbtO5WuyvqZxfMDvgLQu+Moh65EXpHFPXwQOidUNTDC6F3QVEPP4Q+O4p6+CL0WVHUwx+hz4iiHiUg9JlQ1KMUnHufy1jtwbILcXcbOCH0maTZ0Qi6qV4YQQcfhD6btEkj/dLbVvxQz7a04UWzBH6A0GeVFrrX4+FprVdruFQG8iL0maUthT58EXoHFPrwROidUOjDC6F3Q6EPH4TeFYU+8iP07ij0kRehLwCFPnIi9IWg0EcuhL4gaaF7/bnSfa1Xbm6F6yP0BbGh5vpT1m8ZiYdbYGhtMazRpDXa/lGPaXvB7IBPEPoiWK25/mzBrzRLDL3FjVDeu7OhPWt5iPxWD2lE5HE79PTOKOqRG6F3RFEPD5T3Tijq4YWe3gVFPfwQ+uwo6uGL8j4rinr4o6fPiKIeJSD0mVDUoxSU97lUao+R34gbWsEJoc8kPepOT4enE73a9ILZAT9G6LNJ2/Sg0eF2VgM19sLQWeRH6LNKq3Svmd5231Va2pxLYyEvQp8dhT58EXoHJwr9xrtNiIPQO3lX6APZEHpH6VF33m1APITeFWfkIT9CDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBWPJuQTGs1tK7DbiVZN4tKMc/3g0ozFbrzO/4L5Kk/6n/4f3RM/gP+jeS/rf+W+b3rb0/eFkI/bF1Gnk3AddGOXuMbXogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwTzj3cDUCarNJYkbbVKa+/W4JoIPT6wWnMNtZI01FCyjZ7So3ercC2U93jHGj1roX+bRmmU7jTSWkM19mID75bhOgg9jthEE43SLG13z9NKI60lVVoS+34g9GixgRo9Hm/Dp61+SZIqTb3bh2sg9Gibavtx6z1ttJAkTejr+4DQo63S0JItP7y+kiQNVHs3EJcj9Gjbhbq2ybvXN/v/V94NxOUIPdo2l88CpeM4fU/YUENtj3fBWaWBpPXbnvj9qwNV0vvf3XvUXNJ6vw2PXqKn7wWbaKmxlrY/rGYDa+y3GtVq9NvmbzvgrLalXjTRWEt7tfH7+aQn3WmU7o//TOhPWc+5eX2QeOwfqpW09G7FD1v+qkGSprtPoEovajTY/7RR0sv+Xy+q968O9VtJ4zPf4VlJSb/f5tmth5KSdxtKerg3oJxHh0O/1OTwCZJqvb5FOylJlZKSxmo0b4dWy3M/r4b7fzTen/SHyycR+vaDbfrOs6GqNJIkDSVJz5ql1Z+fp7VJ0lTS6KhsX6s+8xDc7qScdZp5f1ZcA6HvvrHeIr4L/VN6av/Ydtvjle7ebamfefjNKk0kvZ2Xh85jR173VYfda7t++/0ZdW9/CjYnXv/rITobaC5pq1HicF5P0NN332p3gG1/KG7xyZ73dwfhbKihztkb36jSViPG1PcHoe+8QzG/6+c/9se1JLW38lu//Zco20QTIt83lPf9sevR34V73/+vPvz2LvRfnoRjE8211h2R7xdC3x/f69FrSeuvttOt0lzrd3v80QOEvifetug//OB0/z/RQF/281Zpqafjc/NsaEPvz4nLEfq++KpH/6z/34fePhy8s6GWWqSHdy9PNBE6j9D3xTe26G2gsQ7FvQ308uHnzyciL1Un9g2gcwh9X9Q6NXLudP+/66/f9vqPj8t8G2ip9YnISzVDb/uAQ3a9sO/RPxb3J/v/d0fux+9O52kkLezjCbq1xAk6fUDo++F7W/StUfY2Pv65zTWRtDz5Lhy66wXK+344fdbdZ8foF5J2ka/UqDWMxpovd9XRz/cCPX0/DCVtPmzR7w6wfTgwl56s0sQaSWP9Oprq64tcE/peIPT9MNNGT+9fTGt7aJ2m2/7Jgz1prI3eXSMnmfcHwe1Z8m5BMazWUqv9yHT0iCX+nLWxTQ8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRDMP94NKExtybsJwG3R0wPB0LE5238Bj2nm3ZIMn3WpWtIqjbxbEhs9PRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtC7soF3CxAPoXdkE716twHx/OPdgKisUqPauxWIiNA7sIGmmh6ebvWUZt5tQhyEPjubqNGfbfmFZmnj3SZEQuizelfUb/SQVt5tQjSEPhuKepSB0GdCUY9ScMgul622rWccn4cbQp9JWuhej4entV6t4dQceCD02aRtmulOf3bcTfViY+9WIR5Cn1XapJF+6W1rfqhnW9rQu1WIhdBnR6EPX4TeAYU+PBF6JxT68ELoHVHowwOhd0Whj/wIvTsKfeRF6ItAoY98CH0hKPSRC6EvSNqkUSv2Qz0Te1wfoS+I1fbaGm2/0n1aeLcJ/UPoC2FDe9ZSbzvwtnpIo7T2bhX6iPH0RbBGk9Zw20c9pu0FswO+QOjdWa25/hyiW2lGD49borx3RVGP/OjpHVHUwwOhd0JRDy+U9w4o6uGJnj47inr4IvRZUdTDH+V9NhT1KAM9fSYU9SgFPX0uK7XvaFOr8m4QoiL0maRVutfscJebSkubM2IeHgh9RulRd3o6PJ3o1aYXzA74EUKfVdqmB430tvtuoMZerL5olsA3EfrsKPThi9C7oNCHH0LvhEIfXgi9Iwp9eCD0zij0kRuhd0ehj7wIfREo9JEPoS8GhT7yIPQFSVvNWje7GGjKzS5wfYS+IDZR+2YXC252gVtgaG0hrFLTCvxGD2l1weyATxH6AthAU/3Zft/qKc2824T+IvTubKKmdXmNhWZpc8HsgL8g9K4o6pEfoXdDUQ8fhN4JRT28EHoHFPXwROgzo6iHN0KfFUU9/BH6bCjqUQZOw83EGr20Iv+Y7og8fBD6XJ7UDvmYoTTwQugzSZs00q/DXW6GeralDS+aJfAjhD6jtNC9Hg9Pa71aw6UykBuhzypt00x3rUJ/qhcKfeRF6LOj0IcvQu+CQh9+LHm3IDAban505H5W8pVyrNqdVvTzQ422VC1plUY/nsNuaW04pekiiYfrQ2O9tp4uNfRu0actXe7+cfEclhfMYfe/xntZdPvRozPyrFKl3bZxpVlaXzi7TNLCVq2z8Wu92qMe0/aimeKbbKjhvuYaapWeLpxd4XoU+qOTXL+xfWyVdvvPt1p5/KlIW83sqVXoTzW2ogv9Hpq0BkF9sungvZ5cT4925KVRMv3St/pIq+1VL6pUaaxGL+Z0rXn26PtKs2RHdx14p5T15Foft2cPPSspqT7rdxv9VqPB/lmtFyUlvby9kr3tAzVHLzReLTnZut5v06s+PffS1pNLHz3q6ffO3q9rE000SrO37ee02t9RrtLS5/BZ2r47Q38i7mnnrrz15FJ92qb/Fhuo0ePxtlna2i+9Sqo0lcOlLazRpLU34kmzruzQs+VZv/Yvu/+e9dvrMi4uUuJ6cjHvUuPqJdquQP5rea9Grydf320e/M5duqk+OnT3ct4GStYWflHe3+Dtll+8T8byvrT15BqP/pX356o0tHSiz9mV14OchbUN7VlLve2422qW7hltX4iC1pNrCVve77+s2ibvjsq+7ROolOmgWXeL+pbr/4kq5aBYMevJ9cQN/UYFHBKzWvNWO9aadbOH//mJtcUrYj25rrjl/W64y9rv7zRFfQmsttq+KtHd15Pr631Pb5Waoxf2J+imJ1tpeCJk1f7/Ny4ve1HUd5rVmmistdYa2Fyr0/vhvdeTW+h96LWVDjtbFtr8OWMvbU4e06/2U92wz+1LUd9dNtBcY6119zZez6Z6OX1Gnt96cjPehw+uftjlwyG7/YGY+Tnj1zTc/+Nm47g03B/s2T1+a+q9xM5u+cVn5F2hDVf5djTQiz4cFFR1/txvv57c8hFhm77WVqP0cNYY7N0Z1Tc7McQavejP5bGedJceL5gdfuZZlbb6dfxiWuv87+LG68lt9T701qjW6Lzy2SpNpI+rw5VaUturpoft+LVG6YHt+Pxsd3Lz04llf2apftv15PZ6vk1vc1UanRctG2guaavRLa7LYnNNWk8fu9lH9MKul/7x3vjbric59Dj0NtBSA92f3Zs2qrTV6EZjpR81aBX2E9tS2Huw/YVWLviWb7ueZNDb8t4qLVVpeO6pFTbR5JZfZdqkXxod9gMP1NiLdfAUzs57uxDGD916Pcmhr6GvtNxvoc3PGf5oE8211t1tv8q0Snd6PKxwlZZ2VutwdT/8nvOsJ7fW19BP9ZBmetRu+ONfWKW51udu+18mzXTf2p6cqNvXYAkl53pyS30N/SgtpDTTStL06zLaKi31lI62/W14u4tVdb3Qt+T3uNqH+MG3m3s9uZ2ehv5Qfj1oK+n58yLahlpqkR7evTw52td+/fZR6HvZrRnfjqrPenIbPQ39m7TRg6TdQZYTbKDnE1+lVN3+9EoKfQ9psftT+70e2nM9ub4eH7LbSQtbaKyxTT8eIrOBllqf+CqlWg9nzPzStm30q3Ue/kCNjQs+D//ynVf/ov9H0v/S/3/hfC47Pv6kqaRaZ1/d3ns9ubbe3dbKGk2lZK1XBnrV4NRhFpurOjm2qta0PYcMbQ4x4u7y21pdpRW79WGT7j78ZHcC1YcTp0pZT67G++T/qw/JeFZ6PyhkPwjn3SWLNf9iNi+ZW93ZYTjf+pQXXwL7Su2o9FsfBstoot9KSnp+93pB68mVPr93A676Vdaavn1xqlXvQn64TnnSq6aHV5svZ/Z827aebH/xl8a8+BMWEvpTsVejV033TxrVqg6vF7aeXOHTezfgih9l+eGl+ugjvp589fTDacikmv0frd1j3sVrrf71Gyog9Gm3c/e3XtWoUaNXzTXYD8JO7VaWuZ5c9ujdNn23WaX54Zos0kYPxe7W+8mnK2Kb/qhFtWpJa63SVrKhxvvdlet+7lXZf2pCXwobtO5eK0m9u3tteaGPqfeH7LrCJmpae/BXOu+iH8C39fzknG6wypaaHyK/0a/U2bHaKB89vbP+F/UoDaF3RVGP/Cjv3VDUwwc9vQuKevgh9A4o6uGJ8j4zinp4o6fPiKIeJSD02VDUowyEPhN7bl31Xppx1Xt4YZs+l4ejO6U11nBNPPgg9JmkbZrpvnVFtalebHzBDIEfIvQZpXUa7a/PK0lDPduyi5dQRrcR+szSk+5ahX6tVwp95EXos6PQhy9C74JCH34IvRsKffgg9I4o9OGB0Duj0EduhL4AFPrIidAXgUIf+RD6YlDoIw9CX5aN2kNtB6LIx9UR+mLY0J611FvfvtUs3afLbw8NvMPQ2kJEuV01/BH6Alituf5sva8169Md7FAayntnJ4t6Io8boqd3RVGP/Ai9G4p6+KC8d0FRDz/09A4o6uGJ0GdGUQ9vlPcZUdSjBPT02VDUowyEPhNbqm495WYXcEN5n0ka6bE1mIabXcANoc8mzXSvxeHpVK828W4TIiL0GaVN+qWR3m5bOdDcllZ5twrREPrM0irdtQr9Wi8U+siL0Dug0IcnQu/iZKFfXzRL4EyE3s2HQn9pcwp93B6hd/Wu0J/o1abebULfEXpnHwr9xl4o9HFLhL4A7wr9ikIft0ToC0Ghj1wIfTFOFvqcuoOrI/RlOS7qt2IcHq6OUXbFsKHmrZF4G83S4oLZAZ+gpy+CDazRayvyj7on8rgNevoC2FhN6xJaKz2kzQWzA75E6J1R1CM3yntHFPXwQE/vhqIePgi9C4p6+KG8z46iHr7o6TOjqIc3Qp8RRT1KQHmfCUU9SkFPn0uj9nXwHtKTd4MQFT19JulBD63hM3N7tuEFswN+jNBnk550pz83sxrrxRrvNiEiQp9R2qaZ7vV2n9qBpvbKpbGQG6HPLK3TqFXoD7Wk0EdehN7BiUJ/7N0mxEHoXXwo9LksFrIh9G72hT6QGaF3xdF65EfogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QDKEHgiH0QDCEHgiG0APBEHogGEIPBEPogWAIPRAMoQeCIfRAMIQeCIbQA8EQeiAYQg8EQ+iBYCz9cIJk5/6kW1Pkfm+rJUmbtCmpVTeaU6WX8j5f+evIpVN8Nocz/e0Nv7/gypyiW61libBEzo/8N8v7bn3cvn+dLBGWyN+mOO0boe/Wx+3718kSYYn8bYrP/HPer5W6tcL2G0vEf4oyW/u5M3v6bn3cvn+dLJGSpiiztV85K/Td+rh9/zpZIiVNUWZrv3ZG6Lv1cfv+dbJESpqizNb+zV9D362P2/evkyVS0hRlt/ZzZ27Td+vj9v3rZImUMEXZrf3KWaHv1sft+9fJEilhirJb+7UzQt+tj9v3r5MlUsIUZbf2b879vY583L5/nSyREqboVms/+uYou6593L5/nSwRlsjfpvjoR0Nrv3/woCtTdKu1LBGWyHmtPfajATfnN6WbU3SrtSwRlsj3ttKvOOCmK+VO38s5lghL5Gtnh75bH7fvXydLpKQpymzt5xhwU+x7lzlFt1obd4l8hQE3hb53mVN0q7Vxl8jXGHBT5HuXOUW3Wht3ifwNA24KfO8yp+hWa1kin2PATXHvXeYU3WotS+QrDLgp7L3LnKJbrWWJfI0BN0W9d5lTdKu1LJG/Off3OvJx+/51skRKmKJbrf3o7KvhnveG319wZU7RrdayRFgi39mTz80uCnvvMqfoVmtZIl/jZhdFvXeZU3SrtSyRv+FmFwW/d5lTdKu1cZfI5zj3vtj3LnOKbrU27hL5CufeF/reZU7RrdbGXSJf49z7It+7zCm61dq4S+RvOPe+wPcuc4putZYl8jnOvS/uvcucolutZYl8hXPvC3vvMqfoVmtZIl/j3Pui3rvMKbrVWpbI35z7ewB64kfXvQfQXYQeCIbQA8EQeiAYQg8EQ+iBYAg9EAyhB4Ih9EAwhB4IhtADwRB6IBhCDwRD6IFgCD0QzP8FAHYLeMW9HHAAAAAldEVYdGRhdGU6Y3JlYXRlADIwMTctMTItMjJUMTI6MDA6MDYrMDA6MDDDbIibAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE3LTEyLTIyVDEyOjAwOjA2KzAwOjAwsjEwJwAAABR0RVh0cGRmOlZlcnNpb24AUERGLTEuNSAFXAs5AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%%tikz -p circuitikz,calc\n",
    "\\usetikzlibrary{calc,patterns,decorations.pathmorphing,decorations.markings,circuits}\n",
    "     \\tikzstyle{spring}=[thick,decorate,decoration={zigzag,pre length=0.3cm,post length=0.3cm,segment length=6}]\n",
    "     \\tikzstyle{dampener}=[thick,decoration={markings, mark connection node=dmp,mark=at position 0.5 with {\n",
    "     \\node (dmp) [thick,inner sep=0pt,transform shape,rotate=-90,minimum width=15pt,minimum height=3pt,draw=none] {};\n",
    "     \\draw [thick] ($(dmp.north east)+(2pt,0)$) -- (dmp.south east) -- (dmp.south west) -- ($(dmp.north west)+(2pt,0)$); \\draw [thick] ($(dmp.north)+(0,-5pt)$) -- ($(dmp.north)+(0,5pt)$);}}, decorate]\n",
    "     \\tikzstyle{ground}=[fill,pattern=north east lines,draw=none,minimum width=4cm,minimum height=0.3cm]\n",
    "\n",
    "     \\begin{scope}\n",
    "     \\node at (0,0) [draw,rectangle, minimum width=2cm,minimum height=1cm,anchor=south,,transform shape](m1) {$m1$};\n",
    "      \\draw [very thick, -latex](m1.north) -- +(0,1);\n",
    "     \\node at (0,-3) [draw,rectangle, minimum width=2cm,minimum height=1cm,anchor=south,,transform shape](m2) {$m2$};\n",
    "     %\\draw [very thick, -latex](m2.north) -- +(0,1);\n",
    "     \\draw [spring] (-0.5,-2) -- (-0.5,0) node[midway,left=0.3cm] {k1};\n",
    "      \\draw [dampener,label=D1,] (0.5,-2) -- (0.5,0)node[midway,right=0.4cm] {d1};\n",
    "      \\node (ground1) at (0,-5)  [ground, anchor=north] {};\n",
    "%        \\draw [ground] (-1,-5) -- (1,-5);\n",
    "      \\draw [spring] (-0.5,-5) -- (-0.5,-3)node[midway,left=0.3cm] {k2};\n",
    " %     \\draw [spring] (-0.5,-5) -- (-0.5,-3)node[draw=none,midway,left=0.3cm] {k2};\n",
    "       \\draw [dampener] (0.5,-5) -- (0.5,-3)node[midway,right=0.4cm] {d2};\n",
    "\n",
    "      \\end{scope}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Trolley problems\n",
    "\n",
    "via https://tex.stackexchange.com/a/342365"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABOQAAAGpCAQAAACHNiGCAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAAAmJLR0QA/4ePzL8AAAAJcEhZcwAAASwAAAEsAHOI6VIAAAAHdElNRQfhDBYMAAimj/aTAAAsl0lEQVR42u3dzY4ryZne8ScGZylgkN0rr2bAbl1AD4+1NSCwdAc8I68NsKydvRiwRhcgJe2NdxoW7L26CF+ATGIuQFO09u4uwhcgFWFg9uFFZvIzszIimWRGJP+/RPepKjKqgm8xiy/jI19jBQAAgBj9VdcdAAAAQDMkcgAAAJEikQMAAIgUiRwAAECkSOQAAAAiRSIHAAAQKRI5AACASJHIAQAARIpEDgAAIFIkcgAAAJEikQMAAIgUiRwAAECkSOQAAAAiRSIHAAAQKRI5AACASH2qu4MZaNDWD7Orrh8uAABAf9Qmcppo2tpPM10/XAAAgP6oT+T2TsfTRhVfL79lJAAAALTIOZGzJ6NpxpZ/veqW4qsAAABoh+Nmh0vTOAAAALTNLZGbHX/qn8YxHgcAANC2BpcfIY0DAAAIgXciRxoHAAAQBp9dq3JJ1kjjAAAAbsNrRK55GsfWBwAAgLZ5JHKkcQAAACFxnlqtS+O4EAkAAMBtOY7IkcYBAACExi2Ry6utksYBAACEw2ONHGkcAABASJwTOdI4AACAsFBrFQAAIFLUWgUAAIgUtVYBAAAiRa1VAACASHkmcqRxAAAAoWgwtdpkiwNbHwAAANrmlchVja7Vj9MBAACgbR6JHGW6AAAAQuKcyJHGAQAAhMUxkSONAwAACI1bIjfN/iGNAwAACIfHGjnSOAAAgJA4J3KkcQAAAGFxTORI4wAAAELjlsjNjj+lvgMAAED3GlR2II0DAAAIgXciRxoHAAAQhk9+d69P1kjjAAAAbqNRrVX/NI6tDwAAAG1rsdYqaRwAAMAtOU+tUqQLAAAgLNRaBQAAiBS1VgEAACJFrVUAAIBIUWsVAAAgUtRaBQAAiBS1VgEAACJFrVUAAIBIUWsVAAAgUp6JHGkcAABAKBpMrTbZ4sDWBwAAgLZ5JXJVo2v143QAAABom0ciR5kuAACAkDgncqRxAAAAYXFM5EjjAAAAQuOWyE2zf0jjAAAAwuGxRo40DgAAICTOiRxpHAAAQFgcEznSOAAAgNC4JXKz40+p7wAAANC9BpUdSOMAAABC4J3IkcYBAACE4ZPf3euTNdI4AACA22hUa9U/jWPrAwAAQNtarLVKGgcAAHBLzlOrFOkCAAAIC7VWAQAAIkWtVQAAgEhRaxUAACBS1FoFAACIFLVWAQAAIkWtVQAAgEhRaxXAjhkZa5Zd9wKAG7M01oy67gW6Ra1VAACASHkmcqRxAAAAoWgwtdpkiwNbHwAAANrmlchVja7Vj9MBAACgbR6JHGW6AAAAQuKcyJHGAQAAhMUxkSONAwAACI1bIjfN/iGNAwAACIfHGjnSOAAAgJA4J3KkcQAAAGFxTORI4wAAAELzyeles+NPqe8A9NRWK6277gQAR2tJ2647gW65JXJHSOOAvrJrPXTdBwCu7FPXPUD3vBM50jgA98gMNNDWMl4JROF+zljPWqvVyRppHID+MqmWGmtpXk1ycktiXszL6VcBdOmezthGtVarkzX3TREAcCmTmsFNfs5cif3GPupZQ81PbpxqrLGGXccCCJ0ZmumNftJdnbEt1loljQNwY0O9mfm1kzkz0cA+7j4dn/VBEltEgFqJUvNmJtf+Mfd2xrZYa5U0DkAHJldP5qa7nftlP2UkaW3ZOQi4GGhulmZ8+Tf6wJ2dsdRaBRC/KyZzZqKtXUmSSbKXgKNbR5Kk1Vmr4ZVfqoB4jfRyvZE53zPWJGZkUpOaUazr5qi1CqAfsmRu1Pr3HWmRfzRWIu0+K26VTl8qEjPVUlefQAIiNtDcvJnpFZInjzPWJCbVm6ZKNNRS7+blNqtu20WtVQD9MdHSvLb8Tn+l5/yjbIzt+ejWYX4fSZKZmLlZ6l0TRfreHrihgVK9tb5lyfmMNYmWGugb+2Af7YO+0kpjvZrotkFQaxVAvww1N+8mbevPsX3OVtOYgUaSFidra07X22w101cnYwAAqiSa6s0s23v75XHGTrWxX4rb7dY+aKNEy9imWKm1CqB/Ek31al7NpMU/ydm7+6PVcKfrbeyzfbKrPi2jBm5ipLl5N/NWR8Pqz9iJkpPFGDNJSWzLIqi1CuCASaK6vtJXH9461Fxzs9Dq7F15E9nLQu0KOeCWzDCiafzvPrw10UQTs9FCz3bTwk+rOWPNUIlGGh79FVloLml0mvWEjVqrAHbMSMuu+9C6scaXp3N5gns62na0Qg64NbNU+9t7ujXQVFOz1kqLS8prOZyxG22VHL8Ns9sYJxM9S3SRxgGI0lhzvZvXxmvnyi8i2rsrUgFBGGqqV/Nm5mbccHlE7Rlrt/qsJ305vDn/69DGeOANeY7IkcYBiNhQQ03NViuttbY+I2kl14uruoYcgFYMNNFEMmuttJbf6lOHM9ZuzqZQy6Zjg9dgarXJFgdrSOgABCLRWEOtzcg+XfR9WCEH3MJQQ600NCuvN1/nas5Yk21zuPSn3JxXIleVjNWP0wFAELZaaaVVe4upbWTv3oGorLVqLbWqO2PnSrQ+nmyNgUciR5kuABG7bPn0oqhwUzBzDcTEKnAdWy20vmCDkvcZa6Yaa62H+Na8OidypHHAnfiL/tR1F5x9p68d7rXWQovLxuDs2qw00mRXw2Ger6ZhYhVd+99677oLjr7S3znca6uFVpeOdPuesWaiVAs9xpfGOSdypHHA3fiTfei6C65qL77Q3jWppC+aa2yWWksaa6FnTcWIHLo3jWVNl8PljZ4vT+F2PM5YM9FcTzaqq8ftuSVy+QAlaRyASGy10KLNFzi71Rcz0EgDrTWzW/MqSbG8hAKBW+u5lQt377ifsWaiVF+KBNIkGlxyBbvb81gjRxoHIAobzdp9SSjYTVGA2ww01Gk5bgD+tlro+Tqpk8sZayaa6uHg5480jmvDg3MiRxoHIAILPbc9SmYSpRpoYQ9fBCb5TwPQ3EbPem77TZfPGWsmmujzUQ+GPb0gMGkcgOA9a9bSerhjqSaShvt38/n1pmZMrAKNrU5SrfY4n7Fmorlmmh4lLBM9dh0aP26J3MkCQOo7AD210Sy2d6OSpK2eW9vUcC7b7Xb4EpBdb6piabQZaiRpaMZcYw5XttA6yjN2ddW3QY5nrJkqlU4vU6K+l+iSSOOA/rIbXVbroBszza560YCtNtrs36Wbicba6Mv5zzTpwYtCohcjSat4dgEjNlca0bqu6yZxkuMZa4ZKy5vHtdVB+ivfBqRxAALyrG/s05Wv/bTQZn99KZNqrpU+l43/2Sdrzg7SOKCw0Rf7cPUlCU5nrF2XnK3GmvgWhnmOyNUna6RxAG7myy0u32mfTKpXk02TjiU9RjkOAnRtbb+5xY+5tzPWa0SueRoXX4YLIHy3ugq7fdKDtpK2+mK/6feLAnAtt6uacF9nbIu1VknjAPSV3SjSq74D9+iezlhqrQIAAETKcWqVNA4AACA0bokctVYBAACC47HZgTQOAAAgJM6JHGkcAABAWBwTOdI4AACA0FBrFcCOSTTUNrYCNcC9MkMlWt/u+mwIkXeJLtI4oMeGWlbVHwQQnFRLDbvuBLpFrVUAAIBIeSZypHEAAAChaDC12mSLA1sfAAAA2uaVyFWNrtWP0wEAAKBtHokcZboAAABC4pzIkcYBAACExTGRI40DAAAIjVsiN83+IY0DAAAIh8caOdI4AACAkDgncqRxAAAAYXFM5EjjAAAAQvPJ6V6z40+p7wD01EYzbbruBABHC605Y++dWyJ3hDQO6Cu70VPXfQDgyj533QN0z7tEF2kcAABAGDxH5OqTNdI4AACA22hUa9U/jWPrAwAAQNtarLVKGgcAAHBLzlOrFOkCAAAIC7VWAQAAIkWtVQAAgEhRaxUAACBS1FoFAACIFLVWAQAAIkWtVQA7JtFQW7vuuh8AXJihEq3ttut+oEveJbpI44AeG2qptOtOAHCUaqlh151At6i1CgAAEClqrQIAAESKWqsAAACRotYqAABApKi1CgAAEClqrQIAAESKWqsAAACRotYqAABApKi1CgAAECnHzQ6kcThkBpp03Qdcxd9I+taM7KrrjqBNZqJB133AVXwr6T+I8/WuUWsVTQyKdZPoob/ViBeGnhlr1HUXcDW/1L/vugvoErVWAQAAIkWtVQAAgEi1VGuVNO5OzexT111Amzhfe21lH7ruAtpkUha5oNHUapMtDmx9AAAAaJtXIlf1br1+nA4AAABt80jkKNMFAAAQEudEjjQOAAAgLI6JHGkcAABAaNwSuXxfDGkcAABAODzWyJHGAQAAhMQ5kSONAwAACItjIkcaBwAAEBq3RG52/Cn1HQAAALrXoLIDaRwAAEAIvBM50jgAAIAwfPK7e32yRhoHAABwG16JXPM0zpp+J3Qm0VAjSdJQG/vYdX8AAMA98EjkLknjun6YVzdRuvt40XVnAADAfWix1uodp3GyM2us0bMkad11bwAAwH2g1mqbtpKkVdfdAAAA94Faq20aSdpaRuQAAMBNUGu1NSbRUEysAgCAm6HWanuyXaskcgAA4EaotdqeoSRWyAEAgJuh1mp7RpJkSeQAAMCNUGu1JfkKOdI4AABwM54lukjjKh2skDMDDRibAwAA1+Y5IkcaVylfIWdGZqkXjZWadzPtulMAAKDPvEfkpCZbHPpea1XFiNxAUz1mY3FmrtQk9qnrjgEAgL7yGpGrSsbqx+n6Ll8ht9VEn3dTqjNJUzPqum8AAKCvPBI5ynR9IEvXtnqw2+JLdqONpEnXXQMAAH3lnMiRxn0oWyH3tE/jJEkbSeOuuwYAAPrKMZEjjauRXUNucfLVRJLMoOvOAQCAfnJL5PLdl6Rx5SqvIZeN05HIAQCAq/BYI0caV6m0yqoZ5h9QfRUAAFyFcyJHGveB8iqr2Ve3J+vmAAAAWuKYyJHGfWgolVRyKE/vAAAAWuKWyM2OP6W+w4mRyhK2bL8qiRwAALgSzxJdEmlchfMVcomkrRaNvhsAAEAt70SONM5RdiHgZ1bIAQCAa/FM5KqTtbtO486SNZNoLGlNpVUAAHA9jWqtVidr7psieuV8HdxUibZ67LpjAACgz1qstXq3aZz0rKPL/pqJptrqwXIFOQAAcEUt1lq92zROdqWFxmaUP+6J5qRxAADg+j653Y1aqzUetdGLeZY01FAzzdjkAITIJBrmtViGrGIFED+3RI5aqzXsVk/mWWMletaKJA4I1kTp7mPOVADRcxyRk0jj6tjN6YWTAYTGzjSTzFwTUQcZQA9QaxXA/cnG4qi7AiB61FoFcH9GkrZsSAIQP2qtArgzJtFQjMcB6AVqrQK4N9muVcbjAPQAtVYB3JuhJEbkAPRCS7VWSeMARIMVcgB6o8HUapMtDmx9ABCGfIUcaRyAXvC4jlz16Fr9OB0ABKJ0hZwZSVpzMW8AsfFI5CjTVRKTUf6i4GNlWZsDdOdshZxJNdZaW83Nk1103T0A8OGcyJHGlRoV5cu8kMgB3RlJUvF2yiRaaqXPklJJqUjkAETFMZEjjQPQB8fXkDNDvWhmn3dFuwAgMm6JXD7qRBoHIHIHK+RMohd9yXevDiQxWg4gOh5r5EjjPvDAujcgCrsVcibRUrPdRUieJW312HX3AMCPcyJHGgegB/IVcibRUs/2ufiyXbA6DkCMHK8jRxoHIH7FCjkz0que9mkcAMTKLZGbHX9KfQcAUcpWyGVbHFgOAaAHGlR2II0DEKlshdyz1poba16M/3UgASAo3okcaRyAaGVVVp/sg77STGMtzbzrLgHAJTwTuepkjTQOQNgOq6zarX3STNLEcP04ABHzSuTqkzX3TREAcGOnVVaz1b/j4lMzNEnXXQQAPx6JHGkcgKidVFm125Pb5/llgQEgGtRaBXAvjqqsSuZohM4MNLDrBt8VADpErdWrMmk2BmAfuu4JcO+Oq6we2OT/jrkkMID4UGv1uobi8gZAGE5XyBUpXJHITfSl6y4CgC+PNXKkcQAiliVyB6NudqOVlK2LM6nWTKwCiA+1VgHch0TS5iRZ+6KlUjPQQAnjcQBi5JjIkcYBiNyTNjqprmq3+mzGGmphWR8HIEpuiRy1VgFEzm70VPr1BZscAMSLWqsArsy3DJZJTHr1FgPvFiMz9WwxNmPPFhOqvwLw47xGrkAahz4wAw20PV4vZYZKJK2PLxObX7Riy0L4pkyqidnYmUeTqaZm7TXZmWpiVnbl0WKukVl4/VZTDczCblzvbhKlSszq7LLDH7fY6huPPt0NzligSku1VknjEBMz0VJjLc0yK8lkEpOad6UaKdW7mReFmszILPWqicZamjff0RVIkhloImnqXvzKDDSVlHq0GGkiyWOEzYw18mwx1VCJTwtNNVAin1G8VIm8RwrvAWdsP5iRWTLmfAW25lAqq3T3mZWVLbmX/fiW7N+6nxXfoTT/YFRx+7IqLnEfGuUfpN32o3Hv35RYaSqrpZWGelWq5OB3+pp/9Fr8ZjXQu6zGXff96rFp/feql/yDuXOLpW8v9Jp/MHW8f6K3/MOJc4v3/MORY4vB7sOhY4virHrXoPXfaxbT5S2fSy32njO2KjL5K1DX/fD4TVpZLdt/ht/34TUiVzW6Vj9OBwRkqpndKrs07MiM9KIn+7SbnllIGpqxSZXooZissxutJU267npszGhXkt5x9Vc+ViZJU+NU+dRM8xqq7uN+011NVddxv1TF/VzX+80P2rop7uc37ncPOGP7ZaQ3MzfDy7/Rx8zAjK7/U0LgkchRpgt9YAYa2uwiFNmL+Ytmh2ur8pU1U430dLT2Zi1Rp8PbYdrjkJ6Y4yTGIWkyh5OXTilQPnVbtHCY+synbjNOU58H6ag0Mg7pxEE6Ko2ZftrjjO2liV6vO81q0nw6/vX0zZpJzIt5cXubGAfnRI40Dj0x3lXbzE7kZ3t0bbH8HdxQX04Wqde8szNDVuScMqkO/1gOHXZ9To9ajBxiuh8rk9zG/Y7Tw6nDe/bj1G1S9xJgThPK2nE/c5pQeu7z7bWWz1iTmJFJTWpG7qswcRUjLc2ry9scf2auxH5jH/Ws4dnZNNFYY91fIkcah94Y7qptZi/5p7spixeLTcnXK3YsmsRMtWQa51i+zUHax7hm6vNgrKy44ltNCrQbK9vufkbNeNlurGzj3KIYK1vnlxOuH/cr0tFVnoLUj/sV6egib8GWh70Wz1iTmFRvmirRUEu992tcJkpDzc27Sdv9PZiJBvZx9+np28HifO6P2sWJ6f6Ts9sqvl52S/n9Yj/Y7NB1Txr0fZItklYiK6uXit/pyW80X7r+UvLd5lrK6k3RLiU/ejwt/l4Ptzlovv/4gxYH2xzctjwcbnNw2fJwsM1h7Lbl4XCbg9uWh8NtDm5bHg63OWi4/7jF32vEmx3aO2OV6FUvu00SiZayenfdjBLiEelmh/Ljpb2tKXrbbXp5OY+PrPLtMX05PNbIMRqHPrDP+QRMMS5zaiRJZ1cky+5d9h5uq5m+ojbAqd02h62eJD0pi/oHU58nY2XF++kPtjzsx8rsTPtRvI/G/XZjZXZht07jfsVY2bNd2f2430dTn8VtM7u2buN+xW3PdmPXeQu2PORaPGOn2tjdBKzd2gdtlGjJFGsQxnox7ya9fHuCmWibPR9MopFOngX5XyCfa04Gz/mCwKRx6Jnsz8XJ6ZxfTPT8JM9O/rN0rVitw5P9zD6d2Up2a2Z5YpLqc9ndD9aVPdmtZDdmlk9IzvVQ0eJoItauzLMmylKgx9IW+6nbR0myC7PSSNnUZ2n5roOp2+xnzMxYQ0kDk9ryFkU6WiR9s3w1zshM7HNpiyId3eTfcaaJEkljM/K6xHHfXX7GTrQ+ielMcyWayOdy1QGKZiL+b2puTzTV1Gy0kN/Fuo+Ndr/3sRKdPguq35THq3aIMj0fmPSbVN1/vevhx/YPpla77skFj+G17Dejcfnj0rs+HIxXqkgnrk4eR0u/192ZcRCzj6c+dy12UTycBi1tcTZdWzf1eT5dWzf1ed7nj6/2VjZdmz+nrN6LSb2TFmd91iT/4K2132vEU6vHv4uzrzqfsfmk9ftJ9KOOy+HSp94db0qbTbYWk/G75/3ReVf2tdiP+jscXRD4uFH1tyv/etcP9grhI5HrvC8NH4HfepvshfWjtVf9SuTaOw4ieZACnf0RPUioDtKjXQr0VtKi9Ltpmn9QknSXf7fzBLLuu+0SyJeSFqXf7aP1fuXfbdei3SPiZ+jlZ2yeMp/EIO649DqRy453vWjabB1j+cpmWfVshZz1r7V6SZkus+x6/LF137rdrXeP/KuuO3CxquH1j9bb5AP0ZkgVR0fPR9f7+mjqcz8Re7AGajf1OSiZ+kx3LQ4uOrGb+hya6XF119Op253Kqc/Tqdudp3zC5mzqs2TXbeZRb5Kk6Wl119Op24Of8Xq9X0uULj5j7dZ81lhlly5xrp+Lm0s01lgyW620Kdmb/JFsle7xOdrDFXLyH5ErPii5Z8UtnSertzhqRuR6e8Q7Ipe9kz15n1c+0ZJ/9XX/WcX3i/Z9/cHjaPObnY28VU19Vo+8lY/UfTTyVjXuVz3yVjX1Wb3TdvfTT6Y+q0feqn76buq2ukWbR8TP0LbP2KPv6lh2LbzjDkbksuNNc/9xufzsKvs70LPSbZ4jcnVjbh/dAgRkJGl7NrJW/q4/GzMp3smP72CH6v/Vjy18l8XJBVortjxUjpWpastD5ViZqsb9Trc5HLUo3fJQOVamqi0PZ9scDpWO+51tczhu0WZxoe/0dYvfrQtXOGNNoomkFVtKArbQWqsmcyDFNpjSy0T37Dfulcg1T+OsMVa9C56kb/W3Tvfr2yP/Sn/XdRcukZ/i538cyk/y7Kv7fVCR73Bz8H35jszLlU59HlwSpKTJPgUa724/uCRISYti6nNiFrvbS6dud8qmPkunbg9+RrZcYmLyyZ6P0lHJbs2TXrLva/IE9yAdLdlja7fle3WbMcu4i1Vd6YydK9FaX7p+dCix1aLiL4Kr8iR/JGmdn4EDyWuqNlgeidwlaVz+b4t/mMJgUpdKjf175GakuFf9+a23SbQ/9cclt8NHkQJNzbPdfjxWJp2kQCu7/XisLG9xMu5XUs3huMV+3C/N0qeTK9Sdtzgc98sSgX06WnqRkZJxv101B55Rta5wxpqpxlrroTRRj00srzDf6b/U3ufyFC5TkuSfrJCb7qq1xK52ljlfI1d84ez2iq+f3lJ+n9gPdq123ZOLfm9O623yHXBLKysN9VZxmYp+rZG76u/1eO2ZSwWH4/s4VXA4rPhQcxmT/Hd/cNkQlwoOx/dxqeBwfJ+PL2PSeswjv/zIFc7Yiey+ykOsR68qO1hZzVus7lCy/jFf3ZpfGkhvtzj3bnFQaxX3aCBpc7buIhtRKbvo77NGJjWpXvSFHasXO6jyUDdWljuo8lA3Vrb7GUWLpGbqVpJ0UuXh46nbosVhlYeDag6VLY6rPBxUc7heqHuj5TPWTDTX077KAzq30ZO+so8tjMR9JNvHupayHcu9OfdcM35Zaq1+EB1G5KI6NFBaejnXSVXVTQ2V7i8zWfFMiHa84+Bx3OT3ut9zWj9WdhBfK6ulS7VTq4Nxv/3O8Zp337t7Fv++143X7Mb9HFscjPsVLVq75K/TY4v2GdruGauJ3vfPOCXxVlvtyYjc/Br7hvMLQB9eyTLN/nrsPov09avksbo+US6t7lB+z9gPErmuexLKQSLn+XNejz6tjdxBCpQdc4cW70dfqH1MB1Of2TGtbXH6sjSpbTE++cKNLnwReyLXaiwmx9OtGpdd2jmOI/pE7r08QW/ppy11cDlgpXrNk7tBPunek4lV56lVJlUBtOj0grk17PHGhoptDictDqdeN/V7je3xfT6eus1arI6WSldsczhqsThafM02h5szE030+Wi6dcgFgTux0Uzf2KcrTm9+0bPGZmlSk5o3SQ92rUdJqUm1rNjBHiXHXaukcQDas9v1KcntD+pu12fWwmFt0+5SJ1LpJUFKFJc6kVSfKub3Gufr6VxbFJc6cUpH0S4z0VwzTY9enib1byTQso1m9W97LmW3ejQzjTTQRg/ZXxn7bFYaKdHn/qRxroncyTvT5kW6gD4yQ40kDQ+udIY6RQrkMFaWK1Igh7Gy3c/ILnXieDGDg0udPLuNlR1c6mTmtg3m4FInbHO4MTNVKpVcMorfwy1tNXM+gy9mN6cXGDn/Svycp1b3SOOAgkmNNVavGkpK9GKssb2rrHsVu6lPt7EyHUx9Oo9j7aY+nUdc8qlPj7EyO9Na5dUcqsy0UXk1B1yRGe52Cp9gJ/rNbDXTN7dL4+6FdyJHGgfs2Sdrzo5YLs7ZMTvT2vPCnzNtXMfKck/aeq6FeZTj1O3Bz/BIR3fr/ZjOuzG7LjlXjTUsA7odu7YeZwpctVRrlTQOgLdHef1Rt1vzKK/RE7s1j34F8uzGfPGbILcr7xYL88g2BwDt8EzkMv5bHPJaqwCw4z+l5Z/++K9avEmL3q3SAdAVr6nVqmSsfpwOAAAAbfNI5CjTBQAAEBLnRI40DgAAICyOiRxpHAAAQGjcErn8AoqkcQAAAOHwWCNHGgcAABAS50SONA4AACAsjteRI41rhmv8AwCA63EbkTupjEZ9BwAAgO5511oljQMAAAiDdyJHGgcAABAGz1qr9ckaaRwAAMBtNKq16p/GsfUBAACgbS3WWiWNAwAAuCXnqVWKdAEAAISFWqsAAACRotYqAABApKi12o6lsc5H2nVnAQBAP1BrFQAAIFKOiRxpHAAAQGjcdq1Sa7XK6kZtAAAAznhWdpBI4w7ZFWkZAADoCrVWAQAAIkWtVQAAgEhRaxUAACBS1FoFAACIFLVWAQAAIkWtVQAAgEhRaxUAACBS1FoFAACIFLVWAQAAIkWtVQAAgEi5JXLUWgUAAAiOd4ku0jgAAIAwUGsVAAAgUp6JHGkcAABAKBpMrTbZ4sDWBwAAgLZ5JXJVo2v143QAAABom0ciR5kuAACAkDgncqRxAAAAYXFM5EjjAAAAQuOWyE2zf0jjAAAAwuGxRo40DgAAICTOiRxpHAAAQFgcEznSOAAAgNC4JXKz40+p7wAAANC9BpUdSOMAAABC4J3IkcYBAACE4ZPf3euTNdI4AACA22hUa9U/jWPrAwAAQNtarLVKGgcAAHBLzlOrFOkCAAAIC7VWAQAAIkWtVQAAgEhRaxUAACBS1FoFAACIFLVWAQAAIkWtVQAAgEhRaxUAACBS1FoFAACIlGciRxoHAAAQigZTq022OLD1AQAAoG1eiVzV6Fr9OB0AAADa5pHIUaYLAAAgJM6JHGkcAABAWBwTOdI4AACA0LglctPsH9I4AACAcHiskSONAwAACIlzIkcaBwAAEBbHRI40DgAAIDRuidzs+FPqOwAAAHTvk38T0jjsJGbkcK+f6Kf5R/+qH7ru8l0h8jjkdr5K3+0++kH/2nWn74pv5JOuO4wQeCdypHE4MNGk6y4AcDTUsusuAGibZyJXn6yRxgEAANxGo1qr/mkcWx8AAADa5jEiRxqHnbUePrj1O/1n/Zvik0Qj/VwDSVv9SSutD+/5e33PGpwW+UT+f3zwfTZdPxC07OmD1VQ/0d/rl/tPhxrp55Kktf6klbb7m/6P/htrLVv0E/0n/bv9p0QeDdmaQ6msUru749ntFV8/vaX8Phz9O5TuPxnaP9hTf7bTw7u/Kum6x305jiP/x7PI/0DkOc4ODfRWfJLY39g/nz1vfm+Hhw0mXfe4L4eG9ZEfEHkOh6P+Dqns/gWiqnn1Nz78vOsHy3H9Q6P9J1Nb5Y+HLw0vXfe5HweR52hyaFl8OLI/VDxrTt58DW/fyz4eevWM/LsGXfSTI/yj/g4H7/KrGld/2+OvdP1gOa59KCneY5aNCB07+AM17rrf8R9EnqPJsf/7/puaZ80f9qNDr133ug9HEfmEyHNcfNTfoYXRuOJrXT9YjmsfxbNlaF38rmjGO80bR/43+8gzvXrHh4bFh39weNb8eZ9QpF33PPaDyHO0eVBrFa0xI02zj37ndP//qHH2QaJ5132Pm2/k/1H5dWGJ/H3Lf/tT/cLhzl/vn11TM+y665Ej8mgRtVbRnt0fp3/r2OB3xVa6keHCwpdoHvmxGXfdeXTDpBpK0kD/4NjiF8X7BfEG4BJEHu2i1ipaYkYaSNJQqXObg3eabqWDUMKMs8iPPCL/7f4lhETuXuXn3O/0tXOTf8ieatKQkaELXBb5QdfdR2i8LgicIY1DqfyP02+8Gv298tcD0onm8hD6Rf4fi2ak0HfJJNkTYOw0uVf4ev8s43nTkBlkkZ80jTx/K3HCO5EjjUOFkSQlXn+cds0kOZbzxrk88q7TqkfN3Aupo1/ydODnns12TxZG5JrKQ/hds2ZEHmc8EznSOJQr3t/75wS7FxLSiUaIPBppOCD7NSO5lyLyaFmDqdUmWxzY+tB7+d8Y33eZ0i+KZfe8z2wmj5vvyIr2Y6dE/h4Nsv99692QkdwLjaULI88ZiyNeiVzV6Fr9OB3uRZO/MPxVakOTFdC8Et+7Js8a/7cMOJJIF0Y+adAYPeaRyNVNnnIpEgAAgFtyTuRI4+Bi26DNputO9wKRx23wrOnKuusOIFCOiRxpHGrkecSfvBv+WLww8PrQDJFHY6sGbf65+IC8opnV7n+educ4kccRt0Quv6g0aRyq2HWWUPj/edq14I9TI0QejeS//u+bNlzbJoPA2J1v/8u7IZFHOY81cqRx+NBCktb60bPZ7l1mk7eokPLIrfUXz2a7kZVF1w8AHcjPN9+R3H8phoA5X5vKz7d/9mxG5FHFOZEjjUON/H3mf/dq9GPxV21tmeBrKv/D/l+9Gv3I+/u7VozkPnu+AfifxQeM4zZURH7hGfnd+U0ihxOOiRxpHGrlf15m+hePRr/iXebliDyaWEnSVr/yaPJ9UXh7axnHbW4lSRv92qPJ98Vb3q3ljMUJt0Rudvwp9R1wzm6K58mvnN9p/rbIIranzzG4O4y8KyIPPRUjQ67r5P6yf4bxrLnEUzEaSuTRhtoMy6SaamafzpuQxqHKVKnDvX7UzxpdNAPViDwA9En9/Gb9iNxf5//lSONQb6bf1t7nR/2SZKJ1RB4A7oytOZTKKt19ln9wdi9bd0tcLeLqbUgRUVp8OLI/2Gq/sUnR5D2Ex9eX50gW+T+7Rf4uIsJz5KMWGuzPP509V/Zf/8EO942mfY7I7Vsk9vclMS9+G388jHwAveU33n1Ezo9GtVarx9xcN0WE3SKu3oYVkeKag9JKP6sYHfpRD/r1fkwo6f7x9eU5UkT+pxVrb04ifwcR4TnycQu7Oazb+eVodWtxR6vf6mf7Taqr/ex9HyNy6xbSVr/cRX5/t+wOT4eRD6C3/MZvHxEXLdZaDTNsPPW6iEhmq1/rs54OLnz5o/5JX/Sz/UbJbf136kdEbpLGfbXb2KZf1kX+LiLCc8S9hSQt9FM96p/0o/YpxVSfD5P/mUb3EpEbtHgqAptF/jCN+15P+ny2syGuxxdXb8OMiBvnzQ7tFemKq0VcvQ0lImas+eE7fWkkaXNaC2qhcXiPL+7nyGnkEw3LIn9HEQmjRdi91Tea7xO0Shs9ankfEblVCzM4j/xIVdcE6rq3l7WIq7ehRMTYsttP3bzWalwt4uptOBGxC31zXC9gpdVxMrHVF9K49lucRn57HnndV0RCaBF6b+3GPuxHhyo86zNpXNstyiK/Io27ixb+36ma24icwgtCmIHuewu/72TGGhfJ2pG1VprpPbzH17fnyJm1hvcdEZ4j1V83A0001uD0O2irlRZ2cX8RuVULM9Cbymy14u1uP1s4z7E4jch5JHIhBSGkQN9PiwsSiuIt5kZrLew2zMfXt+eIHjQqErcs8tdPnsOOCM8R9xb5ObvWwq6JyHVb7GL+qNFuYcRKK7sOsbf8xm8ZEbdE7lPdHQohBSGsQN9Li7h6S0Ss0clSGyJy+xZx9ZaIdB6R59B7y2+8i4jUuVmt1bhaxNVbIkJEQmgRV2+JCBEJoUVcveXxhZjG3bzWahwt4uotESEiIbSIq7dEhIiE0CKu3vL4wkzjvK4jF0IQwgx031vE1VsiQkRCaBFXb4kIEeHxhRgRN96JXFxh46lHRIgIESEi3beIq7dEhMcXUkTqeCZycYWNpx4RISJEhIh03yKu3hIRHl9IEanXYGq1ybBgLC3i6i0RISIhtIirt0SEiITQIq7e8vi6i4gLr0SuKj+szzTjahFXb4kIEQmhRVy9JSJEJIQWcfWWxxdCRMp4JHIhDlaGPSDa1xZx9ZaIEJEQWsTVWyJCRHh8IUTEjXMiF2IQwgj0vbWIq7dEhIiE0CKu3hIRIsLjCzEiVRwTuRCDEGag+94irt4SESISQou4ektEiAiPL8SIVHNL5Ki2GkFviQgRCaFFXL0lIkQkhBZx9ZbHF1Ya57VGLqQghBnovreIq7dEhIiE0CKu3hIRIsLjCzEiH3NO5EIKQpiB7nuLuHpLRIhICC3i6i0RISI8vhAjUscxkQspCGEGuu8t4uotESEiIbSIq7dEhIjw+EKMSD23RG4WThDCDHTfW8TVWyJCREJoEVdviQgR4fGFGBEXDSo7xBU2nnpEhIgQESLSfYu4ektEeHwhRMSNdyIXV9h46hERIkJEiEj3LeLqLRHh8YUUkTqf3O7m/uN8H1KYLeLqLREhIiG0iKu3RISIhNAirt7y+LqISL1GtVbjChtPPSJCRIgIEem+RVy9JSI8vhAi4qLFWqthho2nHhEhIkSEiHTfIq7eEhEeXwgRceM8terbwb60iKu3RISIhNAirt4SESISQou4esvju01E3BI6aq0G+7PDbBFXb4kIEQmhRVy9JSJEhMcXYkSqUWs10J8dZou4ektEiEgILeLqLREhIjy+ECPyEWqtBvmzw2wRV2+JCBEJoUVcvSUiRITHF2JEPkat1QB/dpgt4uotESEiIbSIq7dEhIjw+EKMSJ2iXfUdUk11VqSrmGw9+3r1LXG1iKu3RISIhNAirt4SESISQou4esvju3VEpi5pnWsiBwAAgJuqT+TqLz+yaqkvf63/13U4AAAA+qR2RA4AAABh8irRBQAAgHCQyAEAAESKRA4AACBSJHIAAACRIpEDAACIFIkcAABApEjkAAAAIkUiBwAAECkSOQAAgEiRyAEAAESKRA4AACBSJHIAAACRIpEDAACIFIkcAABApEjkAAAAIvX/AbXW+zSkOzzWAAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDE3LTEyLTIyVDEyOjAwOjA4KzAwOjAwk1PzxgAAACV0RVh0ZGF0ZTptb2RpZnkAMjAxNy0xMi0yMlQxMjowMDowOCswMDowMOIOS3oAAAAUdEVYdHBkZjpWZXJzaW9uAFBERi0xLjUgBVwLOQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%%tikz\n",
    "\\usetikzlibrary{calc,patterns,decorations.pathmorphing,decorations.markings}\n",
    "\n",
    "\\tikzstyle{spring}=[thick,decorate,decoration={zigzag,pre length=0.3cm,post length=0.3cm,segment length=6}]\n",
    "\\tikzstyle{damper}=[thick,decoration={markings,  \n",
    "  mark connection node=dmp,\n",
    "  mark=at position 0.5 with \n",
    "  {\n",
    "    \\node (dmp) [thick,inner sep=0pt,transform shape,rotate=-90,minimum width=15pt,minimum height=3pt,draw=none] {};\n",
    "    \\draw [thick] ($(dmp.north east)+(2pt,0)$) -- (dmp.south east) -- (dmp.south west) -- ($(dmp.north west)+(2pt,0)$);\n",
    "    \\draw [thick] ($(dmp.north)+(0,-5pt)$) -- ($(dmp.north)+(0,5pt)$);\n",
    "  }\n",
    "}, decorate]\n",
    "\\tikzstyle{ground}=[fill,pattern=north east lines,draw=none,minimum width=0.75cm,minimum height=0.3cm,inner sep=0pt,outer sep=0pt]\n",
    "\n",
    "\\node [draw, outer sep=0pt, thick] (M) [minimum width=2cm, minimum height=1.5cm] {$m_1$};\n",
    "\\node [draw, outer sep=0pt, thick] (M2) [minimum width=2cm, minimum height=1.5cm, xshift = 4cm] {$m_2$};\n",
    "\\draw [thick, fill=white] (M2.south west) ++(0.2cm,-0.125cm) circle (0.125cm)  (M2.south east) ++(-0.2cm,-0.125cm) circle (0.125cm);\n",
    "\n",
    "\\node (ground) [ground,anchor=north,yshift=-0.2cm,minimum width=10cm,xshift=2.03cm] at (M.south) {};\n",
    "\\draw (ground.north west) -- (ground.north east) -- (ground.south east) -- (ground.south west);\n",
    "\n",
    "\\node (fill) [ground,xshift=-0.15cm,minimum height = 0.3cm, minimum width = 0.3cm] at (ground.west) {};\n",
    "\\draw (fill.north west) -- (fill.south west) -- (fill.south east);\n",
    "\n",
    "\\draw [spring] (M.east) -- (M2.west) node (k) [midway,above] {$k$};\n",
    "\\draw [thick, fill=white] (M.south west) ++(0.2cm,-0.125cm) circle (0.125cm)  (M.south east) ++(-0.2cm,-0.125cm) circle (0.125cm);\n",
    "\n",
    "\\node (wall) [ground, rotate=-90, minimum width=3cm,anchor=south east] at (fill.north west) {};\n",
    "\\draw (wall.north east) -- (wall.north west) -- (wall.south west) -- (wall.south east);\n",
    "\n",
    "\\draw [damper] (wall.15) -- ($(M.north west)!(wall.15)!(M.south west)$) node [midway,yshift=0.5cm] {$b$};\n",
    "\n",
    "\\draw [-latex,ultra thick] (M2.east) -- +(1cm,0cm) node [right] (u) {$u$};\n",
    "\n",
    "\\draw [-latex,ultra thick] (M.north east) ++(0cm, 0.5cm) -- +(1cm,0cm) node [right] (y1) {$y_1$};\n",
    "\\draw [dashed] (M.north east) -- +(0cm,0.8cm);\n",
    "\n",
    "\\draw [-latex,ultra thick] (M2.north east) ++(0cm, 0.5cm) -- +(1cm,0cm) node [right] (y2) {$y_2$};\n",
    "\\draw [dashed] (M2.north east) -- +(0cm,0.8cm);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Additional Mass Damper Diagram\n",
    "\n",
    "via https://tex.stackexchange.com/a/13952"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABPkAAAJzCAQAAADXKw4OAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAAAmJLR0QA/4ePzL8AAAAJcEhZcwAAASwAAAEsAHOI6VIAAAAHdElNRQfhDBYMAAimj/aTAAAyfUlEQVR42u3dv24rXbgf5ncdfI2LBOBxkQAucsCNBGlzuAt3TkFdgnbOFUhtUkmXQF2CBLi3pUuQmtO52HKdIiJswCncbMWwCwMJMCn4R6RIajjkombN8HkG2NqiZsjFRYnz43pnZqUqAADot79puwEAAJyayAcA0HsiHwBA74l8AAC9J/IBAPSeyAcA0HsiHwBA74l8AAC9J/IBAPSeyAcA0HsiHwBA74l8AAC9J/IBAPSeyAcA0HsiHwBA74l8AAC9J/IBAPSeyAcA0HsiX8elQRq03QYAoHQiX9eNYtx2EwCA0ol8XTeOUdtNAABK91fbDeBIAh8AUCtVbbeAo6Q/MahS260AAMqmsNtpaRyDiORoPgDgSyJft41X/gUA2EHk67bZkXwiHwDwJcfydVgaxJ/5f39U07ZbAwCUyyhfl10u/2ecDwD4gsjXZeMt/wMA2KCw22HpT3xMtva31Xvb7QEASmWUr7PSZazOrnt58B0BAL0n8nXX+IvvAABWKOx21MrZugtKuwDADkb5uupyj1sAACJC5OuuzUKuyAcA7KCw20lpGG9bbnZBZgBgK6N83XTV4FYA4OwZ5euk9BbDLTdPqx9ttwwAKJFRvg5K462BL2KYHM8HAGwh8nXR1QE/AQDOmMJu5+w4dWPBKRwAwAajfN3z9UjeTdvNAwDKY5SvY9Ig3tbm1v3sPX6YhQMAWGeUr2tuvgx8EQPjfADAZ0b5OqV2jC/COB8AsMEoX7dc1QY+43wAwAajfJ2y4xLMnxnnAwDWGOXrkDTZK/AZ5wMAPjHK1xl7Hce3YJwPAFhhlK876s7VXWWcDwBY8VfbDaCBu7Xv/kX885Xv/k38Y9vNAwBKpbDbWWmyNpJ3V9223SIAoFQKuwAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAQD0nsgHANB7Ih8AQO+JfAAAvSfyAWSVqlSlatftu38y//8kPadR288B6B+RDyCjbZFu9fYq1WwxinEM2n4WQP+IfADZ7Ap29YFv8ycAOYl8AJkIfEC5RD6ALAQ+oGQiH0BGAh9QJpEPIJsMgW/c9nMA+knkA8jk+MC363xfgGOJfABZCHxAyUQ+gBMQ+ICyiHwA2R0R+F7abjvQTyIfQGaHBz5n7wKn8lfbDQDol/rCrcAHfD+jfAAZCXxAmUQ+gGwEPqBUCrsAmZhyDSiXUT6ALAQ+oGQiH0BGAh9QJpEPIBtz7AKlEvkAMjHlGlAukQ8gC4EPKJnIB3ACAh9QFpEPIDtz7AKlEfkAMjPHLlAekQ8gK4EPKJHIB3ACza/DJ/ABpyTyAWS06ySM+rE/gFMS+QCyMekaUKq/2m4AQF8IfHyf9Nx2C77da3XbdhO6TeQDyELg41uZmo+GFHYBMhL4gDIZ5QPIJkPgM3ZDE9OYtt2Eb+CvIguRDyATs+zy7Z7O4fg2fxd5KOwCZCHwASUT+QBOQOADyiLyAWR3ROB7abvtQD+JfJ1V3VZpZTmDozmgK8yyC5TH6RsAWdUXbgU+4PsZ5QPISOADyiTyAWQj8AGlUtgFyMSUa0C5jPIBZCHwASUT+QAyEviAMol8ANmYYxcolcgHkIkp14ByiXwAWQh8QMlEPoATEPiAsoh8ANmZYxcojcgHkJk5doHyiHwAWQl8QIlEPoATaH4dPoEPOCWRDyCjXSdh1I/9AZySyAeQjUnXgFKJfACZCHxAuUQ+gCwEPqBkIh9ARgIfUCaRDyCbDIFv3PZzAPrpr7oV0jCGbTcSOmxaTdtuAt/FLLtAuWojX1zFTduNhA67i9u2m8D3EPiAkinsApyAwAeUReQDyO6IwPfSdtuBfqov7C5ctN1U6JjnthtAW46ZZTf5vQFOYu/IV/nkCY248sa5qi/c7n/UH0AuCrsAGZUY+NI4DRpucdl8i8atar7F1am3SIOmrUqD5LI6dITIB5BNkYFvEI8xabjFfbNrNaRhPKZmjzGOx2aBLF3GfbN4lW7iPo0addZN3KdmFyabxGPTeAztEPkAMqk/Uq+Vku5NDOKqUViaxCBuGoWl+4i4ahSWJhEx2T8spUFM5o+z/xY388fZd4th3MSg0RbjuIqBS5nRDSIfQBZlzrGbRvNAsneQSeO4arjFZYwjmoSldBOjiEZh6SaGETFsMJY4iUFEjBuMJc4C5WWDeDxrTbN4DC0R+QAyKivwxUdwG6V949Viiz3DUvqIenuGpfQR9fYMS2m43GLPscRlcN17LHEeXCP2HkucB9fVPoOCiXwA2ZQ2x266Wrm/m32iz0qM2Tcs3axMy7lfWJqNvy3+v4+P+913LHGyssUeYTet3u9eY4lp9X6bjCVCS0Q+gExKm3ItrcejPcJSWo9He4SllfG3iL3C0sr4W8ReYWll/C1ir7HEteC631jizdp88vuMJa4G10bHJUI7RD6ALEoLfDE7cSMi4iXeIyL2OIljstxifg+1YWkx/vY0/1ofliaftqgJSyvBdbFFzVjiSnBdPkbNFsNPW9TG42VwfZ/3lZM4KJ7IB3AC7Qe+5YkbEddxN//fl0FmJcb82nOLxfjbNK7jISJqw9Jy/O21+rVnWFqMv71Uv+I1IurHEhfB9an6FdOIqB9LXITIu7iex+O6scRFCx6WWziJg9JVNUtMZv+pW89isawv8y+TttthafHVr/b/yfI/z1HFOEsbnuf/mVRRRfyef3fzxRYr68Qg3ubfXe1c/2OdyypiEH/m342/2GJlnRgufzDaucXHOsMqYjz//58Y7txibZ2V7wY7t7hcXSeu5t+9fdFTN6vrLPaT8dzC79hZvMO00b99XOpXEPksloOW+ZezeEO2bH3tDwh886CWIfJ9BJdZ1KmPPssY83v+/WXtFp+iznoQ2rrF/fw/99vvYcsWa8F15R4ed27x+9MWj+vfb6y/EW4/P+aWLT6F2/p4fLLfsrN4h6n7LbE06scvVhD5LJaDlvmXs3hDtmx55Q8KfLki30oouVze9ilu7dxi+ehfR5/18bf5bb+/3GIjdtaNJS5j59vKFl+OJW7Gzhgut9g6lrgZO2O0bOfWscTN2LnPWOKJfs/O4h3m8ytkOWypPXgkTWZHWZjuG5qZ/3HdVbdtt4TvNXvlm1yHb3WL9Bzj+Lfx58hGDOZHzL1UF8tHGcTb/Bi3i+rl8wbpfn4c30N1vbxtGG/z//6sXje2eJ4fx7fyO57G8RwREe/xs5pubPF73qrb6m5522U8zrf4Ub1/Wn8Qv+fH8f2qnpa33syPo5tWPzYeYetzXOzHVnvj6+e43OKp+rWxxdbnmB7j8nNvrGxxilM7Nnq/v+Z/O1tePxqpzdZG+SyWg5b5l7P4DG7ZeN2rfW///JPlyFqeZW2U6nPpduUnO0apdhdeN8ff5rfvLLzuevTdY4m7Hn33WOKuR989lrj90b8aS9yxxRdjictHP8VyFu8w238PLAf24xcriHwWy0HL/MtZvCFbPr3qBwe+zJFvMxLtOIlj5+07Cq/rJ258+snWsLQ7RO06iWNb4Xj+kx0ncew+uWNnpN15vOKukzh2n9yxM6BOojrhchbvMNt71nJgP36xgshnsRy0zL+cxRuyZe01PyLwZY18b5vHlW2PPrtH/3aFoq9Ou9h+EsdXRxJuv7evjiTcPpr31ZGE207i+PpIwm2P//WRhNvubWX07zTLWbzDbPv9sBzcj1+sIPJZLAct8y9n8YZsWXnFjwp8+c7Y3dnKjehVe0LERvTZPf42//lG9Pr69IZt0WtX4Xh3m78+X3hb4fXr84W3ncSx3GLr+cLbnmXd+cJH/76dxTvMV6+TpXE/frGCyGexHLTMv5zFG7Ll4xXf8Xuw9/m7J498mxcY+fJM3m0Br/YiJhsF1rqrAn4eS9xdOF5u8Sng1V8VcOOSMrVXBfwc8OrO5N0MeKc8k/fr16Bfy/orZzmyH79YQeSzWA5a5l/O4g3ZsnjFd/wWNL1gywkj3+cy7l6XKl4LS1+Pv823WCu8flU4Xm6xfuHoPS5uvD6WWH+9vs+F17rguhkj99ji01jiKa/XV9eWPi31vw2WBv34xQoin8Vy0DL/chZvyJbZK77jd6D5FfpOGvk+zbKxz6wcK2Nu9eNv8y2WYWmfWTnWx9zqCsfzLVbGEveZlWM93tZfaLqK9ZM19pmVYz2snnZWjvmXs3iHOWU/ntNSv4LIZ7EctMy/nMUbsuWL34BDJl07deT7iD6LUPK7ZouPkb09Y8xH4bWucLzcYnnP9aNp8y2WI3tfXwJ6ZYtF4fV+3/G3j7bsE1yrWBnZm3x9CehMv19n8Q6z+M1oux1dX+pXEPksloOW+ZezeEO27Hz9Dwl8J498K2Fp70fcOJN4WLvF77Vva49mWxk/nC1ve2yxfjbsW22bPp8/WxsilsfvLZbH2i3Gn2440XvAae+9rGXfV8vy9fI3bV8KGqB/6mba2PzJciqkl/gOt7E608VDVf+o12vf3W3OrbHlMda3eP969er90xa3e2xxt3bDddSopvHwRRu3bfG69hjve2zxEk9rW9zVbQHfQ+QDyOzwwPddU1uuhaU9YkxENV3ZYrpPjKleVuLVa7XPFk8rgfeletpji7v4mAjuaY/gGtVtfITVu81p5La4W4nHD3tE3fVAXRtc4buIfABZ7Y51pQS+iLWwVDv+tlhvGZb2jTEf0WffeWCvt/yv7jFm9gqua/e85/jbyujjdL/5bFfGEl+qh322gO8g8gFkVFu4LSLwRcQiLO01/haxEn32Gn+bbzG7730Kx7MtFuOH+xSOZ1ssxhL3G3+LlcLr3uNv1cN89HHfGPoxlrhvDIVvIPIBZNOhwLcISw1CybzwunfwmY8l7j/+FjEbS9yrcLx0G+/7jr+tbNFs/O029iwcL13HvoVj+CZ/td0AgL5oXrhtL/BFRMRtvDeKMRHXcbnvaNryMYZNjmar3tNtROMtGrWpmqa7ZqfJVK/pLhqVaKuXdOfEDcoi8gFk0bnAt3mObP0WzcbfIqJhpIyIfcvGK1s0Pl5u32L2yhaNS7TNt4DTUtgFyKg7gQ84LyIfQDYZAt+47ecA9JPIB5DJ8YHv4zQPgLxEPoAsBD6gZCIfwAkIfEBZRD6A7IqfYxc4OyIfQGblz7ELnB+RDyArgQ8okcgHcALNr8Mn8AGnJPIBZLTrJIz6sT+AUxL5ALLp3qRrwLkQ+QAyEfiAcol8AFkIfEDJRD6AjAQ+oEwiH0A2GQLfuO3nAPSTyAeQiVl2gXKJfABZCHxAyUQ+gBMQ+ICyiHwA2R0R+F7abjvQTyIfQGZm2QXK81fbDQDol/rCrcAHfD+jfAAZCXxAmUQ+gGwEPqBUCrsAmZhyDSiXUT6ALAQ+oGQiH0BGAh9QJpEPIBtz7AKlEvkAMjHlGlAukQ8gC4EPKJnIB3ACAh9QFpEPIDtz7AKlEfkAMjPHLlAel2IGyMoMHECJjPIBZCTwAWUS+QCyEfiAUinsAmRiyjWgXEb5ALIQ+ICSiXwAGQl8QJlEPoBszLELlErkA8jElGtAuUQ+gCwEPqBkIh/ACQh8QFlEPoDszLELlEbkA8jMHLtAeUQ+gKwEPqBEIh/ACTS/Dp/AB5ySyAeQ0a6TMOrH/gBOSeQDyMaka0CpRD6ATAQ+oFwiH0AWAh9QMpEPICOBDyiTyAeQTYbAN277OQD9JPIBZGKWXaBcIh9AFgIfUDKRD+AEBD6gLCIfQHZHBL6XttsO9JPIB5CZWXaB8vzVdgMA+qW+cCvwAd/PKB9ARgIfUCaRDyAbgQ8olcIuQCamXAPKZZQPIAuBDyiZyAeQkcAHlEnkA8jGHLtAqUQ+gExMuQaUS+QDyELgA0om8gGcgMAHlEXkA8jOHLtAaUQ+gMzMsQuUR+QDyErgA0ok8gGcQPPr8Al8wCmJfAAZ7ToJo37sD+CUzLF7kDSOiNfqPSIiDWIUEVF9Oug6jWIQEdNq+sX9jOIyIiKm8bJ7vTSeX5z1PZ6qaRrEZfXQdI3Gz3AYw3ivXrc8o/nzXt46e/6f1oXzZNI1oFRG+RpJ96lKVariOZ5jFJGG6T5+x1VcxWOq0n0azNe7TG8xicuYxFt6S5db72u2TsRrRFzFW7pPwy1r3aQ/cRnTeImXiLhP9/Ecl83WOOB5XsVzXMZzep49ozRIk/QnJjGOSfxZeZ7j9By/4you43nX84TzIfABzW3b+59EVbPEZPafuvXOY4mrmMT9/JtxjOItLuc/GcTvqOJ3FVXEfTzGYLlFFdVirbV7+h2jT/38Z/WW+a2fbov7qOK5yRoHPMtxvMWgiriZ3VOM4ndMls9osnyek/gd4/mtw/iz7Xme8zL/Mmm7HZZvfsWrfW/f/pN4jmrxd2Wx7F7mX87iHWb+5ag9W8lLVHG/2Md+Qz9+sYLIt6tPLuNtLbRdRhVV3MQk7tfW/h1V/Pl0D8Oo4vd6QIrnqOJPDFduGcWWP+f4/fFrX7/GQc/vOa6qqCLGUUUV43hb3QHFaP7sJ+u/ovEcPf6DPKgfZ1/O4g3Z8vGK7/g92DPwiXyWfZf5l7N4h5l/6e0eJqqo4s/H0MqpFoXdw93E3erxa9VTRERcxWV1vbbeS0QM0vpU6ZcRMYrHNFq57SEiBnH1aa3NI/yeGq3RWBrGaH4k4Gyo+THuVo9TnD/nmxjH7dpRfa9hQnjOXoaSrr8iCpCGaby2h4qINErjNF4c2LO8dZDGm+vS2CBu4i1NPvdvTiLf4QZbT5EYxt32tetunUfG0dp9rX8/89pojeYul1f/n0W+h/XnOf/DHsWv9dM4wh88Z88su/SD47lbMot996c6tk/kO9yusbTt0yWtB6KHeI2Il2pz3dVP+NOIuNrI/NOVR65fo7nRMjLO2vI5wi6C4HTL7dOAsyXw0Q9pHDfxs7qOuxjHY0QaxXNE/KguqtvqZ9zFVTxHRKRJTOKu+lH9qq7jZwziUejLYBBX8ZYe0wnG+0W+7Ko9Yk/1Xv2sUnVRs9psdG2W+W8WL381XRl1q1+juZfZvc4vvvL0aSxvEV4/hco0jGEcO74IPSLw0Vk3cVe9x/xwnTSOx7itPg7keYqIUbpMkxjExWLgoprOrz5BHrNx06u8ZV7X5WtZGsQoxrG1LFpN03Xcx/IIvxTxFC+rca5+jeaWW88C5GaAHUdsXoVwvrbIBxFxVOB7cTQfbUrDGM0HJBbHc9+uH8+dIiJuIuLi0/HcY7+7WQ3jPu7TQzxVL8ffWYRRvlalq/QYf+I+IqZxt23Mr3qIX2uh6zLuF8dW7LvGwWYh9PMFpmdjf7sK0kedOAJ9YZZdOszx3GW5iuf0liY5TpAxyteSNIr7GMV7XH89Jlc9xVO6jFGMYjQ/2WMcz/GzyRoHajaaN5uPxLF8sEfhVuCjYL04njtNehVBh3ETN2kaT/F0zExXIl8r0iieYxDvcbHPi1c9zUbP0igu4yoGMUpX60Gxfo3GLZwfybfxg+1jf1cxCGN8EAIfnfcy3590+3juUQ+LzB/B7+WwUq/CbgvSIB5jEBG/vg586T7drH5fvVa3cRHvEYsJ1erXONhXo3m7xv7mbwGuz8T5EvjouuphHvIcz12mYdzEc6rSY7ppurcV+dpwOfs09FVKT1VEDDcHpqvXuI2PYybq1zhUgyP50iAuY1nWTYP4/W09CUWpP1JP4KMzHM9dtsuYxO9Upef9g5/I14Zdf0ibLrfc+hSrn7vq14g0SFeNR97GEfG+MQq5/XPc7LT8RSH50h8956n5qRkCHwVzPHfZ3uMpbuOi2usAsRnH8rVh+/m0l/G+5SeX8fmYvM/HS9SskQbxFoOIdFttnxdki/nnuM1fo61h9dORHZex9+NA/wh89IHjuQv2Gi+HHc1nlK8Nsxdq7Wi7NIir5SGzkUbLUbrN+fYuY/0Pq26Nq3mQvGlw6ZZmR/INIuJ1dvRHutzyczgb5tilJxzPXaKXuI0f1c/q9rD9rFG+RtIoBjGch7VxuoppvFev809Di1MqJvESEa/V+/a1I6qHNIqrGKb76nq+zTAe4yFe4ioibuJ2Ef8i4j2e06+PwfJ0FTdxt/Zi162x+MkgRrHvL8n2M7J2HcPxFON4n/fPJH61/SpBW0y5Rm8cdzx36yPXtbNbFWWPv/vXeIqnYwvnIl8zk5VP4KO4j4iXuIjZ/IMLN7OrksfLjrUjorpOr3ETV2kcTxExiHHcVk8R6TomcZNu4q66jYiI17iIcTyn13lwG8cgflWrUax2jeop3cVNTOO1QeQbRsR04/iA2dWYNobuZxE2TSLisu4sZOgvgY8ecTx3KabxtOUqiAcR+RrZ/rmhetn+iearTxnVQzyk0fyM26fFaF/1kJ5iFNP5i3sXr9V7PMVTmk1jM43rT3+A9WtERHUbtxHpvsETvY3pxhGCUb2m65UJ2VZ/cp0e4jKm8XPjCk5wpgQ+usvx3IV4iJcqY4AW+VpUvcbr509D1fvHH1O1+r+t43P1a6wY7v/Jq5rG7dbbd17euXp1PSb4YI5dOs3x3G17j4dcY3sfnL5xPkZCGXwPc+zScU2P546V47lv6+6cGtO4rv62us1/yRujfGciXcZU0RW+g8BH5zmeuy0vn07QzErkOxdXcdScu0Azza/DJ/BRDMdzt+GkcS9C5DsTaRzDSuSDb7DrJIz6sT8oheO5v900rk9/DKTIdwbSIB6jU9cogq4y6RrQ2PX3DMo4feMc3Metoyvg9AQ+oLnvqsIZ5TsDlfkw4BsIfEDJjPIBZCTwAWUS+QCyyRD4XIgZOAmRDyATs+wC5RL5ALIQ+ICSiXwAJyDwAWUR+QCyOyLwmZIeOAmRDyAzs+wC5XFdPoCs6gu3Ah/w/YzyAWQk8AFlEvkAshH4gFIp7AJkYso1oFxG+QCyEPiAkol8ABkJfECZRD6AbMyxC5RK5APIxJRrQLlEPoAsBD6gZCIfwAkIfEBZRD6A7MyxC5RG5APIzBy7QHlEPoCsBD6gRCIfwAk0vw6fwAecksgHkNGukzDqx/4ATknkA8jGpGtAqUQ+gEwEPqBcIh9AFgIfUDKRDyAjgQ8ok8gHkE2GwDdu+zkA/STyAWRill2gXCIfQBYCH1AykQ/gBAQ+oCwiH0B2RwS+l7bbDvSTyAeQmVl2gfL81XYDAPqlvnAr8AHfzygfQEYCH1AmkQ8gG4EPKJXCLkAmplwDymWUDyALgQ8omcgHkJHAB5RJ5APIxhy7QKlEPoBMTLkGlEvkA8hC4ANKJvIBnIDAB5RF5APIzhy7QGlEPoDMzLELlEfkA8hK4ANKJPIBnEDz6/AJfMApiXwAGe06CaN+7A/glEQ+gGxMugaUSuQDyETgA8ol8gFkIfABJRP5ADIS+IAyiXwA2WQIfOO2nwPQTyIfQCZm2QXK9VfbDQDoB4Evp/TcYOWn6qHt9kL5RD6AExD4jtSkwP3admOhCxR2AbI7IvC9tN12oJ+M8gFkdswsu40Kmv03jeleawG1RD6ArOoLt/sf9Xf2nqrbtpsAfaGwC5CRwAeUSeQDyEbgA0qlsAuQiSnXgHIZ5QPIQuADSibyAWQk8AFlEvkAsjHHLlAqkQ8gE1OuAeVy+gZAFgLfd0hXcRkRUV203RLoGpEP4AQEvhMZKn3DYRR2AbIzxy5QGpEPILNj5thtu+1AXynsAmRlBg6gREb5ADIS+IAyiXwA2Qh8QKkUdgEyMeUaUC6jfABZCHxAyUQ+gIwEPqBMIh9ANubYBUol8gFkYso1oFwiH0AWAh9QMpEP4AQEPqAsIh9AdubYBUoj8gFkZo5doDwiH0BWAh9QIpEP4ASaX4dP4ANOSeQDyGjXSRj1Y38Ap2SOXYBsMk26dpf+tP1MSpee224BdI3IB5BJtll2/77tZ9IBZimBhhR2AbLIFvgATkDkA8hI4APKpLALkE2GwDeOiIvq7C/IXHdSi6DstB+aMsoHkIlZdoFyiXwAWQh8QMlEPoATEPiAsoh8ANkdEfjO/ig+4DREPoDMzLILlMcZuwBZ1RduBT7g+xnlA8hI4APKJPIBZCPwAaVS2AXIxJRrQLmM8gFkIfABJRP5ADIS+IAyiXwA2WSaYxcgO5EPIBNTrgHl2vv0DW9EAF8R+ICSGeUDOAGBDyiLyAeQnTl2gdLUF3bbefv5F/HPIyLiX8W/b+Xxu+J/i7+LiH8X/7rthhTuJiIi/k38YwuPbQd+ho6ZYzc9t916oJ9qI1/10sYuK03mke9fVnaYX0ij+LuI+L+q27ZbUrY0i3z/qJ/4DscEvrbbDvSXwi7ACTS/Dp/At4/qtkpV0lfQnMgHkNGukzDqx/4ATknkA8jGpGtAqUQ+gEwEPqBcIh9AFgIfUDKRDyAjgQ8ok8gHkE2GwDdu+zkA/STyAWRill2gXPWzbwCwB4HvBEZpssdaLy7aD/VEPoATEPiyGO9Z6Bb5oJbCLkB2RwQ+4QU4CZEPIDOz7ALlUdgFyKq+cCvw1dMXkJtRPoCMBD6gTCIfQDYCH1AqhV2ATEy5BpTLKB9AFgIfUDKRDyAjgQ8ok8gHkI05doFSiXwAmZhyDSiXyAeQhcAHlEzkAzgBgQ8oi8gHkJ05doHSiHwAmZljFyiPSzEDZCXw8Y1GabLHWv9d/M/z//3X+D/jP7Xc5v8h/m7+v/8n/l3rrTkjIh/ACTS/Dp/AxwHGjS/r87+23WTaorALkNGukzDqx/4ATknkA8jGpGtAqUQ+gEwEPr5PlXYv8SveV1Z9jbu4qFKV4m/jVzzEdO1nP7+6pxxL3K615iXuZo8ZP+L6U2te4scX93TRdp93nWP5ALIQ+ChDmsTN8puXuK6Woap6j6d4ikjjuI9hRESM4ne6qE54aaB0H1fLb57iulrGv2oaDxER6TLuYxAREeP4nS6q17Z7sK+M8gFkJPDRrjReBr73uK0uqunmOtVL/Iy75bf3aXCy1lwuA997/Kp+Ve+b61RP8SOe5t8M4r6dfjsHIh9ANhkCX9PzL2FF+ohML/Gzutu1XvVe3cbFvKg6PFXMWmnNU/yonnatV71Xv5bF6P0uO8MBRD6ATMyyS+tu5gXb1+3je6uql7iYx6zLdHmS1iwKti/bx/fWWvMUvxbPIY2+rb/OisgHkIXAR9tWirrX+6xfTZfl3RMUd9NVzILk+56teflozWn76VyJfAAnIPDRgkVUutv3FIjqbj6r8yAyl1PTxz1e1403Lt3NS82jdLPnFjQg8gFkd0TgO+G5k/RbGi2LurcNNrueF3dzH0U6nhd1n3Yfw/dZ9TEeeJpC85kT+QAyM8surViEtrsmGy0ulRLDzEfQLe6tWWte5ufujk53FvH5EvkAstod6wQ+TmoRspqOFC+KwHnH+Wb39t74KnunaQ1RfuR7VM8nm5s08bmRU6s/Uk/g40RmxdDXunNjNywiYsZRvjSc39veRd1TtoaZ0iPfICbpd5L1yeMm3tLV8XcDuwh8tGW5p2x8NGj1Pt8m5/Fzw/nXxjNpVK+LUziydg9RbOSrblfm5BvFczrhtcE5AxfLWRwHcZ+eXfOJ06gv3Ap8nJl9z9U9dhv2UGjki6ju4sf8kNKIiKt4U+LlUNVL9SPulh8ixvFbiZf8zLELlKzYyBdRvVfXcbEcFFbi5SjVbfxcOapEiZcTEfhgxfCbtmEPBUe+iIjqpfqpxEse1bT6pcTLaZljl7ZUx5z2MPuty3lNyMU7bePWpOE88invZld45ItQ4iUnJV5OyZRrtGpWxxg3fU9bzq+bMfJV03lka/4R5gStYaYDkS8iIlZPOR+EXTTHeF37fRr5fSIPgY+WHXpFu0Ov5/e1WQBtfoHn07SG6ETkS5fxOz7G9V7iZ6OpZGBFGqbneFweKTKNX9XF3rM/QgMCH99uEZIaHaechvP1m180+WuLe2tUl0ujg68uSK3CI9+nHfR7XFcXmX8pORtpkCbxtvL59y5+7j/3IzRhjl2+3/KKduM0abDZ/WIu3MzNeZnXUy5Tk+v93c+/em8+gYIj35Yd9I/q4Yg75KxtjBb/qG59iuQ0zLFLS67nX2/2Laemm8XEaJG5fla9L2fX3fu0yzSZl3WnYW9/AsVGvm3lXDtoDqOcy3cS+GhL9fIRs/ZZPw2X+9nr/HvY6m4+Zj3YszWjU7aGQiOfci75KOfSjubX4RP4yOBuMV1ZeqwbW0ujeFwUdU/0rni9LO7WXhshjeN58RwqhzecRHGRTzmXnJRz+W67TsKoH/uD41Xvy+LuZbx9dRRdmsTveRn1Y5vcrZkuRx1v4ovJFNIgTeJ5Hj9fnaJ5KoVFPuVc8lHO5fuZdI22rRR3B/G4fawvjdLqvvaEZdRlcTdiGM/bx/rSeGXPf7L4SVGRTzmXfJRzaYPARwmq25VZqy7jT3pOy9M50jhN0u/l+F7ENC5O+95YXSwj6Gyqy8dFa9IgXaZJeovn5Z7/Nez3T6iQikIaxM3atXvu4s7oXr30HOOIeKku2m5JWdJlTFZmaXyJa6N7nF6ewJeeYxwXjmU6HyeaO/6/j/8j/r52rX8V/zr+yzc8yf8x/vf4n/Zozb88yaNP7QFm/mq7ARFbdtC3Uj6HSsO4Xxndm8at0T2+kxE+Gno+/i4O9g/xD20//W9pjb+xiCgg8n3aQb/HrZM1OJTRYtq0PbrtDnQ71jdmf642R3bHO34y3nuLfxL/LP6b+Kef1vp/40/8if8Y/1+Wx9j/9n8S/yz+Nv7bjXv8D/Gf4/8+yWMfMsdvj7Uc+dIw3la+nYbD6znGc6xeflRxDOiIHFdxrJ3ZZXZM3dNHHS3/Y9TfngYxjlEsJ2I73WMXc+xaMVrvjzSKiTLcoRzLt+7TKN97PDjZHyhbqr4n8LWxRZmtPV+tn7FbvVYXy4s1RgzjMT2n4VF3ydmq3qvb+Lkc2h/ETXo70aHRAJmUFYgEvv5qPfJFRFQP8WPlJO5xvNVfpxu22/gQ8ZwefYgASlVWIBL4+qyIyLcxOjO7TvflEXfIWfv0IeIyfqdJ220C2Ee3IlSJjy3w7VJI5ItQ4iUnJV6ge1LVPFqVuUWJraWgyBehxEtOSrxAl3ycT9k0WpW5RVmtpbjIp8RLXkq8QDeUHd8Evj4oLvJFKPGSkxIvUL6PENM8WpW5RUmtZabIyBehxEtOSrxAyco82aFvJ3NQbORT4iUvJV6gTGUEIoGv/wqOfBFKvOSkxAuUp4xAJPCdg8IjX8T2Em/bbaKrlHiBEnUjvgl83fZX2w3Y0+pRfO8H30uB5qNM7x/TXH/jY49m/Vq9HHtPHbN+TOh7v36jgO7pRnzrauBzMsdCByJfuomblZ30U9xW07bblNFzRES8xEULjz2JWeA8oz+HNI5JjJbfvsbt2QVeoDBlBSKBr78Kj3yfdtDTuLaD5lBpEJO4Wn77HnfV3RF3B5BBWYFI4OuzgiPfxg76obptu01016fR4oe4rRR0geJ0K0KV+NgC3y7FRr6el3P5Vsq5QDd0K0KV+NgC325FRj7lXPJRzgW6oi6sNI837WxRcmvPWXGRTzmXnJRzga5ofrGRsrcoq7UUF/mUc5c9sTrOudv/Mvs3Pe9zn1UbZwW3SDkX6I4yi6F9KfYSUVTkU85dM4r954X4pw3WPRPKuUCXdCu+CXzdVEjkU84lJ+VcoEu6Fd+6FvhSJQTOFBH5lHPJRzkX6Ja2A9GptyijtRQQ+ZRztzu34+7yUM4FuqaMQCTwnYOWI18axuPK+N57/FSA4wjPaye9/PLxAeiGbsQ3ga/b/qbdh6+m8TOelt8O4i3dtNwjdFj1M+7i40PDc5qkwRF3B/AtuhHfuhr4HMm30HLki6im1a+4iMWxe4OYpN/J+accqLpd+xBxE2/p6oi7Azi5sgKRwNdfrUe+iIjqpfqxMjoziud0b3SGw2x8iLhPz2mfaxwCtKCsQCTw9VkRkS9iY3TmSomXw336EDGO30q8QDd0K0KV+NgC3y7FRD4lXvJS4gW6JlXNo1WZW5TYWgqKfBFKvOSkxAt0ycf5pU2jVZlblNVaCj2DOQ1jEpfLb3t9dbX5C/Aabcw2Mr8iYt//MNIkrlYuBXQXdy4FBJSm7Pgm8PVBkZEvIiKN4z6Gy297O4dCCS9A//80Nj5E3FYPbbcJ4EOZx7715dg+ZkpIHLsbtz4608uZUkt4Ac7jz+PTh4iXuK1e224TQETbgajvgc8cuwslJI6vmtf7Em8JL8D5/DEo8QJwrkpIHHVN7HWJ98hzku/i7yPi38ZRF7TpU3/WUeIF4Dx1IPJFbJZ4q+u2W1SG9BzjiHipLtpuSZco8QJwfgq7SMt2aRjDlcD3HnbQHCwNYrwS+CJel5dxAYDe+qvtBtRz/BX5pMuYrI3wXVcCHwBnoPDIpwRHPmkY9/Fx7OQ0bqunI+4OADqk4MJuGqbHeF4Gvve4ri4EPg6TBmkSbyuB7y5+CnwAnI9iR/mUc8lHOReAc1dk5FPOJR/lXAAosLCrnEs+yrkAMFPYKJ9yLvko5wLAQkGRTzmXfJRzAWBVIYVd5VzyUc4FgM+KGOVTziUf5VwA2NR65FPOJR/lXADYruXIl4bxuDZ7rhEZjvEYo5XvBD4AmGv5WL5qGj/ibvntIN7SJA2OuEPOWvUzbuPjoIDH9JiGR9wdAPRG66dvVO/VbfyMl+UNN/E7XbbdKrqquosf8bD89jJ+p0nbbQKA9rUe+SIiqtfqIq6XozPDeEzPRmc4TPVeXcdFLI4IHcRNekvjo+4SADqviMgXEVE9rJV4x0q8HK56WSvxDuNZiReA81ZM5NtR4jU6w4GUeAHgQ0GRL2JLidfoDAdT4gWAhcIiX8RGidfoDEdQ4gWAiCIj30aJ1+gMR1HiBYAiI1+EEi85KfECcO6KjXwRSrzkpMQLwDkrOvIp8ZKXEi8A56rwyBexLPEuDOM5PbbdJrpqXuJdzOQ8+xAxOuouAaADOhD5ItIoVqdgm66M00BDaRCX8VHQfY+H6vWIuwOATvir7QbUSYO4iZvlt+/xUN223Sa6K93ETXzM6vIQt9X7EXcHAB1ReORLVzFZ2UE/xW01PeLuOGtpHJP4KOK+xm31csTdAUCHFBz50igm8XGqxjSu7aA5VBrEJK6W377HXXV3xN0BQMcUGvmUc8lJOReAc1dk5FPOJR/lXAAoMPIp55KPci4AzBQV+ZRzyUk5FwAWCop8yrnko5wLAKsKiXzKueSjnAsAnxUQ+ZRzyUk5FwA2tR75lHPJRzkXALZreY7dNFwbkZnGncDHEW5WAl/Eg8AHADMtR75qWv2Iu1gU3obxO03S4Ki75IxVF/ErPj403KfnNDri7gCgN1qOfBER1W38jKfltzfxlq6OuDvOWvUUP+PjZI2xDxEAEFFE5IuoptWvuFiOzgyMznC46r26jR/xUdL1IQIAyoh8ERHVy1qJ1+gMR6imayVeHyIAOHvFRL6IrSXeYdttoqu2lHiN9QFwtoqKfFtKvCIfB9so8fptAuBsFRb5IuYl3qfj7weWJV4AOHMFRr6IiHB1PrKpfIAA4OyVGvkAAMhG5AMA6D2RDwCg90Q+AIDeE/kAAHpP5AMA6D2RDwCg90Q+AIDeE/kAAHpP5AMA6D2RDwCg90Q+AIDeE/kAAHpP5AMA6D2RDwCg90Q+AIDeE/kAAHpP5AMA6D2RDwCg90Q+AIDeE/kAAHpP5AMA6D2RDwCg90Q+AIDeE/kAAHpP5AMA6D2RDwCg90Q+AIDeE/kAAHpP5AMA6D2RDwCg90Q+AIDeS1UJd/G9XqqL73/QNImbtp94Q3fV7fc/aHqOcdtPvJkqtd0CAKj3V4b7eIlY2U2/fPrpuOb279/iNcNzbm66bEN5PbJ9i2kr/bT66pTWI1/fDgAFyzJE93Enn0c8Fj/ZdXs7W7SpzB7ZtUUJPVVWj5TYUwBQL8OxfGXvmr/6yfcru0dKijFl9kiJPQUA+8hR2I2IbTvB5rvN79yiDWX3yOYW7QWbUntk1xYAULpMZ+x2Zdcs8DXdoqSeKqNHBD4AuihL5OvKrrnEGFNGjwh8h27RufPVAThTGSJfN3bNZYzSdKNHBD6BD4C+OcGlmEvcNQt83SpgdqNHBD4AuiN75Ctx11xmjCmzR0rsqTJ7ROADoEsyR74Sd81lxpgye6TEniqzR0roKQDYX9ahihJ3zWXunMvskRJ7qsweKbOnAGC3jKN8Je6ay9w5l9kjJfZUmT1SZk8BwFeyRb4Sd81l7pzL7JESe6rMHimzpwDga5kiX91usHkY+44t2lRmj5TVUyX3SFk9BQB1skS+5hc/bneLNpXZIyWefVpmj5TYUwBQL0Pk61qhssQLDZfRIyWF4zJ7pMSeAoB9ZDuWryu7ZoGv6RYl9VQZPSLwAdA9mSJfV3bNJcaYMnpE4Dt8CwAoX5bI15Vdc4kxpoweEfgO3cKxfQB0Q4bI141dcxmjNN3oEYFP4AOgb7LPsVvmrlng61YBsxs9IvAB0B3ZI1+Ju+YyY0yZPVJiT5XZIwIfAF2SOfKVuGsuM8aU2SMl9lSZPVJCTwHA/rIOVZS4ay5z51xmj5TYU2X2SJk9BQC7ZRzlK3HXXObOucweKbGnyuyRMnsKAL7yV567+RgsbH4kVBtbtKnMHikxxpTZIyX2FADUyTrHbpm75rJ2zmX2SIkxpsweKbGnAKBexjl2y9w1l3WmZdk9UlKMKbNHSuwpANhHpsJud45LKznwldaH7QWbUntk1xYAUDpz7H6bsnukC2OhZfSIwAdAF5lj95uU3SMC36FbuCAzAN1gjt1v1Y0eEfgEPgD6xhy736gbPVJiT5XZIwIfAN1hjt1v040eKbGnyuwRgQ+ALjHH7jfpRo+U2FNl9kgJPQUA+8s6VFHirrnMnXOZPVJiT5XZI2X2FADspjoFANB7Jzh9AwCAsoh8AAC9J/IBAPSeyAcA0HsiHwBA74l8AAC9J/IBAPSeyAcA0HsiHwBA74l8AAC9J/IBAPSeyAcA0HsiHwBA74l8AAC9J/IBAPTe/w9vnJJV9MUJigAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAxNy0xMi0yMlQxMjowMDowOCswMDowMJNT88YAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMTctMTItMjJUMTI6MDA6MDgrMDA6MDDiDkt6AAAAFHRFWHRwZGY6VmVyc2lvbgBQREYtMS41IAVcCzkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%%tikz --no-wrap\n",
    "\\usetikzlibrary{calc,patterns,decorations.pathmorphing,decorations.markings}\n",
    "\\begin{tikzpicture}[every node/.style={draw,outer sep=0pt,thick}]\n",
    "\\tikzstyle{spring}=[thick,decorate,decoration={zigzag,pre length=0.3cm,post length=0.3cm,segment length=6}]\n",
    "\\tikzstyle{damper}=[thick,decoration={markings,  \n",
    "  mark connection node=dmp,\n",
    "  mark=at position 0.5 with \n",
    "  {\n",
    "    \\node (dmp) [thick,inner sep=0pt,transform shape,rotate=-90,minimum width=15pt,minimum height=3pt,draw=none] {};\n",
    "    \\draw [thick] ($(dmp.north east)+(2pt,0)$) -- (dmp.south east) -- (dmp.south west) -- ($(dmp.north west)+(2pt,0)$);\n",
    "    \\draw [thick] ($(dmp.north)+(0,-5pt)$) -- ($(dmp.north)+(0,5pt)$);\n",
    "  }\n",
    "}, decorate]\n",
    "\\tikzstyle{ground}=[fill,pattern=north east lines,draw=none,minimum width=0.75cm,minimum height=0.3cm]\n",
    "\n",
    "\n",
    "\\node (M) [minimum width=3.5cm,minimum height=2cm] {mass, $m$};\n",
    "\n",
    "\\node (ground1) at (M.south) [ground,yshift=-1.5cm,xshift=-1.25cm,anchor=north] {};\n",
    "\\draw (ground1.north west) -- (ground1.north east);\n",
    "\\draw [spring] (ground1.north) -- ($(M.south east)!(ground1.north)!(M.south west)$);\n",
    "\n",
    "\\node (ground2) at (M.south) [ground,yshift=-1.5cm,anchor=north] {};\n",
    "\\draw (ground2.north west) -- (ground2.north east);\n",
    "\\draw [damper] (ground2.north) -- ($(M.south east)!(ground2.north)!(M.south west)$);\n",
    "\n",
    "\\node (ground3) at (M.south) [ground,yshift=-1.5cm,xshift=1.25cm,anchor=north] {};\n",
    "\\draw (ground3.north west) -- (ground3.north east);\n",
    "\\draw [spring] (ground3.north) -- ($(M.south east)!(ground3.north)!(M.south west)$);\n",
    "\n",
    "\\draw [-latex,ultra thick] (M.north) ++(0,0.2cm) -- +(0,1cm);\n",
    "\n",
    "\\begin{scope}[xshift=7cm]\n",
    "\\node (M) [minimum width=1cm, minimum height=2.5cm] {$m$};\n",
    "\n",
    "\\node (ground) [ground,anchor=north,yshift=-0.25cm,minimum width=1.5cm] at (M.south) {};\n",
    "\\draw (ground.north east) -- (ground.north west);\n",
    "\\draw [thick] (M.south west) ++ (0.2cm,-0.125cm) circle (0.125cm)  (M.south east) ++ (-0.2cm,-0.125cm) circle (0.125cm);\n",
    "\n",
    "\\node (wall) [ground, rotate=-90, minimum width=3cm,yshift=-3cm] {};\n",
    "\\draw (wall.north east) -- (wall.north west);\n",
    "\n",
    "\\draw [spring] (wall.170) -- ($(M.north west)!(wall.170)!(M.south west)$);\n",
    "\\draw [damper] (wall.10) -- ($(M.north west)!(wall.10)!(M.south west)$);\n",
    "\n",
    "\\draw [-latex,ultra thick] (M.east) ++ (0.2cm,0) -- +(1cm,0);\n",
    "\\end{scope}\n",
    "\\end{tikzpicture}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Structural Analysis\n",
    "\n",
    "This could be useful? https://github.com/hackl/TikZ-StructuralAnalysis -> https://ctan.org/pkg/stanli\n",
    "\n",
    "Following example from: https://tex.stackexchange.com/a/394933/151162"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuMAAAHECAQAAAA5New8AAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAAAmJLR0QA/4ePzL8AAAAJcEhZcwAAASwAAAEsAHOI6VIAAAAHdElNRQfhDBYMASZjQsodAAAucElEQVR42u2dz48zXXbXvzWTmTBRIKmMQESEUeSXGcGGTT1IgQ0bt0DKJhs/f4Jb7GDVz5/Qz4awde/YthdsMtnY0rACgtpSlqDosSIRaSabNrwBJmEyc1m47G67ft1bdatu3fLnczTz9uO2q25dV399fO655yRGAAAQL98IPQAAAOgCMg4AEDXIOABA1CDjAABRg4wDAEQNMg4AEDXIOABA1CDjAABRg4wDAEQNMg4AEDXIOABA1CDjAABRg4wDAEQNMg4AEDXIOABA1CDjAABRg4wDAETNL4UeAEB7ksf8h63Zhh4LxE6y1EyStDdPocfiBjIOMfNw/gkZh64sNJckbRWZjBNUAQCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhBxgEAogYZBwCIGmQcACBqkHEAgKhpKePJKjHJqvYZWfKavCbZ1aOz5DG3WehLBwCYAr/U7mXmPpGWicx95TN2yZ022iR3Zvfu4Zke8p+22oe+eACA+GkdVDH3etKyziM3O91J2lx75AAA4I8OsXGEHAAgPI0ynqTVv0PIAQBC0+yNbxByAIDx0izjGUIOADBemmX8HiEHABgvjTJunnwKub4f+oIBAKaFRd64eUqklTbJnTlUPsc2j/xz6AsGAJgWVtt/PAr5j0JfMADAtLDMG/cWWnkQAAB4xHozvh+PXH8S+oIBAKaFQ00VT0IOAAAecSqNhZADAIwNxwqHCDkAwLhwLlSLkAMAjIkW9cYRcgCA8dCqbQRCDgAwFtp2/0HIAQBGQUsZR8gBAMZBaxlHyAEAxkCjjCfZRUvkCxByiI9kbv3UvaHtN9Qyjrup2Ru/7m1/AUIOcZE8a+Hy9NDjhTEzlrvJJqiCkMN0uNeTpFQrpZIOulfxvs20UCbp4HpwuDHGcjeZBlOmV70qq33OUkYvSmufs5LRSvPzA/OmM2NYk51/eHR+5Wv967SS0Sb09WFDmjb5D87ve+i7qbn7j0UDNvsytvrXPX4iAViSZDreq5XfMvUp9BghFsLfTRb1xr0K+e/2ezkAVpwWprZVTzAHbWv+LAHeCH43WbWN8CjkP+zvUgCsOd7JO0P0G7oT/G6yzBvPO2l2XuzUv8Mfhy4kS81KHv7nyeO7f23NtuEwR/+p6VlNv4fISeZ6nzD4D07/vbib9uap4TDB7yb77j++hBygC7PSRoC/o985/3xoatxtEcuE2+BQejf99sWjd/WHGMPdZNmLU/IXWgFoj/mkpk0UT41OREUsM3m4ehYyP3HMTk2etu03u6B3k4OMI+QwCup97UZfXBWxzCTV48WzvuI75Q3QdLdEcTc5yThCDuExT7X+eLMvXhXLXFw+wkb8W8Dsa/3xZl98FHeTo4wj5DACqj0kC1+8Mpb5QBDlJvnc8neSxnI3Ocs4Qg6hMU+VG5vtffEL/ymZJc+akZlyi9T44/a+eOC7qV33H7JWICxPpRkGNpHMUyxTej5XKkqPNS8s/mhhiqy1rHi8mVHcTW27/yDkEJLPpTJu44u/xTLf7txMB6V13lOSaWZs/qghQsw22akYW2jOF5cc76YkVab58fleVbFDIZkWRbMojYX5Ma1KHswsXpflPywKjz5UvCLVg14pkzVl07LkQYtyay53k1I96lUbrfISXM+aeRt/pxc7Czkyjvmxd3fSyV6sXveQ/3BVj1MlHwJa5n9yX0S1w4lbXqPwvVmIrP3dpFQvej7rYKqNTJN22lvjEmeyrPHkWeyEQJhtIe3Q5itwdZbvoTQ8eNBn/YZVjBTi5vruWVulCNrfTQ/am4+n55mDudNeqS9dbM5UWSHkMEoKf3hWryrP8p2XxTLNk/lktqzs3ATrhn+XY383LZVeNXz7LClVjbra0yzjO4QcRsnlH5rV8mZllu/WKscFJovZXdwTB5sFbfu7KcmUaq7niwePZ7Dv5VlDs4zf+RRy/aqPQQNIZn/h89ild1XUvzAHw9afW2dd8XM19nfTXgddyb3Pb3nN3X8OXoUcnwf88fbHZuU9aQSVoWG0uMu49d1kDvqgT/r4/rFcLb1s0rfp/uNTyH/gY9AAkt7/sdkuQtpVhoYbxOzP3vLecuuOw91k9ubzldwvJNnfubXYdf/xJ+S/72PQAFLeGuuI1Z/SGCpDw4g5SeoAd1NyXNy02e5vgWVNFV9Crj/0MWiAHCcZb+6ZCDfNkHfTSql2l2GW9liXxvIm5AD+OPpPa8tYN5FxqMHsjpHq/tdZkgcttJO3MiUONVXMIbnTRqtE1bUGbGqtAPjC7JO9ZvVfa5NjqaJUWR6NnCWPuQe1p6Y4XLDVsj5a7eNuSpZ61Fr3Ht0Jxy2rqV5ktKx9Ts0WfTbjY35NKxmLchAnu9xyvbI8x6PYjH8TpoUqa+vkz+h+Ny2bztFi3M4v6CTkyDjm17TUa+/nQMZvxJTKqsBahzMs9fpWSEupn7O5d/8hRg5jYsuCJfjCHLTrMxycLPWgu3ex93lF3XxH2nT/QchhNFzt5QToRo8bFJOllvpw8TGR+dn+0677D4udMBqsivsDWNFfc5BkqZU+6yF5/+BS9z6O3bb7D0ION0Jy7NeSJQv6/0Bbkgc9SiUhlKE245dDaAWmTvKYmMToJU8xe05MYpJN6FFBfCSZHst/48fFTUyXF6faKNN93dfaJNNGUu6RJ3Od/gzuaGALANCd5u4/q+rf4ZEDAISmOaiyRMgBAMZLs4w/IeQAAOOluW3EvU8h1/dDXzAAwLSwSDg094m0TGQqMxwd0g/p/gMA4BWrvHGPQv6j0BcMADAtbNtG+AqteKkgAAAAJ6x3cfrxyPUnoS8YAGBauLSN8CPkAADgEaeaKgg5AMDYcCyN5VHIN4k88QtJbzUFfq739QVM/q+fvXvs5/qr/Kf/p1/UHvmn+mnFb77WX/saPgCMhp35FHoI7jhXOByhR365TPvN0ud8Z5CRAAAMTotCtSMUcgCAm6Vd2wiEHABgJLRtG4GQAwCMgpYyjpADAIyD1jKOkMOZr/NMoL+QJP0s/291ng/AWImy4WSjjCepOVT9rpOQ13T/STKl+rv6TUnf0d/Td/Ubkr6jv6lflfQt/Q1J0rclfVPfrshMgSH5W/l/v1v6250Okg5518F9/t9d9X0FAC40e+Ob5K4nIa/EvT9dkub9EjP9mn5dP5D0d/Rbvc8e2FBRZT6RpIN2Okn88WfR3A/AjWYZz0IIuSvmoOMf/1Xv8nfyfrwW6Xv6QV+jAGdSza8fyreFHd/PXf7/B+G/A1TQ2FI5WWqlne7q/oSSlZZ6qhbyt+bL2o+jpXKJvKeiN1EMHIMy+zxQsxPyDjdPo4x7FvLf17/JHwgq4zWjzDTTTLMyPxFGzFan4MxOB0IzcEtYyLhnIT8xShkvjBhRr+f/6ldCD6GGfS7r0lHm92YfekgA/rGScY9C/l/OkekIZLwwekS9nh/rr/Q/9Zf6ib6nV0kzzUIPqYSdDrnfvtcevx3ix1LGvQn57+k/5D9GJ+NXV4KoN7PVQXsd9LX+Qj+RzvN0/E42JpF/k/atiLZDZFjLuB8hT+bjWOL0TTLTTJnS/H9Qxl577SRt33vAyVHMT5I+nsXmvfbvYu2EY2DEOMi4DyGfqoxfXeVR1GfKRiFI4+SgnXbH/y+/n5Kj7378/5lSSZnSwCM+Rtr32uvgvrsBoB+cZLy7kN+GjF9dcaaZ5iMKIIyPo6RLWxuvNzmK+cmLT6WAQa2t3oQdjx0C4SjjXYX89mT8fOWpMs010zyoR+mPr/XtvCyCb7ba62An6W/kwZks99pDRd5PoZiDdvjrMBTOMt5NyG9Xxt/NwTHYkgUOEfjjT/UT/Zqkf9TDsXd5RH3fThKTTKlOW73eZH44DsfIuq5WBAB80kLGuwg5Mv6eJMsj6FPJdPmP+on+Wv9b380Xe/1ylPS99t3unHwH78l7H9pvP4VhyIgBb7SS8fZCjoyXkxwDLlNZEj3Gunf6Zf1Ec6mH/J2jj77Tzk88Opf2uY7R9iFzZXbnpdMDsg7taCnjbYUcGW8imeeCPo0l0ZOg78w+mWnWS6b9VjvttfMdiU5OH0DHcMxQ35iO3jqxdXCgtYy3E3Jk3JZ8SXQ8EfSfyejbHV6/1057bU8eZzLLl3tnHq9wlwt6T3dWkil9F4wZ4p3Z54umOx34e4EqOsh4GyFHxt1JMi00H1G45X/oVyoaRNhxDIds33zNJM0/rvx9aF19aPRFcvxucRL1/t+j7TkPhvRGONNcqDar+2rnKuTIeFuSVIs8ZXEc/Fh/rp/qN/XbrY9w0PZSzvPrPBY5yDwtPb4L6/Q/Ke9i7Kf/9slpXyy++o3TLOOvuvMn5Mh4V5JMc81HlNvyx/pv+kv9Vgc/ulTO86ud663MQXd6iqLXkcz1tnTabxjmlNzYOZcHYsNGxuVPyJFxPySp5sq06MU3b1N89qCt/lj/S3+/Q2CkRs4lKcnOPnp3Se85il5Ncqq8k+bLvX1x0I6o+q1gEVTRRt6EXGtk3CfJTAtlWvRy8D/Td5xj4DtttdMu397U7jvDIY+d19wfyaliTXf/djtcyKX0St6yd3wu9l5zEvUt+S9TxKb7j08h/6F+N/8nMu6R5Bho6WOJ7c/1Y31D/9jxVeujX91xg9O2Sc7PGS/dUzQP2mrfdLb+SbI8HbPPJdO3WpOI+iSw6/7jU8hPIOPeSWZ53NzOp/u5vulw8D/Sq77l6C/utdVWW3PouGPVQs7P+S7dN1HttB06gl57TSdvvS9f/SjqZL9EjG33H19C/gd44/3jlKL4M33L4dB7/VDf0j9xFMqdtlqf7p3k+L2hNznPzzHrXOYgD+6MaW/lecG0P18933pUXUAYxod99x8vQs4S53A4pSh+rV/o160PvdMf6E/1Tx3THw+5b557fLmct/Ext9q9fSw0zIKPKPpxQXQ7Bv/84tr62hf7xpaIegy4dP/xIOTI+PAkmeZaWPluf6af6vvWB15rq/+qf+G8xLrTWuu3r+9JlnvnrjJ70Fo7rW19Ri9R9G2+/DpKPzWZ56GXfhZL88AL6Yzjw637T2chR8ZDkWRaWCYo/mf9TN+z3NZz0Fpbs26xxHol5h3kfKe1m6ecR9GzDoJ+Ksw12nu4516xTo0+oG9cu/90FHJkPCzJXAstLITyoB/pp/qHlsJ80FpPZpfM8mx2eyEuiHlrOX+3nOowH7NOqZGSQ7w+JO9E3UfWfXEODsdl0rHPw1Rx7/7TSciR8TGQLLSwCoQc9O/1f/S3LYV5r7WezD73/F2ksUTMW8v5VtvisSzm5NSbqa1/HoWc59faZwPwt7wXlkgHo033nw5CjoyPhSTVQgsrqd1rrf+kf2YZNDlLcpLm6Y+2wlgq5nkapWs1mb3W2pl1q3npVlkyIjnPr/hN1P3vKm3VkA9cadf9p7WQI+PjIplpYbn8udNaf6R/aRlfX2t7Wnx03GlaIeYt5fy4CNpKQpLsHEFvw3ZsyYrWV516LyB8mhEkvSfadv9pKeTI+BixXv48aK3PmlnG16W11m8ecTLX3LoKTKWYtypAsNe6i4ect/Jo1wx7l/vnkcm5dFFt0nfw5RRNb9ljFS5p3/2nlZAj4+PFevlzq89m6xBfX2v99l47ZMzUi7lryKaQt95ihk4Ji20WRCOW8/zqs1zQZ54zXzq2zYZu3X9aCDkyPnYs5Xmvz+YpSbXQ3ErM91q/37DjS8xbFe692FPacpbmeX6Lu38evZxLF40+/MbTkfRWdOv+4yzkyHgMWC5/HvSkJ7N3iK/v9aSnd9/N/In50Td3SXW8+mBpOVNtc93PtSCnECfOtx35TmZE0q1pLlS7uuxtf/VbRyFHxuPBUp6f9GR2UjLT0lKS81ecz+NNzFs0vPMi5ufl0DbpipPwzt/NxCwPuvjtInuSdDLTyzENJqNV7e8zvepVWe1zljJ6UWqk+fnBedOZsTGYZnrUl4YnbbQ83w2PerU47PkV5/NkFuc52bMWtWNOtdTKahxH+6LH+jvYYbaWWllfx+WMeBrDeEyp5lrqUZtWM1JlL3rWoxZTm61OM934hJU8CjkyHqdprueGp3zRg9L82QsrEf2ix9MrzuexF/MvetSscdSPerG+SG9i3kHOX/WsZdN1xWk9SfpKj5pPc8Yc5tbiKT6F/PfODyDjkZlmjeL8+iatSrVslH4jo1XxvnEQ8+drr7503Es9W/vmHsW8g5x/0aMW1x9y07FeJH2TS/pkZ61mPq2e5E/I//v5n8h4hKZUD41/es9v761SPVj4w5syMbYW89dmv9wo981tZcOrmBtppoXTN4O3mZlcqKXkffEt6a/a6FEPtyPptk/zJ+QnQ8ajNS21aXjKy3th1sxC/EtCLEYOYr5p9svPY7GVU89ibqS0lZxPONRyNT9HSX9u8YFX/R5u9KjltPXG/ol+hPzfnn+c9LRO3zTXquEpV8Js8YrSEIuRtZhb+uWhxby1nE881FJ4zxfeJX2ii6MuT/Ug5CxxTsk0s8hLWb0XVqtXVPrVyqyWTl+0tJO60GJ+zqlxDSfcQKil8M77lvRJRdLdntxZyJHxqZlSLS0SEi8SBLVs/HOsCLHk52v+Y36t8utLjhhczI2U6cFqQfjyGm8k1FKYK7+x9AkIuevTOwo5Mj5N06IxWv7l0sfuEmLx75e3EPNe5NMxRfJtNIs+RjN+Oy+P2u8RuLbX0NfgZR6cX9BJyJHx6ZpmFtHy5dUrOoRYrJZaj3659b3mJOYOHxKOM9km1HL0zSfgWbaes7keWkj6JvTIvVx9i5d0EHJkfNqmtHEp8kth/2aHEIvlR8HF5iSLq3AR81V/nnCrUMuLHm4rbl4yb6nmetTK4kPe6DH0aL1ccasXtRZyZPwWrFGYN9fvfrcQi6Vfblz8cicx7y3Icp4d11DLlz4/XmIyzTSvXRy1SlMdu7V9WUshR8ZvxRq37xel3C7EUnnfWFV/MfriFnpwEPNNv0GNVqGWm1wErZzBY77L5eLoJHSo/QtbCTkyfkvWKMzPRYmxCLFs6u4dLSzCEK/VQZrKK7ETc6dIfMtZdQ+1vNxagqLFLOaR9NDj8HQ1HV7aQsiR8Vuzxu37qxIpbw6x1Eu5nV++cvVTrcW85yDLeZbcQi2vWt3O5qHbsm4vdhZyZPw2rSG3vMQ3tgixbOrvISu/vIXvbPkh0XuQ5Twal8JfRkbPeiDQMi1rfkLtjegq5Mj47VqtLFeEOdwXS69eb1PKa9NmMVALi0XZQYIs5/HYF/4yOn5j4G9wItb8hBefQo6M37IpbZLy0lc1hVg2TXFfq/ovLTIWLHeUHtMdB/F/lTmkSh7nfHXL2eZTMZsneBRyZPzWTakea35dIaeNIZbGOLeFX+688Gk5tpNZ1Eb3Nsdu3Y+MNgRaYrbmJ5wbsNU8x1rIkXGscb9nlZTXfwBYLVk25pe/ui985ke2icQ71XrxMM8uFdaPM09GS5Rm8xSfQv6vzg8g4zdtmtXKXkXMu3HDv42UNzekaxnPtmyS4bin1MNMuyUoIubRmd2T/An51+d/IuM3b5rX+sZ9Snlz/ZdWC59GDm2lB/TL8wrnLoGWV/aBxmO2T/Ml5Mg4dmENUv5cLslepDxtlNtWC5/50e2CLC9DbwZX5pRtTq55FGb/RD9CTlAFK1hDVvmqVynvaeHT8uinMwy8vJhnm9u/4KarJ47fXJ7qQchZ4sTKzX2DkJEfKbfITu8ktJa10Z9D/EVo4VSlhQotIzW3J3cWcmQcq7aGrPKHilf5kfJ5Y43ETrFsK+930KXPd2NzC7S8kJw4NnN9ekchR8axOmuIV79UdgPyI+U2C58d7lvNLIIsrVMeO8+9S5V16pqPytxf0EnIkXGsyRryw+saSNSLsFWMu3stl8Yz2NRW34Sqg+2Ynkhy4iiszUs6CDkyjtlYrSR/qa05XieRlYGZq6M0L0129Jitdn5+ab+42nn+U6cl0BvuBjoOa/ei1kKOjGO2VrtB6Ln67muQ8hfbO69x4bNz8MOqZ9FApbVKx+eWa06meTBr+7KWQq6ZHnNjkQRrtJqFx1rPumGPqLUANy58dpdymzyW3po3W44RMR+5BR8AhtWbFpUBjtrahrUS/GrfSleZj5h77Rls8ssDZJcX5tO+QgtiPux7E3oAGNZktYuetSJaK+Vf7IVGs9pl1w5bhN6dw2bf53NocVSGmI/Pgg8Aw2xMWWWk+ktD64g6Kd/Ye7iNy6c+pNxu6TNIdvnVu2GfnIiY9/9+hB4AhtmaHiol7rlhnWZR40E6yO8QUm7VVtpDVN7DOF0yzQct0XtrFnwAGGZvNUuXDcmEtRuLXl2ytNtVS3e+Upvs8oBZLBcza5+cSJ55P+9B6AFgmJvVLnnWeqi1+SuNreAujpTVZrD4kvLUIg7duqCu53cFMQ85+6EHgGGuVr/k2fDaeY0wrlxCIg3JiJ6k3Cq73Nu5PLwziHmIeQ89AAxrY5q3XfKsjbFb7vN8N4p6KfcU9NCsMXM74J7PwmgR86FnPPQAMKyt1US7G/xqpTWRZ+t9nvmx5rXLfB0rsFyMeWlRGz34wufFeO3FnKqJXeY69AAwrL1pVrPLsyHQUOtJP7uJSoPAepNyy06iIxJEJzGnBG7bWQ49AAzrZlpU+uSNWeFa1oRXrPd5no81lJQ3F9QN0oSidsT2qYmIufv8hh4AhnW1Tkueda912OeZH61VO7pW12xTUHckC5/vxoyY9zOzoQeAYT6sdsmzYQmtJjTjtM8zP1p9DyOPC5EWCYmjyWG5mG17Md/QA9RqTkMPAMN8WY2ANgZIfO3zNGrsYeRZWq0WPkcnhU5i/jy+D6NxWfABYJg/q/GrN01SViu+zqmDDVLutNnI4mxNBXVHlcNy8X7ZijmVWermMfQAMMyvVS55vjZLce0+T2eftqEdndNmI4uzNW/fH1UOy7uR21dNJMu8fAZDDwDDfFtNVrhF9knNPs8W23lqs0ocNxt1PNvRPGbMeB67m5iP8gMp2NyFHgCG9WGVYvxiIwA1+zxbxJlrfXxvOz3PZ0sbc1hexhtrdhBzclneZi30ADCsL6sIabzaxFhrPPpWwlu719Nxs5HV6Jv6CY0wh+Xd+DOLuuunuSOXBRnHpmyaV4iBlU9dU8ewVe5HzWYjD63gSs9Xv3w4ail36gEavCtS4JkKPQAM69OUVkjxi91SWWV4pZ1PXrfo6VT13PqMTTksI5dyI+vN/DecyxJ8ABjWt1UGV6wErKZ9nON2/fx49ZuNeliAbGwKPX4pTy3K9Z6u5QZzWYIPAMP6N2UVkeJnu1BGpQ9t6dMXjrfwVfXc+oxNpW5HL+VOWeY3VjEx+AAwbAhTWvHVvHGrfv56zz557Z5T57Jc1nPwGLuU59Vk7HJZbmYrf/ABYNhQVhnntsze9u6T12V5e09EzM/ZVLU8CinPA0V2y583EDEPPgAMG84qferGrfoNr2/vk899luWyPuskpNwhl2WUdWU8zkPoAWDYkFaZD26xVT8/gmeffPhExPNZ66X8IQ7hs25MEcXVtJyD0APAsKGtsuqKpUfdg08+eCJift56KY/Ih1XakCW/Cj3CXq8+9AAwbHjTrOJP3mqrvlEvPnldIqJjf1Cn805GyvNclvKrmXR8PPgAMCyMddmqb1STxGi9ZFpyzPnQiYj5eSck5aVb+V9Dj6nnKw49AAwLZd226teGQjosTtYmIvboU05LyvPlz/M7Gno0PV9r6AFgWDjrulW/xn/uUIS2pixXD2W0Ls48NSk/7f6c+L7O4APAsLDWbat+bz55dVku71XKr848MSk3mv5+zuADwLDQVhnlts5uqPHJO2SZ1AhqbxnljWc+XlV0Uj5tCz4ADAtvlVv1X2zlqsYn77A0WXPUnjbsvzs3Uh6NBR8Aho3DKrbqOyQQ1nQc6hCbrUlE7HRcq3M3Snnodw0zQsYx7GwV23peHYQ87Rpprziu16ZyjueeyMb9KVvwAWDYeKwyR8RBqrruEa04anU/z56KaF2cvUnKJ725ZvwWfAAYNi6rqHDiIMI1kfZOC5OVHxAD+OSNUt5LwwvM8r0JPQAMG5spK5VLp6XKCsntuIGnJqO8161B5/Mj5aO04APAsPGZ0tIouXXeilFNAa2Oy4I1G/Ytexl1PH+9lK+mn6U9Pgs+AAwbo1UERix7BZ2PUe47W1Y3rzludRpir1uDziNAykdlwQeAYWO1UrF0DF5URNqtq5tXHre6WG7PW4PeXVm1lJNVPqgFHwCGjde0LH3YKcWuco9o55zrcFuD3o2guuEFUj6YBR8Aho3ZKqogOjUhqMxc6RzLDrk16HxtdVJOVvkgFnwAGDZuq/CmHSVYD6UPO8XaK48cMA3RqKEiI1I+xDsQegAYNnarzFtxikFXBlc6L0qG3Rp0HkOdlJOK2O/shx4AhsVgpSLlsE3fqKa6uYdEwbBbg/IxZDVt6JDyPmc+9AAwLA6ryFtxDBhULEt6iGSH3hqUj2JeK+W9Nr24XQs+AAyLxSryVhyDIhWLph2LZ52PXZUE2GMnT4dRHEeClPue8dADwLB4rGqbvuNR0oqcbw9SW7M1aLA4uVHjBiFSEf3OdugBYFhMpln3bfpGNcEVD55qzdagQeuD1+TQkFXud6ZDDwDD4rKKhUrn6HY/xbPOR6/yyQfKJz/PVf0GIVIR/cxz6AFgWHzmI2+l0rP35DNXbg0aqO7KeRzVYR4jssr9zHHoAWBYjFaxncdZkiqyS7wEV2rCGoPUQnw3DrLK+53f0APAsDitIiji7ElXFs/yEvyo9PgHTEI8j4RUxL7mNvQAMCxW85O3UtMD1FO4oTKoMfgiY0NWOYuebec19AAwLF6rSB10F/KqrTueMksqCwF4qOniPJZ5ZR7N4HH7qVjwAWBYzFZRvdBZyCuj7S2OVDHOquj0oEmI+Wjqsspp0Ow+n6EHgGGxW6lAthHycp/ZOSe98vhVdVcGajRxNZplTSoii55ucxl6ABgWv5Vu028j5OW+vbdc78riXEEyuBuyytm0bz+ToQeAYVOw0koprXpu+slJrzn+SJIQ89HUpyKy6Gk3i6EHgGHTsNK8lVYBER+t4xpGWpWEGCSUUZu/wk5PmxkMPQAMm4p5FHJPOek1ZxhNEuL5iusWPYmU189e6AFg2HTMo5B7ykmvOUNVOdlBq65cjOihdtEz0KhisOADwLApmUchT/3UUqw9Q1Xzt0DZ2w31V1j0rJq30APAsGmZVyHvNW/FaGxJiEa1nUUpb1s1Z6EHgGFTM39C3n/eSm0lxGDbcFj0dJyv0APAsOmZVyHvOW/FqHIHaaAkxPN1Vy96vrDoeTFXoQeAYVM0r0Lee95KTRJi0CyR2u1B1ER8m6fQA8CwaZpXIe89b6V2eTFA1RWLURmxPSi3xAgA+iDJtFF69eBOd+bQ4lipNsoKD7c8WuVZ5noujFiStvro8zyOo5pppXnFLw/6ZJ5CjWwsfCP0AACmitnpTrurBzO9JFmLYx10p3Xh4UybNkerPMtWX5WcRZrrSzJ3PpyvUe3NXclMHkm1SjbhxjYSQn8dwLApW2n2d+tMk/7zVoxqag8GDK7k46pe9PQaYIrNCKoA9EppOOSgO7NrdbSlViUP3/sNLFSGMYIGV6Qk1YOWpWGfnfkQblzBCf05gmFTN88e+QB5K0aVVVcC9Au6Gld5TcSbTkDEGwfoHc8eednSqfRk7j2Puvw80ifzuaeJsh3ZXA8X3xa8X3tcIOMAA+BZyAfJW5GSVM+lwZW17kMGVyQpWeox/5A56IPZhx1NWMhUARgAcyjJtUjb5pkMk7cimYO5U5nnvWiXb+N1bE/6Kh/b022LON44wGD49cilZKVl4cEOx6s8T1U2ueeF1VZjm2mpz6G/GYQGGQcYDO9CXpa30oeQVwVXnvTp1iV0DCDjAAPiXcjLPOUehFxKHvVQ8vBO9/7PBW4QGwcYEL8xcslsdafryHCql2TZ6nB1Z/qkjyp63pk2/s8FbuCNAwxMkmqlxdWDe31oG56oyFvpIXKdzPSssg8cgitBwRsHGBhzMB91LbEzbZK01eHKPXxp1YNHvjcfSjNXltoksx6mCqxAxgECYO4LQp7pufXRDuaDir53D0JeE1x5SRYtDgceQMYBgmDuC5nf82TV6lCn4w0l5OtS7z/Vc/Lo/2zQDLFxgECUxrQ7bSsvzSTvZaN6kupRZR8RO3289c04w4M3DhCI0pj2sov/XOqRL7v4+NVjN/e6J7gyDpBxgGCYQ4kUdgqEmPuSJchehFwyTwRXxgFBFYCglNYR/NBlS03p3s6eagBWBlcCVya/LfDGAYJidvpYeLBTiSvzpKJkL5NV24TG2nMdzL0+lfwiaNu3WwNvHCA4Jf5zx+KrpR6590K257NlelZZ3vhn88n5YOAM3jhAcMxTwaNN9dzFey71yLP2W4wazrbTh9JWzA9JT2eE9+CNA4yCkmTBjt5zstCqEHXvzSOXkgeVLW3u9ZHiWf2CjAOMhBIh77gwWbp82qeQl7d9O+hT+MrkUwYZBxgNyYvf7UCVQt5bcdnKyuREyXsEGQcYDaX7OjsKYKmQ91KR/HzG8srkpCD2BjIOMCKSVF8Kotux5GwAIS9v+0aUvCfIVAEYEeagu5J9nZ22t5ud31YVFmfc6quS/Z0zWkz0A944wMgo8Z47+86+m8dZnXNVur+TKLl38MYBRobZlWSRd/SdK5vH9egdm3uVLc8+JJ0y4qEI3jjACCnZhdmhzVt+zMGavb07Z3kKIo2YvYI3DjBCzFOhUmGHNm/5MQdr9vbunLvSKHmmDeVs/YGMA4wS88lnm7f8mAfdDdUj6O2cpQ3mKGfrEYIqAKMleda1z+qh4Gzp0mOvoZXKjfpr3ZNL3h1kHGC0+G/zlh83hJCX55ITJfcAMg4wYkqF3IPgBhHymZ5LllgP+mi2fZ53+hAbBxgx/tu85ce9L42R97rsaPalkflUm+ShxeHgDN44wMjpYzuQVNGsotcNQVJllPxJn4iStwUZBxg9yVybq4c6Z5FLpUWshhDyqij5xy79jm4ZgioAo8dsC/shZ/LQ7b4kqbHXWivnq/lQmkv+QvfOdiDjABFQ0uZt4SPvuiRGPoSQ73VX0vSNKHlLCKoAREJJdslHs251qKbjDhBaqaxL7iGh8tZAxgGiodAd6KAPPuLJJULuJfbeeN6ybqG9tpmbJgRVAOLh41XyYSov1QJLQiudK7hYnXddWuUl00vfYZ1pgYwDRIPZ6+PVQ5m8VCYx9wU5zQYR8l1plHymDcud9iDjABFhtoWlzqWnZcG7QEJ+MB8L1Rx7r4U+LYiNA0RGScGsDz4WJEs3/g8Up66IkrPcaQUyDhAZJXJ70Fc+xDaokJc3mNiau/7PHTsEVQAio6TOSlrY5dn2yGWhlUHqglc0mKD6oQXIOEB0mF1hV2fmpwmDOehO1ymMy8TDjlHLc19mzKxpv2wDMg4QIWZdWBZ88LMkaA6FtMYBhdzcv7uunYiMW0FsHCBSko0uk/K87b0sjVMPttyYLPWo1NfWplsAGQeIlCTVi2YXD3nbexlYyDNthigHMBWQcYBoSTK9XD20Nh9bHars2EUh/0yseowQGweIlpKlTi91D/Nj3xVi5J7i7+AXvHGAqOmr7qEkJQs9Fx7suWMnuIOMA0ROX3UPpdJGbwj56CCoAhA7PdU9lCTzVJL056GlM/gEGQeInP7qHkoIeQwg4wDR02PdQ1UKORXBRwOxcYBJ0Ffdw/zoxRj5II3ewAZkHGAS9Ff3MD9+oI6d0AwyDjARSjbs7MwHj8cvE3KPHxTQFmLjABOhv7qH+fGLHTvTIfoDQRPIOMBk6K/uYX78opBnfiqdQxcIqgBMiv7qHubHL4ZWaLUWGGQcYFL0WfcwP36x0Rsls4JCUAVgUpS0fZjJWw55RaM3SmYFBW8cYHKUZHl7K5clSclML4UitlRaCQYyDjBBChFsz6mBJcmNZJEHg6AKwAQx91eBj1Reu2maXaGOS6oNG/TDgIwDTJP7qwj5wm/82mwLWeoeKyuCC8g4wCQxu0IO+WMya3WoqjM8Fc4wYztQCIiNA0yWQg751tx5PkMxi9z7OaAJvHGA6XIdWJn73JwvSeZe1xkw88RrFB6awRsHmDAl3TQ9lq+V2A40BpBxgElTqEPuteqhVCHkZJEPCEEVgGlzr8v2yl6rHkqSORSCN9IqWbQ6GLQAbxxg4iTzQhXCO7P1fA62AwUEbxxg4phtITFw5TstsKTWOduBBgMZB5g+n6/2dM7kObAimXWJkHv/uIAyCKoA3ABJpperh7wWy8rP8liopbjTHW3e+gZvHOAGKNnT2YOnbD6VdAcii7x38MYBboTk5SotcG0+tjxU3VmuExzpDtQ7eOMAt8LHQrEsj+0kztwXmkosezkPnMEbB7gZkoerpc2DPph9y4NVn6XYRq6XSDycQMYBbohCsSzvezolssiHhqAKwC1xHVjxvqdTksxOd1fnoRZ5jyDjADeEORSyux/62KRjdroujkUt8t5AxgFuCrMuJAX2sknHPBWEPPO/6QgkZBzg9vh0XSxLvWSSmM+FD4xlHyEcYIkT4OYYolhWfqbrXPWDvmJXp2/wxgFujiGKZeXcXXj+bM3vBbxxgJtkmD2dF8mHiHhP4I0D3CbXGSuLfho9nEvYPpkPiHg/4I0D3CiFeoS9xa2TpVLzuftxoBxkHOBmKezp3Jq70GMCdwiqANwu1z005xSxihFkHOBmMftCxspDMmt1KAgIQRWAm6ZQH5zASnQg4wA3TZLqy1U1QorKRgZBFYCbpqRYFo2QIwMZB7hxzFqX3nfaT40V6AuCKgA3T0lgpacaK9AHeOMAN485FDJWqEQYEcg4AMh81qX33UtXIOgHgioAICnJ9HLxQC/tlqEP8MYBQJLZXQVWUq1CjwnswBsHAElSkupFl3s4ySCPAmQcAHIKXYHo1RMFBFUAIMdsySCPEbxxADhDBnmM4I0DwBkyyGMEGQeAd5BBHh8EVQDgAjLIYwNvHAAuIIM8NvDGAaBA8oUM8nhAxgGgABnkMUFQBQAKmK2eLh4gg3zE4I0DQAlkkMcD3jgAlGAO+nT1EImHIwUZB4BSzBMZ5HFAUAUAKkhm+nLxABnkowRvHAAqMHsyyGMAbxwAaiCDfPwg4wBQQyGDfG++Cj0muISgCgDUUMggn7HQOTbwxgGglkIGOTs6RwbeOADUUqhBnpJBPi7wxgGgkcJC5wezCz0mOIE3DgDN3F/9G398RCDjANCI2V7t6Jwny9BjghPIOADYcO2Pz1odBXoAGQcACy52dG51Zz51OBh4hSVOALAiTzw86JN56n408AcyDgCWJAtl+kzO+NhAxgEAoobYOABA1CDjAABRg4wDAEQNMg4AEDXIOABA1CDjAABRg4wDAEQNMg4AEDXIOABA1CDjAABRg4wDAEQNMg4AEDXIOABA1CDjAABRg4wDAEQNMg4AEDXIOABA1Px/QCw+nbduclYAAAAldEVYdGRhdGU6Y3JlYXRlADIwMTctMTItMjJUMTI6MDE6MzgrMDA6MDDyHp8bAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE3LTEyLTIyVDEyOjAxOjM4KzAwOjAwg0MnpwAAABR0RVh0cGRmOlZlcnNpb24AUERGLTEuNSAFXAs5AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%%tikz -p stanli\n",
    "\\usetikzlibrary{decorations.pathreplacing}\n",
    "\n",
    "    %the points\n",
    "    \\point{begin}{0}{0};\n",
    "    \\point{middle}{2.5}{0};\n",
    "    \\point{end}{5}{0};\n",
    "    \\beam{2}{begin}{end};\n",
    "   \\support{3}{begin}[-90];\n",
    "    %the load\n",
    "    \\load{1}{middle}[90];\n",
    "    \\load{1}{end}[90];\n",
    "    %the inscription of the load\n",
    "    \\notation{1}{middle}{$F_1$};\n",
    "    \\notation{1}{end}{$F_2$};\n",
    "     \\foreach [evaluate={\\in=180-\\b*2}] \\b in {5,10,...,30}{\n",
    "      \\draw[-, ultra thick] (begin) to[out=0,in=\\in] (-\\b:5-\\b*0.01);\n",
    "% quater circles with sortened radius\n",
    "%      \\draw[red] (begin) -- (5-\\b*0.01,0) arc (0:-90:5-\\b*0.01) -- cycle;\n",
    "    }\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "hide_input": false,
  "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.5.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}