{
"cells": [
{
"cell_type": "markdown",
"source": [
"# Cantilever with a Tip Force\n",
"\n",
"This example shows how to predict the behavior of a cantilever beam that is subjected\n",
"to a constant shear force at the tip.\n",
"\n",
""
],
"metadata": {}
},
{
"outputs": [],
"cell_type": "code",
"source": [
"using GXBeam, LinearAlgebra\n",
"\n",
"L = 1\n",
"EI = 1e6\n",
"\n",
"# shear force (applied at end)\n",
"λ = 0:0.5:16\n",
"p = EI/L^2\n",
"P = λ*p\n",
"\n",
"# create points\n",
"nelem = 16\n",
"x = range(0, L, length=nelem+1)\n",
"y = zero(x)\n",
"z = zero(x)\n",
"points = [[x[i],y[i],z[i]] for i = 1:length(x)]\n",
"\n",
"# index of endpoints of each beam element\n",
"start = 1:nelem\n",
"stop = 2:nelem+1\n",
"\n",
"# compliance matrix for each beam element\n",
"compliance = fill(Diagonal([0, 0, 0, 0, 1/EI, 0]), nelem)\n",
"\n",
"# create assembly of interconnected nonlinear beams\n",
"assembly = Assembly(points, start, stop, compliance=compliance)\n",
"\n",
"# pre-initialize system storage\n",
"system = StaticSystem(assembly)\n",
"\n",
"# run an analysis for each prescribed tip load\n",
"states = Vector{AssemblyState{Float64}}(undef, length(P))\n",
"for i = 1:length(P)\n",
"\n",
" # create dictionary of prescribed conditions\n",
" prescribed_conditions = Dict(\n",
" # fixed left side\n",
" 1 => PrescribedConditions(ux=0, uy=0, uz=0, theta_x=0, theta_y=0,\n",
" theta_z=0),\n",
" # shear force on right tip\n",
" nelem+1 => PrescribedConditions(Fz = P[i])\n",
" )\n",
"\n",
" # perform a static analysis\n",
" _, states[i], converged = static_analysis!(system, assembly;\n",
" prescribed_conditions=prescribed_conditions)\n",
"\n",
"end\n",
"\n",
"nothing #hide"
],
"metadata": {},
"execution_count": 1
},
{
"cell_type": "markdown",
"source": [
"## Comparison with Analytical Results"
],
"metadata": {}
},
{
"cell_type": "markdown",
"source": [
"The analytical solution to this problem has been presented by several authors. Here\n",
"we follow the solution by H. J. Barten in \"On the Deflection of a Cantilever Beam\",\n",
"after incorporating the corrections they submitted for finding the tip angle."
],
"metadata": {}
},
{
"outputs": [],
"cell_type": "code",
"source": [
"import Elliptic\n",
"\n",
"δ = range(pi/4, pi/2, length=10^5)[2:end-1]\n",
"\n",
"k = @. cos(pi/4)/sin(δ)\n",
"λ_a = @. (Elliptic.F(pi/2, k^2) - Elliptic.F(δ, k^2))^2\n",
"\n",
"θ_a = @. 2*(pi/4 - acos(k))\n",
"\n",
"ξ_a = @. sqrt(2*sin(θ_a)/λ_a) .- 1\n",
"\n",
"η_a = @. 1-2/sqrt(λ_a)*(Elliptic.E(pi/2, k^2) - Elliptic.E(δ, k^2))\n",
"\n",
"nothing #hide"
],
"metadata": {},
"execution_count": 2
},
{
"cell_type": "markdown",
"source": [
"Plotting the results reveals that the analytical and computational solutions show\n",
"excellent agreement."
],
"metadata": {}
},
{
"outputs": [],
"cell_type": "code",
"source": [
"using Plots\n",
"pyplot()\n",
"nothing #hide"
],
"metadata": {},
"execution_count": 3
},
{
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/plain": "Plot{Plots.PyPlotBackend() n=11}",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAYAAAByNR6YAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACZE0lEQVR4nOzdd1yVdf/H8dc5DNmyREQEXAiogOA29zY3WrfZXY6yeZtmqZUprjKzUu/ql6VmZXangpZW7pV7kLkAJwKKsmWPc871+wM9iXAUlCmf5+PBg3HOdV2fQwZvvt/v9fmqFEVREEIIIYQQZUZd2QUIIYQQQjxuJGAJIYQQQpQxCVhCCCGEEGVMApYQQgghRBmTgCWEEEIIUcYqNGBNnDgRDw8PVCoVZ86cKfJ4Tk4OQ4cOxdPTE39/f/r160dUVFRFliiEEEII8cgqNGCNGDGC/fv34+7ubvA5EyZMIDIykpMnTzJw4EAmTJhQgRUKIYQQQjy6Cg1YXbp0wdXV1eDjZmZmDBgwAJVKBUD79u25fPmywednZWURFhZGVlZWmdcqhBBCCPGwqvQarKVLlzJo0CCDj0dERBAYGEi/fv0YPHgwP/30UwVWJ4QQQghRPOPKLsCQDz74gAsXLvDVV1898LmLFy8mICCgAqoSQgghhHiwKhmwFi1aRGhoKDt27MDCwqKyyxFCCCGEKJUqF7A+/fRTfvrpJ3bs2IGtrW1llyOEEEIIUWoVugbrtddew9XVldjYWHr16kWTJk30jw0YMID9+/czZcoUUlNT6d69O/7+/rRr164iSxRCCCGEeGQqRVGUyi7iYYWFhREYGMiJEydkDZYQQgghqowqfRehEEIIIUR1JAFLCCGEEKKMScASQgghhChjErCEEEIIIcqYBCwhhBBCiDImAUsIIYQQooxJwBJCCCGEKGMSsIQQQgghypgELCGEEEKIMiYBSwghhBCijEnAEkIIIYQoYxKwhBBCCCHKmAQsIYQQQogyJgFLCCGEEKKMScASQgghhChjErCEEEIIIcqYBCwhhBBCiDImAUsIIYQQooxJwBJCCCGEKGMSsIQQQgghypgELCGEEEKIMiYBSwghhBCijEnAEkIIIYQoYxKwhBBCCCHKmAQsIYQQQogyJgFLCCGEEKKMScASQgghhChjErCEEEIIIcpYqQJWVFQUO3bsIDk5ubzqEUIIIYSo9gwGrClTpjBp0iT95xs2bKBZs2b06dOHpk2bcuLEiYqoTwghhBCi2jEYsDZs2EDr1q31n7/77rsMGDCAU6dO0bZtW2bMmFEhBQohhBBCVDcGA1ZcXBxubm4AXLp0icjISGbMmEGLFi2YOHEix48fr7AihRBCCCGqE4MBq3bt2sTHxwOwfft27O3tCQwMBMDU1JTs7OyKqVAIIYQQopoxNvRAly5dmDlzJjdv3mTRokUMHTpU/1hkZKR+dEsIIYQQQhRmcATrs88+w9nZmenTp+Pm5sb8+fP1j/3www907ty5QgoUQgghhKhuDI5g1a9fn127dhX72NatWzE3Ny+3ooQQQgghqjODI1g9evQgIiKi2Mdu3LhB3759y60oIYQQQojqzGDA2rNnD2lpacU+lpaWxr59+8qtKCGEEEKI6uy+ndxVKlWxXz948CBOTk7lUpAQQgghRHVXaA3Whx9+yIcffggUhKvu3bujVhfOYLm5uWg0Gl599dWKq1IIIYQQohopFLA6duzIlClTUBSFOXPmMGrUKFxdXQsdYGpqire3N4MGDarQQoUQQgghqotCAatr16507doVKBjBevHFF3FxcamUwoQQQgghqiuDbRpmzZpVkXUIIYQQQjw2DAYsnU7H8uXLWb9+PbGxseTk5BR6XKVScenSpXIvUAghhBCiujEYsKZNm8Ynn3xCp06d6Ny5M6amphVZlxBCCCFEtWUwYP34448EBwczc+bMiqxHCCGEEKLaM9gHKycnh06dOlVkLUIIIYQQjwWDAWv06NFs2rSpTC82ceJEPDw8UKlUnDlzxuDzLly4QMeOHfH09KRt27acO3euTOsQQgghhChPBqcI27dvz4wZM7h58ya9e/fG1ta2yHOGDx9eqouNGDGCqVOn8sQTT9z3eS+99BITJkxgzJgxrF+/nvHjx3Po0KFSXUsIIYQQorKoFEVRinvg3g7uRQ5UqdBqtQ91UQ8PDzZv3kyLFi2KPBYfH4+npyeJiYkYGxujKAr16tXj8OHDeHh4FHpuWFgYgYGBnDhxgoCAgIeqRQghhBCirBkcwbpy5UpF1qEXExODi4sLxsYFpalUKtzc3IiOji4SsO6YNGkStra2jBo1ilGjRlVgtUIIIYQQRRkMWO7u7hVZRyH3bjJtYJBNb/HixTKCJYQQQogq4/7zgMCWLVuYO3cuEyZMIDo6GoB9+/Zx/fr1cimoQYMGxMbGotFogIJwFRMTg5ubW7lcTwghhBCirBkMWAkJCXTs2JEnn3ySZcuWsWLFChITEwFYuXIl8+fPL5eCnJycaNWqFatXrwYgJCQEDw8Pg9ODQgghhBBVjcGANWnSJBITEzl9+jRRUVGFpul69erFzp07S32x1157DVdXV2JjY+nVqxdNmjTRPzZgwACOHz8OwLJly1i2bBmenp4sWLCAFStWlPpaQgghhBCVxeAarN9++41vvvkGHx+fIncL3pnGK60vvviCL774otjHfv/9d/3HzZo1k7YMQgghhKi2DI5gaTQaLC0ti30sJSVF9iYUQgghhDDAYMBq164dK1euLPax//3vf7KNjhBCCCGEAQYD1rx589i8eTNdunThiy++QKVSsXHjRkaOHMmvv/7K7NmzK7JOIYQQQohSCw0Nxa9VIOYWlvi1CiQ0NPShjy8NgwGrQ4cO7N69G5VKxZQpU1AUhfnz5xMXF8fOnTul75QQQgghylVZhKOgoCBO5zuQMzCY0/kOBAUFlfg89x5fGga3yrlbdnY2KSkp2NraYmFhUaoLlCfZKkcIIYSomkJDQ5k9dz7nIyPwbObFrPffK9UexnfCjap5bxSf3qjObUc5u52QkJASn8evVSCn8x1QJm4GlQoUBdWSJ2lklMLnvx0lUwNZGsjMVwre3/lcA1kahfXj23LL1BHeKDheedGkxPUbvIvwbubm5pibm5f4pEIIIYSovso0HA0M5vS57QQFBd03HOkUhYx89G/Tg+cXHH87HCm93kC15Elee28+R+oPIUMDGfkFx6TfOU6jkJ53V1A6FwFDggvCFRScp3kfLv0STP8t/3RIUKvAwhgsje9+ryI9JhIGj/rn+FK4b8A6fPgw69atIyYmhpycnEKPqVQqfvnll1JfUAghhBDlp6LCkaIoZGrgVh6k5UFavkJaXsHnk2fMB5/C4YglTzJu+ny+Nh98O0QpZGj+CUdZmnsKiSw+HN34JZiQKB1WxmBlosLaBGqbQn1LsDJWY2UCViYFQemTxl5cP7u94Pp3RrDObqOZlzfbRhnrA1Uto6Lb9AH4eXtx+txdx5eCwYC1ZMkSJk+ejJOTE40bN5a2DEIIIUQ5q4hwpCgF02Fp+XfCkXLXxzDjvXlFw9HiJ/n31PnM0g3SPzc9H3SGFhldKhqOaN6H9F+DsTSGuuZgbXI7DBkXBCJrE5U+HFmbwAv/9eLiuaLhyLeFNyefLtlUnfu89wq+H0ueRGneB9XZbSjndvBBaCgNrB4cmGa9X/h4Xny7RNeF+6zBatCgAcOGDWPx4sWo1Q/csrBSyBosIYQQVUVFrTnSKQUjRal5kJILqXnK7fcwY2hr4tT/rBlCUWDxk5jlJlN33uGCEHW/YATwum1BOOo96Z+vbV+M+tdgXt+Vho0p2JhAbVPVXR+DjakKGxOwMYWu7Vtzppi1T75mKZw8cbx03w+fXoXCUWhoKMOGDSvV93XOvA+IjAinmZc3s95/76GPz87KLPFxBgOWg4MDa9eupWfPniU+WUWTgCWEEKIsVPSC7FztP6EoJVchNQ9e6teGWBxQ7g5HS57EKjeZpgsP6wPVrTwwmI8MhCOjX4OZ/qehcPTPx9am0DqwdbELw6tjOKpMBgPWK6+8gqWlJYsWLaromkpMApYQQohH9Sh3q2l1BeGoU7vWnNcWDiUseRIHTTI9vjhCUi4k5Sgk50JSbjHrjeC+I0cTdqRhawp2tcDWVHX7feHPa5tCoISjKsNgwMrNzWX8+PHodDp69eqFra1tkeeUJt2XBwlYQgghHnX0ydCt/G6qFKatPUJSDvqAlHQ7ICXf/jgl9/ZokoFwpPo1mB4ht3AwA4daqtvvwb7W3SFJha0p9O/cmrMSjh4bBgPWX3/9xbBhw4iOji7+QJWqyCbQFU0ClhBCVH+PEpDuN/o0eOgwEnMgPhsSchT9+4R7vnbwqdoog4OLhCN+Ccb4i1Tsb4eiOwHJvhY4mKkKvnb7/duD23BFkXAk/mEwYAUEBKDRaPjwww/x9PQs9i5Cd3f3ci/wfiRgCSFE9Vaa6TlFUUjPhxtZcCNb4UYWTB7YmuvFLOo2yk5G++7hIterZQROZlDHHJzMVNQxgz9ebkuiceFzqJY8SXPTFE6FHSv29n2Dr0PCkbjNYMCysLAgNDSUfv36VXRNJSYBSwghKld5Tc/V1SUzfPlRbmQp3MhG/77I2iUDU3PGm4L56q90nMyhjhk4mReEKSuTov2OJByJ8mCwD5avry/x8fEVWYsQQohq5EE9l3K1Ctcz4VqWwrVMuJZZEJLi7gpNZwx02r7xSzD7b+hwNlfRxAaeqKvG2QLqWahwNgfn2++7fnlPI8jbvZKa+3gz3qtkLYaGDx9OSEhIQTjaHFwQjkoZru6cp7LXJouqw+AIVlhYGOPHj2fhwoV0794dY+MS7apToWQESwghHt7Djj4pSsGdc+3btOaC9p62Ard7LlkHHyah8AYgWBiDiwU4m6twvv1+/fg23DRyLHSO0qxdKqvRJyHKmsHU1LVrV/Lz8+nXrx9qtbrIXoQqlYpbt26Ve4FCCCHKnqHRp/Xr19P9yeFEZ8LVdIXoDIXYe0ahrmXdnqorZisTWvQh/9dgJjZXU99SRX1LqG+hwtWqoM/SvdNz3T+YUWyn7VmhoSV6HWU1+iREWTM4ghUcHPzAhX2zZs0ql6JKSkawhBA11cOOPuXrCqbtenZszWVd0dEndVYyuvf+WRxuqgZXSwqFpfqW4Gqp4p0hbbj8iHfO3XktsnZJPG4MBqzqQAKWEKImut+ddwMGDyMqA66kK0SlK1zNgOgMhegMuJqhcD3r9jYp91kc/tPZdNyswN1KRR1zUBv4Y1um54QwrEQrAGNiYjh48CCZmSXfg0cIIUTxQkND8WsViLmFJX6tAgkt4XQYFKx/mjFnfkG4mrgZek8qeO/di6ffmo/5txq812kYsEXLxIM6fr6k41omNLKG55uq+eoJI7b0M8KzmReqc9sLRq4KTqxfHD6ikZq2TmrqWqgMhiv4Z3rOt1YKZpuD8a2VIuFKiNvuu3L966+/Zvbs2cTFxaFSqTh27BgBAQEMGzaMbt268cYbb1RUnUII8Vh40J13ANkahUtpcDld4XKawuV0br9XiEqH7PDi1z4pm4JZ0cWIRtbQyEZFfQswUhcfkD4Mfu+R1j7dIXfOCVE8gyNYixcv5j//+Q/PPvssf/zxB3fPJHbr1o21a9dWSIFCCPE4mT236OiTyqcXE96ZT6/fNLitycfiWw0tQzQM2ablnWM6dlzTAdCrvpoP2qjxaFr86FMLH2/GNVPTzUWNm5XKYLgCGX0SorwZHMH673//y/vvv8+MGTOKbInTrFkzIiMjy704IYSoakqzuFyrU4jJhAu3FM7fUriQBmfCIwq2Zbmn71Pyr8HY1oLRTdR41lbRtDY0tlFRt5g1UG5zZPRJiKrOYMC6du0aHTt2LPYxExMTWY8lhKhxDE3v/W/denx6DCc8ReFcqkJ4qsK5lIJAlXv771NjFTS2AUvXZqSf3Q73NMb0be7N+l4l6zcorQmEqPoM/t/s7u7O0aNH6dGjR5HHjhw5gqenZ7kWJoQQVc2su6f3VKqC7uGLn+Rfb30A7w0GwMkcvG1VPOGs5gUv8KytwrO2CncrMFarCDV+tL5Pd8jokxBVm8E1WC+++CLz5s1jxYoVpKWlAZCfn89vv/3Gxx9/zEsvvVRhRQohRFko6d17ORqFsESFVed1TDmspd8fGtx/yufMuQgUn95FFpcb34zgz0FGJP7bmJvPmrBnoDH/94QRb7Qwon8DNY1tVBjfXg8la5+EqBnu2wdr4sSJfPHFF6hUKnQ6HWp1QR579dVXWbp0aYUVaYj0wRJClJSh3lH/t3o9Lk8M41SSwukUhVPJCudv3e4VRUF7gxb2KrxtVfw4pi3X1I/eWFMI8fh7YKPRy5cvs2PHDhITE7G3t6dnz540bdq0ouq7LwlYQoiSauEfyDlN4XDE4ichKxneO4xdLfC1V9HSToWvg4qWdgXBysrknwXm0lhTCFFSD1xR2ahRIyZMmFARtQghhEGluXsvOUfhr6SCab6wJIUTCQoXDPSOMtkUzOVRxtS3LLpP3r1kcbkQoqQMBqyff/6Z6Oho3n777SKPLVq0CHd3d0aOHFmuxQkhBNy/OecT/YcRlqhwIvF2oEpUiMooOM7SGFo5qOjfQE1OYy9iz20vWJh+1917Pj7euFrdP1jdTRaXCyFKwmDAWrBgAWPHji32MXNzcxYsWCABSwhRIWYbuHvv6bfmo0kaBEBtUwhwUBHUUE2go4oARxVNbP7pZN51Xtn0jhJCiJIweBfh+fPnadGiRbGP+fj4cP78+XIrSggh8nUFU3tfntMWNOcs5u49JS6CtT2NuPS0MSnPGbNroDGL2hsxqomaZraFO5nL3XtCiIcRGhpKYCt/LC3MS3WcwYBlZmbGzZs3i30sLi4OY+OSNcQTQgi4f4sERVGISlf4+ZKONw9p6fSrBptVGlpv1PDGQR216jeDs8VvDTOykZpGNqoHrp+CgpB1Muw42VmZnAw7LuFKiBrg7oAU2Mq/VJur31meYJVynSmt3Ut1XYMBq2vXrixYsKBIx/bMzEwWLlxIt27dSnUhIUTNdeeH1Ol8B3IGBnM634GgoCBeWLyep3dqcF2joeH/NPxrl5YNUToaWMIHbdQcGGxE2hhjVn88A85tR7XkSdi+uGCa79wOZs14r7JfmhCiHD1KOLpz/N0BySrlOkFBQcWeR1EUFE0+uqwMtKmJ5CdcY97MGXRxd2T1kFZMCPAo1bUNtmmIiIigQ4cO1KpVixEjRuDi4sL169dZv349eXl5HDx4EC8vr1JdrKxJmwYhqoeWrQI5m1+0RYIqK5mOS4/whLOKjnVVtHdS4WRe/EhUaGhowd17EeEFd++9/56MQAlRTkJDQ5k/dw4RkZF4NWvGe+/PLPXNHY96jjvhqIt7Hbq62bM3Ool9VxNZt/p7hvbrjZKXg5Kfiy634L2Sl4OSd/f7XLq8+g61tdmsHhqASqVCURRGbzzBLZ2aba8PR3fPceh0hWrw/GIHb3Voog9Xrou3lLj++/bBunjxIrNmzWLXrl0kJSXh4OBAr169mDVrFk2aNCnxRcqLBCwhqqaUXIV9cQr7bijsv6Fw9OnaBS0Sek/650nbF1NrczA5WbKvqRBlqeyDTTL7riYQEhJS4vOErF/PiJEj6eLuSFc3B304WvPZhwzp0gElL7sgGOXloORmo+Te+TwbJTcHXW42vRZ+h52RwuqhgYXCUWpOPr+P6mD44mo1KlMzVKZmNFmwlrfaNyk0+vR1WBSfHLlM7JrFqE3NUJnUQlXr9nvTWrePrYXa1IwOQ57GOjOR1UMKAlqZBayqTgKWEBXjQT2o0vMU9t9U2HVdYff1glYJCtDAEjo7q9j7eluuqx2lA7oQD1AZ4UjRaW+HmiyUnCza9uyHzV2hoiDYhJFuYsnuRe8WhKE7wSgvRx+I7nyuy81mwA/7sDMzuW840gebWuYFQce04GNVLXPUtcxoMH4GU9o1KhqOjkaR9PfBguffPofa5J9gpDI20T8/sJU/VinXWT2k1T91/BJGlr0rx8P+KvH3s7O7I93cHJizL6LE/y0kYAkh7svQFjOzvl5Hvv8wdl9XOJqgoFXAxQK6u6jo4aKmez0VDW1Uhc8hHdDFY6wiw5Gi0xWMAuVkoeTcfp+bRYchT2OTk3pPODpBmqoWO2ZMQMnNun1MVkEoyskqmBq7y73TYlAQbBYdvsjl4LEFgaaW+V2ByAy1qfk/YamWOXV7j2RK24ZFw9Gxq6SnJBYcqza67/fjUcPR3d/TOwFpT3QSf15NLNXPntDQUD6YN5fwiAgys7JLdAyUYIpw1apVnD9/npycnCKP//rrryW+UHmQgCVE+fNrFcjpYtZPkZVMnbmH6e6ionu9glDVtLbhbuiyhkpUZRUVjhSdrmAEKCcTJTsTXXam/uOOz7xA7bz0ouFIMWbrGyNvH5eln1IrjuFwdInoz98pCEBmFqjMLFDXskBldvvzWha335vTYfBIrDMSCtdRymBTVcLRnfPcCUjeXl689/7MCvnZYzBgHTt2jK5du+Lu7s758+fx9fXl1q1bREVF4erqSpMmTdi1a1e5F3g/ErCEKD9xWQrbYhXG+VqjGxxcZP2U6eZgcjIzStQeQYjyVGHhSFEKFlPfE4x02Zl0/NfYYsPRLZ0R2/4zvOD52ZkouVn/tBu5h8FwdOQy177/+HYoMr/rveVdgckcdS0L2vbog/WtuEcKR2U16lOdw1FZMBiwunfvjpubGytXrsTExITjx48TEBDAwYMHeeaZZ1i2bBl9+/at6HoLkYAlRNnJ1SocuKGwNVZhS6yOU8kFXzdf0I5sM0d4Q9ZPiaqntGuOFK0WXXYGuqx0dFkZKNkZdBj6L2xybxU/cjRxxO0glYUuOwN02mLrMBSO9IupzSxQm1uiMrdCbWZZEIrMLVGbW+k/bt2uPdaP0ahPdQ5HZcFgwLK3t+enn36iT58+GBkZsX//fjp27AjAt99+y9KlS/nrr5L9By8vErCEeLD7LVC/kaWwOVphU7SOHdcUsjRQ1xz6uqro66qmd30Vf/6xQdZPiXJV2hEoRadDyclEl5VB2649Ct3lpQ9HajN2vv/S7TB1V6DKzSpyvvuOHP3wSUEwMrMoCENFPrZEbW5B6/adsE59fMKReHQG27GrVCpMTU1RqVQ4OTlx9epVfcBydXWVrXKEqAYMbZL8r4/XcanJUI4lKKhV0NFJxfut1PRroMbXHtSqolvMzJn3AZGbgwvWT0m4Erc96vTc+rX/Y+TTo+ji7siU1u7sjb5GUFAQP8yYyJP+zdBl3LodjgoCki47HSU7Uz/NFnH+Am91aKKfqlapVHRzd2TR4Uvo8nJQW9li7NQAtYUVagvrgmBkYY3awgrV7a957e7F3ujrvNjKXR+O9kQn4dO8OXZP/adEr2PGzJkEBQUx+pewwuHos69L/L248//aB/Pm8snxgnAU+tnXpf5/TTYkrxoMjmB17tyZ559/nhdeeIFhw4YRExPDqlWrMDEx4dVXXyUhIYFTp05VdL2FyAiWEPdnaIG6OiuZ4cuPMNhdTf8GKhzNZB1VTVRejSB/XvZfBndujy4zDW1mGrrMNHRZaegyCr/XZqYx4Pu9xd7OfytPx7ZJT6G2tCl4s7BGbWFdMMVm8U9I6jjsmaIjWDJyJKoAgwHrhx9+4OrVq8yYMYPw8HD69OnD9evXAbC0tGT9+vX06dOnQou9lwQsIYqXma/wR4zC0z7FL1CXBp/iQWuXFEUpuGst4xbajFR0Gano0lPRZty6/T6VHu99hq1ae/9GkCamGFnaoLawQW11z3tLG5z7Pl387fzHr5bolngJR6KqKnEfrIyMDA4dOkR2djbt27fHycmpvGt7IAlYQvwjLU/ht2iF9Vd0/BGjkK0Fsw/bkWMuC9QfRw8z+qRo8vWBqV3vJ7HJTi5615tWzZYX+6HNSIX8vMInUKkKRpOsbDGytsX9zc94q7hGkMeiuHU9uuC5pmb3ramsbueXcCSqGoNrsO5lZWVF7969H+liFy5c4PnnnycxMRFbW1tWrVqFj49Pkedt3bqVd999F51OR35+Pm+//TbPP//8I11biMdRaq7Cr1cLQtW2awq5WmhbR0VwoJqghmr+dphRsAZryZOFFqjPKuWGqaJsleXU3N3rln5a+jGD2vmjTUtGl5aMNj3ln48zbqFkZ+jPEXnxUrFrlz45chmL1j30Ieru92pL60LNIb2Xby127ZK3tzfGdiX7I/y998tm7ZKsORJVTaERrNLuUl3af9A9evTgueeeY8yYMaxfv55PPvmEQ4cOFXqOoig4Ojqye/dufH19iYqKwsvLi4SEBKytrQs9V0awxOOuuDsA+w0exq9XFX68qGNrrEK+DjrVVRHUUMVwDzXu1qoi55AGn1VHiXsuafLRpqcWhKO73nRpyXR7ewG2qnyDU3MqcyuMbOwxsrFDbWOPkbUdRtZ2qK1qo7a2xcjKlvYDh2OddlPWLglRTgoFLLVaXfIDVSq02uL7gRQnPj4eT09PEhMTMTY2RlEU6tWrx+HDh/Hw8NA/707A2rBhA126dOHUqVP079+fK1euYGpqWuicErDE46zIFjVnt6Oc206t134m128YHZxUjGqsYnhDNfUtZZF6RXnU0aeAVv5Yp1wrtq3A9mlj0N5KRJeWjC4zrfCBKnVBOLKxp+F7y3irXeNityHJuJWKyrRWiV6HhCMhyk+hKcIrV66U24ViYmJwcXHB2LjgkiqVCjc3N6KjowsFLJVKxdq1axk+fDiWlpakpKQQGhpaJFzdbdKkSdja2jJq1ChGjRpVbq9BiIo0e+78gnB1+w5ApdcbsPhJ7HZ+xIEPRtLIRkJVRSs6NXedoKCgfxaG63To0lPQ3kpCm5qI9lbiPe+TiDhzxmBbAZWJKbUat8SotkPBiJONfcFIVG171JY2+uk57x/3G5ia8ypRuAJpCSBEeSsUsNzd3cv1YvduqVHc+nqNRsOHH37IL7/8QqdOnTh27BhDhw7l9OnT2NvbF3vexYsXywiWeGxcTVf47oKO0+ciUIYEFyxOh4L3LfqQujlYwtVDetTRp3mzg+ni7qhfkP1iK3dGbzzB7DdeocPpjWhvJRXu9G1kXBCWbB0xqu2ISYOmeHrsYm90UtGeSy1a4Dh+ZonqKIt1SyDhSIjy9MBF7ufPn+fo0aPExcVRr1492rRpQ7NmzUp9oQYNGhAbG4tGo9FPEcbExODm5lboeSdPnuT69et06tQJgDZt2uDi4sLff/9N9+7dS31dIaqDbI3CxiiFled17LymYGEMtu7NSD23vWDk6s4dgGe30czbu7LLrZYeOPqkKOgyb6FNjkeTEo82+SbalHg0yfEF71NuEnHu3H0XhhvVdrwdphwwsq1ze9Sp8NKLmVqHKtOQUghRfgy2acjIyGDChAmsXbsWnU6HiYkJ+fn5qNVqRo4cyTfffIOVlVWpLtatWzfGjBmjX+S+aNEiDh8+XOg5N2/epGnTphw7doxmzZpx8eJF2rZty+nTp6lfv36h58oaLFGdKYrCiUSFlZEKP13SkZoHnZ1VjPVUM7KRim2bZIuaslTQDqCYtU+KMVteHog2JQElP1f/fJVpLYzs6mJk54SxvRNGdk50fX1G0T3rSrkwHGTdkhA1gmLAmDFjFGtra2XFihVKSkqKoiiKkpKSoixfvlyxtrZWxowZY+hQgyIiIpT27dsrTZs2VQIDA5UzZ87oH+vfv79y7NgxRVEUZc2aNUqLFi0UX19fpWXLlspPP/1U7PlOnDihAMqJEydKXYsQ5S0kJETx9Q9QzMwtFF//ACUkJERRFEVJzdUpS09rlJbr8xS+zlNcVucp7xzVKJEpumLP4dcqUDEzt1D8WgUqoaGhFf0yqoyQkBAlwN9PsTA3UwL8/fTfz7vptFolPzleybnwt5JxaIuSuvlbJXHVB8qNTyYq5iZGyvudmykxb/TVv73fuZlibmqspIT8n5K2O0TJOvmnkht9XtGkpyo6XfH/PQCls7uj8n7nZkpnd0cFqNH/XYQQxTM4gmVjY8NHH33EK6+8UuSxL7/8knfeeYdbt26Vb/p7ABnBElVVkTsAz21HObudnnPXcqj+UHK1MMRdxQteBRsqG6tlTdX9GNqSZfWcqQzwcUebFIcmMQ5N0k3Q5uuPM7J1xMihHsYO9egxcym1NZky+iSEqBAGA5aTkxPff/89/fr1K/LYli1bePbZZ0lMTCz3Au9HApaoqgztAWick8z764/yQjM1LjWstUJpFpgrt9dDaeKvoYmPpdPoF6mtySra9ylXw7ZJT2HsUA9jRxeMHZwxcnAu+Ni+LioT00LXL4u2BEIIURIGF7mPHTuW//u//6Nv376F7v5TFIUvv/ySsWPHVkiBQlRHEZERKAODi9wBaLQ5mJkBRvc99nFkaIH5+rU/M/iJtuTfjEYTH1sQqBKukR8f+0/XcZWK8zHXeat9MYvLj1/FeXrJFofLwnAhREUyOIL10Ucf8fnnn1OrVi0GDRqEk5MT8fHxbNq0idzcXP7zn/8U6mk1efLkCi0cZARLVC2KonDwpsKnp3WEvtAWrGQPwDsC/P2wTr1eZIH53ZsCqy1tMK5TH2MnV/17EydXjB3q0bpdu0fer04IISqSwYBVnl3dy4oELFEV5OsUQq4UBKtjCQrNakO3G7+w7I2Rj80dgCWd3lM0+QUjUHFXC0alblwlP+4qjWZ9y1sdmhTtPH70Cslnj2Nc1w0jq9r3vb5M7wkhqhODU4Q6na4i6xCi2knNVfgmQsfSszpiM6Gni4rf+hrRr4EKtWoEfVxDCvYA3BxcsAdgNQ0Dhqb31q3+noGtm5N//Qr51y6Tf/0y+TdjQKsBQG1jh0ldd8yatcKz4c5im2t6+/hQq3HLB9Yg03tCiOrG4AhWdSAjWKIyxGUpfHZax/+F68jTwjNNVExuYYSvw+O5aD3A3x/r1KL9o/QbC5vUwrieB6b1G2Hi0hATl0YYO7thZGmjP4eMQAkhahqDI1gnTpwgNTWVnj17ApCSksLUqVMJDw+nV69ezJw5s1TTiEJUJ6GhocyeO5/zkRF4NvNi1vvv0arXMD4+pWPleR211PCf5momNlfjbPH4BCtFqyH/RjT5MefJi75AXuwFIs4Wv3feJ0evUPedbzCu46LfI88QGYESQtQ0BgPW5MmT6dmzpz5gTZo0iY0bN9K7d28WLVqEkZER77//foUVKkRFKdTDamAwp89tJygoCPUrP+PQYRizAtS86qOmtmn1CFaG1k8pOi2amzHkxVwgL+YC+dHnybt+GfLzQKXGuG4DTBs0xdOjgcHpPZO6DUpch+x7J4SoSQxOETo6OvLDDz/Qv39/srOzcXR05PPPP2fs2LF88cUXLFmyhPPnz1d0vYXIFKEoD4Z6WLmQzIVTx7Ewrh7BCu5eP+VIVzcHfYPOleMG0cdBjZKXA4CxkyumDZpi0sATU7emmNRvjLqWeaFzyPSeEEKUnMERrKysLCwsLAA4cOAAubm5DBkyBABfX19iY2MrpkIhKpihHlbJm4OrTbhS8vPIi73InGlT6OLmoF8/9WIrd0ZvPMGSnccZsXIpJg2aYuraBLW5pcFzyfSeEEKUnsGA1ahRI/744w+6du3Kjz/+SGBgIPb29gDEx8djY2Nj6FAhqqXIVIXgMC15dZrB2e3Q641/elid3UYzb+/KLtEgbWoiuVHh5F05R15UOHmxl0Cbz/momOLXTx2/inWPESU+v0zvCSFE6RgMWG+++SYvvPACK1asIDk5mR9++EH/2J49e/D19a2QAoUob5fSFOaEaVl9UaG+Bbw05d2CHlZLnizUw2pWaGiF11bc+qlhw4ahTb5J7sVT5F46Te6l02iTbgBgZF8XUw9vbAO7Y9rQB689I9gbHVd0/ZSXV4W/FiGEqEkMBqxx48bRpEkTjh07RkBAAN27d9c/5uDgwBtvvFEhBQpRXuKzFWaH6fg6XEcdc1jSQc2LXmpqGVWNHlZ3r58q6D91jaCgIJY/1YW+9cxBpcKkXkPMm7fDtHFLanl4Y1TbodA5ZsycRVBQEKN/CSu8fuqzkm0vI4QQ4uFIHyxR42TmF/Sx+uiUDiMVvOuv5j/N1ZhXofVVmuSbtGnfEZvcW0U2OE43tuDgLz9Tq1EL1BbWDzxXaGgoH8ybS3hEwfqpO6NgQgghyk+hEaywsDC8vb0xNzcnLCzsgQdLqBHViVansOq8wswTWhJy4HUfNe+1UuNgVvnBSpeTSe6FU+REhpEbGYYm4RqRV64aXD9l3qJDic8t66eEEKLiFQpYrVu35vDhw7Rt25bWrVvrf7DfS1GUStt/UIgHKa5JqHmbobx9RMvZFBjVWMX81kY0tCnfYHW//fsUrZa86EhyI8PIiQwj72oE6HQYOdTDrFkrbAaOxWv3dfZG35D1U0IIUQ0VCli7d+/Gx8dH/7EQ1Y2hJqG8/DNdnxzO0aFq2tQp/x0IDO3ft+ajYPrUNSUn4gRKdgYqcyvMmvphG/QaZs1aYezooj/HjFnBsn5KCCGqKVmDJR4rhpqEuqtTuHL6mMFR2bIW2Mofq5Ti9+/bPvU5zHzaYubdGtMGnqiMDG8zI+unhBCiejJ4F2FmZiZhYWHcuHEDlUqFs7MzgYGBmJubV2R9QpTKeQNNQm9uDq6QcKVoNeRePEXEuXNMadeo6PqpY1epO+W/JT6frJ8SQojqqUjAyszM5M033+T7778nLy+POwNcKpWKWrVqMXbsWD7++GN9l3chqorTyQrqehXfJFTJzyMnMozsv/eTffYwSlYGjR2s2Xu1mP37vGX9lBBC1ASFApZGo6FXr14cO3aMoKAg+vTpQ4MGDVAUhdjYWLZs2cKyZcs4deoUe/fuRa0u/7UsQjxIlkZhTpiOT07pqDtsOlmfPF3uTUJ1uTnkhB8j++/95Jw7hpKbhbFTA6w6DcTc7wlmdf6bESNGyPopIYSoqZS7rFixQjExMVG2bt2qGPLHH38oJiYmyqpVqww+p6KcOHFCAZQTJ05UdimikmyN0SqNfspTaq3IU+ac0Cg5Gp0SEhKi+LUKVMzMLRS/VoFKaGhoqc8bEhKiBPj7KRbmZkqAv58SEhKi6PJylayTfyqJK+cqsW8NVmLe6Kvc+OgV5daW1UpeXFSx5whs5a9YmJspga38H6oOIYQQ1VOhRe79+/enfv36LF++/L6h7IUXXuD69ev8/vvv5Z3/7ksWuddc8dkKkw9pWXNJoXs9FV89YYSnbdmssbr7DsCubvbsjU5i39VEvh7ajv7utpi4NsGiVVfMfTthXMflwScUQghR4xSaIjx16hTjx49/4EH9+/dn4sSJ5VaUEIYoisL3FxQmH9aiAr7tasTzTVVluoB93tw5dHF3ZPWQVqhUKl5s5c7ojSf48nQcY75ah0ndBmV2LSGEEI+nQouokpKScHF58F/k9erVIykpqdyKEqI4cVkKQ7ZpGbNXy4AGKiJGGjPGU11m4UqbkUr67hAizp6hq5tDkTsAz1+Pl3AlhBCiRAqNYOXl5WFsbLBzg56RkRH5+fnlVpQQd1MUhZ8uKbx+UIupGn7pY8Rg97K5wULRacmJCCPryFayzxwGVDSt58Te6GLuAJQO6kIIIUqoSJr66aef2L9//30Pio6OLreCRM127zY3k6a+yybHIWyIUhjVWMV/OxqVyd6BmqQbZB7ZStbR7WhTEzGp50HtwS9g0boHszrtkA7qQgghHkmhRe6labtQFfYilEXuj5dC29z49EZ1djvKue3Y/OdnVk4eQVDDRxu1UnQ6ciKOk7l/Eznhx1HVMscioBuW7fti0sCz0FSjdFAXQgjxKGSrHFFlGNrmpnmtFM78dbzE57l3k+V33n6bvvXMyDiwGW1iHCauTbB6YhDmAV1Rm5qV4ysSQghRU0mnUFFlnI+MQPHpXWSbm0uR4SU+x51RMKuU60xp7Y5VciwjRz/L/z77gFruXtSZ9BlOU/6LZfu+Eq6EEEKUGwlYokrQKQp27re3ubkzqPoQ29zc3WJhQoAHq4cG0tnNga+u5mP/72nU8vCusA2fhRBC1FwSsESli89WGLBFS1y36XBuO6olT8L2xQXb3ZzbwawZ7z3wHIpWQ+bxnUScPVtsi4WICxfK+2UIIYQQehKwRKXaf0NHq1ANJ5MUtr43gpCQEHxrpWC2ORjfWimEhobed3G5LieL9N0h3Jg7lpTVH9PU2YG90Un6TcqlxYIQQojK8OCmV0KUA0VR+PS0jmlHdXSqq+J/PY2oZ6EC1+EMHz78gcfrcjLJ2PcrGXtC0eVmYxHQDeseI5jVJUxaLAghhKh0MoIlKtytPIURO7S8dUTHmy3V7HzydrgqAV12JmlbfyRuzvOkbV2DeUA3nGd8i/3otzCp58Hw4cMJCQkhy96VT45fJcve9YGjYEIIIYQhoaGh+Af4Y25hXqrj7juClZ+fz3fffcfhw4eJi4ujXr16tG/fnueeew5TU9NHKljUTKeSFIJ2aEjIgQ29jRjqUbKMr8vKIGPfRtL3bkTJz8WqQ3+sez6Fka1jkecOH16yUTAhhBDifu7cmW7f0gmXYY1LdazB327nz5+nWbNmvPzyyxw/fpz8/HyOHz/Oyy+/TLNmzYiMjHzkwkXNsu6yjva/aLAygRPDjIsNV6GhoQS28sfSwpzAVv6ErP2Z9F3riJs7hrSda7Fs25t676/CNujVYsOVEEIIcbe7R6D8A/wJDQ297/MVRSE9L4PotFjemfUudi2daD61Ha4DShewDI5gvfTSS5iamhIREUGTJk30X79w4QKDBw/mlVdeYdeuXaW6mKiZdIrCrBM65v2lY1RjFSu6GGFuXHRK8M5fCl3c6zCltTt7o68x4ul/8fXAAEY8Pxab3qMwqu1QCa9ACCFEdXTvCFTM6TiCgoII/mo+Lbv5kZKTSlJ2Cik5qSTnpJKSnUpKTip5uoL9li9euIhHkNdDtfcx2MndwsKCH374gaCgoCKPrVu3jueff56srKxSX7AsSSf3qi89T+Hfe7T8elXhwzZqpvqpDf5DDWzlj1XKNVYPCdBvsjx64wkyaztz4tSZCq5cCCFEZQsNDWXOvDlERkTSzKsZM2fMLLQERFEUMvIzScpOJik7hcSsgvfJOckkZiXzxbhP0ZrraDG1vf73yumPDqPJzKPt/J7Ym9liZ2aLvbkd9ma19R/bmdlib2bLM32e5qaSSPOp7VCpVOwd/UuJazc4guXi4mLwF6FarcbZ2bkU3yJRU9y9WbNHUy8y+0wjtfkwfu1jxED3+6+3ioiIYEobjyI9rD45fqkiShdCCFGGHhSO7kdRFFav/ZHn/vVv7G6PPkXfHn0aMfdZ6ratT+LtUJWnzSt0rJWJJQ7mdjiY25Mam4zbcM9Cv1fsfZ24vuESO/61/oEjU/NmzSUoKIizC49g27IOjC756zcYsGbNmsX777+Pv78/jRo10n/90qVLzJw5k1mzZpX8KqJGKLRZ88BgIs5uh0VPs/i79Qx0LzoSeoc2PZW031bR2KYWe68m8mIrd/1fGtLDSgghqh9DU3Pr16+nz8C+JGQnkZiVpA9JSfe8T85O4fC7O7BrWYcWt0eP6vdvxJmPDrP3+x1M6jGNlnW8cTC3v/1mh6O5PfbmdpgZ19LXsdt7EzGn46jfv5H+90rq6QS8vEs27XfnzvS58+cSsSGiVN8Dg1OEgwYNIiwsjPj4eFq0aIGTkxPx8fGcOXOGunXrFpqSU6lU/PJLyYfNyopMEVYtxW3WrFryJL61UjgZVnSzZkWrIePPX0nb8iOoVeyxaMqz7y+gs7tj4R5W0mZBCCEq1MOMPml1WlJyUknITmZw1ydJUt/Sh6M7U3PaTA2t5nYudFztWjaFApKDuT2O5vaM8g/CdXiTQovLY3+/xPUNl8jOyi7x67gT9Gxb1iH1dALJp+Mr5PeKwRGsjIwMPD098fT0BCAvLw9bW1ueeOIJANLT08u1MFH9nI+MQBkYXGizZqV5HyI3Bxd5bs75v0hd/yWahGtYdhyAzYDnGG1pg7lPGz6YN5dPjkfg7eVF6GdfS7gSQogKZGj06ZvVy2nTqz0J2UkkZCWSmJVMQlYiCdkF75OzU9AqOgAuX7xcaHH4nam52NALzO48FUdzB+pYOGBvZouJkUmxdXh7exkcfSqpe0egvLy9WB76VYX8XjE4glUdyAhW1aHRKTh7BZJk4ghv3DOCZZbCyRMFI1jajFRubfyGrOM7MW3cEtugVzB1afSAswshhCip0o4+5Wpyic9KJD4rkZuZCbw88AUyTLKKXRgeMLcrAJYmFjia21PHwpE6Fg76tzvBaXiPIVzXxesXhyuKwtmFR3AzceGvE3+V+HVU1uhTWZCtcsQjy8xX+NcuLSm9psP/PV2wSXPzPqjObivYrDk0FEVRyDq6nVu/fIOCgt2oN7Fo2/uhbn0VQghRPEOjT4u/XYpv9wDisxK4mZlQEKgyE7iZmUhq7q1C57h+JRb3Ykafrm24yPcDv6COhT0WJhb3rWP2zNmFFoffCUfLQ78q8WupzNGnslBoBOvTTz9l9OjR1K1bl08//fT+B6pUTJ48udwLvB8Zwap8yTkKA7ZqOZuisK6nEVnHNjJn3gdERoTTzMubWe+/x8BObUhdu5Tci6ewCOxO7aEvYWRtW9mlCyFElVPa0afM/CziMxNvh6YEXhv8Epkm2QZHn8yNzahrWQcnizo4WTpS9857yzo4WThSx8KRtm3aEKOJe6TRpzuvZe78uUSEF4SjmTNmVptwVBYKBSy1Ws3hw4dp27YtavX9b6lXqVRotdpyL/B+JGBVrrgshT6/a7iRDVv6GRNYR0VoaCjz584hIjKSZp7NmDKsN13TIzCydcRu5H8wayb/nYQQojiGpsQ+W7mEFt39uZER/8/IU1YC8ZmJZORn6o83UqnZN24T7kFeRRaGX9twkYSURKxMLR84c1Ddp+aqigpdg3XhwgWef/55EhMTsbW1ZdWqVfj4+BR5Xm5uLlOmTGHr1q2YmprSqlUrVq9eXeR5ErAqz+U0hd6/a8jTwfYBxnjZqgp1Yu/qZs/eq4nsi07i+6kvMXruYtSmZpVdthBClJvSjj4pikJq7i1uZMQTl3GTsQOeI9040+Dok5WpJc6WTjhZOBYZhapr6Yi9uT2tWwfK6FMVUShgGRkZcejQIdq2bVsuF+vRowfPPfccY8aMYf369XzyySccOnSoyPMmT56MVqtlyZIlqFQq/UbT95KAVTnOJCv0+UODlXFBuHK3LvhrqNhO7L+EkWXvyvGwkv+PLYQQ1Y2hUZ/vf/qBNr3bEZcRz43Mm9zIiOdGZkGgupEZT+5dTTL3j/8NDwOjT4mpSViZWj50HTL6VPEMThGWtfj4eDw9PUlMTMTY2BhFUahXrx6HDx/Gw8ND/7zMzEzq169PbGwsVlZW9z2nBKyKdzReR78tWtwsYWt/Y+pa/DPUbGluzpQ27kwI8NB/7euwKD45fpXMEvYsEUKIyvCwXcezNTlcT79Bv859iu37dPedd+bGZtSzqouzpZP+vbOVE/Us6+Js5UTn9k/I6NNjpMLuIoyJicHFxQVj44JLqlQq3NzciI6OLhSwLl26hIODA/PmzWPHjh2Ym5sTHBxMz549DZ570qRJ2NraMmrUKEaNGlXeL6XGOnhTR78/tLS0V/FbXyNsa/0TrrLPHaWxrZl0YhdCVDuG7rwLCQlh2LBhpObe4lr6Da5n3H67/fG19Bsk56QA9+v7dJGv+31CPau6WJta3Xf908wZMx/5zjsouPuupFvSiPJTJGClpaWRnJxcooPt7e1LdbF7/2EVt/wrPz+fy5cv4+Pjw4IFC/j777/p1asX586do06dOsWed/HixTKCVc4O3CgYuQpwUPFbPyOsTAr+WyqafG5t/paMPaG8OaQXY5etZ/QvYYU7sX/2dSVXL4QQhs2eOwf7lk76kaM7W7JMmPYy3+SvJUvzzwi8ba3auFg7U9/KmQBnX1ysnKlvVY9RXpeJOx1fpCmmj483zRyalKiO6t6WQBRWJGD17du3xAeX5i7CBg0aEBsbi0aj0U8RxsTE4ObmVuh57u7uqNVqRo8u2FHRz8+Phg0bcvbsWbp161bi64mys/+Gjv5btAQ6FoxcWd4OV5qE6yR9/yH5169Qe+gEnu86DJs+G6QTuxCiQpVkei9fm09c5k1i0+KISb/OtfQ4rmfEcT3jJmfCzxQZfbLzdSIm9DzPtXyK+lb1cLF2xsWqrsH+T3PKoO8TyOjT46RIwHrvvfdo3Lhxcc99JE5OTvq7AceMGUNISAgeHh6FpgcBHB0d6dmzJ1u3bmXAgAFcvXqVK1eu0KxZszKvSTzYn3EF4apNHRWb7wpXWX/vJ2XNpxhZ2+I06TNMGzQF5IeDEKJiGZrem/rfd3Ft35CY9OvEpl/nRsZN/TYuZka1bgcmZ55wbcehxttJOZ1QZPSpuU9zRvmU7OeZjD6Je1XYIneAyMhIxowZQ1JSEjY2Nnz33Xc0b94cgAEDBjBnzhxat27N5cuXGTduHElJSRgZGTFr1qxi/5HKIvfytS9Ox4AtWto5qdjU1wgLYxWKVkva79+RvnMt5v6dsfvXZNRm9+/oK4QQ91OaBeaKopCck0ps+nVi067z8qAXim2sqc3MZ8iSUbhau9DApj6u1vX0Hzua2xdasiJ33onyUKEBq6xJwCo/R+J19PpdS9s6/4QrbcYtkr//kNyLp6g9aDxW3YbLVjdCiEdiKNysWfsTgT3bcvVWDDHp14hJu8619OvEpF3Xr4lSoWL/+N9wD2pWpLXB9Q2XyC7F3cty550oaxKwhF5oaCiz584nIjICjVMzmo56hxNzRmBpoiIv5gJJK+ei5Odi//y7mDX1q+xyhRCPgZb+Lbmmi79vewN7M9vbo08uuFrXv/3eBRdrZ9q1aVsmrQ2EKGuF1mB9++235bL+SlR9d/6KVDXvjTIwGM5uJ3LBU2xtE0K/ho4kr1mEST0PHMa+j7Fd8XdzCiFqnpJM7+kUHTczE7h6K4arabFcvRVLdFrB27mIcAPtDS7wZd+FuNnUx9rUcE/EsmptIERZq9CtcsqajGCVHb9WgZzOd0CZuBlUKlAUVEueZH7THP7tYYlFYHfs/jUZlYlpZZcqhKgiDE3vBX81D/cOjbmaFkN02jWi066Rd7tjuZlRLRrY1Me9tivuNg2YM+o9UtRpjzQCJdN7oiqSgCUAMLOwJHdgMPSeBICpLp+PTrzKCIsb2PT7N9Z9n5H1VkIIALQ6LXGZ8XTv0NVg9/IeCwffDlGuuN15b+OKk6UjapVafy5ZYC4eVxXWyV1UXel5CirnZnBoNUaHV8ON8zSwt8YkwIWFOguWLh5d2SUKIcrBg6b3FEUhPiuRK6nRXLl19fb7aK7eiiFXm2ewe/n1DZf4ZcT3JapB2huIx5UErBouX6cwYocWXbPusOVTOro50K29O3uvJvL6H38zderUyi5RCFEODPWPemvJdJzauXIl9SpRt2LIzM8CCvbRa1jbjaZ2Dent0ZWGtu6MNdC93Mu7dNtjSf888TiSKcIaTFEUxuzV8tMlhYaftcE5/Rqrhwbof1CO3hhGloMrx8PkThwhHhcanZbY9Gv07NSTFHVakek9bWY+Iz9/noa13Wlk60ZDW3c8ajegrmWdQlN7INN7QtxPoRGs0NDQUh0sf3FUb+8f1/H9BYU13Y144fUInm7jUWiov5u7A58cj6jkKoUQxSnJ3XvpeRlcTrnKxdQrXEq5wsWUKK7ciiZPm0fUxSsGp/dWDFhcohpkek8IwwoFrBEjRpT4QJVKVaq9CEXV8tU5LfNP6vi4nZqhOUf50MaMvVcTebGVu/6v2T3RSXh7lW6oXwhR/gxN7835aj4u7d25lBrFpZQobmTGA2CiNsajthuN7Tzo3bArjW09eN7rCtdP35TpPSHKSaEpwqtXr5bqYHd39zIvqDRkivDh/BatY/A2La/7qPmg1j5SfvyYnbm1GfPl/+js7kg3Nwf2RCfx59VEGeoXoorR6rS0bOXLDSWx2Lv3en08hMZ2HjS2bUgTOw8a2zXEzaY+xurCS25lek+I8lXo/7jKDkyi/J1JVvjXLi2D3FTMM95NyupPsWjdg+dGTca650g+mDeXT45H4O3lRehnX8sPWiEqkUanJTotlvPJl/RvF1OuEBkZaXB6b0PQdyU6t0zvCVG+ZJF7DZKQrdB2o4baprDTeSfZ65dg2a4vtk9NRKVWP/gEQogyUdz6qcFDBxN1K0YfpCKTL3Ep5Qq5txt0ulq74GnfGE/7RgQ//S6JqhTZHkaIKuy+bRpWr17NV199xfnz58nJySnyeFpaWrkVJspWrlZh+HYt2VrYa7+N7HX/xbLTQGyDXpVwJUQFunf9VPTt9VO+k9pj29oJFSoa2NSnmX1jurl1opl9Y5raN8LSxEJ/jlqzke1hhKjiDAas1atX88ILLzBmzBgOHjzIuHHj0Gq1bNq0CVtbW5577rmKrFM8AkVReHm/lmOJCkfrbUH9yxdYdR5M7eGvSHd2ISpAcnYq4UnnCU88z7vT3sauZR396FP9/o0489Fh0rbF88M7K2hi1xALE/P7nk+m94So+gwGrE8++YT333+f6dOn8/XXX/Pqq68SEBBAeno6ffr0wcrK8Oabomr55LSOVecVtjvvxn7rF1h1GUrtYS9JuBLiITyoPUK2JofzSZcKAlXSBcITz3MzKwEAO7PapMQk4RbkWWj9lN3t9VO+Tj4lrkPu3hOiajMYsC5cuECnTp0wMjLCyMhIPx1obW3NtGnTmDRpEm+++WaFFSoezrZYHVOP6FhufxCvXYux7NBfwpUQD8lQe4R3v5hJ7YA6RCRd4MqtaHSKDjOjWjRzaEJ39054O3ri7eCJk4UjrXxOEXM67pHbIwghqjaDAat27drk5uYCUL9+fc6dO0e3bt0A0Gq1JCUlVUiB4uFFpSuM2qXlLbPj9D34MeYBXbEd+bqEKyEeUvDc2di1dCo0vXf6o8MsXbSEEf99Dm8HT4Y3exJvB0/cazfAWG1U5BwzZ8yU9VNC1AAGA1br1q05deoUffv2ZfDgwcyePRudToeJiQkLFiygXbt2FVmnKKUcTcEeg12zTzPxwoeYebfB/pkpqIr5gS+EKF5CVhKnE8I5mxDBmcQIzoafNdge4dsnl5bonLJ+SoiawWDAeuedd/SNR+fMmcPVq1eZPHkyWq2WNm3a8PXXX1dYkaJ0FEXhtQNalOsXWXx1LrUaNcfh+XdQGcne3qJmu9/6Ka1Oy6XUKE4nhHMmIYKzCRH6tVMuVs40d/TCrbE7qacTpPu5EOKBStUHKzc3l9zcXGxsbMqzphKTPljF+zpcx5xdcey99hZWDo7UeX0h6lr3vytJiMedoc7lkz57CxM/a84khJOZn4WJ2hhP+8a0qONNizpeNHf0wsHc7r7nkO7nQoh7SaPRx0RoaCiz584nIjICCxcvtg5pSgNbc5wmfYqRtV1llydEpfNr5UeM9kaR7WV0mRpeXP4f/Jxa4OfUnGYOTahlZGrwPKGhoQXTe+EF03szZ8yUcCWEKOK+c0bbtm1j/fr1xMbGFmk0qlKp2LlzZ7kWJ0rmzl/VRubWKHn52KdcYP/pPLyem8wACVeihsrKz+ZsYiR/3zzD3/FnOXOf9VOLeswu8Xllek8IURIGA9bHH3/MtGnTaNCgAV5eXlhaWlZkXaIU3pzyFgAd6tSiu7sLe68mMmX7Gdwvzifq2XGVXJ0QZeNB/adytXmcSQjnxI1T/HXjFJHJl9AqWmrXssHPqTlujdzKZP2UEEKUhMEpQg8PDwYOHMh///vfKntbv0wRFjA1NqJdfTvWDA3U/+IYteEER6+nkqfRVHZ5QjwyQ2ufFn+7FPvWzoTdOMWZhHDydPnY1qpNgHNL/J1a4Fe3Be42rqhUKlk/JYSoUAZHsJKTkxk+fHiVDVfiHzoFurs7Fpr66OHhyOFrKZVcmRBlY868OdgX039qxpz3eeLDfvjVbcGL/v8m0NmPhrZuqFVF99eU9ghCiIpkMGANGjSI/fv306NHj4qsR5RSRr6ClWPBtOCLrdz1I1i7ohJo6OFe2eUJ8UjiMxM4FneSs+fOFdlext7XiWsbLrJp5GqM1SVrQSLrp4QQFcXgT6UxY8bw6quvkp2dTe/evbG1tS3ynJo8LVdVvLPjJjPae/D2r/t5ZsMJuns4sjsqkQMxyYSGLq/s8oQolVxtHqfiz3L0+l8cjQsj6lYMapWa2q52pBSzfsrb27vE4UoIISqSwTVYanXhIfa7pwoVRUGlUqHVasu3ugeo6Wuw1kdmY/7t2zQzSuWkd28WfLqY8IgIvL28eO99uXVcVC3FLVIfNmwY0WnXOBb3F0euh/F3/BlytXnUMXegrUsr2tQLoLWzH9t/2ybrp4QQ1YrBP/12795dkXWIUrqWqRC39gsG5kXjMukTmrg1ZcSzz1d2WUIU695NkqNvb5LccWpPjH0tMVEb4+fUgvF+o2lbLwCP2g0K/VEn66eEENWNNBqthhRFYd43fzD23FJMRr5J3U59KrskIQxSFIXm/i24oSQWafJpqTFj/Y4N+NdtgZlxrcouVQghyswDFy/s27ePP//8k+TkZOzt7enSpQudO3euiNqEAT/uu8Co8P8juWU/fCVciSooX5vPqYRzHLp2jIOxx4iMjDTY5LN9/cBKrlYIIcqewYCVmZnJsGHD2LFjB8bGxjg4OJCUlIRWq6VXr15s2LABCwuLiqxVAJFx6TT67QNSarvR+rlXK7scIfRSc9I4fP04h64d51jcX2TmZ1HH3IEOrm041rQxydLkUwhRgxgMWNOmTePIkSOsWbOGkSNHYmRkhFarZf369bz00ktMnz6dpUuXVmStNV6eRsvZ5Yvw1qZT/5UPUZkY3i9NiLJkaIH61bQY9scc5dC1Y5xNjERBwduhKU97D6Vj/TY0sWuISqWi4ey6BAUFcXbhkUKL1JeHflXZL00IIcqFwTVY9erVY+bMmbzyyitFHvvyyy+ZO3cucXFx5V7g/dS0NVihq9fT9vhykkbOwq9Th8ouR9QQhjqgd57eB1ULM8yNzQh09qOja1vauwTiYF78/peySbIQoiYxOIKVmppKo0aNin2scePGpKamlldNohjhZy/gf2IVJ72GMVDClahAs+fOwa6YLuqXNpxj3WuhBNbzo5bRg0dTpcmnEKImKbqfxG3e3t589913xT723Xff4ePjU25FicLysrNJW72AqxZu9BwzprLLETVArjaPg7HHWHBoKWfDz2DXsk6RBerJ0Yl0dG1TonAlhBA1jcERrJkzZxIUFERUVBQjR47E2dmZmzdvsnbtWo4ePUpISEhF1lmjHfr2K+rlJGL04n8xN5Nb2UX5yNXkcvj6CfZEH+DQteNka3Jws6lPvYb1SZUF6kIIUSoGA9bQoUPZsGEDs2fP5q233tJ3b/f392fDhg0MGjSoIuussS4d2Efj81vZ1mYi43zcKrsc8ZjJ0+ZzLO4vdl39kwOxR8nW5NDUrhGjm4+gc4P2eNRuQGh+qCxQF0KIUrpvH6zBgwczePBgMjMzSU1NxdbWFktLy4qqrcbLT45Hs2EJh+w78fRT/Sq7HFFN3XsH4Hvvvodbx8bsuvon+2OOkJGfScPa7jzjE0R39040sKlf6Hjpoi6EEKUnndyrKEVR+GvRu2hvRpP96pd0aVS7sksS1dC9dwCmnI4n5XQCPm+0plWP1nR3f4Ie7p1paCujo0IIUZYKjWBNnDiRt956Czc3NyZOnHjfA1UqFUuWLCnX4mqya3t+x+naX6zpNJupEq7EQ1AUhfeCZxS5A/DMR4dRdmfxw2dfFtrvTwghRNkpFLA2bdrE+PHjcXNzY9OmTfc9UAJW+dEk3SBn83L+cOzDhMFtK7scUc3czExg+5W9bLuym/PnzxfZosbO14krGy5JuBJCiHJUKGBduXKl2I9FxVF0Oi6u+pRUtRV1h0/Atpb8EhQPlpGXyd7og2yL2sPJm2eoZWRK5wbtady0MUlyB6AQQlS4B272LCpW6p+bsYo5xZcB81jiLTcUCMM0Og1Hr//Ftiu7OXDtGPnafAKcfXmnwxt0adAeCxMLms9uKHcACiFEJTDYaPTEiRPs3LlT/3lKSgovvvgiTzzxBMHBweh0ugopsCYIDQ0lsJU/lubmdHnqOabfcOH1wYEyhSMIDQ3FP8Afcwtz/AP8CQ0NJTotlq/++o6RG8bzzt55XE2LZZzvKNYNW8GnPefQr1EPLEwKNmK/cwegm4kL1zdcws3EhdDQULkDUAghypnBuwi7dOlCz549mTVrFgDPP/88GzdupHfv3mzZsoVp06bx/vvvV2ix93oc7iK8c5dXZzcHurk7svdqIvuikwgJCZFtRWq4oncAJpByOh6fN1rTsJMnvRt2ZUCjXjS1L35LKyGEEJXHYMBydHTkhx9+oH///mRnZ+Po6Mjnn3/O2LFj+eKLL1iyZAnnz5+v6HoLeRwCVpNGDamnSWfNsED9GplRG45zw7g2Fy9fruzyRCXyb+VPjDZOfwegoiic+egwDthy5uQZ2aJGCCGqMINThFlZWVhYFEwzHDhwgNzcXIYMGQKAr68vsbGxpbrQhQsX6NixI56enrRt25Zz587d9/mzZ89GpVJx5syZUl2nuom6Gk13D8dCd3n18KhD1NWrlVyZqCypOWn8HP4LZ8LPYnvPHoB2vk7EXoqRcCWEEFWcwYDVqFEj/vjjDwB+/PFHAgMDsbe3ByA+Ph4bG5tSXeill15iwoQJnD9/nqlTpzJ+/HiDzw0LC+Pw4cO4uT3+zQ/VKthzNZE7A4mKorArKhG1LL+qURRF4UxCOPMOfMbIDeP45uT31HF3IuV0QqF/G3IHoBBCVA8G7yJ88803eeGFF1ixYgXJycn88MMP+sf27NmDr69viS8SHx9PWFgY27ZtAyAoKIjXX3+dqKgoPDw8Cj03NzeX1157jTVr1tC9e/dSvpzqp32zxvx57gLPbAyju7sDu64mcSAmCY9Gsq6mJsjKz2L7lb38cmELl1KjcLFyZpzfM/Rv1ItdtXbIHYBCCFFNGQxY48aNo0mTJhw7doyAgIBCYcfBwYE33nijxBeJiYnBxcUFY+OCy6lUKtzc3IiOji4SsGbOnMmzzz5Lw4YNS3z+SZMmYWtry6hRoxg1alSJj6tsiiaPrwe1YmVthU9P3uTgtcuoTM0A+HTRokquTpSnSylR/HJhC9uv7CFHm0uH+q15udXztK7nj1pVMLAsewAKIUT1dd8+WF26dKFLly5Fvh4cHFzqC93bcqC4tfWHDh3i2LFjLFiwoFTnXrx4cbVc5J6+OxTz3Az2tH0Ln5xvuBARTjMvT2a9/578En0M3LvJ8rvvvot9a2c2nv+D0wnnsDezY4TXIAY26UNdyzrFnmP48OFyN6kQQlRD9w1YWq2WI0eOEBsbS05OTqHHVCoV//73v0t0kQYNGhAbG4tGo8HY2BhFUYiJiSmyxmrv3r1EREToR69iY2Pp27cvy5cvp3///qV5XVWeJiWB1G0/scJ+EG+MepHRi1+q7JJEGbq7xYLLsMZEn47j6ZFP4/NGa3oN7EPwE1Pp3KAdxmrp9SuEEI8jg20awsLCGD58ODExMcWONqlUKrRabYkv1K1bN8aMGcOYMWNYv349ixYt4vDhw/c9xsPDg82bN9OiRYtiH6/ObRoSV33A9bOneLPD12wfZi1NRR8z/gH+xGiKtlhwVtfh3N9nK7s8IYQQ5czgXYSvvPIK1tbWbNu2jWvXrpGQkFDoLT4+vlQXWrZsGcuWLcPT05MFCxawYsUK/WMDBgzg+PHjD/8qqpmcC3+Tc3If85zG8FFnCVePE52i49C145w9V3yLhSsXpLeZEELUBAbnJ86ePcvatWvp2bNnmVyoWbNmHDp0qNjHfv/992K/HhUVVSbXrkoUrZbkkP/jlGUzrNr0ILCOhKvHQbYmhy2XdrI+cjOx6dep7WpPimyyLIQQNZbBgOXp6Ul6enpF1lIjZB78De2Nq8xq8gkbWptUdjniEaXmpLHh/G9sOP8bGXmZdGnQgXc6vEGk2VlGjBghLRaEEKKGMhiwPvvsM9544w38/Pzw8pK/usuCNuMWqb99zzr73vRu1wxXKxm9qq7iMm6yNuIXfru4HYAnm/TmKa8h1LOqC0CLIC9psSBEJYmOjiYxMbGyyxDVmKOj4yM3OzcYsF5//XVu3LhBixYtcHFxwdbWttDjKpWKv//++5EuXtOkb1tDrkbLV67/5pifweVvogq7mHKFn86FsvvqfqxMLXnGJ4ihngOwNSu6s4G0WBCi4kVHR+Pt7U1WVlZllyKqMQsLC8LDwx8pZBkMWIGBgbL4ugzlJ1wjff9mPqszmv+0tae2qXxvq7JCPayaNePZ/zzPTY9bHI37C2dLJ14LHM+Axr0wNzar7FKFEHdJTEwkKyuL1atX4+3tXdnliGooPDycZ599lsTExPIJWKtWrXrok4qi0jZ/S4qpPbvcB/Opt4xeVWVFe1hd5+3xb9L93QHMeGEq3d07Sf8qIao4b2/vate+RzxeSvSbXlEUrl+/jkajKe96Hku5V86R/fd+5jg8y+z2FpgayehVVTZn7hzsWjrRfGo7XAc0psXU9ti1cCLlj+v0bthVwpUQQogHum/A2rp1K+3bt8fMzIwGDRpw6tQpACZMmMCPP/5YIQVWd4qicOuXb7hs1YhrTbszoqGEq6pKURQOxB7lTPhZ7Ir0sKpDREREJVcohBCiujAYsH766ScGDBiAu7s7S5cuLdTNvXHjxnz77bcVUmB1l/33fvKiwnnXYRwL2hnLurYqSKfo2BdziBf/eJN3987HroEDKacT9P/mpYeVEEKI0jI41zF37lwmTZrEJ598glar5ZVXXtE/1rx5cz777LMKKbA6UzT53Nq0kiN2rbFo5k/XerL2qipRFIXD10+w/O/VXEy5gn/dFizuNY/LZuelh5UQQohHYjBgXb58mQEDBhT7mKWlJbdu3Sq3oh4XmUe2oUm6wYzG7/JdawlXVcnJm2f45uQPnEmMwLeOD0t7fYBf3eYAtApqKT2shBBCPBKDAcvZ2ZmIiIhit8o5deoU7u7u5VpYdafk5ZK2dQ07HLrg2awhbZ0kYFUF4UkXWH5yNcdvnMTTvjEfd59Fm3qtikzdSg8rIYQQj8Lgb/1nnnmG4OBgdu7cqf+aSqXizJkzLFy4kGeffbZCCqyuMg5sRpOewhy7Z5gbaFTZ5dQooaGh+Af4Y25hjn+AP6GhoVxJjWbGvg95ectbJGQlMqfzNL7u9wltXQJkXZwQolKkp6djZWXFCy+88MjnCg4O5q233nrg81atWsX58+f1n//666+8/fbbD33dVatWMWLEiIc+/nFmMGAFBwfTsWNHevfujbOzMwD9+/fHz8+P1q1bM3369AorsrrR5WaTtuNnNjn2or1PfXwd5Bd4RbnTwypGE4fLsMbEaOIICgpiYPAILiRf5p0Ob/Dtk0vp6tZRgpUQolL973//IyAggJCQEDIyMirkmvcGrMGDB/Pxxx9XyLVrGoMBy9TUlF9++YWdO3fy4osv8sILLzBmzBi2bt3KL7/8grGx9AIyJGPvRjTZWSywe5rgABm9qkhz5s3B/q4eVs2ntsO2RR2ytyexetCX9GvUAyO1/DcRQlS+FStWMG3aNDp37szatWuBggDUt29fRo0aRcuWLWndujWXL18G4MaNG3Tv3p3AwECaN2/OxIkTC93hf0fLli05dOiQ/vNly5bx9NNPs3z5co4fP87EiRPx9/fn999/LzIC9e233+Lv768fTImKikKj0dC3b19at25N8+bNGT16tGxFVAIPTEndu3ene/fuFVHLY0GXlU76rvWsdexP7+Z1aWYroyQVKTIiEpdhjQv1sLL3deL6hkuYGJlUcnVCiKoiLi6OuLi4Mj9vvXr1qFev3gOfd/bsWWJiYujXrx8ajYaFCxcybtw4AI4cOcLff/+Nu7s706dP56OPPmLZsmXY2tqyadMmrKys0Gq1DBkyhJCQkCJTdBMnTuSLL76gQ4cOAHzxxRd88cUXdO7cmdWrV/PWW28xcOBAoPCuLXv27GH+/Pn8+eef1KtXTx+ijIyMWLNmDQ4ODiiKwquvvsqXX35ZoinJmuyBASstLY3Y2FhycnKKPCbbEBSVvjuEfI2GT+1HckhGryqMoij8GXsY83pWpJyOp37/RqhUKulhJYQo1rJly5g9e3aZn3fWrFkEBwc/8HkrVqzgueeew8jIiCeffJKXX36Z8PBwAJ544gn9jWQdOnTgv//9LwA6nY5p06axf/9+FEUhPj4ef3//IgHr2WefZdasWcTHxxMeHo5KpaJz584PrOm3337jueee0wdECwsL/XU/++wzfvvtNzQaDbdu3aJLly4l/p7UVAYD1rVr1xg/fjzbt28v8piiKKhUKrRabbkWV91oM26RvncjPzgOYqCPPR7WMnpVEcITz/NF2LecTjhH5+d68GvwWulhJYS4r5deeonBgweX+XlLMnqVn5/P6tWrMTEx4aeffgIgKyuLlStX0rx5c8zM/tlE3sjISL9N3aeffkpSUhJHjhzBzMyMN998s9jBD3Nzc55//nmWL1/OX3/9xeuvv/5Ir2nNmjXs3buXffv2YW1tzdKlS9m3b98jnbMmMBiwnnvuOS5cuMDSpUvx9PTE1NS0IuuqljL2bkCjg6W2QzniL6NX5S0pO4WvT37Plsu7aGzrwaIewbQZ3YrQlk9LDyshxH2VdCqvPPzyyy80atSIw4cP67925swZevbsybx58wwel5KSgrOzM2ZmZty8eZN169bx9NNPF/vc1157jQ4dOpCbm8t3332n/7qNjY3BPpaDBg1i3LhxTJgwAWdnZ/0UYUpKCg4ODlhbW5Oens6qVato1KjRw7z0GsVgwDp69CirV69myJAhFVlPtaXLSidj3y/85DiAJ71saWQjo1flJV+bT0jkZr47/TMmRsZMafsqTzbupV+8Lj2shBBV2YoVKxg9enShr7Vo0QIXFxfS09MNHjdx4kRGjhyJv78/9evXp1evXgaf6+rqir+/P56envqpPijYS3jKlCl8/PHHfPDBB4WO6dKlCzNmzKBPnz6oVCpMTU1Zv349zz33HL/88gs+Pj7Ur1+fzp07c+3atYd89TWHSinuFgSgVatWvPfee1W6v0VYWBiBgYGcOHGi0taDhYaGMn/uHCLCz9GotgXRQz/j2Edj8JTF7eXiyPUw/nv8G65n3GCIZ3/GtXwG61pWlV2WEKKKqAq/F6qCjIwMvLy8+PPPP2nYsGFll1OtlNW/IYNtGhYtWsS8efMK9csQhd3puWSVco0pbRthZ2ZC2tfjOLNrQ2WX9tiJTY/jnT3zmLp7No4WDnzT/zPeaD1BwpUQQtzjq6++wsvLi1dffVXCVSUyOEXYs2dPevTogY+PDy4uLtja2hZ6XKVS8ffff5d3fVXa/Llz6OLmyOohBd3AX2zlzuiNYXwwb65MUT2k0NBQ5sybQ2REJM28mjH9nXfIaJrP/8I3YG9mx+zOU+naQJqECiGEIS+//DIvv/xyZZdR4xkMWG+//TaLFy8mMDBQFrkbcO5cOG+3a1io51I3dwcWHjlXyZVVT3dGBO1bOuEyrDHRp+MY9dS/aDmpHZPGTuQZnyDMjGtVdplCCCHEAxkMWMuXL2fOnDnMmDGjIuupVoyMjdlzNZEXW7nrey7tikrERLrcP5S7u7CrVCrq92/EmY8Oo92dwbjPnqns8oQQQogSM5gETExMaNeuXUXWUu3odFr+jE5i9MYTdHN3ZNfVJA7EJGEio30Ppbgu7Ha+TlzecKmSKxNCCCFKx+Ai9/Hjx7NmzZqKrKXa+U/3QL4a4MdxXV3mH7nKYeqDqy8+LVpWdmnVTmTSRSxdrEk5Ha/fW0u6sAshKptOpyMjIwOdTlfZpYhqxuAIlp2dHWvWrKFTp0707Nmz2EXukydPLu/6qixFp+UFP1d2nYok29oFOoxBdXYbnNvBrKWhlV1etZGjyWXF3z+yPnITvk+3YfcHv0sXdiFEpTt58iSfffYZ69atIzs7G3Nzc0aOHMnkyZPx9/ev7PJEdaAYoFKp7vumVqsNHVphTpw4oQDKiRMnKvzaWX/vV2Le6Kv4v/JfpWGLQMXM3ELxaxWohIaGVngt1VXYjVPKqI0TlF4/jVDWnA1R8rUaJSQkRPEP8FfMzM0U/wB/+X4KIUqlLH4v/Pjjj4qxsbECFHkzNjZW1qxZ88h15uXlKcHBwUqzZs0UHx8fxd/fXxkyZIjy119/KT///LPi5+en5OXlKYqiKBqNRmnXrp3y7bffKoqiKM8//7xSv359xc/PT2nZsqXSuXNnJTw8/JFrEgXKKlsYHMGS4VDDFEUhfec6Iu1aoGvxEpe+eE3aBpRCZn4WX/31Hb9e2IJvHR8+6j6TBjb1AenCLoSoXCdPnuT555/X7/93L41Gw3PPPYe3t/cjjWSNHTuWjIwMDh06hJ2dHQCbNm3i7NmzjB49mvXr1zN79mzmzZvHggULqFOnDmPGjNEfP336dP0egwsXLuT9999n3bp1D12PKHtyu9tDyLt8lryrEXzgNpOpfkYSrkrh8LUTfHL0S9LzMpjUegJDPPujVhlcCiiEEBXqs88+Mxiu7tBoNCxevJhVq1Y91DUuXLjAhg0biImJ0YcrKNgL8I4vv/wSPz8/3N3d+fzzz/nrr7+KPZeiKKSmphY6z9atW5k7dy7Z2dkYGxvz8ccf06VLF27cuMGoUaNIS0sjJyeHnj17smTJElQqFcHBwURGRpKenk5ERAQBAQG88847TJkyhaioKIYOHcqnn376UK+3pioUsMLCwvD29sbc3JywsLAHHlxTtyFI37WOm9ZuXHYO5OlGEq5KIi03nS/CVrLl8i5aO/vzVrtXqWdVt7LLEkIIPZ1OV+JRoLVr17Jy5UrU6tL/gfjXX3/RpEkT7O3tDT7H0dGRjz76iH//+9+sXr0aZ2fnQo8vWLCA5cuXk5CQgJGREfv27QPg8uXLzJ49my1btmBjY8PFixfp2rUrUVFR2NrasmnTJqysrNBqtQwZMoSQkBD9lnjHjx/n+PHjWFlZERAQwPTp0/njjz/QaDQ0bNiQl19+GU9Pz1K/3pqqUMBq3bo1hw8fpm3btrRu3drgyIyiKKhUKrRabYUUWZXk37hKztkjfFz/Dd70NcZYLQHrQQ5fO85Hhz8nT5vH1Pb/YUCjnjLqJ4SocrKyssjOzi7Rc7Ozs8nOzsbS0vKhrnX3z8BLly4RFBREdnY2Xbp04ZtvvgFgw4YNuLq68tdffxXZHPruKcKVK1cyYsQIjh8/zpYtW7h48SJdunQp9PyYmBicnZ2ZNm0a+/fvR1EU4uPj8ff31wesvn37Urt2bQB8fX3x8/OjVq1a1KpVi2bNmnH58mUJWKVQKGDt3r0bHx8f/ceiqPRdIaSbObCnTldWNJOprXvdvdVN02aetHq6HVENbtK2XgBT279OHQuHyi5RCCGKZWFhgbm5eYlClrm5Oebm5g91nVatWnHhwgVSUlKws7OjcePGnDx5klWrVrF582YAfvzxR65cuaLfeHj48OF07Nix2PP961//Yvz48SQkJKAoCv369eP7778v8rx58+aRlJTEkSNHMDMz48033yQnJ0f/uJmZmf5jIyOjIp8/aOpUFFYoIXTt2hUrKyv9xw96q2m0t5LIOrGLL+wG80rLWliayCjM3e5sdROjicNlWGNitXF8/843dEzxZWH3mRKuhBBVmlqtZuTIkSV67lNPPfVQ04MATZs2ZciQIYwfP57U1FT91zMzMwGIi4vjrbfe4rvvvqNOnTosW7aMcePGGQx+O3fuxMHBAQcHB/r06cOWLVs4c+aM/vGjR48CkJKSgrOzM2ZmZty8eVMWxZczGYIphYwDm8lXmbDWvi+vNZdv3b3u3urGdUBjWkxtj31LJ7Z8u1mmBIUQ1cLkyZMxfsB2Z8bGxkyaNOmRrrNq1SpatmxJu3bt8PHxoVOnTuzYsYO3336bF198kf/85z+0bFnQtLp///507NiR9957T3/8ggUL8Pf3x8/Pj7lz5xISEoJaraZp06asXr2aF154AT8/P7y9vVmyZAkAEydO5ODBg/j7+zNu3Dh69er1SK9B3J9KUW63zQYaNmxYql+Ely9fLpeiSurO0OmJEyfKfcG9kp/H9eBn+cmyCxc7v8JXnY3K9XrVkZmFGfWHNcF1QGP912J/v8T1DZfIzirZugYhhHgUZfF7Yc2aNQZbNRgbG/P9998zatSoRy1VVFFllS0KxfQnn3yyUMDauHEjqamp9OjRg7p163Lz5k127dqFnZ0dQ4cOfeiLVkdZJ3ajy0zni3qD+KOljF7dTVEU1kduwrSuOSmn46nfv5F+82vZ6kYIUd0888wz+Pj4sHjxYtauXavv5P7UU08xadIk6eQuSqRQwPr888/1Hy9atAhXV1dOnz5daJuclJQUBgwYgKura4UVWdkURSF930aO2Qfi06Q+XrYy3XVHSk4qCw4t5fD1Ezz16ihWTP0/2epGCFHt+fv7s2rVKlauXKkPWA+75krUTAb/tSxZsoR33323yB6EdnZ2vPPOOyxdurS8a6syci+eQnP9Cp9ZD2GyjF7pHYs7ybjfJhGRdJGPus1k+dtfEhISgpuJC9c3XMLNxIXQ0FCGDRtW2aUKIcRDUavVWFpaSrgSpWZwJV9ycjK3bt0q9rFbt26RkpJSbkVVNRn7fuG6lRu3GvjR00VGrzQ6Dcv//pGfzoXS2tmfdztOwsG8oIuwbHUjhBBC3Cdg9ezZk2nTptGgQYNCLRn27NnD9OnT6dmzZ4UUWNk0iXFknznEUudXmdTSuMbfDXctPY45Bz7hQvJlXm71PE97D5WtboQQQoh7GAxYy5YtY/DgwfTo0YPatWtTp04dEhISuHXrFq1ateKrr2rGupqMP38lx8SKP527sbxxzQ5X26/s5dOj/4edWW2+6PsR3g5NK7skIYQQokoyGLDq1avHsWPH2LJlC0ePHiUuLo569erRtm1b+vXrV5E1VhpdThYZh7fyve0AxrawwMy4ZgasXG0en59Ywa8XttDLowtvtn0FSxOLyi5LCCGEqLLu300N6NevX40JVPfKOrYDXV4OP9gN4Ih3zZwGu55+g1n7FxKVGs1b7V5jYOPeNX6aVAghhHiQB6aG69evc/ToUfbt21fk7XGmKAoZ+zezx7Y9PXycqGtR80LFgdgjvPjHm2TkZfJF348Y1KSPhCshhChDK1euxN/fHw8Pj0Ld4detW8crr7yi/3z48OEcOnSoROdct24dY8eOZejQoXh6euLv70+/fv2IiorSP0dRFDp37syVK1fK6qWIexgcwbp8+TL//ve/OXz4MFDwH+NuKpUKrVZbvtVVorxLp9HcjGaZ+wS+bFEzurbfvVGzo7sTNn2dGDpsGNM7TMTa1KqyyxNCiMfKxo0b+e233zhy5AgAzZs356WXXsLb25uNGzfy3HPPAQV7CaamptKhQ4cSn/epp57CxMSE/v37o1Kp+Pzzz5kwYQLbtm0DCn6HT548mdmzZ7Nq1apyeX01ncGA9eKLLxIdHc3XX3+Nj48PpqamFVlXpcs48BvXLVxRN/LDz+HxH7W5s1GzfUsnXIY1JuV0PLFLjjOn83QJV0IIUcYURWHatGkcOHCAWrVqAeDu7k54eDhNmjThwIEDfPvtt0DBTWejR48GCpp9t2zZkpycHFxdXcnNzeXy5cuMGTOGZcuWkZ+frz/27t/b7du3Z/HixYVqGDRoEC+//DLp6elYW1tXzAuvQQwGrKNHj/Ldd9/VyJ5G2rRksv7ezzKncbxeQ0av5sybg93tjZpVKhX1+zfi7MIjzPtgHkFBQZVdnhBClKm4LIW4rLI/bz0LqFeCJSUHDx4kLi6u0IbL4eHhBAcHs3v3bjp27KgPSHv27OGtt94CCpp9P/PMM9ja2vLuu++ydetWPvroI5YtWwZQ5Ng7li5dyqBBgwp9zcTEhBYtWnDgwIEau9a6PBkMWPXr18fIqGzDxYULF3j++edJTEzE1taWVatW4ePjo388JyeHf/3rX5w7dw4LCwucnZ356quv8PDwKNM6HiTz8FY0KmP+rNeD//N4/EevFEUhPDwc1+FN9WusVCoVti3rELEhopKrE0KIsrcsXMfsMF2Zn3dWgJrgwAf/7jx69Cjjx4/ns88+AyAmJoamTZvi7+/PtGnTCu2AERsbi7Ozs/7zkydPMnnyZABOnDhRaG/EjRs3Ftk944MPPuDChQvFtldydnYmNja2VK9RlIzBgDVv3jwWLFhA586dsbe3L5OLvfTSS0yYMIExY8awfv16xo8fX2TR3oQJEwzOGVcERacl/eDvbKzdhWda2GBq9HgHrFxtHp8e/T9Mnc1JOZ0gGzULIWqEl7zVDHYv+7vD65Wwg01iYiIWFv88ee3atfTr1w8rKyu2bt3KwoUL9Y9ZWFiQnZ2NnV3BjhknT57Uh6oTJ04wZMgQoOCP5XuPXbRoEaGhoezYsaPQ9e7IycnB3Ny8tC9TlIDBgPXdd98RGxuLh4cH/v7+RfYkVKlU/PLLLyW+UHx8PGFhYfqwFBQUxOuvv05UVJR+hMrMzIwBAwbojyluzri85Zw7ipKawPeNB/Cb1+PdmiE+M4EZ+xZw5VY0b059kzmvzpKNmoUQNUI9C1WJw1B58PLyYsWKFUDB7M7SpUvZsWMHR48exdvbGyurf9a++vr6EhERgYuLC9euXUOtVlOvXj0ATp06xYwZMwCKHPvpp5/y008/sWPHjiK/w+8IDw/Hz8+vHF9pzWUwQWRkZNCkSRMCAwMxMjIiPT290FtaWlqpLhQTE4OLiwvGxgWZTqVS4ebmRnR0tMFjipszLs6kSZMYPHgwP/30U6lqKk7G/s2EW3nSzKcpLpaP7+jVyZtnePGPKaTm3OKLPguY/cpM2ahZCCEqyIgRI7CwsMDDw4NRo0bx/fff07RpUzZs2MDQoUOLPPePP/4A4K+//io0JWhjY8OXX34JUOjY2NhYpkyZQmpqKt27d8ff35927doVOu+dtg0tWrQol9dY0xkcwdq9e3eZX+zeHkr3tn642/3mjO+1ePFiAgICHrk+TeJ1ciNO8I3LG/zH5/EcvVIUhdDzv/HFiZX4OvkQ/MTb2JrVBmSjZiGEqCjm5ub8/vvvRb6+adOmIr9/x4wZQ4cOHQgODmbgwIEMHDhQ/9iJEyeKPdbV1fW+v2MBvvrqK/3ieVH2HtjJHW433czIwMrK6qEbTTZo0IDY2Fg0Gg3GxsYoikJMTAxubm5FnvugOePyknHgNzJNrLji0ZknnB+/0as8bT6fHP2SLZd3MdJrMC+3GoOxumbcJSmEENXB2bNni3zN2tqaxYsXc+XKlfuONhV37P24uLgwduzYUtcoSua+wzR79+6lR48emJubY2tri7m5OT179uTPP/8s9YWcnJxo1aoVq1evBiAkJAQPD48idwjemTPevn27wTnj8qBo8kg/sp01Nj2Z4Gvx2HUsT8lJ5c2d77Mr6k/e6ziZ1wPHS7gSQohqolevXmU+lTdx4kTU6sdztqYqMDiCtX37dgYMGICnpyfvvPMOzs7OxMXFsX79enr27Mnvv/9eqH9HSSxbtowxY8bwwQcfYGNjw3fffad/bMCAAbz77rtMmTKFRo0a0b17dwBq1aql73JbnrJPHUSVlcYm97782fjxCleXUqJ4d+988rR5LOk9Hx/HZpVdkhBCCPFYMxiwZsyYwYABA9i4cWOh0ZxZs2YxdOhQZsyYUeqA1axZM4N7Kd2Zi37QnHF5ST+0hTBLH7r5umNp8vgErIOxx5hzYBH1revxQdcPqGtZp7JLEkIIIR57BscGT58+zSuvvFJkqkylUvHKK69w6tSpci+uomgSr5N/4STf2/bllcdkcbuiKPwcvpF3984n0NmP//b+UMKVEEIIUUEMjmBZWVlx7dq1Yh+LjY0t1KOjugoNDWX+3DlEnDuHh60lpuMSaVq7+o9e5Wvz+fTYV/x+aQejm4/gBb/RqFWPR3AUQgghqgODv3UHDx7M9OnT2bp1a6Gvb9u2jffee0/fOba6urO5sVXKNaa0a4SDmTEnFzxLaGhoZZf2UEJDQ/EP8MfcwhwXrwb8+POPvNthEhP8/y3hSgghhKhgBn/zfvzxxzRq1Ij+/ftja2tLs2bNsLW1pX///jRs2JCPP/64Iussc/PnzqGLmyOrhwQwIcCDH4cG0NnNkQ/mza3s0krtTliM0cThMqwxWnMdpxcfIfNkSmWXJoQQQtRIBqcI7ezsOHToEJs3b2b//v2kpKRgb2/PE088wZNPPlntb+08dy6ct9s1LLS5cTd3BxYeOVfJlZXenHlzsGvpRPOp7VCpVNTv34izC48wd/5caRwqhBBCVIL7NhpVq9UMHjyYwYMHV1Q9FcbI2Jg9VxN5sZW7fnPjXVGJmBiXqPdqlRIeHo7r8KaFwqJtyzpEbIio5MqEEEKImql6D0M9gjyNhj+jkxi98QRfh0UxamMYB2KSyNVoKru0ElMUhR/OrMXU2ZyU0/H6FheKopB6OgEvb69KrlAIIURx+vfvz+eff17k635+fmzYsKHU5wsODiYvL0//ub+/P9nZ2Q9dn0qlIiMjo9jHhg8fbrDl0v2sW7eOV155BYCDBw/SqFGjQo8rikLnzp25cuVK6QuuggoFLGtra2xsbEr0Vrt27cqquUz4+LRgfm8/ovJrMf/IVQ5TH1x98WnRsrJLKxGNTsPCI5+z/O8fGTNxPCmnEzi78Aixv1/i7MIjJJ+OZ+aMmZVdphBCiGKMHz+eb7/9ttDXjh8/zo0bNwrtNfggmtuDArNnzy4UsE6ePIm5uXnZFHuXo0ePkpqaSocOHUp97MaNG/WbUYeFhdGqVatCj6tUKiZPnszs2bPLotRKV2g+bMqUKQ/cImb//v3s3Lmz2m8ls+j15/E6/Tuh6XWIeaIPqrPb4NwOZi2t+ncRZuRlMuvPhZyMP8O7HSbRt1F3urt3Yu78uURsiMDL24vloV8xbNiwyi5VCCGqJO2tJLRpyWV+XiMbe4xqOzzweYMHD+bVV1/l77//xs/PD4CVK1fy3HPPYWJiwrFjx5g2bRppaWnodDree+89goKCgIIgsmjRIjZt2kSbNm1IT08HoGPHjqjVarZt20bdunVJT0/HysqKQ4cOMXXqVNLS0lAUhblz5zJkyBCeffZZIiIiyMvLw83NjZUrV+Lk5HTfupctW8bo0aP1n+fm5mJvb09mZiYA7du3p0mTJqxevZpLly4xcOBAzp49i1ar5cCBA/pQeeLECQICAoqcf9CgQbz88sukp6djbW1dgu941VUoYAUHBxt84oEDBwgODmbnzp34+voyc2b1Hh1pqc4gXLElJvUCZn8H08zLm1mhoVU+lMRnJjBtz1ziMxP5uPssApx9gYIhW1nQLoQQJZNx8HfSt/5Y5ue17jua2v3//cDnmZqa8uyzz/Ltt9+yePFicnJy+N///seBAwdITU3lpZde4rfffqNevXokJiYSGBhIp06dcHZ2BgqCzZ49e/TnW7ZsGQcPHizSozI5OZlhw4YRGhpKx44d0el0pKamArB48WIcHR0BWLBgAXPmzCl22vJue/bs4a233tJ/XqtWLQC0Wi379u3DzMyMpKQk/fknT56MWq1mx44ddOzYEVNTU6BgBOupp54qcn4TExNatGjBgQMH6Nev3wO/j1XZA1d0//nnn8yePZvdu3fj5+dHaGiofoivutJlppN7+iBrnJ/n6NERuFtXj9G488mXmb5nDiZqE77ouwCP2m6VXZIQQlRLVh0HYN6ifZmf18jGvsTPHT9+PN26dWPhwoWEhobi7e2Nt7c3v//+O5cvX6Z///765yqKQmRkpD5gjRs3rkTXOHToED4+PnTs2BEouHnN3r6gxh9//JEffviB3NxcsrOz9ee+n9jY2CLPs7OzIz09nYULF/LBBx8wefJkUlNT+eOPP/QtnTZu3KgfwMjJyeHcuXPFjmABODs7ExsbW6LXV5UZDFj79u0jODiYPXv2EBAQwIYNGx6buwmzwnajUxRu+XSvNuHq8LXjBO//GDcbVz7sNgMHc7vKLkkIIaoto9oOJZrKK0/NmzencePGbNq0iZUrVzJ+/HigIEz5+vqyb98+g8c+6m4q+/fv5/PPP+fgwYPUqVOHX3/9lTlz5jzwOAsLC7Kzs7Gz++d3kL29PXv37sXExIQOHTpw69Ytvv76a8aNG4eZmRmKorB161YWLlwIwKlTp3BycqJu3brFXiMnJ6dc1o9VtCJ3Ee7Zs4fu3bvTrVs3MjMz2bRpE8ePH39swhVAwsHt7LJqzWj/kv+lUZl+vbCFd/bOJ8DZlyW950u4EkKIx8T48eP54IMPOHbsmH7KrGPHjly4cIFdu3bpn3fy5MlCi9jvZW1tza1bt4p8vWPHjoSHh3Pw4EEAdDodycnJpKSkYGNjg729PXl5eSxbtqxE9fr6+hIRUbgFkL29Pe+//z7Tpk1DpVKRm5vLqlWr9HcMHj16FG9vb30oNLT+6o7w8HD9urTqrFDA6tatGz179iQ3N5fff/+dI0eO8OSTT1ZWbeUiPy4Kk7gL7KzbkwENqvbolU7R8fVf3/PJ0f9jqGd/5naejrmxWWWXJYQQooz861//IjIykhEjRugDiJ2dHZs2bWLu3Ln4+fnh4+PD9OnT0el0Bs8zZcoUevTogb+/P/Hx8fqv29nZsWHDBt5++218fX1p1aoV+/fvp3///jRp0gQvLy/69u2Lv79/ieodMWIEf/zxR6GvOTg4ULt2bTp16gQUjMD169dPP8q1YcOGQkuLwsLC2Lt3Lx4eHvq3lStXAhAVFQVAixYtSlRPVaZS7jRPAn13dktLywfeJahSqYpNyxUpLCyMwMDAB6bhuyWEfk3cgR2EPvU9we2qbljJ0+bz4aEl7L66n1cDxjLSa3C1v3NTCCHK28P8XhAll56eTocOHThy5AiWlpYlOqZ58+bs3r37gXcoAkyfPp2mTZvqp0srQ1n9Gyq0BmvWrFmPXFhVpmg1pB3bTWjtbjzvXauyyzEoLTed9/Z+QETyRWZ3nkpXt46VXZIQQgiBtbU1ixcv5sqVKyUeZTp79myJz+/i4sLYsWMftrwqpUYFrJyIE9TKTiHarxcNbarOaFBoaChz5s0hMiKSRk0b4zjAldqBjnzWcy4t6kg3diGEEFVHr169yu3cEydOLLdzV7QatVXOzf3bOGvWkD5tmlR2KXqhoaEEBQURo4nDZVhj4pR49i3YyqCcbhKuhBBCiGqqxgQsbWYaqsgjbHbsxVD3qjN6NWfeHOxbOtF8ajtcBzSmxdT22Ld04svPvqjs0oQQQgjxkGpMwMo4vgtFUbAM6I6ZcdUJWJERkdi2rKNfwK5SqbBtWYeI8IgHHCmEEEKIqqrGBKwbB3aww6oNz/hVnR5SOkWHo1sdUk7Hc+dmTkVRSD2dgJe3TA8KIYQQ1dUDt8p5HORdv4xl/EX+av4vXnCoGqNXuZpc5h9ajE0/Z2KXHOPswiPYtqxD6ukEkk/Hszz0q8ouUQghhBAPqUaMYMUf2EGiUW38OrSp7FIASM1J482dMzl87TjLpvyXkJAQ3ExcuL7hEm4mLoRWg02nhRBCCGHYYz+CpWg1ZJ/YxSa7bkz0NK3scohNj2Pa7jlk5GWypNd8vB09oQEMHz68sksTQgghRBl57Eewss8dwzwnlYyWvahtWrnTg2cSInh161RUwP/1XVgQroQQQtRIHh4enDlzptDXunXrxubNm0t9Ln9/f7Kzs8uqtEKCg4Pvuw/i3VQqFRkZGQYfHz58OIcOHSqr0li3bp1+z8ODBw/SqFGjQo8rikLnzp25cuVKmV2zpB77gBX953bOmDViYNvGlVrHvuhDTN75Pm42rnzR9yNcrJ0rtR4hhBCPB41Gw8mTJzE3Ny+X88+ePbvEAet+jh49SmpqKh06dCiDqgps3LhRv89hWFgYrVq1KvS4SqVi8uTJzJ49u8yuWVKPdcDSZaVjdvEo+5y784Rz5Y1erYv4lZl/fsQTrm35pOdsateyqbRahBBCVA9btmwhICAAX19funbtyrlz5/SPqVQqPvnkE7p168Y777yjHzn6/fff8ff317+ZmZnx3Xffleh8H330Ee3ataNhw4Z8++23ALz88ssAdOzYUb+R9LPPPkvr1q3x9fVl4MCBhTaXvp9ly5YxevToUn0PnnnmGdasWaP//M033+TTTz8FID8/nwMHDtC9e3cAg3sHDho0iN9//5309PRSXftRPdZrsBKP7QOdDud23Splo2StTssXYSsJidzMKJ/hTPD/N2rVY51phRCiWkjKTiYpO6XMz+tgboeDuX2Jnz9ixAjMzMz0n1+8eBFAH2R2795Ny5Yt+fHHH3nqqacKTSnm5uayZ88eABYtWgTAgAEDGDBgAACrV6/mk08+Yfjw4SU6n5mZGUeOHCE8PJy2bdvy73//m6+++oply5Zx8OBBrKysAFi8eDGOjo4ALFiwgDlz5vD5558/8LXu2bOHt956q8TfGygITTNnziz0+ZAhQwDYvXs3HTt2xNS0YH11WFgYTz31VJFzmJiY0KJFCw4cOEC/fv1Kdf1H8VgHrLiDOzll1Yqn/Bwr/No5mlzmHviEg9eOMbnNywz17F/hNQghhCjerxe2sur0/8r8vGNa/ouxvqNK/Pz169cX2jS5W7duABw5cgR/f39atmwJwOjRo3nttdeIi4ujXr16AIwbN87geXfv3k1wcDD79u3D2tqaPXv2PPB8d0aXvL29MTY25saNG7i6uhY5948//sgPP/xAbm4u2dnZODuXbMlLbGxsoed27tyZ8PDwYp/7119/YWtrS1xcHJ6eBeuVdTodJ0+e1E8Dbty4UX/HfU5ODufOnSt2BAvA2dmZ2NjYEtVZVh7bgKVJvI7DzXNE+b+Ns0XFjl6l5KTy7p75XE69yvwu79LRtWq0hxBCCFFgcNO+dHJtW+bndTAvm2bWiqIUO/Ny99fujCjd68yZM4wdO5bffvsNFxeXEp/v7pE0IyMjNBpNkefv37+fzz//nIMHD1KnTh1+/fVX5syZU6LXZGFhQXZ2NnZ2Bd+jP//8877P37t3L76+vqjVBTM/p06dom7dutjY2KAoClu3bmXhwoX6x5ycnKhbt26x58rJySm3NWqGPLYB6+K+XShqcwKeKLvFdCURdSuGd/bMI0eTw5Le8/FyaFqh1xdCCPFgDub2pZrKq2gdOnRg/PjxhIeH4+3tzf/+9z9cXV0fOFp07do1hg4dyrfffkvz5s0f+XwA1tbW3Lp1CysrK1JSUrCxscHe3p68vDyWLVtW4tfk6+tLRESEPvQ9SEJCAvb2Bf+NNBoNc+bMITAwEChYMO/t7a0PmYbWX90RHh6On59fiWstC49lwFIUhfwTO9lv15EXG5VvYg0NDWXOvDlERkTi2qgBVn2d8O3uz6c951DPqvgkLYQQQtxPnTp1+OGHHxg9ejRarRZbW1vWrl37wOOWL19OQkICkydP1n9tzpw5DB48+KHOBzBlyhR69OiBubk527Zto0mTJnh5eeHq6krHjh3ZunVric4zYsQI/vjjD3r06FGi5/fs2ZNFixbRt29fPDw8UKlU+hC1YcMG/d2DULD+au/evXh4eOi/NnPmTMaNG0dUVBRAoanYiqBS7myCVw2FhYURGBhYJLmmXzzHrc/fJLTHB0wcbDjRPqrQ0FCCgoKwb+mEbcuCPQVTTiew5uc1jHqq5HPwQgghyoah3wui8qWnp9OhQweOHDmCpaXlI52refPm7N69Gycnpwc+d/r06TRt2pTx48eX6Nxl9W/osbylLXLPDq6ZONLnifIdDpwzbw72LZ1oPrUdrgMa02Jqe+xbOrHwo4Xlel0hhBCiurG2tmbx4sVl0vTz7NmzJQpXAC4uLowdO/aRr1laj13AUjR5WEfs42i9bnjZG5XrtSIiIrFtWUe/SFClUmHbsg4R4RHlel0hhBCiOurVq1eFT9VNnDhRv1C+Ij12Aeta2FEsNRk4tu9VrtfZE32AWs7mpJyO584sq6IopJ5OwMvbq1yvLYQQQoiq7bFb5B69fydp5o15srV7uZw/V5vHlydWsvHCH/QZ9yTr31/N2YVHsG1Zh9TTCSSfjmd56Fflcm0hhBBCVA+P1QiWNiONejHHuNK4BzblsLFzbNp1Xts6jd8v7WBym5dZ+973hISE4GbiwvUNl3AzcSE0NFTf+EwIIYQQNdNjNYJ1Zt9eais6fLp0L9PzKorC1iu7WXxsGQ7m9nzZdyFN7Qt27B4+fDjDhw8v0+sJIYQQonp7rAJW5vGdnLcNIKhp2XTSBUjKTuGTo19yIPYofRt2Z1KbCViYWJTZ+YUQQgjx+HlsAlZa3DXckiO4+MRU1GW0sfOuqD9ZfHwZapWa+V3e5YkG7crkvEIIIYR4vD02Aevkzp04qi3o0qPjI58rLuMmn59Ywf7YI3Rz68TkNi9ja2ZTBlUKIYQQoiZ4PAKWomB5dhd/OXditL3Zg59/293b3DTzasa7775Lrhf8cGYdNqbWzO48lW5uncqxcCGEEEI8jh6LgJV87Tpe2TeI6zupxMfcvc2Ny7DGRJ+O4+mRT9NiUltee+5VnmvxFBYmFbvzthBCiLIRHh5e2SWIaqqs/u08FgErI/IM103q0L1jyxIfc/c2NyqVivr9G3Hmo8Podmfw8mfPl2O1QgghyoujoyMWFhY8++yzlV2KqMYsLCxwdHR8pHNU+4BlamxE0xsnON2kP21NS7Y1Tr42n/DwcFyHNy20zY2drxOXN1wqz3KFEEKUIzc3N8LDw0lMTKzsUkQ15ujoiJub2yOdo9oHrB6dulNbm0HjLj0f+NyYtGtsu7KHPy7txPT2Njf1+zdCpVLJNjdCCPGYcHNze+RfjkI8qgrt5H7hwgU6duyIp6cnbdu25dy5c4/0PIAt+3bQ5MsdvDrmqSKP6RQdkUkX+eHMOl7e8hbPbnqVkMjNtKsfyMI5H5FyOoGzC48Q+/slzi48QvLpeGbOmPlQr+2nn356qOPKmtRRtWoAqeNeVaGOqlADSB33qgp1VIUaQOq4V1Woo9Q1KBWoe/fuyrfffqsoiqKsW7dOad++/SM978SJE0qjZ5orti3qKIAS0LO1suXSLuWrsFXKlJ0zlUHrRitdVg9W+v7vKeW9vR8ou6L2KzmaXP3xISEhin+Av2Jmbqb4B/groaGhD/3aBg0a9NDHliWpo2rVoChSx72qQh1VoQZFkTruVRXqqAo1KIrUca+qUEdpa6iwKcL4+HjCwsLYtm0bAEFBQbz++utERUXh4eFR6ucBZGdn4zqgMfX7N+L0R4c5f/kCHxxaTG1ja1zN6tHeOoAmdT3wsGiAscoIkuBs0hn98R4eHqz4ZkWhc4aFhT3U60tNTX3oY8uS1FG1apA6qmYdVaEGqaNq1lEVapA6qmYdqampZGVlYWFRwt1cyinoFXH8+HHF29u70NfatGmj7N2796GepyiKsnr1aqXL6sFKl9WDlUbPNFdQoYC8yZu8yZu8yZu8yVvZv504caLEuadCF7mr7tnCRlGUR3pe3759YWvBSJT5FHOYUjZ1CiGEEELcy8ur5DfCVVjAatCgAbGxsWg0GoyNjVEUhZiYmCJ3epT0eVBwG+Xo0aMr6iUIIYQQQpRIhd1F6OTkRKtWrVi9ejUAISEheHh4FFlXVdLnCSGEEEJUVSrF0PxbOYiMjGTMmDEkJSVhY2PDd999R/PmzQEYMGAAc+bMoXXr1vd9nhBCCCFEVVehAUsIIYQQoiao0EajQgghhBA1QbUNWKXp9l6eJk6ciIeHByqVijNnzjz4gHKQk5PD0KFD8fT0xN/fn379+hEVFVXhdfTp0wdfX1/8/f3p3LkzJ0+erPAa7jZ79uxK/e/i4eGBl5cX/v7++Pv78/PPP1dKHbm5ubz++us0bdqU5s2bV8omuKmpqfrvg7+/P56enhgbG5OcnFyhdWzdupXAwEBatWpFixYt+O677yr0+nds2bKF1q1b4+vrS/v27fn7778r5Lr3+3lVkT9T71dHRf5MNXStiv6Zer/XXFE/V0vyfa+In6n3q6Mif6ber44S/0z9//buPyqm/P8D+HP6ISlmlB81ldoNJaVJv2iLFM5qU7EiqyO/1u4R2V0cP5a2rJXD2rXsHw6t9eMoscnWxqaTIm2caGoV+VFJFBtFhWw/3t8/nOZrmpnMaLpT+3k9zplz9L533u/XvXd69fK+d+5V+oYOPYyyd3vvbufPn2eVlZXM0tKSXbt2TSMxvHz5kqWmprK2tjbGGGN79uxhU6ZM4TyOuro6yb+TkpKYk5MT5zG0u3r1Kvvwww/ZsGHDNHZcNPmZeNMXX3zBVqxYIfl8VFVVaTgixnbs2MH8/f05HbOtrY0ZGRmxwsJCxhhj5eXlTE9Pj9XX13MaR21tLTM2NmbXr19njDGWlZXFRo8ezcnYneUrLnNqZ3FwmVMVjcV1Tu1sm7nKq2/b71zl1M7i4DKndhaHsjm1VxZYjx49Ynw+nzU3NzPGXifOoUOHsvLyco3F1FP+mDLGWF5eHrO2ttZoDAcPHmTOzs4aGbupqYmNGzeOlZWVafS49ITPRGNjI+Pz+ayhoUGjcXRkZ2fHkpKSOB2zvcBqv2lxYWEhEwqF7NWrV295p3rl5eXJ3EzZ0NBQpRsYdlXHz6amcmpnvyNc/v68bSyucurb4uAir8qLQRM5VdMFlqIxVcmpvfIUYWVlJYRCIXR0Xt/Gi8fjYdiwYbh3756GI+sZdu/ejenTp2tk7Pnz58PCwgIbN27U2OmXyMhIhIaG4r333tPI+G+aN28eHBwcsGTJEtTU1HA+fmlpKYyNjbFlyxa4uLjAy8sLGRkZnMfxptzcXDx58gT+/v6cjsvj8XD8+HHMnDkTlpaW8PT0xKFDh9CnTx9O4xgxYgRqampw6dIlAEBSUhIaGxs1clq/HeXUzmkypwKaz6uUU/+fKjm1VxZYgPJ3e/9fs3XrVty+fRvfffedRsY/fPgwKisrsWXLFqxZs4bz8XNzc5GXl4dly5ZxPnZHFy5cQGFhIfLz82FsbIywsDDOY2hubkZZWRns7Oxw5coV/PzzzwgJCdFIYmp34MABzJ8/X/LHnCstLS2IiYnB77//joqKCmRkZCAsLIzz68D4fD4SExOxbt06ODs7IysrC3Z2dtDV1eU0jo4op8qn6ZwKaDavUk6VplJO7b6Jte7z6NEjNmDAADpF2MGOHTuYs7Oz1Dl7Terbty97/Pgxp2PGxMQwU1NTZmlpySwtLZm2tjYTCoXs9OnTnMbRUVVVFTM0NOR83JqaGqalpcVaWlokba6uriwzM5PzWBh7Pb3ev39/duPGDc7HlndqzsXFhZ07d47zWN7U1NTEBAIBu337NmdjyjtFqImc2tNPEXKdU5XZ5u7Oqx1j0FROfdu+4CqndoxDlZzaK2ew6G7vsn744QfEx8cjPT0dAoGA8/Hr6+tRVVUl+TkpKQnGxsYwMjLiNI5169ahqqoKd+/exd27d2Fubo60tDRMmzaN0zieP3+Op0+fSn6Oj4+Hk5MTpzEArx8n5evri7S0NABARUUFysvLYWNjw3ksAHDixAmMGTNGped5qUv7Y7hu3rwJALhz5w5KS0sxcuRIzmOprq6W/Pvbb7+Fj48Phg8fznkc7SinytJ0TgV6Rl6lnCpNpZza7eVfNykpKWHjxo1jI0aMYM7OzqyoqEgjcSxbtoyZmZkxbW1tNnToUI1cXF5ZWckAsPfff585OjoyR0dH5ubmxmkM9+7dY66ursze3p6NGTOG+fr6MrFYzGkM8mhqZrG0tJSJRCLm4ODA7O3tWUBAgMZmWEtLS9nEiROZvb09c3R0ZCdPntRIHIwx5unpyQ4cOKCx8ePi4iSfUQcHBxYfH6+ROBYvXsxsbGyYtbU1Cw0N5WyGpLN8xWVO7SwOLnOqorG4zqmK4uAyryq737s7pyqKg+uc2tn+UDan0p3cCSGEEELUrFeeIiSEEEII6cmowCKEEEIIUTMqsAghhBBC1IwKLEIIIYQQNaMCixBCCCFEzajAIoQQQghRMyqwCCGEEELUjAosQgghhBA1owKLEEJIp9ra2tDY2NgtfT9//rxb+iVE06jAIoQQohBjDD/++CN4PN47vT8zMxNLlixBSEgIYmJiZJbn5OTg+vXrXQ2TkB6HCixCCCEK7dy5E9OmTYOBgQEA4MyZMxAIBBCJRLC2tkZwcDCam5uRkpICLy8vmfdPmjQJsbGxOHbsGC5cuCCzfOrUqYiNjUVNTU23bwshXKICixBCiFy3bt1CZWUl7OzsJG1isRhhYWEoKChASUkJioqKkJycjIKCAohEIoV9HTx4EB999JHcZatWrUJUVJSaoydEs6jAIoQQIldUVBQ+//xzqTaxWAx7e3sAgK6uLkxNTdHc3AyxWKywwPrpp5/Q0NCA5cuXy11uZmaG2tpa3Lx5U63xE6JJVGAR0k2ioqLA4/EwYcIEucsMDQ05i+Xx48fg8Xg4ePCgxmJQF03Hrez47ce/48vW1paDKLuuqqoKZWVlGDVqlFT7mwXWxYsXUVxcDG9vbxQUFMDJyUmmn2PHjmH37t0oLi5GeHi4wvE+/fRTREdHq3cjCNEgHU0HQMh/XXZ2Ns6dOwcfHx9NhyJlyZIlCk/Z9GS9KW59fX2cO3dOpq03SEhIwJQpU6Ta6uvrUV5ejqVLl0JPTw+mpqY4ffo09PX1cf/+fYwePVqmn5CQEISEhLx1vIkTJ2L27Nl49uwZ+Hy+2raDEE2hAouQbmRgYAB7e3tER0f3uALL3Nwc5ubmmg5DZb0pbi0tLYwbN06tfTLG8O+//0JPT0+t/XZ05swZrFy5UqqtoKAAtra2uHbtmlT7+fPnYWtr26WYtLW1YWdnh7NnzyI4OPid+yGkp6BThIR0s8jISFy4cAFZWVmdrnfq1Ck4OTmhb9++MDExQXh4uNS9hxYsWAB7e3tkZWXByckJBgYGcHNzw9WrV2X62r9/P6ysrNCvXz/4+vrizp07Mut0PNXV3n9aWhocHBygr68PLy8vlJeXo7a2FnPmzMGAAQNgbW2NhIQEmf5yc3Ph4+MDAwMD8Pl8fPLJJ/jnn39Ujr+4uBh+fn4wNjZGv379YGNjg+3btyuMW537Ljc3FwEBARAKhTAwMIBIJMKRI0dktlVd3hb3m7GfPn0ajo6O0NPTQ3JysiTeqVOnYsCAAejfvz/c3d2Rnp4utT2dHRNFGGO4fPmyzIyUWCzG2LFjZdbv7PorVbi7uyMnJ6fL/RDSE1CBRUg38/Pzg6ura6ffkkpOTsbMmTMxcuRIJCUlYdOmTThy5AiCgoKk1nv48CEiIiKwZs0aJCQk4MWLF5gxYwaam5sl6/zxxx9YunQpJk2ahKSkJPj4+Ch1igYAqqursXbtWmzatAlHjx5FeXk55s2bh5CQENjb2yMxMRHOzs4IDQ1FRUWF5H25ubnw9vYGn89HQkIC9u3bh7y8PAQEBKgcf0BAAOrq6vDLL78gNTUVq1ev7vRmlOrcdxUVFfjggw8QGxuLlJQUfPzxx1i8eDEOHz6s1P6Tp6WlRerFGFMpbuD19VArV67EV199hT///BMikQg5OTnw9vbGq1evEBsbi8TERAQGBuLevXsAlD8m8lRXV6OhoQFmZmZS7YoKrIKCAiQnJ8PKykryys/PV3lfmZubo6SkROX3EdIjMUJIt/jmm2+YgYEBY4yxlJQUBoBlZWXJLGOMMScnJ+bm5ib1/ri4OAaAZWZmMsYYCwsLYzwejxUVFUnWSU9PZwBYdna2pM3d3Z15eXlJ9bV+/XoGgP36669y43uz/+vXr0va9uzZwwCwtWvXStrq6uqYtrY227Vrl6RtwoQJzMPDg7W1tUnaioqKGI/HY6mpqUrHX1NTwwCw5ORkufu0u/fdm9ra2lhzczNbunQpGz9+vMLxO4sTgMzryJEjSsfdHjsAdvnyZal1PTw8mJ2dHWtpaZE7vjLHRJHLly8zPp//1m1Ut0OHDjE7OzvOxyWkO9AMFiEc8Pf3x9ixY+V+S6qxsREFBQWYPXu2VHtwcDB0dHSQnZ0taRMKhVKnbdrvT3T//n0AQGtrK65evYoZM2ZI9TVr1iyl4hQKhVLfGhs5ciQAYPLkyZI2gUCAIUOGoLKyEgDw4sUL5OTkIDg4GK2trZKZGhsbG5iamiIvL0/p+I2NjWFpaYn169fj0KFDknZF1LnvAKCurg4RERGwtLSErq4udHV1sW/fPty6davTOBTR19dHXl6e1MvPz0+luAFg0KBBcHNzk/z84sULXLp0CWFhYdDW1pYZV5VjIk9jYyP69u37TtvcFQKBoNseyUMI16jAIoQjkZGRyMzMlPnj+fTpUzDGYGJiItWuo6MDY2Nj1NbWStoEAoHUOn369AEANDU1AQBqamrQ0tKCIUOGSK03dOhQpWJU1L+89vYx6+rq0Nraii+//FJSlLS/qqqqJIWYMvHzeDykpaVh1KhRCA8Ph4WFBZydneXeARxQ774DXl/vFB8fj9WrV+Ps2bPIy8vDokWLpNZRhZaWFlxcXKReRkZGKsUNQOZ41tXVoa2tDUKhUO64qhwTeRhjaG1tfYct7ho9PT3JKVRCejv6FiEhHAkMDIRIJEJ0dDQ8PT0l7QKBADweD48ePZJav6WlBU+ePIGRkZHSYwwePBg6OjoyFzJ37Fud2uPfsGGD3OuHBg0apFJ/NjY2OHHiBJqbm/HXX39hw4YNmD59Oh48eCBzcbs6911TUxNSU1Oxc+dOrFixQtLe1tamUvzKUDXujs8BFAgE0NLSQlVVVaf9v+sxMTAwwLNnz2Ta3/V5hIp0LKZevnwpeSQPIb0dzWARwqHIyEhkZGTg4sWLkjZDQ0OIRCIcP35cat3ExES0tLTIfb6bItra2hg7diySkpKk2n/77beuBd4JAwMDjB8/Hjdu3JCZrXFxcYGVldU79aurq4uJEydi3bp1qK+vl1tMqHPfvXr1Cq2trZKZLQBoaGiQfGNPnboad/s+P3z4sNyZpq4eE6FQiObmZpkvFzDG1Prq6MmTJzIX1hPSW9EMFiEcCgoKwpgxY5CRkSH1P/WoqCgEBQVh7ty5CAsLQ1lZGdavXw9fX194e3urNMbXX3+NwMBALFy4ECEhIbhy5Qri4uLUvCXSduzYAR8fH8yZMwchISEYOHAg7t+/j/T0dCxcuFDpbfj777+xatUqzJkzB9bW1nj27BliYmJgZWUFa2true9R177j8/lwdXXFtm3bJDOB27ZtA5/PV+rWBqrqatzbtm2Dj48PJk+ejGXLlmHgwIHIz8/HoEGDsGjRoi4dE3Nzc+jr66O6uhrDhw9X30a/xYMHDyTX/RHS29EMFiEc4vF4iIyMlGkPCAhAYmIiSkpKEBgYiOjoaISGhuLUqVMqjxEQEIC9e/ciIyMDQUFBSE9PR3x8vBqiV8zDwwMXL15EY2MjFi5cCD8/P2zevBn9+vVT6Q+0iYkJTExMEBMTg2nTpuGzzz6DhYUFzp49K/dibkC9+y4uLg7W1tYICwtDREQEZs2ahfnz56vcjzK6GrenpyeysrLA4/GwYMECzJw5E0lJSbC0tATQtWOipaUFd3f3t14MD7y+IalAIIBIJJK8Nm7cKFmmyixicXExxo8fr/T6hPRkPEZXFBJCCOlg+/btqKqqwq5duzpdb+vWraitrcX3338vd9nDhw+xe/dupca0trZGdna2wov3CelNaAaLEEKIjLlz58o8R1EesVgMZ2dnhcuUvcN7SUkJLC0tqbgi/xlUYBFCCJFhYWEBDw8PmduKdCQWi7F582bJ6UF/f3+pZU5OTkqNd/ToUURERHQpZkJ6EjpFSAghRK6HDx9i06ZN2L9/v9zl9fX1MDMzQ319vcwtHOrr6zF48GA0NDRIfTNTnqamJsyaNQspKSlqvxUEIZpCM1iEEELkMjExgUgkUjiLVVBQAFtbW7lFUUFBAUaNGvXW4goA9u7di+joaCquyH8KFViEEEIUCg8PR35+Pl6+fCmzTCwWSz1+qOMyZa6/KiwshFAoVHgdFyG9FRVYhBBCOrV8+XK5jwsSi8U4c+aM1C0aTp48KVmmTIE1ZMgQmWcyEvJfQNdgEUIIIYSoGc1gEUIIIYSoGRVYhBBCCCFqRgUWIYQQQoiaUYFFCCGEEKJmVGARQgghhKjZ/wFD0GasErIVWgAAAABJRU5ErkJggg==",
"text/html": [
"
"
],
"image/svg+xml": [
"\n",
"\n",
"\n"
]
},
"metadata": {},
"execution_count": 4
}
],
"cell_type": "code",
"source": [
"u = [states[i].points[end].u[1] for i = 1:length(P)]\n",
"θ = [states[i].points[end].theta[2] for i = 1:length(P)]\n",
"w = [states[i].points[end].u[3] for i = 1:length(P)]\n",
"\n",
"# set up the plot\n",
"plot(\n",
" xlim = (0, 16),\n",
" xticks = 0:1:16,\n",
" xlabel = \"Nondimensional Force \\$\\\\left(\\\\frac{PL^2}{EI}\\\\right)\\$\",\n",
" ylim = (0, 1.2),\n",
" yticks = 0.0:0.2:1.2,\n",
" ylabel = \"Nondimensional Tip Displacements\",\n",
" legend=:bottomright,\n",
" grid = false,\n",
" overwrite_figure=false\n",
" )\n",
"\n",
"plot!([0], [0], color=:black, label=\"Analytical\")\n",
"scatter!([0], [0], color=:black, label=\"GXBeam\")\n",
"plot!([0], [0], color=1, label=\"\\$ \\\\theta/(\\\\pi/2) \\$\")\n",
"plot!([0], [0], color=2, label=\"Vertical \\$\\\\left(w/L\\\\right)\\$\")\n",
"plot!([0], [0], color=3, label=\"Horizontal \\$\\\\left(-u/L\\\\right)\\$\")\n",
"\n",
"plot!(λ_a, θ_a*2/pi, color=1, label=\"\")\n",
"scatter!(λ, -4*atan.(θ/4)*2/pi, color=1, label=\"\")\n",
"\n",
"plot!(λ_a, η_a, color=2, label=\"\")\n",
"scatter!(λ, w/L, color=2, label=\"\")\n",
"\n",
"plot!(λ_a, -ξ_a, color=3, label=\"\")\n",
"scatter!(λ, -u/L, color=3, label=\"\")"
],
"metadata": {},
"execution_count": 4
},
{
"cell_type": "markdown",
"source": [
"---\n",
"\n",
"*This notebook was generated using [Literate.jl](https://github.com/fredrikekre/Literate.jl).*"
],
"metadata": {}
}
],
"nbformat_minor": 3,
"metadata": {
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.9.1"
},
"kernelspec": {
"name": "julia-1.9",
"display_name": "Julia 1.9.1",
"language": "julia"
}
},
"nbformat": 4
}