{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# DiscreteDP Example: Option Pricing " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Akihiro Sasaki**\n", "\n", "*Department of Economics, University of Tokyo*\n", "\n", "From Miranda and Fackler, Applied Computational Economics and Finance, 2002, Section7.6.4\n", "\n", "Julia translation of the [python version](http://nbviewer.jupyter.org/github/QuantEcon/QuantEcon.notebooks/blob/master/ddp_ex_MF_7_6_4_py.ipynb)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "using QuantEcon\n", "using PyPlot" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "T = 0.5 # Time expiration (years)\n", "vol = 0.2 # Annual volatility\n", "r = 0.05 # Annual interest rate\n", "strike = 2.1 # Strike price\n", "p0 = 2 # Current price\n", "N = 100 # Number of periods to expiration\n", "\n", "# Time length of a period\n", "tau = T/N\n", "\n", "# Discount factor\n", "beta = exp(-r*tau)\n", "\n", "# Up-jump factor\n", "u = exp(vol*sqrt(tau))\n", "\n", "# Up-jump probability\n", "q = 1/2 + sqrt(tau)*(r - (vol^2)/2)/(2*vol);" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# Possible price values\n", "ps = [u^i for i = -N:N] * p0\n", "\n", "# Number of states\n", "n = length(ps) + 1 # State n: \"the option has been exercised\"\n", "\n", "# Number of actions\n", "m = 2 # 0: hold, 1: exercise\n", "\n", "# Number of feasible state-action pairs\n", "L = n*m - 1; # At state n, there is only one action \"do nothing\"" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [], "source": [ "s_indices = Array{Int}(L+1)\n", "a_indices = Array{Int}(L+1)\n", "\n", "i = 1\n", "for s = 1:n\n", " for a = 0:m-1\n", " s_indices[i] = s\n", " a_indices[i] = a\n", " i += 1\n", " end\n", "end" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# Arrays of state and action indices\n", "s_indices, a_indices = s_indices[1:end-1], a_indices[1:end-1];" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# Reward vector\n", "R = zeros(m, n)\n", "R[1, :] = 0\n", "R[2,1:end-1] = strike - ps\n", "R = vec(R)[1:end-1];" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# Transition probability array\n", "Q = zeros(L, n)\n", "for i = 1:L-1\n", " if a_indices[i] == 0\n", " Q[i, min(s_indices[i]+1, length(ps))] = q\n", " Q[i, max(s_indices[i]-1, 1)] = 1 - q\n", " else\n", " Q[i, n] = 1\n", " end\n", "end\n", "Q[L, n] = 1;" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# Create a DiscreteDP\n", "ddp = DiscreteDP(R, Q, beta, s_indices, a_indices);" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": true }, "outputs": [], "source": [ "v = zeros(n)\n", "sigma = Array{Float64}(n)\n", "max_exercise_price = Array{Float64}(N)\n", "for i = 1:N\n", " v, sigma = QuantEcon.bellman_operator!(ddp::DiscreteDP, v, sigma)\n", " max_exercise_price[i] = ps[Int(sum(sigma))]\n", "end" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA/QAAAGJCAYAAAAkBnhUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3XdYU2f/BvA7YYQdQaYKivC6R9U6cWulrqpVWxwFHHUUtG6ldWtLbWu1ttXq6x5Yd2udP9x7i4p1A05AQNky8/z+sOQ1BpB9CNyf68pVczjjPgn05JvzDJkQQoCIiIiIiIiIdIpc6gBERERERERElH8s6ImIiIiIiIh0EAt6IiIiIiIiIh3Egp6IiIiIiIhIB7GgJyIiIiIiItJBLOiJiIiIiIiIdBALeiIiIiIiIiIdxIKeiIiIiIiISAexoCciIiIiIiLSQSzoiShPZDIZZs+eLXWMIpGRkQGZTIb58+dLHYWIiPKoLF2H8mL27NmQyWTl7tglISwsDDKZDGvXrpU6ik6pVq0avL29pY5Bb2FBT5QPa9euhUwmUz+MjIxQo0YN+Pr6IjIyskD7XLp0ab4vKElJSZg3bx4aNGgAExMTKJVKtGnTBuvXr4cQokA5AGDfvn2l6sNSgwYN4OTklOs5ubm5wc7ODhkZGSWYjIhIGrwOSef06dPo06cP7OzsoFAoUK1aNYwcORKPHj0q8D6Tk5Mxe/ZsHDt2rOiCliBvb2+N38e3fzfpf44dO6b1GllZWaFFixbYtGmT1PFIh+lLHYBIF82dOxfOzs5ISUnBqVOnsGzZMuzbtw/BwcEwMTHJ176WLl0Ka2vrPH/jGRkZiU6dOuHWrVvw8PCAr68vUlJSsGPHDnh5eWHfvn3YtGkT9PT08n1e+/btw2+//Zbth6lXr15BX79k/5cxaNAgTJs2DSdPnkTbtm21fh4WFoazZ8/C19e3xLMREUmJ16GS9csvv+DLL79E9erVMWbMGDg4OODWrVtYuXIltmzZgn379qFVq1b53m9ycjLmzJkDAGjfvr3Gz6ZPn45p06YVRfxipVAosHLlSq3lBXn/S0rVqlXx6tUrGBgYlPixx44di6ZNmwIAYmJisGXLFgwePBixsbHw8fEp8Tyk+/gJmKgAunbtivfffx8AMHz4cFSsWBE//fQT/vrrLwwYMKBYj+3l5YVbt25h165d+Oijj9TLx44di8mTJ+PHH39Eo0aNMHXq1CI9rhTftA8cOBB+fn4ICAjItqDfvHkzhBAYNGhQiWcjIpISr0Ml5/Tp0xg3bhxat26NAwcOaHxhMnr0aLi5uaFfv364efMmLC0ti+y4+vr6OvFltb6+PgYPHix1DKSkpMDQ0BBy+bsbIEvZgqBNmzbo16+f+vno0aNRvXp1BAQElKuCPj/vF+WOryBREejYsSMAIDQ0FEDOfc+ymkqGhYUBeN0X6ebNmzh+/Li6+dXb39C/6dy5czh48CC8vb01PkRl8ff3x3/+8x8sWLAAr169AvC/fmI//vgjFi1ahKpVq8LY2Bjt2rVDcHCweltvb2/89ttvAKDRHCxLdn0Xr169iq5du8LCwgJmZmbo1KkTzp07l+05nz59GhMmTICNjQ1MTU3Rp08fREVF5XiuAODo6Ii2bdti+/btSE9P1/p5QEAAXFxc0Lx5cwCvLw4zZsxA48aNoVQqYWpqinbt2uHEiRO5HgcABg8eDFdXV63l06dPz/YD1bp169CkSRMYGxvDysoKAwcOxNOnT995HCKi4sDrUPFchwBg3rx5kMlkWLdunVbrBxcXF3z//fcIDw/H8uXLNc7FzMwMISEhcHd3h6mpKSpVqoS5c+equySEhYXBxsYGADBnzhz1+WadY3bvoUwmg6+vL7Zt24Y6derA2NgYLVu2xI0bNwAAy5cvh6urK4yMjNC+fXv1+5zl5MmT6N+/P5ycnKBQKODo6Ijx48er36viIIRAhw4dYGNjg+fPn6uXp6WloX79+nBxcUFSUpJ6+dOnTzF06FB114a6deti9erVGvvMar7+xx9/YPr06ahcuTJMTEwQHx8PAIiNjcX48eNRrVo1KBQKVKlSBZ6enoiOjgaQfR/6iIgIDBkyBFWqVIFCoYCDgwN69eql9Rru378fbdq0gampKczNzdG9e3fcvHmzwK+PoaEhLC0ttT5rZGRkYN68eXBxcVF38fjqq6+QmpqqsV5OY0u83d89P38HQgjMnz8fVapUgYmJCTp06JDtOb548QKTJk1C/fr1YWZmBgsLC3Tt2hXXrl3TWC+n9ysoKAgymQyLFi3S2veZM2cgk8mwefPmd72E5V7p/9qPSAc8ePAAAFCxYsV8bbd48WKMGTMGZmZm+PrrrwEAdnZ2Oa7/999/AwA8PT2z/bm+vj4GDhyIOXPm4PTp0+jcubP6Z+vXr0dCQgJ8fHyQkpKCn3/+GR07dsSNGzdgZ2eHkSNH4tmzZwgMDMSGDRvemf3mzZto06YNLCwsMGXKFBgYGGD58uVo3749jh8/ri6ys4wZMwaWlpaYNWsWwsLCsHjxYvj6+mLLli25HmfQoEEYMWIEDh48iB49eqiX37hxA8HBwZg5c6Z6WWxsLFavXo0BAwZgxIgRiI+Px8qVK/HBBx/g0qVLqF+//jvPKy/mzJmDOXPmwMPDA8OHD8fz58+xZMkSnD9/HlevXoWFhUWRHIeIKK94HSqe61BycjIOHz6MNm3awNnZOdt1Pv30U4wYMQJ79uzRaCKfmZmJDz/8EC1atMD333+PAwcOYNasWcjIyMDcuXNhY2ODZcuWYfTo0ejTpw8+/vhjAK/Hj8nNyZMnsXv3bvXdXH9/f/To0QNTpkzB0qVL8cUXX+Dly5f4/vvvMXToUBw5ckS97bZt25CcnIzRo0ejYsWKuHDhAn755Rc8efIE27Zty/0Fz0VWofwmQ0NDWFhYQCaTYfXq1WjQoAFGjRqFnTt3AgBmzZqFmzdv4tixYzA1NQXwujtHixYt1F9c2NjYYP/+/Rg2bBji4+Mxbtw4jWPMmzcPhoaGmDRpElJTU2FoaIjExES0adMGt27dwtChQ9G4cWNER0dj9+7dePLkCaytrbM9h759++LmzZsYM2YMqlWrhufPnyMwMBCPHj1CtWrVAAAbNmyAl5cX3N3dsWDBAiQnJ2PZsmVo3bo1rl69ql4vNwkJCerX68WLFwgICEBwcDBWrVqlsd7w4cOxbt069OvXDxMnTsT58+fh7++vbh1TUHn5O5g5cybmz5+Pbt26oVu3brhy5Qq6dOmCtLQ0jX2FhITgzz//RP/+/eHs7IzIyEgsX74c7dq1wz///INKlSpprP/2+1WrVi24ublh06ZNGD9+vMa6mzZtgrm5OXr16lXgcy03BBHl2Zo1awQAcejQIREVFSUeP34s/vjjD1GxYkVhbGwsnjx5IoQQYtasWSK7P6+s7UNDQ9XL6tatK9q1a5en4/fu3VsAEC9fvsxxnZ07dwoAYsmSJUIIIUJDQwUAjXxCCHH+/HkBQIwfP169zMfHJ9vcQggBQMyaNUsji6GhoXjw4IF62bNnz4S5ublo27at1jl37txZqFQq9fLx48cLPT09ERsbm+s5v3jxQigUCjFgwACN5dOmTRMAxJ07d9TL0tPTRWpqqtb21tbWYsSIERrrARDz5s1TLxs0aJBwcXHROv7XX38t9PT01M/v378v5HK5WLBggcZ6QUFBQk9PT2s5EVFR4nVolkaW4r4OBQUFCQDiyy+/zHEdIYRo0KCBsLKyUj/38vISAMSYMWPUy1QqlejevbswNDQUUVFRQgghoqKitM4rS3bvIQChUCg03r/ly5cLAMLe3l7Ex8erl/v5+Wm918nJyVrH8ff3FzKZTDx8+DDXY2cn6zyze7i7u2usm5Vz48aN4ty5c0JPT0+MGzdOY51hw4YJBwcHER0drbHcw8NDKJVKdf6jR48KAKJ69epa5zRz5kwBQOzcuVMrb9b7n/U7uWbNGiGEEC9fvhQAxA8//JDjuSYkJIgKFSqIzz//XGN5RESEUCqVWsvflpX57YdcLhfffPONxrpZv3fDhw/XWD5p0iQBQBw5ckS9LKffn6pVqwovLy/187z+HTx//lwYGhqK7t27a6z31VdfCQAa+0xJSRGZmZkaxw0NDRUKhULMnTtX69yze7+yfi9u3bqlXpaWliasra01jkU5Y5N7ogLo3LkzbGxs4OjoCA8PD5iZmWHXrl2oXLlysR43ISEBAGBubp7jOlk/y2p2lqV3794a+Zo1a4bmzZtj3759+c6RmZmJ//u//0Pv3r1RvXp19XIHBwcMHDgQp06d0jr+iBEjNJoOtmnTBpmZmXj48GGux7K0tES3bt2we/dudZM8IQT++OMPvP/++6hRo4Z6XX19fRgaGgIAVCoVXrx4gczMTLz//vu4cuVKvs8zO1l3Fvr164fo6Gj1o3LlyqhevTqOHj1aJMchIsoNr0Mlcx3Ky/lm/fzt4wGAr6+v+t9Zd53T0tJw6NCh3E8wF506ddK4E5zVEqFv374aObOWh4SEqJcZGxur/52UlITo6Gi0atUKQghcvXq1QHmMjIwQGBio9fjuu+801hsxYgTc3d0xZswYfPbZZ3BxccG3336r/rkQAjt27EDPnj0hhNC4xrq7uyMuLk7rWu7l5aVxTgCwY8cONGzYEH369NHKmtNUfMbGxjA0NMSxY8fw8uXLbNcJDAxEbGwsBgwYoJFNT08PzZs3z/P1f+bMmerXaMuWLRgwYAC+/vpr/Pzzz+p1sv4mJkyYoLHtxIkTAQB79+7N07Gy866/g0OHDiEtLQ1jxozRWO/t1hHA6wERs/rAZ2ZmIiYmBmZmZqhZs2a2n7uye78++eQTGBkZaYz0f/DgQURHR5eKsRl0AZvcExXAb7/9hho1akBfXx92dnaoWbNmiQzqkXWhTkhIQIUKFbJdJ6cPH//5z3+01q1Rowa2bt2a7xxRUVFITk5GzZo1tX5Wu3ZtqFQqPH78GHXr1lUvd3Jy0lgva+CgnC6cbxo0aBB27dqFv/76CwMHDsSZM2cQFhaGL7/8UmvdNWvWYOHChbhz547GVHbZnX9B3Lt3DyqVCi4uLtn+/F0f+oiIigKvQyVzHXrzfHOTkJCgdb5yuVzjywYA6i+h3+6XnR9vn4dSqQTwetyZ7Ja/eX6PHj3CzJkzsXv3bq3zjouLK1AePT09ja4VuVm1ahVcXFxw7949nDlzRqO4i4qKQmxsLFasWIEVK1Zku/2bffABZNsN4sGDB+jbt28+zuB1YbpgwQJMnDgRdnZ2aNGiBXr06AFPT0/Y29sDeH39B/43XsXb8trdrn79+hqv1yeffIK4uDhMmzYNAwcOhI2NDR4+fAi5XK41to+9vT0qVKjwzpshuXnX30HWvt/+e7WxsdEa9FGlUuHnn3/G0qVLERoaiszMTPXPsuv+k937VaFCBfTs2RMBAQGYN28egNfN7StXrpzja02aWNATFUCzZs3UowtnJ6dvgN/8H11B1K5dG3/++SeuX7+e7ajvAHD9+nUAQJ06dQp1rKKW0/Q1Ig/zFffo0QNKpRIBAQEYOHAgAgICoKenBw8PD4311q5di6FDh+Ljjz/GtGnTYGNjAz09PcyfPx9PnjzJ9Rh5fc9UKhX09PSwf//+bLdhQU9EJYHXofwryHXI1dUV+vr66nPKTmpqKu7cuZPr+1GUcjqPd51fZmYmPvjgA7x48QJTp05FrVq1YGpqiqdPn8Lb2xsqlarYMmc5duyYelC3GzduoGXLluqfZR1/8ODB8PLyynb7t8cXePtub2GMGzcOPXv2xJ9//omDBw9ixowZ8Pf3x5EjR9CoUSN1vg0bNqiL/DcVZkaCTp06Yc+ePbhw4QK6d++uXp7T33Fe5PS3XpjPY2/79ttvMWPGDAwdOhTz5s2DlZUV5HI5xo0bl+3vU07vl6enJ7Zt24YzZ86gfv362L17N7744guOgJ9HLOiJikHWN5ixsbEadzCy+0Y1P/+z7tGjB/z9/bF+/fpsP0hlZmYiICAAlpaWcHNz0/hZ1jfLb7p7965Gs728ZrGxsYGJiQnu3Lmj9bPbt29DLpdr3SkoDIVCgX79+mH9+vWIjIzEtm3b0LFjR60L6vbt21GjRg3s2LFDY3nWQE+5sbS0RGxsrNbyt98zFxcXZGZmwsXFRevOCxFRacHrUNFch0xNTdGhQwccOXIEDx8+RNWqVbXW2bp1K1JTUzUGbgVeF6ghISEaXcPu3r0LAOpzLkzBll83btzA3bt3sW7dOo1BDQMDA0vk+OHh4RgzZgy6dOmiHhjN3d1d/Zra2NjA3NwcmZmZeb7jnx0XFxeN2RPyu+3EiRMxceJE3Lt3D++99x4WLlyIjRs3qlvm2draFipfdrJaFCYmJgIAqlatCpVKhXv37qF27drq9SIjIxEbG6vxe5jd55e0tDSEh4cXKEvWvu/du6fxOScqKkqrVcf27dvRoUMHrQH9YmNjcxx8MDsffvghbGxssGnTJjRv3hzJycn47LPPCpS/POLXHkTFIOt/+m9Ol5aUlIR169ZprWtqapptIZmdVq1aoXPnzlizZg327Nmj9fOvv/4ad+/exZQpU7S+Bf3zzz81plW7cOECzp8/j65du2pkAfDOPHp6eujSpQv++usvjWaDkZGRCAgIQOvWrYt8pPdBgwYhPT0dI0eORFRUVLZzz2d96/zmt8ynT5/GxYsX37l/FxcXxMTEaEzL8vTpU+zevVtjvb59+0Iul2POnDla32YLIRATE5Ov8yIiKg68DhXddWj69OkQQsDb21trerfQ0FBMmTIFDg4OGDlypNa2v/76q/rfQgj8+uuvMDAwQKdOnQBAPQ1eXl//wsjuGimE0Oi7XZw+//xzqFQqrFq1CitWrIC+vj6GDRumzqOnp4e+fftix44d2RbkeZliEHh9nb527Vq2I8HndBc6OTkZKSkpGstcXFxgbm6ublHg7u4OCwsLfPvtt9lOpZvXfNnJ+ltq2LAhAKBbt24AXs9C8aaffvoJADTu4ru4uGhNz7tixYoCt8bp3LkzDAwM8Msvv2i8Xm9nAV6/Z2+/ptu2bcv3NL76+voYMGAAtm7dirVr16J+/frvnO2B/od36ImKQZcuXeDk5IRhw4Zh8uTJ0NPTw+rVq2FjY4NHjx5prNukSRMsW7YM8+fPh6urK2xtbXPtM7R+/Xp06tQJvXr1wsCBA9GmTRukpqZi586dOHbsGD799FNMnjxZaztXV1e0bt0ao0ePRmpqKhYvXoyKFStiypQpGlkAYOzYsXB3d8+2WXuW+fPnIzAwEK1bt8YXX3wBfX19LF++HKmpqfj+++8L8rLlql27dqhSpQr++usvGBsbq6f3eVOPHj2we/dufPzxx+jatStCQkLw+++/o3bt2lrztr5t4MCB+Oqrr9CrVy+MGTMGSUlJWLp0KWrWrKkxn2qNGjUwZ84czJgxAyEhIfjoo49gZmaG0NBQ7Ny5E76+vtkOHENEVJJ4HSq661Dbtm3x448/YsKECWjQoAG8vb3h4OCA27dv47///S9UKhX27dun1b/YyMgIBw4cgJeXF5o3b479+/dj7969+Oqrr9TzzxsbG6NOnTrYsmULatSoASsrK9SrVw/16tUrsvxZatWqBRcXF0yaNAlPnz6FhYUFduzYkaexbHKTkZGBjRs3ZvuzPn36wNTUFGvWrMHevXuxdu1aVKlSBQDwyy+/YPDgwVi2bBm++OILAMB3332Ho0ePonnz5vj8889Rp04dvHjxAleuXMGhQ4fw4sWLd+aZPHkytm/fjv79+2Po0KFo0qQJXrx4gd27d+P3339XF81vunv3Ljp16oRPPvkEderUgb6+Pnbt2oXIyEj175+FhQWWLVuGzz77DI0bN4aHh4f672nv3r1wc3PT+AInJydPnlR/eZCV6/jx4/Dw8ECtWrUAvC7svby8sGLFCsTGxqJdu3a4cOEC1q1bh969e6NDhw7q/Q0fPhyjRo1C37598cEHH+DatWs4ePBgvu6Qv8nGxgaTJk1ST4fYrVs3XL16Ffv379faZ48ePTB37lwMGTIErVq1wo0bN7Bp06YCtWD09PTEkiVLcPToUSxYsKBA2cutEh5Vn0inZU35cfHixXeue/nyZdG8eXNhaGgonJycxE8//ZTtdEERERGie/fuwtzcXADI09RBCQkJYvbs2aJu3brC2NhYmJubCzc3N7F27VqNKUaE+N/ULD/88INYuHChcHR0FAqFQrRp00Zcu3ZNY92MjAwxZswYYWNjI2QymcaUNchmWpQrV64Id3d3YWZmJkxMTESHDh3EmTNn8vSaZU1hcvTo0Xeeb5bJkycLAOKTTz7J9ucqlUrMmzdPODk5CSMjI9G4cWOxb98+rSnpspu2Tggh9u/fL+rUqSMMDAxErVq1xObNm7Wmrcuybds24ebmJkxNTYWpqamoVauW8PX1Fffu3cvz+RAR5RevQ7M01i/J69CJEydEr169hLW1tTAwMBBOTk7i888/F2FhYVrrenl5CVNTU/HgwQPRpUsXYWJiIuzs7MSsWbO0pvk6c+aMaNKkiTA0NNQ4x5ymrfPx8dFY9ubrm935bdu2Tb3sn3/+EZ07dxZmZmbC2tpafP755+LatWsaU7jldOzs5DZtXdbv2ePHj4VSqRQ9e/bU2r5Pnz7C1NRUhISEqJdFRkYKHx8f4ejoKAwMDIS9vb3o1KmTWLFiRa7n9qaYmBjh6+srKleuLAwNDUWVKlWEl5eXejq8t6eti46OFj4+PqJWrVrC1NRUKJVK0bx5c7F161atfR89elS4u7sLpVIpjIyMhIuLi/D29haXLl3K9bXKbto6Q0NDUatWLfHNN9+ItLQ0jfXT09PFnDlzhLOzszAwMBCOjo7Cz89PpKSkaKyXmZkppk6dKqytrYWJiYlwd3cX9+/fz3Haurz8HWRmZoo5c+YIBwcHYWxsLNq3by+Cg4O19pmSkiImTpyoXs/NzU2cPXtWtGvXTuP/I+96v7LUrVtXyOVyjekt6d1kQhRgBAQi0hlhYWFwdnbGDz/8gEmTJkkdh4iIypnyeB3y9vbG9u3b1X2iiejdGjVqBCsrKxw+fFjqKDqFfeiJiIiIiIhIMpcuXUJQUJDGgI2UN+xDT0RERERERCUuODgYly9fxsKFC+Hg4IBPP/1U6kg6h3foiYiIiIiIqMRt374dQ4YMQXp6OjZv3gwjIyOpI+kc9qEnIiKiYuPv74+dO3fi9u3bMDY2RqtWrbBgwQLUrFkzx23Cw8MxceJEXLp0Cffv38fYsWOznTKJiIiovOMdeiIiIio2x48fh4+PD86dO4fAwECkp6ejS5cuSEpKynGb1NRU2NjYYPr06dlOMUVERESv8Q49ERERlZioqCjY2tri+PHjaNu27TvXb9++Pd577z3eoSciIspGuRsUT6VS4dmzZzA3N4dMJpM6DhEREYQQSEhIQKVKlSCXl+3Gc3FxcQAAKyurIt1vamoqUlNT1c9VKhVevHiBihUr8npPRESSE0IgJiYGVlZWGtd6hUIBhUJR4P2Wu4L+2bNncHR0lDoGERGRlsePH6NKlSpSxyg2KpUK48aNg5ubG+rVq1ek+/b398ecOXOKdJ9ERETFbdasWZg9e3aBty93Bb25ubn63y4uLjh16hRMTEwkTEREROVdfHw8HB0dNa5RZZGPjw+Cg4Nx6tSpIt+3n58fJkyYoH4eFxcHJycnPH78GBYWFkV+PCIiovzIutY/evQISqVSvbwwd+eBcljQZzW7s7e3x4MHD/Dtt99iyZIlEqciIiJCmW4a7uvriz179uDEiRPF0gohpyaLFhYWLOiJiKjUUCqVRXpdKtsd9XKxdOlSAMAvv/yCQ4cOSZyGiIiobBJCwNfXF7t27cKRI0fg7OwsdSQiIqIyo9wW9J06dcLo0aMBAEOGDEFsbKzEiYiIiMoeHx8fbNy4EQEBATA3N0dERAQiIiLw6tUr9Tp+fn7w9PTU2C4oKAhBQUFITExEVFQUgoKC8M8//5R0fCIiolKt3E1bFx8fD6VSibi4OOjp6aFhw4Z48OABNmzYgMGDB0sdj4iIyqE3r01lrXl4Tt0I1qxZA29vbwCAt7c3wsLCcOzYsVy3q1q1KsLCwvJ03LL8mhIRke4prutSuetD/yZTU1Ns2rQJCQkJ6Ny5s9RxiIiIypy83DdYu3ZtgbYjIiIq78p1QQ8AzZs3lzoCERERERERUb6V2z702QkNDcXs2bN5V4CIiIiIiIhKvXJ/hz5LYmIimjZtipiYGDg7O8PLy0vqSERERFRIz2JfsQ89ERGVWZLeoff390fTpk1hbm4OW1tb9O7dG3fu3HnndseOHUPjxo2hUCjg6uqabd+7/DIzM8PkyZMBAGPHjsXDhw8LvU8iIiKS1tVHL6WOQEREVGwkLeiPHz8OHx8fnDt3DoGBgUhPT0eXLl2QlJSU4zahoaHo3r07OnTogKCgIIwbNw7Dhw/HwYMHC51n0qRJaNWqFeLj4zFkyBCoVKpC75OIiIikc/khC3oiIiq7StW0dVFRUbC1tcXx48fRtm3bbNeZOnUq9u7di+DgYPUyDw8PxMbG4sCBA+88xrumC7h//z4aNmyI5ORkLF68GF9++WXBT4iIiCgPOMVa0ct6TTt+uw+H/bpKHYeIiMq54rrWl6pB8eLi4gAAVlZWOa5z9uxZrSnm3N3dcfbs2WzXT01NRXx8vMYjN66urli4cCEAYNq0abh9+3Z+ToGIiIhKkXvPExGbnCZ1DCIiomJRagp6lUqFcePGwc3NDfXq1ctxvYiICNjZ2Wkss7OzQ3x8PF69eqW1vr+/P5RKpfrh6Oj4ziwjR46Eu7s7UlJS4O/vn/+TISIiolKDze6JiKisKjUFvY+PD4KDg/HHH38U6X79/PwQFxenfjx+/Pid28hkMqxatQozZszAihUrijQPERERlawLYS+kjkBERFQsSsW0db6+vtizZw9OnDiBKlWq5Lquvb09IiMjNZZFRkbCwsICxsbGWusrFAooFIp8Z6pcuTLmzp2b7+2IiIiodLkYyoKeiIjKJknv0Ash4OtWXztQAAAgAElEQVTri127duHIkSNwdnZ+5zYtW7bE4cOHNZYFBgaiZcuWxRUTGRkZ+PXXX7Nt0k9ERESl242ncUhJz5Q6BhERUZGTtKD38fHBxo0bERAQAHNzc0RERCAiIkKjcPbz84Onp6f6+ahRoxASEoIpU6bg9u3bWLp0KbZu3Yrx48cXW85+/fphzJgxmD59erEdg4iIiIqerbkh0jMFgh7HSh2FiIioyEla0C9btgxxcXFo3749HBwc1I8tW7ao1wkPD8ejR4/Uz52dnbF3714EBgaiYcOGWLhwIVauXAl3d/diy/n5558DABYtWoRjx44V23GIiIioaDVysgQAXGI/eiIiKoMk7UMvhHjnOmvXrtVa1r59e1y9erUYEmWve/fuGD58OFauXAlvb29cv36d8wQTERHpgCZVLRF4PwEXwjjSPRERlT2lZpT70u6nn35CtWrV8PDhQ0yYMEHqOERERJQHjZwqAACuPHyJTNW7byQQERHpEhb0eWRubo5169app7T7+++/pY5ERERE71DDzgLmCn0kpmbgVni81HGIiIiKFAv6fGjbti0mTpwIABg7diwyMjIkTkRERES50ZPL0Lgq+9ETEVHZxII+n+bNm4fBgwdj//790NeXdAgCIiIiyoNmzlYAgIvsR09ERGUMK9J8MjIywoYNG6SOQURERHnUtNrrgv5C2AsIISCTySROREREVDR4h76Qzpw5gydPnkgdg4iIiHLQoIoShnpyRCWk4mFMstRxiIiIigwL+kJYsWIF2rRpg6FDh+ZpCj4iIiIqeUYGemjoqATw+i49ERFRWcGCvhDatm0LQ0NDBAYGYtmyZVLHISIiohxkNbu/GMqCnoiIyg4W9IVQq1YtLFiwAAAwefJk3Lt3T+JERERElJ2sgfF4h56IiMoSFvSF5Ovri44dOyI5ORleXl6cyo6IiKgUalLVEnIZ8DAmGZHxKVLHISIiKhIs6AtJLpdjzZo1sLCwwNmzZ/HDDz9IHYmIiIjeYm5kgNoOFgCAC2x2T0REZQQL+iLg5OSEJUuWAABmzZrFpvdERESl0P/mo2dBT0REZQPnoS8inp6eOHz4MNq0aQNXV1ep4xAREdFbmlWzwprTYbxDT0REZQYL+iIik8mwfv16qWMQERFRDpr+e4f+TmQC4pLToTQxkDgRERFR4bDJfTGJi4vDP//8I3UMIiIi+pe1mQLVbUwhBHDpIe/SExGR7mNBXwyuX7+O+vXr46OPPkJiYqLUcYiIiOhfzf6dj57N7omIqCxgQV8MqlatCplMhgcPHmDy5MlSxyEiIqJ/cT56IiIqS1jQFwOlUok1a9YAAH7//XccOHBA4kREREQEAE3/vUN/40kcktMyJE5DRERUOCzoi0nHjh0xduxYAMCwYcPw4gXvBBAREUmtiqUxKimNkKESCHoUK3UcIiKiQmFBX4z8/f1Rs2ZNPHv2DL6+vlLHISIiKvdkMpl6tPvz7EdPREQ6jgV9MTIxMcH69euhp6eHzZs3Y8eOHVJHIiIiKveymt1fZD96IiLScZyHvpg1a9YMX331FW7duoV27dpJHYeIiKjca/7vHforj14iLUMFQ33e3yAiIt3Egr4EzJ49GzKZDDKZTOooRERE5Z6rrRksTQzwMjkdN57GoUlVS6kjERERFQi/ki4BcrlcXcwLIRASEiJxIiIiovJLJpOx2T0REZUJLOhLUGJiIj799FM0aNCART0REZGEmlevCAC4wIHxiIhIh7GgL0HGxsaIjIxEUlISvL29kZmZKXUkIiKicimrH/3F0BfIVAmJ0xARERUMC/oSpKenh7Vr18LMzAwnT57EokWLpI5ERERULtV2sICZQh8JqRm4FR4vdRwiIqICYUFfwpydndWF/Ndff43g4GCJExEREZU/enIZ3q/2ejA8NrsnIiJdxYJeAsOGDUP37t2RlpaGzz77DGlpaVJHIiIiKneaO7/uR38+NEbiJERERAXDgl4CMpkMK1euRMWKFREUFIR58+ZJHYmIiKjcafZvP/oLoS8gBPvRExGR7mFBLxF7e3v8/vvvcHR0RPv27aWOQ0REVO7Ur6yEkYEcL5PTcf95otRxiIiI8o0FvYT69euH27dvo1OnTlJHISIiKncM9eVo7PS6H/159qMnIiIdxIJeYiYmJup/v3z5UsIkRERE5c//+tGzoCciIt3Dgr6UWLVqFapWrYrDhw9LHYWIiKjc+F8/+hj2oyciIp3Dgr6UuHz5MhISEuDt7Y3Y2Fip4xAREZULjZwqwFBPjsj4VDx6kSx1HCIionxhQV9K/PDDD3BxccGTJ0/w5ZdfSh2HiIioXDAy0ENDRyUA4HwIm90TEZFuYUFfSpiammL9+vWQy+VYv349du3aJXUkIiKiQvP390fTpk1hbm4OW1tb9O7dG3fu3HnndseOHUPjxo2hUCjg6uqKtWvXFlvGrGb37EdPRES6hgV9KdKqVStMmTIFADBy5Eg8f/5c4kRERESFc/z4cfj4+ODcuXMIDAxEeno6unTpgqSkpBy3CQ0NRffu3dGhQwcEBQVh3LhxGD58OA4ePFgsGZv9OzDehbCYYtk/ERFRcZGJcjYCTHx8PJRKJeLi4mBhYSF1HC2pqalo1qwZrl+/jl69emHXrl2QyWRSxyIiomJU2q9NRSkqKgq2trY4fvw42rZtm+06U6dOxd69exEcHKxe5uHhgdjYWBw4cCBPx8nPa5qYmoGGc/4PmSqBM9M6olIF47yfEBERUR4U17Wed+hLGYVCgfXr18PIyAi1atVCZmam1JGIiIiKTFxcHADAysoqx3XOnj2Lzp07ayxzd3fH2bNnc9wmNTUV8fHxGo+8MlPoo16l1x+uLrDZPRERFaO3r1WpqamF2h8L+lKoYcOGCA0NxXfffQd9fX2p4xARERUJlUqFcePGwc3NDfXq1ctxvYiICNjZ2Wkss7OzQ3x8PF69epXtNv7+/lAqleqHo6NjvrL9rx89m90TEVHxcXR01Lhe+fv7F2p/rBZLKXt7e/W/MzIyIJfLIZfz+xciItJdPj4+CA4OxqlTp4p8335+fpgwYYL6eXx8fL6K+ubOFfHfk6Ec6Z6IiIrV48ePNZrcKxSKQu2PFWIpd/v2bbRq1Qq//PKL1FGIiIgKzNfXF3v27MHRo0dRpUqVXNe1t7dHZGSkxrLIyEhYWFjA2Dj7/u0KhQIWFhYaj/xo6mwFmQwIiU7C8/iUfG1LRESUV29fq1jQl3FHjx7FxYsXMW3aNNy+fVvqOERERPkihICvry927dqFI0eOwNnZ+Z3btGzZEocPH9ZYFhgYiJYtWxZXTCiNDVDb/vWXAJy+joiIdAUL+lJu1KhRcHd3R0pKCjw9PZGRkSF1JCIiojzz8fHBxo0bERAQAHNzc0RERCAiIkKjL7yfnx88PT3Vz0eNGoWQkBBMmTIFt2/fxtKlS7F161aMHz++WLM2r/66H/25EPajJyIi3cCCvpSTyWRYtWoVKlSogIsXLxZ60AQiIqKStGzZMsTFxaF9+/ZwcHBQP7Zs2aJeJzw8HI8ePVI/d3Z2xt69exEYGIiGDRti4cKFWLlyJdzd3Ys1a4vqr+ej5x16IiLSFZyHXkds3rwZAwcOhL6+Ps6dO4cmTZpIHYmIiIqIrl6bSrOCvKYvk9LQaF4gAODS9M6wNitcv0YiIqIsnIe+nPPw8ED//v2RkZEBT09PpKRwwB4iIqKiZGlqiFr25gA4Hz0REekGFvQ6QiaTYenSpbC3t4epqSnO3X4idSQiIqIyp3nWfPTsR09ERDqABb0Osba2xuEjR9Hsy2UYsuUOB+0hIiIqYs3Zj56IiHQIC3odU6d2LZgZGUAIYNK2a0hISZc6EhERUZnR7N879LcjEvAiKU3iNERERLljQa+Dpveog8oVjPDk5St86PdfqeMQERGVGdZmCvzH1gwA+9ETEVHpx4JeB5kp9DGsrgGEUOGpoir81++ROhIREVGZkTUf/flQdm0jIqLSjQW9jhrasy1qitcD4/1+OR73H4VLnIiIiKhsaO78bz/6EN6hJyKi0o0FvQ7bNsMTsvgIyIyV6P/9DqhUKqkjERER6bysO/S3IuIRl8yxaoiIqPRiQa/DlOam+P7juhCZGXhp5oxpy7ZLHYmIiEjn2Zobobq1KYQALoTxLj0REZVe+lIHoMLp37kFNh//L66kV8KW+4DH7RA0rlVd6lhEREQ6rXl1K4REJ2HT+YcIj3uV43oVTRXoVt8eMpmsBNMRERG9xoK+DAiY7oV649ciw6Iyvv7rFvbVdOYHCyIiokJoUb0iNl94jGN3onDsTlSu6y7/rAnc69qXUDIiIqL/YUFfBhgpDPHfoa0xalcYbr0U2Hj+ET5rUVXqWERERDqraz0HXHOLQ2R8So7r3HuegLuRiTh5L4oFPRERSYIFfRnR4f068Es1xuy//8G3e2+htas1nK1NpY5FRESkkwz15ZjZs06u6xy8GYGRGy7jHEfDJyIiiXBQvDLEs2U1uLlWxKv0TPScvxWpaRyZl4iIqLg0d7aCTAbcf56IqIRUqeMQEVE5xIK+DJHLZZjpXh0iNRmJRrb47Jt1UkciIiIqsyqYGKKWvQUA4HxojMRpiIioPGJBX8bUdLTFx9UyAQDnk63x57GLEiciIiIqu1r8O2f9uRAW9EREVPIkLehPnDiBnj17olKlSpDJZPjzzz9zXf/YsWOQyWRaj4iIiBJKrBsWjvkEyoRQyPQMMGHrNSQk5TzdDhERERVci+oVAQBnH7CgJyKikidpQZ+UlISGDRvit99+y9d2d+7cQXh4uPpha2tbTAl1k1wux9bJfSBexUFl4YBP5rHpPRERUXHI6kf/ICoJzxNyHhGfiIioOEha0Hft2hXz589Hnz598rWdra0t7O3t1Q+5nD0H3lazWmUMa2ACAPhHVQmr/z4ucSIiIqKyp4KJIWpn9aPnaPdERFTCdLISfu+99+Dg4IAPPvgAp0+fznXd1NRUxMfHazzKi5nDesMuORQyuR4WnnqOpNQMqSMRERGVOVnN7tmPnoiISppOFfQODg74/fffsWPHDuzYsQOOjo5o3749rly5kuM2/v7+UCqV6oejo2MJJpbejq8HoIKhCkkyE3y775bUcYiIiMocDoxHRERS0amCvmbNmhg5ciSaNGmCVq1aYfXq1WjVqhUWLVqU4zZ+fn6Ii4tTPx4/flyCiaVXxa4ilnq1BABsOv8Ix+48lzgRERFR2dKM/eiJiEgiOlXQZ6dZs2a4f/9+jj9XKBSwsLDQeJQ3rVysMcStGgBg5KqTCHvGop6IiKiosB89ERFJRecL+qCgIDg4OEgdo9Sb+mEtGLyKQarcCP38t0kdh4iIqExhP3oiIpKCpAV9YmIigoKCEBQUBAAIDQ1FUFAQHj16BOB1c3lPT0/1+osXL8Zff/2F+/fvIzg4GOPGjcORI0fg4+MjSX5dYmSghzldXSBUmYg2rYavft8hdSQiIqIyI6sf/VkW9EREVIIkLegvXbqERo0aoVGjRgCACRMmoFGjRpg5cyYAIDw8XF3cA0BaWhomTpyI+vXro127drh27RoOHTqETp06SZJf1wz80A0N9cMBAJtuZ+DG/Ufv2IKIiIjyorlzRchkQEhUEp7Hsx89ERGVDJkQQkgdoiTFx8dDqVQiLi6uXPanT3qVggYT1yPTojLM4h/i+q+jIJfrfM8LIiKdVt6vTcVBite0+5KTuPksHsNaO6NupYIfU19PjnY1bKA0NijCdEREJKXiui7pF9meSCeYGhvh10FNMWpXGBItqsJ30WYsnThI6lhEREQ6r0X1irj5LB6rToUWel+936uExR6NiiAVERGVZSzoy6Gubo3Q5shVnEqyw75wY4RGJcDZxlzqWERERDptiFs1hMe9QkJKRoH3kZKeiYthL3HsbhRUKgG5XFaECYmIqKxhQV9OrZnmibazdyAcZpiy4wb+GNESevzQQEREVGBVLE2wdFCTQu0jPVOFhnP+D7HJ6bgdkYA6hWi6T0REZR87T5dTBgb62DqxJ0wN9XAx7CVWnQqROhIREVG5Z6AnR9NqHDGfiIjyhgV9OeZoZYJZPesCABbsu4U9p65InIiIiIhaurye0/7sAxb0RESUOxb05Vz/96vASS8OmZDhy4ArSEzmVDtERERSaln9dUF/PjQGmapyNRkRERHlEwv6ck4mk2HpkNZQpSQg08IBA+avlzoSERFRuVa3kgXMFfpISMnAP8/ipY5DRESlGAt6Qj1XJ3jVMQQAXM+wx4b9pyROREREVH7p68nRzDmrH320xGmIiKg0Y0FPAIB5Iz6GTVIoZHI9zNx7DzFxCVJHIiIiKrfYj56IiPKCBT2p7fjqU4iklxBmtug3d4PUcYiIiMqtFv/2o78Y9hIZmSqJ0xARUWnFgp7UnBxsMK7V6w8QoQZVse3EDYkTERERlU91HCygNDZAYmoGbjyNkzoOERGVUizoScP4AV1RxyAKAPDT6eeIe5UucSIiIqLyRy6Xobkz56MnIqLcsaAnLdtnDIaztSnC41IwZ/dNqeMQERGVS+xHT0RE78KCnrSYGOpj4ScNIZcBO68+xQ+bA6WOREREVO5kFfSXwl4iLYP96ImISBsLespWYydL9HRRAAB+PReDW6FPJU5ERERUvtSwNYeVqSFepWfi+pNYqeMQEVEpxIKecjR/gBvk8c8gMzbHpz/+BZWKdweIiIhKilwuQ4vq//ajZ7N7IiLKhr7UAaj0Mjc1xqJPG2Hs308Qb14VE37ZisVfekgdi4iIqNxoUb0i9t2IwM6rT5GQmlGofb1f1RJd6toXUTIiIioNWNBTrnq1ex8BR4NwPsUBux7qYUDwfTSv5yp1LCIionLBzdUaABAanYQVJ0IKta9Vchkuft0ZVqaGRRGNiIhKARb09E4bvvZGvbErkaZ0hNfSQwhe4gx9fT2pYxEREZV5LjZmWPzpe7gVHl+o/ey+9gzhcSk4FxKDbvUdiigdERFJjQU9vZOhgT7+O7wNPDfdRoqFI8b+uhNLx/WXOhYREVG50LtRZfRuVLlQ+0jNUGHtmTCceRDNgp6IqAzhoHiUJ+2a1EW3Sq8AAIeizXAvMkHiRERERJRXrf6dAu8MB9cjIipTWNBTni2dOAjta9ogLUOFCVuvIT2To94TERHpgubVK0IuA0KikhARlyJ1HCIiKiIs6CnPZDIZFvRtAKWxAW48jcOk1YFSRyIiomISGxuLlStXws/PDy9evAAAXLlyBU+fPpU4GRWE0tgA9SorAQBnQ6IlTkNEREWFBT3li52FEb7q4gwA+PNeKrYeOidxIiIiKmrXr19HjRo1sGDBAvz444+IjY0FAOzcuRN+fn4Sp6OCavlvs3vOaU9EVHawoKd8+6SFK6wSQiCT62Hqrn8Qm5AkdSQiIipCEyZMgLe3N+7duwcjIyP18m7duuHEiRP52teJEyfQs2dPVKpUCTKZDH/++ec7t/ntt99Qu3ZtGBsbo2bNmli/fn2+z4G0tXJ5PQUe+9ETEZUdLOgp32QyGbZN6w+RHAthbod+8zZIHYmIiIrQxYsXMXLkSK3llStXRkRERL72lZSUhIYNG+K3337L0/rLli2Dn58fZs+ejZs3b2LOnDnw8fHB33//na/jkram1SyhL5fhyctXePwiWeo4RERUBFjQU4G4ONrji/df98W7L3fE7zsPS5yIiIiKikKhQHy89rznd+/ehY2NTb721bVrV8yfPx99+vTJ0/obNmzAyJEj8emnn6J69erw8PDAiBEjsGDBgnwdl7SZGOrjPccKAIAzD9iPnoioLGBBTwU21bMHKqeEAgC+O/oM4dEvJU5ERERF4aOPPsLcuXORnp4O4HXLrEePHmHq1Kno27dvsR47NTVVo5k/ABgbG+PChQvqPDltFx8fr/EgbZy+johIWm9fq1JTUwu1Pxb0VCg7ZwwGEqMAUyv0998idRwiIioCCxcuRGJiImxtbfHq1Su0a9cOrq6uMDc3xzfffFOsx3Z3d8fKlStx+fJlCCFw6dIlrFy5Eunp6YiOzvmusr+/P5RKpfrh6OhYrDl1Vcs3+tELISROQ0RU/jg6Ompcr/z9/Qu1P/0iykXllF3FCvi6kxO+OZeMJwaOCPwnEh/UsZM6FhERFYJSqURgYCBOnz6Na9euITExEY0bN0bnzp2L/dgzZsxAREQEWrRoASEE7Ozs4OXlhe+//x5yec73Ifz8/DBhwgT18/j4eBb12WjkVAEKfTmiElLxICoRrrbmUkciIipXHj9+DAsLC/VzhUJRqP2xoKdC+7xPR0QrbmH5iRD47byOxk5tUdGscL+YREQkPTc3N7i5uZXoMY2NjbF69WosX74ckZGRcHBwwIoVK2Bubp5r/32FQlHoD0XlgZGBHt6vZonT92Nw5kEMC3oiohJmYWGhUdAXVqGb3CcmJrLPGmFClxqoaWeO6MQ0+Kw9BZVKJXUkIiIqoLFjx+LXX3/VWv7rr79i3LhxJZLBwMAAVapUgZ6eHv744w/06NEj1zv0lHfq6evusx89EZGuK9CVMTQ0FN27d4epqSmUSiUsLS1haWmJChUqwNLSsqgzkg5Q6Ovhp08bQiZUOPckBVOWbpc6EhERFdCOHTvQunVrreWtWrXC9u35+/97YmIigoKCEBQUBOD1Z4igoCA8evQIwOum8p6enur17969i40bN+LevXu4cOECPDw8EBwcjG+//bYQZ0RvavnvwHhnQ2KgUrEfPRGRLitQk/vBgwdDCIHVq1fDzs4OMpmsqHORDqpbSYnGhuG4nF4Z2x7IMOBWCJrUri51LCIiyqeYmBiYm2s3xbawsMh1YLrsXLp0CR06dFA/z+rn7uXlhbVr1yI8PFxd3ANAZmYmFi5ciDt37sDAwAAdOnTAmTNnUK1atYKdDGlpUFkJM4U+4l6lo/bMAyjpj3HVKppi66iWsDAyKNkDExGVQQUq6K9du4bLly+jZs2aRZ2HdNym6d6oN241MpSO+OzXg7ixZCT09NhEkohIl7i6umL//v3w9fXVWL5//35Ur56/L2rbt2+f62jqa9eu1Xheu3ZtXL16NV/HoPzR15OjRwMH/HHxMVIzSr6L3O2IBBy7E4WPGlYq8WMTEZU1BSromzZtisePH7OgJy1GCkMsH9IKQ7fcQ7K5E0b+uAkrp34mdSwiIsqHCRMmwNfXF1FRUejYsSMA4PDhw1i4cCEWL14scToqCv4f18eXnf+DzBJucr/02AMEnH+Esw+iWdATERWBAhX0K1euxKhRo/D06VPUq1cPBgaaTaYaNGhQJOFIN3Vq1gAdD13B0XgjBD43xZGLN9GxaV2pYxERUR4NHToUqamp+OabbzBv3jwAQLVq1bBs2TKN/u6ku2QyGRyUxiV+3M61bRFw/hFOc0A+IqIiUaCCPioqCg8ePMCQIUPUy2QyGYQQkMlkyMzMLLKApJv+O+Uz1B/zO14pq2HkmtMIblgDCkP2lSMi0hWjR4/G6NGjERUVBWNjY5iZmUkdicqAZs4VoS+X4dGLZDx+kQxHKxOpIxER6bQCdW4eOnQoGjVqhLNnzyIkJAShoaEa/yXS19fDhjEfAmnJSLeojBWnwqSOREREBWBjY8NinoqMmUIfDR0rAADOPMjfAItERKStQHfoHz58iN27d8PV1bWo81AZ8n4dF/w4QB+TdgTj50P30KGmLepVVkodi4iIstG4cWMcPnwYlpaWaNSoUa4z2Fy5cqUEk1FZ4+ZSEZcfvsTp+zH4tKmT1HGIiHRagQr6jh074tq1ayzo6Z36vu+Ew3eisT84AmMDLuPP0S1gYcbmdUREpU2vXr2gUCgAAL1795Y4DZVlrVytseTIfZx5EKPurklERAVToIK+Z8+eGD9+PG7cuIH69etrDYr30UcfFUk40n0ymQzze9fD6TsRCIl5hU/mb8CB70ZKHYuIiN4ya9YsAK/nge/QoQMaNGiAChUqSJyKyqJGThVgZCBHdGIq7kYmoqa9udSRiIh0VoEK+lGjRgEA5s6dq/UzDopHb6topkB/53SsumuAW6ISVu0+jmEftZM6FhERZUNPTw9dunTBrVu3WNBTsVDo66FpNSucvBeNMw+iWdATERVCgQbFU6lUOT5YzFN2ZgztDfvkEMhkcsz7vzBExsRKHYmIiHJQr149DnJLxaqVizUAcPo6IqJCKlBBT1QQ278eAJEYA5hZo983m6WOQ0REOZg/fz4mTZqEPXv2IDw8HPHx8RoPosJyc60IADgfEoOMTJXEaYiIdFeBmtxn19T+TTNnzixQGCrbqthZY2o7Byy4lILHhk74MWA/Jg3sKnUsIiJ6S7du3QC8HhPnzQHLsgYwY2s8Kqy6lZSwMNJHfEoGbjyNQyMnS6kjERHppAIV9Lt27dJ4np6ejtDQUOjr68PFxYUFPeXoi/4fYPu53xBiUA2/nHuJj9tFonplO6ljERHRG44cOcKRx6lY6cllaOlSEQdvRuLMgxgW9EREBVSggv7q1atay+Lj4+Ht7Y0+ffoUOhSVbTtmeaLxtK2AhT2+CwzDci9bfnAkIipF2rdvL3UEKgfcXK1x8GYkTt+Phk8HToVMRFQQRdaH3sLCAnPmzMGMGTOKapdURllamCNg7IfQl8vwf7ejsfvaM6kjERERgKSkJHzxxReoXLkybGxs4OHhgaioKKljURmVNTDepYcvkZLObhxERAVRoDv0OYmLi0NcXFxR7pLKqJa1KmNMx2QsOnQXM/+6ifp2xqjuYCV1LCKicm3GjBlYv349Bg8eDIVCgc2bN2PEiBFaXe2IioKLjSlszRV4npCKsZuvooKJgdSRtHSt54AOtWyljkFElKMCFfRLlizReC6EQHh4ODZs2ICuXTnIGeXNFx1ccODGE9yKTMZH87fg+i8jIZdz4gUiIqns2rULa6e2PkgAACAASURBVNasQf/+/QEAnp6eaNGiBTIyMqCvX6T3AIggk8nQroYNtl1+gv/7J1LqONk6EByBKzM+gL4eP58QUelUoKvzokWLNJ7L5XLY2NjAy8sLfn5+RRKMyj4DPTm+bF4BI3e+RKK5E3wW/YFlEwdKHYuIqNx68uQJ3Nzc1M+bNGkCAwMDPHv2DE5OThImo7JqatdaqGlvjtSM0jd13e/HH3AUfiIq9QpU0IeGhhZ1DiqnPmzVCG0PX8XJJDvse6bAqaA7aP1eTaljERGVSyqVCgYGms2e9fX1OU0dFRtrMwWGt6kudYxs3XgShwM3I3D6fjQLeiIqtdh+iCS3eponFHEPITMwwtAVx5CWniF1JCKickkIgU6dOqFx48bqR3JyMnr27KmxjKg8cPvP60H7Tt2PljgJEVHO8nyH/uOPP8batWthYWGBjz/+ONd1d+7cWehgVH4YGOhjzahO8Fh3HWkWVeD93QYEzBgidSwionJn1qxZWst69eolQRIi6bV2fV3QX3kYi1dpmTA21JM4ERGRtjwX9EqlUj1XuFKpLLZAVD61algTvapcwt/PTXA63hJ7Tl1Fj9aNpI5FRFSuZFfQE5VX1SqaoJLSCM/iUnAx7AXa1rCROhIRkZY8F/Rr1qzJ9t9EReXncQNwzHcpEiyc8cPJ5/igRSYU+vw2nIiIiEqeTCaDm6s1tl1+gtP3o1nQE1GpxD70VGrI5XLsmNoXViYGeBiXgZ8P3ZM6EhEREZVjrdmPnohKuQIV9DExMfDx8UGdOnVgbW0NKysrjQdRQdWo6oBvP24A4PV0MRdCeAElIqL/Z+/Ow6Iq+zeA32cGGPYBZRcEUyGXQNQ0pHKDCP1ZViqmJi6UmlpKaZrlUq+vZa4lLq/hmqmZW66JGyhqboyKu0KCCigqO7LMnN8f5hThAjhwZpj7c13n0jlzlnueAZ75zjnnOUTSaNfwQUF/9mYO7uYXS5yGiKi8Kt227r333sOVK1cwePBgODs7a6+tJ9KF15u74G3/etiQcAO9Z2/HH193h6OdrdSxiIiIyMg42ijg42yDixm5OHz1Drr6ukodiYiojCoV9AcOHMDBgwfh5+en6zxEAIBR7d2x/uAZaKzqoud/fsb+GUOljkREZDRWrFiBsLAwKBSKMvOLi4uxZs0a9O/fX6JkRDUvsJEDLmbk4uCVTBb0RKR3qnTK/fPPP4/CwkJdZyHSqu/igNHtHpzm9qeJB+au3SVxIiIi4zFw4EBkZ2eXm5+bm4uBA3lbUTIuLzeuCwCI53X0RKSHqlTQz58/HxMmTEBsbCzu3LmDnJycMlNFxcXFoVu3bnBzc4MgCNi0adNT19m/fz9atmwJhUKBRo0aYdmyZVV5CWQARr0bCs+iZADArPjbSE2/I3EiIiLjIIriIy+nu379Om9dS0anTYO6MJEJSLlbgJQ7BVLHISIqo0oFvZ2dHXJyctCpUyc4OTnB3t4e9vb2sLOzg729fYW3k5+fDz8/P0RFRVVo+eTkZHTt2hUdO3aESqXCqFGjEBERgd9//70qL4MMwK8T3wNyMyBY2uGdaWuljkNEVKv5+/ujZcuWEAQBnTt3RsuWLbWTn58fXnnlFQQFBUkdk6hGWStM4F/fDgAQf5VH6YlIv1TpGvq+ffvC1NQUP//88zMNihcaGorQ0NAKL79w4UI0aNAAM2fOBAA0adIEBw8exOzZsxESElKlDKTfHO2V+Cq0Ib6My8YtC09Mjt6MyYPflDoWEVGt1L17dwCASqVCSEgIrK2ttc+ZmZnBy8sL77zzjlTxiCTTrqEDjv15D/FXMvFum/pSxyEi0qpSQZ+YmIiEhAT4+PjoOs8THT58uNyRgZCQEIwaNeqx6xQVFaGoqEj7uDKXBJB+6P9/7bH24AKcRX2suSrgw5z7cLI1lzoWEVGtM2nSJACAl5cXwsLCYG7Ov7VEwIP70c/dcxkHLmdi7u7LUsepFiZyAW/514ObnYXUUYioEqpU0Ldu3Rqpqak1XtCnp6fD2dm5zDxnZ2fk5OSgsLAQFhbl/wBNmzYNU6ZMqamIVE3WfDkAb847gKR7wLgNZxAd3pq3SyQiqibh4eHIysrCTz/9hKtXr2LMmDGoU6cOTp48CWdnZ9SrV0/qiEQ1qoWHHawVJsguLMHs3ZekjlNtVKlZWNy/tdQxiKgSqlTQjxw5Eh9//DHGjBmDF154AaampmWe9/X11Uk4XRg/fjwiIyO1j3NycuDh4SFhIqoKGysLLBwQgP/74SD2XriFtcdS0ZunvBERVYvTp08jKCgISqUSf/75J95//33UqVMHGzZsQEpKClasWCF1RKIaZSqX4Yc+/th9LkPqKNWioFiNjQk3cOhKJkrUGpjKqzTMFhFJoEoFfVhYGABg0KBB2nmCIGhHxVWr1bpJ9y8uLi7IyCj7hzQjIwO2traPPDoPAAqFotx9dMkweTvbYMxrPpi6/TzGrzsOL8tivNS8kdSxiIhqndGjR2PAgAGYPn06bGxstPO7dOmCPn36SJiMSDodfZzQ0cdJ6hjVQqMREXvpNu7mFyMhJQttGtSROhIRVVCVvn5LTk4uNyUlJWn/rS4BAQHYs2dPmXkxMTEICAiotn2Sfhn0cgMoclIhyhUYsGAPSkqr58sjIiJjdvz4cQwZMqTc/Hr16iE9PV2CRERUnWQyAe0a1gUAHLx8W+I0RFQZVSroPT09nzhVVF5eHlQqFVQqFYAHXxSoVCqkpKQAeHC6fP/+/bXLDx06FElJSRg7diwuXLiA+fPn45dffsHo0aOr8jLIAMllAhYPfgWa4kLct3HH4Ok/SR2JiKjWUSgUjxxE9tKlS3B0dJQgERFVt1caOwAADl7hrfmIDEmVL5BZuXIlAgMD4ebmhmvXrgEA5syZg82bN1d4G8ePH4e/vz/8/f0BAJGRkfD398fEiRMBAGlpadriHgAaNGiAbdu2ISYmBn5+fpg5cyZ+/PFH3rLOyLzaqhlCXQoAALF3ldh5+LTEiYiIapc33ngDX331FUpKSgA8uKwuJSUFn332GW9bR1RLvdz4wZd1p65nI+d+icRpiKiiqlTQL1iwAJGRkejSpQuysrK018zb2dlhzpw5Fd5Ohw4dIIpiuWnZsmUAgGXLlmH//v3l1klISEBRURGuXr2KAQMGVOUlkIFb8Gk/WGUnQzAxxYifjqLgfrHUkYiIao2ZM2ciLy8PTk5OKCwsRPv27dGoUSPY2Nhg6tSpUscjompQz84CzzlYQa0RcfjqHanjEFEFVamg/+GHH7B48WJMmDABcrlcO79169Y4c+aMzsIRPY5MJsOqUV2huZ+LUhtX9JnKEZeJiHRFqVQiJiYGW7duxffff48RI0Zg+/btiI2NhZWVldTxiKiavPzwtPvLPO2eyFBUaZT75ORk7Wny/6RQKJCfn//MoYgqooVPA/TxPoY1KYCq2Bknr91FS0+OykpEpCuBgYEIDAwEAGRlZUmchoiqW2AjB6w4fA3xvI6eyGBU6Qh9gwYNtAPZ/dPOnTvRpEmTZw5FVFHThvVEKwcNIMjw6brTKCzmqPdERM/q22+/xdq1a7WPe/Xqhbp166JevXo4deqUhMmIqDoFNKwLuUxAUmY+bmQVSh2HiCqgSgV9ZGQkhg8fjrVr10IURRw9ehRTp07F+PHjMXbsWF1nJHosQRCw5MPX4WyrQFJmPr7deUHqSEREBm/hwoXw8PAA8OD2sDExMdixYwdCQ0MxZswYidMRUXWxNTeFn7sSAG9fR2QoqnTKfUREBCwsLPDFF1+goKAAffr0gZubG+bOnYvevXvrOiPREyktTTG9hx/ClxzFskN/wqkkHR++01nqWEREBis9PV1b0G/duhW9evXCa6+9Bi8vL7Rt21bidERUnV5u7IiTKVk4cDkTYS/WlzoOET1FpY/Qi6KIlJQUvPPOO7h8+TLy8vKQnp6O69evY/DgwdWRkeip2ns7orGQDgCYHpuGG7fvSZyIiMhw2dvbIzU1FcCDy+mCgoIAPPgM8PDONkRUO73c6MHAeIeu3oFGI0qchoiepkoFfaNGjbQdvaWlJZycnHQejKiyVn7yNpB7G7C0R4+pa6SOQ0RksN5++2306dMHwcHBuHPnDkJDQwEACQkJaNSokcTpiKg6+de3g5WZHHfzi3EuLUfqOET0FJUu6GUyGRo3bow7d3h/StIvLg72+CK4PkSNGmnm9TF1+TapIxERGaTZs2djxIgRaNq0KWJiYmBtbQ0ASEtLw4cffihxOiKqTqZyGV56ri4A4CBHuyfSe4IoipU+l2bLli2YPn06FixYgObNm1dHrmqTk5MDpVKJ7Oxs2NraSh2HqkHImChclHtBLMxBTGQHeHu6Sh2JiOiJ2DfpHtuUqOqWxidjypZzqF/HEoGN6kodR2+4Ki0wvGMjyGWC1FHIAFVXv1SlQfH69++PgoIC+Pn5wczMDBYWFmWev3v3rk7CEVXFukkD0WLMKmiUbug1YxNOzh0CmaxKN3QgIjJKy5cvh4ODA7p27QoAGDt2LP73v/+hadOmWL16NTw9PSVOSETVqb23IwAg5W4BUo4WSJxGv3g72+D15i5SxyDSqlJBP2fOHF3nINIZW2tLfNfjBUTuTEeWVX0s3n0GQ17zkzoWEZHB+O9//4sFCxYAAA4fPoyoqCjMnj0bW7duxejRo7FhwwaJExJRdXrO0RqL3muFS+m5UkfRG/FXM3Ek6S7iLt9mQU96pVIFvUajwYwZM7B582YUFxejc+fOmDRpUrkj9ERSe6fzS4hJ3ICdacC8+HR0fbEx3O0tpY5FRGQQUlNTtYPfbdq0Ce+88w4++OADBAYGokOHDtKGI6IaEdLMBSHNWLg+1NTNFkeS7uLgZY4rQPqlUuchT506FZ9//jlsbGxQr149zJ07F8OHD6+ubETPJGrkW2jlaY/colKMWXeat14hIqoga2tr7eC3u3btQnBwMADA3NwchYWFUkYjIpLES8/VhalcQMrdAly7ky91HCKtShX0K1aswPz587Fz505s2rQJW7ZswapVq6DRaKorH1GVyWUCZvb0g4WpHIeT7mDE9+ukjkREZBCCg4MRERGBiIgIXLp0CV26dAEAnD17ttLXz8fFxaFbt25wc3ODIAjYtGnTU9dZtWoV/Pz8YGlpCVdXVwwaNIh31yEiSVkpTNCyvj0AII5H6UmPVKqgT0lJ0d6LFgCCgoIgCAJu3ryp82BEuuDlYIW3n3vw/23XTbH3+FlpAxERGYCoqCgEBATg9u3bWL9+PerWfTDK9YkTJ9CnT59KbSs/Px9+fn6Iioqq0PLx8fHo378/Bg8ejLNnz2LdunU4evQo3n///Uq/DiIiXXqlsQMA4ODl2xInIfpbpa6hLy0thbm5eZl5pqamKCkp0WkoIl36qn8wNgxfgEK7Bhiy9BDO+HrD3MxU6lhERHrLzs4O8+bNKzd/ypQpSExMrNS2QkNDyxwMeJrDhw/Dy8sLH330EQCgQYMGGDJkCL799ttK7ZeISNdeaeyIGbsu4dCVOyhVa2Ai512USHqVKuhFUcSAAQOgUCi08+7fv4+hQ4fCyspKO4+j35I+kcvlWDnydby9+ARKbNzQ778r8evkQVLHIiIyGLm5uVi9ejWio6Nx/PhxqNXqattXQEAAPv/8c2zfvh2hoaG4desW1q1bpz3t/3GKiopQVFSkfZyTk1NtGYnIODWvp4SdpSmyCkpw6noWWnnWkToSGaB/908KhaJMfV1ZlfpaKTw8HE5OTlAqldqpX79+cHNzKzOPSN+0btYIPZ97MNbDsfy6+HXvMYkTERHpv7i4OISHh8PV1RUzZsxAx44dceTIkWrdZ2BgIFatWoWwsDCYmZnBxcUFdnZ2Tz1lf9q0aWU+i3h4eFRrTiIyPnKZgMCGD067j7vE6+ipajw8PMr0V9OmTXum7QmiKBrV0N85OTlQKpXIzs6Gra2t1HGoBmk0GrQcMQ9Ztg0hy81Awre9obTmreyISHr61Delp6dj2bJliI6ORk5ODnr16oWFCxfi1KlTaNq06TNtWxAEbNy4Ed27d3/sMufOnUNQUBBGjx6NkJAQpKWlYcyYMXjxxRcRHR392PUedYTew8NDL9qUiGqPNUdTMG7DGbTytMf6Ye2kjkMG5GFfn5qaWqZfqtEj9ESGTCaT4ZfPekAsyILGxhkj/7dL6khERHqlW7du8PHxwenTpzFnzhzcvHkTP/zwQ41mmDZtGtq1a4cxY8bA19cXISEhmD9/PpYsWYK0tLTHrqdQKGBra1tmIiLStZf/GhhPlZqF7EKOI0aV9+++6lmKeYAFPRkZb083fPiiHQDgwG1THL7K2yARET20Y8cODB48GFOmTEHXrl0hl8trPENBQQFMTMoO8fMwh5GdVEhEesjd3hLPOVpBrRH5OZL0Agt6MjqfvdcV77bxgCgCn647hdz7/HaViAgADh48iNzcXLRq1Qpt27bFvHnzkJn5bNeJ5uXlQaVSQaVSAQCSk5OhUqmQkpICABg/fjz69++vXb5bt25Yv349FixYgKSkJMTHx+Ojjz5CmzZt4Obm9kxZiIh04ZVGf92+7gpvX0fSY0FPRmlC16bwqGOBG1mF6Ddzo9RxiIj0wksvvYTFixcjLS0NQ4YMwZo1a+Dm5gaNRoOYmBjk5uZWepvHjx+Hv78//P39AQCRkZHw9/fHxIkTAQBpaWna4h4ABgwYgFmzZmHevHlo3rw5evbsCR8fH95Bh4j0xiuNHQEABy5zYDySHgfFI6MVc+pPRPx8BoIgw4cvyDC2b8Xvk0xEpEv63DddvHgR0dHRWLlyJbKyshAcHIzffvtN6lhPpc9tSkSGLa+oFC2m7EKpRkTsmA7wrGv19JXI6FVXv8SCnoxa50/m4appA4iFOdg7tjMa1nOSOhIRGSFD6JvUajW2bNmCJUuWsKAnIqPXa9FhHE2+C0szOUzlPOm5OnnVtcRPEW1hY24qdZRnwoJeR9jB0z/dy8lFy3FrIdq6om5BCo7PHQpBEKSORURGhn2T7rFNiag6rfrjGiZsTJQ6htGY27sF3mxRT+oYz6S6+iWTpy9CVHvZ29pg2ptN8NnuTNyxrI/PF23EtKFvSx2LiIiIiPRY37ae6OjjhMIStdRRarUlB5Ox6o8UxF66bfAFfXVhQU9Gr3dIIFbHLsQpjQd+vliK3peuwc/bU+pYRERERKTH3OwspI5Q63V5wRWr/kjBgcuZEEWRZ9I+Ai/4IAKw+ssBkGenQlBYYdCifdBojOpKFCIiIiIivdPayx4WpnLczi3ChfTK32nFGLCgJwJgaWGOqH5tIRPVuGPqiJ/+uCZ1JCIiIiIio6YwkeOl5+oAAOIu3ZY4jX5iQU/0l9cD/THxjRcAAP/dfh7JmfkSJyIiIiIiMm6vejsCAOIus6B/FBb0RP/QP8ALgY3q4n6JBv3n70ZRcYnUkYiIiIiIjNYrjR8U9MeS76GguFTiNPqHBT3RP8hkAqa/4wuZugipBSYI/2aV1JGIiIiIiIxWQ0cr1LOzQLFagz+S70odR++woCf6l3r2lujqWgAAOJxrh98OnJQ4ERERERGRcRIEAa96OwDgdfSPwoKe6BHmjuoDm+wkCHJTjF6jQl7BfakjEREREREZpVf/Ou2eBX15LOiJHkEmk2Htp29CU5ANtY0zwv6zUupIRERERERGqV1DB8gE4OrtfNzIKpQ6jl5hQU/0GE0b1sfAF8wBAIlqF6zYHi9xIiIiIiIi46O0NEULDzsAwAEepS+DBT3RE0x5/2045l2FIMgwacdV5BYWSx2JiIiIiMjo8PZ1j8aCnugpfp3QGybFORCt6mLazotSxyEiIiIiMjoPC/qDlzNRqtZInEZ/mEgdgEjfebo6YcXwYPRZ/Ad+/iMFwU2d0dHHSepYRERERERGw7eeErbmJsi5X4qp28+jrpWZ1JEqpTA/r1q2y4KeqALaNXTAoMAGWBKfjMg1J7Dpgxfh6eogdSwiIiIiIqNgIpfhlcaO2HYmDUvj/5Q6TqVpigqqZbss6IkqaOzrPth89BLuFJqixzfrcGzuMKkjEREREREZjbGv+8DRRoGiUrXUUSqtqCAPs6thuyzoiSrI3FSOj1+yx5ex2bhtUR8Tf9yMryLelDoWEREREZFR8KxrhclvNJM6RpXk5ORUS0HPQfGIKqF/1/ZoLlwHACw/ex9nk65LnIiIiIiIiIwVC3qiSlo7aSBkWdchKKzRe9ZWaDQcZZOIiIiIiGoeC3qiSrK2tMDcd1tCLC1GrrUHPvp+ndSRiIiIiIjICLGgJ6qCbu1fRDvrTADAllQT/HEuSeJERERERERkbFjQE1XR8s8HwDwnFYKpOabHpkOtEaWORERERERERoQFPVEVmZmZYtfX/WCtMMGJa/ew+ACP0hMRERERUc1hQU/0DOrXtcbE/2sKAJi56yJiVZclTkRERERERMaCBT3RM+rZ2h1t6lmgRC1i8OI4FNwvljoSEREREREZARb0RM9IEAR8EdIAmsIclNq44N2pK6WORERERERERoAFPZEO+Ho3QF8fOQBAVeSI1buOSJyIiIiIiIhqOxb0RDry32E9USfnCgSZHJ//dgF3c/KljkRERERERLUYC3oiHREEAb+O7wUx/y5Ea0f0/M8qqSMREREREVEtxoKeSIee83DFiDZ2AICrsnqIWr9X4kRERERERFRbsaAn0rEx73WDe+FVAMCyc6XILiiROBEREREREdVGLOiJqsGW/0TAq64lbueXYPKWs1LHISIiIiKiWogFPVE1sLexxKywFpAJwMaEG1gTf1HqSEREREREVMuwoCeqJi3r2yOiXX0AwGfrEnDh2k2JExERERERUW3Cgp6oGo3o2BBC9g0I5jYIm/EbNBqN1JGIiIiIiKiWYEFPVI1srS0xs6cvxNISZFt54NP566WOREREREREtYTkBX1UVBS8vLxgbm6Otm3b4ujRo49ddv/+/RAEodyUnp5eg4mJKuftTi+hjUUGAGB9koBj55IkTkRERERERLWBpAX92rVrERkZiUmTJuHkyZPw8/NDSEgIbt269cT1Ll68iLS0NO3k5ORUQ4mJqmblFwNhmnUNgpkF+kfthlrNU++JiIiIiOjZSFrQz5o1C++//z4GDhyIpk2bYuHChbC0tMSSJUueuJ6TkxNcXFy0k0wm+YkGRE9krjDDokGvQFNciEKbeoiY8bPUkYiIiIiIyMBJVgkXFxfjxIkTCAoK+juMTIagoCAcPnz4ieu2aNECrq6uCA4ORnx8/BOXLSoqQk5OTpmJSAqd2ryAYMdcAMCBLDtczsiVOBERERERERkyyQr6zMxMqNVqODs7l5nv7Oz82GviXV1dsXDhQqxfvx7r16+Hh4cHOnTogJMnTz52P9OmTYNSqdROHh4eOn0dRJXxv7H9EdjADqWigNG/qFDCU++JiIiIiKiKDOpcdR8fHwwZMgStWrVCu3btsGTJErRr1w6zZ89+7Drjx49Hdna2dkpNTa3BxERlyeVyzHq3FZQWpki8kYPZuy5IHYmIiIiIiAyUZAW9g4MD5HI5MjIyyszPyMiAi4tLhbfTpk0bXLly5bHPKxQK2NralpmIpORsa47/dG8OAIjadwXr9x2TOBERERERERkiyQp6MzMztGrVCnv27NHO02g02LNnDwICAiq8HZVKBVdX1+qISFRt/s/XFXY5VyDI5Biz/iyy8wqkjkRERERERAZG0lPuIyMjsXjxYixfvhznz5/HsGHDkJ+fj4EDBwJ4cLp8//79tcvPmTMHmzdvxpUrV5CYmIhRo0Zh7969GD58uFQvgahKBEHAus96QpN/DxprR/T6zyqpIxERERERkYExkXLnYWFhuH37NiZOnIj09HS0aNECO3fu1A6Ul5aWhpSUFO3yxcXF+OSTT3Djxg1YWlrC19cXu3fvRseOHaV6CURV1tjTDUNbWuN/F4GLcMPizbF4/832UsciIiIiIiIDIYiiKEodoibl5ORAqVQiOzub19OTXggYORdpVo2Agns4MqkbXOoqpY5ERDWMfZPusU2JiEifVFe/ZFCj3BPVRuu/7Acx9zZgaY8e/10jdRwiIiIiIjIQLOiJJObmVBfjOrpBFDW4buqOnYlpUkciItKZuLg4dOvWDW5ubhAEAZs2bXri8gMGDIAgCOWmZs2a1VBiIiIiw8GCnkgPDOvxGnq9UAcAMGFjIjLziiRORESkG/n5+fDz80NUVFSFlp87dy7S0tK0U2pqKurUqYOePXtWc1IiIiLDI+mgeET0t//0fgmnf4jHxYxcfL7hDBa91wqCIEgdi4jomYSGhiI0NLTCyyuVSiiVf48lsmnTJty7d097BxwiIiL6G4/QE+kJhYkcs8L8YCITsOtcBsYt2ih1JCIiyUVHRyMoKAienp5PXK6oqAg5OTllJiIiIn3z776qqOjZzsxlQU+kR5q5KeErSwUArLmkhupSylPWICKqvW7evIkdO3YgIiLiqctOmzZNe3RfqVTCw8OjBhISERFVjoeHR5n+atq0ac+0PRb0RHrmpwnhkGelQDCzRN/vd0Ct1kgdiYhIEsuXL4ednR26d+/+1GXHjx+P7Oxs7ZSamloDCYmIiConNTW1TH81fvz4Z9oeC3oiPWNpYY75/V+CpuQ+8q3d8eGctVJHIiKqcaIoYsmSJXjvvfdgZmb21OUVCgVsbW3LTERERPrm332VQqF4pu2xoCfSQyHt/NHRLgsAsDPNHLEnL0iciIioZsXGxuLKlSsYPHiw1FGIiIj0Fgt6Ij0VPT4c5lnJEEzM8P6SgyguKZU6EhFRpeXl5UGlUkGlUgEAkpOToVKpkJLyYIyQ8ePHo3///uXWi46ORtu2bdG8efMazUtERGRIWNAT6SkTExMs/zAYmqJ8FFu74rttp6WORERUacePH4e/vz/8iCdCVwAAIABJREFU/f0BAJGRkfD398fEiRMBAGlpadri/qHs7GysX7+eR+eJiIiegvehJ9JjbX198FG7dMw7kYelR9Px5ovZaF5P+fQViYj0RIcOHSCK4mOfX7ZsWbl5SqUSBQUF1ZiKiIioduAReiI990mPV9HlBReUakSMXqvC/RK11JGIiIiIiEgPsKAn0nOCIOA/3V+Ag7UZLt/KQ8+vlksdiYiIiIiI9AALeiIDUMfKDL28HgyKd7rYEcu2HZQ4ERERERERSY0FPZGBGNsvFE55VyAIMkzZmYTb93KkjkRERERERBJiQU9kQH6d8C7EvEyIVnXRY+pqqeMQEREREZGEWNATGZD6rk745GUniKIG10zcMWttjNSRiIiIiIhIIizoiQzMR71D0aA4GQAw99Bt/JmWKXEiIiIiIiKSAgt6IgO0ftIAIDsNgoUSkT8ffeI9nomIiIiIqHZiQU9kgOra2eLbt5tBLgAnb4v45Xiq1JGIiIiIiKiGsaAnMlBhwQGIfM0HADDpt7O4mJ4rcSIiIiIiIqpJLOiJDNiw9g3xqrcj7pdo0GP2DuTdL5E6EhERERER1RAW9EQGTCYTMKGzO9R5d5ErWKHf7N+kjkRERERERDWEBT2RgfPxdEN3p7sQNWqoss2xaNcpqSMREREREVENYEFPVAvM/nw4bK7FAQC+iUnGhbRsiRMREREREVF1Y0FPVAuYmJhg/dcfoDjlFES5KfpG7UVBcanUsYiIiIiIqBqxoCeqJXy8vTHqJXuU5t7BnVIzDFsaz/vTExERERHVYizoiWqRT0cMQYObeyGqSxGbnIcFsVeljkRERERERNWEBT1RLSKTybD6h/8g1PnBPem/+/0i9l28JXEqIiIiIiKqDizoiWoZd3d3LPy0H95tUx+iCHz0cwIuZ+RKHYuIiIiIiHSMBT1RLTXljWbwd7dFblEpBiw9itu5RVJHIiIiIiIiHWJBT1RLmcoF3Fg7ESX3buJG1n28v+I47peopY5FREREREQ6woKeqJYSBAFfjo3ErV+nQF2YC1VqFkatUUGt4cj3RERERES1AQt6olqsW7du6N89BLc3/hdQl2Dn2XRM3JzI29kREREREdUCLOiJarnZs2fDRZaD21tmAKKIVX+kYO6ey1LHIiIiIiKiZ8SCnqiWs7W1xfLly1F46RDuxCwAAMzZfRnL4pMlTkZERERERM+CBT2REWjfvj1GjRqFvITtUFzeAwCYvOUcVh9NkTgZERERERFVlYnUAYioZkydOhUajQaffz4Oi49lYvGBZHy+8QzM5DK808pd6nhERERERFRJLOiJjISFhQXmzJkDAPi8iyOKSjVYcfgaxvx6CiKAHizqiYiIiIgMCk+5JzJCgiCgSeFZvNm8LjQi8Om6U1j1xzWpYxERERERUSWwoCcyQlOmTEGfd3vjwk+TMaCdFwBgwsZERB/kQHlERERERIaCBT2REerVqxfMzc2x6/ff4XzjAIa2bwgA+HrrOXyz4wI0Gt6nnoiIiIhI37GgJzJCTZo0wTfffAMA+PTTT/BOIznGhPgAABbGXsUn606huFQjZUQiIiIiInoKFvRERmrkyJHo2LEjCgoKMGDAAAx9tQG+6+ELuUzAxoQbCF9yFPfyi6WOSUREREREj8GCnshIyWQyLF26FDY2Njh06BBmzJiBnq09EB3eGpZmchxOuoPu8+Nx5Vau1FGJiIiIiOgRWNATGTFPT0/MnTsXADBp0iSkpaWhg48TNnzYDu72Frh2pwBvRR3C7nMZEiclIiIiIqJ/Y0FPZOQGDBiAoUOHYsuWLXB1dQUAPO9ii83DA9HGqw5yi0oRseI4pu04jxI1r6snIiIiItIXLOiJjJwgCFiwYAGCg4PLzK9rrcBPEW0xMNALALAoNgl9Fh/BzaxCCVISEREREdG/saAnojL+/PNPqFQqAICZiQyTujXD/L4tYa0wwbE/7yFkThw2JdyAKPLWdkREREREUmJBT0RasbGx8PX1Rc+ePZGfn6+d3+UFV2wd+TL8POyQe78Uo9aqMOLnBNzJK5IwLRERERGRcWNBT0Rafn5+UCqVuHLlCsaOHVvmOS8HK6wfGoDIYG/IZQK2nUlD51mx+OV4Ko/WExERERFJgAU9EWnZ2dlh6dKlAID58+dj165dZZ43kcvwUefG2PhhOzzvYoOsghKM/fU03l18BBfTeXs7IiIiIqKaxIKeiMoICgrCiBEjAACDBg3CvXv3yi3j626HLSNfxvjQ52FuKsORpLsInRuHCRvPIJOn4RMRERER1QgW9ERUzrfffgtvb2/cuHEDI0eOfOQypnIZhrRviJjR7RHa3AUaEVj1Rwo6frcfc3ZfQs79khpOTURERERkXFjQE1E5lpaWWLFiBWQyGVatWoVt27Y9dlmPOpZY0K8V1n7wEl6op0RuUSnm7L6Ml7/Zix/2XEZ2IQt7IiIiIqLqIIhGNppVTk4OlEolsrOzYWtrK3UcIr02adIkFBYW4quvvoK5uflTl9doRGxPTMPc3Zdx+VYeAMDKTI5eL3pgUGADeNSxrO7IRAaJfZPusU2JiEifVFe/xIKeiHROrRGx7UwaovZewcWMB4PlyQSgcxNn9GlTH696O0IuEyROSaQ/2DfpHtuUiIj0SXX1SyY62xIR1WqlpaVITExEixYtnrqsXCbgDT83dPN1RdzlTPx4IAkHLmci5lwGYs5lwE1pjrdbuuONFm7wdrapgfRERERERLUPC3oieqrMzEz83//9H86ePYvTp0+jQYMGFVpPEAS093ZEe29HXMrIxZqjqdiQcB03s+9j3r4rmLfvCp53sUHXF1wR1NQZz7vYQBB45J6IiIiIqCJ4yj0RPZVarUbHjh1x4MABvPrqq9i3bx9ksqqNqXm/RI2YcxnYrLqJ2Eu3UKL++09QPTsLdHreCS83dsBLDepCaWmqq5dApNfYN+ke25SIiPQJr6HXEXbwRFWTlJQEX19f5OfnY+bMmYiMjHzmbWYVFGNnYjpizmXg4JVMFJVqtM/JBKB5PSVa1reHr7sSvu5KPOdgDRmvvadaiH2T7rFNiYhIn7Cg1xF28ERVt3jxYnzwwQdQKBQ4efIkmjZtqrNtFxarEX8lE3GXbyP+Siau3s4vt4y1wgTN69nC190Ozesp8ZyDFRo4WMFKwauHyLCxb9I9tikREemT6uqXJL8PfVRUFLy8vGBubo62bdvi6NGjT1x+//79aNmyJRQKBRo1aoRly5bVTFAiQkREBEJDQ1FUVIT33nsPJSW6u8e8hZkcQU2d8dWbzbHnkw44Mr4z5vZugYGBXmjtaQ9zUxnyikpxJOku/heXhI9WJ+D/fjiIZpN+R5upuxG26DDGrT+NRbFXsTMxDSdT7uFGViGK/3HUn4hqXlxcHLp16wY3NzcIgoBNmzY9dZ2ioiJMmDABnp6eUCgU8PLywpIlS2ogLRERkWGR9LDW2rVrERkZiYULF6Jt27aYM2cOQkJCcPHiRTg5OZVbPjk5GV27dsXQoUOxatUq7NmzBxEREXB1dUVISIgEr4DIuAiCgOjoaDRv3hwnT57E1KlTMXny5GrZl4vSHG+2qIc3W9QDAJSqNbhyOw+nU7Nx6noWzqfl4M87BbibX4xbuUW4lVuEP5LvPnJbdazM4GSjgJOtORytFbCzNIXS4h/TPx5bmZnAwlQOczMZzOQyDtJH9Izy8/Ph5+eHQYMG4e23367QOr169UJGRgaio6PRqFEjpKWlQaPhl3NERET/Jukp923btsWLL76IefPmAQA0Gg08PDwwcuRIjBs3rtzyn332GbZt24bExETtvN69eyMrKws7d+6s0D55Ch7Rs/vll18wffp0rFixQqen3VdFdkEJku/k48/MfCRlPvg39V4BbuUU4Vbu/TKD7lWWIAAWpvIHBb6pHOamMliY/f1YYSKHiUyAiVz4618ZTOUC5DIBJrKH///nPAGCIEAQAJkgQMBf/woPviyRCXgwT/bguX8u++C5fy37j3lVf41PXvlpm37S6sJT1n5a7mfZt6HJz8tFjwCfWt83CYKAjRs3onv37o9dZufOnejduzeSkpJQp06dKu+L/T0REemTWncf+uLiYpw4cQLjx4/XzpPJZAgKCsLhw4cfuc7hw4cRFBRUZl5ISAhGjRr12P0UFRWhqKhI+zg7OxvAgwYloqp5/fXXERwcDLlcLvnvkgDgOaUMzylt0Klh2Xvai6KIewUlyMy7j1u5xbiVcx+ZeUXIvV+KnMIS5NwvQU5hKbLvl2gfFxZroNb8/SVA3n0gr4ZfExkfTVEBgAc/s8but99+Q+vWrTF9+nSsXLkSVlZWeOONN/D111/DwsLiseuxvyciIn32sD962D89pFAooFAoqrxdyQr6zMxMqNVqODs7l5nv7OyMCxcuPHKd9PT0Ry6fk5ODwsLCR3b006ZNw5QpU8rN9/DweIb0REREupebmwulUil1DEklJSXh4MGDMDc3x8aNG5GZmYkPP/wQd+7cwdKlSx+7Hvt7IiIyBPXr1y/zeNKkSc90CWutHxp6/PjxZW6vlZWVBU9PT6SkpBj9hyZ9kZOTAw8PD6SmpvK0SD3B90S/8P3QP7p+T0RRRG5uLtzc3HSQzrBpNBoIgoBVq1Zp++lZs2ahR48emD9//mOP0rO/r178O6R7bFPdYnvqHttUt7Kzs1G/fn1cu3YNdnZ22vnPcnQekLCgd3BwgFwuR0ZGRpn5GRkZcHFxeeQ6Li4uj1ze1tb2sR38405hUCqV/MHUM7a2tnxP9AzfE/3C90P/6PI9YdH5gKurK+rVq1emPZo0aQJRFHH9+nU0btz4keuxv68Z/Duke2xT3WJ76h7bVLfs7Oxqx23rzMzM0KpVK+zZs0c7T6PRYM+ePQgICHjkOgEBAWWWB4CYmJjHLk9ERESGJTAwEDdv3kRe3t+jV1y6dAkymQzu7u4SJiMiItI/kt6HPjIyEosXL8by5ctx/vx5DBs2DPn5+Rg4cCCAB6fP9e/fX7v80KFDkZSUhLFjx+LChQuYP38+fvnlF4wePVqql0BERERPkJeXB5VKBZVKBeDBLWhVKhVSUlIAlO/r+/Tpg7p162LgwIE4d+4c4uLiMGbMGAwaNOiJg+IREREZI/nk6rqJdAU0b94cdnZ2mDp1KmbMmAEAWLVqFXx8fAAAP/30E65du4YBAwYAAOzt7REYGIj58+fj66+/xtmzZzFz5kz06NGjUvuVy+Xo0KEDTExq/RACBoPvif7he6Jf+H7oH74nFXPo0CEEBARg0aJFAIDff/8dixYtwr1799C9e/dyfb2ZmRm6dOmCjRs34ssvv8SOHTvw9ttv47vvvoOpqWml9s33SLfYnrrHNtUttqfusU11qzraU9L70BMRERERERFR1Uh6yj0RERERERERVQ0LeiIiIiIiIiIDxIKeiIiIiIiIyACxoCciIiIiIiIyQEZX0EdFRcHLywvm5uZo27Ytjh49KnUkoxUXF4du3brBzc0NgiBg06ZNUkcyatOmTcOLL74IGxsbODk5oXv37rh48aLUsYzaggUL4OvrC1tbW9ja2iIgIAA7duyQOhb95ZtvvoEgCBg1apTUUYxSZfvz/fv3o2XLllAoFGjUqBGWLVtWM0ENSGXaNC0tDX369IG3tzdkMhl/Dx6jMm26YcMGBAcHw9HRUfs3//fff6/BtPqvMu158OBBBAYGom7durCwsMDzzz+P2bNn12Baw1DV2ig+Ph4mJiZo0aJFNSc0LJVpz/3790MQhHJTenp6pfZpVAX92rVrERkZiUmTJuHkyZPw8/NDSEgIbt26JXU0o5Sfnw8/Pz9ERUVJHYUAxMbGYvjw4Thy5AhiYmJQUlKC1157Dfn5+VJHM1ru7u745ptvcOLECRw/fhydOnXCm2++ibNnz0odzegdO3YMixYtgq+vr9RRjFJl+/Pk5GR07doVHTt2hEqlwqhRoxAREcFi6R8q26ZFRUVwdHTEF198AT8/vxpOaxgq26ZxcXEIDg7G9u3bceLECXTs2BHdunVDQkJCDSfXT5VtTysrK4wYMQJxcXE4f/48vvjiC3zxxRf43//+V8PJ9VdVa6OsrCz0798fnTt3rqGkhqGq7Xnx4kWkpaVpJycnp8rtWDQibdq0EYcPH659rFarRTc3N3HatGkSpiJRFEUA4saNG6WOQf9w69YtEYAYGxsrdRT6B3t7e/HHH3+UOoZRy83NFRs3bizGxMSI7du3Fz/++GOpIxmdyvbnY8eOFZs1a1ZmXlhYmBgSElKtOQ3Js3xG4u/Bo+nic2fTpk3FKVOmVEc8g6OL9nzrrbfEfv36VUc8g1TVNg0LCxO/+OILcdKkSaKfn191xzQYlW3Pffv2iQDEe/fuPdN+jeYIfXFxMU6cOIGgoCDtPJlMhqCgIBw+fFjCZET6KTs7GwBQp04diZMQAKjVaqxZswb5+fkICAiQOo5RGz58OLp27VqmP6GaU5X+/PDhw+Xer5CQEPb/f+FnJN3TRZtqNBrk5uayH4Zu2jMhIQGHDh1C+/btqyumQalqmy5duhRJSUmYNGlSTcQ0GM/yM9qiRQu4uroiODgY8fHxld63SaXXMFCZmZlQq9VwdnYuM9/Z2RkXLlyQKBWRftJoNBg1ahQCAwPRvHlzqeMYtTNnziAgIAD379+HtbU1Nm7ciKZNm0ody2itWbMGJ0+exLFjx6SOYrSq0p+np6c/cvmcnBwUFhbCwsKi2vIaAn5G0j1dtOmMGTOQl5eHXr16VUdEg/Is7enu7o7bt2+jtLQUkydPRkRERHVGNRhVadPLly9j3LhxOHDgAExMjKaMrJCqtKerqysWLlyI1q1bo6ioCD/++CM6dOiAP/74Ay1btqzwvvlOEFE5w4cPR2JiIg4ePCh1FKPn4+MDlUqF7Oxs/PrrrwgPD0dsbCyLegmkpqbi448/RkxMDMzNzaWOQ0S12M8//4wpU6Zg8+bNlb+elso4cOAA8vLycOTIEYwbNw6NGjXCu+++K3Usg6NWq9GnTx9MmTIF3t7eUsepFXx8fODj46N93K5dO1y9ehWzZ8/GypUrK7wdoynoHRwcIJfLkZGRUWZ+RkYGXFxcJEpFpH9GjBiBrVu3Ii4uDu7u7lLHMXpmZmZo1KgRAKBVq1Y4duwY5s6di0WLFkmczPicOHECt27dKvOtuVqtRlxcHObNm4eioiLI5XIJExqHqvTnLi4uj1ze1tbW6I/OA/yMVB2epU3XrFmDiIgIrFu3jpf2/OVZ2rNBgwYAgBdeeAEZGRmYPHkyC3pUvk1zc3Nx/PhxJCQkYMSIEQAenNEpiiJMTEywa9cudOrUqUay6yNd/R1t06ZNpQ+oGc019GZmZmjVqhX27NmjnafRaLBnzx5ej0oEQBRFjBgxAhs3bsTevXu1HSDpF41Gg6KiIqljGKXOnTvjzJkzUKlU2ql169bo27cvVCoVi/kaUpX+PCAgoMzyABATE8P+/y/8jKR7VW3T1atXY+DAgVi9ejW6du1aE1ENgq5+RtmH/q2ybWpra1uuDxw6dKj2TMK2bdvWZHy9o6ufUZVKBVdX18rt/JmG1DMwa9asERUKhbhs2TLx3Llz4gcffCDa2dmJ6enpUkczSrm5uWJCQoKYkJAgAhBnzZolJiQkiNeuXZM6mlEaNmyYqFQqxf3794tpaWnaqaCgQOpoRmvcuHFibGysmJycLJ4+fVocN26cKAiCuGvXLqmj0V84urc0ntafjxs3Tnzvvfe0yyclJYmWlpbimDFjxPPnz4tRUVGiXC4Xd+7cKdVL0DuVbVNRFLV9eKtWrcQ+ffqICQkJ4tmzZ6WIr5cq26arVq0STUxMxKioqDL9cFZWllQvQa9Utj3nzZsn/vbbb+KlS5fES5cuiT/++KNoY2MjTpgwQaqXoHeq8nv/TxzlvqzKtufs2bPFTZs2iZcvXxbPnDkjfvzxx6JMJhN3795dqf0aVUEviqL4ww8/iPXr1xfNzMzENm3aiEeOHJE6ktF6eKuGf0/h4eFSRzNKj3ovAIhLly6VOprRGjRokOjp6SmamZmJjo6OYufOnVnM6xkW9NJ5Un8eHh4utm/fvszy+/btE1u0aCGamZmJzz33HP+2PUJl2/RRfYanp2fNhtZzlWnT9u3b83PRU1SmPb///nuxWbNmoqWlpWhrayv6+/uL8+fPF9VqtQTJ9Vdlf+//iQV9eZVpz2+//VZs2LChaG5uLtapU0fs0KGDuHfv3krvUxBFUazcMX0iIiIiIiIikprRXENPREREREREVJuwoCciIiIiIiIyQCzoiYiIiIiIiAwQC3oiIiIiIiIiA8SCnoiIiIiIiMgAsaAnIiIiIiIiMkAs6ImIiIiIiIgMEAt6IiIiIiIiIgPEgp6I9MbkyZPRokULqWMQERHpzP79+yEIArKysqSOUqsJgoBNmzZJHYOoxrGgJzJAhw8fhlwuR9euXWt0vxUtuCdPngxBECAIAkxMTODl5YXRo0cjLy/viet9+umn2LNnj67iEhERVauHfd3jpsmTJ6Ndu3ZIS0uDUqmUJJ8uityH/frrr79e7rnvvvsOgiCgQ4cOkmR7KC0tDaGhoQCAP//8E4IgQKVS6Wz7RPrKROoARFR50dHRGDlyJKKjo3Hz5k24ublJHamcZs2aYffu3SgtLUV8fDwGDRqEwsJCLFy4sNyyoihCrVbD2toa1tbWEqQlIiKqvLS0NO3/165di4kTJ+LixYvaedbW1jAzM4OLi4sU8XTK1dUV+/btw/Xr1+Hu7q6dv2TJEtSvX1+yXMXFxbWmjYmqgkfoiQxMXl4e1q5di2HDhqFr165YtmxZmefv3buHvn37wtHRERYWFmjcuDGWLl0K4EGnN2LECLi6usLc3Byenp6YNm2adt2srCxERETA0dERtra26NSpE06dOgUAWLZsGaZMmYJTp05pjzz8e9//ZGJiAhcXF7i7uyMsLAz9+vXD5s2bAfx9+uGOHTvQqlUrKBQKHDx48JFnACxZsgTNmjWDQqGAq6srRowYUaG8RERE1c3FxUU7KZVKCIJQZp61tXW5U+6XLVsGOzs7bN26FT4+PrC0tESPHj1QUFCA5cuXw8vLC/b29vjoo4+gVqu1+yoqKsKnn36KevXqwcrKCm3btsX+/fsfm83LywsA8NZbb0EQBO1jAFiwYAEaNmwIMzMz+Pj4YOXKlU99rU5OTnjttdewfPly7bxDhw4hMzOz3BmDx44dQ3BwMBwcHKBUKtG+fXucPHnyqdkGDBiA7t27l9nWqFGjyhz979ChA0aMGIFRo0bBwcEBISEhAMoe8W/QoAEAwN/fX3v2QFxcHExNTZGenl5u+6+88spTXz+RvmJBT2RgfvnlFzz//PPw8fFBv379sGTJEoiiqH3+yy+/xLlz57Bjxw6cP38eCxYsgIODAwDg+++/x2+//YZffvkFFy9exKpVq8p08D179sStW7ewY8cOnDhxAi1btkTnzp1x9+5dhIWF4ZNPPkGzZs2QlpaGtLQ0hIWFVTi3ubk5iouLy8wbN24cvvnmG5w/fx6+vr7l1lmwYAGGDx+ODz74AImJidi2bRu8vb0rlJeIiEhfFRQU4Pvvv8eaNWuwc+dO7N+/H2+99Ra2b9+O7du3Y+XKlVi0aBF+/fVX7TojRozA4cOHsWbNGpw+fRo9e/bE66+/jsuXLz9yH8eOHQMALF26FGlpadrHGzduxMcff4xPPvkEiYmJGDJkCAYOHIh9+/Y9NfegQYPKfJm/ZMkS9O3bF2ZmZmWWy83NRXh4OA4ePIgjR46gcePG6NKlC3Jzc5+YraKWL18OMzMzxMfHP/LMv6NHjwIAdu/ejbS0NGzYsAGvvvoqnnvuuTJfXpSUlGDVqlUYNGhQpfZPpFdEIjIo7dq1E+fMmSOKoiiWlJSIDg4O4r59+7TPd+vWTRw4cOAj1x05cqTYqVMnUaPRlHvuwIEDoq2trXj//v0y8xs2bCguWrRIFEVRnDRpkujn5/fUjP9e7vjx46KDg4PYo0cPURRFcd++fSIAcdOmTU9cz83NTZwwYcIj91GRvERERDVl6dKlolKpLDf/YZ9379497XIAxCtXrmiXGTJkiGhpaSnm5uZq54WEhIhDhgwRRVEUr127JsrlcvHGjRtltt25c2dx/Pjxj80EQNy4cWOZee3atRPff//9MvN69uwpdunS5bHbedg/FxcXi05OTmJsbKyYl5cn2tjYiKdOnRI//vhjsX379o9dX61WizY2NuKWLVuemC08PFx88803y8z797bbt28v+vv7P/G1JicniwDEhISEMst8++23YpMmTbSP169fL1pbW4t5eXmPzU6k73iEnsiAXLx4EUePHsW7774L4MFp7WFhYYiOjtYuM2zYMKxZswYtWrTA2LFjcejQIe1zAwYMgEqlgo+PDz766CPs2rVL+9ypU6eQl5eHunXraq9lt7a2RnJyMq5evVrprGfOnIG1tTUsLCzQpk0bBAQEYN68eWWWad269WPXv3XrFm7evInOnTs/8nld5yUiIqoplpaWaNiwofaxs7MzvLy8yowj4+zsjFu3bgF40Keq1Wp4e3uX6fNiY2Mr3eedP38egYGBZeYFBgbi/PnzT13X1NQU/fr1w9KlS7Fu3Tp4e3s/8gy7jIwMvP/++2jcuDGUSiVsbW2Rl5eHlJSUSmV9nFatWlVpvQEDBuDKlSs4cuQIgAeXP/Tq1QtWVlY6yUUkBQ6KR2RAoqOjUVpaWmYQPFEUoVAoMG/ePCiVSoSGhuLatWvYvn07YmJi0LlzZwwfPhwzZsxAy5YtkZycjB07dmD37t3o1asXgoKC8OuvvyIvLw+urq6PvB7Pzs6u0ll9fHzw22+/wcTEBG5ubuVOxwPwxA7UwsLiidvXdV4iIqKaYmpqWuaxIAiPnKfRaAAeibQWAAAEnklEQVQ86PPkcjlOnDgBuVxeZrmaHkx20KBBaNu2LRITEx97qnp4eDju3LmDuXPnwtPTEwqFAgEBAeUuvfs3mUxW5jJC4MFp8f9W1QLcyckJ3bp1w9KlS9GgQQPs2LHjieMQEBkCFvREBqK0tBQrVqzAzJkz8dprr5V5rnv37li9ejWGDh0KAHB0dER4eDjCw8PxyiuvYMyYMZgxYwYAwNbWFmFhYQgLC0OPHj3w+uuv4+7du2jZsiXS09O1t5l7FDMzszID9DyJmZkZGjVqVOXXa2NjAy8vL+zZswcdO3Ys93xF8hIREdUG/v7+UKvVuHXrVqUGcDM1NS3Xbzdp0gTx8fEIDw/XzouPj0fTpk0rtM1mzZqhWbNmOH36NPr06fPIZeLj4zF//nx06dIFAJCamorMzMynZnN0dERiYmKZeSqVqtyXHU/z8CDCoz6zRERE4N1334W7uzsaNmxY7mwFIkPDgp7IQGzduhX37t3D4MGDy93L9p133kF0dDSGDh2KiRMnolWrVmjWrBmKioqwdetWNGnSBAAwa9YsuLq6wt/fHzKZDOvWrYOLiwvs7OwQFBSEgIAAdO/eHdOnT4e3tzdu3ryJbdu24a233kLr1q3h5eWF5ORkqFQquLu7w8bGBgqFotpe8+TJkzF06FA4OTkhNDQUubm5iI+Px8iRIyuUl4iIqDbw9vZG37590b9/f8ycORP+/v64ffs29uzZA19f33KjzD/08IvxwMBAKBQK2NvbY8yYMejVqxf8/f0RFBSELVu2YMOGDdi9e3eF8+zduxclJSWPPSOucePGWLlyJVq3bo2cnByMGTOm3Jl3j8rWqVMnfPfdd1ixYgUCAgLw008/ITExEf7+/hVvLDw4Em9hYYGdO3fC3d0d5ubm2s9OISEhsLW1/f/27t/V2DCO4/jHjyxmlLIyKINYlSIGxWjxK4vJopTBIpMoDMogG2aDxX+AiX+A2SCZdZ7h9Kjneeqc8/zqnFvv13rXdX/u8dN9fa9L7XZbrVbrt9YFviJm6AGDmEwmisViv5R56bXQ73Y77fd72Ww2NRoNBQIBRSIRWSwWLRYLSa9/vTudjkKhkMLhsI7Ho1arlcxms0wmk1arlSKRiEqlkrxer7LZrE6nk1wu1+M9yWRS0WhUDodD8/n8v35zoVBQv9/XaDSS3+9XKpV6nOb7kbwAADyL6XSqfD6vWq0mn8+nTCaj7Xb75h3wvV5P6/VaHo/nUYozmYwGg4G63a78fr/G47Gm0+kPV8O9x263vzneNplMdLlcFAwGlcvlVK1W5XQ6382WSCTUbDZVr9cVDod1u92Uz+c/nOs7q9Wq4XCo8Xgst9utdDr9eGY2m1UsFnW/3/9obeCrMb38PKgCAAAAAE+qXC7rfD5ruVx+dhTgr7HlHgAAAMDTu16vOhwOms1mlHk8DQo9AAAAgKeXTqe12WxUqVQUj8c/Ow7wT7DlHgAAAAAAA+JQPAAAAAAADIhCDwAAAACAAVHoAQAAAAAwIAo9AAAAAAAGRKEHAAAAAMCAKPQAAAAAABgQhR4AAAAAAAOi0AMAAAAAYEDfAMpykGuvKnA8AAAAAElFTkSuQmCC", "text/plain": [ "PyPlot.Figure(PyObject )" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, axes = subplots(1, 2, figsize=(12, 4))\n", "\n", "\n", "axes[1][:plot]([0, strike], [strike, 0], \"k--\")\n", "axes[1][:plot](ps, v[1:n-1])\n", "axes[1][:set_xlim](0, strike*2)\n", "axes[1][:set_xticks](linspace(0, 4, 5))\n", "axes[1][:set_ylim](0, strike)\n", "axes[1][:set_yticks](linspace(0, 2, 5))\n", "axes[1][:set_xlabel](\"Asset Price\")\n", "axes[1][:set_ylabel](\"Premium\")\n", "axes[1][:set_title](\"Put Option Value\")\n", "\n", "\n", "axes[2][:plot](linspace(0, T, N), max_exercise_price)\n", "axes[2][:set_xlim](0, T)\n", "axes[2][:set_ylim](1.6, strike)\n", "axes[2][:set_xlabel](\"Time to Maturity\")\n", "axes[2][:set_ylabel](\"Asset Price\")\n", "axes[2][:set_title](\"Put Option Optimal Exercise Boundary\")\n", "axes[2][:tick_params](right=\"on\")" ] } ], "metadata": { "kernelspec": { "display_name": "Julia 0.5.1", "language": "julia", "name": "julia-0.5" }, "language_info": { "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", "version": "0.5.1" } }, "nbformat": 4, "nbformat_minor": 2 }