{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "cab47263",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-06-19T08:14:07.882964Z",
     "iopub.status.busy": "2026-06-19T08:14:07.882726Z",
     "iopub.status.idle": "2026-06-19T08:14:13.372363Z",
     "shell.execute_reply": "2026-06-19T08:14:13.371449Z"
    },
    "tags": [
     "remove-cell"
    ]
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "An NVIDIA GPU may be present on this machine, but a CUDA-enabled jaxlib is not installed. Falling back to cpu.\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import brainmass\n",
    "import brainstate\n",
    "import brainunit as u\n",
    "import jax.numpy as jnp\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "brainstate.random.seed(0)\n",
    "brainstate.environ.set(dt=0.1 * u.ms)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e97cabec",
   "metadata": {},
   "source": [
    "# Epileptor (Seizure Dynamics)\n",
    "\n",
    "The **Epileptor** is a six-dimensional model of epileptic seizure dynamics. It couples a fast oscillatory subsystem $(x_1, y_1)$, a second (spike-wave) subsystem $(x_2, y_2, g)$, and a **slow permittivity variable** $z$ that acts on a much slower timescale. The slow variable autonomously ramps seizure-like discharges *on* and *off*, reproducing the characteristic onset/offset structure of focal seizures. The excitability parameter $x_0$ sets the epileptogenicity of the node. The recorded signal is the LFP proxy $\\text{lfp} = x_2 - x_1$.\n",
    "\n",
    "**Reference:** Jirsa, Stacey, Quilichini, Ivanov & Bernard (2014), *On the nature of seizure dynamics*, Brain 137(8):2210-2230."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1ed74451",
   "metadata": {},
   "source": [
    "## Build the model\n",
    "\n",
    "`x0 = -1.6` makes the node epileptogenic (it will seize autonomously). A bounded but long run is needed for the slow `z` variable to drive a full onset/offset cycle."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "badb82e7",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-06-19T08:14:13.375721Z",
     "iopub.status.busy": "2026-06-19T08:14:13.375283Z",
     "iopub.status.idle": "2026-06-19T08:14:13.400643Z",
     "shell.execute_reply": "2026-06-19T08:14:13.399827Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "EpileptorStep(\n",
       "  in_size=(1,),\n",
       "  out_size=(1,),\n",
       "  a=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(1., dtype=float32)\n",
       "  ),\n",
       "  b=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(3., dtype=float32)\n",
       "  ),\n",
       "  c=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(1., dtype=float32)\n",
       "  ),\n",
       "  d=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(5., dtype=float32)\n",
       "  ),\n",
       "  r=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(0.00035, dtype=float32)\n",
       "  ),\n",
       "  x0=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(-1.6, dtype=float32)\n",
       "  ),\n",
       "  Iext=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(3.1, dtype=float32)\n",
       "  ),\n",
       "  slope=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(0., dtype=float32)\n",
       "  ),\n",
       "  Iext2=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(0.45, dtype=float32)\n",
       "  ),\n",
       "  tau=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(10., dtype=float32)\n",
       "  ),\n",
       "  aa=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(6., dtype=float32)\n",
       "  ),\n",
       "  bb=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(2., dtype=float32)\n",
       "  ),\n",
       "  Kvf=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(0., dtype=float32)\n",
       "  ),\n",
       "  Kf=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(0., dtype=float32)\n",
       "  ),\n",
       "  Ks=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(0., dtype=float32)\n",
       "  ),\n",
       "  tt=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(1., dtype=float32)\n",
       "  ),\n",
       "  modification=Const(\n",
       "    fit=False,\n",
       "    t=IdentityT(),\n",
       "    reg=None,\n",
       "    val=Array(0., dtype=float32)\n",
       "  ),\n",
       "  init_x1=Constant(value=-1.5),\n",
       "  init_y1=Constant(value=-10.0),\n",
       "  init_z=Constant(value=3.5),\n",
       "  init_x2=Constant(value=-1.0),\n",
       "  init_y2=Constant(value=0.0),\n",
       "  init_g=Constant(value=0.0),\n",
       "  method=exp_euler\n",
       ")"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "node = brainmass.EpileptorStep(in_size=1, x0=-1.6)\n",
    "node"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "92692db0",
   "metadata": {},
   "source": [
    "## Run a simulation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "aa8ad6c4",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-06-19T08:14:13.403214Z",
     "iopub.status.busy": "2026-06-19T08:14:13.402881Z",
     "iopub.status.idle": "2026-06-19T08:14:13.664195Z",
     "shell.execute_reply": "2026-06-19T08:14:13.663358Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(30000, 1)"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sim = brainmass.Simulator(node, dt=0.1 * u.ms)\n",
    "res = sim.run(3000. * u.ms,\n",
    "              monitors={'lfp': lambda m: m.lfp(), 'z': 'z'})\n",
    "res['lfp'].shape"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d48833df",
   "metadata": {},
   "source": [
    "## Visualize\n",
    "\n",
    "Top: the LFP proxy shows quiescence, an abrupt seizure onset (fast discharges), and offset. Bottom: the slow permittivity variable `z` is the hidden driver — its slow ramp times the onset and offset."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "eca79d46",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-06-19T08:14:13.666426Z",
     "iopub.status.busy": "2026-06-19T08:14:13.666166Z",
     "iopub.status.idle": "2026-06-19T08:14:13.840520Z",
     "shell.execute_reply": "2026-06-19T08:14:13.839490Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3kAAAJOCAYAAAAK+M50AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAA8jlJREFUeJzs3Xd8U1X/B/DPTdKku6V0UShlb5C9ZAqKDBUQRUUBRUXcoqg8j4Kbx8c9cf5AcQEO9HEwXKgsAZnK3rstdNHd5P7+SJPmJjc3J6vp+Lxfr77a3Jx778loe7/5nvM9kizLMoiIiIiIiKhO0IW6A0RERERERBQ4DPKIiIiIiIjqEAZ5REREREREdQiDPCIiIiIiojqEQR4REREREVEdwiCPiIiIiIioDmGQR0REREREVIcwyCMiIiIiIqpDGOQRERERERHVIQzyiIjqgcceewySJCm2NWvWDFOnTg1Nh6hWOnbsGMLDw7FmzZpQd6VaLV++HNHR0cjKygp1V4iIhDDIIyKqYRYuXAhJktx+rV+/PtRdFPb999/jscceq9ZzDhkyBJ06ddJsYwt61b7eeustezvH7TqdDmlpabjkkkvw66+/BvlR1ExPPPEE+vTpgwsvvNCr/SwWCxYuXIjLL78c6enpiIqKQqdOnfDUU0+hpKQkSL31bOXKlZg2bRo6deoEvV6PZs2aqba79NJL0apVK8ybN696O0hE5CNDqDtARETqnnjiCTRv3txle6tWrbw+1iOPPIKHH344EN3yyvfff4833nij2gM9UfPnz0d0dLRiW58+fRS3L774YkyePBmyLOPQoUN48803cdFFF+G7777DyJEjq7O7IZWVlYUPPvgAH3zwgdf7FhUV4cYbb0Tfvn1x2223ITk5GevWrcPcuXPx008/4eeff3bJNFeHTz75BIsXL0b37t2Rlpam2Xb69Ol44IEH8PjjjyMmJqaaekhE5BsGeURENdTIkSPRs2fPgBzLYDDAYKgbf/ItFgvKysoQHh7u97EmTJiAxMREzTZt2rTB9ddfb789btw4dOnSBS+//HLQgrySkhIYjUbodDVnwM1HH30Eg8GAyy67zOt9jUYj1qxZg/79+9u33XLLLWjWrJk90Bs+fHgguyvkmWeewbvvvouwsDCMGTMGO3fudNv2yiuvxF133YWlS5fipptuqsZeEhF5r+b89yAiIq8cPnwYkiTh+eefx0svvYSMjAxERERg8ODBLheranPy1OTm5uLee+9Feno6TCYTWrVqhWeffRYWi8Xr806dOhVvvPEGAOWwR5vCwkLcf//99nO1bdsWzz//PGRZVvRJkiTceeed+Pjjj9GxY0eYTCYsX77cp+csEDp37ozExEQcOnRIs51t2OjmzZvRv39/REREoHnz5orhoADw66+/QpIkfPbZZ3jkkUfQuHFjREZGIj8/HwCwdOlS9OjRAxEREUhMTMT111+PEydO2PefO3cudDodfvrpJ8Vxb731VhiNRmzbtg3nz59HVFQU7rnnHpd+Hj9+HHq93uNQxGXLlqFPnz6KzOeuXbsQERGByZMnK9r+8ccf0Ov1eOihhwBYgzzHAM9m3Lhx9uMEys8//wydToc5c+Yotn/yySeQJAnz58+3b0tLS0NYWJjQcZOTk9GlSxd8/fXXAesrEVGw1I2PdYmI6qC8vDxkZ2crtkmShIYNGyq2ffjhhygoKMAdd9yBkpISvPLKK7jooouwY8cOpKSkCJ+vqKgIgwcPxokTJzB9+nQ0bdoUa9euxezZs3Hq1Cm8/PLLXp13+vTpOHnyJFatWoVFixYp9pVlGZdffjl++eUXTJs2DV27dsWKFSswa9YsnDhxAi+99JKi/c8//4wlS5bgzjvvRGJiotu5U946d+6c4rZer0eDBg0098nJyUFOTo7QsNmcnByMGjUKV199Na699losWbIEM2bMgNFodMkGPfnkkzAajXjggQdQWloKo9GIhQsX4sYbb0SvXr0wb948nDlzBq+88grWrFmDLVu2ID4+Ho888gj+97//Ydq0adixYwdiYmKwYsUKvPvuu3jyySdxwQUXALAGVIsXL8aLL74IvV5vP++nn34KWZYxadIkt4+jvLwcGzduxIwZMxTb27dvjyeffBKzZs3ChAkTcPnll6OwsBBTp05Fu3bt8MQTT2g+P6dPnwYAj9lUb1x00UW4/fbbMW/ePIwdOxbdu3fHqVOncNddd2H48OG47bbbfD52jx49sGzZsoD1lYgoaGQiIqpRFixYIANQ/TKZTPZ2hw4dkgHIERER8vHjx+3bN2zYIAOQ77vvPvu2uXPnys5/8jMyMuQpU6bYbz/55JNyVFSUvHfvXkW7hx9+WNbr9fLRo0e9Pu8dd9zhcl5ZluVly5bJAOSnnnpKsX3ChAmyJEny/v377dsAyDqdTv777781nzebwYMHyx07dtRsY3s+nL8yMjIU7QDI06ZNk7OysuTMzEx5w4YN8rBhw2QA8gsvvOCxH87tSktL5a5du8rJyclyWVmZLMuy/Msvv8gA5BYtWshFRUX2tmVlZXJycrLcqVMnubi42L7922+/lQHIc+bMsW/bsWOHbDQa5ZtvvlnOycmRGzduLPfs2VMuLy+3t1mxYoUMQP7hhx8U/ezSpYs8ePBgzceyf/9+GYD82muvudxnNpvlAQMGyCkpKXJ2drZ8xx13yAaDQd64caPmMWVZlocPHy7HxsbKOTk5Htt6o7CwUG7VqpXcsWNHuaSkRB49erQcGxsrHzlyxO0+o0ePdnn9nT3zzDMyAPnMmTMB7S8RUaBxuCYRUQ31xhtvYNWqVYqvH374waXd2LFj0bhxY/vt3r17o0+fPvj++++9Ot/SpUsxcOBANGjQANnZ2fav4cOHw2w247fffgvYeb///nvo9Xrcfffdiu33338/ZFl2eZyDBw9Ghw4dvHo8Ir744gvF8/vxxx+7tHn//feRlJSE5ORk9OnTB2vWrMHMmTNx7733ejy+wWDA9OnT7beNRiOmT5+OzMxMbN68WdF2ypQpiIiIsN/etGkTMjMzcfvttyvmH44ePRrt2rXDd999Z9/WqVMnPP7443jvvfcwYsQIZGdn44MPPlDMwxw+fDjS0tIUj3Hnzp3Yvn27Ys6hmrNnzwKAapZTp9Nh4cKFOH/+PEaOHIk333wTs2fP9jif9JlnnsGPP/6I//znP4iPj9ds663IyEgsXLgQu3btwqBBg/Ddd9/hpZdeQtOmTf06ru3xO2fYiYhqGg7XJCKqoXr37i1UeKV169Yu29q0aYMlS5Z4db59+/Zh+/btSEpKUr0/MzMzYOc9cuQI0tLSXKoUtm/f3n6/I7Uqo4EwaNAgj0MFr7jiCtx5552QJAkxMTHo2LEjoqKihI6flpbm0rZNmzYArHMb+/bta9/u/Bhtz0Hbtm1djtuuXTv88ccfim2zZs3CZ599hj///BPPPPOMS1Cs0+kwadIkzJ8/H0VFRYiMjMTHH3+M8PBwXHXVVUKPR3aaL2nTsmVLPPbYY5g1axY6deqERx99VPM4ixcvxiOPPIJp06a5DAFVc+7cOZSVldlvR0REIC4uTnOfCy+8EDNmzMAbb7yBESNGBKRYiu3xh6ISKBGRNxjkERERAGvVyosvvhgPPvig6v224CQUHDNc1a1JkybVUvnR38d48OBB7Nu3DwCwY8cO1TaTJ0/Gc889h2XLluHaa6/FJ598gjFjxngMmGzzQHNycty2WblyJQDg5MmTOHv2LFJTU1XbrVq1CpMnT8bo0aNditC4M378eKxevdp+e8qUKVi4cKHmPqWlpfb1DA8cOGAPbP1he/yBnENIRBQMDPKIiGo524W9o71793pdnKRly5Y4f/68cEAjcl53GY+MjAz8+OOPKCgoUGTzdu/ebb+/Ljh58iQKCwsV2by9e/cCgMfXx/Yc7NmzBxdddJHivj179iieI4vFgqlTpyI2Nhb33nsvnnnmGUyYMAHjx49X7NepUyd069YNH3/8MZo0aYKjR4/itdde8/g4mjZtioiICLcVRd966y2sWrUKTz/9NObNm4fp06erVqHcsGEDxo0bh549e2LJkiXCy3q88MILigDT05p2gLXq6K5du/D888/joYcewsMPP4xXX31V6HzuHDp0CImJiW6z3URENQXn5BER1XLLli1TlNT/888/sWHDBq/XcLv66quxbt06rFixwuW+3NxcVFRUeH1eW3CTm5ur2HfUqFEwm814/fXXFdtfeuklSJJUZxYZr6iowNtvv22/XVZWhrfffhtJSUno0aOH5r49e/ZEcnIy3nrrLZSWltq3//DDD9i1axdGjx5t3/biiy9i7dq1eOedd/Dkk0+if//+mDFjhurcsRtuuAErV67Eyy+/jIYNGwo912FhYejZsyc2bdrkct+hQ4cwa9YsXHnllfjXv/6F559/Ht988w0+/PBDRTtbn5s1a4Zvv/3Wq8xljx49MHz4cPuXp/mZGzZswPPPP497770X999/P2bNmoXXX39dkQ30xebNm9GvXz+/jkFEVB2YySMiqqF++OEHe2bLUf/+/dGiRQv77VatWmHAgAGYMWMGSktL7Rfv7oZdujNr1ix88803GDNmDKZOnYoePXqgsLAQO3bswOeff47Dhw8rhqmJnNcWyNx9990YMWIE9Ho9rrnmGlx22WUYOnQo/v3vf+Pw4cO44IILsHLlSnz99de499570bJlS2+fLoWsrCw89dRTLtubN2+uuVRAoKWlpeHZZ5/F4cOH0aZNGyxevBhbt27FO++843F9trCwMDz77LO48cYbMXjwYFx77bX2JRSaNWuG++67D4A1eHr00UcxdepU+0LlCxcuRNeuXXH77be7zJG87rrr8OCDD+Krr77CjBkzhNeJu+KKK/Dvf/8b+fn5iI2NBWCdo3bTTTchIiLCvv7c9OnT8cUXX+Cee+6xF3spKCjAiBEjkJOTg1mzZimKxgDWLHKggqeSkhJMmTIFrVu3xtNPPw0AePzxx/G///0PN954I3bs2GH/8GH79u345ptvAAD79+9HXl6e/X1zwQUXKBZ+z8zMxPbt23HHHXcEpJ9EREEV2uKeRETkTGsJBQDyggULZFmuWsrgueeek1944QU5PT1dNplM8sCBA+Vt27YpjimyhIIsy3JBQYE8e/ZsuVWrVrLRaJQTExPl/v37y88//7y95L83562oqJDvuusuOSkpSZYkSdGHgoIC+b777pPT0tLksLAwuXXr1vJzzz0nWywWxTEAyHfccYfw82dbukDta9iwYYrnIysrS/NY3p7buR8dO3aUN23aJPfr108ODw+XMzIy5Ndff13RzraEwtKlS1WPs3jxYrlbt26yyWSSExIS5EmTJtmXrqioqJB79eolN2nSRM7NzVXs98orr8gA5MWLF7scc9SoUTIAee3atcKP58yZM7LBYJAXLVrkco4vvvhC0fbo0aNybGysPGrUKFmWq94z7r6c34f+uO+++2S9Xi9v2LBBsX3Tpk2ywWCQZ8yYYd+m9bvm3Kf58+fLkZGRcn5+fsD6SkQULJIsuymVRURENdrhw4fRvHlzPPfcc3jggQfq/HlrmyFDhiA7Oxs7d+4MdVdcjBs3Djt27MD+/fu92m/atGnYu3cvfv/99yD1rObq1q0bhgwZgpdeeinUXSEi8ohz8oiIiOqRU6dO4bvvvsMNN9zg9b5z587Fxo0bsWbNmiD0rOZavnw59u3bh9mzZ4e6K0REQjgnj4iIqB44dOgQ1qxZg/feew9hYWGKRdpFNW3aFCUlJUHoXc126aWX4vz586HuBhGRMGbyiIiI6oHVq1fjhhtuwKFDh/DBBx+4XceOiIhqP87JIyIiIiIiqkOYySMiIiIiIqpDGOQRERERERHVISy8osJiseDkyZOIiYmBJEmh7g4RERERERFkWUZBQQHS0tKg07nP1zHIU3Hy5Emkp6eHuhtEREREREQujh07hiZNmri9n0GeipiYGADWJy82NjbEvSEiIiIiIgLy8/ORnp5uj1fcYZCnwjZEMzY2lkEeERERERHVKJ6mlLHwChERERERUR3CII+IiIiIiKgOYZBHRERERERUhzDIIyIiIiIiqkMY5BEREREREdUhDPKIiIiIiIjqEAZ5tczqvVl489f9kGU51F0hIiIiIqIaiOvk1TJT/u9PAED71FgMbZcc4t4QEREREVFNw0xeLbXrdH6ou0BERERERDUQg7xaqrTcEuouEBERERFRDcQgj4iIiIiIqA5hkEdERERERFSHMMirpVhbk4iIiIiI1DDIq0W4bAIREREREXnCIK8WqbAwyCMiIiIiIm0M8mqRCrNDkMesHhERERERqWCQV4tUWLhsAhERERERaWOQV4soMnlEREREREQqGOTVIuUOmTwzh2sSEREREZGKOhfkzZs3D7169UJMTAySk5MxduxY7NmzJ9TdCgizQ+EVxnhERERERKSmzgV5q1evxh133IH169dj1apVKC8vxyWXXILCwsJQd81vjsM1mckjIiIiIiI1hlB3INCWL1+uuL1w4UIkJydj8+bNGDRoUIh6FRjl5qrhmozxiIiIiIhITZ0L8pzl5eUBABISEty2KS0tRWlpqf12fn5+0PvlC8fhmhaumUdERERERCrq3HBNRxaLBffeey8uvPBCdOrUyW27efPmIS4uzv6Vnp5ejb0UV+4wXJMxHhERERERqanTQd4dd9yBnTt34rPPPtNsN3v2bOTl5dm/jh07Vk099I7jOnkWjtckIiIiIiIVdXa45p133olvv/0Wv/32G5o0aaLZ1mQywWQyVVPPfFfhOFyTQR4REREREamoc0GeLMu466678NVXX+HXX39F8+bNQ92lgKkwM8gjIiIiIiJtdS7Iu+OOO/DJJ5/g66+/RkxMDE6fPg0AiIuLQ0RERIh75x/lcM0QdoSoFln592k0bRiJdqmxoe4KERERUbWoc3Py5s+fj7y8PAwZMgSNGjWyfy1evDjUXfObcjF0RnlEnuw4nodbF23GpS//HuquEBEREVWbOpfJq8vBj2K4pkWjIREBAA6fLQx1F4iIiIiqXZ3L5NVljoVXzHU4mCUKlDB91Z+4uvwBEBEREZEjBnm1SIWZSygQeUMnVf3suM4kERERUV3GIK8WqVDMyQthR4hqCUmqivLKzRzjTERERPUDg7xahIuhE/mOQR4RERHVFwzyahHlOnkh7AhRLeE4D6+MQR4RERHVEwzyahHH4ZoWRnlEHjkuO8I5eURERFRfMMirRRRBHodrEnlU7hjkVTCTR0RERPUDg7xahNU1ibxjdpjHWsHFJYmIiKieYJBXi5gtnJNH5A3HIZqckkdERET1BYO8WkS5hAKjPCJPzBziTERERPUQg7xaxHG4ppmpPCKP+DtDRERE9RGDvFqkpNzhgpXXq0QelZuZySMiIqL6h0FeLVJSbrb/zOGaRJ5xHisRERHVRwzyapGSCscgL4QdIQqxXafysfVYrsd25RYO1yQiIqL6xxDqDpA4x+GaHHpG9ZUsyxj5yu8AgE2PDEditMltW7OZxYqIiIio/mEmrxZxHK7JII/qK8cqs8fOFWm2dVwMnZk8IiIiqi8Y5NUiikwe1/yieqrCMTvnsa1j9jtIHSIiIiKqYRjk1SLF5RX2n5nJo/rKcZ6dpyGYXCePiIiI6iMGebVITmG5/WdesFJ9VV7hWEzFQ1szh2sSERFR/cMgrxbJKSqz/8zrVaqvHOfklXuI8iosLFZERERE9Q+DvFpClmWcKyxT3Caqj8oqxJdFqOBwTSIiIqqHGOTVEsXlZpRWsIgEkWPgZvYQuCkKr7BYEREREdUTDPJqCccsHsCsBNVfysBNPJPnKSAkIiIiqiu4GHot4Vh0BWARCQq9I2cL8drP+zF9UAu0TomptvOWmb0YrsnF0Imojtl7pgBbjubAbAHS4sPRLjUWqXHhoe4WEdUwdTLI++233/Dcc89h8+bNOHXqFL766iuMHTs21N3yy7kiZSaP16sUajOXbMPmIzn4bvsp7Hry0mo7b4UXFTMdC694qsRJRFSTlVaY8eDn2/H11pMu9zWOj0CvZg3Qq3kCejdLQKvkaEiSFIJeElFNUSeDvMLCQlxwwQW46aabMH78+FB3JyByKodr6iTrfDwO16RQ23osF4B1vmh1cqyo6XlOHodrElHd8Oiynfh660noJKBfy4YIN+hx5FwRDmUX4kRuMU5sLcayygCwQWQYujSJxwVN4tC58ntyLLN9RPVJnQzyRo4ciZEjR4a6GwFlm5PXMNqErIJSj0FeSbkZB7LOo7DUjAaRYWjSIBIRRn11dJXqiTC9JDxs+NM/j+Ln3Zl47dpuCA/z733ozdp3ZV7M3yMiqqm2HsvFkk3HIUnAwht7Y1CbJPt9haUV+OtoDjYezsHGQ+ew5VgOcorKsXpvFlbvzbK3S4k1oW1qLFonR1u/UqLRKikGcZFhoXhIRBRkdTLI81ZpaSlKS0vtt/Pz80PYG3W2NfIaRhmRVVDqdrhmZkEJXly5F8u2nkBJedUFrk4CWifHoEuTOHRr2gDdM+LROjkGeh2Hc5BvwvQ6xXtMy+wvdwAAPlx3GLcOaunXeb1Z+85xuQUZ2m1P55Xg+x2ncFXPJogJ50UPEdUcH6w9DAAY17WxIsADgCiTAQNbJ2Fga+v2sgoL/jmVjx3Hc7HteB52HM/DvswCnMkvxZn8LPzmEPgBQGK0CU0TIpCeEIn0BpFIT4io/B6JRnHhMOhZo4+oNmKQB2DevHl4/PHHQ90NTbZMXmK0CUCB6sXtPyfzMXXBn8gssAasDSLDEBcRhnOFZcgvqcCeMwXYc6YASzcfBwBEmwzomh6P7k3j0S2jAbqnN+AneiTM4MMHBDlF5Z4beeA4XNNxOKYa5Zp62se97t31OJhdiK3HcvHqtd386iMRUaCUlJuxfOdpAMDk/s08tjcadOiaHo+u6fG4oXJbYWkFdp/Ox74z57Ev0/q1/0wBTuaVIPt8KbLPl+Kvo7kux5Ik63VHSqwJKTHhSI4Nt/5c+T05JhwNo41oEGn0e5QGEQUWgzwAs2fPxsyZM+238/PzkZ6eHsIeucqtvDhOiDICcF0n73ReiT3Aa50cjWfGd0bPjAb2iddn8kuw/Xgeth3LxV9Hc7DtWC7Ol1bgj/3Z+GN/tv04LZOi0L1pA3RJj0eb5Gi0SYlBg8pzEjnS+TCp39Me324/idN5Jbh5YAu3bRyHa3rM5JnFs34HswsBAKv+OeOhl0RE1WfT4RwUl5uRHGPCBU3ifDpGlMmAHhkJ6JGRoNheUFKOw9lFOJZThGPnrN+PnivG8XNFOJ5TjDKzBVkFpcgqKMVOaI9yijTqkRBlREKUNehrGGVEg8rbCVFGxEWEITY8DLERhsrvYYgJNyCMmUKioGCQB8BkMsFkMoW6G5pyi62ZPFuQ5zgXSZZlPPjFdmQWlKJtSgyWzuiHWKfhZimx4bi4Qzgu7pBi33/P6QL8dTQHfx3NwZajuTiUXYgDWdYvW7YPABKjjWidHIM2KdFonhiFjIZRSE+IRJMGEfzkrh7zpXCbp8Dwzk+2AACGtE1Cq2T1ZRmU1TW1z6cYrilYeKWcZTiJqAaxfRA7sHVSwCtmxoSHoXOTOHRWCR4tFhlnC8twJr+k8qsUZ/JLkFng+HMpcgrLUGGRUVRmRlFZMY7nFHvVh0ij3h78xYSHITbcgNjKgDAm3IAokwHRJgMijXpEm6y3o0x663dj1W2TgdcjRI4Y5NUSzpk8xwvW5TtP47e9WTAadJh/fXeXAE+NXiehQ1osOqTF4vq+GQCsQ0K3VAZ9u04VYO+ZAhzPKUb2+TJknz+LdQfPKo4hSUBqbDjSEyLRNCESjeMj0CguHKmVX41iIxAbYWAZ5zrKl9dVdJdj54rdB3mKZRG0AzJvhmva1NTKtfkl5Xhp1V5c0bUxuqbHh7o7RFRNdp7IAwD0bNagWs+r00lIijEhKcaETo3dZxBlWUZ+SQVyCstwrqgM585bv+cUluGcw1d+STnyiytQUFKO/JIKnC+tAIDK4NCM036WQwjTS4g0GioDQb3iZ1swGGnSIzLMGjBGGPWIrPyKMFZuC7NtM9jvZ6aRaqs6GeSdP38e+/fvt98+dOgQtm7dioSEBDRt2jSEPfOdLchr4DRcU5ZlvPLTPgDA9EEt0CIp2udzJEQZMax9Coa1T7FvKyytwIGs89h75jz2nSnAkbNFOHKuCEfPFqKwzIxTeSU4lVeCPw+dUz1mRJgejeLCkRIbjkZx4UiKNSExyoTEGCMaRpnQMNqIpGgTGkQZ+Ye0lvEldBcNDLWyad5U1yytEB+uGUqv/rQPBSXl+PfoDm7bvPHzfixYcxgL1hzG4f+Mdttu+/Fc/L4vG9MHtWDBBKJaTpZl/H3SGuR1TIsNcW/USZKEuAhrDYBmiBLer8JswfnSCuQXV1QGgOXW7yUVlT9bvxeVVaCw1IzzpRUoKqvA+VIzCu0/V9gLgJWbZeQVlyOv2P+5345swaMiMAyrCgKrgkWDQ5DoEDga9YgMUwaOtvZGA/9GU/DUySBv06ZNGDp0qP22bb7dlClTsHDhwhD1yj+2P1oJkbYgz3rB+uueLOw+XYAoox43D3A/j8lXUSYDujSJR5cm8YrtsizjXGEZjp4rsn6dLcLJvBKczivGqTzr0I6conIUl5txMLvQPt9JS4PIMDSMNiEx2mj9Xjmev0GkEfGRYWgQWfVzQpQRkUY9s4Qh5MucPNFaLVoBWYVinTzt4zjOyRMdrlndKswWvLhqLwBgUp8MNEtUv0j655TYx9yXv74GABAepse0Ac0D00kiColTedb/pXqdhDYp6qMbaiuDXof4SCPiI/2b919htqCwzFwZDFqDwKJSawBYWFZ1u7C0orKdGcVlFdbv5WZ7JtG+rcyMonKz/UPEYAWPgLWAmVqQqNhmDxKtgWN0ZZbSNmTV/rMtW2nUw2TQ8fqI6maQN2TIkBp7QeeL8spPuwCgQZR1KKYtgfHxhqMAgGt7N63WypiSJKFhtAkNo03o1lR9CElJuRln8kvsQd+pvBJkFZTi7PnSyiGg1u/nCkthka2VF3OKyrE/U6wPRr3OHvzZg8CoMMRHGtEg0vbd+rMtWIyLCOOyEQEi+v/DcX06STD/p/XrW+5wPE9r3xVV/t4ArsWK3J5brFnAlDhkGys0hp96+ydtl2BQWJ2WbDyG6HADRnVuFOquENUKe88UAABaJEZxDrwbBr0OcRE6xEUE7hpIlmWUmS0oLjOj0CEAtAeBlUGl2yCxMlB02Va5j21ESoVFRkFJBQpKKgCUanfKCwadhEij3h70RTsGg7afK+93HNoaad9mvT/SZEC00TrMlaOtap86GeTVNUVlZvvPtvl2tkzar3usEdHEXjWrGihgzSRkNLQWatFitsjILSpD9vkynD1fiqzzpThbGQTmFJUjt6gMOUVlyC0qR05RGXIKy1FmtqDMbEFmQal9yQgRkgTERYRVBX+VnyImVAaH1qpgVT/bgkf+cXMlGiqXOwQuok+jVjxjdlxCQSNyM1tkFDr87ogu3F7dSsqr+uhLdtSdmvZRxsncYjz4xXYAwIFnRvHDFiIBx84VAYDbDD8FhyRJMBmsxVziIwN//LIKS2XGUCUI1AgSCyvvLyyzDlm1Dlut/LmsauhqhcU6TzK/pMJDT8QZDTrEhluL48SEG6xfJtvPVYVzYsINTu2q2rM4TvVikFcLlFZYLwJ1EmCqHL9tkWV8t+MUKiwyOqbFonUtHsah11VlBQHPj0OWrVW8FIGfLRgsLK/c7rCtqBw5hWUoKK2ALFvnN+YWleOQF32MMRkqs4HKADAh0oj4yu0JkdbhpUkxJiREGqGr4xex4vPrqoIrrTliotl3x8BOa1jn+VLlP7dQzMk7lF2Io+eKMNhp8WJHjkGeViDqaTH3ms621idgHZ2g1/GfPZEnRyuDvKYJQYg0KGSMBh2MBh3iENgRWNYPNytQVGoNCK2BoNk+f9EeEFZuK6yc72gLFm3zHgsr9y8qNdunPZRVWCpHYZV56IV71kCxqoKq44ft1uurqlFY8fZRWGGICOP0HF8wyKsFSis/mTEZ9PbAwWyRsfJv6+Kol1+QFrK+hYIkSfbhBE28KDZWbrY4ZAMrg7/KLKHtdm6RtQpYblE5zhWVIa+4HLIMFJRWoKC0AkfV68u40OskJEYbkRwTjqQYE5IrK5TZvifFhCM5xrqgbG2deK0T7Ha5w3BEvcYfacf4RnO4pllsuKZzkCca4wUyFhz6/K8AgC9v74/uboY1OxaHKdeYZOhtvwLx/7Ck3Iy1B7LRr0UiIoz+BWWOWUqtDKw3/j6ZB4NOh7aptfdDLiItR84yyCNxep1UGUQFLngsq7DYg0Tb0NL84nIUlJZX3S6p+rnA/nPVNtv/Y2ugWIrs894NTTUadPaA0Fawz3pNFW6vAGu7voqLCGNAWIlBXi1guwg0GnT2C6WiMjM2VFa0HNY+OWR9q03C9Dr7HwNRZot1wrUtO3iu0ClILKwaSnquqKpUtNkiV64jpP2HTJKAlJhwNGkQgcYNItA4PgJNGkSicYMIpDeIQNOEyBpbIVF0fp1jpUytv7uOmTytrJVj4RWtUCGnUPlpozmE83S3Hs11G+SVm8Xm5HmbiRR9fbQ8umwnlm4+jjFdGuH167q7bbf3TAF2ny7Q/MDJcXhmRQDWIiwuM2P0q38AAPY8dSmHAVGdxEwehZo16+hfgRyzRcZ5RTBoLWTjMhrLaRRWblEZys0yyiosQtdUgLVeQ1KMCSmxJjRpEGm/vrL/HF9/1nhmkFcL2IZrmgw6e3XCCosMWGQ0jo9ASz+WTSBtep2EhCijfX1CEeVmC86eL0NmgbXQTGZBaeV3623btsyCUpRVWHA6vwSn80uw6UiOy7HC9BKaJ0ahVXI0WiVFo2VyNDqmxaJFYnTIh4OKflDmWOFSK1ARTe6UCw7XPJmrXJA3lEsoaJ27vEJsSQhv5xTq9f6/P5ZuPg4A+Hb7Kbx+nft2Y99Yg6IyM4x6HS7tlKraxvHtqpWxFFVUVpWpzSsuR3JM9f7T/vGfM3jlp314/qoLmEmkoJBl2T4nL51BHtViep2EuMgwrwsEyrJ1bn1OYZn9A3fH6yvlNVYp8oqtNRtO5BbjRG4x/jqaq3rc5BgTWqdEo3VyTNX35Gj7MmV1BYO8WsCWyTOF6VwKMwxqk8S0dA0TptfZF4TXIssyss+XWf8Y5RTjeE4RTuQW43iO9fbRc0UoLjdj7xnrOoWOYkwGdEmPQ9f0ePRsloC+zRv6PZzOW6JFQhTDKzWu7R0DIa14TJHJ02h3Kq9EcTuUBXcNGgG5Mgh2fwyvg7xq/LtgKw617Xiu2yDPsfdaGUtRjn/3PL22b/yyH5uP5OCt63sEbHj0zR9uAgA89MV2LLvjwoAc0+ZEbjGOni1Cv5YNA3pcql3OV5b8B4C0eO3/J0R1kSRJiDZZK3+KlBcsrTAj+3wZMisrutuurY7nFFd+FaGwzGz/oH3N/rOK/ZsmRKJb03h0S49H94wG6JgWV6uLhDHIqwXU5uTZ8CKg9pIkyT58tGt6vMv9FouMk3nF2J953v61L/M8/jmZj4LSCqzZf7byD9QBGA069GmegIs7pOCyLmnV8mmU6J+9Mh8WJNdqVS44XHN/pjIw9rTcgjdkWYZFhuYff8fhp1pDbssFM53eJr9C8Y9J6zl2fGwVAcjkOQa9ngLg51bsAQCs+Ps0LgvwHOYjZz2vAeqty1/7A2cLy7D41r7o0yJ4f+NlWUZphaXeDF2qbWzFiiIqF9ImIm0mgx6N461DMrup3C/L1ik4B7ML7ddVe88UYN+Z89YP1yrXfv5660kA1vWbB7dJwtB2yRjSJrlalyoLBP7VqAXUhmvadG8aX/0domqh00mVY8gjMaRt1bzLCrMFe8+cx9ZjudhyNAdr9mfjZF4Jft+Xjd/3ZeOpb3fh4g4puHVQC1ygEjwGjOhwzQqxzJtoAOi4pIhWRc4dJ/IAADHhBhSUVAR0Tt51727AmfwSrLhvkNvlNRyLi2hl8hRBnlaQJBCkKgLLEAR5WsGWY/IuEIVXHB+raJYzEBlEZ4EYeursbOXF/a97s4SCvKNnizD5/zbgxgubY0r/ZsLnmfvN3/ho/RGsvG8QWiV7P+T0u+2n8L9tJ/H81Rcg2sTLiUCzVTH0ZroAEbknSRLiI43o3tToMk8+r7gc24/nYsvRXGw9louNh88hp6gcy7aexLKtJ/He5J4Y3iElRD33Df8q1wL24ZoG5XDNpBgTGsdHhKpbFCIGvQ4d0mLRIS0W1/VpClmWsT/zPH7Zk4mvt57E3yfz8d2OU/huxykMa5eMuZd1RNOGgZ/PITpc0/YhBaAdqDjGYFrBW7EiyFNvk1dUbg/yemY0wC97soTn/IlYd9A6xGPXqXx0aRKv2sYxuBXN5GkFoiKBjGPwFIpMnlbw5hjEB2J+pONzpXVex/dSIIrROCsPQBEZt8euEDv2az/vw+GzRZj7zd9eBXkfrjsCAHj95/14+Rq1z7213fHJXwCAxqsi8OiYDl7vT9psmbzEaAZ5RMEWFxGGga2TMLC1dcmjCrMFfx3Nxc+7M/HH/iz0b1X7Rs4xyKsF4iLC0LdFAtqkxCiKXXRLj+d8PIIkSWidEoPWKTG4dVBL/HMyH+//cQjLtp7AT7szsfbAWTx2eQdM7NXU47Hmfr0Tv+/Pxtd3XIgYDyWYRd95tsVZAe3hlaIX/sUO68q522fVrjMwW2S0SYlGRsMoAFnC6/B54ngcrSyO49IIYRpFUMocCq9odVEkyHMMOGryPIJAvBaOT4dZI0MnK9oFPusWzKmeZYIBpGN22xf+ZlZP5RV7bkReO1dorSTITB5R9TPodejdPAG9mycAaBfq7vikZtZmJ4W+LRris1v74YkrOimyJ93clGSn+q1DWixeuPoCrLpvEPq2SEBxuRkPfbEDz6/Y4/Hi+oN1R3AwqxBfVFZVDARFJk/j/KKXmYrhmmrHkWUsXGtd6v6yLmn2D0YCdYHveBitpQDKBLMwikyeVnVNgcDIsVJnTQvyLIrhlQE4nkUsk2fxYlhnaYUZs7/cgeU7T4t3JIhRnugHH/5W2vU35g7EHEtyVTVcU3zZHyIiGwZ5tYxjZbgLmsSFsCdU07VIisYnN/fFvcNbAwBe/2U/3v/jkNC+gVqsGlBmtLQOKztc+GtddOYWlzscz7Xh+oPnsPNEPsLDdJjUN8NeZTJQD0m04IdjcKvVTrTwishlvOOcs5oV4jln3gKRyRMr5KIIyj2c9/PNx/Hpn0dx20ebhfuhtaajv0SDL39Xy/D39QhGhpQ4XJOI/MMgr5aJMRkwqnMqejdLQM9mCaHuDtVwOp2Ee4e3weyR1qEGz3y/C9uP51ZrHxyDHa1MomjWIrvAYTFUlV3e+/0gAODK7k2QEGW0ZzkCNVxTtJ+OmbyABHkCF/KOp6lp192BnpMnGrwplubwEJCdrcyceEPkoVSYLVi+8xQyC0o8N3Y8tmA7fzN5/hYlCuSHQlTl7HkO1yQi3zHIq2UkScKbk3pgyW39ArbeE9V9tw5qgTFdGsEiA7O/3BGQ5QREj1BaLhbEyIqf1duVmy2KC2XnVvszz+On3ZmQJGDagOYAEIThmg7H0bi2LhVcOqLMcR1BjWGMIkVDFBUnQ7AwoFYgKisC0EAvoSA2J8/TW8CXUEnkkXy28Rhu++gvXPbaH94dW/Bp8regTCjXkCT3bFVWG0ZzuCYReY9RAlE9IEkSHru8I6JNBvx9Mh8/7c6stnM7fsovuhi6u2v2nSfylIurOx3w/T+sWbzh7VPQIikaQFUV0GAM19RSqsjkuW/nWEFRKzDzPpNXs67cZVnsfQBYFy9/6PPtHqqxev++ClQ2193x3flx1xkAwJn8Ug8tnYn1198ho/7W72L9r+CwZZYbMpNHRD5gkEdUTyRGm3B93wwAwII12nPzRK6FRa/rHA+lmclzuEutVUm5GY/97x+3x84qKMUXf50AYM1c2lTNyQt84RWtDIrjUg9awZticXc/+6gMlLWPFcjF4UV4MyfvuRV7sHjTMWw+muO2jVnwsZoVQZ5AR70kckxfn+rqitMZo9VM5+yZPAZ5ROQ9BnlE9cikPtZlFNYdPIvMfPfzgwJaTMJT9GZv5n7OlizLeGDpNmw7lou4iDCM6pzq0m7R+iMoq7Cga3o8emZUVZ61TVcK2Jw8wSv2PMcCMRr7VCiGHfreL0C8guWuU/m44ImVeGv1Ac3jeZuh0Qp6RTNqjkGv1jpxjtle0YI+ngL9YMVUWlVYtQh/MOFnx/3O5Pm3O6mQZdke5HFOHhH5gkEeUT2SnhCJbk3jIcvWteRsLBYZp/Oqgj6h+V+C5xQdQqho53TV/vZvB/Ht9lMw6CS8dX0PNE+MsvahslmF2YLFG48CAG4e2FyxfqTt50DNURM9Tn5JVZCnlbnSCm4diayJKTrv7cVVe1FQUoH//LDb4zEDRXS9usLSCvvPWstAiAaNFi8yecHKnPm6xEA1xXh+z+njeq2BV1BaYV8nsSGXUCAiHzDII6pnBrZOAgBsOHjOvu3d3w+i77yf7LcDec0mOnfKMXvo2OxQdiFeXLkXAPD4FR3Rr2VDxXqRAPD7vmycyS9FQpQRl3RIVdxnCxT8GZ2YW1SG5TtP48jZQuFqjXkelnqwH0MwMBN5SUQrWJYKruEXyEt30fdBQUmFUDtF1lIwyPNmyG4g5++Va1XU0eqDaDs/+8pMXs1zrnI+XqRRjwijPsS9IaLayBDqDhBR9erbPAGvAvjz0DmsO3AW32w7gU//PKZo4xxEOXOcb+Zs/cGz+PTPo2iaEInbh7QSnpPnLuP39Hf/oMxswcDWibiut3W4qeTUbskma//Hdm3sUnXW3+GaW47m4MaFG5FbVA6DTsKTYzvZ79M6ZEGJWJAnmunUOTwsi0VWLZsveizHPc+XViDaFJh/BVpvG9G+iQbHjtlA7WBQ7LyAMmi3yP6vP2dT0xcLZyKu5jlbyOUTiMg/DPKI6pku6fEAgNP5Jbj23fWqbRwv+krKzThytggtkqIQptfh0z+P4pFlO1WH3P34zxlM/2iz/b7jOcXo3DjOfr/WNbasMqxu16l8/LgrEzoJeOzyjvZhYbbvsgzkFZXbqxde3auJymOpHK7pQyov+3wpbqoM8Ix6HcrMFvzrqx32+7WCBscqoFpTskTn0TkOqbPIMnQq+ROz4Pw+x3PuOpWPXn6suSkaPItm1PKLxYa5igZvjv3zZj23crMFel1gMijlXszJcxyqXHuGa/rZAXJhr6zJ5ROIyEccrklUz0SbDGgcH6HZxpbJ++dkPgY/9wtGvPwbRrz8G7LPl+KlVXtVL77PFZZh1ufbYLbI6NXMWvjkqy0ncCi70N7G3cX42fOl2HWqwH7bdmH+7m/WJRFGdm6ElpVLIgBVF5UWWcYvezJRbpbROjka7VJjXY7tz3DNp7/bhZyicrRLjcH6fw1Dwyij5rprZ/JLcDK3GIBzpkkjCHH4WbSoi7tmykDZ/bESHS4cw/Tu/w2IzLUSDZ4dW2ntIz7MVXROnsPPHvoa6LX8bLz5gEFRDVR0CQUfuqp4zvwO0hjlBZp9jTxm8ojIRwzyiOqh1inRmvdLknW44bQPNtrX9jqYVYjnlu+BKUz9z8brP++3B0Sf3NIXA1snAgBW/eNQ4MXN1ehdn27BLR9ust+WYc3QfbvjFADg1oEtFO1tmQcZwIq/TwMARnRUzsWz0TkEhGr2nC7AxS+udtl+IOs8lm21Lsnw7JVdrPP9OqYo2jge89Wf9qHvvJ9w4bM/Y+GaQ4oLe63heuKFV9TPq+xP1c+igVRpufuhtyJZOsegROtSX7QAii/DNbUSZWYvMnmiQ0C95U3A6E2hGPV9vP+gwKdF4GvYOox1zTkGeUTkJwZ5RPVQs4ZRmvdLkoRXf9qHU3klaJoQiQVTewEAlm09gUaxrlnAvKJyfFZZ3fJfo9ojTK/D4DbWAi+nHZZqcHdduPbAWcVti0XGdztOoazCgnapMejSJE5xvy1wKykzY/XeLABaQV7V0E417/9xEPsyz7tsf2f1QcgyMLx9Mi6oHOLaI0M5rNF2obvqnzN4cdVeyLL1PM/8sBvZ56sWvtYqDKKoOik4j859kCdWjOSMw2tSphEhiVzGC66QoQgKRANQ7SGnjj9rZPIsYucFnNfUC1wQ482hHGu0+NIH0eC0wsdiMDaOzyWHawaebbhmAtfIIyIf1dkg74033kCzZs0QHh6OPn364M8//wx1l4hqjEZx4Zr35xeX4+MN1qDtscs7YEjbJDSOj0BphQWHzxa6tF+y6RiKysxolxpjz+B1TItzaed4AZpTWIYSN1kkGcCXfx0HAIzv3thl2KDt5poD2SgqMyMtLhydGrsO1bS2VZ+Tt2jdYTy3YrfqBfj50gp8vc2axZs+uKV9u+P8QsB6QV5utuCxb/4GANwysDk6NY5FWYUF31dmIQHtddKUc9XcNlNcSbsLVkSzQEUOhXM0K54KBAyiWSqL4Hp1+cIFa8QCMtGlG6x9DH0mT5EZFYyeRIfCKs5j8f48jhyzoozxAs9WeIWZPCLyVZ0M8hYvXoyZM2di7ty5+Ouvv3DBBRdgxIgRyMzMDHXXiGqEVA9B3nfbT9mDtqFtkyFJEvq0sGaxMgtKXdp/tcUaEE3qm2G/YGyjMiTUdjG+ZOMxdH9qFXo//SN2HM9zaXc6vwSbjuRAkoArujZ2ud92DttQ0sFtk9xeqOrdDNd89Ou/8cYvB7D3TIHLPst3nkZJuQUtEqMUC6s3S4xUtLPIMr7dfhIncouRFGPCzIvbYli7lMr7qtqpXXiXVVjwzm8HsONE1ePXmjOmzOSptxENaBSBpZ/RjOjuirmHGkFPTpHYIvKiwbFFMIPofH8gM3nePMVmX4InH+YS+hukMZMXXFXDNVl4hYh8E/Dqmt26dRP+VPCvv/4K9OkBAC+++CJuueUW3HjjjQCAt956C9999x3+7//+Dw8//HBQzklUmzSK0y68YsumXNAk3v773C41xm37f07lw6CTMKZzI/u2hCgjIsL0KC53zBpZLww//+s4ZBnIL6nAk9/943K8n3dZP5Dp0iQeKbGuAanzn5j+LRPd9s221IC7a99tKkHmssqgdVw3ZRbRZFBWW7TIwII1hwEAU/plIMKoxwXprhlMtblgP+06g2e+Vy5GrnaBnltUhhdX7VUEo+6CH9EKlr6uHefpWFqHEg3KciovbgHtIafKOXlij9XTnDzlEFDNpl4RLaACOL22gsGT4/GFM3l+LuugDBIZ5QVaNodrEpGfAh7kjR07NtCH9EpZWRk2b96M2bNn27fpdDoMHz4c69atU92ntLQUpaVV2Yn8/Pyg95MolBJ9uHBonew+yAOAZolRaOAwtEiSJKTEmnD4bJF9m23UouOF7J+HqhZlt7HNkRvWLln1XM7r+PVr2dBtv3xZQmHtgWwA6lnEmHCDfcHuA1nnsf14Hgw6CddWruHXItE1g2kbrinLMs7klyIx2oiTeSUu7dT6+O7vB/HhuiOKbWqB2Y//nMHDX4ot7+A4dNLd8+IYbGlRzsnTGjYpFoiccziv6FBS0XUIPWfyqp6YQFbX9Gb6m1Zg647ovE5/z6PYn5m8oDpXOVwzkZk8IvJRwIO8uXPnBvqQXsnOzobZbEZKirIKXkpKCnbv3q26z7x58/D4449XR/eIaoTYiDDN+23Xf44Xb2kel11w3ZYcE64I8mxBgOjl5UVugjzHU7VLjVEsCeBML9mWUPAimyIDbVNi0LRhpMt9KbHhKCixBqHfbD0JABjQOtG+nlWj+HBIkvLCu8Iio8JswZ2fbMHyv0+jQ6NY3DSgucuxbft8teU4lmw8jlFdGmHHCdcPnWzX13+fzMPfJ/MxqnMj3OxQnRRwH9DIsuxU0VP9OXjwi+3qdzj5dY/YMHjRQin5lQE0oD1cUzm00v15RYNLQFmEJpBBnjdDPy3KqNlrosNvFZlQHx5rRZACYrK+X2wfdjCTR0S+4mLoAGbPno2ZM2fab+fn5yM9PT2EPSIKrjiPQZ71os0xyEuJ1f5EWW3IVkOnCxTna8H4yDDkOszBcuZuiKhjv3o3117Iu2oJBcd+eL4ovai9eoAZEVY1ZHNP5RDKy7qk2beZDHokRZsUcxfNFhmr/jmD5ZXLPfxzKh+f/nnU5dhmWcYf+7Ixc8k2yDKw7uBZl2Iv1sciY83+bEz+vz9htsj2IjXKNq59//Kv45jz9d84X1oVSLl7LhyXvtByz2dbhdqJZt5El5QQHXLqzRIK245VDd0NZNzivFafTu0TEfv9jvuJdcKb4jJq7XyZlykaZJP38ksqUF45nJaFV4jIV0EN8sxmM1566SUsWbIER48eRVmZcvjPuXOuw7T8lZiYCL1ejzNnlBcoZ86cQWqqeol1k8kEk4lDIqj+CNPrEGnUK6osOlK75ouLCIPRoENZhfjYs5hw5Z8Y20Wr7UI+o2EUcoty3e7vPCxTbXu8h4BV55DJ23wkB4CMbukNNPcBrPMR1TgGeTbtGykreyZEGRVBXrlZRtZ5ZcEaa1+UzBYZ/1m+S3HR7FiYxaaswoKHv9xuv9Bef9D1b6ntwj2zoAS/7snCwNaJmLlkm+s5na7QS8rN+GHnKZd2IrQu9nOKHObaaS6NIJZ5c7xHc06eY9DkIZj551RV1jSgmTwoA1Kdxhw2syJ4FTu+4z7CwzX9DNIc135kjBdYtixelFGPcJW/N0REIoJaXfPxxx/Hiy++iIkTJyIvLw8zZ87E+PHjodPp8NhjjwXlnEajET169MBPP/1k32axWPDTTz+hX79+QTknUW2klc2rulCsuhiVJAmJGp8qq8Vj0SblOWyHtR09I8F1OKSnY3rbyHb3z7szceX8tbhy/jos3XzM42H1brItYQbX7c5diI9UPm6zxWJ/7G1T3M9tzMwvwc4T+ZAk4Olxndy2+2HnKRw7V4yEKCMGtFIvOmORZWQVlGLMq3/gwc+3Y/yba920U97+eusJ3LfYNRh0lltUhvd+P+hyTjW/78vCI8t2emzn3B/RjJ9owRetTJ5zRjOwSyiIH1e5jINYJxTBq2hg6MN53O3PhdED62zlB0IcqklE/ghqJu/jjz/Gu+++i9GjR+Oxxx7Dtddei5YtW6JLly5Yv3497r777qCcd+bMmZgyZQp69uyJ3r174+WXX0ZhYaG92iYRWYO8UyrFPwD3F20x4WGAm33URLvN5FlvexwCKpDJ80QtWHv9l/0e93N3Br3O9bMxlyAvQnlxVmGpmgfXLDESB7PP24djObLNCQvT6XChRsXQ3/dZC8OM79YYDaKM+GN/tksbswy89ONee0bR3WttsVjn/3z651G0So7GF5tPuD3vV1uO481fDmBI2ySczi/F/7adVJ5TJYIpLjPjPz84VRH1c2kEWZbx1He7VPdxbevYP/cpRMe5gJ766C3RIaiAbxk2x+BVtNeKhd8F93F7TsZ4AXWWyycQUQAENcg7ffo0OnfuDACIjo5GXp512NGYMWPw6KOPBu28EydORFZWFubMmYPTp0+ja9euWL58uUsxFqL6TGsYkFrhFcA1aPMkxqRs73wxmODjRYxjvzyFe44BYbem8dh5Ig/HzhX7dF6gat09RX+ceuGcJTVbqgbsGfQ6NE2IxIEs10Xl7Rf4EtCkQQTC9JJ6MFg5ZLZBlBEtk1yreQJAUWkFPt9snavXu1kC/jysPjy+zGzBjQv+tC8lkRyj/prsOV2AWUu3o8Ii26ufOpMrH8MrP+7FjhN5eGBEW9y/ZBt2n1auRagWP8myjP8s342DDs+Lu2GYh7ILcTyn6jUUXidPo1250xjSoM3J8xTk+bC8hS9LYlj8DNL8zQSSe2fP24I8ZvKIyHdBHa7ZpEkTnDplndvRsmVLrFy5EgCwcePGoM+Bu/POO3HkyBGUlpZiw4YN6NOnT1DPR1TbhIe5//W3XWg6xzPOc+wcqWXdokxuMnmVtxOitOfTuT2X4rzabR3nQ905tBU6OM2fc3sON8cVyeRFmpQBdLlZrgqc4X4xelt2RII1GHTXzvGiOj1BverpucIylFVYoNdJuGlAM9U2ALB6T5ZirUC1xe4B4M1f93ssXCLLMv7vj0N49ef9+GVPFu79bKtLgAe4zhs7nVeCN389gLdXiw3/LHMKyLTmoWWfdyyA4z6T55y5C+gSCl7Ms/Mm62ejmB/nwzw+X4ZbOlbXZIgXWLb3rFbVYCIiT4Ia5I0bN84+N+6uu+7Co48+itatW2Py5Mm46aabgnlqIvJAK5PnbqhaTLh3QZlzIGm/aK38Hh9pdBtMaQVvjtUJPS3EbBs62SAyDIPbJKFlsnrmS/T8BpXhn85bnJ9bs8VivxCWJMntxZvtebed291wrQqHdklujmW7iNdLElpprHF4Ot86jFOj4CMA4Icd1sqgdw5t5bZNUZkZr/28z37bbcavchmHw9mFKCk345FlO/Hcij0u7dTeh7/szsSlL//ucjw1WQWlmPZB1dISFRqpPOdTBXYJBfHjmhXVNcWO78vQS38LpygzeT4cgNyyB3kxzOQRke+COlzzP//5j/3niRMnomnTpli3bh1at26Nyy67LJinJiIPwg0+DNc0aWTyVLYZDc5BXuXxK28bdBIiwtxX+RQ5l6dM3uA2SXjiio7o3DgOBr0OzRpGeXUuZ3qV8ZrOfXCuwOk4J0+CVmCmvO1u0Xp7MAgJCVFGl3X5AIeLeMl95tDWNwC4ID0eW47mum1XZragY1osxnZLczun8fsdp1BultE4PgKN4yPcDhE1W2Tcv2QbvtxyAu0bxWLXKde1AAHXCpo3vL8Baw+cdWmnFgwWl5kxc8lWZTvBpRacz+3seE6Rxr2uHM8reyhO60tBE1+WQ7D4mckrN/u3P7nHTB4RBUK1rpPXr18/VrgkqiG0hmta5KogwpHa8gE2asGWUa88h/PFoCTB7VIOWrGbu4IsanQ6CZP7NbPfTnIz58z1/OrnENnq/Nw6Zk0kyX1lU8fgDbAuxaDGMZNn0OsQHxGGHKf1Bi0OQWW0yYAoox6FKs+zbQhjYrQJEWF6FJe7D7j7NG+IdI2KqLYL/3HdGiO3uMxtkPfLniz8tjcLANwGeIA1cD2QdR7bj+eiU1qcaoAHuGaSlmw8hqe/34W8YuVz4m646fnSCkx8e53Tud0HLte+u97tfWq8yeR5M7TTRnRtPHf7+JKJ8+WcJCa7wDonj0EeEfkj6EHevn378MsvvyAzMxMWp/kQc+bMCfbpicgNzeGabq7fTBqBoRrnTN6PuzKx+cg5xUVvhFG9H1qBnDeFV5y5y6K5nkR9c5MG2ss+AGrDNZVz8pznKtrY5jnZHp+74bFmS1UAB1gL4jgHec4X4UkxJhSedc1A2QJQnQQ0bhCB/W6GWALW5SNMBj3iIsJcAihHDaONLstIOLIFeJ4cyDqP51bsQXG5GRekx7ttZ5Fl7D1TgLX7szGkbTIe/GK7aju1wGTj4XO46q11Ltu1YhhvC/c4ntfTOna+FFHZeizX/rNoUs3iwxBPRxUO40pZeCWwsguZySMi/wU1yHv33XcxY8YMJCYmIjU1VXHRJkkSgzyiENIO8pRzw2ycM3OOVDN5Btf2V85fZx+GKEFCZJj3f4Z0ir8l3u2bKJjJc+fOi1rhrdUHFNuc++D83Fpk2V4ARpIkt8NenYO3KDcBsNmifH2ijK7Hc24T6yZ7aHbIHjaMMkJrcQlbhrFhlFEzyJMApMS6HyJq0zYlBnvOuBZmsfnftpP2zOI2h0DG2aHsQrywci+Ky81osf6I23aOwUhBSTlGvvK7okqnu7Y2247l4uUf97o9vpq/T+YpsqOeAqIDWVVBtkiy7GCWMij3rViLL4VX/KvOSe5lF9iCPM7JIyLfBTXIe+qpp/D000/joYceCuZpiMgHQksoOG3XyuSpDW80qQR5AJBdWSIckkYmz+2ZnOfkeRflxQouA+HuqNEmA0Z0TMGKv8+4bev8uM2yWCavKjCzHjFCJXhzbGcLdtWOZ3YacqsWCAJVQYFO57qIuzN7QZhoIw5muy4BYaPTSWgQ6fkCtUVSlGaQpzV01NHHG47afz6osjSFjSwDmw6fw7oDZ9G0YaTbAA8AHAeerD2Qjfd/P4SfdmcK9cfm75N5GP3qHy590OK4EL1I8OU81Fm48IqfQZpjdU1m8gKntMJsX7ORmTwi8kdQg7ycnBxcddVVwTwFEfnIqLbgWyV3F21amTz19u4DScAa8ES6Ha6psZ+3YzQdaBWPUZ5DY7ioU1jn3DZM71pwxv6MSq5LLNhUOGfy3LRzHvKn9hyazcpMnrs1DquW5pMQ7ybbZ2PrV7yHAE4C0EBgeQxP2b5AT/vKLSrHpPc2oLTCoppldlRmtmDl36eRnhCJuz/dqliKQcupvGJsP56HWUu3uSywDmjPZcsvUZ9X6amfjkSzcp9tPFa1jw8DNn1ZtoE8s30AFqaX3M7dJSISEdQg76qrrsLKlStx2223BfM0ROQDtfXedJL1wrpquKYyeDF5WXglzOA5GnOX7dNaGsHTsgla3GXRvOH81Dn3xnmZBYtiTp774ZoWpxRqpJvsW4VZ+fqoHc95rUNP55Qk8Uyeu8DcsaFIJi8pxqRaGdQm0MU9TuUVo7RyIXnbgvLuPPvDbqw7eBZGg85jW8AaXBWUVmDo87+ipNx9e+dHZLHI2H4iDz/sPOW6TqDn07r0TfQp+9+2k16dxxmHawbHyVxrdjk1LlyxVAwRkbeCGuS1atUKjz76KNavX4/OnTsjLEx5AXH33XcH8/REpMGgksnTSRIssuy+8IqH7IfLOVQCSUeSJMHgZXbQuqPjMbzbVatCqJtTqNznnMlT3u/83O44kYchbZPsbd0tX2EP3ipvuwumLE4BnNrQ2wqnoZ9us4IO7dwFlTa2x+0pyNNJYsF0mF5SrQzqzKjXuWSs1DSOj8CJXPdDML0JGtcdtFbyFAnwnvjfP/j0z6MY1bmRZoAHVGXaZFnGjI/+wvK/T7ttq5XJk2UZfx46h2veca70KfYYWyZF4UDl0FafMnmOQR6XQw8YW5DXOD4ixD0hotouqEHeO++8g+joaKxevRqrV69W3CdJEoM8ohBSW9RbZ0vluaEV5KkFRXoPn0RLsF7oCx9Q5S5vs3oB+XRccr6p3KAW3L728377ru4ynM4X9c7DPm0qnIqqqD2HFpciLm4yefb5fUCE0VNQbv0e4aFYjgRJKJiWICHK5FoZ1FmDqDCcyfc8XDItPlwzyAtW1f//W3MIAPDFX8c9tj2dV4JPNhxFr+YJmgEeoMyQmS0yft2TiS1Hc/H2bwcU69S520dLekJkVZDny5w8RXVN7/cndbb3bxqDPCLyU1CDvEOHDgXz8ETkB7UAzHmTS4ZKKzOnklLTe0izSZL7Y2rt6U91TQAe14PzdFyd052uz5P28g/ugjfnwivuAmDnAE7teBVVk+3ctgGUwzq1ivE4ns9TJs/6GD2/MJIkNkcyIcokFOR5mitYEwqE/Hf5Huv6gb8e8NjWIss4k1+CgpIK/LEvC4/97x+BfZxuW2T8cyofR84W4UDWeZzMLcYvezIVz6cvT8vrv1TVYeVi6IHDTB4RBUq1LoZORDWHWiDiHJQ5Z6i8HVkpEoCpDRv15ri+5OWMBp3nIE9zTqA2rQymBMltARtbcsb2+Ny3UwaDWkVEbD1x18aWkNFJktthpFUHs1X99BwMilQ91VpOwlFDN4vCO2vgYU5hTVjAO7OgRLhtUZkZl778m8dMpyPb0ElZljFlwUahNQl9GW7pWMW0BjytdYZtDUYGeUTkr6AGeTNnzlTdLkkSwsPD0apVK1xxxRVISEgIZjeISIVeJYBwzlC57KORyfNtuKb7gCdY1TUBsSyTFk/n15pnqJ3Jq1wM3cNxnCtnaq9faMsKqreRHYq9hAsEb4DneY2295GnuXRay0k4Eq2I6qnYS03I5HkTEO05U+B1YPrnoXOY/P6fGNkpVXjRecenpcJswRd/HUf2+TL8vDsTJeVmZJ8vRYNII3o1S8DIzqno3zJRub9XPSQtu0/nAwBap0SHuCdEVNsFNcjbsmUL/vrrL5jNZrRt2xYAsHfvXuj1erRr1w5vvvkm7r//fvzxxx/o0KFDMLtCRE7CVAIw5+DF+bZWJk8t8PEY5EnuM3mi1TV9CfjcBTzOfXPH/+Ga6vfvPJFf2UZ7uKZz5Uytx2MPBN1l8hzW0wv3UFjH07GqGlq/mcI8BHmS+4IwjtwtOeHM3YLvNjUh41QhUEDGxpfM4xebjyOzoBQfrDsivI9FlrF2fzaW/30aKbHheG7FHpc2Z/JLsft0ARapLTbvoUAM4P16lvXRucIy+zDatqmxIe4NEdV2QQ3ybFm6BQsWIDbW+gcrLy8PN998MwYMGIBbbrkF1113He677z6sWLEimF0hIidqAZjzhZhzC0+ZPtfjebgfGnPyBDN5viynIBTkeXGf8/OmPQRVQpinYKryu6fCK7YnQmi4pps+5VYOBdRJInPytINP53OaDHoUwHWtOMd2ImsvelyyoVKUh3Y1YbhmeZD7UFgmtoC8IxnWwkC2iqLessjWYG7niXw0S4zEfYu3YceJXMW8v7iIMFx2QSNcfkFj9Mho4PEDoPpo1ynrhzwZDSOFs9dERO4E9a/Ic889h1WrVtkDPACIi4vDY489hksuuQT33HMP5syZg0suuSSY3SAiFWqBiKegTHuumUp7gaAw2EMn1fgyD1B5Us2bmgVqJMlzYFNVNVO9nUjhlapjeQ4Ebef0FPza+iWyNIa1nadMruQ5Kwj36wU68zRXsCbwJpPnC18CWYtsHRrqKxkylm46jge/2I6kGBOyClyL5OQVl+Oj9Ufx0fqjAIClt/VDr2acquHon5PWIK9DI2bxiMh/PixQJS4vLw+ZmZku27OyspCfb/1jFh8fj7KysmB2g4hUqM2vc81QOe3jZUTl8dN6yf28M81cmJ9Dv0SyR1od8DRc01P3PJ9fO2MmsoSC8kieAzidJHkMfqvmCmq3s73sIsN1RYI80bUNPWUia4IKN0sfBIpP2UpZRlp8uM/ntFiA9YesWUC1AE/NVW+tw7u/HfTcsB6xZfLaM8gjogAIapB3xRVX4KabbsJXX32F48eP4/jx4/jqq68wbdo0jB07FgDw559/ok2bNsHsBhGpUMuyuAzXdLqttcacWuDlKRiTIKnODfS0r78DvUQyed5U13RuqxXbSBBfq89jZg2eM2ai8+gkSSAor+Q5k2f97qFZ5TqJnv8NifZLax1HR0JBfqVAjyqsCPJwTV+CPH97JMO3YdNPf78Lk95bj6Iy90N665N/GOQRUQAFNch7++23MWzYMFxzzTXIyMhARkYGrrnmGgwbNgxvvfUWAKBdu3Z47733gtkNIlKhdvHq6TLN6+GaApkcrYqdWvtV/RycOXmi51e/7TnoEjm+6PBJ7aBRfOik5zbW754zfrZzeoryxIZrigZaou8FU5j46x/ouVEVlpo4XFOGP92SZdnntfLW7D+LDnNWYMvRHN87UAeUVViwP/M8AKB9o5gQ94aI6oKgzsmLjo7Gu+++i5deegkHD1qHZbRo0QLR0VWlgbt27RrMLhCRG+qZN6fbTvd7W3hF5OLcXRutXRWLoXvVIytPwQzgXXVNl3017xPJItrOI9rOc1DpOasqEpQLFl6xBZ8C/TcJVToVe5VF3wsRYXoUlIhlj6JNBuQLtI0y6oWKnpQHe7imD8GWLPu3vIQs+3ZeR+PeXItrezfFY5d3gMnTeo110L7MAlRYZMSEG7hGHhEFRFAzeTbR0dHo0qULunTpogjwiCh03AyS1LypmclTuUskGHLbxMOQR63zeiK0ULfm/tptRYIuzXMLBma282ier/K7p/mUOkkSznCKFl4RyeT6u5yFI9EPIUSrdQKASXCen8h6f9XBp+Gasn+VRy2yHJDlKT798yjaPrIc//fHIf8PVsvsOmUtfNO+USyXmyCigAj4f6Xx48dj4cKFiI2Nxfjx4zXbfvnll4E+PREJUrsg9lhdUzOYcL1PJDPk7oJGNMjy5XLI/0so7ShPc/kHgaPbsj0eE44CGTPbxbfAyEmhoAwQGa5p65fnYZ0iwzVF53uJXht7M1xXdP6eN0NAg8mX6p0WWfYrEycjsAvNP/HtP3ji239w/8VtMPXCZogJ117/sC74q3K4KitrElGgBDzIi4uLs1+0xcXFBfrwRBQgQhkl54IiXl7HerzIl5RtrurRBEs3Hxc4sn9hmkjGx5tP010Lr2hl8jwft6zC4vE41vNWnk8jOCstNwsdSycyJ89e9VNsrqBI0Oiuuqqybx6bVLYLbDAIiFX/BLwr5hJMvmbU/InRrHP6Aj8M9YVVe/HCqr2INOoxa0RbjO7cCMmxvlcBram2HcvFsi0nAAAXtUsOcW+IqK4IeJC3YMEC+89vvvkmLBYLoqKiAACHDx/GsmXL0L59e4wYMSLQpwYAPP300/juu++wdetWGI1G5ObmBuU8RLWdN8MGbTQv2lWHa3p3jvjIMIftYkMefRnaJBKsah9WeUHr7RIKnpSZBYM8geGaxYJBHiCeyRPJ0AkdD96/RzweUKiZ+AskGuR5O181WM6Xel+pcvdp39fIA4CDWYU4mFXo1zG0FJWZ8fj//sHj//sHAJAUY8KIjino3zIR3Zs2QHKMSbhibajJsoy84nLsPl2Abcdy8du+LKzZb11+onezBAxolRjiHhJRXRHUSQRXXHEFxo8fj9tuuw25ubno27cvwsLCkJ2djRdffBEzZswI+DnLyspw1VVXoV+/fnj//fcDfnyiukK18IqHfbSGa6pV1xMp9uHYwvFCTXTIoy/X1oG+IHc+mr/VNW2ZPElgCQJAO1Cylez3FHBZM3melmzwfD7H+z2+/pLYayH6egUlkyeYoaspQV59kFVQqlhY3VFKrAldmsSjbUoMmjSIQGK0CbERYYiLCENshAFRJgMMOsmeubb9Xnj7YZEsy6iwyCguN6Ok3IziMjMKS804X1qB/OJy5BSVIb+kArlFZcjML8XJvGKcyivBqdxi1QI9Y7o0wtPjOteaYJWIar6gBnl//fUXXnrpJQDA559/jpSUFGzZsgVffPEF5syZE5Qg7/HHHwcALFy4MODHJqpL1C4lPBUU0boO8qVqoCQpjyl+ke4QDHp9VjHeHNdlPUHNAFX8yCLDXQGxdeQ8PbWSJDZvz5t+eZ5SKAX0olY44efFBb1oJo8xXs1wJr8Uq/45g1X/nPF63yijHkaDDpIkKX6HbSNRLbKM8goLys3WOYz+FKtpHB+BDmmx6NeiIYa2S0bzxCifj0VEpCaoQV5RURFiYqzrvaxcuRLjx4+HTqdD3759ceTIkWCemog8UC284nSZ7M36b+U+FHxwvizXCwZvivt8Ga4pNCfP/X2e5i9pLqTuRXc9LkEgmDEDPFfXtA6dFBuG6fn5s2VHPDcTW2YjwJk8oVZW4kEeo7zarrDMLLQMhjNJsi7LEWk0ICbcgGiTAQ2ijIiLCEN8RBiSY0xIjQtHWnwEUuPC0SguHJHGmlGNlYjqrqD+lWnVqhWWLVuGcePGYcWKFbjvvvsAAJmZmYiNrTkVpEpLS1FaWmq/nZ+fH8LeEFUPtWtSj9kejft8C/KUAZHjBb/wnDxfzqmyU7vUGOG5Sc5BnvPhvJy66JZ40CVwLIFGwhk60eGaHs4nwXPwKXI+b9t5U0BIOMgTP2Sd8c4NPXDros2h7gYax0egQVQYEqNNSIg0Qq+TEGUywGjQIS4iDOFheoQ7VD+1yNbhlgadDhZZhk6SoNcBkUaDfX6fQWcdvuw4pFivs1aDNeisy42YDDpEGPUw6nUcZklENU5Qg7w5c+bguuuuw3333Ydhw4ahX79+AKxZvW7dugkf5+GHH8azzz6r2WbXrl1o166dT/2cN2+efZgnUX2hnslzvu2c2dPK5AVguKbDhZLaHD/H/fyh9thbJkc7BXniJ/Em4+lN30WGWAJigZJYpdNA9Us88yYHMJMXjMtskcXagfo5XLO6Kl0Oa5eMMRc0Qpcm8UiNDUekUc/MKRGRB0EN8iZMmIABAwbg1KlTuOCCC+zbhw0bhnHjxgkf5/7778fUqVM127Ro0cLXbmL27NmYOXOm/XZ+fj7S09N9Ph5RbSByieR8HaUVBNiKhfjTB8eL+Zyico1+OQzr9OFaz5fHrn087WGuyvsCNyfPm2OKVbH0FAiKDde03StSeMWXKq/u24k1LK8Q/0BCZL4jAFT48CGHKEnyb4mDYAlWmBVl1OPFiV0xuE0SwgUXoyciIqWgDwpPTU1FamqqYlvv3r29OkZSUhKSkpIC2S0Fk8kEk8kUtOMT1UiqwzXFhgeq8WW4pvP5RC+oJcXP3l9q+lJZ1JEM7fGamuvkeXEe8SUUPB/L03NrkcWH64oO6/REkkSLxogGu2LnLa0Qn3flaeF3m6Jy75cuEO6DTvIpU+7JPcNa45Wf9vm8f6ArijaOj8CS2/qhcXxEQI9LRFQf1bmZv0ePHsW5c+dw9OhRmM1mbN26FYB1fmB0dHRoO0dUgwgVH3G+rbFLhY+V5nzJyvmdyRMIcL05rDcZT28O7LHwSuV3kUDJ0+tdbrYEcE6epOifjfO8RwmSWCbPcxPFeT0p9SLrLPrBgzfZQW9ZP8gI/PEHt03yOci76cLmAR2i+snNfdCfa8QREQVMnQvy5syZgw8++MB+2zb375dffsGQIUNC1Cuimkft+ky0hL4aX8qJOx9OZG6Z836+XGeqXbf7c73qaS6j6H0ubQWDrkCsNVduloWDSo9xjxeZPLGhpoHN5Fm8GPvoae1AG18/5BA7tm9FjTzxJxM357IO+PtkXkD68ee/hlXb/D4iovrCixpjtcPChQshy7LLFwM8IiW1anAuAYhzdkvjotCbC2fHwzseMRiLWatRLTrj8tDdn8Sluqbz86TxlzWQ2Q/b6yVyTE+xilgmT1J8d98vW3v1/RX9CuCcPNHCMd6MLBbN5JmDFIgBVeu0uTOwdSJiw73/zNbfgpD+DtdsGGXE34+PYIBHRBQEdS6TR0Ri1C7PXC7Kne7Xuij0pTCE83A94eGaDj0L2HBND7e9Or4f+3rD9noEIpNnkWXhap6iwaDLdpV2QksoeGzhXUutyq3ODKKFV4KYyRPh6++fL/47oQsA/4K8bk3jsXR6PxgEq5cSEZF3+NeVqJ4SGSbnfN3oeFF404XNFff5kslzJlx4xTEw9KXwiso+/pRkdw2GA5vxdHteL4I8jy1kL6pmejhU1Tp5ypbO2UQJYs97wDN5XrwG4pm8EAd5Xrbv0CjW56zy1T2t1ad9zQRe0ysdX9zWnwEeEVEQ8S8sUT0lks1yTg84XtT1b9lQcV+1Dtd0e8N3rlkm921dgl+ntppBXkCDAfHhmp6CKRkCT6XwEgr2KE99u8PhhIZrCr7IooG6NwFZbcnkefv7982dF/qUiZs9smo9Wl8+GLm2d1P858ouXDyciCjIGOQR1VPqwzWVW12uWzWGVvo6JcnxnMIXfpLqj16c1PM2rwqkeLFOXiCJVrsEPAdTsix7sf6dYL889EmCJJwpEyF6JG/iIb1g4ZVAZvL+b2pPr9ob9TqvgzyDXudxnqaa6YNb2n/29qUL00uYN76z9yclIiKvMcgjqqdEPsV3vnDUmgvnayZPeXwxgV6fy3puP4ZruhQY8bMzouet/C40XNPjnDzx4M1TMO7uXufdrJk8oVSeENH3hVfVNQXXyQtkkNc6Ocar9g+PbFctc/L+Paq94ra3v4frZg/zqj0REfmOQR5RPaV2feZckML5ulGz8EoA+iFeeMVxn8BEVN4Eap4uqIMRhKqpWgxdpPCK9v0yPD+X9mqegv1yOZzKWoSBXCcvEEsohDkFddUVsDtyDqK1XruHR7ZD65QYr37/fpw5yONxnfXMaICbByrn4XrzPn9kdHskRpvET0hERH5hkEdUT6ldoDknI1wyeRoXdb6tk+db/kyxGLov+wtuEz6ejxlJf4kOnwQ8X5B7U3FSdNF0l+0q7QKRhfR0Xmfuhha3SYl2WRfPnwyvr1wznu77YJ8z6MWvX6vKTKE3H5Bc2inVdakQwd1HdkrFzQNbCJ+LiIj8xyCPiNxyvu53vPgMRJFI54tE0WP6kv3zti9aZLgfxgpoBy4BLK7pUF3T/2OJdEv0fPbg0+NcRSmgWU/RgMxdJs8iu97nrgDk57f1E+5Xcox3GSyX502jrbdzGtukRNt/Ft316p5NMKV/M5ftonNoB7dJEjsREREFDIM8onpKdbim06W+cyXIQGc1JNHxes77OR/D6/O67hTI4inVX3hFYLimhwtykUxeVebQUyZP/X7VTJ7AfyHhuZp+/kdTew7cved7NksQPu5nt/YVbvvI6PYuwZe7QDg5xoQJPZoAcP3ddefTW6r6Ihpg3za4JcJUol2RGO/Fqy+wL7lARETVh0EeUT0lcoHnUlzT4S9GIDNS3vI3iBLZ3avqml4MrwskWx8DMSdP6HyS8runc3ks5AL15zk8zGnIpPBcTf8epAzPy2N4q2dGA7RIivbcEEB6QgRuHtjC9f3jpg+r7huMmPAwAGK/j62To9HQYV6caJCXEhuuul1k//Hdm3C5BCKiEGCQR1RPiVzfuVbXDHAffL4sd5yTF4ILSKcL6pBdwkqKbwCATo1jVZt6npMncjqxoNLtnDyXOV2SattbneZvCQd5/r4QKlGevwG7N5+FhBv0led06oOb9o6VP0XO0zU9XnlcgYf28/2DEWUyqN5XXQWGiIjIewzyiOopteDI+ULf+bbjBW8gEnkuc/J82C9Q15ku8+y8OG51Ze5czlv53fFi+/0pvTTbuiO0rIDwHEBrA9GMnzNf31uiCSOt4zu/D6ozCWWbX+ccPLkLpoyGqn/hnobbDm+fjEcv66A8rsCD08pCetp94Y3q70UiIgo+BnlE9ZTaBZrnpQEc2/of5rlMyRM8pr/X3UJl+70J8nzvil/Ulipw1xdPgag3L6en7Kl9uKaHdpKk3l9PHzYA1gXAZ49s59ozlQNOHyxW2VEt6PH3tU2JFS+6Ypu75jonz7XtG9d1V8yT8/Ty3XdxG8RWDu3UOq43tDJ5HRrFYkjbZP9OQEREPmOQR1RP+VSwJMjhjHgmT1L92a9zexHkBHreloip/ZtheHvlRbNaJs/dSySyTp4oycN/Drfr5Dm3g6T6+jln09Rem2fGd0bThEjFNp1K0Pj45R1xQ98M7Y5Uun1IS9eKsipP3KNjOrhsU9MiKQqPXd5RqO3cyzrYK1iKVGsd3aWR0HFt1Kpwehpu+db1PTTv19rf26qfREQUWAzyiOot7y/CHK/pgjJcU/Cg1XH9qBXQOmcxq2O4pmrlR5UCJ+4uvAOSyZO1z2E/l+h2SXxuqOt7xXWb2mNMjQt36a/ac/nW9d0xvnsTj31OjDZi2oDmnjsNYMHUXkiOUS9a4qxvi4b2wMgliBZ4jjy9fikq/dB67v97ZRdc2ilV85hawT6LrRARhRaDPKJ6SuTi2jWYCXgvFMGU6BBQx3186ZLqEEHnNkG6RvUvOPac4XHXbY+ZPK8WQ9e+333hFafbAscC1AMY2X4EZb/UAj3nrJLaI23SIFL1PufjeZPx9aYwieNxnffz96345e390SDK6LLd38Ipavvb1sR76NK2fh2biIj8o14yi4jqPN+Co8BGPjWp8EooTRvQHO//cUizjQz3VReVz4ePmTwPfbS2sYVWnubk2YZrepqTJ0HtnegyJw8yJOc+qnRYggSdBJidjiX2gYbtu/YHG96837zJZkWZ9FX7edjto2l9xDsBoHvTBqrbtYI8kXX31Po597IOSIox2Zd2ICKi0GAmj6ie8uVTfGXhlcD0w7EbFh+O6UvgKTK8UqtFIJcIzGgYiWt7N/XYTvX5Vilw4i5AcHzI8yd1x5R+ynlq3mTyHI/Vu1mCx7lbNs4LiEtQD5qce2IN1JyzcerDV13fD7LQe93WxPmoLlk1L35vIsP0nhsBuKFvBjIaRlWdA86Ptcqbk7pjQOtE4T5o0QomwwX6rva86nUSAzwiohqAQR5RPaV6ce10oe9p6JrfffBwfrf7+ZnJq2nJP9GEj3MztQDXXdDreEGeEhfukmUSyuRVNnIc/njXsFZo3yhG2QeVLrx1fXcMdaq26O61c83kqTx2ySW3B7XEYJheB73Am8T2mFyWDXFqd51AQA4Aj13WQXWIpJonx3ZSnlNjrmqgPlyxnsf1eXngkjYY0jYJozp7Luyi9rRy7TwiopqBwzWJ6ilfMmDKuCAASyhIvg0ArY4F0KvrWlUtS6XaDmqFRlQaujmUp0BSLYv6yOj2eOq7XQ59sB2r6mAW2fX1UHt9Lu3UCJuPnHNp59zygvR4oeqaMSaD6jqOjsdrGGVEv5YNUV7h+XjNHDJpymNW/fzRtD7o2yJBtZ2zqReKFWdRoxUoCa1n6CAm3P2/ebX3xJ0XtRY+tlo/WXCFiKhmYCaPqJ4SGSbnuk9wM3miF7Ai68J5y/XUGvOVAjleE+LFR1yDqcrvAplNZYEbz+ebPrgFbh6ovr6cY3/VKl/q3PxncX7/SE6FUmaNaIv3p/R02c85wL20Yyou6ZiqMrRS+fjX/2sYIo0Gt/2x2TbnEkQY1YcnOj5vA1onwqAP/r9NtUqi9p8Fj/Hqtd1wUbtkLJnez22bYBReEcmaEhFR8DGTR0RuBTqY8UR0Tp7fhVeq4Tr0ixn9cSi7EA8s3abZThacMwaoFQGpLHDisM39EgraxxYZKmtr4hiYybLskr0RGTJqbafs+80Dm8Nk0Lvs7/y+eOsG9TmA1sIrVfvaFgv39PzGRdasOWTO/XV8/J3SYoWO0TEtFpdfkObVebyl9uGEp4CaiIiqB/8cE9VTasOqvFoQXLDtha0aur1Pcsq8mAWjPGUQEJiIzXmIoL8JiR4ZDTChh+u6a2qEqj+qtFPbz92hlBf0np9ntUBNrdiJNcPouV/WPricxCkr61321GXhckn98dfEeWK2OYBqwymdnycZMn6bNRRfzOiHFknRQscXecT+Pi2qy1XUwOeaiKg+YpBHVE+pXYqJlE331qKb+mDLoxe76YNTxkIwyHO8CA7WNaXWYf15llyKiojOyVNdMsC1P6KZPLV++NIHi+x+XTfnnri2U2be3D4Nsvo702V4r6T+2lRXdimycshnQ4GCKx/f3AejOqfi89v6u9yntjZf04aR6JEhNh8QEAtsgxH8Oq9JSEREoVGngrzDhw9j2rRpaN68OSIiItCyZUvMnTsXZWVloe4aUY2jXl3T6bZGOCMa6Oh0ktsqg84l730arinYD8X+antV89BUR2LXxbLrnDwvhss5XtCrLy7uOZPpLtBy7r/7tfo83Fbdy3pekfmaOklSbedPMBMbIT6Uc+lt/TC8fTI+uaWvx7atk6Px5qQeaJsa47GtL29NkYccjHjMl2VQiIgo8OrUnLzdu3fDYrHg7bffRqtWrbBz507ccsstKCwsxPPPPx/q7hHVKKrl97246AvGfD3xyoGBvzoN9nIRWkSCEFltHYHKDY5Pm7vhcp4u6H19PdX6ZetCvNNcN9U5eYqsrHonZdldYOp0Xjd91Hp+PT0vw9sn46oeTdAlPV67IYCOaXF4b0ovj+0AL99fQVo/MtCZvJhwA+K8CIqJiCh46lQm79JLL8WCBQtwySWXoEWLFrj88svxwAMP4Msvvwx114hqnBozdcahH77MYfMlGFOdyyZZ13JT6ZYLbxYO90RtTptou6oFvKv64364pqT6s+PxhTrhpEVSlNvhmv8e1R59mifg1Wu7ue2T4yFtAZdWdUmt7ZIk4f+m9kJsuAGvXNPV5biOPripN5o1jMRnt7qvPglYhx8+d9UFuKFvhmY7b3nzrvVlGLXIr0Ug/wb898ou2Pjv4RyuSURUQ9SpTJ6avLw8JCSIz2Mgqi/Ui3a4zgWqzj6kJ0SicXwETuQWa+/n5mfh87rZ3jDa5MPR/Ce+Tp56MBVlrPpTLrJOmVoTl+Gaqn2o8u1dA3AmvwRtUmJwrlA5JN7Wz+TYcCzWKOEvSda5Zs77OZ9bBvDwyHb4zw+7cZPG+nM6CejfMhHb5l7iMagd3CYJv84a6vZYWvsGgjeHDdYQyEA+tgijHuFh6stQEBFR9avTQd7+/fvx2muveRyqWVpaitLSUvvt/Pz8YHeNKOSqY0Fxj31QucjUWrxZbb9AXac6L8ytddxAXnNLktjcKLWFw23PQ3pCJO4d3lp4qFyUyfU5dl1YXPsYnRrHoVPjuMq+Kbl7PGrVMGPDw7Bu9kUwGdwHCBYZmD6oBS7tmIoMh6DQme0Zqs6htr4K9u+ft4umExFR3VIrhms+/PDDkCRJ82v37t2KfU6cOIFLL70UV111FW655RbN48+bNw9xcXH2r/T09GA+HKIawd2QxWrn03wjh58DFeRJ4uX8AxnlWYM8zw8iOdakWajk3uFtcKNGlgsAZo9sh5sHNEebFNdiH7aHZJtHd1G7FNc2bh63c0Dh7rlzHXZobdcoLgIJHipSSpKEZolRTmv0ObWpFf/RKnk1/9X7N5wv2b8HL23r9T6tkq1LOvRt4X6pFCIiqn61IpN3//33Y+rUqZptWrRoYf/55MmTGDp0KPr374933nnH4/Fnz56NmTNn2m/n5+cz0KM6T+0as01KDE7lldhvBzsXICF484283V99oKC6QGZJJEiaj+fNSd3x3Y5TuG1wS/zrqx2K+7wtnDF9cEv3d1Y+pN8eHIqTucVol+q66La7YMPsMjfOzSm8zBZ6Oq9rcFl7BHu4prdT49Y+fBHS4iO8Ps+3dw1AYWlFyIY6ExGRuloR5CUlJSEpKUmo7YkTJzB06FD06NEDCxYsgE5ggSSTyQSTif+gqH5xzIgsu+NCGHQSPt5wtJr74Nu8P8dMkb/D3q7v2xQfrT+KGUNaIre4XNE3d/yJ8dSWKtAaXjiqcyOM6tzI2tbpvmBkXmPDwxCbWjXsc9G03rjh/T8193GcE+hNv9w1E90/0um8NXHRc3c89XRYu2T8tDtT+HhPju2ER5fttN8WLYDy+4NDkVdc7lOABwDhYZyLR0RUE9WmwS0enThxAkOGDEHTpk3x/PPPIysrC6dPn8bp06dD3TWiGsfxEjAtPtw+v6p6++B6ISoyn0rRxM/r+qfGdsbep0aiRVK08KECkckb3j4ZAHDThc2Dsl6Zt9xlVAe2TnJooy7KZMCcMR3st0XnxLk7nmgBoIs7pCA1NtzhvEKnDZgFU8WWS1Dj6Tl6d3JPr453fZ+m+G3WUAxvn4JBbZLQWDBoS0+IDMnvPhERBVetyOSJWrVqFfbv34/9+/ejSRNlKfZAljwnqgscfyPczqGqhl+bjIZRXu8T6MXQjQZd5XFFh2v6cFInb07qgf2Z59G+UQyKysxC+7hU1wxgUCPyWmu16dWsqopxoGMtdwGoXifh2QldMOX/rJnG6s7kDW2X7PO+Bg+RvUiVVEeSJKFpw0i8N8W74JCIiOqmOpXJmzp1KmRZVv0iIiWzQ6RSNbRL/HclPMz/Px+SBPRr2RBPj+uEpbdpr1em3E+7PL7P/XHqmzsdGrkWLvGW0aBDh7TYyuJRYvu4DNcMYDjl719J5dqFbs4ha992pzr/hCdGaxeA8dXA1omK296sJ2cy1Kl/1UREVA34n4OonnL88MP99abr1fWsEW1xacdUDGmbjPf9zBrYTjupT4Y9EyRy6RusfI1OMHi8fWiroJ1Xk3N1TT+eCOcPv0Q+DNNq4dgXd4+nqcbyB1q0MqeBrrT6wU29/T+IijcndcdDl7az39YLdPbNSd2RGG3Eh0HqExER1V0M8ojqqViHNdXU1k1z546hrfDWDT2g10kY1t61zL5XfLwo93u4ppudRIOE8DA9bh6gvVxBIPoTrP0A14BNJFmmFQgqi+Goi4sIw48zB3k8q6+PSytYvqFvhtAxOqbFYf6k7r51QENMeBimD2qB/i0bYnSXRkLDMUd1boSN/x6OPlyegIiIvFSn5uQRkbjwMD1+mzUUkgSE6a2f97hew9eAiiAqFAGFD10UGf5XnY9cdNilczt/hmv6OnTSHdFiOI3ivK/iKLrMhtaz8eTYTli0/ojYcYI0t0+nk/DJLX292qc2LOxOREQ1DzN5RPVY04aRSE+oGkLneqEf3MlQvgYpwouW+3NcD4f19dpbLZgSnZ7lshh6IAuv+Lm/Yqirxmvi2M5dYNnPKXOlFYCKDBP1Vk2odkpEROQPBnlEZOfLwuT+UF2UXOAC2995WG6Ha4YocykanARy8W+X11oglae1dIQy2HJ/DJGH2r9VIj71MuMlemyx44gfKDnGusZqIAoRERERBQr/KxGRXXUXovX5mtzPi3l3jzNYGUJPRGMK5377V3jF6bbvh3KhFSQ53qV1zn4tq7J5ohWSAzW00ZuA7aOb+2B4+xR8MaN/QM5NREQUCJyTR0R2zlUMa+rqIyJFPnw6boiG6YkGJ6M7N8JXW0447qnarlFcOE7llWgey6Xwip/r5DnSejS+vHbV/T68sGUiRndphHYpnpfKaJMSw7XpiIioxmEmj4jsqnu4ptowRZF4RzEcMEjDNWtivYth7ZOx4MZe9tvu+vjDPQODcn7Rd4fWc+fNvEebQCw+DwBD2yYBACb1aarZTqeT8MZ13XHXsNaBOTEREVE1Y5BHRHYzBres1vP5vnSAfxHY4DbWi32jXvkn0N+lGYJNkiR0aRxXddtNu/hII9o3itU8VlunLJVIgC+cydN4fZSBvdizrNU3b7J8r1/XHe9N7olHx3QQ34mIiKgWYpBHRHatU2Kw+8lL7beDPUzO12BNsZcPfby0Uyo+vKk3/nh4qPvj1lCiz5mnVpP6NMWDl7a13xYarqnxZIsWw/HlOdYqTOPNyx9lMmB4hxSEh+l96AUREVHtwSCPiBSq8wI4EIuA+xKHSpKEQW2SkBwT7r4/HpdQ8K3z7oqIjO/eWGh/ZTDle1hq0Otw+5BWDv3yvE9A5uT5MFzTZHD/r0q0KAsREVF9wiCPiEJG7RpfpKpl8Cpfiq31pmwZGC9e3dVejl/zvI4BruCSBiKqr7qmeOGVmRe3Qce0WEzqm+G2TXyk0dvuERER1XmsrklEbgW7EIvPi1cHKcaricVWnIkGuF4HeQIZMdGsmfji7toN7x7WGnd7KH7SNT0e9wxrjaYJkWInJSIiqgcY5BFRyARiuGYg+bvIugh/w2YphOMvxIdrVm+0fN/Fbar1fERERDUdh2sSUcj4uoSCY5NATsnyZiihr/ztr2i/ghFoaXXdcZkDneB/lhZJUX71h4iIiNQxk0dEbtXUmhb+LqEQynP4e9jg9cvzcbXeD45Dez0Nw9029xKUVpgRGx4m3D8iIiISxyCPiELG13hFdM6XP4Jf2sWVSEytyGJqtRN8AJFGPYrKzOjfsqHHtlpzNC2Wqp/1Hk4eFxEGgAEeERFRsDDIIyK3gp3IUx2uKbBfdcz5ClTC7M9/D8Oxc8W4cv7ayuP6d2DH3bUyZqJnWXHvIPy06wwm9mrqsa1oZldXHVE4ERERucUgj4jcCvaluq/HV66TF5xQNFCBZHJMuMt6fO7P6ZlOcN6g6LOSnhCJqRc2F2qrFZ+2S41B72YJSBJYBoKIiIiCi0EeEbkV7ExeIOaXBW3eoKeuBSEC9vahaD19jeMjsP14nl/9cRYf4X6IpU4nYclt/QJ6PiIiIvINq2sSUbW4Y2hLAMCNFzZDWlw42qbE+Dy3rjqWUAiWXs0S/NrfMajVGq75+BUdMbpzI3xySx+/zgcAz17ZGZd2TMU1vT0P6SQiIqLQYyaPiKrFA5e0xdiujdEyKRr/HtUeOklSz+QJRHC1aU6ezepZQ7DjRB5GdU716zh6h8g42uT+T3hyTDjemNTdr3PZTOzVVGjOHhEREdUMDPKIqFpIkoTWKTEAAJ2fQVo1rKAgMFrTu05kNIxCRkPtdeFeu7Ybbnh/A+aM6eC2jdGgw+vXdUO52YIGUUav+kBERET1A4M8Iqp1qmNYZXWsxeesb4uG2PXEpTDotUfSj+mSVk09IiIiotqozs3Ju/zyy9G0aVOEh4ejUaNGuOGGG3Dy5MlQd4uoVqrXi6EH/QzqPAV4RERERJ7UuauJoUOHYsmSJdizZw+++OILHDhwABMmTAh1t4gogEQXBPfrHFzqjYiIiGqpOjdc87777rP/nJGRgYcffhhjx45FeXk5wsLcl/8motqjeubkMcojIiKi2qnOZfIcnTt3Dh9//DH69+/PAI+oDnEcrikHaUwpM3lERERUW9XJIO+hhx5CVFQUGjZsiKNHj+Lrr7/WbF9aWor8/HzFFxGRFjnoS8UTERER+aZWBHkPP/wwpMo1tdx97d69295+1qxZ2LJlC1auXAm9Xo/Jkydrfto/b948xMXF2b/S09Or42ER1Xj1LZBh9o6IiIjqgloxJ+/+++/H1KlTNdu0aNHC/nNiYiISExPRpk0btG/fHunp6Vi/fj369eunuu/s2bMxc+ZM++38/HwGekS1RKgKr1zSIRVvrz6IxGiuVUdEREQ1S60I8pKSkpCUlOTTvhaLBYB1SKY7JpMJJpPJp+MTUf3UI6MBVt43CI3iwkPdFSIiIiKFWhHkidqwYQM2btyIAQMGoEGDBjhw4AAeffRRtGzZ0m0Wj4jIV21SYkLdBSIiIiIXtWJOnqjIyEh8+eWXGDZsGNq2bYtp06ahS5cuWL16NTN1RD6ob8sIxIZXVeHVc4IeERER1VJ1KpPXuXNn/Pzzz6HuBlGtd//FbbBo/RHMvKRNtZ+7Y1osth3LFd8hgJPyGkQZ8db13WE06GDQ16nPwIiIiKgeqVNBHhEFxl3DWuPOi1op1qNz58WrL8DMJdvwyjVdA3Luh0e2Q2x4GC67oJHYDgFOuF3aSfC8RERERDUUgzwiUiUS4AHA+O5NMKZLGoyGwGS+YsPD8PDIduI71K9VHoiIiIg84ngkIvJboAI8IiIiIvIfr8yIiIiIiIjqEAZ5RFSrJcWwci4RERGRI87JI6Ja6a3re+B4ThE6NY4LdVeIiIiIahQGeURUK13aKTXUXSAiIiKqkThck4iIiIiIqA5hkEdERERERFSHMMgjIiIiIiKqQxjkERERERER1SEM8oiIiIiIiOoQBnlERERERER1CIM8IiIiIiKiOoTr5KmQZRkAkJ+fH+KeEBERERERWdniE1u84g6DPBUFBQUAgPT09BD3hIiIiIiISKmgoABxcXFu75dkT2FgPWSxWHDy5EnExMRAkqRQdweANWpPT0/HsWPHEBsbG+ruUJDwda4/+FrXH3yt6w++1vUDX+f6oya+1rIso6CgAGlpadDp3M+8YyZPhU6nQ5MmTULdDVWxsbE15k1GwcPXuf7ga11/8LWuP/ha1w98neuPmvZaa2XwbFh4hYiIiIiIqA5hkEdERERERFSHMMirJUwmE+bOnQuTyRTqrlAQ8XWuP/ha1x98resPvtb1A1/n+qM2v9YsvEJERERERFSHMJNHRERERERUhzDIIyIiIiIiqkMY5BEREREREdUhDPKIiIh89Ouvv0KSJOTm5oa6K0RERHYsvEJERCRoyJAh6Nq1K15++WUAQFlZGc6dO4eUlBRIkhTazhEREVUyhLoDREREtZXRaERqamqou0FERKTA4ZpEREQCpk6ditWrV+OVV16BJEmQJAkLFy5UDNdcuHAh4uPj8e2336Jt27aIjIzEhAkTUFRUhA8++ADNmjVDgwYNcPfdd8NsNtuPXVpaigceeACNGzdGVFQU+vTpg19//TU0D5SIiGo9ZvKIiIgEvPLKK9i7dy86deqEJ554AgDw999/u7QrKirCq6++is8++wwFBQUYP348xo0bh/j4eHz//fc4ePAgrrzySlx44YWYOHEiAODOO+/EP//8g88++wxpaWn46quvcOmll2LHjh1o3bp1tT5OIiKq/RjkERERCYiLi4PRaERkZKR9iObu3btd2pWXl2P+/Plo2bIlAGDChAlYtGgRzpw5g+joaHTo0AFDhw7FL7/8gokTJ+Lo0aNYsGABjh49irS0NADAAw88gOXLl2PBggV45plnqu9BEhFRncAgj4iIKIAiIyPtAR4ApKSkoFmzZoiOjlZsy8zMBADs2LEDZrMZbdq0URyntLQUDRs2rJ5OExFRncIgj4iIKIDCwsIUtyVJUt1msVgAAOfPn4der8fmzZuh1+sV7RwDQyIiIlEM8oiIiAQZjUZFwZRA6NatG8xmMzIzMzFw4MCAHpuIiOonVtckIiIS1KxZM2zYsAGHDx9Gdna2PRvnjzZt2mDSpEmYPHkyvvzySxw6dAh//vkn5s2bh++++y4AvSYiovqGQR4REZGgBx54AHq9Hh06dEBSUhKOHj0akOMuWLAAkydPxv3334+2bdti7Nix2LhxI5o2bRqQ4xMRUf0iybIsh7oTREREREREFBjM5BEREREREdUhDPKIiIiIiIjqEAZ5REREREREdQiDPCIiIiIiojqEQR4REREREVEdwiCPiIiIiIioDmGQR0REREREVIcwyCMiIiIiIqpDGOQRERERERHVIQzyiIiIiIiI6hAGeURERERERHUIgzwiIiIiIqI6hEEeERERERFRHcIgj4iIiIiIqA5hkEdERERERFSHMMgjIiIiIiKqQxjkERERERER1SEM8oioTmjWrBmmTp0a6m7UKlOnTkWzZs38Po4kSXjssce83u+xxx6DJEl+n78mWbhwISRJwuHDh73ed8iQIejUqZPHdocPH4YkSVi4cKH3HQyCY8eOITw8HGvWrAnK8QP1Pq0pvPlb5c/7KRD69u2LBx98MCTnJiL/MMgjohptx44dmDBhAjIyMhAeHo7GjRvj4osvxmuvvRbqrtU5RUVFeOyxx/Drr7+63Pf999/7FMh565lnnsGyZcuCfh4KnCeeeAJ9+vTBhRdeGOquUIA99NBDeOONN3D69OlQd4WIvMQgj4hqrLVr16Jnz57Ytm0bbrnlFrz++uu4+eabodPp8Morr4S6e7Xeu+++iz179thvFxUV4fHHH3cb5D3++OOqxykuLsYjjzzi9fkfeeQRFBcXK7bV9iDvhhtuQHFxMTIyMkLdlWqRlZWFDz74ALfddlvQzuH8Pq1PQv1+uuKKKxAbG4s333wzJOcnIt8ZQt0BIiJ3nn76acTFxWHjxo2Ij49X3JeZmRmaTtUAhYWFiIqK8vs4YWFhAegNEB4e7tN+BoMBBkPd+Ddke030ej30en2ou1NtPvroIxgMBlx22WVBO0eg3qeiZFlGSUkJIiIiqvW8jkL9fioqKkJkZCR0Oh0mTJiADz/8EI8//nidG15NVJcxk0dENdaBAwfQsWNHlwAPAJKTkz3uf/DgQVx11VVISEhAZGQk+vbti++++85+vyzLSExMxMyZM+3bLBYL4uPjodfrkZuba9/+7LPPwmAw4Pz5827PZ5s/89tvv2H69Olo2LAhYmNjMXnyZOTk5Li0/+GHHzBw4EBERUUhJiYGo0ePxt9//61oM3XqVERHR+PAgQMYNWoUYmJiMGnSJADWuXB33nknli5dig4dOiAiIgL9+vXDjh07AABvv/02WrVqhfDwcAwZMsRlXo/jXKfDhw8jKSkJAOwXc7a5dlOnTsUbb7xhP6fty8ZxTt7nn38OSZKwevVql8f79ttvQ5Ik7Ny5E4DrnDxJklBYWIgPPvjAfo6pU6fil19+gSRJ+Oqrr1yO+cknn0CSJKxbt071Ndm0aRMkScIHH3zgct+KFSsgSRK+/fZbAMCRI0dw++23o23btoiIiEDDhg1x1VVXuTxvttd59erVuP3225GcnIwmTZoo7nPc5+uvv8bo0aORlpYGk8mEli1b4sknn4TZbFbt8+bNm9G/f39ERESgefPmeOutt1TbOdu9ezcmTJiAhIQEhIeHo2fPnvjmm2887jdkyBDF6+r45Wne37Jly9CnTx9ER0crtu/btw9XXnklUlNTER4ejiZNmuCaa65BXl6eot1HH32EHj16ICIiAgkJCbjmmmtw7NgxRRvnOXki/XU331Pt9WnWrBnGjBmDFStWoGfPnoiIiMDbb78NAMjNzcW9996L9PR0mEwmtGrVCs8++ywsFouHZ9X69+Wpp55CkyZNEBkZiaFDh7r8fjv2SeT9NGbMGLRo0UL1fP369UPPnj0V20SeX9tc0M2bN2PQoEGIjIzEv/71L/v9F198MY4cOYKtW7d6fMxEVHPUjY9QiahOysjIwLp167Bz506hghSOzpw5g/79+6OoqAh33303GjZsiA8++ACXX345Pv/8c4wbNw6SJOHCCy/Eb7/9Zt9v+/btyMvLg06nw5o1azB69GgAwO+//45u3bq5XMyqufPOOxEfH4/HHnsMe/bswfz583HkyBH8+uuv9gvPRYsWYcqUKRgxYgSeffZZFBUVYf78+RgwYAC2bNmiuKitqKjAiBEjMGDAADz//POIjIy03/f777/jm2++wR133AEAmDdvHsaMGYMHH3wQb775Jm6//Xbk5OTgv//9L2666Sb8/PPPqn1OSkrC/PnzMWPGDIwbNw7jx48HAHTp0gWFhYU4efIkVq1ahUWLFmk+9tGjRyM6OhpLlizB4MGDFfctXrwYHTt2dPtaLlq0CDfffDN69+6NW2+9FQDQsmVL9O3bF+np6fj4448xbtw4xT4ff/wxWrZsiX79+qkes2fPnmjRogWWLFmCKVOmuPSnQYMGGDFiBABg48aNWLt2La655ho0adIEhw8fxvz58zFkyBD8888/iucdAG6//XYkJSVhzpw5KCwsdPucLFy4ENHR0Zg5cyaio6Px888/Y86cOcjPz8dzzz2naJuTk4NRo0bh6quvxrXXXoslS5ZgxowZMBqNuOmmm9ye4++//8aFF16Ixo0b4+GHH0ZUVBSWLFmCsWPH4osvvnB53hz9+9//xs0336zY9tFHH2HFihWaH6aUl5dj48aNmDFjhmJ7WVkZRowYgdLSUtx1111ITU3FiRMn8O233yI3NxdxcXEArJn6Rx99FFdffTVuvvlmZGVl4bXXXsOgQYOwZcsW1Q93/Omvlj179uDaa6/F9OnTccstt6Bt27YoKirC4MGDceLECUyfPh1NmzbF2rVrMXv2bJw6dQovv/yy5jHnzJmDp556CqNGjcKoUaPw119/4ZJLLkFZWZlqe5H308SJEzF58mRs3LgRvXr1sm8/cuQI1q9fr3g/efP8nj17FiNHjsQ111yD66+/HikpKfb7evToAQBYs2YNunXr5umpJKKaQiYiqqFWrlwp6/V6Wa/Xy/369ZMffPBBecWKFXJZWZlL24yMDHnKlCn22/fee68MQP7999/t2woKCuTmzZvLzZo1k81msyzLsvzcc8/Jer1ezs/Pl2VZll999VU5IyND7t27t/zQQw/JsizLZrNZjo+Pl++77z7N/i5YsEAGIPfo0UPRx//+978yAPnrr7+29yM+Pl6+5ZZbFPufPn1ajouLU2yfMmWKDEB++OGHXc4HQDaZTPKhQ4fs295++20ZgJyammp/TLIsy7Nnz5YBKNpOmTJFzsjIsN/OysqSAchz5851Odcdd9whu/uX4bzPtddeKycnJ8sVFRX2badOnZJ1Op38xBNP2LfNnTvX5ZhRUVGK19Gx/yaTSc7NzbVvy8zMlA0Gg2p/nfcNCwuTz507Z99WWloqx8fHyzfddJN9W1FRkcu+69atkwHIH374oX2b7XUeMGCA4jE63uf4PKsdd/r06XJkZKRcUlJi3zZ48GAZgPzCCy8o+tm1a1c5OTnZ/p46dOiQDEBesGCBvd2wYcPkzp07K45nsVjk/v37y61bt9Z6elysWbNGDgsLUzw3avbv3y8DkF977TXF9i1btsgA5KVLl7rd9/Dhw7Jer5effvppxfYdO3bIBoNBsd35fSrSX7X3liyrvz4ZGRkyAHn58uWKtk8++aQcFRUl7927V7H94YcflvV6vXz06FG3fcrMzJSNRqM8evRo2WKx2Lf/61//kgEo3uPevJ/y8vJkk8kk33///Yp2//3vf2VJkuQjR47Isuzd82t737311ltuH4/RaJRnzJjh9n4iqnk4XJOIaqyLL74Y69atw+WXX45t27bhv//9L0aMGIHGjRt7HIb2/fffo3fv3hgwYIB9W3R0NG699VYcPnwY//zzDwBg4MCBMJvNWLt2LQBrZmzgwIEYOHAgfv/9dwDAzp07kZubi4EDBwr1+9Zbb1XMI5oxYwYMBgO+//57AMCqVauQm5uLa6+9FtnZ2fYvvV6PPn364JdffnE5pnO2xGbYsGGKrF+fPn0AAFdeeSViYmJcth88eFDoMfhj4sSJyMzMVBRw+fzzz2GxWDBx4kSfjjl58mSUlpbi888/t29bvHgxKioqcP3113vsT3l5Ob788kv7tpUrVyI3N1fRH8c5WOXl5Th79ixatWqF+Ph4/PXXXy7HveWWW4TmSzket6CgANnZ2Rg4cCCKioqwe/duRVuDwYDp06fbbxuNRkyfPh2ZmZnYvHmz6vHPnTuHn3/+GVdffbX9+NnZ2Th79ixGjBiBffv24cSJEx77CQCnT5/GhAkT0LVrV4/FNs6ePQsAaNCggWK7LVO3YsUKFBUVqe775ZdfwmKx4Oqrr1b8DqSmpqJ169aqvwP+9ldL8+bN7Rldm6VLl2LgwIFo0KCBoo/Dhw+H2WxWjABw9uOPP6KsrAx33XWXYtjovffe63YfkfdTbGwsRo4ciSVLlkCWZfv2xYsXo2/fvmjatCkA759fk8mEG2+80e15bc8BEdUeDPKIqEbr1asXvvzyS+Tk5ODPP//E7NmzUVBQgAkTJtgDNTVHjhxB27ZtXba3b9/efj8AdO/eHZGRkfaAzhbkDRo0CJs2bUJJSYn9PseAUUvr1q0Vt6Ojo9GoUSP7vJp9+/YBAC666CIkJSUpvlauXOlSVMZgMNjn6DizXdTZ2C6w09PTVberzQ0MtEsvvRRxcXFYvHixfdvixYvRtWtXtGnTxqdjtmvXDr169cLHH39s3/bxxx+jb9++aNWqlea+F1xwAdq1a+fSn8TERFx00UX2bcXFxZgzZ459/lViYiKSkpKQm5vrMpcMsAYGIv7++2+MGzcOcXFxiI2NRVJSkj0wdT5uWlqaS1Ed23Pmbq20/fv3Q5ZlPProoy7vp7lz5wIQK1RUUVGBq6++GmazGV9++SVMJpPQ43MMNgDr8zJz5ky89957SExMxIgRI/DGG28oHuu+ffsgyzJat27t0uddu3YFtb9q1F7Lffv2Yfny5S79Gz58OADt59T298X5b0FSUpJLUKzVBzUTJ07EsWPH7PNQDxw4gM2bNys+sPD2+W3cuDGMRqPbc8qyzKIrRLUM5+QRUa1gNBrRq1cv9OrVC23atMGNN96IpUuX2i9ifRUWFoY+ffrgt99+w/79+3H69GkMHDgQKSkpKC8vx4YNG/D777+jXbt29sIk/rIVbVi0aBFSU1Nd7neuOGkymaDTqX8m5+6Tf3fbnS/Ig8FkMmHs2LH46quv8Oabb+LMmTNYs2YNnnnmGb+OO3nyZNxzzz04fvw4SktLsX79erz++utC+06cOBFPP/00srOzERMTg2+++QbXXnut4rm+6667sGDBAtx7773o168f4uLiIEkSrrnmGtVCGyLVF3NzczF48GDExsbiiSeeQMuWLREeHo6//voLDz30kFABD09sx3jggQdcslE2ngJhAJg1axbWrVuHH3/80e2HCo4aNmwIQP2DgxdeeAFTp07F119/jZUrV+Luu+/GvHnzsH79ejRp0gQWiwWSJOGHH35Qfa+KzH3V6q+7gMRdsRu119JiseDiiy92uxi4rx9YuCNazfOyyy5DZGQklixZgv79+2PJkiXQ6XS46qqr7G28fX49nTs3NxeJiYlC/SOimoFBHhHVOrYKcqdOnXLbJiMjQ3VtLdvwOMd1pwYOHIhnn30WP/74IxITE9GuXTtIkoSOHTvi999/x++//44xY8YI92/fvn0YOnSo/fb58+dx6tQpjBo1CoC1mAhgrRBqywrUBFqf1Hv7Kf7EiRPxwQcf4KeffsKuXbsgy7LQUE2t81xzzTWYOXMmPv30UxQXFyMsLEx4+OfEiRPx+OOP44svvkBKSgry8/NxzTXXKNp8/vnnmDJlCl544QX7tpKSEkWVVW/9+uuvOHv2LL788ksMGjTIvv3QoUOq7U+ePOmyRMbevXsBQDEs15Gt2mJYWJjP76fPPvsML7/8Ml5++WWXgjnuNG3aFBEREW4fS+fOndG5c2c88sgjWLt2LS688EK89dZbeOqpp9CyZUvIsozmzZv7FCx56q8tW5abm6soMGLLsIlo2bIlzp8/79Nzavv7sm/fPkU1zKysLL+z6VFRURgzZgyWLl2KF198EYsXL8bAgQORlpam6Ls/z6+jEydOoKyszD4KgohqBw7XJKIa65dfflHNPNnmtqkNx7QZNWoU/vzzT0Vp/cLCQrzzzjto1qwZOnToYN8+cOBAlJaW4uWXX8aAAQPsgcbAgQOxaNEinDx5Ung+HgC88847KC8vt9+eP38+KioqMHLkSADAiBEjEBsbi2eeeUbRziYrK0v4XIFkqx6pFtTYgg7RgGf48OFISEjA4sWLsXjxYvTu3VtoOFpUVJTbcyQmJmLkyJH46KOP8PHHH+PSSy8Vzi60b98enTt3tvenUaNGiqALsGY/nd9vr732mtvsjwhbFsXxuGVlZW7nj1VUVNjL99vavv3220hKSrJXOXSWnJyMIUOG4O2331b94MPT+2nnzp24+eabcf311+Oee+7x+JhswsLC0LNnT2zatEmxPT8/HxUVFYptnTt3hk6nQ2lpKQBg/Pjx0Ov1ePzxx12ec1mW7fP9fO2v7YMUx3lztuU5RF199dVYt24dVqxY4XJfbm6uy2N0NHz4cISFheG1115TPD5PFTlFTZw4ESdPnsR7772Hbdu2uXzY4c/z68w2F7R///7+d5yIqg0zeURUY911110oKirCuHHj0K5dO5SVlWHt2rVYvHgxmjVrplko4OGHH8ann36KkSNH4u6770ZCQgI++OADHDp0CF988YVi+GO/fv1gMBiwZ88ee+l+ABg0aBDmz58PAF4FeWVlZRg2bBiuvvpq7NmzB2+++SYGDBiAyy+/HIC1eML8+fNxww03oHv37rjmmmuQlJSEo0eP4rvvvsOFF14oPAwxkCIiItChQwcsXrwYbdq0QUJCAjp16oROnTrZA4y7774bI0aMgF6vd8mEOQoLC8P48ePx2WefobCwEM8//7xQH3r06IEff/wRL774ItLS0tC8eXN70RjAOmRzwoQJAIAnn3zSq8c3ceJEzJkzB+Hh4Zg2bZrLENgxY8Zg0aJFiIuLQ4cOHexDAW3DEn3Rv39/NGjQAFOmTMHdd98NSZKwaNEit8Nm09LS8Oyzz+Lw4cNo06YNFi9ejK1bt+Kdd97RXBT8jTfewIABA9C5c2fccsstaNGiBc6cOYN169bh+PHj2LZtm9t9bb9HgwYNwkcffeTSf3frsgHAFVdcgX//+9/Iz89HbGwsAODnn3/GnXfeiauuugpt2rRBRUUFFi1aBL1ejyuvvBKANQh76qmnMHv2bBw+fBhjx45FTEwMDh06hK+++gq33norHnjgAZ/7e8kll6Bp06aYNm0aZs2aBb1ej//7v/+z/56JmDVrFr755huMGTMGU6dORY8ePVBYWIgdO3bg888/x+HDh91+yJCUlIQHHnjAvqTJqFGjsGXLFvzwww8BGfZoWzPzgQceUDyvNv48v85WrVqFpk2bcvkEotqmmqt5EhEJ++GHH+SbbrpJbteunRwdHS0bjUa5VatW8l133SWfOXNG0dZ5CQVZluUDBw7IEyZMkOPj4+Xw8HC5d+/e8rfffqt6rl69eskA5A0bNti3HT9+XAYgp6enC/XXVu589erV8q233io3aNBAjo6OlidNmiSfPXvWpf0vv/wijxgxQo6Li5PDw8Plli1bylOnTpU3bdpkbzNlyhQ5KipK9XwA5DvuuEOxzVZe/7nnnnM5F5zK2quVpl+7dq3co0cP2Wg0KpZGqKiokO+66y45KSlJliRJUZ4eUF92YdWqVTIAWZIk+dixYy73q5W53717tzxo0CA5IiLCpdS8LFuXFGjQoIEcFxcnFxcXqz4v7uzbt08GIAOQ//jjD5f7c3Jy5BtvvFFOTEyUo6Oj5REjRsi7d+92eW/ZXueNGze6HEOtRP+aNWvkvn37yhEREXJaWpp9KRAA8i+//GJvN3jwYLljx47ypk2b5H79+snh4eFyRkaG/PrrryvOobaEgixb3++TJ0+WU1NT5bCwMLlx48bymDFj5M8//1zzebEtIaD25XwOZ2fOnJENBoO8aNEi+7aDBw/KN910k9yyZUs5PDxcTkhIkIcOHSr/+OOPLvt/8cUX8oABA+SoqCg5KipKbteunXzHHXfIe/bssbdxfp+K9nfz5s1ynz59ZKPRKDdt2lR+8cUX3S6hMHr0aNXHV1BQIM+ePVtu1aqVbDQa5cTERLl///7y888/r7qUiyOz2Sw//vjjcqNGjeSIiAh5yJAh8s6dO/1+P9lMmjRJBiAPHz7cbR9Enl/b+87dY2jUqJH8yCOPaD5WIqp5JFmuhln4RET1wMKFC3HjjTdi48aN9nmDFFgVFRVIS0vDZZddhvfffz/U3SEA06ZNw969e+1VaKnuWLZsGa677jocOHAAjRo1CnV3iMgLnJNHRES1xrJly5CVlYXJkyeHuitUae7cudi4cSPWrFkT6q5QgD377LO48847GeAR1UKck0dERDXehg0bsH37djz55JPo1q2bcAVICr6mTZuipKQk1N2gIHAsXEVEtQszeUREVOPNnz8fM2bMQHJyMj788MNQd4eIiKhG45w8IiIiIiKiOoSZPCIiIiIiojqEQR4REREREVEdwsIrKiwWC06ePImYmBhIkhTq7hAREREREUGWZRQUFCAtLQ06nft8HYM8FSdPnkR6enqou0FEREREROTi2LFjaNKkidv7GeSpiImJAWB98mJjY0PcGyIiIiIiIiA/Px/p6en2eMUdBnkqbEM0Y2NjGeQREREREVGN4mlKGQuvEBERERER1SEM8oiIiIiIiOoQBnlERERERER1CIM8IiIiIiKiOoRBHhERERERUR3CII+IiIiIiKgO4RIKRERERCEmyzJ2nSrA6r1ZyD5fimYNIzGmSxoaRBlD3TUiqoUY5BERERGF0J+HzuGZ73dh67FcxfZnl+/B45d3xJU9moSmY0RUazHIIyIiIgqBQ9mFmPf9Lqz85wwAwGjQYXCbJGQkROKP/dnYfboA9y/dhhO5xbjrolYeFz8mIrJhkEdERERUjUrKzZj/6wHM//UAyswW6HUSrumVjnuHt0FSjAkAYLbIeGnVXrz+y368uGovDHoJtw9pFeKeE1FtwSCPiIiIqJqsPZCNR77aiYPZhQCAwW2S8Mjo9midEqNop9dJeGBEW8RGGPDM97vx3+V7kBobjvHdOXSTiDxjkEdEREQUZGfPl+Lp73fhy79OAACSYkx47LKOGNU5VXMY5q2DWiKroBTv/n4ID36+HU0aRKJ384Tq6jYR1VJcQoGIiIgoSGRZxpKNxzDsxdX48q8TkCRgcr8M/HT/YIzu0khont3ske0xunMjVFhkzPhoM07kFldDz4moNmMmj4iIiCgIjp0rwuwvd+CP/dkAgPaNYjFvfGd0TY/36jg6nYTnruqCQ9mF+OdUPqYv2oSl0/sjwqgPQq+JqC5gJo+IiIgogCwWGYvWH8GlL/+GP/ZnIzxMh3+Naof/3Xmh1wGeTaTRgHcm90BClBE7T+TjoS+2Q5blwHaciOoMBnlEREREAXLsXBEmvbcBjy7bicIyM3o1a4Af7hmEWwe1hEHv32VXkwaReOO67tDrJHyz7STe+e1ggHpNRHUNgzwiIiIiP1ksMj5YexgjXv4N6w6eRUSYHnMv64DFt/ZD88SogJ2nX8uGmHtZBwDAf5bvxq97MgN2bCKqOxjkEREREfnhyNlCXPPuesz95m8UlZnRp3kClt87EDde2Bw6XeAXML+hbwYm9kyHLAN3fboFhyqXYyAismGQR0REROQDWZbx8YYjuPTl3/HnoXOINOrx5BUd8ektfZHRMHDZO2eSJOGJsR3RvWk8CkoqcMuHm1BQUh608xFR7cMgj4iIiMhLWQWlmPbBJvz7q50oLjejX4uGWHHvINzQr1lQsnfOTAY93rq+B1JiTdifeR73Ld4Gi4WFWIjIikEeERERkRdW/n0aI17+DT/vzoTRoMMjo9vj45v7ID0hslr7kRwbjreu7wGjXocfd53Byz/tq9bzE1HNxSCPiIiISMD50go89Pl23LpoM84VlqFdagz+d+cA3DywRbVk79R0a9oAT4/rBAB49ad9WL7zVEj6QUQ1CxdDJyIiIvJg85Ec3Ld4K46eK4IkAbcObIGZl7SByRD6Bcmv6pmOv0/mY+Haw5i5ZBuaJ0ajbWpMqLtFRCHETB4RERGRG+VmC15YuQdXvbUWR88VoXF8BD69pS9mj2pfIwI8m3+Pbo9+LRqiqMyMWz7chNyislB3iYhCiEEeERERkYoDWedx5fy1eO3n/bDIwPhujfHDvQPRt0XDUHfNRZhehzcmdUeTBhE4eq4Id326BRVmS6i7RUQhUmOCvP/85z+QJAn33nuvZrulS5eiXbt2CA8PR+fOnfH9998r7pdlGXPmzEGjRo0QERGB4cOHY98+TkQmIiIiMbIsY+mmY7jstT+w/Xge4iLC8MZ13fHixK6IDQ8LdffcSogy4t3JPRERpsfv+7Lx3xV7Qt0lIgqRGhHkbdy4EW+//Ta6dOmi2W7t2rW49tprMW3aNGzZsgVjx47F2LFjsXPnTnub//73v3j11Vfx1ltvYcOGDYiKisKIESNQUlIS7IdBREREtVxBSTnuXbwVsz7fjqIyM/q3tC6NMLpLo1B3TUj7RrF4/qoLAADv/HYQy7acCHGPiCgUJFmWQ7qoyvnz59G9e3e8+eabeOqpp9C1a1e8/PLLqm0nTpyIwsJCfPvtt/Ztffv2RdeuXfHWW29BlmWkpaXh/vvvxwMPPAAAyMvLQ0pKChYuXIhrrrlGqE/5+fmIi4tDXl4eYmNj/X6MREREVPNtP56Luz7dgiNni6DXSZh5cRvcNrgl9CGqnOmP51bsxhu/HIDJoMPnt/VH5yZxoe4SEQWAaJwS8kzeHXfcgdGjR2P48OEe265bt86l3YgRI7Bu3ToAwKFDh3D69GlFm7i4OPTp08fehoiIiMiRxSLj3d8O4sr5a3HkrLW4ypLpfXHH0Fa1MsADgPsvbouL2iWjtMKC6Ys2Ift8aai7RETVKKRLKHz22Wf466+/sHHjRqH2p0+fRkpKimJbSkoKTp8+bb/fts1dGzWlpaUoLa3645efny/UHyIiIqrdss+X4v4l27B6bxYAYGSnVPxnfBfERdbcuXcidDoJL1/TFWPfWIODWYW4/aO/8NHNfWA0hPzzfSKqBiH7TT927BjuuecefPzxxwgPDw9VNwAA8+bNQ1xcnP0rPT09pP0hIiKi4FuzPxsjX/kdq/dmwWTQ4elxnfDmpO61PsCziQ0Pw7uTeyLGZMCfh8/hiW//DnWXiKiahCzI27x5MzIzM9G9e3cYDAYYDAasXr0ar776KgwGA8xms8s+qampOHPmjGLbmTNnkJqaar/fts1dGzWzZ89GXl6e/evYsWP+PjwiIiKqocwWGS+t2ovr39+ArIJStE6Oxjd3DsCkPhmQpNo5PNOdlknReOXarpAk4KP1R/HR+iOh7hIRVYOQBXnDhg3Djh07sHXrVvtXz549MWnSJGzduhV6vesCo/369cNPP/2k2LZq1Sr069cPANC8eXOkpqYq2uTn52PDhg32NmpMJhNiY2MVX0RERFT3ZBWUYvL/bcArP+2DLAPX9ErHN3cOQNvUmFB3LWguapeCBy5pCwCY+83fWLs/O8Q9IqJgC9mcvJiYGHTq1EmxLSoqCg0bNrRvnzx5Mho3box58+YBAO655x4MHjwYL7zwAkaPHo3PPvsMmzZtwjvvvAMA9nX2nnrqKbRu3RrNmzfHo48+irS0NIwdO7ZaHx8RERHVLOsPnsVdn25BVkEpIsL0eGZ8J4zr1iTU3aoWtw9piX1nCrBs60nM+PgvLLvjQjRPjAp1t4goSEJaeMWTo0ePQqerSjb2798fn3zyCR555BH861//QuvWrbFs2TJFsPjggw+isLAQt956K3JzczFgwAAsX7485PP+iOj/27vvuCrr/o/jr8PeIIIMcQKiqLgHmi1xlJXmSM3EyqzM9va+73Z32rjrbppmjoZSemtbG+ZIBSfurSiiIIqy9znX7w+KX6SWCXgx3s/H4zweec73XH6uvh643uf6DhERc9hsBtNWHuQ/P+zFZkB4Iw/eG9OZ8IC6e/fujywWC1OHRXHkdD6JyZmMn7OBxff0rjPzD0WkItP3yauJtE+eiIhI3XA6r5iHPttSvnrmsM4hvDCkLW5ONfp77mpzMqeIwe+s5nhWIZeF+THntm442GvFTZHaotbskyciIiJSHbanZHHdW/+/euYrw6P4z00d6m3AA/D3dGbmuG64Odmz+sApXvhml9kliUg1UMgTERGROmfR5hSGv7+W41mFtPBz58t7e3NTV22RBBAZ7MUbIzsCMDf+CB/HHza1HhGpegp5IiIiUmeUWG08+9VOHv58K0WlNvq2bsSX9/amdaCmX/zegLaBPD6wbMXNZ7/exer9WnFTpC5RyBMREZE64VRuEbfMXMectYcBuL9vOB/EdsXLRYuLnMvEK0IZ2qkxVpvBPZ9u4tDJXLNLEpEqopAnIiIitd62lExueHs165JO4+HswPSxXXi4Xyvs7OrW5uZVyWKx8NLQ9nRu6kN2YSnj524kK7/E7LJEpAoo5ImIiEit9r9NKQx/P57jWYW09Hfni0m9GNA20OyyagUXR3umj+1KYx9Xkk7lcc+8TZRYbWaXJSKVpJAnIiIitdJv8+8eWbCV4lIbMW0a8cWk3oQ1qj/731WFshU3u+LmZM+aAxk89/VOtMOWSO2mkCciIiK1zqncIsb8bv7dA33DmTFW8+8uVpsgL94c1QmLBT5JSGb2msNmlyQilaCQJyIiIrXKtpRMrn97Net/nX/3QWxXHtL8u0rrFxnAkwNbA/DCt7v4YWeayRWJyMVSyBMREZFaY8HGowx/P57U8vl3vekXGWB2WXXGnZe3ZHT3phgGPBC3he0pWWaXJCIXQSFPREREarwSq41nvtzBYwu3/Tr/LoAvJ/UmrJGH2aXVKRaLhecHt6VPuB8FJVbGz93A8cwCs8sSkb9JIU9ERERqtJM5RYz5YB1z448A8FBMK2aM7YKn5t9VC0d7O94d05mIAE/Sc4q4fc4Gcgq1tYJIbaKQJyIiIjXWlqO/zr87fBpPZwdmxnblgZhwzb+rZl4ujsy6rRv+ns7sScvh3nmJlGprBZFaQyFPREREaqTPNx7lpunxpGUXEurvzhf39iZG8+8umcY+rnw4risujnas3HeSZ77S1goitYVCnoiIiNQoxaU2nvpiB4//Ov+uf2QAX0zqTai/5t9dalEhPuVbK3y6LpkPVyeZXZKIXACFPBEREakxTuYUMWZmAh8nHMFigYf7teL9WzT/zkwD2gbyz2vbAPDv73bzvbZWEKnxFPJERESkRkhMPsP1b69mw+Ez5fPv7u+r+Xc1wfjLWjC2Z7Nft1ZIZFtKptklicifUMgTERER032+8SgjpyeQll1IWCMPvry3N33baP5dTWGxWHjm+kiujPCnsMTG+LkbSTmTb3ZZInIeCnkiIiJimlKrjee/3lU2/85qY0Dbsvl3LTX/rsZxsLfj7dGdaB3oyclft1bIKtDWCiI1kUKeiIiImCIzv5jb5mxg1pqyxTweimnFtDFd8HB2MLkyOR9PF0dm3dqNAC9n9p3I5c6PNlJUajW7LBH5A4U8ERERueT2n8hhyLtr+GX/Kdyc7Hn/ls7a/66WCPZxZfat3fFwdmBd0mke+XwrNpu2VhCpSRTyRERE5JL6dlsqN7yzhsMZ+TT2ceV/E3sxsF2Q2WXJ3xAZ7MX0sV1wsLPwzbZUpi7dY3ZJIvI7CnkiIiJySVhtBlOX7GHSvM0UlFjpFdqQr+7tTZsgL7NLk4vQO8yPV0dEATBj1SFmr9EeeiI1hQa9i4iISLU7nVfMA3GJ/LL/FAB3XNaCyde2wV7DM2u1GzuFkJpVyCtL9/L8N7sI9HLhmva6KytiNoU8ERERqVZbjmYy6dPNHMsswMXRjleGd+CGDsFmlyVVZOIVoaRmFvJxwhEe+GwLfp7OdGvua3ZZIvWahmuKiIhItTAMg4/iD3PT9HiOZRbQvKEbX0zqrYBXx1gsFp69oS39IgMoLrVxx9yNHEjPMbsskXpNIU9ERESqXFZ+CRM/2czTX+6kuNRG/8gAvr7vMloHav5dXWRvZ+GtUZ3o1NSHrIISxs3aQHp2odllidRbCnkiIiJSpbanZHHDu6tZujMNBzsL/7y2DdPHdsHTxdHs0qQauTrZ8+G4brTwc+dYZgG3zt5AblGp2WWJ1Eumhrxp06YRFRWFl5cXXl5eREdHs2TJkvO2v/LKK7FYLGc9Bg0aVN7m1ltvPev1gQMHXorTERERqdfyikp55ssdDJ22hiO/bo/wxaTeTLi8JRaLFlipD3zdnZh7W3f8PJzYlZrNxE82UWK1mV2WSL1j6sIrISEhTJ06lfDwcAzDYO7cuQwePJjExETatm17VvtFixZRXFxc/ueMjAw6dOjAiBEjKrQbOHAgs2fPLv+zs7Nz9Z2EiIhIPZddWMKCjSl8knCEpFN5AAxsG8jLw6LwdtPdu/qmaUM3Zt3ajVEzEvhl/ymeWLiN10Z00Eb3IpeQqSHv+uuvr/Dnf//730ybNo2EhIRzhjxf34orNcXFxeHm5nZWyHN2diYwMLDqCxYREZFyhmGwav8pnvlyB4cz8oGyOzlvjupIn3B/k6sTM0WF+PDumM7cMXcjixKP4evuxD8HtdEdXZFLpMbMybNarcTFxZGXl0d0dPQFvefDDz9k1KhRuLu7V3h+xYoVNGrUiIiICCZOnEhGRkZ1lCwiIlJvFZZYefjzrYybtZ7DGfkEernw3A1tWfbwFQp4AsBVEY14ZVjZZukzVycxfdUhkysSqT9M3ydv+/btREdHU1hYiIeHB4sXLyYyMvIv37d+/Xp27NjBhx9+WOH5gQMHMnToUFq0aMHBgwf5xz/+wTXXXEN8fDz29vbnPFZRURFFRUXlf87Ozq7cSYmIiNRhJ3OKuPPjjSQmZ2JvZ2Fsz2ZMvDKUAC8Xs0uTGmZYlxBO5xXz7+92M3XJHnzdnbipaxOzyxKp8yyGYRhmFlBcXExycjJZWVksXLiQmTNnsnLlyr8MenfddRfx8fFs27btT9sdOnSI0NBQfvrpJ/r27XvONs8++yzPPffcWc9nZWXh5aWlnkVERH6zJy2b8XM2ciyzAG9XR94b05neYX5mlyU13JQlu5m+8hD2dhbev6UL/SIDzC5JpFbKzs7G29v7L3OK6cM1nZycCAsLo0uXLkyZMoUOHTrw5ptv/ul78vLyiIuLY/z48X95/JYtW+Ln58eBAwfO22by5MlkZWWVP44ePfq3z0NERKSu+2FnGkPfW8uxzAJa+Lmz+J5eCnhyQZ4c2JoRXUKw2gzunbeZ9UmnzS5JpE4zPeT9kc1mqzB08lwWLFhAUVERt9xyy18eLyUlhYyMDIKCgs7bxtnZuXwbh98eIiIiUsYwDN75eT93fbKJ/GIrvcMasmhiL1r6e5hdmtQSFouFKUPbE9OmEUWlNsbP3cDuVE2PEakupoa8yZMns2rVKg4fPsz27duZPHkyK1asYMyYMQDExsYyefLks9734YcfMmTIEBo2bFjh+dzcXB577DESEhI4fPgwy5YtY/DgwYSFhTFgwIBLck4iIiJ1SV5RKfd8upnXftiHYUBsdDPm3NadBu5OZpcmtYyDvR1vj+5Mt+YNyCksJXbWeo6ezje7LJE6ydSFV9LT04mNjSU1NRVvb2+ioqL4/vvv6devHwDJycnY2VXMoXv37mX16tX88MMPZx3P3t6ebdu2MXfuXDIzMwkODqZ///688MIL2itPRETkbzp4MpcJH23k0Mk8HO0tPHdDO27u0dTssqQWc3WyZ2ZsN0bOiGdPWg5jP1zHwom98PPQdZpIVTJ94ZWa6EInNIqIiNRVn284yrNf7yS/2EqAlzPvjelMl2a+f/1GkQtwIruQYdPWknKmgHaNvZg/oSeeLo5mlyVS49WahVdERESk5igssfLCN7t4/H/byC+20qOFL9/e30cBT6pUgJcLH4/vQUN3J3Ycy+bOjzZRWGI1uyyROkMhT0RERABYc+AUV7y6nA9XJwEw6apQ5k/oqaF0Ui1a+Lkz57buuDvZE38ogwfiEim12swuS6ROUMgTERGp59YePMWYmQnc8uE6TmQXEeTtwszYrjw2oDV2dhazy5M6rH2INx/EdsXJ3o7vd57g8YXbsNk0k0ikskxdeEVERETMk1NYwr+/3c1nG4/y2wz94V1CeHFIO1wc7c0tTuqNXmF+vHNzJyZ+uplFicfwcHHguRvaYrHoCwaRi6WQJyIiUg8dSM/lzo82cuhUHgAjuzZh4pWhNPdzN7kyqY/6tw3kPyM68NDnW/go/gieLg48NqC12WWJ1FoKeSIiIvXMqn0nmTRvMzmFpQR7u/DfUZ3o3kILq4i5hnRqTG5RKf/6YgfvLj+Ih7MjE68MNbsskVpJIU9ERKQe+Tj+MM9+vQurzaBrswa8P7aLFlaRGuOWns3ILSpl6pI9vLx0Dx4uDozt2czsskRqHYU8ERGReqDEauO5r3fySUIyUDb37t83tsPZQXPvpGa5+4pQcgpLeHf5QZ7+cgcezvbc2CnE7LJEahWFPBERkTruZE4Rd3+yiU1HzmCxwGMDIph4RagWtpAa69H+EeQWljI3/giPLtiGu5MD/dsGml2WSK2hLRRERETqsA2HTzPorV/YdOQMXi4OzBjblXuuDFPAkxrNYrHwzPVtGdY5BKvN4N55iazef8rsskRqDYU8ERGROqjEamPGqoOMmpFAek4RYY08WDypN/0iA8wuTeSC2NlZeHlYewa2DaTYamPCRxvZdOSM2WWJ1AoKeSIiInXMvhM59H9jFS99twerzeCGDsF8Oak3of4eZpcm8rc42Nvx5uiO9An3o6DEym2z17PzeJbZZYnUeBbD+G37U/lNdnY23t7eZGVl4eXlZXY5IiIiF2TpjjTe/nk/O49nA9DQ3YmH+7fi5u5NNTxTarX84lJiP1zPxiNnaODmSNyd0UQEeppdlsgld6E5RXfyREREarnswhIeW7CVuz/ZVB7wLm/lz9IHL2dMj2YKeFLruTk5MOu2bnQI8eZMfgljZiZwID3X7LJEaiyFPBERkVosM7+YsTPXsWBTCgA3dQ1h479i+Oj27vh7av87qTu8XBz56PYeRAZ5cSq3mJs/SODwqTyzyxKpkRTyREREaqmjp/MZNm0tW1Oy8HV3YuHd0bwyvIM2N5c6y9vNkU/u6EFEgCfpOUXc/EECR0/nm12WSI1zwXPyOnXqdMHDPTZv3lyposymOXkiIlLTbTpyhjs/2khGXjFB3i7Mvb07rQI0R0nqh5M5RYyaEc/Bk3mENHDls7uiaezjanZZItXuQnPKBW+GPmTIkKqoS0RERCpp3rpknvpyB1abQbvGXsyM7Uagt4vZZYlcMv6ezsyb0JOR0+M5nJHPzR8k8Pld0QR46XMgAlpd85x0J09ERGqi/OJSpi7Zw0fxRwAY1D6IV4ZH4e58wd/ZitQpxzMLGDkjnqOnCwj1dyfuzmjNRZU6TatrioiI1CG7U7Pp9/qq8oD3UEwr3rm5kwKe1GvBPq7Mu6Mnwd4uHDyZxy0z13E6r9jsskRMd1Ehz2q18tprr9G9e3cCAwPx9fWt8BAREZGqUWK1MfOXQ1z39mqOZRYQ7O3CrFu78kBMuLZGEAGa+Loxb0JPAryc2Xsih1tmriMzX0FP6reLCnnPPfccr7/+OiNHjiQrK4uHH36YoUOHYmdnx7PPPlvFJYqIiNQ/Wfkl3D8/kahnf+DFb3djtRlcGeHPl/dextWtA8wuT6RGae7nzrwJPfHzcGZXajZjP1xPVn6J2WWJmOai5uSFhoby1ltvMWjQIDw9PdmyZUv5cwkJCcybN686ar1kNCdPRETMtHLfSR7+bAsZvw4783V34qGYcG7pqY3NRf7MvhM5jJqRwOm8Yto39uaT8T3wdnM0uyyRKlOtc/LS0tJo3749AB4eHmRlZQFw3XXX8e23317MIUVERARYuiOVO+ZuICOvmJZ+7sTd2ZON/4xhbHRzBTyRv9AqwJN5E3rg6+7E9mNZjPkwQUM3pV66qJAXEhJCamoqUHZX74cffgBgw4YNODtrRSMREZGLMWt1EhM/3UyJ1WBQ+yCWPng5PVs2xM5O4U7kQrUO9GL+hJ40dHdix7FsxmiOntRDFxXybrzxRpYtWwbAfffdx1NPPUV4eDixsbHcfvvtVVqgiIhIXVditTF50Tae/2YXhgFjejTlzVEdcXLQItgiFyMi0JP5d/bEz8OJncezufmDdZzRqptSj1TJPnnx8fHEx8cTHh7O9ddfXxV1mUpz8kRE5FI5k1fMpHmbWXswA4sFHhsQwcQrQjU0U6QK7D+Rw+gP1nEqt4g2QV58ekfZUE6R2upCc4o2Qz8HhTwREbkUftiZxuRF28nIK8bNyZ53b+7MVa0bmV2WSJ1yID2HUTPKgl7rQE/mTeipoCe1VrVvhr5//35mzJjBiy++yPPPP1/hcaGmTZtGVFQUXl5eeHl5ER0dzZIlS87bfs6cOVgslgoPFxeXCm0Mw+Dpp58mKCgIV1dXYmJi2L9//8WepoiISJUrKLbyyOdbufPjTWULrPi7s+ieXgp4ItUgrJEncXf2xN/TmT1pOdz8QQIZuUVmlyVSrRwu5k0ffPABEydOxM/Pj8DAwApDSiwWC08//fQFHSckJISpU6cSHh6OYRjMnTuXwYMHk5iYSNu2bc/5Hi8vL/bu3Vvh7/u9V155hbfeeou5c+fSokULnnrqKQYMGMCuXbvOCoQiIiKX0u7UbD5OOMLyPemkZhUCMLJrE168sR2O9pp/J1Jdwhp5MH9CT0Z/kMCetBzGzFzHp3f0oKGHFgyUuumihms2a9aMe+65hyeeeKLKC/L19eXVV19l/PjxZ702Z84cHnzwQTIzM8/5XsMwCA4O5pFHHuHRRx8FICsri4CAAObMmcOoUaMuqAYN1xQRkapktRl8FH+YV7/fS36xFYAGbo68e3NneoX5mVydSP1x8GQuo2ckkJ5TRESAJ59O6IGfgp7UItU6XPPMmTOMGDHioos7F6vVSlxcHHl5eURHR5+3XW5uLs2aNaNJkyYMHjyYnTt3lr+WlJREWloaMTEx5c95e3vTo0cP4uPjz3vMoqIisrOzKzxqksISK7fP2UDCoQyzSxERkb/JZjN44n/beO7rXeQXW+kQ4s0HsV1Z9fhVCngil1iovwdxd/YkwMuZvSdyGDk9nhPZhWaXJVLlLirkjRgxonxvvMravn07Hh4eODs7c/fdd7N48WIiIyPP2TYiIoJZs2bx5Zdf8sknn2Cz2ejVqxcpKSlA2SbtAAEBARXeFxAQUP7auUyZMgVvb+/yR5MmTark3KrK9JWH+HlPOmNmrmP6yoNorRwRkdqhuNTGA59tYeGmFOztLDxzfSSf3x1Nv8gAPF0czS5PpF5q6e9B3J3RBHu7cPBkHjdNjyflTL7ZZYlUqYsarjllyhRef/11Bg0aRPv27XF0rPiL6v7777/gYxUXF5OcnExWVhYLFy5k5syZrFy58rxB7/dKSkpo06YNo0eP5oUXXmDt2rX07t2b48ePExQUVN7upptuwmKx8Nlnn53zOEVFRRQV/f8E3OzsbJo0aVJjhmvmF5fyz8U7WJx4DIDBHYN5dXgH7Z8kIlKDnc4r5vY5G9hyNBNHewuvjejA4I6NzS5LRH519HQ+N89M4OjpAoK9XZg3oSfN/dzNLkvkT1XrFgotWrQ4/wEtFg4dOvR3D1kuJiaG0NBQpk+ffkHtR4wYgYODA/Pnz+fQoUOEhoaSmJhIx44dy9tcccUVdOzYkTfffPOCjlkT5+QZhsEn65J57qudlNoMeoU25P2xXfDSN8EiIjXOnrRsJny0kaOnC/B2deTt0Z24vJW/2WWJyB+kZRVy88wEDp3Mo5GnM5/e0YPwAE+zyxI5r2qdk5eUlHTeR2UCHoDNZqtwV+3PWK1Wtm/fXn7XrkWLFgQGBrJs2bLyNtnZ2axbt+5P5/nVBhaLhbE9m/Hhrd1wd7Jn7cEMbno/nrQsjSMXEakpDMPgrWX7GfTWao6eLqCprxv/m9hLAU+khgr0duGzO6NpHehJek4RI2cksPN4ltlliVSaqeP9Jk+ezKpVqzh8+DDbt29n8uTJrFixgjFjxgAQGxvL5MmTy9s///zz/PDDDxw6dIjNmzdzyy23cOTIEe644w6gLAg9+OCDvPjii3z11Vds376d2NhYgoODGTJkiBmnWOWuaOXPZ3dFl+/1MmzaWg6dzDW7LBGRei85I59Bb63m9R/3YbUZ9IsM4MtJvQlr5GF2aSLyJ/w9nZk/oSftG3tzOq+Y0TMS2HI00+yyRCrlovbJe/jhh8/5/G+bk4eFhTF48GB8fX3/9Djp6enExsaSmpqKt7c3UVFRfP/99/Tr1w+A5ORk7Oz+P4eeOXOGCRMmkJaWRoMGDejSpQtr166tMH/v8ccfJy8vjzvvvJPMzEwuu+wyli5dWqf2yGvX2JtFE3sRO2s9SafyGP5+PHNu60ZUiI/ZpYmI1CtWm8GMVYdYuS+drUezKCix4upoz+MDI7it9/mnNohIzdLA3YlPJ/Tgttkb2HTkDLfMXMesW7vRvcWfX8uK1FQXNSfvqquuYvPmzVitViIiIgDYt28f9vb2tG7dmr1792KxWFi9evUFLaBS09TEOXnnciq3iFtnr2fHsWzcneyZEduV3lqOW0TkkjiZU8S98zazLul0+XORQV7MHNeVYB9XEysTkYuVV1TKHXM3En8oA1dHez6I7cpl4bq2kpqjWufkDR48mJiYGI4fP86mTZvYtGkTKSkp9OvXj9GjR3Ps2DEuv/xyHnrooYs+Aflrfh5lwwt6hTYkr9jKbbM38O22VLPLEhGp807lFnHzBwmsSzqNi6MdTwxszTf3Xca391+mgCdSi7k7OzD7tm5c0cqfghIrt8/dwM97TphdlsjfdlF38ho3bsyPP/541l26nTt30r9/f44dO8bmzZvp378/p06dqrJiL5XacifvN0WlVh76bAvfbU/DYoHnB7djbM9mZpclIlIn7T+Rw10fb+LQqTwCvVyYN6EHLf01706kLikqtXLfvER+2HUCBzsL/7lJW6BIzVCtd/KysrJIT08/6/mTJ0+SnZ0NgI+PD8XFxRdzePmbnB3seXt0Z27u0RTDgKe+2MGbP+3XpukiIlVs05EzjJgez6FTeb/uq6WAJ1IXOTvY8+6YzgzuGEypzeDBz7bwccIRs8sSuWAXPVzz9ttvZ/HixaSkpJCSksLixYsZP358+SqW69evp1WrVlVZq/wJezsL/x7Sjvv7hgPwxk/7eOarndhsCnoiIpVVYrXx35/2MXpGApn5JXRq6sNX912mgCdShzna2/HGTR0Z27NZ+Zfo7/ysL9Gldrio4Zq5ubk89NBDfPTRR5SWlgLg4ODAuHHjeOONN3B3d2fLli0AFTYlry1q23DNP5q79jDPfr0Tw4DrooJ4/aaOODmYuluGiEitdSQjj38u3sHqA2XTD/pHBvDGyI64O1/UAtUiUssYhsHrP+7j7Z8PAHDHZS3456A2WCwWkyuT+uhCc8pFhbzf5Obmlm9+3rJlSzw86sY3mrU95AF8tfU4j3y+hRKrQZ9wP96/pYsuSERELlB+cSnrkk6zIyWLWWuSOJNfgqO9halDoxjaubEu7kTqoQ9XJ/HCN7sAGNElhClD2+Ngry/R5dK6JCGvrqoLIQ9g1b6T3P3JJvKLrXQI8Wb2bd3xdXcyuywRkRrtQHoOt8/ZSPLp/PLnokK8eXlYFG2Cau/vBBGpvIWbUnjif9uw2gz6Rwbw1uhOuDjam12W1CNVHvKGDh3KnDlz8PLyYujQoX/adtGiRX+v2hqmroQ8gC1HM7lt9nrO5JfQ0t+dj8f3oLGW9xYROac1B04xfu4GCktseDo7cGXrRnRu6sPo7k11ISciAPywM4175ydSXGqjV2hDZsR2xUOjpeQSqfLVNb29vcuHp3h7e//pQ2qOjk18WHB3L4K9XTh0Mo9h761l/4kcs8sSEalxvt+Zxu1zygJe9xa+fHt/H94e3YnberdQwBORcv3bBjLntm64O9mz9mAGYz5I4HSeVpSXmuWihmsWFBRgs9lwd3cH4PDhw3zxxRe0adOGAQMGVHmRl1pdupP3m9SsAsZ+uJ4D6bl4uzoy69ZudGnWwOyyRERqhHnrkvnH4u0A9G3diPdu6Yyzg4KdiJzf1qOZ3PrraKmwRh58PL47Qd4aLSXVq1r3yRs8eDAff/wxAJmZmfTs2ZP//Oc/DBkyhGnTpl1cxVKtgrxdWXBXNJ2a+pBVUMKYmQks33v2XociIvVJRm4RT/5vW3nAG9OjKdNu6aKAJyJ/qUMTHxbcHU2glwsH0nMZ+t5a9mm0lNQQFxXyNm/eTJ8+fQBYuHAhAQEBHDlyhI8++oi33nqrSguUqtPA3YlP7+jBFa38KSyxMWHuRhYnpphdloiIKVbuO8mQ99YQt+EoAPdfHcaLQ9ppyxkRuWBhjTxZODGaUH93UrMKGT5tLRsOnza7LJGLC3n5+fl4enoC8MMPPzB06FDs7Ozo2bMnR44cqdICpWq5OTkwc1xXhnQMptRm8NBnW/lwdZLZZYmIXBJZBSVsOnKGl77bzbhZ6zl6uoAgbxfmT+jJw/0jtDWCiPxtIQ3cWHh3Lzo39SG7sJQxM9exdEeq2WVJPXdRIS8sLIwvvviCo0eP8v3339O/f38A0tPT68wctrrM0d6O12/qyG29mwPwwje7eGXpHrSbhojUZd9sO063F39i2LS1zFhVtsfryK5NWHxPb6JDG5pcnYjUZmWjpXoS0yaA4lIbEz/dzMfxh80uS+qxiwp5Tz/9NI8++ijNmzenR48eREdHA2V39Tp16lSlBUr1sLOz8PR1kTw2IAKA91Yc5Mn/bafUajO5MhGRqjdj1UHunZdIsdWGj5sjvUIb8sbIDkwd1p5AbxezyxOROsDVyZ73b+nM6O5NMQx46sudvPq9vkQXc1z0ZuhpaWmkpqbSoUMH7OzKsuL69evx8vKidevWVVrkpVYXV9f8M3Hry1aVsxloY08RqVNsNoNXf9jLtBUHARjdvQnP3aB5dyJSfQzD4K1lB3jjp30AjOgSwktD2+Nor587UnlVvhl6fVLfQh7A0h1p3B9XtrFnjxa+fDCuK14ujmaXJSJy0UqtNp79eiefJCQD8NiACCZdFWZyVSJSX/z+S/SrIvx5d0xn3Jy0abpUTrVuoSB1z8B2gcy9rTsezg6sSzrNyOkJpOcUml2WiMhF2Xcih+Hvx5cHvKlD2yvgicglNap7U2aM7YqLox3L955k9AfryMgtMrssqScU8qRcdGhD4u7siZ+HE7tTsxk+LZ4jGXlmlyUicsEMw2DW6iSGT1vLlqOZuDra896Yzozq3tTs0kSkHoqJDGDehJ40cHNk69FMhk1bS9IpXVtJ9dNwzXOoj8M1f+/wqTzGzlrH0dMF+Hk48+G4rnRo4mN2WSIi55SZX8zKfSc5lVtM/MEMftp9AoB2jb2YGdtNC6uIiOkOnsxl3Kz1pJwpoIGbIzNiu9Ktua/ZZUktpDl5lVDfQx5AenYh42ZvYHdqNi6Odvx3ZCcGtgs0uywRkQoOncxlzMx1pGb9//Byi6Vs/t3tvVtoESkRqTFO5hRxx9wNbE3Jwsnejtdu6sANHYLNLktqGYW8SlDIK5NbVMq98zazYu9JLBb4xzVtuKNPC20WLCI1wsGTuYyakcDJnCI8nR24snUjgn1cuK59MO1DvM0uT0TkLAXFVh78LJHvd5aNOHhsQAT3XBmqayu5YAp5laCQ9//+uDrdmB5Nee6GtjhoGWARMdHetBxiZ63jRHYRrQM9+eSOHvh5OJtdlojIX7LaDKZ8t5uZq5MAGNm1CS/e2E5bLMgF0eqaUiUc7O14YXA7/jWoDRYLfLoumfFzN5JTWGJ2aSJSDxmGwReJx7j+7dWcyC6iVYAHnyrgiUgtYm9n4V/XRfLC4LbYWeCzjUe5dfZ6sgp0bSVVRyFP/pLFYuGOPi15/5YuuDras3LfSUa8H8/xzAKzSxOReiS3qJR/fbGDBz/bQrHVxuWt/Jk3oScNFfBEpBYaG92cmeO64uZkz5oDGQyftpaUM/lmlyV1hEKeXLABbQP57K6e+Hs6sycthyHvrmF7SpbZZYlIHWcYBpuTz3DFK8v5dF3Z0PHbejdnZmxX3cETkVrt6tYBfH5XNAFezuxPz2XIu2vZlpJpdllSB2hO3jloTt6fO5ZZwO2zN7D3RA6ujva8NboT/SIDzC5LROqY77an8tr3ezmWWUBRqQ2AIG8XnrymNYM7Nja5OhGRqpOaVcBtszewJy0HF0c7/jOiI4OigswuS2ogLbxSCQp5fy2nsIR7Pt3ML/tPYbHA5GtaM6FPS60OJSKVZhgG7/x8gP/8uK/8OXs7C/3aBPDy8Ci8XR1NrE5EpHrkFJZw3/xEVuw9CcBDMa24v2+Yrq2kglqx8Mq0adOIiorCy8sLLy8voqOjWbJkyXnbf/DBB/Tp04cGDRrQoEEDYmJiWL9+fYU2t956KxaLpcJj4MCB1X0q9Y6niyOzb+3GmB5NMQx46bs9PLpgG0WlVrNLE5FazGYzeHnp3vKANygqiFWPXcXO5wbw/tguCngiUmd5ujjy4bhujL+sBQBv/LSPe+cnUlCsayv5+0wNeSEhIUydOpVNmzaxceNGrr76agYPHszOnTvP2X7FihWMHj2a5cuXEx8fT5MmTejfvz/Hjh2r0G7gwIGkpqaWP+bPn38pTqfecbC348Uh7XjuhrbY21n43+YURs9IID2n8K/fLCLyB/nFpYyfu4H3Vx4E4KnrInn35s40beimTc1FpF6wt7Pw1HWRTB3aHgc7C99uS2XkjHjSsnRtJX9PjRuu6evry6uvvsr48eP/sq3VaqVBgwa88847xMbGAmV38jIzM/niiy8uugYN1/z7ftl/kkmfbia7sJRgbxdmxHalXWNtRiwiF+bgyVxiP1zPscwCXBzteP6GdtzUrYnZZYmImCbhUAYTP9nEmfwSAryc+SC2K1EhPmaXJSarFcM1f89qtRIXF0deXh7R0dEX9J78/HxKSkrw9fWt8PyKFSto1KgRERERTJw4kYyMjD89TlFREdnZ2RUe8vf0Cffni0m9aenvzvGsQka8H8+S7almlyUiNZxhGLy9bD8D3ljFscwC/D2dmTWumwKeiNR7PVs25MtJlxHeyIMT2UWMeD+er7ceN7ssqSVMv5O3fft2oqOjKSwsxMPDg3nz5nHttdde0Hvvuecevv/+e3bu3ImLiwsAcXFxuLm50aJFCw4ePMg//vEPPDw8iI+Px97+3MN9nn32WZ577rmzntedvL8vq6Bs0vCqfZo0LCLntuNYFu/8fICTuUWkZRVy7Nc9N7s1b8B7Y7rg76ltEUREfpNTWML98xNZ/uuCLPf3DefBvuHY2enaqj6qNatrFhcXk5ycTFZWFgsXLmTmzJmsXLmSyMjIP33f1KlTeeWVV1ixYgVRUVHnbXfo0CFCQ0P56aef6Nu37znbFBUVUVRUVP7n7OxsmjRpopB3kUqtNl76bg+z1iQBZQsnvDa8A65OmlMjUt99vuEoT3+1g8ISW/lzTg52/GtQG2Kjm5tXmIhIDWa1GUxdspsPfim7trq2fSCvDu+Au7ODyZXJpVZrQt4fxcTEEBoayvTp08/b5rXXXuPFF1/kp59+omvXrn95TH9/f1588UXuuuuuC6pBc/KqxmcbkvnXFzsosRq0a+zF9LFdaezjanZZImKSxYkpPPz5VgwDerb05ZaezfB1cyKqiQ8eulAREflLn288yj8Xb6fEahAR4MmM2C40a+hudllyCdW6OXm/sdlsFe6q/dErr7zCCy+8wNKlSy8o4KWkpJCRkUFQkDaUvNRGdmvKp3f0xNfdiR3Hsrn+7dWsPXjK7LJE5BIzDIPPNiTz6IJtGAaMi27G/Ak9uS4qmF5hfgp4IiIX6KauTYi7syf+ns7sPZHD9W+vZsXedLPLkhrI1JA3efJkVq1axeHDh9m+fTuTJ09mxYoVjBkzBoDY2FgmT55c3v7ll1/mqaeeYtasWTRv3py0tDTS0tLIzc0FIDc3l8cee4yEhAQOHz7MsmXLGDx4MGFhYQwYMMCUc6zvurfw5at7e9M22IvTecWM/XA9M385RA27gSwi1eRMXjGPLNjKE//bjtVmMKxzCM9c31bzdEVELlKXZr58c99ldGrqQ3ZhKbfN2cB7Kw7o2koqMDXkpaenExsbS0REBH379mXDhg18//339OvXD4Dk5GRSU/9/hcZp06ZRXFzM8OHDCQoKKn+89tprANjb27Nt2zZuuOEGWrVqxfjx4+nSpQu//PILzs6ayG+WkAZu/G9iL4Z2aozVZvDit7t58LMt2txTpA4qsdoosdooKrUSfzCDq/+zgkWbj2GxwAN9w3l5WHstFiAiUkkBXi7E3dmTUd2aYBjwytK93DsvkbyiUrNLkxqixs3Jqwk0J696GIbB3LWHefHb3ZTaDNoEeTFjbBea+LqZXZqIVFKp1cZTX+4kbkMyf/yt0ryhG09e05qB7TRsXkSkqs1bl8wzX+3QPL16otYuvFITKORVr3WHMpg0bzOncovxcXPk7dGd6BPub3ZZInKRbDaDRxZsZXHisQrPuzvZc0WEPy8MbkdDD42mEBGpLpuOnObuTzZzMqcILxcH3hrdiSsjGpldllQDhbxKUMirfqlZBdz9yWa2Hs3EzgKPD2zNXZe31DwdkVrGMAye+3oXc9YexsHOwivDo+jbOgAAL1cHfaZFRC6RE9mF3P3JJhKTM7FY4NH+EUy8IlRD5OsYhbxKUMi7NIpKrTzz5U7iNhwFyvZ8eXlYFJ4ujiZXJiIXIqewhDvmbmRd0mkA3hjZgRs7hZhclYhI/VVUauXZr3Yyf33ZtVXf1o14/aaOeLvp2qquUMirBIW8S+v3Y8lb+rnz3i2daR2o/+8iNdmXW47x7Fc7OZNfgrODHc9c35abezQ1uywREQHi1ifz9Fc7KS61EdLAlWljutA+xNvssqQKKORVgkLepbc5+Qz3frqZ41mFuDja8eKQ9gzvojsCIjXBjmNZvP3zfgpKbDjYWcguKGHjkTMABHq5MCO2C1EhPuYWKSIiFew4lsU9n24m+XQ+TvZ2PHNDJDd3b6ph9LWcQl4lKOSZ43ReMQ99toWV+04CMLJrE54b3BYXR3uTKxOpv7alZDJm5jpyCs9elnt4lxCmDG2Po72pu/GIiMh5ZBWU8MjnW/lp9wkAhnZqzIs3tsPNycHkyuRiKeRVgkKeeWw2g/dWHOD1H/dhM6BNkBfvjelMCz8tBSxyqR0+lcewaWvJyCsm1N+dO/q0xN5iwdnRjrbBXoQ18jS7RBER+QuGYTB91SFe/X4vVlvZNgvv3dKZUH8Ps0uTi6CQVwkKeeZbc+AUD8Qlciq3GE9nB14dEaU9tkQuoTUHTvHogq2kZhXSrrEXcXdG4+Gsb35FRGqrhEMZ3Dc/kZM5RXg4O/DysCgGRenaqrZRyKsEhbya4UR2IffNS2T94bKV+27v3YInr2mNk4OGholUl8ISK+8tP8D7Kw9RbC2bsL/onl408nQxuzQREamk9OxC7p2fyPpfV0UeF92Myde20dSYWkQhrxIU8mqOUquNV3/Yy/SVhwDo2MSHt0d3oomvm8mVidQNhmGw9mAGaVmFFJZaWbgphcTkTABi2gTwnxEdtPS2iEgd8sdrq8ggL965uRMtNXyzVlDIqwSFvJrnx10neOTzLWQXluLp7MCUYe25LirY7LJEarX0nEKmfLeHxYnHKjzv4ezAM9dHcmOnxjhoURURkTpp+Z50HlmwldN5xbg52fPikHYM7ayVzWs6hbxKUMirmVLO5PNA3BY2/bp0++juTXj6ura4OmmIgcjflXQqj5HT40nPKQIgKsQbPw9nwgM8GNO9GU0b6m65iEhddyK7kAfiEkk4VDZ8c1jnEJ4f3BZ3zcGusRTyKkEhr+Yqtdr470/7eXfFAQwDwht58PbNnbR5usjfsOnIGe6fn8ixzAJa+Lnz1HVtuLp1gNlliYiICaw2g3eXH+C/P5WtbN7Sz523b+5E22Btnl4TKeRVgkJezbf2wCke/GwL6TlFODvY8a/rIrmlhzb4FPkzVptB3IZknv1qJyVWg2YN3Vh4dy/8PZ3NLk1EREy27lAGD8RtIS27ECcHO/41qA1jezbTtVUNo5BXCQp5tUNGbhGPLNjKir1lm6cPbBvIy8OitEiEyK/Ssgo5kV2Ig72F/GIr7684yLI96QBcGeHPK8OjtGqmiIiUO51XzGMLtpb/rhjQNoCpQ6No4O5kcmXyG4W8SlDIqz1sNoNZa5J4eekeSqwGjX1cef2mDvRo2dDs0kRMtXxPOnd/somiUluF550c7Lj/6jAmXN4SZwfNZxURkYoMw2DWmsNMXbKbEqtBI09nXhvRgctb+ZtdmqCQVykKebXPtpRM7pufyJGMfCwWuPPyljzcr5UuYqVemrU6iee/2VX+50aezjjYWYgM9uaBvuG0D9E8CxER+XM7jmVxf1wih07mAXBb7+Y8MbC19tQzmUJeJSjk1U65RaU8//VOPt+YApTt+/LfUR1pFeBpcmUil86n647wry92YBgwtFNjnh/SDg+tkiYiIhehoNjKS9/t5uOEIwC0CvDgzVGdaBOk62OzKORVgkJe7bZ0RxqTF23jTH4JTg52PDmwNbf2ao6dnSYOS911PLOAd5cf4NN1yQCM6taEKUPba8K8iIhU2vI96Ty2cCuncotxsrfj8YER3N67ha6tTKCQVwkKebVfenYhj/9vW/miLH3C/Xh1eAcCvbXIhNR+Wfkl/HLgJI72dng4O5CRV8wL3+zi5K973t1zZSiP9I/AXr98RUSkipzKLeLJ/23jp91li7L0Cm3If27qQJC3q8mV1S8KeZWgkFc3GIbBJwlH+Pd3uyksseHt6shLN7ZnUFSQ2aWJXLQD6bmM/XAdqVmFZ70W0sCVp66LpH9kgO7giYhIlTMMg/nrj/LCN7soKLHi5eLAC0PacUOHYP3euUQU8ipBIa9uOZCey0OfbWH7sSwABncM5tnr22o5YKl19qRlM3J6AlkFJQB0aupDfpEVFyd7OoZ4M+mqMBp56W61iIhUr0Mny66ttqaUXVsNbBvIize2w89D+65WN4W8SlDIq3tKrDbeWrafd5cfwGaAn4czL93Yjv5tA80uTeSCbDpymomfbCY9p4gOId68O6YzIQ3czC5LRETqqRKrjWkrDvLWsv2U2gx83Z14YXA7jZiqZgp5laCQV3clJp/hsYXbOJCeC+iuntR8OYUlzFh1iFmrk8grttLS353FE3vj7eZodmkiIiLsPJ7Fowu2sTs1G4DrooJ4fnA7fHVtVS0U8ipBIa9uKyyx8t+f9jNj1UHd1ZMaxTAMvtp6nJM5RXg4O2AzYP765PKhxr1CG/L+2C54uSjgiYhIzVFcauOdn/fz7oqDWG0Gfh5OvDikPQPb6dqqqinkVYJCXv2gu3pSk+QWlTLlu93lWyD8npuTPU9dF8mwziE4OdiZUJ2IiMhf256SxSMLtrDvhK6tqotCXiUo5NUf57qr9+8b2zFAd/XkEioothI7ax0bDp8BILyRB0193TCAFn7u3NyjKaH+HuYWKSIicgGKSq28+dN+3l/527WVE09f35bro4K0AmcVUMirBIW8+uePd/WuaRfIsze0JUArFUo1S88p5N5PE1l/+DSeLg68PCyKa9oF6hehiIjUaonJZ3h84Tb2/3ptdXXrRrwwpB2NfbSvXmVcaE4xddzPtGnTiIqKwsvLCy8vL6Kjo1myZMmfvmfBggW0bt0aFxcX2rdvz3fffVfhdcMwePrppwkKCsLV1ZWYmBj2799fnachdUCnpg345r7LmHhlKPZ2FpbsSCPmPyv5JOEINpu+B5GqV1hi5f2VBxk2bS3rD5/G2cGOD2K7cm17fdMpIiK1X6emDfj2/j48FNMKJ3s7ft6TTv/XVzJnTRJWXVtVO1Pv5H399dfY29sTHh6OYRjMnTuXV199lcTERNq2bXtW+7Vr13L55ZczZcoUrrvuOubNm8fLL7/M5s2badeuHQAvv/wyU6ZMYe7cubRo0YKnnnqK7du3s2vXLlxcLuyujO7k1W+7jmczedG28r1fujZrwJSh7QkP8DS5MqnNsgpK+G57KiVWGw52dny55Rjrkk4DEOztwtzbu+vfmIiI1EkH0nN48n/b2XikbFpCp6Y+vDwsilb6vfe31drhmr6+vrz66quMHz/+rNdGjhxJXl4e33zzTflzPXv2pGPHjrz//vsYhkFwcDCPPPIIjz76KABZWVkEBAQwZ84cRo0adUE1KOSJ1WbwUfxhXvt+L3nFVhztLUy8Mox7rgzFxdHe7PKkFrHaDH7Ymcaby/azJy2nwmsOdhYe6R/B6O5N8HHTpHQREam7bDaDT9cn8/KSPeQWlZZdW10RyqSrw3B20LXVhaoVwzV/z2q1EhcXR15eHtHR0edsEx8fT0xMTIXnBgwYQHx8PABJSUmkpaVVaOPt7U2PHj3K25xLUVER2dnZFR5Sv9nbWbitdwt+fPgKYto0osRq8Nay/Vz71i8kHMowuzypJQzD4NEFW5n46ebygBfTJoC+rRsxLroZ/5vYi4lXhirgiYhInWdnZ2Fsz2b8+PDl9IsMKLu2+vkA1/z3F37Zf9Ls8uocB7ML2L59O9HR0RQWFuLh4cHixYuJjIw8Z9u0tDQCAgIqPBcQEEBaWlr56789d7425zJlyhSee+65ypyG1FHBPq58ENuVJTvSeOarnRw6mceoGQnc2Kkxk69tTSNPLcwi55aRW8S/v93N4sRjAIzoEsIdfVoSEaihKSIiUn8FebsyY2wXlu5I4+mvdnLoVB5jP1zPoPZB/Ou6NgR5a2GWqmB6yIuIiGDLli1kZWWxcOFCxo0bx8qVK88b9KrD5MmTefjhh8v/nJ2dTZMmTS7Z3y81m8Vi4dr2QfQO8+PlpXuYvz6ZxYnH+GnXCR7s14px0c1wsK8xN8XlEikqtfLQZ1vYejSLlv7uhPp70MLPHW9XR0qsNqavOlS+WuuLQ9pxS89mJlcsIiJSM1gsFq5pH8Rl4X688eN+5sYf5tvtqSzfm84DfcO5/bIWOOraqlJMD3lOTk6EhYUB0KVLFzZs2MCbb77J9OnTz2obGBjIiRMnKjx34sQJAgMDy1//7bmgoKAKbTp27HjeGpydnXF2dq7sqUgd5+3qyEs3tmdk1yY8/eUOtqZk8cI3u1iw8SjP3dCWHi0bml2iXCKFJVbGz93AmgNlQ3ePZRbwy/5TZ7Vr4ObIm6M6cXkr/0tdooiISI3n6eLI09dHMrxLCE99uYNNR84wZckeFm5K4YUh7eipa6uLVuMiss1mo6io6JyvRUdHs2zZsgrP/fjjj+Vz+Fq0aEFgYGCFNtnZ2axbt+688/xE/q4OTXxYfE9vpgxtj4+bI3vSchg5I4GHPttCenah2eVJJS3anMIdczfy8tI9fLPtOIdP5VXYRqPUauO++YmsOZCBu5M9TwxszcvD2nPn5S0Z0DaAy8L8uCzMj1t7NeeLSb0V8ERERP5CZLAXC+6K5pXhUfi6O7E/PZdRMxJ4MC6RE7q2uiimrq45efJkrrnmGpo2bUpOTk75lgjff/89/fr1IzY2lsaNGzNlyhSgbAuFK664gqlTpzJo0CDi4uJ46aWXztpCYerUqRW2UNi2bZu2UJBqcSavmFd/2Mv89ckYBng4O/BA33DG9WqOk0ON+w5F/sK321K5b/5m/rh9j7ODHf6ezni7OnIiu4hTuUU4Odjx0e3d9S2jiIhIFcrML+bV7/cy79drKzcne+65MpQ7+rTUCufUki0Uxo8fz7Jly0hNTcXb25uoqCieeOIJ+vXrB8CVV15J8+bNmTNnTvl7FixYwL/+9S8OHz5MeHg4r7zyCtdee23564Zh8MwzzzBjxgwyMzO57LLLeO+992jVqtUF16WQJ3/XtpRMnvpiR/nees0aujH5mjYMaBugja0vwum8YvKKSglp4Hre/3+GYXAiu4gAL+cq+X+86cgZbv4ggaJSG80buhEd6seu41nsTsuhuNRWoa2nswOvj+xIv8iA8xxNREREKmPr0Uye/XonicmZADT2cWXyta0Z1D6oXl9b1YqQV1Mp5MnFsNkMFm5K4dUf9nIyp2zIcY8Wvjx1XSTtGnubXF3tsfN4FqNnJJBdWIqvuxMdQrzp2KQBbYI8adzAFX9PZ4pKbLz+4z4WJx7Dz8OJ6FA/eoc2pHeYH0183f7W33cqt4g3f9rPos0p5BVbuSrCn5njumFvV/YLpMRqIzWzkJO5hWQXlOLp4kDbYG9cnfRtooiISHUyDIOvth7n5SV7OJ5VNmyza7MGPH19JFEhPuYWZxKFvEpQyJPKyC0q5f0VB/ngl0MUldqwWGBopxAeHxhBgJe2XPgzRzLyGDYtnlO5556XeyGa+LpyWZgfvUL96Nbc97x3+hKTz7Bq3ym+2nqMgyfzAIhu2ZAPb+2Km5Ppa1KJiIjIrwqKrXzwyyGmrThIQYkVgGGd6+e1lUJeJSjkSVU4llnAq0v38MWW4wC4Otpz1xUtmdCnJe7OChG/Z7UZbD+WxX3zN3P0dAFtgrz46PbuHMssYOvRTLYczeTgyVyOnSkgI68YezsLrQI8ue/qMBq6O7HmwCnWHMxgy9FMrH+YUOfn4Uz7xl5EBnsR6O2Kt6sjR07l8cZP+8rn3vm4OfLikHYMbBuo7TBERERqqLSsQl75fg+LNpftQevqaM/tlzXnzstD8XZ1NLm6S0MhrxIU8qQqJSaf4cVvd7PpyBkA/DycmHRVGDf3aIqzQ/0d8lditeFgZ+F0XjHj525ky9FMoGw+44K7o8+70bzNZmBnd+6x+LlFpaxPymDNgQzWHsxgb1r2WYuo/F7rQE+u7xDMsM4hBHrXr28CRUREaqstRzN58ZtdbPz12srHzZF7rgwlNrp5nV+cRSGvEhTypKoZhsG321N57fu9HM7IB8omED8QE87QTo3r1d2j4lIbD3++hW+2peJgZ6H0dymsX2QAz97QlsY+rlXydxUUW9mVms2OY1nsO5FDek4R2QUleLo40ifcj5t7NNVmqyIiIrWQYRj8tDudV5buYX96LgBB3i48FNOKoZ3r7rWVQl4lKORJdSmx2liwMYU3l+3jRHbZvLNQf3ce7R/BwHaBdXq1qPTsQu75dDOJ5xhS2dLPnVdHdKBLswYmVSciIiK1kdVmsGhzCm/8uK98cZawRh480q8VA9oGnnf0T22lkFcJCnlS3QpLrHwUf5j3VhwkM78EgPaNvXmgbzh92zSqc2Gv1GpjzMx1rEs6DZQNWX18YGv6hPvh4mBPA3cnkysUERGR2qywxMonCUd4Z/mB8mur1oGePNA3vE6FPYW8SlDIk0slu7CEmasOMXN1EvnFZatFRQZ5cd/VYbXmB1KJ1canCUdwdrSnQ4gPrQI8KgyR2HU8m/vjEjmQnouHswPzJ/QkMtirfIsCERERkary27XV7DWHySkqBSAiwJP7+4ZzTbvacW31ZxTyKkEhTy61U7lFfPDLIT6OP1Ie9loFeHDv1eEMah9UYwNRRm4RTy7azo+7TpQ/5+poT2gjdwK9XLBYLCQczCCnqBQXRzveGd2ZGG0gLiIiItUsK7+ED9ckMXt1UnnYaxXgwX1Xh3NtDb62+isKeZWgkCdmOZNXzKw1Scz53bdPLf3dufvyUG7oGFyjVozKyi8hdtY6tqZkAWVbFRSWWMn9te7fa9fYiw/Hdat3e9mIiIiIubLyS5i1JolZa5LIKSy7Rmnq68YdfVowoksTXJ1qzrXVhVDIqwSFPDFbVkEJc9ce5sPVSWQVlI0r9/NwIja6OWN6NKWhh7NptVltBqv2neTlpXvYk5aDu5M9s27tRo+WDbHaDA6dzCX5dD5p2WWTnxv7uNI7zE+rWIqIiIhpsgpKmLPmMHPWJnHm1zl7DdwcGRvdnHHRzUy9tvo7FPIqQSFPaorcolLmrTvCnDWHy1eMcnawY2jnEMZf1pywRp7V8vduTj7Dv7/dTaC3Cx1CvGkd6EUjL2dsNpi6dA+r9p0EwNvVkY9u706HJj7VUoeIiIhIVcovLmXhphRm/pJE8umyba2cHewY1iWE2OhmtA6s2df+CnmVoJAnNU2J1caSHWnM/OUQ234dHgnQK7Qho7s3pX/bgCrZWL3EauO5r3eyYGMKRaW2P217Q4dgJl0VRkRg9QRNERERkepitRks3ZHGjFUHy6eeAHRt1oBbejbjmvaBVXJtVdUU8ipBIU9qKsMw2HjkDDN/OcQPu07w26fX192J4V1CGNWtCS39PS76+C99t5sZqw4BEOjlwoiuIRxIz+VAei6ncouwWCyE+rvzYEwreof5VcUpiYiIiJjGMAzWJZ3mo/jD/LDzBKW/7uXr6+7EiK4hjOxauWurqqaQVwkKeVIbpJzJ5/MNR/ls49HyjdUBOjTx4fqoIAZFBRHk7XpBx9p6NJPXf9zHyl+HYb42ogNDOzWu9csMi4iIiFyo9OxC4jYcZf76ZFJ/nSYD8M7NnbguKtjEyv6fQl4lKORJbVJqtbF870ni1iezfG86tt99ors1b0DfNgFcFuZHZJDXWaFtT1o2n204yucbjpL369YN91wZyuMDW1/KUxARERGpMUqtNn7ek07chqPEH8wgfvLV+Lg5mV0WoJBXKQp5Ulul5xSyZHsa32w7zobDZyq81tDdiagQbyICvQj0cqbYauOtZQfKtzxoHejJayM60K6xtxmli4iIiNQ42YUleLk4ml1GOYW8SlDIk7ogNauA73eksfrAKeIPZpTfqfuj5g3duPfqcK6LCqpR+/CJiIiISEUKeZWgkCd1TYnVxraUTHan5rA3LYfT+cVgQGSwF7f0bIa3a835hkpEREREzu1Cc4rDJaxJREziaG9Hl2a+dGnma3YpIiIiIlLN7MwuQERERERERKqOQp6IiIiIiEgdopAnIiIiIiJShyjkiYiIiIiI1CEKeSIiIiIiInWIQp6IiIiIiEgdopAnIiIiIiJSh2ifvHP4bX/47OxskysREREREREp81s++S2vnI9C3jnk5OQA0KRJE5MrERERERERqSgnJwdvb+/zvm4x/ioG1kM2m43jx4/j6emJxWIxuxygLLU3adKEo0eP4uXlZXY5Uk3Uz/WH+rr+UF/XH+rr+kH9XH/UxL42DIOcnByCg4Oxszv/zDvdyTsHOzs7QkJCzC7jnLy8vGrMPzKpPurn+kN9XX+or+sP9XX9oH6uP2paX//ZHbzfaOEVERERERGROkQhT0REREREpA5RyKslnJ2deeaZZ3B2dja7FKlG6uf6Q31df6iv6w/1df2gfq4/anNfa+EVERERERGROkR38kREREREROoQhTwREREREZE6RCFPRERERESkDlHIqwXeffddmjdvjouLCz169GD9+vVmlyR/w7PPPovFYqnwaN26dfnrhYWFTJo0iYYNG+Lh4cGwYcM4ceJEhWMkJyczaNAg3NzcaNSoEY899hilpaWX+lTkD1atWsX1119PcHAwFouFL774osLrhmHw9NNPExQUhKurKzExMezfv79Cm9OnTzNmzBi8vLzw8fFh/Pjx5ObmVmizbds2+vTpg4uLC02aNOGVV16p7lOTP/irvr711lvP+pwPHDiwQhv1dc03ZcoUunXrhqenJ40aNWLIkCHs3bu3Qpuq+pm9YsUKOnfujLOzM2FhYcyZM6e6T09+50L6+sorrzzrc3333XdXaKO+rvmmTZtGVFRU+V530dHRLFmypPz1OvuZNqRGi4uLM5ycnIxZs2YZO3fuNCZMmGD4+PgYJ06cMLs0uUDPPPOM0bZtWyM1NbX8cfLkyfLX7777bqNJkybGsmXLjI0bNxo9e/Y0evXqVf56aWmp0a5dOyMmJsZITEw0vvvuO8PPz8+YPHmyGacjv/Pdd98Z//znP41FixYZgLF48eIKr0+dOtXw9vY2vvjiC2Pr1q3GDTfcYLRo0cIoKCgobzNw4ECjQ4cORkJCgvHLL78YYWFhxujRo8tfz8rKMgICAowxY8YYO3bsMObPn2+4uroa06dPv1SnKcZf9/W4ceOMgQMHVvicnz59ukIb9XXNN2DAAGP27NnGjh07jC1bthjXXnut0bRpUyM3N7e8TVX8zD506JDh5uZmPPzww8auXbuMt99+27C3tzeWLl16Sc+3PruQvr7iiiuMCRMmVPhcZ2Vllb+uvq4dvvrqK+Pbb7819u3bZ+zdu9f4xz/+YTg6Oho7duwwDKPufqYV8mq47t27G5MmTSr/s9VqNYKDg40pU6aYWJX8Hc8884zRoUOHc76WmZlpODo6GgsWLCh/bvfu3QZgxMfHG4ZRdnFpZ2dnpKWllbeZNm2a4eXlZRQVFVVr7XLh/njhb7PZjMDAQOPVV18tfy4zM9NwdnY25s+fbxiGYezatcsAjA0bNpS3WbJkiWGxWIxjx44ZhmEY7733ntGgQYMKff3EE08YERER1XxGcj7nC3mDBw8+73vU17VTenq6ARgrV640DKPqfmY//vjjRtu2bSv8XSNHjjQGDBhQ3ack5/HHvjaMspD3wAMPnPc96uvaq0GDBsbMmTPr9GdawzVrsOLiYjZt2kRMTEz5c3Z2dsTExBAfH29iZfJ37d+/n+DgYFq2bMmYMWNITk4GYNOmTZSUlFTo49atW9O0adPyPo6Pj6d9+/YEBASUtxkwYADZ2dns3Lnz0p6IXLCkpCTS0tIq9K23tzc9evSo0Lc+Pj507dq1vE1MTAx2dnasW7euvM3ll1+Ok5NTeZsBAwawd+9ezpw5c4nORi7EihUraNSoEREREUycOJGMjIzy19TXtVNWVhYAvr6+QNX9zI6Pj69wjN/a6He7ef7Y17/59NNP8fPzo127dkyePJn8/Pzy19TXtY/VaiUuLo68vDyio6Pr9GfawbS/Wf7SqVOnsFqtFf5RAQQEBLBnzx6TqpK/q0ePHsyZM4eIiAhSU1N57rnn6NOnDzt27CAtLQ0nJyd8fHwqvCcgIIC0tDQA0tLSzvlv4LfXpGb6rW/O1Xe/79tGjRpVeN3BwQFfX98KbVq0aHHWMX57rUGDBtVSv/w9AwcOZOjQobRo0YKDBw/yj3/8g2uuuYb4+Hjs7e3V17WQzWbjwQcfpHfv3rRr1w6gyn5mn69NdnY2BQUFuLq6VscpyXmcq68Bbr75Zpo1a0ZwcDDbtm3jiSeeYO/evSxatAhQX9cm27dvJzo6msLCQjw8PFi8eDGRkZFs2bKlzn6mFfJEqtk111xT/t9RUVH06NGDZs2a8fnnn+uHu0gdMWrUqPL/bt++PVFRUYSGhrJixQr69u1rYmVysSZNmsSOHTtYvXq12aVINTtfX995553l/92+fXuCgoLo27cvBw8eJDQ09FKXKZUQERHBli1byMrKYuHChYwbN46VK1eaXVa10nDNGszPzw97e/uzVvg5ceIEgYGBJlUlleXj40OrVq04cOAAgYGBFBcXk5mZWaHN7/s4MDDwnP8GfntNaqbf+ubPPr+BgYGkp6dXeL20tJTTp0+r/2u5li1b4ufnx4EDBwD1dW1z77338s0337B8+XJCQkLKn6+qn9nna+Pl5aUv/y6x8/X1ufTo0QOgwudafV07ODk5ERYWRpcuXZgyZQodOnTgzTffrNOfaYW8GszJyYkuXbqwbNmy8udsNhvLli0jOjraxMqkMnJzczl48CBBQUF06dIFR0fHCn28d+9ekpOTy/s4Ojqa7du3V7hA/PHHH/Hy8iIyMvKS1y8XpkWLFgQGBlbo2+zsbNatW1ehbzMzM9m0aVN5m59//hmbzVZ+MREdHc2qVasoKSkpb/Pjjz8SERGh4Xs1WEpKChkZGQQFBQHq69rCMAzuvfdeFi9ezM8//3zW8Nmq+pkdHR1d4Ri/tdHv9kvnr/r6XLZs2QJQ4XOtvq6dbDYbRUVFdfszbdqSL3JB4uLiDGdnZ2POnDnGrl27jDvvvNPw8fGpsMKP1GyPPPKIsWLFCiMpKclYs2aNERMTY/j5+Rnp6emGYZQt3du0aVPj559/NjZu3GhER0cb0dHR5e//bene/v37G1u2bDGWLl1q+Pv7awuFGiAnJ8dITEw0EhMTDcB4/fXXjcTEROPIkSOGYZRtoeDj42N8+eWXxrZt24zBgwefcwuFTp06GevWrTNWr15thIeHV1hWPzMz0wgICDDGjh1r7Nixw4iLizPc3Ny0rP4l9md9nZOTYzz66KNGfHy8kZSUZPz0009G586djfDwcKOwsLD8GOrrmm/ixImGt7e3sWLFigrL5ufn55e3qYqf2b8tt/7YY48Zu3fvNt59913Tl1uvb/6qrw8cOGA8//zzxsaNG42kpCTjyy+/NFq2bGlcfvnl5cdQX9cOTz75pLFy5UojKSnJ2LZtm/Hkk08aFovF+OGHHwzDqLufaYW8WuDtt982mjZtajg5ORndu3c3EhISzC5J/oaRI0caQUFBhpOTk9G4cWNj5MiRxoEDB8pfLygoMO655x6jQYMGhpubm3HjjTcaqampFY5x+PBh45prrjFcXV0NPz8/45FHHjFKSkou9anIHyxfvtwAznqMGzfOMIyybRSeeuopIyAgwHB2djb69u1r7N27t8IxMjIyjNGjRxseHh6Gl5eXcdtttxk5OTkV2mzdutW47LLLDGdnZ6Nx48bG1KlTL9Upyq/+rK/z8/ON/v37G/7+/oajo6PRrFkzY8KECWd9Gae+rvnO1ceAMXv27PI2VfUze/ny5UbHjh0NJycno2XLlhX+Dql+f9XXycnJxuWXX274+voazs7ORlhYmPHYY49V2CfPMNTXtcHtt99uNGvWzHBycjL8/f2Nvn37lgc8w6i7n2mLYRjGpbtvKCIiIiIiItVJc/JERERERETqEIU8ERERERGROkQhT0REREREpA5RyBMREREREalDFPJERERERETqEIU8ERERERGROkQhT0REREREpA5RyBMREREREalDFPJEREQu0ooVK7BYLGRmZppdioiISDmLYRiG2UWIiIjUBldeeSUdO3bkv//9LwDFxcWcPn2agIAALBaLucWJiIj8ysHsAkRERGorJycnAgMDzS5DRESkAg3XFBERuQC33norK1eu5M0338RisWCxWJgzZ06F4Zpz5szBx8eHb775hoiICNzc3Bg+fDj5+fnMnTuX5s2b06BBA+6//36sVmv5sYuKinj00Udp3Lgx7u7u9OjRgxUrVphzoiIiUuvpTp6IiMgFePPNN9m3bx/t2rXj+eefB2Dnzp1ntcvPz+ett94iLi6OnJwchg4dyo033oiPjw/fffcdhw4dYtiwYfTu3ZuRI0cCcO+997Jr1y7i4uIIDg5m8eLFDBw4kO3btxMeHn5Jz1NERGo/hTwREZEL4O3tjZOTE25ubuVDNPfs2XNWu5KSEqZNm0ZoaCgAw4cP5+OPP+bEiRN4eHgQGRnJVVddxfLlyxk5ciTJycnMnj2b5ORkgoODAXj00UdZunQps2fP5qWXXrp0JykiInWCQp6IiEgVcnNzKw94AAEBATRv3hwPD48Kz6WnpwOwfft2rFYrrVq1qnCcoqIiGjZseGmKFhGROkUhT0REpAo5OjpW+LPFYjnnczabDYDc3Fzs7e3ZtGkT9vb2Fdr9PhiKiIhcKIU8ERGRC+Tk5FRhwZSq0KlTJ6xWK+np6fTp06dKjy0iIvWTVtcUERG5QM2bN2fdunUcPnyYU6dOld+Nq4xWrVoxZswYYmNjWbRoEUlJSaxfv54pU6bw7bffVkHVIiJS3yjkiYiIXKBHH30Ue3t7IiMj8ff3Jzk5uUqOO3v2bGJjY3nkkUeIiIhgyJAhbNiwgaZNm1bJ8UVEpH6xGIZhmF2EiIiIiIiIVA3dyRMREREREalDFPJERERERETqEIU8ERERERGROkQhT0REREREpA5RyBMREREREalDFPJERERERETqEIU8ERERERGROkQhT0REREREpA5RyBMREREREalDFPJERERERETqEIU8ERERERGROkQhT0REREREpA75P1CXO+97p0FhAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 900x600 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, axes = plt.subplots(2, 1, figsize=(9, 6), sharex=True)\n",
    "brainmass.viz.plot_timeseries(res['lfp'], ts=res['ts'], ax=axes[0])\n",
    "axes[0].set_title('Epileptor LFP proxy (x2 - x1)')\n",
    "brainmass.viz.plot_timeseries(res['z'], ts=res['ts'], ax=axes[1])\n",
    "axes[1].set_title('Slow permittivity variable z (seizure driver)')\n",
    "plt.tight_layout()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "151783f6",
   "metadata": {},
   "source": [
    "## Try it: vary the epileptogenicity `x0`\n",
    "\n",
    "`x0` controls whether the node is healthy or epileptogenic. Around `x0 ~ -2.0` lies the boundary: more negative values keep the node healthy (no large excursions), less negative values make it seize. We report the peak of the fast variable `x1` as a simple seizure indicator."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "48750c96",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-06-19T08:14:13.842823Z",
     "iopub.status.busy": "2026-06-19T08:14:13.842635Z",
     "iopub.status.idle": "2026-06-19T08:14:10.794996Z",
     "shell.execute_reply": "2026-06-19T08:14:10.794358Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "x0 = -2.4  ->  max(x1) = -1.514  (healthy)\n",
      "x0 = -1.8  ->  max(x1) = +2.016  (seizing)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "x0 = -1.6  ->  max(x1) = +2.047  (seizing)\n"
     ]
    }
   ],
   "source": [
    "for x0 in [-2.4, -1.8, -1.6]:\n",
    "    m = brainmass.EpileptorStep(in_size=1, x0=x0)\n",
    "    r = brainmass.Simulator(m, dt=0.1 * u.ms).run(\n",
    "        3000. * u.ms, monitors={'x1': 'x1'})\n",
    "    x1max = float(np.max(u.get_magnitude(r['x1'])))\n",
    "    state = 'seizing' if x1max > 0 else 'healthy'\n",
    "    print(f'x0 = {x0:.1f}  ->  max(x1) = {x1max:+.3f}  ({state})')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.13.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
