Lightcurve.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "try:\n",
    "    from PyFRS import cgs,fs,rs\n",
    "except ImportError:\n",
    "    import cgs,fs,rs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "Z={ 'z':         0.1,              # Redshift\n",
    "    'k':         0,                # Medium density profile index\n",
    "    'SSC':       'No',             # SSC\n",
    "    'XIC':       'No',             # EIC\n",
    "    'E0':        1.e53,            # Isotropic-equivalent energy in erg\n",
    "    'Gamma0':    20.,               # Initial bulk Lorentz factor\n",
    "    'theta_j':   5.,               # Jet opening angle\n",
    "    'theta_obs': 2.8,              # Viewing angle in deg\n",
    "    'n18':       1.,               # CNM density at 10^18cm\n",
    "    'p':         2.2,              # Electron energy distribution index\n",
    "    'epsilon_e': 0.1,              # Fraction of the shock energy into electrons\n",
    "    'epsilon_B': 0.01 }            # Fraction of the shock energy in magnetic field    \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "t=np.geomspace(1.e3,1.e8,500)\n",
    "nu=np.empty(t.shape)\n",
    "nu[:]=1.e15"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Forward-Shock"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "Fnu=fs.FS_flux(t,nu,**Z)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x125c612d0>]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh0AAAGYCAYAAADm9VARAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABBpElEQVR4nO3dd3xT5eIG8OckadKd0k3bMFpGWWW0UIoDHFxE5ToQlN0yRES9V+5VQa97FNSfXierQAtluHBcFFHcMrrYo0AZbWmhkyadSZuc3x9otYLYQpM34/l+PuePhiTn8YC8D2e8ryTLsgwiIiIiK1OIDkBERESugaWDiIiIbIKlg4iIiGyCpYOIiIhsgqWDiIiIbIKlg4iIiGyCpYOIiIhsgqWDiIiIbEIlOsCvLBYLiouL4ePjA0mSRMchIiKiVpBlGdXV1QgLC4NCcelzGXZTOoqLi6HT6UTHICIiostQWFiIiIiIS77HbkqHj48PgPOhfX192/W7DQYDdDqdVb6bfsPjbBs8zrbB42w7PNa2Ya3j/Ov3/jqOX4rdlI5fL6n4+vpa7Q+dNb+bfsPjbBs8zrbB42w7PNa2Ya3j3JpbI3gjKREREdkESwcRERHZhEuUDo1Gg6effhoajUZ0FKfG42wbPM62weNsOzzWtmEPx1mSZVkWtvffMRgM0Gq10Ov1vKZHRETkINoyfrvEmQ4iIiISj6WDiIiIbIKlg4iIiGyCpYOIiIhsgqWDiIiIbIKlg4iIiGyCpYOIiIhsgqWDiIiIbMJuFnwj11Td0Ii80hrkldagxNCAsmojKusaYbHIkHF+3jqNSokOnmr4e7nB30uD8A4e6BrghTA/d6iU7M1ERI6CpYNsSl/fiJ+OlWFbXgV2HC/HqYq6y/4uN6WEnqE+6BeuRd9wLQbqOiA61AcKxV+vdEhERLbH0kFWZ2wyY8vBEny2pxg/Hi2DyWxp8evBPhp0C/ZGRAcPBPlo4O+lgUohQZIACUB9oxmVtY04V2tCeY0RhefqkF9RB2OTBQeKDDhQZABQCADQerghvqs/EqICcF3PYHQJ9LL9fzAREV0U114hq9HXNWJtZj5St51CabWx+fXuwd64tkcQhkUFILZzB/h5qtv83RaLjNPn6nGgWI/9RXrsP63HroJzqDOZW7yvZ4gPRvUJwai+oejd0ReSxLMgRETtqS3jN0sHtbuGRjNW/HwSi78/jhpjEwAgxFeDcbE6jOkfhp6hPlbZb6PZggNFeuw4UYFteeXIOFGJJstvf7z7hvviiZt7IyEqwCr7JyJyRSwdJIQsy/hsbzEWbs7FGX0DACA61Af3XhuJW2PCoFbZ9qZPfV0jvsktwZaDZ/HdkTKYmixQSMCLd/TDhCGdbJqFiMhZsXSQzZ3VN+A/n+zH1sOlAIBwPw88Mqon/t4/zC5u7DxXa8Lzmw5h4+4iAMCjN/XEnOFRvNxCRHSF2jJ+80ZSumKb9hVjwcb9qG5oglqpwAPXd8O910bC3U0pOlqzDl5q/N/4/gjRumPx98fx8pdHUFljwuM397KLUkRE5ApYOuiymZoseOmLw0jdfgoA0D9Ci1fG9UePEOvcs3GlJEnCYzdFI8BLjRc+P4yUn0+iss6ERWNj4Mb5PoiIrI6lgy5LWbURs9dkY1dBFQDg/hFRmDeyh0NM1jXzmkh08FTj0Y/2YeOuIhjqG/H2xEF2dWaGiMgZ2f8IQXbneFkN7ly8DbsKquDjrsLyqXF49KZohygcvxobG4Glk2OhUSmw9XAppq7IhKGhUXQsIiKn5jijBNmFnPxKjF28HYWV9egc4IlP516Fkb1DRMe6LDf2DsGaGfHwcVch81QlJi3PwLlak+hYREROi6WDWm1bXjkmpWSgqq4R/XV++GjOMEQGeYuOdUWGdPXHe/cmIMBLjf1FekxYvhNlv5vIjIiI2g9LB7XKj0fLMD01Cw2NFozoGYQNs4Yi0FsjOla76B3miw33DkWwjwa5Z6tx97IdOPvLPCNERNR+WDroL/14tAwzV2fD2GTBjb2CsXRKLDzUznXTZfcQH7w/OwHhfh44UVaL8Ut3oLDy8hejIyKiC7F00CXtKjiH2WtyYGqy4G+9Q/DupFhoVM5VOH7VJdAL780eik7+niiorMPdS3fgZHmt6FhERE6DpYP+VF5pNaanZqG+0YxrewTh7YmDbD6Vua1FdPDE+7MTEBXkhWJ9A8Yv3YFjJdWiYxEROQXnHkHosp3R12Pqiszmm0YXT3L+wvGrUK07NtybgOhQH5RVG3HPsp04yuJBRHTFXGMUoTbR1zdi2spMFOsbEBnohVWJg+Glca155IJ8NFg/ayj6hPmiotaECSweRERXjKWDWjBbZPxjw24cLalBsI8GadOHwN9LLTqWEB281Fg7M75F8ThylsWDiOhysXRQC69sOYLvj5RBo1JgZeJg6Pw9RUcSys/zfPHoG36+eExczuJBRHS5WDqo2Wd7i7Hkh+MAgJfvikHfcK3gRPbBz1ON9BksHkREV4qlgwAAB4r0ePTDvQCA+4ZH4bYB4YIT2Zc/Fo8JLB5ERG3G0kHQ1zVi9pqc5tlGHxnVU3Qku+TnqcbaGUPRN9wXlb8Uj9yzBtGxiIgcBkuHi5NlGY98uBdFVecXcHvjnoFQKiTRseyW1tOtRfGYuDyDxYOIqJVYOlxc2vZT+OpQCdRKBd6ZOAhaDzfRkezer8WjX7i2uXgcPsPiQUT0V1g6XNj+03q89EUuAODxm6N542gbaD3dkD4jvrl4TErJ4MylRER/gaXDRRkaGjF33S6YzBaM6hOCacO6iI7kcH4tHs2XWlIycKKsRnQsIiK7xdLhop7+9CAKKusQ7ueBl8f2hyTxPo7LofV0w5rp8c1Tpk9cnoH8Ci4SR0R0MSwdLujzfWfw8e4iKCTgzQkDofXkfRxXooOXGukz49E92BtnDQ2YuDwDp8/ViY5FRGR3WDpcTKmhAU98sh8AMPe6bojt3EFwIucQ6K3B2lnxiAz0QlFVPSYuz8AZfb3oWEREdsWqpePHH3/EmDFjEBYWBkmS8Mknn1hzd/QXzj8euw9VdY3oF67FQzd0Fx3JqQT7uGPdrKHo5O+Jgso6TFyegVJDg+hYRER2w6qlo7a2Fv3798c777xjzd1QK6VnFOCHo+fXVXn97v5wU/JEV3sL1bpj3ax4hPt54GR5LSamZKC8xig6FhGRXbDqqDN69Gi88MILuOOOO6y5G2qF/IpavPT5YQDAYzdFo1uwj+BEziuigyfWzxqKUF935JXWYHJKBs7VmkTHIiISzu7+qWswGFpsRiP/lXilZFnG/I/2o77RjITIACTy8Vir6xTgifX3DkWQjwa5Z6sxZWUG9PWNomMREV0xo9F4wVjdWnZXOnQ6HbRabfOWnJwsOpLDW59ZiB0nKuDhpsSisTFQcJpzm+ga6IV1M+MR4KXGgSIDpq7MRHUDiwcRObbk5OQW47ROp2v1Z+2udBQWFkKv1zdvCxYsEB3JoZ3R1+OlL85fVvn3qJ7oFOApOJFr6R7ig7Wz4tHB0w17C6uQuCoLtcYm0bGIiC7bggULWozThYWFrf6s3ZUOX1/fFptGoxEdyWHJsownPj6AGmMTBnby42UVQaJDfbFmRjx83VXIyT+H6alZqDeZRcciIrosGo3mgrG6teyudFD7+WxvMb7NLYVaqcDLY2O4eqxAfcO1WD0jHt4aFTJOVuLeNdkwNrF4EJFrsWrpqKmpwZ49e7Bnzx4AwMmTJ7Fnzx4UFBRYc7cEoLzGiGc+OwgAePD6bugewqdVRBug80Pa9MHwVCvx07FyPLBuNxrNFtGxiIhsxqqlIzs7GwMHDsTAgQMBAPPmzcPAgQPx1FNPWXO3BOD5TYdwrq4R0aE+uG9ElOg49IvYzv5ImRoHtUqBrw+VYN77e2G2yKJjERHZhMqaXz5ixAjIMv9CtbWfj5Xj0z3FUEjAy3fFcBIwOzOsWyCWTB6E2Wty8L+9xfBwU2DhnXyqiIicH0cjJ9PQaMaTnx4AAExN6IKYCD+xgeiiro8OwRv3DIRCAt7PPo3nNh1iQScip8fS4WSW/nACJ8trEeSjwby/9RAdhy7h5n4d8cpd/QEAqdtP4ZUtRwQnIiKyLpYOJ3KqvBbvfJ8HAHjq1t7wdeeS9fZubGwEXri9LwDg3e+P4+1vjwlORERkPSwdTkKWZTz56QGYmiy4pnsgbo3pKDoStdLkoZ3xxM29AACvfnUUK34+KTgREZF1sHQ4iU37zuCnY+VQqxR47ra+kCTelOhIZl0biYdvPH857PlNh7Aug4+VE5HzYelwAoaGRjy/6RAA4P4RUega6CU4EV2Oh27ohtnDIwEAT3yyHx/vPi04ERFR+2LpcAKvfXUUpdVGdA30wn3DOSeHo5IkCfNvisbUhM6QZeDfH+zDlwfOiI5FRNRuWDoc3MFiPVbvOAUAeP62vnB3U4oNRFdEkiQ8M6YP7oqNgNki48H1u/HdkVLRsYiI2gVLhwOTZRnPfHYQFhm4NaYjru4eKDoStQOFQsKisTG4JaYjGs0y7luTgx3HK0THIiK6YiwdDuyzvcXIOnUOHm5KPP7L0w/kHJQKCf+9ewBu7BUMY5MFM9KykJN/TnQsIqIrwtLhoGqNTXjpi8MAgLnXRSHMz0NwImpvbkoF3p44CFd3C0SdyYzEVZk4UKQXHYuI6LKxdDiot7/LQ4nBiE7+nph5TaToOGQl7m5KLJsai8FdOqC6oQlTV2biWEm16FhERJeFpcMBnSyvxYqfzk8g9eStvXnzqJPzVKuwInEwYiK0qKw1YVJKBk6V14qORUTUZiwdDuj5TYdgMlswvEcQbuwVLDoO2YCvuxvSkoagZ4gPSquNmJSSgaKqetGxiIjahKXDwXybW4Jvc0uhUkh4akxvzjzqQjp4qZE+Mx6RgV4oqqrH5JQMlFUbRcciImo1lg4HYmwy47n/nZ95dPrVXREV5C04EdlakI8G6TPjEe7ngZPltZiyIgNVdSbRsYiIWoWlw4Gs/PkUTlXUIchHgwev7yY6DgkS5ueBtTPjEeSjQe7ZaiSuykKNsUl0LCKiv8TS4SDO6hvw1i/Lns+/KRo+XLbepXUJ9EL6jHj4ebphT2EVZqVlo6HRLDoWEdElsXQ4iJe35KLOZMagTn64Y2C46DhkB3qG+iAtaQi8NSrsOFGB+9fuQqPZIjoWEdGfYulwAHsLq7BxVxEA4OkxfaBQ8OZROq+/zg8rpsVBo1Lg29xSPPzeHpgtsuhYREQXxdJh52RZbl62/s6B4eiv8xMbiOxOfGQAlk6JhZtSwqZ9Z/D4xv2QZRYPIrI/LB127ov9Z5Gdfw7ubgo8clNP0XHITo3oGYw37hkIhQS8l12I5zcdZvEgIrvD0mHHGhrNSN58fn2V2ddGoaOW66vQn7u5X0csGhsDAFi57ST+u/WY4ERERC2xdNixVdtO4fS5eoT4ajB7ONdXob82Lk6HZ8b0BgC88c0xLP/xhOBERES/YemwU2XVRrzzXR4A4NFR0fBUqwQnIkeReFVXPDLq/KW4F784jHUZBYITERGdx9Jhp177+ghqjE2IidDyEVlqs/tHROG+4VEAgCc+2Y9P9xQJTkRExNJhlw6fMeC9rEIA51eR5SOy1FaSJOGxm3pi8tBOkGVg3vt78fWhEtGxiMjFsXTYGVmW8cLnh2CRgVv6dcTgLv6iI5GDkiQJz/29L+4cGA6zRcbcdbuwLa9cdCwicmEsHXZm6+FSbMurgFqpwPzR0aLjkINTKCS8fFcMRvUJganJglmrs5GTf050LCJyUSwddsTUZMFLX5x/RHb61V2h8/cUnIicgUqpwJsTBuKa7oGoM5mRtCoTh4oNomMRkQti6bAja3bm42R5LQK91Zh7XZToOORENCollk6JRVznDjA0NGHqygwcL6sRHYuIXAxLh504V2vCG1uPAgD+9beeXEWW2p2nWoWVSYPRN9wX5TUmTE7JwOlzdaJjEZELYemwE298cwyGhiZEh/pgfJxOdBxyUr7ubkhLGoJuwd44o2/ApJQMlBoaRMciIhfB0mEHTpbXIn1nPgDgP7f0hpKPyJIVBXhrkD4jHjp/D+RX1GHKikycqzWJjkVELoClww68/GUumiwyRvQMwtXdA0XHIRcQqnXH2hlDEeKrwZGSaiSuykR1Q6PoWETk5Fg6BMvJr8TmA2ehkIAFo3uJjkMupFOAJ9JnxMPfS429p/WYkZaNepNZdCwicmIsHQLJsowXPz//iOy4WB16hvoITkSupnuID1ZPHwIfjQqZJysxZ20OTE0W0bGIyEmxdAi0+cBZ7CqogoebEvP+1kN0HHJRfcO1WJU0GB5uSnx/pAz/fG83mswsHkTU/lg6BDE1WbDoy1wAwKxrIxHi6y44EbmyuC7+WDY1FmqlAl/sP4v5G/fDYpFFxyIiJ8PSIUj6znzkV9Qh0FuD2ddGio5DhGu6B+HNCQOhVEj4MOc0ntt0CLLM4kFE7YelQwB9fSPe/PYYAGDeyB7w0qgEJyI676a+oXh1XAwAIHX7Kbz29VHBiYjImbB0CPDu93moqmtEt2BvjI+LEB2HqIU7Bkbg+dv7AgDe+jYPS344LjgRETkLlg4bO32uDqu2nQIALBgdDZWSvwVkf6YM7dy8yvHCzbnNk9cREV0Jjng29uqWIzA1WZAQGYDro4NFxyH6U/cNj8ID13UDADz56QF8vPu04ERE5OhYOmxo/2k9PtlTDAB4/OZekCROd0727V9/64HEYV0gy8C/P9iHLQfPio5ERA7MJqXjnXfeQZcuXeDu7o74+HhkZmbaYrd2RZZlvPjFIQDA7QPC0C9CKzgR0V+TJAlP3dobd8VGwGyR8eC63fjpWJnoWETkoKxeOt577z3MmzcPTz/9NHbt2oX+/ftj1KhRKC0ttfau7cq3uaXYeaISapUC/x7VU3QcolZTKCQsvLMfRvcNhclswb2rc5B9qlJ0LCJyQFYvHa+99hpmzZqFpKQk9O7dG0uWLIGnpydWrlxp7V3bjSazBcmbz08ElnRVF0R08BSciKhtVEoF3rhnIIb3CEJ9oxlJqVk4UKQXHYuIHIxVS4fJZEJOTg5uvPHG33aoUODGG2/Ejh07LvoZg8HQYjMajdaMaBPvZ59GXmkNOni64f4R3UTHIbosapUCSybHYkgXf1Q3NGHaykwcL6sRHYuIbMxoNF4wVreWVUtHeXk5zGYzQkJCWrweEhKCs2cvfkOaTqeDVqtt3pKTk60Z0epqjE3NEyw9dEN3aD3cBCciunweaiVSEuPQN9wXFbUmTE7JwOlzdaJjEZENJScntxindTpdqz9rd0+vFBYWQq/XN28LFiwQHemKLPvxBMprjOgS4IlJ8Z1FxyG6Yr7ublg9PR7dgr1xRt+AySkZKK1uEB2LiGxkwYIFLcbpwsLCVn/WqqUjMDAQSqUSJSUlLV4vKSlBaGjoRT/j6+vbYtNoNNaMaFUlhgYs//EEAODRm6KhVtldxyO6LP5eaqTPiEdEBw+cqqjD1BWZqKoziY5FRDag0WguGKtby6qjoFqtRmxsLL755pvm1ywWC7755hskJCRYc9d24bWvjqK+0YxBnfwwuu/FSxaRowrVumPtzHgE+2iQe7YaiauyUGtsEh2LiOyY1f/pPW/ePCxfvhxpaWk4fPgw5syZg9raWiQlJVl710IdOVuND3LOn3J64hZOBEbOqXOAF9JnxsPP0w17Cqswa3U2GhrNomMRkZ2yeum4++678eqrr+Kpp57CgAEDsGfPHnz55ZcX3FzqbJI3H4ZFBkb3DUVsZ3/RcYispkeID9KShsBLrcT24xV4YN1uNJotomMRkR2SZFmWRYcAzj8qq9Vqodfr23R9yB5tyyvHpJQMqBQSvp43HF0DvURHIrK6HccrkLgqE8YmC24fEIbXxg+AQsEzfETOri3jN+9sbGcWi4zkzYcBAJOHdmbhIJeREBWAxZMHQaWQ8MmeYjz12QHYyb9piMhOsHS0s037z+BAkQHeGhUevJ4TgZFruT46BK/fPQCSBKTvLMDLW46IjkREdoSlox2Zmix49Ze/ZGdfG4kAb8d93Jfoco3pH4aX7ugHAFj8/XG8+32e4EREZC9YOtrR2ox8FFTWIchHgxnXdBUdh0iYCUM64fGbowEAL395BGt2nBIbiIjsAktHO6luaMRb357/F90/b+wOT7VKcCIise69Nqr5EuOTnx7Ex7tPC05ERKKxdLSTpT+cQGWtCZFBXrg7rvXz0BM5s3kjeyBxWBcAwL8/2IevDl58zSUicg0sHe2gxNCAlJ9/me58VDRUSh5WIgCQJAlP3dobYwdFwGyR8cC63diWVy46FhEJwtGxHfx36zE0NFoQ27kDRvVx7knPiNpKoZCwaGw/3NQnFCazBbNWZyMn/5zoWEQkAEvHFcorrcH72eenO18wOprTnRNdhEqpwBsTBuCa7oGoM5mRtCoTh88YRMciIhtj6bhCL3+ZC7NFxsjeIYjrwunOif6MRqXE0imxiOvcAYaGJkxZkYmT5bWiYxGRDbF0XIHsU5X46lAJFBLw6KieouMQ2T1PtQorEgejd0dflNcYMTklA0VV9aJjEZGNsHRcJlmWkbw5FwAwPk6H7iE+ghMROQathxtWzxiCyCAvFFXVY0pKBsqqjaJjEZENsHRcpq8OlSAn/xzc3RR4eGQP0XGIHEqgtwbpM+IR7ueBE+W1mLoyE/q6RtGxiMjKWDouQ5PZgpe/PH+WY8bVXRHi6y44EZHjCfPzQPrMeAR6a3D4jAFJqZmoMzWJjkVEVsTScRk+yDmN42W16ODphtnDo0THIXJYXQO9kD5zCLQebthVUIXZa3JgbDKLjkVEVsLS0UZ1pia8/vVRAMCD13eHr7ub4EREji061BerkgbDU63ET8fK8dD63WgyW0THIiIrYOloo5U/n0RptRE6fw9MGtpJdBwipzCoUwekTI2DWqXAloMlePTDfbBYZNGxiKidsXS0QUWNEUt+OD/d+b//1hMalVJwIiLnMaxbIN6ZOAhKhYSNu4vw7P8OQpZZPIicCUtHG7z1bR5qjE3oG+6LMTFhouMQOZ2RvUPw2vj+kCQgbUc+/u+ro6IjEVE7YulopYKKOqzNyAcAzL+pFxQKTndOZA23DQjH87f1BQC8/V0elv5wXHAiImovLB2t9OpXR9BolnFN90Bc3T1QdBwipzZ5aGc8dlM0ACB5c25z4Scix8bS0Qr7T+vx2d5iSBIwf3S06DhELmHOiCjcP+L8I+n/+eQAPt1TJDgREV0plo6/IMsyFn55GABw+4Bw9AnTCk5E5DoeGdUTU4Z2hiwD/3p/L745XCI6EhFdAZaOv/DjsXJsy6uAWqnAPE53TmRTkiTh2b/3wR0Dw9FkkTFn7S5sP14uOhYRXSaWjkuwWGQs/GVRtykJnaHz9xSciMj1KBQSXrkrBiN7h8DUZMGstGzsLjgnOhYRXQaWjkv4dG8RDp8xwMddhQeu6yY6DpHLUikVeGvCQFzVLQC1JjMSV2XhyNlq0bGIqI1YOv5EQ6MZr245P0fAnBFR6OClFpyIyLW5uymxbEocBnbyg76+EZNXZOBUea3oWETUBiwdfyJ9Zz6KquoR6uuO6Vd1FR2HiAB4aVRITRyC6FAflFUbMSklA2f09aJjEVErsXRchL6+EW9/lwcAmDeyB9zdON05kb3QerphzYx4dA30QlFVPSanZKCixig6FhG1AkvHRSz+/jiq6hrRI8QbY2MjRMchoj8I8tEgfWY8wrTuOF5Wi6krM6GvbxQdi4j+AkvHHxRX1WPVtpMAgEdHRUPJ6c6J7FK4nwfWzIxHgJcaB4sNmJGahXqTWXQsIroElo4/+O/WozA2WTCkiz9u6BUsOg4RXUJUkDdWzxgCH3cVsvPPYXZ6DoxNLB5E9oql43eOllTjw5zTAID5N0dDkniWg8je9QnTIjVpMDzclPjxaBn+uWEPmswW0bGI6CJYOn5n0eZcWGRgdN9QDOrUQXQcImql2M7+WDY1FmqlApsPnMX8jfthsciiYxHRH7B0/CLjRAW+yS2FUiHhkVE9Rcchoja6pnsQ3pwwEEqFhA9zTuO5TYcgyyweRPaEpQPnF3VL/mW68wlDdIgM8haciIgux019Q/Hy2BgAQOr2U3h96zHBiYjo91g6AGw+cBZ7CqvgqVbioRu6i45DRFdgbGwEnrutDwDgzW+OIeWnE4ITEdGvXL50NJoteGXLEQDAzGsiEezjLjgREV2pqQldmi+TvvD5YWzILBCciIgAlg5syCrEyfJaBHqrce+1kaLjEFE7uX9EFGb/8v/0go/34397iwUnIiKXLh21xia88cs134du6A5vjUpwIiJqL5IkYf7oaEyM7wRZBh5+bw++yy0VHYvIpbl06Vj+0wmU1xjRJcATE4Z0Eh2HiNqZJEl4/ra+uG1AGJosMu5Lz8HOExWiYxG5LJctHWXVRiz78fwNZo+Mioab0mUPBZFTUyokvDquP26IDoaxyYKZadnYd7pKdCwil+SyI+2b3xxDncmM/jo/3NwvVHQcIrIiN6UC70wahITIANQYmzB1ZSaOllSLjkXkclyydJwsr8X6X+5mXzCa050TuQJ3NyWWT4tDf50fquoaMTklAwUVdaJjEbkUlywdr245giaLjOt6BmFoZIDoOERkI94aFdKSBqNniA9Kq42YtGInzuobRMcichlWKx0vvvgihg0bBk9PT/j5+VlrN222p7AKn+8/A0kCHhsdLToOEdmYn6caa2YMQecATxRW1mPyigxU1ppExyJyCVYrHSaTCePGjcOcOXOstYs2k2UZyV8cBgCMHRSB6FBfwYmISIRgX3ekz4hHqK878kprMG1lJgwNjaJjETk9q5WOZ599Fg8//DD69etnrV202XdHSpFxshIalQLzRvYQHYeIBNL5eyJ9Zjz8vdTYX6THzNRs1JvMomMROTW7u6fDYDC02IxGY7t8r9kiY9Hm89OdJ17VBWF+Hu3yvUTkuLoFe2P19CHw0aiQeaoSc9bmwNRkER2LyK4ZjcYLxurWsrvSodPpoNVqm7fk5OR2+d6Nu07jSEk1tB5uuH94t3b5TiJyfH3DtViZNBjubgp8f6QMD7+3B2aLLDoWkd1KTk5uMU7rdLpWf7ZNpWP+/PmQJOmSW25ubpv/A36vsLAQer2+eVuwYMEVfR8ANDSa8drXRwEAD1zXDVpPtyv+TiJyHoO7+GPplDi4KSV8vv8MHt+4H7LM4kF0MQsWLGgxThcWFrb6s21abORf//oXEhMTL/meyMgrWzTN19cXvr7te4Nn6vZTOKNvQLifB6YkdG7X7yYi5zC8RxDevGcg5q7bhfeyC+HtrsJ/bunFeXyI/kCj0UCj0VzWZ9tUOoKCghAUFHRZOxKlqs6Ed7/LAwDMG9kD7m5KwYmIyF6N7tcRi8bG4JEP92HFzyfh6+6Gf9zYXXQsIqdhtWVVCwoKUFlZiYKCApjNZuzZswcA0K1bN3h7e1trtxc4WlIDpUJCdKgPbh8YbrP9EpFjGhenQ3VDE57bdAivbz0KH3cVpl/dVXQsIqcgyVa6cJmYmIi0tLQLXv/uu+8wYsSIC143GAzQarXQ6/XtfnnF0NCIUoMR3YJtV3aIyLG9sfUYXt96/l6wl8fGYPzg1t8sR+RK2jJ+W610tJU1SwcRUVvJsowXPz+MlJ9PQiEBb08chJv7dRQdi8jutGX8trtHZomI7IEkSXjill64Z7AOFhn4x4bd+P5IqehYRA6NpYOI6E9IkoQX7+iHW2I6otEs4770HGSdqhQdi8hhsXQQEV2CUiHh9fEDMKJnEBoaLZi1OhvHy2pExyJySCwdRER/Qa1SYPGkWAzQ+aGqrhFJq7JQUdM+SzQQuRKWDiKiVvBQK5EyLQ46fw8UVNZh1upsNDRygTiitmDpICJqpUBvDVYlDoGvuwq7Cqow7/09sHCdFqJWY+kgImqDbsHeWDb1/DotX+w/i0VfXtl6U0SuhKWDiKiNhkYG4OW7YgAAS388gfSd+YITETkGlg4iostwx8AIzBvZAwDw1KcH8F0u5/Ag+issHUREl+nB67vhrtgIWGTggXW7cLBYLzoSkV1j6SAiukySJOGlO/phWFQAak1mTE/Nwhl9vehYRHaLpYOI6AqoVQosnhyL7sHeKDEYkbQqC9UNjaJjEdkllg4ioiuk9XDDqqTBCPLRIPdsNeau241Gs0V0LCK7w9JBRNQOIjp4YsW0OHi4KfHj0TI89ekB2Mki3kR2g6WDiKidxET44c0JAyFJwPrMQiz54YToSER2haWDiKgdjewdgqdv7Q0AWPRlLjbtKxaciMh+sHQQEbWzxKu6IumqLgCAee/vRfapSrGBiOwESwcRkRX855beGNk7BKYmC2atzsap8lrRkYiEY+kgIrICpULCG/cMQP8ILc7VNSJxVSYqa02iYxEJxdJBRGQlnmoVUqYNRrifB05V1OHe1dloaDSLjkUkDEsHEZEVBflokJo0GD7uKmTnn8O/P9gLi4WP0pJrYukgIrKy7iE+WDo5FiqFhE37zuDVr46IjkQkBEsHEZENDOsWiIVjYwAA735/HOszCwQnIrI9lg4iIhu5KzYC/7ihOwDgP58cwI9HywQnIrItlg4iIhv6543dcefAcJgtMu5fuwuHzxhERyKyGZYOIiIbkiQJC8fGYGikP2qMTZiemoUSQ4PoWEQ2wdJBRGRjapUCSyfHISrIC2f0DUhalYUaY5PoWERWx9JBRCSA1tMNqUlDEOitxqEzBjy4bheazBbRsYisiqWDiEgQnb8nUqYNhrubAt8dKcMz/zsIWeYcHuS8WDqIiAQaoPPDf+8eCEkC0ncWYPlPJ0RHIrIalg4iIsFu6huKJ27uBQB46YtcfL7vjOBERNbB0kFEZAdmXN0VicO6AAAefn8PMk9Wig1EZAUsHUREdkCSJDx5a2+M6hMCU5MFs1ZnI6+0WnQsonbF0kFEZCeUCglv3DMQAzv5QV/fiGkrs1BazTk8yHmwdBAR2RF3NyVSpsahS4AniqrqMT01C7Wcw4OcBEsHEZGdCfDWIG36EAR4qXGgyIAHOIcHOQmWDiIiO9Q5wAsrEn+bw+PJTw9wDg9yeCwdRER2aoDOD29NGASFBKzPLMTb3+aJjkR0RVg6iIjs2MjeIXj2730AAP/39VF8lHNacCKiy8fSQURk56YkdMF9w6MAAI99tA8/HysXnIjo8rB0EBE5gEdH9cTf+4ehySLjvvQcHD5jEB2JqM1YOoiIHIBCIeGVcTEYGumPGmMTklZlobiqXnQsojZh6SAichAalRJLp8She7A3zhoakLQqC/r6RtGxiFqNpYOIyIFoPdyQOn0Ign00OFJSjfvW5MDUxDk8yDFYrXScOnUKM2bMQNeuXeHh4YGoqCg8/fTTMJlM1tolEZFLCPfzwKqkwfBSK7HjRAUe+2gf5/Agh2C10pGbmwuLxYKlS5fi4MGDeP3117FkyRI8/vjj1tolEZHL6BOmxeLJsVApJHy8uwivfnVEdCSivyTJNqzHr7zyChYvXowTJ05c8GsGgwFarRZ6vR6+vr62ikRE5NA+yC7EIx/uAwC8eEdfTIrvLDgRuZq2jN8qG2UCAOj1evj7+1/yPQZDy8fANBoNNBqNNWMRETmscXE6FFXV479bj+HJTw4g1NcdN/QKER2LnJjRaITRaGz++Y/j9qXY7EbSvLw8vPXWW5g9e/Yl36fT6aDVapu35ORkGyUkInJM/7ihO8bHRcAiAw+s2429hVWiI5ETS05ObjFO63S6Vn+2zZdX5s+fj0WLFl3yPYcPH0Z0dHTzz0VFRRg+fDhGjBiBlJSUi37m19MzhYWFLU7P8EwHEdFfazRbMCMtGz8eLUOgtxofzRmGzgFeomORE7rYmQ6dTteqyyttLh1lZWWoqKi45HsiIyOhVqsBAMXFxRgxYgSGDh2K1NRUKBQXP7nCezqIiK5MjbEJdy/dgYPFBnQJ8MSHc4Yh0Jv/aCPrasv4bdUbSYuKinDdddchNjYW6enpUCqVf/pelg4ioitXWt2AsYu3o7CyHjERWqyfNRReGpvevkcupi3jt9Xu6SgqKsKIESPQqVMnvPrqqygrK8PZs2dx9uxZa+2SiMjlBfu4Iy1pCPy91Nh3Wo/71+5Co5mTh5F9sFrp+Prrr5GXl4dvvvkGERER6NixY/NGRETWExnkjRXT4uDhpsQPR8s4eRjZDauVjsTERMiyfNGNiIisa2CnDnh30iAoFRI27irCy1s4eRiJx7VXiIic1HXRwVh4Zz8AwOLvjyN120nBicjVsXQQETmxcXE6PDKqJwDg2U2H8Pm+M4ITkStj6SAicnL3j4jC1ITOkGXg4ff2YMfxS097QGQtLB1ERE5OkiQ8PaYPRvcNhclswb2rs3H4TOunriZqLywdREQuQKmQ8PrdAzCkqz+qjU2YtjITp8/ViY5FLoalg4jIRbi7KbF8ahx6hvigtNqIqSszca7WJDoWuRCWDiIiF6L1cEPq9MHoqHXHibJazEjLQr3JLDoWuQiWDiIiF9NR64HV04dA6+GGXQVVeHD9LjRx1lKyAZYOIiIX1D3EByumxUGjUmDr4VL855MDnLyRrI6lg4jIRcV18cebEwZCIQEbsgrx+tZjoiORk2PpICJyYaP6hOL52/sCAN785hjWZuQLTkTOjKWDiMjFTYrvjIdu6A4AePKTA9hykKuBk3WwdBARER6+sTvuGayDRQYeWr8bWacqRUciJ8TSQUREkCQJL9zeFzf2CoaxyYLpqVmctZTaHUsHEREBAFRKBd6aMAiDu3RAdUMTpq7MREEFZy2l9sPSQUREzTzUSqRMHYzoUB+UVRsxZWUGyqqNomORk2DpICKiFrSebkibPgQRHTyQX1GHaSszYWhoFB2LnABLBxERXSDE1x3pM+IR6K3GoTMGzErLRkMjp0unK8PSQUREF9Ul0AupSUPgo1Eh42QlHly/m9Ol0xVh6SAioj/VN1yL5dPioFYp8PWhEjz+8X5Ol06XjaWDiIguaWhkAN76Zbr097NPY9GXR0RHIgfF0kFERH9pVJ9QLLwzBgCw5IfjWP7jCcGJyBGxdBARUauMH6zD/NHRAIAXvziMD3NOC05Ejoalg4iIWm32tZGYdU1XAMBjH+3D1kMlghORI2HpICKiVpMkCY/f3AtjB0XAbJExd90uZJ7kOi3UOiwdRETUJpIkYdHYfs3rtMxIy8LBYr3oWOQAWDqIiKjNVEoF3p44CEO6+J9fp2VFJk6U1YiORXaOpYOIiC6Lu5sSKYlx6BPmi4paEyanZKCoql50LLJjLB1ERHTZfN3dsHr6EEQFeaFY34ApKVwgjv4cSwcREV2RAG8N0mfGI9zPAyfKazF1ZSb09Vwgji7E0kFERFeso9YD6TPjEeitweEzBkxPzUKdqUl0LLIzLB1ERNQuugZ6Yc2MIfB1VyEn/xxmr8mBsYkr09JvWDqIiKjd9Oroi9TpQ+CpVuKnY+X454Y9XJmWmrF0EBFRuxrUqQOWTYmDWqnA5gNnsWDjflgsXJmWWDqIiMgKru4eiDcnDIRSIeGDnNN44fPDkGUWD1fH0kFERFZxU99QvDz2/Mq0K7edxJvf5AlORKKxdBARkdWMjY3AM2N6AwBe33oUK38+KTgRicTSQUREVpV4VVfMG9kDAPDcpkP4ILtQcCIShaWDiIis7sHru2Hm1V0BAI99tA+b958RnIhEYOkgIiKrkyQJT9zSC3fH6WCRgYc27MZ3uaWiY5GNsXQQEZFNSJKEl+7sh1tjOqLRLGN2eg625ZWLjkU2xNJBREQ2o1RIeP3uARjZOwSmJgtmpmUj82Sl6FhkIywdRERkU25KBd6eOBDDewShvtGM6alZ2FNYJToW2QBLBxER2ZxGpcTSKbFIiAxAjbEJU1dk4GCxXnQssjKWDiIiEsLdTYmUaXGI69wBhoYmTFmRiaMl1aJjkRVZtXT8/e9/R6dOneDu7o6OHTtiypQpKC4utuYuiYjIgXhpVFiZNBgxEVpU1powKSUDJ8trRcciK7Fq6bjuuuvw/vvv48iRI/joo49w/Phx3HXXXdbcJRERORhfdzesnj4E0aE+KKs2YuLynSisrBMdi6xAkm24As9nn32G22+/HUajEW5ubi1+zWAwQKvVQq/Xw9fX11aRiIjITpTXGHHPsp3IK62Bzt8D789OQEeth+hY9BfaMn7b7J6OyspKrF27FsOGDbugcPyewWBosRmNRltFJCIigQK9NVg7Mx6dAzxRWFmPScszUFrdIDoW/YHRaLxgrG4tq5eOxx57DF5eXggICEBBQQE+/fTTS75fp9NBq9U2b8nJydaOSEREdiLE1x3rZg1FuJ8HTpTXYnJKBiprTaJj0e8kJye3GKd1Ol2rP9vmyyvz58/HokWLLvmew4cPIzo6GgBQXl6OyspK5Ofn49lnn4VWq8WmTZsgSVKLz/x6eqawsLDF6RmNRgONRtOWiERE5ODyK2oxfukOlBiM6BPmi3Uzh0Lr+ednycl2jEZji6sQBoMBOp2uVZdX2lw6ysrKUFFRccn3REZGQq1WX/D66dOnodPpsH37diQkJLT4Nd7TQUREv5dXWoN7lu1AeY0JA3R+SJ8ZD2+NSnQs+oO2jN9t/t0LCgpCUFDQZQWzWCwAwPs0iIjoL3UL9kb6zHjcs2wn9hRWYfqqLKROHwxPNYuHo7LaPR0ZGRl4++23sWfPHuTn5+Pbb7/FhAkTEBUVdcFZDiIioouJDvXFmunx8NGokHmqEjNSs1FvMouORZfJaqXD09MTGzduxA033ICePXtixowZiImJwQ8//MB7NIiIqNX6RWiROn0IvNRK7DhRgRlpWSweDsqm83RcCu/pICKiS8nJr8TUFZmoNZlxdbdApEyLg7ubUnQsl2eX83QQERFdidjO/kibPgSeaiV+zivHrNXZaGjkGQ9HwtJBREQOI66LP1KTzhePn46xeDgalg4iInIoQ7r6Y1Xi4ObiMXtNDouHg2DpICIihxMfGYCViYPh4abED0fLcF86i4cjYOkgIiKHNPSX4uHupsD3R8owJz0HxiYWD3vG0kFERA4rIeq34vHdkTLMSd/F4mHHWDqIiMihDYsKxMppg6FRKfBtbinmrmXxsFcsHURE5PCGdQvEysTzxWPr4VLMXbsbpiaL6Fj0BywdRETkFK7qFogV034tHiWYu24Xi4edYekgIiKncXX3QCyfGge1SoGvD5XgARYPu8LSQURETuXaHkHNxeOrQyW4fy2farEXLB1EROR0hvcIQsrUuOZ7PO5dzXk87AFLBxEROaVrewRh1e8mEJuZls3VaQVj6SAiIqc1rFsgUpMGNy8Sl5SaiVpjk+hYLoulg4iInFp8ZADWzBgCH40KO09UYtrKTFQ3NIqO5ZJYOoiIyOnFdvbHmpnx8HVXITv/HKasyIS+nsXD1lg6iIjIJQzQ+WHdrKHw83TDnsIqTE7JQFWdSXQsl8LSQURELqNvuBbrZw2Fv5ca+4v0mLA8A5W1LB62wtJBREQupVdHX2y4dygCvTU4fMaACct2oqzaKDqWS2DpICIil9MjxAfvzR6KEF8NjpRU455lO1BqaBAdy+mxdBARkUuKCvLGe/cmIEzrjuNltbh72U6c0deLjuXUWDqIiMhldQn0wnuzExDu54GT5bUYv3QHCivrRMdyWiwdRETk0nT+nnj/vgR0DvBEYWU97lqyHXml1aJjOSWWDiIicnnhfh74YHYCeoR4o8RgxPilO3GwWC86ltNh6SAiIgIQ7OuODfcmoF+4FpW1JkxYthM5+edEx3IqLB1ERES/8PdSY+2seAzu0gGGhiZMWZGB7XnlomM5DZYOIiKi3/F1d0Pa9CG4pnsg6kxmJKZm4dvcEtGxnAJLBxER0R94qlVImRaHkb1DYGqy4N7VOdi0r1h0LIfH0kFERHQRGpUS704ahNsGhKHJIuOh9bvxfnah6FgOjaWDiIjoT7gpFXht/ABMGKKDRQYe/XAfUredFB3LYbF0EBERXYJSIeGlO/ph5tVdAQDP/O8Q3vkuT3Aqx8TSQURE9BckScITt/TCP27oDgB4ZcsRLNycC1mWBSdzLCwdRERErSBJEh4e2QOP3xwNAFjyw3Es2LgfZguLR2uxdBAREbXBvddGYeGd/aCQgA1ZhZi7dheMTWbRsRwCSwcREVEb3TOkE96ZOAhqpQJfHjyL6alZqDE2iY5l91g6iIiILsPofh2RmjQYXmoltuVVYOLynaisNYmOZddYOoiIiC7TsG6BWH/vUPh7qbHvtB53LdmOoqp60bHsFksHERHRFYiJ8MP7sxMQpnXHibJa3LV4O/JKa0THskssHURERFeoW7A3PpwzDFFBXjijb8C4Jduxt7BKdCy7w9JBRETUDsL8PPDBfcPQP0KLc3WNmLh8J7ZxhdoWWDqIiIjaib+XGmtnDcVV3QJQazIjaVUWvjxwRnQsu8HSQURE1I68NSqsTByM0X1DYTJbcP/aXdiQWSA6ll1g6SAiImpnGpUSb08c1LxQ3PyN+7H4++MuP206SwcREZEV/LpQ3P0jogAAi77MxUtfHIbFhadNZ+kgIiKyEkmS8OhN0Xji5l4AgOU/ncS/PtgLU5NFcDIxbFI6jEYjBgwYAEmSsGfPHlvskoiIyG7MujYSr47rD6VCwse7izAjzTWnTbdJ6Xj00UcRFhZmi10RERHZpbtiI7BiWhw81Ur8dKwc9yzbgdLqBtGxbMrqpWPz5s346quv8Oqrr1p7V0RERHZtRM9grJ81FAFeahwoMmDs4u04WV4rOpbNWLV0lJSUYNasWVizZg08PT1b9RmDwdBiMxqN1oxIRERkU/11fvhozjB08vdEYWU9xi7ejj0ONHup0Wi8YKxuLauVDlmWkZiYiPvuuw9xcXGt/pxOp4NWq23ekpOTrRWRiIhIiC6BXvhozjD0C9eistaECct24rvcUtGxWiU5ObnFOK3T6Vr92TaXjvnz50OSpEtuubm5eOutt1BdXY0FCxa06fsLCwuh1+ubt7Z+noiIyBEE+Wiw4d6huLZHEOobzZi5OhvvZxeKjvWXFixY0GKcLixsfWZJbuNMJWVlZaioqLjkeyIjIzF+/Hj873//gyRJza+bzWYolUpMmjQJaWlpLT5jMBig1Wqh1+vh6+vblkhEREQOq9FswWMf7cPGXUUAgH+N7IEHru/WYvy0Z20Zv9tcOlqroKCgxXWe4uJijBo1Ch9++CHi4+MRERFx2aGJiIiciSzLeGXLEbz7/XEAwOShnfDs3/tCqbD/4tGW8VtlrRCdOnVq8bO3tzcAICoq6oLCQURE5Mp+nUQsxNcdz/zvINJ3FqDUYMSbEwbC3U0pOl674YykREREdmLasC54Z+IgqFUKfHWoBJNTMlBVZxIdq93YrHR06dIFsixjwIABttolERGRw7m5X0esmT4EPu4qZOefw9jF21FYWSc6VrvgmQ4iIiI7Ex8ZgA/vG4aOWnccL6vFHe9uw14Hmsvjz7B0EBER2aGeoT74+P6r0KujL8prTLhn2U58fahEdKwrwtJBRERkp0K17vjgvgQM/2Uuj9lrspG2/ZToWJeNpYOIiMiOeWtUSJkWhwlDdLDIwNOfHcQLmw7BYrHKjBdWxdJBRERk59yUCrx0Rz88MqonACDl55OYu24XGhrNgpO1DUsHERGRA5AkCXOv64Y37hkAtVKBzQfOYuLynaiocZyFUVk6iIiIHMhtA8KxZsYQaD3csKugCncu3o6T5bWiY7UKSwcREZGDiY8MwEdzhkHn74H8ijrc+e42ZJ+qFB3rL7F0EBEROaBuwd7YOOcq9Nf54VxdIyamZGDTvmLRsS6JpYOIiMhBBflosGHWUIzsHQJTkwUPrNuNJT8ch5XWcr1iLB1EREQOzEOtxJLJsUi6qgsAYOHmXPznkwNoMlvEBrsIlg4iIiIHp1RIeHpMHzx1a29IErA2owAzV2ejuqFRdLQWWDqIiIicxPSru2LJ5Fi4uynw/ZEyjFuyA8VV9aJjNWPpICIiciKj+oTi/dkJCPLRIPdsNW57Zxv2na4SHQsASwcREZHTiYnwwydzr0J0qA/Kqo0Yv3QHthw8KzoWSwcREZEzCvfzaF4srqHRgvvSc7D8xxNCn2xh6SAiInJSPu5uWDEtDpOHdoIsA0t+OI7KWpOwPCpheyYiIiKrUykVeP62vuga6I0BOi0CvDXisgjbMxEREdmEJEmYcXVX0TF4eYWIiIhsg6WDiIiIbIKlg4iIiGyCpYOIiIhsgqWDiIiIbIKlg4iIiGyCpYOIiIhsgqWDiIiIbIKlg4iIiGzCJUqH0WjEM888A6PRKDqKU+Nxtg0eZ9vgcbYdHmvbsIfjLMkil5v7HYPBAK1WC71eD19fX4f5bvoNj7Nt8DjbBo+z7fBY24a1jnNbvtclznQQERGReCwdREREZBN2s8rsr1d5DAZDu3/3r99pje+m3/A42waPs23wONsOj7VtWOs4//p9rblbw27u6Th9+jR0Op3oGERERHQZCgsLERERccn32E3psFgsKC4uho+PDyRJEh2HiIiIWkGWZVRXVyMsLAwKxaXv2rCb0kFERETOjTeSEhERkU2wdBAREZFNsHQQERGRTTh16Vi8eDFiYmLg6+sLX19fJCQkYPPmzaJjObWFCxdCkiT885//FB3F6TzzzDOQJKnFFh0dLTqWUyoqKsLkyZMREBAADw8P9OvXD9nZ2aJjOZUuXbpc8OdZkiTMnTtXdDSnYjab8eSTT6Jr167w8PBAVFQUnn/++VY93moNdjNPhzVERERg4cKF6N69O2RZRlpaGm677Tbs3r0bffr0ER3P6WRlZWHp0qWIiYkRHcVp9enTB1u3bm3+WaVy6v+FhTh37hyuuuoqXHfdddi8eTOCgoJw7NgxdOjQQXQ0p5KVlQWz2dz884EDBzBy5EiMGzdOYCrns2jRIixevBhpaWno06cPsrOzkZSUBK1Wi4ceesjmeZz6b6wxY8a0+PnFF1/E4sWLsXPnTpaOdlZTU4NJkyZh+fLleOGFF0THcVoqlQqhoaGiYzi1RYsWQafTYdWqVc2vde3aVWAi5xQUFNTi54ULFyIqKgrDhw8XlMg5bd++HbfddhtuueUWAOfPMK1fvx6ZmZlC8jj15ZXfM5vN2LBhA2pra5GQkCA6jtOZO3cubrnlFtx4442iozi1Y8eOISwsDJGRkZg0aRIKCgpER3I6n332GeLi4jBu3DgEBwdj4MCBWL58uehYTs1kMiE9PR3Tp0/nPE3tbNiwYfjmm29w9OhRAMDevXvx888/Y/To0ULyOPWZDgDYv38/EhIS0NDQAG9vb3z88cfo3bu36FhOZcOGDdi1axeysrJER3Fq8fHxSE1NRc+ePXHmzBk8++yzuOaaa3DgwAH4+PiIjuc0Tpw4gcWLF2PevHl4/PHHkZWVhYceeghqtRrTpk0THc8pffLJJ6iqqkJiYqLoKE5n/vz5MBgMiI6OhlKphNlsxosvvohJkyaJCSQ7OaPRKB87dkzOzs6W58+fLwcGBsoHDx4UHctpFBQUyMHBwfLevXubXxs+fLj8j3/8Q1woF3Hu3DnZ19dXTklJER3Fqbi5uckJCQktXnvwwQfloUOHCkrk/P72t7/Jt956q+gYTmn9+vVyRESEvH79ennfvn3y6tWrZX9/fzk1NVVIHqc/06FWq9GtWzcAQGxsLLKysvDGG29g6dKlgpM5h5ycHJSWlmLQoEHNr5nNZvz44494++23YTQaoVQqBSZ0Xn5+fujRowfy8vJER3EqHTt2vOBsaK9evfDRRx8JSuTc8vPzsXXrVmzcuFF0FKf0yCOPYP78+bjnnnsAAP369UN+fj6Sk5OFnLlz+tLxRxaLBUajUXQMp3HDDTdg//79LV5LSkpCdHQ0HnvsMRYOK6qpqcHx48cxZcoU0VGcylVXXYUjR460eO3o0aPo3LmzoETObdWqVQgODm6+0ZHaV11d3QXroSiVSlgsFiF5nLp0LFiwAKNHj0anTp1QXV2NdevW4fvvv8eWLVtER3MaPj4+6Nu3b4vXvLy8EBAQcMHrdGX+/e9/Y8yYMejcuTOKi4vx9NNPQ6lUYsKECaKjOZWHH34Yw4YNw0svvYTx48cjMzMTy5Ytw7Jly0RHczoWiwWrVq3CtGnT+Pi3lYwZMwYvvvgiOnXqhD59+mD37t147bXXMH36dCF5nPp3ubS0FFOnTsWZM2eg1WoRExODLVu2YOTIkaKjEbXZ6dOnMWHCBFRUVCAoKAhXX301du7cecGjh3RlBg8ejI8//hgLFizAc889h65du+K///2vuBvvnNjWrVtRUFAgbAB0BW+99RaefPJJ3H///SgtLUVYWBhmz56Np556SkgerjJLRERENuEy83QQERGRWCwdREREZBMsHURERGQTLB1ERERkEywdREREZBMsHURERGQTLB1ERERkEywdREREZBMsHURERGQTLB1ERERkEywdREREZBMsHURERGQT/w8fzgMKfuzueQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(np.log10(t),np.log10(Fnu))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saving figure lc.png\n"
     ]
    }
   ],
   "source": [
    "fig,ax=plt.subplots(1,1)\n",
    "ax.plot(t,Fnu)\n",
    "ax.set(xscale='log', xlabel=r'$t$ (s)',\n",
    "       yscale='log', ylabel=r'$F_\\nu$[$10^{9}$ Hz] (mJy)')\n",
    "\n",
    "fig.tight_layout()\n",
    "print(\"Saving figure lc.png\")\n",
    "fig.savefig(\"lc.png\")\n",
    "plt.close(fig)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Reverse-Shock"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "Fnu=rs.RS_flux(t,nu,**Z)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(-10.0, 10.0)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAGdCAYAAADnrPLBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAvgklEQVR4nO3deXwU9eH/8ffm2ghkw5UDZCFcEm6QIwStQEHBUistP6vIfXkU0BK1Er8ICNJAq996YPEAExUQtVU8qiggwtdyJJwCAhIEEo4EQchCkE3Izu+PPly7EsKV3dnZvJ6PxzwezGdnZt870s6bmdlZm2EYhgAAACwqzOwAAAAAV4MyAwAALI0yAwAALI0yAwAALI0yAwAALI0yAwAALI0yAwAALI0yAwAALI0yAwAALI0yAwAALM2vZWb16tW67bbbVL9+fdlsNi1ZssTndcMwNGXKFNWrV0/XXHON+vTpoz179lx0uy+88IKSkpIUHR2tlJQUZWdn++kTAACAYOfXMlNcXKz27dvrhRdeKPf1v/zlL3ruuef04osvav369apevbr69u2rs2fPXnCbb731ltLS0jR16lRt2rRJ7du3V9++fXX06FF/fQwAABDEbIH6oUmbzab33ntPAwYMkPSfszL169fXQw89pIcffliSVFRUpISEBGVlZemuu+4qdzspKSnq0qWL5syZI0nyeDxyOp2aMGGCJk2aFIiPAgAAgkiEWW+8b98+FRQUqE+fPt6x2NhYpaSkaO3ateWWmZKSEm3cuFHp6enesbCwMPXp00dr164t9308Ho/279+vyMhI2Ww277jdbpfdbq/ETwQAAK6U2+2W2+32zhuGodLSUiUlJSksrOILSaaVmYKCAklSQkKCz3hCQoL3tZ87duyYysrKyl1n165d5a5z+PBhNW3atBISAwCAQMvPz1eDBg0qXMa0MhMoMTExkqSvv/7a+2epcs7MuFwuOZ1O5efny+FwXNW2cGHs58BgPwcG+zlw2NeBUVn7+ednZk6dOqVWrVr5HLsvxLQyk5iYKEkqLCxUvXr1vOOFhYXq0KFDuevUrVtX4eHhKiws9BkvLCz0bu/nfry0dO211/rtL7PD4eB/KAHAfg4M9nNgsJ8Dh30dGJW9n10ulyT53CJyIaY9Z6Zx48ZKTEzUihUrvGMul0vr169XampquetERUWpU6dOPut4PB6tWLHigusAAIDQ5tczM6dPn1Zubq53ft++fdqyZYtq166thg0b6o9//KOefPJJNW/eXI0bN9bjjz+u+vXre7/xJEm9e/fWb3/7W40fP16SlJaWpuHDh6tz587q2rWrnnnmGRUXF2vkyJH+/CgAACBI+bXMbNiwQb169fLOp6WlSZKGDx+urKws/elPf1JxcbHuuecenTx5UjfeeKOWLl2q6Oho7zp79+7VsWPHvPN33nmnvvvuO02ZMkUFBQXq0KGDli5det5NwYFgt9s1depUvhXlZ+znwGA/Bwb7OXDY14ERDPs5YM+ZMYvL5VJsbKyKioq4ZgoAgEVczvGb32YCAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWZnqZSUpKks1mO28aN25cuctnZWWdt2x0dHSAUwMAgGARYXaAnJwclZWVeee3b9+um2++WXfccccF13E4HNq9e7d33maz+TUjAAAIXqaXmbi4OJ/5WbNmqWnTpurRo8cF17HZbEpMTPR3NAAAYAGmX2b6byUlJVqwYIFGjRpV4dmW06dPq1GjRnI6nbr99tu1Y8eOi27b5XL5TG63uzKjAwCAq+B2u887Vl+qoCozS5Ys0cmTJzVixIgLLtOiRQu9+uqrev/997VgwQJ5PB51795dBw8erHDbTqdTsbGx3ikjI6OS0wMAgCuVkZHhc5x2Op2XvK7NMAzDj9kuS9++fRUVFaUPP/zwktcpLS1Vy5YtNWjQIM2YMeO8110ul2JjY5Wfny+Hw+Edt9vtstvtlZIbAABcHbfb7XPVxOVyyel0qqioyOf4XR7T75n50YEDB7R8+XK9++67l7VeZGSkOnbsqNzc3AqXczgcF90ZAADAHFdzkiFoLjNlZmYqPj5e/fv3v6z1ysrKtG3bNtWrV89PyQAAQDALijLj8XiUmZmp4cOHKyLC92TRsGHDlJ6e7p2fPn26PvvsM3377bfatGmThgwZogMHDmjMmDGBjg0AAIJAUFxmWr58ufLy8jRq1KjzXsvLy1NY2E+d68SJExo7dqwKCgpUq1YtderUSWvWrFGrVq0CGRkAAASJoLoB2B9+vAH4Um4gAgAAweFyjt9BcZkJAADgSlFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApVFmAACApZleZqZNmyabzeYzJScnV7jOO++8o+TkZEVHR6tt27b6+OOPA5QWAAAEG9PLjCS1bt1aR44c8U5ffvnlBZdds2aNBg0apNGjR2vz5s0aMGCABgwYoO3btwcwMQAACBZBUWYiIiKUmJjonerWrXvBZZ999ln169dPjzzyiFq2bKkZM2bo+uuv15w5cwKYGAAABIugKDN79uxR/fr11aRJEw0ePFh5eXkXXHbt2rXq06ePz1jfvn21du3aCt/D5XL5TG63u1KyAwCAq+d2u887Vl8q08tMSkqKsrKytHTpUs2dO1f79u3TL37xC506darc5QsKCpSQkOAzlpCQoIKCggrfx+l0KjY21jtlZGRU2mcAAABXJyMjw+c47XQ6L3ndCD/muiS33nqr98/t2rVTSkqKGjVqpLffflujR4+utPfJz8+Xw+Hwztvt9krbNgAAuDrp6elKS0vzzrtcrksuNKaXmZ+rWbOmrrvuOuXm5pb7emJiogoLC33GCgsLlZiYWOF2HQ6HT5kBAADBw263X/GJBtMvM/3c6dOntXfvXtWrV6/c11NTU7VixQqfsWXLlik1NTUQ8QAAQJAxvcw8/PDDWrVqlfbv3681a9bot7/9rcLDwzVo0CBJ0rBhw5Senu5d/sEHH9TSpUv19NNPa9euXZo2bZo2bNig8ePHm/URAACAiUy/zHTw4EENGjRIx48fV1xcnG688UatW7dOcXFxkqS8vDyFhf3Uubp3765FixZp8uTJeuyxx9S8eXMtWbJEbdq0MesjAAAAE9kMwzDMDuFPLpdLsbGxKioq4p4ZAAAs4nKO36ZfZgIAALgalBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBplBkAAGBpppeZjIwMdenSRTExMYqPj9eAAQO0e/fuCtfJysqSzWbzmaKjowOUGAAABBPTy8yqVas0btw4rVu3TsuWLVNpaaluueUWFRcXV7iew+HQkSNHvNOBAwcClBgAAASTCLMDLF261Gc+KytL8fHx2rhxo2666aYLrmez2ZSYmHjJ7+NyuXzm7Xa77Hb75YUFAAB+4Xa75Xa7vfM/P25XxPQzMz9XVFQkSapdu3aFy50+fVqNGjWS0+nU7bffrh07dlS4vNPpVGxsrHfKyMiotMwAAODqZGRk+BynnU7nJa9rMwzD8GO2y+LxePSb3/xGJ0+e1JdffnnB5dauXas9e/aoXbt2Kioq0lNPPaXVq1drx44datCggc+yLpdLsbGxys/Pl8Ph8I5zZgYAgOBR3pkZp9OpoqIin+N3eYKqzNx///365JNP9OWXX55XSipSWlqqli1batCgQZoxY4bPaz+WmUvZGQAAIDhczvHb9HtmfjR+/Hh99NFHWr169WUVGUmKjIxUx44dlZub66d0AAAgWJl+z4xhGBo/frzee+89ff7552rcuPFlb6OsrEzbtm1TvXr1/JAQAAAEM9PPzIwbN06LFi3S+++/r5iYGBUUFEiSYmNjdc0110iShg0bpmuvvdZ70+706dPVrVs3NWvWTCdPntRf//pXHThwQGPGjDHtcwAAAHOYXmbmzp0rSerZs6fPeGZmpkaMGCFJysvLU1jYTyeRTpw4obFjx6qgoEC1atVSp06dtGbNGrVq1SpQsQEAQJAIqhuA/YEbgAEAsJ7LOX6bfs8MAADA1aDMAAAAS6PMIORtzjuhmf/6WsdPuy++MADAcky/ARjwp3c3HdSkf25TSZlHDWpV0/DuSWZHAgBUMsoMQpJhGHr6s280Z+VPD1IsOecxMREAwF+4zISQYxiGnvjwa2+RqVkt0uREAAB/oswgpBiGob98ultZa/ZLkjJ+11a/bBFvbigAgF9RZhBS5nyeq7lf7JUkzRjQRoO6NjQ5EQDA3ygzCBnz/u9bPb3sG0nS//yqpYZ2a2RyIgBAIFBmEBIWrj+gJ/+1U5I0sc91GntTE5MTAQAChTIDy/vnxoOavGS7JOm+Hk31QO9mJicCAAQSZQaW9q+vjuiRf2yVYUgjuifp0X4tZLPZzI4FAAggygwsa8XOQj24eLM8hvT7zg005detKDIAUAVRZmBJX+45pvsXbtI5j6HftK+vjN+1U1gYRQYAqiLKDCwnZ//3Gvv6BpWc8+iWVgl6+vftFU6RAYAqizIDS9maf1IjM3P0Q2mZelwXp+fv7qjIcP4aA0BVxlEAlrHziEvDXs3Wafc5pTSurReHdJI9ItzsWAAAk1FmYAm5R09r6Pz1KvqhVB0b1tT8EV10TRRFBgBAmYEF5B0/oyHz1uvY6RK1ru9Q1siuqmHnB98BAP9BmUFQO3zyB909b50KXGfVPL6G3hidothr+BVsAMBPKDMIWt+dcmvIvPU6eOIHJdWppoVjUlS7etQVb8+QUYnpAADBgjKDoHSiuERD5q3Xt8eKdW3Na7RwbDfFO6LNjgUACEKUGQQd19lSDXs1W7sLTyk+xq6FY1J0bc1rzI4FAAhSlBkElWL3OY3MzNG2Q0WqXT1KC8ekKKlu9avbKM/TA4CQRplB0DhbWqaxr2/QxgMn5IiO0Buju6p5QozZsQAAQY4yg6BQcs6j+xds1Jq9x1U9Klyvjeqq1vVjzY4FALAAygxMd67MowcXb9bK3d8pOjJM80d0UceGtcyOBQCwCMoMTOXxGHrkH1/pk+0FigoP08tDO6tbkzpmxwIAWAhlBqYxDEP/s2S73tt8SOFhNs25u6Nuui7O7FgAAIuhzMAUhmFoxkc79WZ2nmw26W93dtAtrRPNjgUAsCDKDEzx9Gff6NV/75MkzR7YTr9pX9/kRAAAq6LMIOBeWJmrOStzJUkzbm+t33d2mpwIAGBllBkE1Pwv9+mvn+6WJKXfmqyhqUnmBgIAWB5lBgHzZnaeZnz0tSTpwd7NdW+PpiYnAgCEAsoMAuK9zQf12HvbJEn33tREf+zT3OREAIBQERRl5oUXXlBSUpKio6OVkpKi7OzsCpd/5513lJycrOjoaLVt21Yff/xxgJLiSnyy7YgefucrGYY0LLWRJt2aLJuNH0wCAFQO08vMW2+9pbS0NE2dOlWbNm1S+/bt1bdvXx09erTc5desWaNBgwZp9OjR2rx5swYMGKABAwZo+/btAU6OS7Fy11E9sHizyjyG/l+nBpp2W2uKDACgUtkMwzDMDJCSkqIuXbpozpw5kiSPxyOn06kJEyZo0qRJ5y1/5513qri4WB999JF3rFu3burQoYNefPHF85Z3uVyKjY1Vfn6+HA6Hd9xut8tut/vhE+FHa3KPaURWjkrOefTrdvX07F0dFR4W+CKT9vYWvbvpkNJvTeY+HQAIUm63W2632zvvcrnkdDpVVFTkc/wuj6lnZkpKSrRx40b16dPHOxYWFqY+ffpo7dq15a6zdu1an+UlqW/fvhdc/kdOp1OxsbHeKSMj4+o/AC5o44HvNeb1DSo559HNrRL0tzs7mFJkAADWkJGR4XOcdjov/bEdEX7MdVHHjh1TWVmZEhISfMYTEhK0a9euctcpKCgod/mCgoIK36u8MzPwj20HizTi1RydKSnTL5rX1Zy7Oyoy3PQrmgCAIJaenq60tDTv/I9nZi6FqWUmkBwOx0VPU+Hq7Spwaeir63XKfU5dG9fWy0M7yx4RbmommzgjBADB7mpu/zD1n8t169ZVeHi4CgsLfcYLCwuVmFj+7/QkJiZe1vIInG+/O60h87J18kyp2jtr6tURXXRNlLlFBgAQ+kwtM1FRUerUqZNWrFjhHfN4PFqxYoVSU1PLXSc1NdVneUlatmzZBZdHYOR/f0aD563XsdNutarn0Osju6qGvcqc+AMAmMj0o01aWpqGDx+uzp07q2vXrnrmmWdUXFyskSNHSpKGDRuma6+91nvD7oMPPqgePXro6aefVv/+/bV48WJt2LBBL7/8spkfo0orKDqru+et05Gis2oWX0NvjO6q2GqRZscCAFQRppeZO++8U999952mTJmigoICdejQQUuXLvXe5JuXl6ewsJ9OIHXv3l2LFi3S5MmT9dhjj6l58+ZasmSJ2rRpY9ZHqNK+O+XW3fPWKf/7H9SoTjUtHJOiOjW4uRoAEDimP2fG3358zsylfE8dl+fkmRLd9fI67So4pfqx0Xr7vlQ1qFXN7FjneejtrfrnpoM8ZwYALORyjt98XxZX5NTZUg1/NVu7Ck4pLsauhWO7BWWRAQCEPsoMLtuZknMalZWjrQeLVKtapBaOSVHjutXNjgUAqKIoM7gsZ0vLdM/rG5Wz/4RioiP0xugUXZcQY3YsAEAVRpnBJSs559G4hZv0Ze4xVYsKV9bIrmpzbazZsQAAVRxlBpfkXJlHE9/aohW7jsoeEab5w7uoU6NaZscCAIAyg4vzeAz96Z9f6V/bjigy3KaXhnZSatM6ZscCAEASZQYXYRiGHn9/u97ddEjhYTbNuft69WwRb3YsAAC8KDO4IMMwNPNfO7VwfZ5sNul/f99efVvzG1gAgOBCmcEF/W3ZN5r35T5J0qzftdXtHa41OdHVCemnQwJAFUaZQbn+/kWunvs8V5I07bZWurNLQ5MTAQBQPsoMzpP57336y9LdkqRH+yVrxA2NTU50dWw2sxMAAPyJMgMfb+Xk6YkPv5YkPfDLZrq/J79lBAAIbpQZeL2/5ZAmvbtNkjTmxsaaePN1JicCAODiKDOQJC3dXqC0t7fKMKTBKQ31P/1bysb1GQCABVBmoJW7j2rCm5tU5jE08PoGmnF7G4oMAMAyKDNV3Nq9x3XfGxtVWmaof9t6mj2wrcLCKDIAAOugzFRhGw+c0OjXcuQ+51Hv5Hj97c4OigjnrwQAwFo4clVR2w8VaURmts6UlOnGZnX1wuDrFRXBXwcAgPVw9KqCvik8paHz1+vU2XPqklRLLw/rpOjIcLNjAQBwRSgzVcy+Y8W6+5X1OnGmVO0bxOrVEV1ULSrC7FgAAFwxykwVcvDEGQ1+ZZ2OnXYrOTFGr43qqpjoSLNjAQBwVSgzVUSh66zufmW9DhedVdO46lowJkU1q0WZHQsAgKtGmakCjp126+5X1inv+zNqWLuaFo7ppro17GbHAgCgUlBmQtzJMyUaOj9be78rVr3YaC0ck6LE2GizYwEAUGkoMyHs1NlSDc/M0c4jLtWtYdfCMSly1q5mdizTGIbZCQAA/kCZCVE/lJRpdNYGbc0/qVrVIrVwTIqaxNUwOxYAAJWOMhOCzpaW6Z43Nih7//eKsUfo9VEpapEYY3YsAAD8gjITYkrLPBq/aJP+b88xVYsKV9aoLmrbINbsWKbil6YAILRRZkJImcfQxLe2aPnOo4qKCNO8YZ3VqVFts2MBAOBXlJkQ4fEYevSfX+mjr44oMtyml4Z0Uvdmdc2OBQCA31FmQoBhGJr6wQ79Y+NBhYfZ9PygjuqVHG92LAAAAoIyY3GGYSjjk116Y90B2WzS03e0V7829cyOBQBAwFBmLO6Z5Xv08upvJUl//m1bDeh4rcmJAAAILMqMhb20aq+eXbFHkjTl1600qGtDkxMBABB4lBmLen3tfmV8skuS9EjfFhp1Y2OTEwEAYA7Tysz+/fs1evRoNW7cWNdcc42aNm2qqVOnqqSkpML1evbsKZvN5jPdd999AUodHN7OydeU93dIksb3aqZxvZqZnAgAAPNEmPXGu3btksfj0UsvvaRmzZpp+/btGjt2rIqLi/XUU09VuO7YsWM1ffp073y1alXn94Y+2HpYj777lSRp1A2N9dAt15mcCAAAc5lWZvr166d+/fp555s0aaLdu3dr7ty5Fy0z1apVU2Jior8jBp3PdhRo4ltbZBjS3SkN9fivW8pm4/m2AICqLajumSkqKlLt2hd/Yu3ChQtVt25dtWnTRunp6Tpz5sxF13G5XD6T2+2ujMgBs+qb7zR+0WaVeQz9ruO1evL2NhQZAEDIcLvd5x2rL1XQlJnc3Fw9//zzuvfeeytc7u6779aCBQu0cuVKpaen64033tCQIUMuun2n06nY2FjvlJGRUVnR/W7dt8d1z+sbVFLm0a/aJuov/6+dwsIoMgCA0JGRkeFznHY6nZe8rs0wDKMyw0yaNEmzZ8+ucJmdO3cqOTnZO3/o0CH16NFDPXv21Lx58y7r/T7//HP17t1bubm5atq06Xmvu1wuxcbGKj8/Xw6Hwztut9tlt9sv673MsCnvhIbOW6/ikjL9MjleLw7ppKiIoOmglvDIO1v1zsaDerRfsu7vef7fEQCA+dxut89VE5fLJafTqaKiIp/jd3kq/Z6Zhx56SCNGjKhwmSZNmnj/fPjwYfXq1Uvdu3fXyy+/fNnvl5KSIkkXLDM/cjgcF90ZwWbH4SKNeDVbxSVluqFZHf198PUUmatgqFJ7OwCgEl3NSYZKLzNxcXGKi4u7pGUPHTqkXr16qVOnTsrMzFRY2OUfqLds2SJJqlcvtB7hv6fwlIbOz5br7Dl1blRLrwzrrOjIcLNjAQAQdEz7Z/6hQ4fUs2dPNWzYUE899ZS+++47FRQUqKCgwGeZ5ORkZWdnS5L27t2rGTNmaOPGjdq/f78++OADDRs2TDfddJPatWtn1kepdPuPFWvwvPX6vrhE7RrE6tWRXVQtyrQvngEAENRMO0IuW7ZMubm5ys3NVYMGDXxe+/E2ntLSUu3evdv7baWoqCgtX75czzzzjIqLi+V0OjVw4EBNnjw54Pn95eCJMxo8b72OnnKrRUKMXhvZVY7oSLNjWRpf+gKA0GZamRkxYsRF761JSkrSf9+f7HQ6tWrVKj8nM89R11kNmbdeh07+oCZ1q2vBmBTVqh5ldiwAAIIad5MGieOn3Ro8b732Hz8jZ+1rtHBsiuJigv/bVgAAmI0yEwSKzpRq6Pxs7Tl6WomOaC0a0031Yq8xOxYAAJZAmTHZafc5Dc/M1tdHXKpbI0oLx6bIWbvq/NYUAABXizJjoh9KyjQ6K0db8k+qZrVIvTE6RU3japgdCwAAS6HMmMR9rkz3Ltio9fu+V4w9Qq+P6qqW9az1UD8AAIIBZcYEpWUeTVi0Wau/+U7XRIYrc2QXtWtQ0+xYAABYEmUmwMo8htLe3qrPvi5UVESY5g3vrM5JF/+lcAAAUD7KTAB5PIbS3/1KH249rMhwm14ccr1uaFbX7FgAAFgaZSZADMPQEx/u0NsbDirMJj17V0f9MjnB7FgAAFgeZSYADMPQrKW79NraA7LZpKfuaK9ftQ2tH8YEAMAslJkAeG5Frl5a9a0k6ckBbfS76xtcZA0AAHCpKDN+9srqb/W35d9Ikh7/dSsNTmlkciIAAEILZcaP3lh3QDM/3ilJeviW6zT6xsYmJ6ra/us3SwEAIYQy4yf/2HhQjy/ZLkn6Q8+mGv/L5iYnAgAgNFFm/ODDrYf1p39slSSN6J6kR/q2MDlR1WaTzewIAAA/osxUsuVfF2riW1vkMaRBXZ2aelsr2WwcTAEA8BfKTCX6vz3f6Q8LN+mcx9CADvX15IC2FBkAAPyMMlNJ1n97XGNf36CSMo/6tU7UU3e0V3gYRQYAAH+jzFSCLfknNSorR2dLPerZIk7PDeqoiHB2LQAAgcAR9yp9fdilYfPXq7ikTKlN6ujFIZ0UFcFuBQAgUDjqXoXco6c0dP56uc6eU6dGtTRveGdFR4abHQsAgCqFMnOFDhwv1uB563W8uERtrnUoc2QXVbdHmB0LAIAqhzJzhbbkn9TRU25dl1BDr49KkSM60uxIAABUSZxKuEK3d7hW9ohwXd+opmpXjzI7DgAAVRZl5ir0a5NodgQAAKo8LjMBAABLo8wAAABLo8wAAABLo8wAAABLo8wAAABLo8wAAABLo8wg5Nn48XIACGmUGQAAYGmUGQAAYGmUGQAAYGmmlpmkpCTZbDafadasWRWuc/bsWY0bN0516tRRjRo1NHDgQBUWFgYoMQAACDamn5mZPn26jhw54p0mTJhQ4fITJ07Uhx9+qHfeeUerVq3S4cOH9bvf/S5AaQEAQLAx/YcmY2JilJh4aT/YWFRUpPnz52vRokX65S9/KUnKzMxUy5YttW7dOnXr1s2fUQEAQBAy/czMrFmzVKdOHXXs2FF//etfde7cuQsuu3HjRpWWlqpPnz7eseTkZDVs2FBr166t8H1cLpfP5Ha7K+0zAACAq+N2u887Vl8qU8vMAw88oMWLF2vlypW699579ec//1l/+tOfLrh8QUGBoqKiVLNmTZ/xhIQEFRQUVPheTqdTsbGx3ikjI6MyPgIAAKgEGRkZPsdpp9N5yetW+mWmSZMmafbs2RUus3PnTiUnJystLc071q5dO0VFRenee+9VRkaG7HZ7pebKz8+Xw+Hwzlf29gEAwJVLT0/36QUul+uSC02ll5mHHnpII0aMqHCZJk2alDuekpKic+fOaf/+/WrRosV5rycmJqqkpEQnT570OTtTWFh40ftuHA6HT5kBAADBw263X/GJhkovM3FxcYqLi7uidbds2aKwsDDFx8eX+3qnTp0UGRmpFStWaODAgZKk3bt3Ky8vT6mpqVecGQAAWJdp32Zau3at1q9fr169eikmJkZr167VxIkTNWTIENWqVUuSdOjQIfXu3Vuvv/66unbtqtjYWI0ePVppaWmqXbu2HA6HJkyYoNTUVL7JBABAFWVambHb7Vq8eLGmTZsmt9utxo0ba+LEiT7Xy0pLS7V7926dOXPGO/a3v/1NYWFhGjhwoNxut/r27au///3vZnwEAAAQBEwrM9dff73WrVtX4TJJSUkyDMNnLDo6Wi+88IJeeOEFf8YDAAAWYfpzZoBA+XkxBgCEBsoMAACwNMoMQp7NZnYCAIA/UWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWYAAIClUWZQZRiG2QkAAP5AmQEAAJZGmUEVYDM7AADAjygzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzAADA0igzqDIMswMAAPzCtDLzxRdfyGazlTvl5ORccL2ePXuet/x9990XwOSwGpvN7AQAAH+KMOuNu3fvriNHjviMPf7441qxYoU6d+5c4bpjx47V9OnTvfPVqlXzS0YAABD8TCszUVFRSkxM9M6Xlpbq/fff14QJE2S7yD+lq1Wr5rMuAACouoLmnpkPPvhAx48f18iRIy+67MKFC1W3bl21adNG6enpOnPmzEXXcblcPpPb7a6M2AAAoBK43e7zjtWXKmjKzPz589W3b181aNCgwuXuvvtuLViwQCtXrlR6erreeOMNDRky5KLbdzqdio2N9U4ZGRmVFR0AAFyljIwMn+O00+m85HUr/TLTpEmTNHv27AqX2blzp5KTk73zBw8e1Keffqq33377otu/5557vH9u27at6tWrp969e2vv3r1q2rTpBdfLz8+Xw+Hwztvt9ou+FwAACIz09HSlpaV5510u1yUXmkovMw899JBGjBhR4TJNmjTxmc/MzFSdOnX0m9/85rLfLyUlRZKUm5tbYZlxOBw+ZQYAAAQPu91+xScaKr3MxMXFKS4u7pKXNwxDmZmZGjZsmCIjIy/7/bZs2SJJqlev3mWvCwAArM/0e2Y+//xz7du3T2PGjDnvtUOHDik5OVnZ2dmSpL1792rGjBnauHGj9u/frw8++EDDhg3TTTfdpHbt2gU6OgAACAKmfTX7R/Pnz1f37t197qH5UWlpqXbv3u39tlJUVJSWL1+uZ555RsXFxXI6nRo4cKAmT54c6NgAACBImF5mFi1adMHXkpKSZBg/PYTe6XRq1apVgYgFAAAswvTLTAAAAFeDMgMAACyNMgMAACyNMgMAACyNMgMAACyNMoMq47++GAcACCGUGYQ8m9kBAAB+RZkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRpkBAACWRplBlWHIMDsCAMAPKDMIeTab2QkAAP5EmQEAAJZGmQEAAJZGmQEAAJZGmQEAAJZGmQEAAJZGmQEAAJZGmQEAAJZGmQEAAJZGmQEAAJZGmQEAAJZGmbkKbrdb06ZNk9vtNjtKSGM/Bwb7OTDYz4HDvg6MYNjPNsMwQvrX91wul2JjY1VUVCSHw2GZbeMnV7ufJy/ZpgXr8vTHPs31xz7X+SFhaODvc2CwnwOHfR0Y/trPl7Ndv52ZmTlzprp3765q1aqpZs2a5S6Tl5en/v37q1q1aoqPj9cjjzyic+fOVbjd77//XoMHD5bD4VDNmjU1evRonT592g+fAAAAWIHfykxJSYnuuOMO3X///eW+XlZWpv79+6ukpERr1qzRa6+9pqysLE2ZMqXC7Q4ePFg7duzQsmXL9NFHH2n16tW65557/PERAACABUT4a8NPPPGEJCkrK6vc1z/77DN9/fXXWr58uRISEtShQwfNmDFDjz76qKZNm6aoqKjz1tm5c6eWLl2qnJwcde7cWZL0/PPP61e/+pWeeuop1a9f/7x1fryKdujQIblcLu+43W6X3W6/qs/44/b+e7uofFe7n91nTsvjPqMfTp/mv1UF+PscGOznwGFfB0Zl7We32+1z382pU6ck/XQcr5DhZ5mZmUZsbOx5448//rjRvn17n7Fvv/3WkGRs2rSp3G3Nnz/fqFmzps9YaWmpER4ebrz77rvlrpOfn29IYmJiYmJiYrLglJ+ff9Gu4bczMxdTUFCghIQEn7Ef5wsKCi64Tnx8vM9YRESEateufcF16tevr7179yoyMlI2m807XhlnZgAAQOX4+ZkZwzBUWlpa7lWXn7usMjNp0iTNnj27wmV27typ5OTky9msX4WFhalJkyZmxwAAAH5yWWXmoYce0ogRIypc5lKLQ2JiorKzs33GCgsLva9daJ2jR4/6jJ07d07ff//9BdcBAACh7bLKTFxcnOLi4irljVNTUzVz5kwdPXrUe+lo2bJlcjgcatWq1QXXOXnypDZu3KhOnTpJkj7//HN5PB6lpKRUSi4AAGAtfvtqdl5enrZs2aK8vDyVlZVpy5Yt2rJli/eZMLfccotatWqloUOHauvWrfr00081efJkjRs3znsvS3Z2tpKTk3Xo0CFJUsuWLdWvXz+NHTtW2dnZ+ve//63x48frrrvuuqRragAAIPT47QnAI0aM0GuvvXbe+MqVK9WzZ09J0oEDB3T//ffriy++UPXq1TV8+HDNmjVLERH/OWH0xRdfqFevXtq3b5+SkpIk/eeheePHj9eHH36osLAwDRw4UM8995xq1Kjhj48BAACC3UW/7wQff//73422bdsaMTExRkxMjNGtWzfj448/NjtWyMvIyDAkGQ8++KDZUULO1KlTz/sqZIsWLcyOFZIOHjxoDB482Khdu7YRHR1ttGnTxsjJyTE7Vkhp1KhRuV/v/cMf/mB2tJBy7tw5Y/LkyUZSUpIRHR1tNGnSxJg+fbrh8XhMyWPaV7OtqkGDBpo1a5aaN28uwzD02muv6fbbb9fmzZvVunVrs+OFpJycHL300ktq166d2VFCVuvWrbV8+XLv/I9nR1F5Tpw4oRtuuEG9evXSJ598ori4OO3Zs0e1atUyO1pIycnJUVlZmXd++/btuvnmm3XHHXeYmCr0zJ49W3PnztVrr72m1q1ba8OGDRo5cqRiY2P1wAMPBDwP/491mW677Taf+ZkzZ2ru3Llat24dZcYPTp8+rcGDB+uVV17Rk08+aXackBUREcE3Av1s9uzZcjqdyszM9I41btzYxESh6edfUpk1a5aaNm2qHj16mJQoNK1Zs0a33367+vfvL0lKSkrSm2++ed63lAPFbzcAVwVlZWVavHixiouLlZqaanackDRu3Dj1799fffr0MTtKSNuzZ4/q16+vJk2aaPDgwcrLyzM7Usj54IMP1LlzZ91xxx2Kj49Xx44d9corr5gdK6SVlJRowYIFGjVqlM9DU3H1unfvrhUrVuibb76RJG3dulVffvmlbr31VlPycGbmCmzbtk2pqak6e/asatSooffee++CXyfHlVu8eLE2bdqknJwcs6OEtJSUFGVlZalFixY6cuSInnjiCf3iF7/Q9u3bFRMTY3a8kPHtt99q7ty5SktL02OPPaacnBw98MADioqK0vDhw82OF5KWLFmikydPXvT5aLh8kyZNksvlUnJyssLDw1VWVqaZM2dq8ODB5gQy5U4di3O73caePXuMDRs2GJMmTTLq1q1r7Nixw+xYISUvL8+Ij483tm7d6h3r0aMHNwAHwIkTJwyHw2HMmzfP7CghJTIy0khNTfUZmzBhgtGtWzeTEoW+W265xfj1r39tdoyQ9OabbxoNGjQw3nzzTeOrr74yXn/9daN27dpGVlaWKXk4M3MFoqKi1KxZM0lSp06dlJOTo2effVYvvfSSyclCx8aNG3X06FFdf/313rGysjKtXr1ac+bMkdvtVnh4uIkJQ1fNmjV13XXXKTc31+woIaVevXrnncFt2bKl/vnPf5qUKLQdOHBAy5cv17vvvmt2lJD0yCOPaNKkSbrrrrskSW3bttWBAweUkZFhyplGykwl8Hg8Pj+OhavXu3dvbdu2zWds5MiRSk5O1qOPPkqR8aPTp09r7969Gjp0qNlRQsoNN9yg3bt3+4x98803atSokUmJQltmZqbi4+O9N6iicp05c0ZhYb633YaHh8vj8ZiShzJzmdLT03XrrbeqYcOGOnXqlBYtWqQvvvhCn376qdnRQkpMTIzatGnjM1a9enXVqVPnvHFcnYcffli33XabGjVqpMOHD2vq1KkKDw/XoEGDzI4WUiZOnKju3bvrz3/+s37/+98rOztbL7/8sl5++WWzo4Ucj8ejzMxMDR8+nMcM+Mltt92mmTNnqmHDhmrdurU2b96s//3f/9WoUaNMycN/5ct09OhRDRs2TEeOHFFsbKzatWunTz/9VDfffLPZ0YArcvDgQQ0aNEjHjx9XXFycbrzxRq1bt67SfocN/9GlSxe99957Sk9P1/Tp09W4cWM988wz5t0wGcKWL1+uvLw80w6sVcHzzz+vxx9/XH/4wx909OhR1a9fX/fee6+mTJliSh6//ZwBAABAIPCcGQAAYGmUGQAAYGmUGQAAYGmUGQAAYGmUGQAAYGmUGQAAYGmUGQAAYGmUGQAAYGmUGQAAYGmUGQAAYGmUGQAAYGn/H82BluwkNa5nAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(np.log10(t),np.log10(Fnu))\n",
    "plt.ylim(-10,10.)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "base",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.3"
  },
  "vscode": {
   "interpreter": {
    "hash": "5d5727baf96e20f3ff6d9eb9958f1aee38dd4ea4d6cb478932055b20d56a57e6"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}