{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "\n", "How to read this lecture…\n", "\n", "- If using QuantEcon lectures for the first time on a computer, execute ] add InstantiateFromURL inside of a notebook or the REPL \n", "- For some notebooks, enable content with “Trust” on the command tab of Jupyter lab \n", "- Code should execute sequentially if run in a Jupyter notebook \n", "- Please direct feedback to [contact@quantecon.org](mailto:contact@quantecon.org\") or [discourse forum](http://discourse.quantecon.org/) " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Job Search II: Search and Separation\n", "\n", "\n", "" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Contents\n", "\n", "- [Job Search II: Search and Separation](#Job-Search-II:-Search-and-Separation) \n", " - [Overview](#Overview) \n", " - [The Model](#The-Model) \n", " - [Solving the Model using Dynamic Programming](#Solving-the-Model-using-Dynamic-Programming) \n", " - [Implementation](#Implementation) \n", " - [The Reservation Wage](#The-Reservation-Wage) \n", " - [Exercises](#Exercises) \n", " - [Solutions](#Solutions) " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Overview\n", "\n", "Previously [we looked](mccall_model.ipynb#) at the McCall job search model [[McC70]](zreferences.ipynb#mccall1970) as a way of understanding unemployment and worker decisions\n", "\n", "One unrealistic feature of the model is that every job is permanent\n", "\n", "In this lecture we extend the McCall model by introducing job separation\n", "\n", "Once separation enters the picture, the agent comes to view\n", "\n", "- the loss of a job as a capital loss, and \n", "- a spell of unemployment as an *investment* in searching for an acceptable job " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Setup" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "using InstantiateFromURL\n", "activate_github(\"QuantEcon/QuantEconLecturePackages\", tag = \"v0.9.0\") # activate the QuantEcon environment\n", "\n", "using LinearAlgebra, Statistics, Compat # load common packages" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## The Model\n", "\n", "The model concerns the life of an infinitely lived worker and\n", "\n", "- the opportunities he or she (let’s say he to save one character) has to work at different wages \n", "- exogenous events that destroy his current job \n", "- his decision making process while unemployed \n", "\n", "\n", "The worker can be in one of two states: employed or unemployed\n", "\n", "He wants to maximize\n", "\n", "\n", "\n", "\n", "\n", "
\n", "$$\n", "{\\mathbb E} \\sum_{t=0}^\\infty \\beta^t u(Y_t)\n", "$$\n", "\n", "(1)\n", "
\n", "\n", "The only difference from the [baseline model](mccall_model.ipynb#) is that\n", "we’ve added some flexibility over preferences by introducing a utility function $ u $\n", "\n", "It satisfies $ u'> 0 $ and $ u'' < 0 $" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Timing and Decisions\n", "\n", "Here’s what happens at the start of a given period in our model with search and separation\n", "\n", "If currently *employed*, the worker consumes his wage $ w $, receiving utility $ u(w) $\n", "\n", "If currently *unemployed*, he\n", "\n", "- receives and consumes unemployment compensation $ c $ \n", "- receives an offer to start work *next period* at a wage $ w' $ drawn from a known distribution $ p_1, \\ldots, p_n $ \n", "\n", "\n", "He can either accept or reject the offer\n", "\n", "If he accepts the offer, he enters next period employed with wage $ w' $\n", "\n", "If he rejects the offer, he enters next period unemployed\n", "\n", "When employed, the agent faces a constant probability $ \\alpha $ of becoming unemployed at the end of the period\n", "\n", "(Note: we do not allow for job search while employed—this topic is taken\n", "up in a [later lecture](jv.ipynb#))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Solving the Model using Dynamic Programming\n", "\n", "Let\n", "\n", "- $ V(w) $ be the total lifetime value accruing to a worker who enters the current period *employed* with wage $ w $ \n", "- $ U $ be the total lifetime value accruing to a worker who is *unemployed* this period \n", "\n", "\n", "Here *value* means the value of the objective function [(1)](#equation-objective) when the worker makes optimal decisions at all future points in time\n", "\n", "Suppose for now that the worker can calculate the function $ V $ and the constant $ U $ and use them in his decision making\n", "\n", "Then $ V $ and $ U $ should satisfy\n", "\n", "\n", "\n", "\n", "\n", "
\n", "$$\n", "V(w) = u(w) + \\beta [(1-\\alpha)V(w) + \\alpha U ]\n", "$$\n", "\n", "(2)\n", "
\n", "\n", "and\n", "\n", "\n", "\n", "\n", "\n", "
\n", "$$\n", "U = u(c) + \\beta \\sum_i \\max \\left\\{ U, V(w_i) \\right\\} p_i\n", "$$\n", "\n", "(3)\n", "
\n", "\n", "Let’s interpret these two equations in light of the fact that today’s tomorrow is tomorrow’s today\n", "\n", "- The left hand sides of equations [(2)](#equation-bell1-mccall) and [(3)](#equation-bell2-mccall) are the values of a worker in a particular situation *today* \n", "- The right hand sides of the equations are the discounted (by $ \\beta $) expected values of the possible situations that worker can be in *tomorrow* \n", "- But *tomorrow* the worker can be in only one of the situations whose values *today* are on the left sides of our two equations \n", "\n", "\n", "Equation [(3)](#equation-bell2-mccall) incorporates the fact that a currently unemployed worker will maximize his own welfare\n", "\n", "In particular, if his next period wage offer is $ w' $, he will choose to remain unemployed unless $ U < V(w') $\n", "\n", "Equations [(2)](#equation-bell1-mccall) and [(3)](#equation-bell2-mccall) are the Bellman equations\n", "for this model\n", "\n", "Equations [(2)](#equation-bell1-mccall) and [(3)](#equation-bell2-mccall) provide enough information to solve out for both $ V $ and $ U $\n", "\n", "Before discussing this, however, let’s make a small extension to the model" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Stochastic Offers\n", "\n", "Let’s suppose now that unemployed workers don’t always receive job offers\n", "\n", "Instead, let’s suppose that unemployed workers only receive an offer with probability $ \\gamma $\n", "\n", "If our worker does receive an offer, the wage offer is drawn from $ p $ as before\n", "\n", "He either accepts or rejects the offer\n", "\n", "Otherwise the model is the same\n", "\n", "With some thought, you will be able to convince yourself that $ V $ and $ U $ should now satisfy\n", "\n", "\n", "\n", "\n", "\n", "
\n", "$$\n", "V(w) = u(w) + \\beta [(1-\\alpha)V(w) + \\alpha U ]\n", "$$\n", "\n", "(4)\n", "
\n", "\n", "and\n", "\n", "\n", "\n", "\n", "\n", "
\n", "$$\n", "U = u(c) +\n", " \\beta (1 - \\gamma) U\n", " + \\beta \\gamma \\sum_i \\max \\left\\{ U, V(w_i) \\right\\} p_i\n", "$$\n", "\n", "(5)\n", "
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Solving the Bellman Equations\n", "\n", "We’ll use the same iterative approach to solving the Bellman equations that we\n", "adopted in the [first job search lecture](mccall_model.ipynb#)\n", "\n", "Here this amounts to\n", "\n", "1. make guesses for $ U $ and $ V $ \n", "1. plug these guesses into the right hand sides of [(4)](#equation-bell01-mccall) and [(5)](#equation-bell02-mccall) \n", "1. update the left hand sides from this rule and then repeat \n", "\n", "\n", "In other words, we are iterating using the rules\n", "\n", "\n", "\n", "\n", "\n", "
\n", "$$\n", "V_{n+1} (w_i) = u(w_i) + \\beta [(1-\\alpha)V_n (w_i) + \\alpha U_n ]\n", "$$\n", "\n", "(6)\n", "
\n", "\n", "and\n", "\n", "\n", "\n", "\n", "\n", "
\n", "$$\n", "U_{n+1} = u(c) +\n", " \\beta (1 - \\gamma) U_n\n", " + \\beta \\gamma \\sum_i \\max \\{ U_n, V_n(w_i) \\} p_i\n", "$$\n", "\n", "(7)\n", "
\n", "\n", "starting from some initial conditions $ U_0, V_0 $\n", "\n", "Formally, we can define a “Bellman operator” T which maps:\n", "\n", "\n", "\n", "\n", "\n", "
\n", "$$\n", "TV(\\cdot) = u(\\cdot) + \\beta (1-\\alpha)V(\\cdot) + \\alpha U\n", "$$\n", "\n", "(8)\n", "
\n", "\n", "In which case we are searching for a fixed point\n", "\n", "\n", "\n", "\n", "\n", "
\n", "$$\n", "TV^{*} = V^*\n", "$$\n", "\n", "(9)\n", "
\n", "\n", "As before, the system always converges to the true solutions—in this case,\n", "the $ V $ and $ U $ that solve [(4)](#equation-bell01-mccall) and [(5)](#equation-bell02-mccall)\n", "\n", "A proof can be obtained via the Banach contraction mapping theorem" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Implementation\n", "\n", "Let’s implement this iterative process" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "solve_mccall_model (generic function with 1 method)" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "using Distributions, LinearAlgebra, Compat, Expectations, Parameters, NLsolve, Plots\n", "\n", "function solve_mccall_model(mcm; U_iv = 1.0, V_iv = ones(length(mcm.w)), tol = 1e-5, iter = 2_000)\n", " # α, β, σ, c, γ, w = mcm.α, mcm.β, mcm.σ, mcm.c, mcm.γ, mcm.w\n", " @unpack α, β, σ, c, γ, w, dist, u = mcm\n", "\n", " # parameter validation\n", " @assert c > 0.0\n", " @assert minimum(w) > 0.0 # perhaps not strictly necessary, but useful here\n", "\n", " # necessary objects\n", " u_w = u.(w, σ)\n", " u_c = u(c, σ)\n", " E = expectation(dist) # expectation operator for wage distribution\n", "\n", " # Bellman operator T. Fixed point is x* s.t. T(x*) = x*\n", " function T(x)\n", " V = x[1:end-1]\n", " U = x[end]\n", " [u_w + β * ((1 - α) * V .+ α * U); u_c + β * (1 - γ) * U + β * γ * E * max.(U, V)]\n", " end\n", "\n", " # value function iteration\n", " x_iv = [V_iv; U_iv] # initial x val\n", " xstar = fixedpoint(T, x_iv, iterations = iter, xtol = tol).zero\n", " V = xstar[1:end-1]\n", " U = xstar[end]\n", "\n", " # compute the reservation wage\n", " wbarindex = searchsortedfirst(V .- U, 0.0)\n", " if wbarindex >= length(w) # if this is true, you never want to accept\n", " w_bar = Inf\n", " else\n", " w_bar = w[wbarindex] # otherwise, return the number\n", " end\n", "\n", " return (V = V, U = U, w_bar = w_bar) # return a NamedTuple, so we can select values by name\n", "end" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The approach is to iterate until successive iterates are closer together than some small tolerance level\n", "\n", "We then return the current iterate as an approximate solution\n", "\n", "Let’s plot the approximate solutions $ U $ and $ V $ to see what they look like\n", "\n", "We’ll use the default parameterizations found in the code above" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "#5 (generic function with 2 methods)" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# a default utility function\n", "u(c, σ) = (c^(1 - σ) - 1) / (1 - σ)\n", "\n", "# model constructor\n", "McCallModel = @with_kw (α = 0.2,\n", " β = 0.98, # discount rate\n", " γ = 0.7,\n", " c = 6.0, # unemployment compensation\n", " σ = 2.0,\n", " u = u, # utility function\n", " w = range(10, 20, length = 60), # wage values\n", " dist = BetaBinomial(59, 600, 400)) # distribution over wage values" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deWBU5bk/8Oc9J5nsJCSEbCSBhACyy6KWRG2KFbcQsHC1FRV79V7k0iq2Xq1bRKiQVJFaUYsUqP1RbMWtFotiAohS17AYSFhky76vM8ks531+fwzkxrCFSTJnZs7389fMO+fkPHNycr55z/IewcwEAABgVIreBQAAAOgJQQgAAIaGIAQAAENDEAIAgKEhCAEAwNAQhAAAYGgIQgAAMDQEIQAAGBqCEAAADA1BCAAAhubpQXj06FGz2dzDiR0OB0aMczNN06SUeldhLFJKTdP0rsJYmNnhcOhdheHY7Xb3LMjTg/Duu+/eu3dvDyduaGjATtnNWlpabDab3lUYi8VisVgseldhLDabraWlRe8qjEVK2dDQ4J5leXoQAgAA9CsEIQAAGFpvg7CoqCgkJKTzbXp6ujhjwYIF3SbeuXPnxIkTw8LCJk6c+MknnzgbGxsbs7KyIiMjZ86c2djY2Mt6AAAALkmvgrC5uXn+/PmdpyuYuaSkpKysrLW1tbW1ddWqVd2mnzdv3uOPP97Q0PDYY4/NmzfP2Zibm5ucnFxZWZmUlJSXl9ebegAAAC6Vn8tzMvP8+fMfffTRuXPnOluqq6ttNlt2dvahQ4emT5++Zs2awMDArrMMGDCgubm5ra2ttbU1NDTU2fjOO++89957AQEBixYtys7OXr58eddZpJRlZWVHjhzptvSwsLDo6OhujZqm4Wo6N9PO0LsQA9E0jZmxzt0J27n7Oa+O7vk6VxRFCOHasoTL9xusWLGirq7uueeeE+L0D9m7d+9DDz20cuXKpKSkxYsX22y2TZs2dZ3l66+/njp1qvP1V199NWXKFCIKDQ2tra0NCgpqb2+PiYnpdmnWuHHjamtrAwICui39lltuycnJ6dZYV1c3cOBAVVVd+0bggqampqCgoLN/QdB/zGYzM3f+KwluYLPZLBZLRESE3oUYiJSyvr7+7A7P+YSHh7u8I3IxCLdv375kyZJt27b5+/t3BmFXlZWVY8aM6Xbx6/Tp06dMmfLggw++8MILhYWFH3/8MRGFhITU19cHBgZaLJbo6Ohudw2mp6fn5eWlp6f3pKqampqoqCgEoTs1NjYGBQV16/pDv2pra2PmsLAwvQsxEKvVajabIyMj9S7EQKSUtbW1MTExbliWi4dG8/Pzd+7caTKZnG+FELt27QoODu7o6Jg2bRoRmUyms8P5iy++2LhxY2xs7COPPDJ06FBnY3x8fGlpaVpaWnl5eUJCgovfAwAAwCUuXiyzbNkyPoOImDkjI8NsNs+ePbu4uNhmsy1dunTWrFnOiXfs2OF8MX78+D/96U9tbW2vv/76hAkTnI1ZWVnr1q1j5nXr1mVnZ/f2CwEAgDfTmL6plauK5CmLm27w68vFZGRk5OTkZGVlJSQkNDY25ubmOtszMzOdL9atW/fBBx/ExcVt3rx57dq1zsannnpq//79iYmJBw4ceOKJJ/qwHgAA8CJHW3hNsbxnp/Z0oSyo4M/qXL+c85L0wWI6TxAKIRYuXLhw4cLzTTBq1KjPPvus26cRERFbtmzpfRkAAOCNjrfyzkreVcU17afDIjFUXB0jxvi7aaxRN+UtAABAV9XttLNS7qiQpWcukYwOpGvilGtiRcoAIaWsrXXTQxQQhAAA4D7NNv60indWcknz6cOF4SaRHkPXximXRbh8K2CvIAgBAKDf2STtreftFfx5DTskE5FJoStilMx4MSmK/BQ9AvAMBCEAAPQXybSnTu6opM9rZIdGROSniKnR4ofxypXRIsAz7vpGEAIAQN872sI7KvmTSm60MhEJQaMHimtjxdVxSpi/3sV9H4IQAAD6TF0H767m/Ao+1nL6UpchIeLqWJEZr8QF61vaeSEIAQCgtywO2l3N+eXyQOPpG+bCTXRNrJIZL9LC9Tz/1xMIQgAAcJHzFOD2Sv53NdskEVGASldGi8x45fJBQvX0BDwNQQgA4N1uvvnmiooKNy/UIalDow6NtDNDpvgrIkilQD/6juivvfjJI0eOfOONN/qkyB5CEAIAeLcDBw68/PLLcXFxehfSBw4fPvzss8+6eaEIQgAArzd69OjOR/p4NUVx00Db31uo+xcJAADgORCEAABgaAhCAAAwNAQhAAAYGoIQAAAMDUEIAACGhiAEAABDQxACAIChIQgBAKCPLV++fObMmV1brrvuupdfflmvei4MQQgAAH3stttu++ijj1paWpxvGxsbd+3adeutt+pb1fkgCAEAoI+lpKSMHz/+/fffd77dsmXLVVddFRsbq29V54OxRgEAfE2rnWvb3bSsMJOIDjxH++2337558+Y77riDiN599925c+e6qaBLhyAEAPBukru3fFlLq77V3LP06QnKg2PPcXBx7ty5TzzxRGtrq7+//0cfffTiiy+6px4XIAgBALySXdIXNfLDMjY7un8U6kfDwtz0VNxzdgeJKDExcfLkyVu2bAkLC5swYUJ8fLx76nEBghAAwMuUm3lrGeeXy1b7uSe4crC4crDq3qLO4bbbbnvzzTcjIyPnzJmjdy0XgiAEAPAODslf1NK/SuX+BnY+Fn74ADFjiNjvr3dl5zFnzpzf/OY3wcHBOTk5etdyIQhCAABPV9dBH5bJj8q4wcpEFKjSNXHihiFKWrggIjcdA710sbGxV155pcViGTJkiN61XAiCEADAQzHRvnr+sIz/XS01JiJKCKbrEpQZiSLM32Pj73s+/vhjvUu4OAQhAIDHabXzx+X0r1JZaWEi8lcoPUbcmKiMjfSO/PMuCEIAAA9ypJk/KJWfVLJNEhHFBNENQ5QfD1HCTXpX5rsQhAAA+uu8F2JvPRORIJoYJWYMEdNiFAWdwH6GIAQA0FN1O/2rVG4rly02IqIwf/pxgnJjoogNRgC6CYIQAEAHzLSvgf95Sn5Z47wVgtLCxc1JytUxZFIRgW6FIAQAcKt2jQrK5T9PcZn59IUwGbEiK+n0vRA+Qwhx+m7Hc731KAhCAAA3qbTQR+Vya6lssxMRRQaIHyWIWxKVqPOMUgbugSAEAOhfzFRYJ98vpcI66ewUjR0obk4SP4hRcBDUEyAIAQD6S4dG+eXy/ZOy3EJEFKDSD+OUW5KUoWH9vGCpSbu1n5dxmlD8hL9339uBIAQA6HtVFv6wnD8s5VY7E1FkgJgxRGQlK2FuGRfUenR/68533LEkosARl4deO9s9y+onCEIAgL60v4H/cfL/rgUdM1BkJStXDRZuPQqq+omAIDcty89Tx/zuMQQhAEAfsEv6pFL+4xQfazl9Leg1sWLmUDWlv4+CnktA6riA1HE6LLiL8PDwxsbGgQMHElF9fX1ERIS+9VwAghAAoFeabfyvUv6glButp4+C3pgobkw0+qBoM2bM+M1vfrNixQpN0x5++OFbbrlF74rOC0EIAOCicjNvOcUflUurRkSUMkDcOET8KEExKXpX5gFWr179wAMPpKamKopyyy23vPjii3pXdF4IQgCAS+O8HeLdk7yvgZlJEF05WGQnK+PwaIguBg0atHHjRr2r6BEEIQBAT9kl7aiU756Qp9qIiIJUum6IkpWkxAXrXRn0AoIQAODi2hzi42Pyn6dOnwgcFEhZycqMBBHiJQ/IhQtAEAIAXEiVhd/+TnxYGiBVSUQpA0R2snJNLPnh8Ui+AkEIAHBuxU389nH5RQ1rkiSLH0Qrs4aK8TgR6HMQhAAA38NEX9Tw28dlcRMTkUkVP4ql6VEd44a46xZ1cC8EIQDAaTaNCyr43ZNcbmYiCvMXNyWKW5KVIHaYzR76CCHoPQQhAABZHPRxuXz7hKzvICIaHEg3Jio3Jp6+FsbqpvGrQR8IQgAwtLoOeveE9mEZd2hERKkDxK1DRXosHpBkIAhCADCo0jZ+67jcWckOJiFo0iBx61BlQhQC0HAQhABgON+18Hsn5c5KlkyCaGq0+GmqkhbuxRG4du1a5/DW3q6ystL9C0UQAoCBHGzkzcflV7WnHxBxbZy4LUVJCPHiCCSijRs3vv322+3t7XoX0jfuuusuNy8RQQgAvk8yfVbNbx6Tx1uZiEL96aZEZWayjzwgIj09PT09Xe8qvBiCEAB8mYNpe4V86/jpOyKiAik7WbkhUQlS9a4MPAaCEAB8k11Sfrn82zFZ10FEFBNEM5OVG4YIE64Hhe9DEAKAr2nXaFuZfOs4N1iZiJJDxa3DlGvjBBIQzglBCAC+o8VG752UW0rZbGciGhku5qYoV0QLgQiE80MQAoAvaLLR28e1f5Wevi9+QpT4jxQFA2RDTyAIAcC7NdnovZPy/ZPSqhERTYwS89KUkd58UyC4GYIQALxVTTu/e4I/LGebxkL4wn3xoAsEIQB4n0oL/e07bUcla0yKoGtixX+kKsmhiEBwBYIQALxJhYX/9p3cWckak58ipsfRXO8fGgb0pfRy/qKiopCQkM636enp4owFCxZ0m1ic5aKzAAA4VbfT6oNy4adaQQULITLjxep05cFxKlIQeqlXPcLm5ub58+dbLBbnW2YuKSkpKysLDw8nIj+/7j+8tbW183VeXp7NZrvoLAAAVRZ+6wRvK5POXmBmHP00VY0L1rss8BWu9wiZef78+Y8++mhnS3V1tc1my87OjouLmzdvXktLS7dZQs84fvz47t27ly1bdtFZAMDIysz8u33af+3StpZKQXRjovLHDOWhcUhB6Euu98Byc3NTU1PnzJnT2VJVVTV16tSVK1cmJSUtXrz4gQce2LRp09kz2my2e++9d/369X5+fhedpamp6brrrlPV7sMC/uQnP8nLy+vWWFdX53A4zp4Y+k9zc3NgYGBAQIDehRiI2Wxm5s4jMb6qsl28VxHw73o/yeSv0DWDbDMTHJEmyS1U7fZ/mG02m8Visdvt7l6wgUkp6+vrez59RESEyzsiwcwuzLZ9+/YlS5Zs27bN399fiHP8kMrKyjFjxjQ0NJw97/LlyysqKv7whz/0ZJZp06bl5ORcddVV3Sb29/cPDAzs1lhTUxMVFYUgdKfGxsagoKCzfxfQf9ra2pg5LCxM70L6S6WF/36cdlRIjclP0PVDxJxhYlCgnicCrVar2WyOjIzUsQajkVLW1tbGxMT0cHpFcf0Ap4s9wvz8/J07d5pMpx9hIoTYtWtXcHBwR0fHtGnTiMhkMp0znDVNe/XVV/Pz851vCwsLLzyLECI0NNR5BvGilDNc+1LgAqxz91MUhZl9cp3XtPObx3lbmfOKUOXaWPppqhIXrP+1MNjOdeG2de7iMpYtW8ZnEBEzZ2RkmM3m2bNnFxcX22y2pUuXzpo1yznxjh07OmcsKChITEwcPny48+35ZgEAQ6ntoJcOyP/+VG4tlUKIHyeIP16tPjRO9YQUBJ/Xl1dpZmRk5OTkZGVlNTc333TTTZ0HPzMzMzuPnW7YsOH666+/6CwAYBANVt58nLeWSrskVdD0eHF7qhKL/AM3cvEcodukp6fn5eX18OHLOEfofjhH6H4+c46wxUZvHZdbSqVVI0H0w3jl9lQR75ERiHOE7nep5wh7A/ftAYC7dWj0z5Ny8wk220+PETpvuJIywBMjEIwAQQgA7tOh0fsn5TsnuNXORDQ1WrljuEhFBIKuEIQA4A4Opn+dkn8/JptsRETjI8W8NOWyCEQg6A9BCAD9SzLtqOSNR7SaDiKiURHizjQ8Mhc8CIIQAPrR3nped0geb2UiSgyhnw1X0mMUgRAET4IgBIB+UdzEfz4sDzQyEUUH0n+kKNcPURREIHgeBCEA9LHvWnjDYbm3noloYIC4LUXMSFT8EIHgqRCEANBnKi38/47yrirJTKH+dOtQJStZCcSdveDZEIQA0AeabfTGd3JrqXQwmRTKGqrMGaaE+utdFkAPIAgBoFe+d3c8UUascs9IZTDGGgLvgSAEABc5JH9Yxm98d/rWwCsHi7tHqIkhepcFcIkQhABwyZhpdw3/+bBWaSEiuixC3DMSd8eDt0IQAsClKW7idYdkSRMTUWKouDtNuXIwIhC8GIIQAHqqwsJ/OcKfVUtmGmCi21KUm5MUFSEIXg5BCAAX12rnt0/weyekXVKASlnJytxhSjD2H+ATsCEDwIXYNP7HKd58nM12VgRdP0TcMVyNDNC7LIC+gyAEgHNjpl1V8s9HuKadiWhKtLhnhJoUqndZAH0NQQgA51DSxGtL5KFmJqKUAeLnI5QJUTgZCL4JQQgA31PXQa8f0XZUMjNFBtBPUzFYNvg4BCEAnNah0dvH5Vsn2KYxrogB48A2DgAkmT4skxuPcrONBdGP4sWdacqgQHQDwRAQhABGt7+BXyvRTrQSEY0dKP5zlDJ8ACIQDARBCGBclRZef1j+u5qJKCaIfj5SnRaDCATDQRACGFHX04GBKs0eqswZJkwYJAYMCUEIYCzMtL2SNxyWjVYWgjJilZ+PVKLx1CQwMAQhgIEUN/GaYnm0hYnosghx3yglLRy9QDA6BCGAITRYecNh6bw7cFAg3TNSvTpGCIQgAIIQwOc5JH9QyhuPSouD/ATdmKTcNUIJVPUuC8BjIAgBfNm+el5Top1qIyKaGi3+a5QSG4xuIMD3IAgBfFOlhV8/wp9WSSKKDxb3jVKmRCMCAc4BQQjga2yS3jwm3z7BNo2D/ej2VCUrSfhhtFCA80AQAviUL2r4tRKtup2EoOkJyt1pYmAAIhDgQhCEAD6iysJrSuRXtUxEw8LE/aOVyyIQgQAXhyAE8Ho2SZuPnR4mJsSPfjZcuSUJD04C6CkEIYB3+7KWXyuRVRYWgjLjxc9HqhEmvWsC8CoIQgBvVdNBfzyofYljoQC9gyAE8D4OyZuPyTeOSatGIf7ijuHipkQFI2YDuAZBCOBlipvF2iOmCqskoqnRYuFoPEEXoFcQhABeo9XOrx/hLSdUZkocQAsuUycNQgQC9BaCEMALOJ+d9KdDWouN/BWRPUSbd5mKxwcC9AkEIYCnO9HKqw/KkiYmokmDxJ1JjphAiRQE6CsIQgDPZZO06ah89yQ7JEcG0H2j1IxY0dbGzHpXBuBDEIQAHmpPPb9yUFZaWBDdkqTcmaYE4+8VoB/gDwvA4zgvivmwTDLT0DBaNEYdiefIA/QbBCGAB+l6UYxJoZ+kKP+RqvghBAH6E4IQwFOUW+ilIq2okYlocrRy/2VKTJDeNQEYAIIQQH8OyW+f4L8dY5vGESa67zL1mlh0AwHcBEEIoLMjzfyHA/J4KwtBP04QPx+phvrrXROAkSAIAXRjk7TpO/n2cSmZYoPF/4xWJkahIwjgbghCAH0UNfBLB2W5mVVBWUnKXSOUQFXvmgAMCUEI4G5mB204LDvvjvjlGDUNd0cA6AdBCOBW/67mV4plo5VNqrg9Rdw6DI9PAtAZghDATZpt/Gqx/LSKiWhspFg0WkkIQQYC6A9BCOAOn1bJV4u52cYBKt2eovwkRUEGAngIBCFA/2qy0SsHtd3VpzuCvxyjxgXrXRMAdIEgBOhHn1bJlw/KVjsF+9E9I5QZQxSBniCAh0EQAvSL+g566aD2dS0T0ZRo8T+jlUGByEAAT4QgBOh7BRW8plgzOyjMX9w7SvlRPCIQwHMhCAH6UqOVVx+UX9QwEV05WPzPaGVgAFIQwKMhCAH6zKdV8pVi2WKjEH8xP03ckKjoXREAXByCEKAPNNv4lYPys2omokmDxC/G4IwggNdAEAL01le18qUD3GBlXBoK4I0QhACua7PTq8XazkomoolR4pdj1ehAvWsCgEuEIARw0Z56/n2RVt9BgSrdM0K5MREdQQCvhCAEuGTO5wi+dUwy0chw8dB4JT4YGQjgrXp7VVtRUVFISEjn2/T0dHHGggULuk0szkJEjY2NWVlZkZGRM2fObGxs7GU9AP3tcDP/8jPH5mNSVcRPU5W8K1WkIIBX61WPsLm5ef78+RaLxfmWmUtKSsrKysLDw4nIz6/7D29tbe18nZeXZ7PZiCg3Nzc5OXnz5s2/+tWv8vLyli9f3puSAPqPxvTOcbnxO3ZISgqlh8YpqQMQgQBez/UeITPPnz//0Ucf7Wyprq622WzZ2dlxcXHz5s1raWnpNkvoGcePH9+9e/eyZcuI6J133lm0aFFAQMCiRYvefvttl+sB6FelbfzwF9qfj0iNOStJWfUDP6QggG8QzOzanCtWrKirq3vuueeEOP1D9u7d+9BDD61cuTIpKWnx4sU2m23Tpk1nz2iz2a6++ur169ePHj2aiEJDQ2tra4OCgtrb22NiYrrF5/jx46WUYWFh3X7Iddddt3jx4m6NdXV1AwcOVFXVtW8ELmhqagoKCgoICNC7kH7ETB9WqptO+ttZRAfI+9PsowZIHesxm83MHBoaqmMNRmO1Wi0Wy8CBA/UuxECklPX19dHR0T2cPjQ01GQyubYsFw+Nbt++fevWrdu2bevaOHHixIKCAufrFStWjBkz5pzzPv/881dccYUzBYmImZ0nC5lZ07RuEwcFBd1xxx1jx47t1j548OCu5yadLBZLSEgIgtCdbDZbUFBQYKDP3jTQaOXVh9TCOiaFro/jn49QglSdLzFjZmY+e/uH/uM80YN17k5Syvb29p6v87NPxvWci3Pm5+fv3LmzM36FELt27QoODu7o6Jg2bRoRmUymc/YSNE179dVX8/PzO1vi4+NLS0vT0tLKy8sTEhK61+fnN3ny5PT09J5U5VwogtCdAs7Qu5B+UVjHvy+SDVYOM4mFo5Vr4jxiyDS73c7MvrrOPZbD4cA6dycp5flypM+5+Ie9bNkyPoOImDkjI8NsNs+ePbu4uNhmsy1dunTWrFnOiXfs2NE5Y0FBQWJi4vDhwztbsrKy1q1bx8zr1q3Lzs52/asA9B2bpDXF8ulCrcHKE6PE6nQ/D0lBAOhzfXmQJyMjIycnJysrq7m5+aabbvrDH/7gbM/MzOw8E7lhw4brr7++61xPPfXUHXfckZiYOGnSpL/85S99WA+Aa4408/P7tXIL+Sv0s+HKT1IUXBUD4MP6IAg7Q04IsXDhwoULF55vAiLauHFjt08jIiK2bNnS+zIAeo+Z3j8l1x+SDqbEEPr1eDUFl4YC+DqMLANwWl0Hr/xWftvAQlD2UOWu4cKkIgUBfB+CEICI6PMafrFIa7VTZIB4cJxyeRQiEMAoEIRgdDaNNxzmf5ZKZro8Siweh2fKAxgLghAMrbSNf7dfHm9lXBcDYFgIQjCuggp++aBm1WhIiHh4vILrYgCMCUEIRmS28+qDclcVE1FmvFg4Wg3EMAwARoUgBMM50MjP7Zd1HRzqT4vGqOkx6AgCGBqCEAyEid48Jv96VGpMoweKX49Xo312kFQA6CkEIRhFs41fKOJvaqUQlJWk/Hyk8MOgaQCAIASD2N/Az++XDVYON4nFY8XkaGQgAJyGIAQfJ5ne+E7+7ZiUTOMixa/Hq5F4hAAAdIEgBF/WZKPn92t761kR9NNU5fZUBfcJAkA3CELwWd/UyheKuNnGkQH08AR17EBkIACcA4IQfJBk2nhUvnlMMtHkQWLxOCXchBQEgHNDEIKvabLR7/Zp+xtYETRvuDJ3mCIQggBwfghC8ClFDfy7/VqDlcJN4tfjlYl4iAQAXAyCEHzE6WfqHmaHpLGR4uHxSiQeIgEAPYAgBF9gtvOLB3h3ddeb5ZGCANAjCELwekeaecU+WdPOYf700Dh1SjQiEAAuAYIQvNu/SuVrJdIuaWSEeGQCxg4FgEuGIARvZZP08kGZXy6JKCtZ+flIxQ9dQQC4dAhC8Eq1HbR8r3akmU2quP8y5boEZCAAuAhBCN7n61pe+a1stXNcMD02URkahhQEANchCMGbMNNbx+XrRyQTTY0WvxqvhmATBoDewV4EvIbFQS98q31ew0LQnGHKXWm4QwIA+gCCELzDsRZ+dq9W3U4DTPTweBVDxgBAX0EQghf4pIpfLNKsGqWFi99MxD0SANCXEITg0STTX47Kt45LZsqMVxaNUUx4tjwA9CkEIXiuVjvn7ZN761kVNC9NmZOCDASAvocgBA91opV+u1dWWXiAiR6ZoI6PxElBAOgXCELwRLuq+MUirUOjlDB6/HJ1cBBSEAD6C4IQPEvXk4I/jFN+MRYnBQGgfyEIwYOYHZS3TyusYz9F/OcocUsSMhAA+h2CEDxFuYWWFmrlZg430aMTlLE4KQgAboEgBI+wp55z92pmBw0Noycn+Q3GnYIA4C4IQtDf1lL5arHUmKZGi1+PV4OxVQKAG2GXA3qyS3rloNxWLoWgOSkYPhQAdIAgBN202Gj5Xq2okU2q+MUY8cM4XBoDADpAEII+jrfysj2ypp2jAunxiUpaOLqCAKAPBCHoYHe1XPmttGp0WYR47HI1wqR3QQBgYAhCcLfNx+XrhyUTXZcgFo5W/XFAFAB0hSAE93EwvVSk5VewEPTTFOVnw5GBAKA/BCG4Saudl++V3zZwoEq/Gq9eNRgnBQHAIyAIwR0qLPzMN1q5hSID6MlJ6vABSEEA8BQIQuh3e+t5xV7N7KBhYeKpSeogjBoDAJ4EQQj966MyfqVYOiRNixEPjVMDVL0LAgD4PgQh9Bcm+lOJ9t5JJqK5KcqdwxWBA6IA4HkQhNAvrBo9/63272r2V2jRGPVH8chAAPBQCELoe612WlqoFTdxqD89NlEdhwcqAYAHQxBCH6u00NPfaBUWjgminElqYihSEAA8GoIQ+lJxEy/bo7XYaES4eHISxk4DAC+AIIQ+82mVfKGIbRpdNVj8ejwuEAUA74AghL7xj5NybYlkoqwk5b7L8FhBAPAaCELoLY3pxSK5rVwqgu4bqWQlYwRRAPAmCELoFasUz+8X+5pkoEq/Hq9eiRFEAcDbIAjBdc02euZb/5PtYmAg5UxS8XBdAPBGCEJwUZWFc76RJ9uU+FBaOlVNCEEKAoBXQhCCK44085JCrdlGKaEyZ5KIQQoCgCjPwhsAABa/SURBVNdCEMIl21vPz+7R2jWaGCUWDrOHm3CfBAB4MQQhXJr8Cn7pgHRImp6gLBotWptZ74oAAHoFQQiXYPMx+fpRyUxzUpS78DQJAPAJCELoEcn0ykFtaxkrghaOVm5IxM2CAOAjEIRwcXZJz+/XPqtmk0KPTFSviEZPEAB8B4IQLqJDo2f3aHvqOcSPnpykjhmIFAQAn4IghAtptfOSb+ShZh4YIJ6epKQMQAoCgK/p7ZmeoqKikJCQzrfp6enijAULFnSb2OFwLFy4MDo6Oj09vby8vCezgI5q2vnhz7VDzRwTRLlXqEhBAPBJveoRNjc3z58/32KxON8yc0lJSVlZWXh4OBH5+XX/4atWrWppaTl58mROTk5OTs7atWsvOgvopdRMT30t6zooOVQ8M0WNDNC7IACA/uF6j5CZ58+f/+ijj3a2VFdX22y27OzsuLi4efPmtbS0dJvlr3/968MPPxwcHPzUU0/df//9PZkFdHG4mR/90lHXwWMHitwrFKQgAPgw13tgubm5qampc+bM6WypqqqaOnXqypUrk5KSFi9e/MADD2zatKnrLCdPnty0aVNmZmZKSsr69et7MovFYsnNzY2Nje229MmTJ99+++3dGltaWvz8/FQVA530yp4G5YUSP5sUU6O0Xwy3Oyyi+fwTt7S02Gw2q9XqvvoMz2w2M7OUUu9CDMRqtVosFuxb3ElK2dLSEhgY2MPpg4OD/f39XVuWi0G4ffv2rVu3btu2rWvjxIkTCwoKnK9XrFgxZsyYbnO1tLQw84EDB1avXn3fffd9/vnnF51FVdXk5ORhw4Z1a09MTDz7O/v7+/v7+2Nj7Y3dNeLFEsXBdF28/O8RQhGmC0/vf4Z7ygM6cwYB69ydpJTYzt3sUte5orh+gNPFIMzPz9+5c6fJdHovKYTYtWtXcHBwR0fHtGnTiMhkMgUEdD+gFh0d/eCDD8bFxS1atGjVqlVEVFhYeOFZAgICbr/99vT09J5UFRQUFBwcjCB02bZyfumwxgrd1uOBY6xWa1BQUM//a4Pek1Iyc3BwsN6FGIiqqljnbialdO7S3bAsFyN02bJlfAYRMXNGRobZbJ49e3ZxcbHNZlu6dOmsWbOcE+/YscP5YsaMGRs2bLBarWvWrJkyZQoRnW8WcL9/npJ/KNIk05wU5e40DJ8GAEbRlwNlZWRk5OTkZGVlJSQkNDY25ubmOtszMzOdL5YvX15QUBATE5Ofn7927doLzAJutvmY/GOxJEH3jlTuTsPwaQBgIMLZpfNY6enpeXl5PTw0WlNTExUVhUOjl4SZ1h2W756QiqBfjFGvS7i0nmBjYyMOjbpZW1sbM4eFheldiIFYrVaz2RwZGal3IQYipaytrY2JiXHDsnDfnqFJppcOaNvK2U/QwxPUaTE4HgoAhoMgNC4H0/P75adVHKjSY5erl0chBQHAiBCEBmXVaPk++U2tDPWnnEnqqAikIAAYFILQiDo0eqZQ+7aBI0z0zBR1WBhSEACMC0FoOGY7LymUxU0cFUi/neqXgDujAMDYEITGYnZQTqE81MTRgfTbqX5xSEEAMDwEoYG02jnnG3mkmQcH0rNX+MUE6V0QAIAHQBAaRZONnvxanmjlhGD67VS/KNz4BwBARAhCg2i08uNfaaVmSgwVy/BwQQCALhCEvq+2gx7/Squ0UEoYLZ2iDrjI8yQAAIwFQejjKi30+FeO2g4aES6WTFZD8RgZAIDvQxD6sgoLP/aVVt9BoweKnElqMH7bAABnwa7RZ3Wm4LhI8dQkNRBDkQMAnAseuOObOlNwzECkIADAhSAIfVC5+f9S8OnJSEEAgAvBoVFf40zBBitSEACgRxCEPqUzBcdGihwcEQUA6AEcGvUd5RZCCgIAXCoEoY8ot9BjXzqQggAAlwqHRn1BuYV+86XWaKWJUeLJSaoJ/94AAPQYgtDrOceOQQoCALgGe03vVttBT3zlcI4d8/jlSEEAgEuGHqEXq+vgx76SNR10WYRYgjslAABcgh6Et3I+X7DKwiPDxZIpSEEAABehR+iVmm302JeOMjOlDBA5k5UgpCAAgKsQhN6nxUbOp+wODaNlU9QwPFkJAKAXEIRexuygnG+0k208JEQsRQoCAPQazhF6E7ODnvxaO9rCQ0LE8ivUCDxrHgCg1xCEXqNDoyXfaEeaOT5Y/HYqUhAAoG/g0Kh3sGm8bI8sbuLBgfTbqWpkgN4FAQD4CvQIvYBkWvmt3FfP4SZ6ZqrfoEC9CwIA8CEIQk/HRCu/1T6r5hA/emaymhCsd0EAAL7Ftw6NOuyyrUmovnNXHTOtP6wVVnCMHz16mTpUEbJN75q+j83NLG3S0aF3IQbCZjMzS6HpXYiBsM1GFovEGIZuJKUkh809y/KpIOTyI0178oUQehfSZ0rbOKGd7iIaGSHCyqhB73rOZrfbHapqVrCDcB+Hw0FENj+f+uP1cFJKTdMa/HG7kvswsxx7NSUkumFZPvW3xKq/GhbhM0H4XQsfk6wE0IQoJcJTzws6rFbFz0/1oV6455N2OxGp2Cm7k6ZJh0MNwFVq7sPMmr+bLo73qSBUhqSFT7jKN3bK75+Sa4qlIujh8WparOdGe2NjY1BQUGCgpwa1L2pra2PmsLAwvQsxEKvVajabB0ZG6l2IgUgpa2tr3bMsHNHyRB+WyddKpCB6YKya4cEpCADgAxCEHufTKvnyQclMC0arP4pHCgIA9C8EoWfZW8/P75eS6a405aZEpCAAQL9DEHqQI8382z2ag+mWJGVuCn41AADugL2tp6iw8JJCrUOjH8Yp/3UZfi8AAG6CHa5HqO+gJ7+WzTaaGi0eHKfgkCgAgNsgCPXXaqcnv9Fq2nlkuPjfCaqKGAQAcCMEoc6sGi0t1ErbeGiYeHqyGugL90ACAHgTBKGeHEzL98niJo4NFksmK6EYKgQAwO18amQZ78JEz+/XvqnlcJN4epISGYBDogAAOkCPUDdriuWnVRziL56ZrCSEIAUBAPSBINTH5mPyn6ekSRVPTFRSBiAFAQB0gyDUwSdV/PpRKYgeGqeMjUQKAgDoCUHobvsbeNW3GjPdO0pJj0EKAgDoDEHoVida+dm90i5pzjBlZjJWPgCA/rAvdp/qdnrqG2m28zWx4q40rHkAAI+A3bGbtNrp6W+0RiuPjxQPjlMFjokCAHgGBKE72CQtLdTKzDw0TDx2ueqPtQ4A4DGwS+53kum5/VpxE8cE0TOTlRCMYQAA4EmwV+53rxzU/l3NYf60ZIrfwAC9qwEAgO9Dj7B/vXVcbi1jk0JPTVITgvWuBgAAzoIg7EefVfOfD0tB9KvxyqgIXB4DAOCJEIT95Ugzv/CtxkT3jFSmxWA9AwB4KOyg+0V1Oy0p1KwazRiizB6KlQwA4LlwsUzfa7Vzztdas40mDxL3j0YKAgB4NOym+5iDacVeWW6h5FDx8HhFxZlBAADPhh5hX2KmF4u0/Q0cGUBPT1ZD8MR5AACPhx5hX9p4VG6v4CCVciapgwL1rgYAAHoAQdhn8iv4b8ekIujXE1Q8axcAwFsgCPvGtw380gFJRPePVq+IRgoCAHiN3gZhUVFRSEhI59v09HRxxoIFC7pN7HA4Fi5cGB0dnZ6eXl5e7mxsbGzMysqKjIycOXNmY2NjL+vRRaWFlu/VHJJvHarcMAQpCADgTXoVhM3NzfPnz7dYLM63zFxSUlJWVtba2tra2rpq1apu069ataqlpeXkyZPTpk3LyclxNubm5iYnJ1dWViYlJeXl5fWmHl1YHLRsj6PVTlOixd0j0MMGAPAygpldm5OZb7311jvuuGPu3LnOH1JVVZWWljZy5MhDhw5Nnz59zZo1gwcP7jrLpEmT1q9fP2HChNbW1sOHD0+ePJmIRo4c+d57740aNaqkpCQ7O/vQoUNdZ5kyZcq1116bmprabenDhw/PyMjo1lhbWxsZGamqqmvfyAWSaUWRUlgvhoXR0olaoPuW7CmampoCAwMDA3FpkPu0tbURUWhoqN6FGIjVarVYLAMHDtS7EAORUtbV1XULkQsICAhweefv+u0Tubm5qampc+bM6WypqqqaOnXqypUrk5KSFi9e/MADD2zatKnrLCdPnty0aVNmZmZKSsr69eudjeXl5cnJyUTk7Bd2W4qmaQcOHKitre3WzsxXXnllt0a73W6326WULn+pS7XuO79v6ijCxA+NtKmS7e5bsqew2+2qqrrznw9wOBzMbLfb9S7EQOxn6F2IgUgpL2md+/v7uzsIt2/fvnXr1m3btnVtnDhxYkFBgfP1ihUrxowZ022ulpYWZj5w4MDq1avvu+++zz//nIiYWQjhfKFpWrdZgoODn3zyyfT09J5UZbVaw8PD3bZT/ucpWVArgwPE01OV1HCDdomklEFBQegRupOqqswcFhamdyEGYrVa/fz8wsPD9S7EQKSUNpvNPevcxXNa+fn5O3fuNJlMzgwTQnz66aeFhYW7d+92TmAymQICuj98Lzo6+sEHH4yLi1u0aFFRUZGzMT4+vrS0lIjKy8sTEhJc/B5ut7ee15ZIIegXY8TIcFwgAwDgrVwMwmXLlvEZRMTMGRkZZrN59uzZxcXFNptt6dKls2bNck68Y8cO54sZM2Zs2LDBarWuWbNmypQpzsasrKx169Yx87p167Kzs3v7hdyi3EIr9moa020pyg/jcIEMAIAX68udeEZGRk5OTlZWVkJCQmNjY25urrM9MzPT+WL58uUFBQUxMTH5+flr1651Nj711FP79+9PTEw8cODAE0880Yf19JNWOz/zjWZ20LQY8bNUpCAAgHfrg7FGO687FUIsXLhw4cKF55sgNja222lFIoqIiNiyZUvvy3APh+QVe2WFhVMHiMXjVIFjogAAXg4dmkvzajHvb+CoQHpqkmrAmyUAAHwPgvASbDklPyyTASo9ebka2f1KIAAA8EoIwp462Hj6MtFfjlVTMaY2AICvQBD2SE07P7tXOph+MlS5JhYpCADgOxCEF2eTtHyvbLbx5VHizjSsMQAAn4Ld+kUw0++L5NEWTggRj0xUFfQGAQB8C4LwIt48Lj+plEEqPTZRDemDm00AAMCzIAgvpLCONx6VgujhCWoSxvoHAPBFCMLzKrdQ3j5NMt2ZpkzFQ+cBAHwUgvDcLA76baHD7KBpMcqcYVhLAAA+C7v4c2Ci3+3XSs2UMkAsHqdgHDUAAB+GIDyHvx6VX9dyuEk8fjnGUQMA8HEIwu6+rOW/fScVQb8arwzG42YBAHwdgvB7Kiy8cr/GRPNHqJdH4ZAoAIDvQxD+nw6Nnt0jzQ66arCYlYwUBAAwBATh/3nloHayjROCCQ8aBAAwDgThae+d5IIKDlLp8Ul+wRhBBgDAMBCERETFTbzhkCYEPTBOTQzRuxoAAHAjBCE1WnnFXulgmjtMSY/BIVEAAGMxehA6JK/YJxusPCFK3DHc6GsDAMCAjL7rf+0QH2zk6ED63/F4xBIAgBEZ+rKQggr+4JQ0qeKxy5UBJr2rAQAAPRi3R3iyjV8+qBHRgsvE8AHoDAIAGJRBg7BDo9x90qrRD+OUHycYdCUAAAAZNghfPqiVtnFyqFg0xqBrAAAAnIwYAx+U8vYKDlTpkQlKAB4uAQBgbIYLwuOt/KcSjYgWjlYSQ3FqEADA6IwVhGYHPbtX2iTdnKRkxhvruwMAwDkZKAyY6fdFWpWF08LFf4400BcHAIALMFAevHuS/13NIX70v+MVfwN9bwAAuBCjBMKhZn79sCYEPTBWjQ3GqUEAADjNEEHYaue8fdLBdOtQ5QcYVhsAALrw/SBkouf2y5p2Hj1Q3Jnm+98XAAAuie8Hw1vHZGEdh5vof8erKnqDAADwfT4ehIebeeN3LAT9cqwaFah3NQAA4Hl8OQjNDsrbpzkkz05WrohGZxAAAM7Bl4Pw5YOyup3SwsWdI3z5awIAQG/4bEJ8cEp+UilD/MUjE1Q/9AYBAOA8fDMIT7bxnw5JIrr/MiUmSO9qAADAg/lgEDqfNWiTdEOicm0cOoMAAHAhPhiEfyzWSts4KZTuxYCiAABwMX56F9DHPqumj8vZpIqHx6t41iAAAFyUTwVhdYd46bAkEv89Shkapnc1AADgDXzn4KFD8ivfBVoclBGrXD8EpwYBAKBHfKdHqArxgyiHbBaLRiMFAQCgp3wnCIWgGbH2uZeJQH8EIQAA9JTvHBp1whN3AQDgkiA3AADA0BCEAABgaD4VhCtWrKirq9O7CmP585//vG/fPr2rMJZt27Zt2bJF7yqMpbi4+LXXXtO7CmNpbW195pln3LMsnwrCt956q6WlRe8qjKWgoODYsWN6V2Es+/bt+/rrr/WuwlhKS0s/+ugjvaswFovF8ve//909y/KpIAQAALhUCEIAADA0BCEAABiaYGa9a7iQG264QVXVqKionky8devWq6++OiQkpL+rgk6ff/75kCFDhgwZonchBlJcXKxp2tixY/UuxECqqqqOHj2akZGhdyEG0tHRkZ+ff/PNN/dw+gceeGDy5MmuLcvTg7CysnLbtm16VwEAAB4tMzMzMTHRtXk9PQgBAAD6Fc4RAgCAoSEIAQDA0BCEAABgaAhCAAAwNO8LQk3TRo0a1bWlsbExKysrMjJy5syZjY2N3aa/8KfQE2ev8/fee2/s2LERERHXXHPN4cOHu02fnp4uzliwYIEbK/UpZ6/2C69YbOq9d/Y6F2fp+ik29V46557E/bt0LwvC3//+99OmTTt06FDXxtzc3OTk5MrKyqSkpLy8vG6zXPhTuKiz1/mpU6fmzZv32muvVVZWzpw585577uk6PTOXlJSUlZW1tra2trauWrXK7SX7grNX+0VXLDb1Xjrn7qW1iyeffPKRRx7p/Aibei+db0+iwy6dvUpBQcH777/frewRI0YUFxczc3Fx8YgRI7rNcuFP4aLOXufbt2+/9957na9ramqioqK6Tl9ZWRkaGjp58uTQ0NDs7Ozq6mq3lusrzl7tF12x2NR76Zy7l0779++fPn263W7vbMGm3kvn25O4f5fuZUHo1G1LDQkJsVgszGyxWMLCwrpNfOFPoYfOuXdwOBwLFixYuHBh18Y9e/ZkZmbu2bOnvr7+rrvuuv32291Vow/qutovumKxqfeJc27qVqv1iiuuOHDgQNdGbOp9pduexP27dF8IwuDg4Pb2dmY2m83BwcHdJr7wp9BDZ+8dtm3bdvnllz/yyCNd/0fupqKiYuDAgf1cmi87X+/knCsWm3qfOOc6f/bZZxctWnSBubCpu+zsPYn7d+l+fXOAVVfx8fGlpaVpaWnl5eUJCQmX9Cm4gJkfe+yxzz777I033hgxYkS3TwsLCzs6OqZNm0ZEJpMpICBAjxp90EVXLDb1fqJp2quvvpqfn9+tHZt6L51vT+L+XbqXXSxzTllZWevWrWPmdevWZWdnOxt37NhxgU+hN3bv3v3OO+/84x//iI+Pb2tra2trc7Y717nZbJ49e3ZxcbHNZlu6dOmsWbP0rNWHnG/FYlPvbwUFBYmJicOHD+9swabeJ863J9Fhl94n/Uo361Z2Y2PjTTfdlJCQkJWV1dTU1G2ac34Kl6rrOl+2bNk5tyLnCynl6tWrU1NTBw0adNdddzU3N+tTsU/outrPt2Kxqfets/eKP/vZz5YsWXL2NNjUe+l8exL379Ix6DYAABiaLxwaBQAAcBmCEAAADA1BCAAAhoYgBAAAQ0MQAgCAoSEIAQDA0BCEAABgaAhCAAAwNAQhAAAY2v8HEkarlhv3C7wAAAAASUVORK5CYII=" }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# plots setting\n", "gr(fmt=:png)\n", "\n", "mcm = McCallModel()\n", "@unpack V, U = solve_mccall_model(mcm)\n", "U_vec = fill(U, length(mcm.w))\n", "\n", "plot(mcm.w, [V U_vec], lw = 2, α = 0.7, label = [\"V\" \"U\"])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The value $ V $ is increasing because higher $ w $ generates a higher wage flow conditional on staying employed\n", "\n", "At this point, it’s natural to ask how the model would respond if we perturbed the parameters\n", "\n", "These calculations, called comparative statics, are performed in the next section" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## The Reservation Wage\n", "\n", "Once $ V $ and $ U $ are known, the agent can use them to make decisions in the face of a given wage offer\n", "\n", "If $ V(w) > U $, then working at wage $ w $ is preferred to unemployment\n", "\n", "If $ V(w) < U $, then remaining unemployed will generate greater lifetime value\n", "\n", "Suppose in particular that $ V $ crosses $ U $ (as it does in the preceding figure)\n", "\n", "Then, since $ V $ is increasing, there is a unique smallest $ w $ in the set of possible wages such that $ V(w) \\geq U $\n", "\n", "We denote this wage $ \\bar w $ and call it the reservation wage\n", "\n", "Optimal behavior for the worker is characterized by $ \\bar w $\n", "\n", "- if the wage offer $ w $ in hand is greater than or equal to $ \\bar w $, then the worker accepts \n", "- if the wage offer $ w $ in hand is less than $ \\bar w $, then the worker rejects \n", "\n", "\n", "If $ V(w) < U $ for all $ w $, then the function returns np.inf\n", "\n", "Let’s use it to look at how the reservation wage varies with parameters\n", "\n", "In each instance below we’ll show you a figure and then ask you to reproduce it in the exercises" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### The Reservation Wage and Unemployment Compensation\n", "\n", "First, let’s look at how $ \\bar w $ varies with unemployment compensation\n", "\n", "In the figure below, we use the default parameters in the McCallModel type, apart from\n", "c (which takes the values given on the horizontal axis)\n", "\n", "\n", "\n", " \n", "As expected, higher unemployment compensation causes the worker to hold out for higher wages\n", "\n", "In effect, the cost of continuing job search is reduced" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### The Reservation Wage and Discounting\n", "\n", "Next let’s investigate how $ \\bar w $ varies with the discount rate\n", "\n", "The next figure plots the reservation wage associated with different values of\n", "$ \\beta $\n", "\n", "\n", "\n", " \n", "Again, the results are intuitive: More patient workers will hold out for higher wages" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### The Reservation Wage and Job Destruction\n", "\n", "Finally, let’s look at how $ \\bar w $ varies with the job separation rate $ \\alpha $\n", "\n", "Higher $ \\alpha $ translates to a greater chance that a worker will face termination in each period once employed\n", "\n", "\n", "\n", " \n", "Once more, the results are in line with our intuition\n", "\n", "If the separation rate is high, then the benefit of holding out for a higher wage falls\n", "\n", "Hence the reservation wage is lower" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Exercises" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Exercise 1\n", "\n", "Reproduce all the reservation wage figures shown above" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Exercise 2\n", "\n", "Plot the reservation wage against the job offer rate $ \\gamma $\n", "\n", "Use" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.05:0.0375:0.95" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "γ_vals = range(0.05, 0.95, length = 25)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Interpret your results" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Solutions" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Exercise 1\n", "\n", "Using the solve_mccall_model function mentioned earlier in the lecture,\n", "we can create an array for reservation wages for different values of $ c $,\n", "$ \\beta $ and $ \\alpha $ and plot the results like so" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3daXzU1P4G8N/JtDNdBlraTlu2Fgp0A0uhoMhmuXoVkIJcFURZ5KLIRQS5CIrIIkVAXBDQP1wXkEVAUVnE5YqAIOBF9rUtW6EspS3d9+nknP+L0VrKtFQ6TaaT5/vhxUwmk/yShjyT5OSECSEIAABAqyS1CwAAAFATghAAADQNQQgAAJqGIAQAAE1DEAIAgKYhCAEAQNMQhAAAoGkIQgAA0DQEIQAAaBqCEAAANM2hg/DcuXMFBQXKz7esrEz5mWqcLMucc7Wr0Bxs6srDpu6AHDoIn3nmmYMHDyo/36ysLGypCisqKioqKlK7Cm3hnGdlZaldhebk5+eXlpaqXQXcRNEglGU5PDy84pBu3bqxP4wZM0bJYgAAAIjIRbE5LVq0aO3atUlJSeVDhBCJiYlXrlzx8vIiIhcX5YoBAACwUi57oqKiWrVqFRcXVz4kLS3NbDYPGDAgKSnp/vvv//DDD93c3Cp+pbi4+OOPP962bVulSUVERDzyyCN1V2phYaG7u7skOfR5YydTWFjIGFO7Cm3hnBcWFqpyGV7LCgsLy8rKZFlWuxBn4+bmdsdHU8oFYa9evSoNuX79eufOnd99992goKCJEydOmDBh3bp1lcbhnN96uY5zXqePURR/qLtZQCXWtY11riRs56rAandATOG/B2O255iamtq2bdtKl+5jY2NnzZoVGxurUHF/SEtLM5lMOCJUUn5+PmPMaDSqXYiGcM4zMjICAgLULkRbcnJyDAaDu7u72oXAn9Tc1x8+fHjfvn3W13q93mAwqFgMAABok5pBWFhYOHDgwISEBLPZHB8fX6eX/QAAAGxSMwi7d+8+c+bMuLi4pk2bZmdnv/nmmyoWAwAAjqNYweZESt+xUPECIWNs7NixY8eOVbgGAABwWBYufroq1pzjL7bTdTIp0Zgct+4BAIBDEIJ2Xxerz8ppxUREv2UIBCEAAGhFQo5YkcQTcgQRNTeyJ1ux7oEKXbxDEAIAgJouF4i158We65yIfN3oiRDpwWaSpGAHGwhCAABQR2YJrb/Af7zCuSA3HfULlgaFSO46pctAEAIAgNJKZPo6mX+VzM2cXCT2YDP2VGvJW69OMQhCAABQjrVR6GfneI6ZiKiziT0brmvsoWZJCEIAAFCCELQ3ja86y1OLiIjCvNg/w6TIRur3to8gBACAOpeYI5b/0Si0mSd7qrVyjUJvC0EIAAB16Eqh+Ozc741CvfRsUAh7OEjSqX8c+CcEIQAA1ImsUlp3nm+7wmVVG4XeFoIQAOrQ/v37x4wZo3YVDoRzzhhz+sdQC0GFFiqyWDvVZO4uZHSh9xi9V4PvhoWFrV+/vo4LvAmCEADqUF5ensFgWLp0qdqFQP1w5syZuXPnKjxTBCEA1K0GDRp06NBB7SqgflDlieiO0mgHAABAFQhCAADQNAQhAABoGoIQAAA0DUEIAACahiAEAABNQxACAICmIQgBAEDTEIQAALfhUD2iTZ8+PTAw0F5Ts/ui2bc8ZSAIAQDqk6VLl548ebKWE3nggQesL15++eVaV3QTu5SnMAQhAEB9kpmZ6efnV8uJbN++3fpi/vz5ta7oJnYpT2EIQgBQWmoRXcgTyvwrkSvPvWvXrt999x0RvfzyywEBAUIIIURAQMCJEyeqqXnmzJkmk6lr167nzp2zDtmyZUt0dLS3t3fjxo3ffvtt60DG2OrVqyudG7Q55vLlyxs3buzn57d48eJK87I5vtUjjzxCRNHR0XTzWc3y14yx9evXt2/f3tfX97333iOiGzduPPbYYyaTKSws7KuvvrI5kaysrOHDhzdu3LhJkyYjRozIysqqamrlbH6l4pStbp27A0Kn2wCgtIUnZOuTyhWw4B5dhPdNl8H69ev3008/9e3bd8eOHe7u7qdOnRJCuLu7t2vXrprpcM5TU1NnzJgxfvx4a45Onz592LBhEydOPH78eJcuXV566SXrmPv37y8/3rKyOeakSZN2795tMBief/758ePH33Z8q02bNjHGjh49Wk2pKSkpR48e3blzZ79+/V588cWJEyeGhISsX7/+119/HTBgQL9+/W6dyIsvvqjX6y9cuEBE48ePnzRp0ooVK2xOrfqv3DrlW+duMBiqKV4VTAiFNsc7EBsbO2vWrNjYWIXnm5aWZjKZVOkEXbPy8/MZY0ajUe1CNIRznpGRERAQUKdz2bZt24IFC7Zt21Zx4Hsn+YU8hfY8E++SWja4KQiPHz/+1FNP7dmzp0uXLgMGDGjcuLEQ4vz580uWLKlqIoyx8+fPh4SE3Lhxo1WrVrm5uUTEOT9w4MCpU6d27dq1atUq676UMZaenm4ymSp+3eaYAwYM4JwPHz584MCBLi4utx2/YjHl8yr/qOLAvLy8Bg0alA/08/NLSkry9fUloqysLC8vL51OV2kifn5+p0+f9vf3J6K0tLSoqKi0tDSbUysvo5qvVBrt1rlX8/c6duzY8OHDjx07Vs04docjQgBQ2ovt1PyVedddd+Xl5X3xxRc9e/a8//77lyxZIoSodExWFUmSXF1dra8HDRqk1+uHDBkyb968VatWlY9TKQWrGnPTpk0//vjjp59++uGHH1b6oVDVlKuSn59f8a01t8rJslz+sz41NdXNzc3Dw+PWiVQ8uSrLf55QrjS1mnzlDuauLhz0AIC2MMYefvjh+Pj42NjYbt26/e9//zt69Oh9991X/bc+/fRTi8Xy7rvv9urVyzpk27Zt06ZN69ev3w8//EBEFoulqu/aHLNly5YtW7acMWPGoUOHajL+rQwGw86dO4UQ1T/3+IEHHnjrrbdkWf7111+7dOlSPrWysrLycfr06TNt2rSSkpLi4uJp06b17du3+rVR869UNXeHgiAEAM3p16/f5cuXY2NjPTw82rVr16VLF71eX/1XXF1dGzduvGfPnnfffdc6ZO7cubGxsXfddVdmZuZDDz00atSoqr5rc8wpU6bce++9sbGxb731Vk3Gv9WcOXMeffTRqKio6s9vL1my5MSJE/7+/kOHDv30008bNmxIRH379m3VqlX5OO+9915xcXGLFi1CQkLMZnOldjE21fArNufuaHCN0AZcI1QerhEqT8VrhABVwTVCAAA12exmxZGPFsAuEIQAAL9D5mkTzv4BAICmIQgBAEDTEIQAAKBpCEIAANA0BCEAAGgaghAAADQNQQgAAJqG+wgBoA65urqeP3++hl1aQw1xQfIfdzxKjCRGNjoCqJ9u3Lih/EwRhABQh2JjYxcsWHDt2jW1C3EUxcXFLi4u5Y+w+EsEUXKe2J8h8i1ERM086J4Aydfhnu5XK61bt46Li1N4pghCAKhbjz32mNolOJCcnByDweDu7v5Xv3g0U6xIki/kkx9RRyM9HarrbHKa40CVIQgBABxaSgF9ekY+kCGIyM+NDQ5hDzaTJISg/SAIAQAc1I0S+vwC//EK54KMrvRoC6l/C0mPNo72hiAEAHA4BWX01UW+5ZIwy8KF0YPNpaGtmZceh4F1AkEIAOBALIJ+usLXnBO5ZsEYdQ+URrRhgR6IwDqEIAQAcAhC0N40vvKsuF4kiKi9LxsZKrVqiAiscwhCAAD1HcsUK87w83mCiJp70pOtpe6BuBioEAQhAICaLhfS2nN8z3VORL5u9ESIhEahCkMQAgCo40aJ+PyCsDYKddPRwBbSYyFoFKoCBCEAgNJKZPo6mX+VzM2cXCT2YDM2tLXkpVe7LK1CEAIAKMfCxe5r7IsUS66ZrI1Ch7dhjdEoVFUIQgAAJVgbhX58yi29lOl0FO7N/hkmRXgjAtWHIAQAqHMJOWJFEk/IEWVlFNSAhoehUagDQRACANShSwXi0zP8YIYgIl83GtDC8mBzF08PpKADQRACANSJGyXis3N8xzXBBXm40KMtpQHBUnG+LDHseB2Lcn8PWZbbtm2bmJhYafjJkyfvueeewsJCxSoBAKhThWXiq4ti80Vu5uTC6OFgaXDI7z2FFqtdG9xKoSBctGjR2rVrk5KSKg3Pzc19+umni4qKlCkDAKBOWbj46ar47BzPMRMRdTaxZ8N1jT3ULguqpVAQRkVFtWrVqtJzh4UQTz/99CuvvPL4448rUwYAQB2xNgpddZanFhERhXmxf4ZJkY3QKLQeYEII5WbGbprd/Pnzb9y48fbbb1caXq5z586enp4mk6nS8E6dOo0ePbru6szIyPD19ZUkXM1WTkFBAWPM09NT7UI0hHOemZl56/8vuANn86U1yS5n8iUiauIuHg8q6+LHbY6Zm5trMBjc3NyULdD5GY1GV1fXO/uuatdsd+7c+cMPP2zbtq2acfR6fZcuXcLDwysNDwoKcnd3r7va3Nzc3N3dEYRKslgsjLE6/bNCJZxz66audiH129UiWp8s7UsnImrkRo8G8z7NmI4ZqhrfbDbr9XqsdrurzR5btSDcvn37rl279Prf+xRijP3yyy/du3evOI6rq2vv3r1jY2MVrs36ew1BqKSysjLGGH4mK4lzjkOT2sgqpXXn+bYrXBbk4Ur9gqVBIZK77jbfKikpwWp3NKoF4Zw5c+bMmWN9XdWpUQAAB1Qi09ZL/IsLvFgmF4n1bsaebMUaGXA5sL7C7SwAADV1a6PQZ8KlJugptJ5TNAirOuzD4SAAOL7fMsQnifxakSCiMC82Mkxqi0ahTgFHhAAAt5GQIz5J4kk5goiaebIRoVIXf0Sg80AQAgBU51AGn3NUWLhoZGBPtpb+3pTpEILOBUEIAFClpFwx/xi3cOobJI0Mldxu1ygU6iMEIQCAbZcL6fVDcolMvZqwMeESw4Ggk0IQAgDYcKNEzDrE88vobhOb0E6HFHRiuGccAKCyPDPNOMTTi0W4N5vSXoeLgs4NQQgAcJNSmeYckS8XiBYN2MyOOgOuCzo7nBoFAPiThYt5x0RCjgj0YK/HSMY77MYZ6hMEIQDA7wTRuyf4oQzhpWezYnQ+VXadDU4Fp0YBAH73SRL/5brwcKHXY6SmeJquZiAIAQCIiNad55svchdGU6N1rRqieYyGIAgBAOj7y3ztOS4xeqm9LtoXKagtCEIA0Lr/pYtlCZwxej5S1y0AKag5CEIA0LTjWWLBMZkLGtZaerAZUlCLEIQAoF1nc8Wcw3IZp4eDpMdDsD/UKPzhAUCjUoto9mG5WKb7GrPnIrAz1C787QFAi7JK6bUDlhwzRfuyF+9CH2qahhvqAUBzCsvErEM8vYTCvNm0DjoXxKC24YgQALTFzGn2EZ6cL4KMNLMjHjEIOCIEAC2xcDHvKD+dLfzc2KwYXQN0JQo4IgQA7RCCPjjND2aIhnqK7ySZ3NQuCBwDghAAtGL5Gf7TVWHQ0fQOumaeuDAIv0MQAoAmfHGBb7rIXRi9Gq0L90YKwp9wjRAAnJxF0MZkvuYcZ0STonQd/ZCCcBMEIQA4LSFob5pYdZanFgki+leE1D0QKQiVIQgBwDmdzBYrkviZXEFEzY3s6VDpbhNSEGxAEAKAs7lSKD47J/Zc50TUUE+DQ6SHgyR0HgNVQRACgPPIKqV15/m2K1wW5KajfsHS4y0lD+znoFrYQADAGZTItPUS/+ICL5bJRWK9m7EnW7FGBhwGwu0hCAGgfrNw8dNVsfa8yC4VRBTty54Jl4KNiECoKQQhANRjv2WITxL5tSJBRGHebGSo1LYRIhD+GgQhANRLSbliRRI/lS2IqKkHDW0jdQuQGEIQ/joEIQDUM1eLaM1ZvjeNC0EN9TQwWHqkpYRHKcEdQxACQL2RZ6b15/l3l9EoFOwJWxAA1AOlMn2ZzDdd5CUy6Rj1aS4NQaNQsBMEIQA4OgsX846JQxmciLr4sxGhEp4dAXaEIAQAhyaI3jvJD2UILz29Gq2LRKNQsDcbj2ESQixdurRHjx6BgYGJiYlTp07dunWr8pUBABDRJ4nyrlThrqNZMUhBqBM2gnDlypUvvPDC3XffnZaWRkR5eXlxcXEbN25UvDYA0Lr15/nmS8KF0dQOutYNkYJQJ2wE4cKFC6dMmfLOO+9Y377//vvDhg2bN2+esoUBgNb9cEV8do5LjCZFSR18kYJQV2wE4ZkzZ+67777yt4yxgQMHnj59WsGqAEDr9qeLpadlxmhspNQ90MaeCsBebGxewcHBGRkZFYdcvHgxKChIqZIAQOuOZ4kFxzkXNLS19FAzpCDULRtb2IQJE+bPn3/+/HkiKi4u3rp168yZMydPnqx4bQCgRRfzae4R2SyLvs3ZoBCkINQ5G7dPjBkzxmg0Dho0yNPT86GHHmrbtu2aNWv69++vfHEAoDWpRTT9oKXQQj0D2ZhIndrlgCbYCELG2LBhw4YNG6Z8NQCgZVml9NoBS46Zon3Zi3fhkfKgENxQDwAOobBMvH6Yp5dQqBeb1kHninOioBQb2xq7haura+PGjTt16vR///d/ZrNZ+SoBwLmZOcUf4RfyRHNPmhWjc8M5UVCQjSPCX375ZciQIcOHD+/Tp48kSd9///26des+/PDDlJSUWbNmpaSkzJ8/X/lCAcBZcUFvH5dPZQs/N/Z6J10DV7ULAo2xEYTvvvvuiBEj5syZY33btWtXIcSKFStWr14dHBzcr1+/uXPnShJOWwCAHQhBS07Jv6aJhnqaHSOZ3NQuCLTHRp7t3Lmze/fuFYd069bthx9+IKLOnTsXFRVdvXpVoeoAwNmtOMN/uioMOnqtg665Ee1jQAU2gjAwMPD48eMVh5w4ccLb25uILl++TEQ+Pj7KFAcAzu3LZL7xIneR2NT2UoQ3UhDUYePU6Lhx41555ZXAwMCHH36YMfbtt9/Gx8fPnj07OTl57NixvXv39vT0VL5QAHAyO6+JVWc5I/r3XVKMCSkIqrERhGPHjmWMTZ06dcSIEURkMpni4+PHjx+/du1af3//999/X/EiAcDZ/JYhFp2UhaBnw6UegUhBUJPtG+rHjh07duzYzMxMi8Xi7+/PGCOioUOHDh06VPEKAcDZJOaIBcdkWdCQVlL/YLS8A5VVtwn6+voGBAQwxsxm88WLF5UqCQCc2cV88fphuVSmPs2lJ1sjBUF9tnuWyc3NTU1NLX/722+/jRs3Li8vT6mqAMA5XS8SMw/xgjLq4s/GRCAFwSHYCMJ169YNHTqUc14+RKfTjR49WsGqAMAJ5ZrFrMM8q1RE+bAp7XUSrgyCY7Dxiyw+Pv7ZZ5/NyckJDQ1NT09PS0uLiopCEAJAbVwupOkH+dVC0caLvdYRXYmCA7GxMV64cKFv375eXl6xsbHHjh3z9/efPHnypEmTaj8zWZbDw8MrDvnhhx8iIyO9vb0jIyN//PHH2s8CABxNVim9f4qP22tJzhdNPWhmR8kdXYmCI7ERhD4+Punp6UQUFRW1e/duImrevPmBAwdqOadFixZ17do1KSmpfAjn/KmnnlqyZElWVtbs2bNHjhxZy1kAgEMplmnNWf7cL5b/XuESY32bszfv0XnpcUoUHIuNa4SdO3deuHBhx44do6OjR40a9fzzz69bty4gIKCWc4qKimrVqlVcXFz5EIvFsnr16r/97W8FBQUGg8HaeQ0AOAELFz9dFWvPi+xSQUTRvuyZcCkYPaiBQ7IRhG+88UavXr22bNkyc+bMTp06BQYGenh4rF27tpZz6tWrV6Uher2+b9++BQUFDRs2ZIzt2bOn0gj5+fnPPvus0WisNDw2NnbatGm1rKcaWVlZjDF0LK6kgoICxlhJSYnahWgI5zwrK0uns/9pysPZus8u6q+XMCJq3YA/0dwc4cWphG7gz0uUm5trMBjc3NC5uJ01bNhQr9ff2XeZEOLWoaWlpcXFxdZDtIyMDKPR6O7uXqsay+fHbMyxsLBw0aJFGzdurHQCtkePHk888UTHjh0rje/j4xMSEmKXemxKS0szmUwIQiXl5+czxm790QN1h3OekZFR+5M9FZ3Jo5Vn6XQOEVETD3oyRHT9vUMO+F1OTo7BYLDXHhXK6XS6O95p2zgitFgsBoPBYDBY35pMpjsvrVoXL1784IMP3nrrLU9Pz1GjRs2dO7fSCDqdrm3btvfee28dFVAVV1dXV1dXBKGSXF1drY+AVrsQDeGcWzd1u0ztahGtOcv3pnEhyMtAA4OlR1pKLojAW7j+Qe1C4E829vX+/v5Dhw79/PPPc3Jy6nTeTZo0+eSTT3bt2iWE+Pzzzzt06FCnswOAupBnpg8T+PN7LHuuc4NEj4VIH/VweSwEKQj1hu2nT2zfvt3arWjPnj3j4uL69evXunVru89br9dv3Ljx3//+d3Jycnh4+PLly+0+CwCoOyUybb3ENyTzIgsxol5N2NOhko8BAQj1jO1rhESUn5+/e/fuHTt2bN++/dixYxEREadPn1a4uNjY2FmzZsXGxio8X1wjVB6uESqvNtcIZUHbrvB150XWH41CR4VJLRogAm8P1wgdkO2+Rq3YH4goLS1NqZIAwKEdzRQfJ/JLBYKIQr3YyDCpXSNEINRjNoJw5syZ27dv379/v4uLS/fu3Z944omPP/4YF/AAIClHfJLEE3IEETXxYMNDJTQKBSdgIwhnz57t4uIyatSoV199NSgoSPmaAMABncgSsw5zsyy89DSklfRQM+aCbrPBKdgIwrVr1/788887d+78z3/+ExkZ2atXr9jY2J49e/r7+ytfHwA4ggv5NOeIbJbp702lZyPQWSg4FRtBOGTIkCFDhhDRtWvXdu3atXXr1sGDB3POq2pWAwDOLbVIzDrEiyzUs7H0QjscBoKzsd1YRpblw4cPb9u27ccff9y3b5+np+dDDz2kcGUA4AhyzTTrEM8uFe192YtIQXBGNoLw8ccf3759e3Z2dkhISFxc3GuvvdazZ8877sMNAOqvwjIx4xC/ViTaeLHXOuAhguCcbARhWlraK6+8EhcXFx4ejgZhAJpl5hR/hF/IE0082MyOOjdcFwQnZSMIrc8gBAAt44LeOc5PZQtfN4rvJHnhlBA4r+puqAcAbRKC3j8l70sTDVwpPkbn744zQ+DMcMofACpbcUbedlUYdDS9o645nqYLzg5HhABwk6+S+caLwkViU6OlCG+kIDg/HBECwJ92XuMrz3LGaFwki/FDCoIm2D4iLCwsvHz5cqWB4eHhdV8PAKjmtwyx6CQXgp4Jk+5vil/JoBU2gnD9+vUjRowwm82VhqNnGQAnlpQr3jomy4KeaCUNaIEUBA2xsblPnTp1yJAh+fn54mbKFwcAyrhUIGYdkktk6t2MPdUaKQjaYmOLz83NHTNmDB6RCqAR14vEjIO8oIy6+LN/ReK2edAcG0H4wAMPHD9+XPlSAEB5uWYx65CcVSqifNiU9jr0JQoaZOMa4UsvvTR+/Pji4uLu3bt7enqWD0djGQAnUyyzeYf41SJq2YC9iq5EQatsBOE999xDRPv37680HJcJAZyJRdDis27nS0RjD5rdSeeJm4pBq2xs+wg8AKdn7Ur0ZK6ucUOa09nFG12JgobhVAiA5ghBH5yW96YJDxeaFaPzd1O7IABV2Q7Cb775pkePHn5+fj4+Pj169Pj2228VLgsA6s7qc/zHK0Iv0aTQ4pYN0DwGtM5GEG7YsOEf//hHjx49Nm3aZE3EAQMGfP3118oXBwB2920K33CBu0js5fYstIGsdjkA6rNxjXDu3LlTpkx54403rG+7devGOX/jjTf+8Y9/KFsbANiThYtvUmjFGc4YjW8rdTaJjAy1awJwADaOCM+cOdO9e/eKQ+67776kpCSlSgIAOxOC9qaJ5/fy5UmyEDQqTNerCc6IAvzOxhFhcHDwqVOn+vTpUz7k5MmTwcHBClYFAHZzOlssT+JJuYKImnrQ02G6Lv5IQYA/2QjC0aNHT58+PSAg4OGHHyaib7/9dvbs2fHx8YrXBgC1crWI1pzle9O4ENRQTwODpUdaSi4IQYCb2QjC8ePHWyyWiRMnDh8+nIh8fX1nzJgxfvx4xWsDgDuUZ6b15/l3l7ksyE1H/YKlx1tKHrhlHsAWG/8zJEl66aWXJk2alJGRQUQmk4kx/IYEqB9KZNp6iW9I5kUWYkS9mrCnQyUfA/4LA1Spyp+IjDF/f38lSwGA2pAFbbvC154X2aWCiKJ92agwqQVuEwS4nd+DMDExkYhCQ0MlSbK+vhU63QZwWEczxceJ/FKBIKJQLzYyTGrXCBEIUCO/B2FERAQR5efnG41G6+tboQ9SAAd0JlesOMNPZgkiauLBhrVh3QIkXM0AqLnfg7BiyCHwAOqFq4VizTlhbRTawJX+0UIa0ELCo5QA/iob1wgZY6mpqYGBgeVDDh8+/MADD2RlZSlYGABUKc9MGy/xTcncIsigo7hg6bEWzNMVh4EAd+LPICwoKLhy5Yr19dmzZ3Nycso/2r17d1lZmdKlAcAtzLL4+qL4OpkXyyQxeqiZ9GRrycegdlkA9dmfQfjzzz/HxcVZX/fs2bPiSDqdbuLEiYrWBQC3kAUtOM73pwsiutvEng7TNfdUuyaA+u/PIOzXr5/16uCtp0YBQHVC0Aen+P500cCVTY2W7vLBiVAA+7BxYV0IUSkFzWbzxYsXFaoIAGxZeZZvu8oNOpreESkIYE+2b6jPzc1NTU0tf/vbb7+NGzcuLy9PqaoA4CZbU/hXydxFYlOjpQhvpCCAPdkIwnXr1g0dOpRzXj5Ep9ONHj1awaoA4E87r/EPEzhjNK6tFOOHFASwMxunRuPj45999tmcnJzQ0ND09PS0tLSoqCgEIYAqDmSIRSe5IBoVprsfDxEEqAM2gvDChQt9+/b18vKKjY09duyYv7//5MmTJ02apHxxABqXlCsWHJNlQYNDpAHBSEGAOmEjCH18fNLT04koKipq9+7dRNS8efMDBw4oXRqAtl0qEK8f4iUy9W4uDW2DDmMA6oqN/12dO3deuHDh4cOHo6Ojv/jii7S0tHXr1gUEBChfHIBm3QynnxgAAB9vSURBVCihWYfk/DJxjz8bE4EUBKhDNhrLvPHGG7169dqyZcvMmTM7deoUGBjo4eGxdu1a5YsD0KZcs3jtgHyjhO7yYVOiJB3OiQLUJRtB2K5duytXrhQXFzPG1qxZs3DhQqPR6O7urnxxABpULNOsQ/xqEbVoQNM66PQ6tQsCcHY2TrksWbKkqKjI29vb+tZkMiEFAZRhETTviHwuTzT2YLNjdJ5VPjkbAOzGRhC+9NJLjRs3HjZs2O7du/FIJgDFCKJ3jvMjmcLHQPGddI0MOCUKoAQbQZiamrp48eKUlJT77rsvPDz87bfftjYiBYA69VEi33Ode7qymR11ATgLA6AU27dPjB49eteuXZcuXRo5cuTKlSubNWs2aNAg5YsD0I7VZ/k3l7heoukdpJCGOBYEUE51zbKDgoIeeeSRxx9/3MvLa8OGDYrVBKA1310WX1zgEqNJUVLbRkhBAEXZDsJz587NnTs3Ojo6IiJi9erV48aNO3v2rMKVAWjE7uti2WnZ2pVo1wDcMgigNBuN0jp16nTo0CFfX9/BgwcvW7bsnnvuYQw/UQHqxNFM8d4JWRD9M1T396b4jwagAhtBGBISMmPGjN69e+v1euULAtCOM7nijSNyGafHQqSBLZCCAOqwcR7m0qVLhYWFSEGAOnWtSMw+LJfI1KuJNLw1zogCqMbGfz9/f/8dO3YoXwqAdtwoEdMP8lwz3W1iE9pJuPgAoCIbp0ZnzJgxduzYKVOmxMXFmUym8uHh4eEKFgbgtPLLaMYhnl4swr3ZlPY6dCUKoC4bQXj33XcT0cGDB996662Kw2vfy4wsy23btk1MTCwfsnnz5mnTpl25ciUqKurjjz8ODQ2t5SwAHFyJTLMOyZcLRIsGNCtGZ0BXogBqs3FqVFShlnNatGhR165dk5KSyoekpKQMHTr0o48+Sk1N7d+//8iRI2s5CwAHZxE094h8JlcEerDX0ZUogGNQ7hJ9VFTU9OnTKw65cOHCE088ce+997q7u48YMaJiRgI4H0H07nH5SKbw0rNZMTofdCUK4BjYrYd6Qohly5atXbv27NmzP//888qVK7t169avXz/7zI/ZmKMsy+PGjZMk6YMPPqg4PDo6+uzZszpd5ZNHffr0Wbx4sV3qsSkjI8PX11eS0JBPOQUFBYwxT09PtQupQ5+lGH5IdXXXiVcjilt4crXLIc55ZmZmxXYAoIDc3FyDweDm5qZ2Ic7G29vbYDDc2XdtnJpZuXLlCy+8MGHChD179hBRXl5eXFzc119/PXDgwFqVWYWffvppypQpDz744Jw5cyp95OXltX79+u7du1cartfr6/TJUJxzk8mEIFSSm5sbY8xoNKpdSF1Zd17syOQeBprRUYr29VC7HCIizjkRIQgV5urqajAY8Gw7u6vNHttGEC5cuHDKlClz58599913iej999/Pz8+fN2+e3YNQCPHqq6/u3bt3/fr1NpvJMMYaNGjQqFEj+873tqQ/KDxfLZMkiTHmrOv8+8t8/QWhk9jk9rqOJgc6I4rtXHnYvTggG3+MM2fO3HfffeVvGWMDBw48ffq03ee9b9++jRs3btmypUmTJgUFBQUFBXafBYDq/pculiVwxuj5SF23AAdKQQCwsnFEGBwcnJGRUXHIxYsXg4KC7D7vn3/+OSkpqeIBH54DDE7meJZYcEzmgkaESg82QwoCOCIbR4QTJkyYP3/++fPniai4uHjr1q0zZ86cPHmyXeZXMeqmTZtm3zs0ABzK2Vwx57BcxunhIOmxljgVBuCgbBwRjhkzxmg0Dho0yNPT86GHHmrbtu2aNWv69++vfHEA9VdqEc0+LBfLFNtYei4CKQjguGwEIWNs2LBhw4YNU74aAOeQWUKvHbDkmCnal024S8IpUQBHhp4tAOyssEy8fpinl1CYN5vWQeeCGARwbDbO2CQkJHTu3Hn16tVENHnyZKPR2KtXr8uXLyteG0D9Y+Y0+whPzhdBRprZUXJDV6IADs9GEL7wwgsBAQF9+vRJSUlZtmzZl19+qdfrJ06cqHxxAPWLhYt5R+XT2cLPjWbFuDRwxcEgQD1g49To/v37P/74Yz8/vwULFvTp06d37965ubnPP/+88sUB1CNC0Pun+MEM0VBP8Z10JnShBVBP2DgidHFxYYwR0d69e6131ru7u5vNZqVLA6hXlp/h268Jdx29HqNr5oljQYB6w0YQ3nPPPd98801CQsKOHTsef/zxsrKy9evXt2/fXvniAOqLLy7wTRe5C6NXonWtGyIFAeoTG0G4YMGCn376KTIycvDgwf7+/s8///yOHTsqPaQXAMptv8rXnOOMaFKU1NEPKQhQz9i4RhgVFXX58uW0tLQmTZoQ0Ztvvrl06dJbn4UEAES0P10sOcWFoNERUvdA3DgPUP/Yvo/QxcWladOm1tfKP/wBoL44kSUWHOeyoKFtpLggpCBAvWTjv64QYunSpT169AgMDExMTJw6derWrVuVrwzAwV3MpzeOyGZZ9G3OBocgBQHqKxv/e60P5r377rvT0tLojwfzbty4UfHaABxXapGYcUgutFDPxtKYSFw4AKjHbASh9cG877zzjvXt+++/P2zYsHnz5ilbGIDjyiql6Qfl7FLR3pe92A5diQLUb2o+mBegPiq00MxDlrRiCvNmr3XQueKcKEA9Z+M/sWIP5gWod8ycZh+WL+ZTcyOb2VGHrkQBnIDSD+YFqL+4oHeO89PZwteNZnWUGriqXRAA2AMezAtQI9auRPel8QauFB+j83fHlUEAJ4EH8wLUyIoz8rarwqCj6R11zY1IQQDnYbuv0XXr1ilfCoDD+iqZb7woXCQ2tb0U4Y0UBHAqNoLQ399/x44dypcC4Jh2XuMrz3LG6IW2UowJjUQBnI2NU6MzZswYO3bslClT4uLiTCZT+fDw8HAFCwNwCL9liEUnuRD0TJj0tyY4FgRwQjaC8O677yaigwcPVnrihBBCoaIAHENSrnjrmCwLeqKVNKAFjgUBnJONIETgARDRxXwx65BcIlPvZuyp1khBAKeF/94ANlwvEjMP8YIy6uLP/oWuRAGcmu3HMAFoWa5ZzDokZ5VSlA+b0l6HvkQBnBuOCAFuUmShmYf41SJq2YC9iq5EATQA/8sB/mQRNO+ofD5PNPag2Z10njhjAqABCEKA33FBbx+Tj2YKHwPN6ezirVe7IABQBH7xAhARCUEfnJb3pglPVzYrRvJ3U7sgAFAKjggBiIhWneU/XhF6iWZ0kFo2QPMYAA3BESEAfZvCv0zmLhKbGi1FNkIKAmgLjghB63aliv8kcMbo+UipkwkpCKA5CELQtKOZ4r0TsiAaGSo90BQpCKBFCELQrqQc8cYR2SLo8RBpILoSBdAq/OcHjUopoNcP8xKZ/taEDUNXogAahsYyoEU3SsSsQ3J+Gd3jz8a30zGcEwXQMPwQBs0RRG8e4xkl1K4RmxIl6ZCCANqGI0LQnG1XRGKO8DHQax11ejxYAkDzcEQI2pJfJlaelYloVJiErkQBgBCEoDUrkniemaJ9Wc/G2PgBgAhBCJpyJlf8dFW4MHouAls+APwOuwPQClnQB6e5IHosRGrmiRYyAPA7BCFoxTeX+IU84e/OHm2JzR4A/oQ9AmhCdqlYf0EQ0ZgIyQ0tRQGgAgQhaMKHibywTNwbwDqjW20AuBmCEJzf0Uyx57ow6OiZcBwMAkBlCEJwcmWclp2WiWhIKzx3HgBsQBCCk/symV8touaeNCAYJ0UBwAYEITiz1CLxZbJgjMZE6lwkBCEA2IAgBGf2UaIwy6JXYxblgxQEANsQhOC09qWJAxnc05WNCMV2DgBVwg4CnFOpTJ8kykQ0vDXzMeBwEACqhCAE5/TZOTm9hFo3ZL2bYyMHgOpgHwFOKKWAvkkhRjQmQkITGQCoHoIQnI0QtCxBtnDRN0gK80YMAsBtIAjB2exIFSeyhLeehrbG5g0At4c9BTiVgjJakSQT0T/DJKOr2tUAQH2gQhDKshweHl6TgQB/1aqzPNdMbRuxWDyAHgBqRumdxaJFi7p27ZqUlHTbgQB/1dlc8d8rXMdoTITEcHEQAGrGReH5RUVFtWrVKi4u7rYDiUgIkZmZmZqaWmm4m5ubl5dX3RXJOeec19304Vacc8ZYbVa7IFp6mstcPNZSCvIUnAs7lueU+B/ULkRbsNrriCTd+XEdE0KF/QVjNuZ768AOHTokJye7ula+1NO7d++333677srLyMjw9fWtzWqFv6qgoIAx5unpecdT+PG66+pLBh89X9C+2CAhBW+Pc56ZmWkymdQuRFtyc3MNBoObG56EYmfe3t4Gg+HOvqv0EeFf4uXltWnTptjYWOVnbTKZEIRK8vDwYIwZjcY7+3qOmb49aTEY6MUOuiD/BvatzVlxziVJCggIULsQbTEYDAaDwd3dXe1C4E/Y14Mz+CSRF1ooxiR18ce1QQD4axCEUO+dzBK7rnO9RGPCkYIA8JchCKF+swj6vwQuBA1uJQV6IAgB4C9TJwhtttBRpdkO1HdfJ/PLBaKJBxvYAr/qAOBOYN8B9VhGCW24wIlobKTkim0ZAO4Idh5Qj/0nQS6RqWcga++Lk6IAcIcQhFBfHboh9qcLDxf6Z5hO7VoAoB5DEEK9VCrT/52SiWhoG8kXtyYDQC0gCKH+EUTvneTpJRTSgB7GA+gBoHYcumcZAJuWJfA917mnK5t4lw4PoAeAWsKvaahnPjvHv0vheh17LVpqgc7UAKDWEIRQn3yXwtef5xKjSXdJ7XxwMAgAdoAghHpjdypflsAZo3FtdV0DkIIAYB8IQqgfjmeJ905yQTSijfT3pkhBALAbBCHUA2dzRfxhuYxTXLD0aEtstABgT9ingKO7ViRePyyXyNSrifRsGLZYALAz7FbAoWWW0PSDPNdMnU1sQjuJ4ZwoANgb7iMEx5VfRtMPWtKLKcyLTWmv0yEFAaAO4IgQHFSpTPGH5cuFFGxkM2MkN/QnCgB1A0eE4IgsXMw7yhNyhJ8bzYrRNXBVuyAAcF4IQnA4QtD7p8WhG8JLz+Z01vmhT20AqEs4NQoO55MzfPtV7uFCs2Kkph5qVwMAzg5BCI5l/Xm++SJ3YTQ1Wte6IZrHAECdQxCCA/nhivjsHJcYTYqSovHQeQBQBK4RgqM4kCktTZIZo7GRUvdA/EQDAIVgdwMO4WQOW3ha4oKGtpYeaobNEgCUgz0OqC85X7xz2qWMU9/mbFAItkkAUBR2OqCy1CKacVAutFB3fz4mErfNA4DSEISgpqxSeu2AJcdMUY3E8+EczWMAQHloLAOqKSwTsw7x9BIK9WIvhVtc0KM2AKgBQQjqsAh6/TBPzhdBRpoVI1GJ2gUBgFbh1CioY9NFnpAjTG70eoyugSuOBQFANQhCUEFGCX1+nhPR8211fm5IQQBQE4IQVPCfBLlEpp6BLMYPKQgAKkMQgtIO3RD704W7jv4ZhpslAEB9CEJQlJnT0tMyEQ1tI/ni+UoA4AAQhKCoz8/ztGJq0YA9HIRtDwAcAnZGoJxrRWLjRc6Ino+UdLg4CACOAUEIyvm/07yM04PNpHBvxCAAOAoEIShkV6o4likauNLwNtjqAMCBYJcESiiy0PIkmYhGhkkN9WpXAwBQAYIQlLDmHM8qpQhv9kATbHIA4FiwV4I6dzGfvkvhOkZjIiR0rA0AjgZBCHVLEH1wWpYFxQVLIQ0RgwDgcBCEULf+e0Uk5ohGBvZECFIQABwRghDqUH4ZrT4rE9HocOaJR0wAgENCEEIdWp4k55kp2pd1D8SWBgAOCrsnqCsJOWL7NeEq0ZhIdK4NAI4LQQh1Qha09DQXgh5rKTX1ULsaAICqIQihTmy5xJPzRWMP9lhLXBoEAIeGIAT7yyoV689zIno2nOnRuzYAODYEIdjfhwm8yEJdA1hnEzYwAHB02E+BnR3JFHvThEFHo8LRRgYA6gEEIdhTGaf/nJaJ6KnWOn88gB4A6gMEIdjThgv8ahEFGSkuSO1SAABqBkEIdpNaJL66KBijMRE6FwltZACgfkAQgt18lCjMsvhbY3aXD1IQAOoNBCHYx57r/EAGN7rSyDC0kQGA+gRBCHZQKtOKM4KIRoTqvPAAegCoVxCEYAefnePpxaKNF3uoGU6KAkA9gyCE2kopoG9SBCMaE4EWMgBQ/yAIoVaEoGUJsoWLfsFSqBdyEADqH0WDUJbl8PDwikOys7Pj4uJ8fHz69++fnZ2tZDFgF9uviRNZwltPT7XGjyoAqJeU23ktWrSoa9euSUlJFQe++eabwcHBqampQUFBCxYsUKwYsIv8MvHpGZmIRoVLni5qVwMAcEeU23tFRUW1atUqLi6u4sCNGzdu3rzZYDCMGzduwIAB8+bNq/ipxWI5cuQI57zSpPz8/MLCwioNTMylz5Ptc2qusFDv4VHGGE703UZmicguYXc1oi6N5NLSWk2qtLSUMebq6mqn0uD2OOdms7m0ln85+IusK1yScAbFzlxdXe94rTIhhH2ruc382E1zNBqNGRkZ7u7uxcXFAQEBeXl5FUeOiYkpKSlp2LBhpYn06NHjlVdeqTTwQKZuYZJ9Wu6Xlpbq9XoEYU3oSMyPLm3qUdutqKCggDHm6elpl6qgJjjnmZmZJpNJ7UK0JTc312AwuLmhK147MxqNev0dRoDK57OEENa8EULIslzp0wYNGrzzzjuxsbE1mVQXI833sU+oZ2YWNmrUAD/ZasLoSq0b2uEJ9K6urowxo9FY+0lBDXHOZVn28fFRuxBtkSTJYDC4u7urXQj8SeUgbNKkyeXLl9u0aXP16tWmTZvWZlJeeor2tc8xXJpFNvkyCfcCAABogMoHPXFxccuXLxdCLF++fMCAAeoWAwAAGqRyEM6YMeP48ePNmzc/derUa6+9pm4x5dLS0hS+dAoFBQUFBQVqV6EtnPP09HS1q9CcnJyc4uJitauAmygdhJUCxtvb+9tvv71y5cqWLVu8vLwULqYqPXv2zMnJUbsKbVm0aNH777+vdhXakpmZWcML8GBHkyZN2rp1q9pVwE3QHgQAADQNQQgAAJqGIAQAAE1z6H6xzGbz999/n5ycrPx8161bh5u7lXTixAlJklasWKF2IRqSl5dXWlqKda6wixcv/vLLL7feNg21dP/99wcFBd3Zd5XuWeYv2bdv37Jly3Q6pZ94XlhY6OHhgZ5llGTtd8pgMKhdiIYIIYqKivCDT2ElJSUuLi4uLg59EFIfvfDCCx07dryz7zp0EAIAANQ1XCMEAABNQxACAICmIQgBAEDTEIQAAKBpCELavHlzu3btvL29e/bseebMmUqfduvWjf1hzJgxqlTofKpfq9nZ2XFxcT4+Pv3798/OzlalQufDblHxU2zn9iXLcnh4eMUh1W/V2ObVpfUgTElJGTp06EcffZSamtq/f/+RI0dW/FQIkZiYeOXKlfz8/Pz8/Pfee0+tOp3Jbdfqm2++GRwcnJqaGhQUtGDBAlWKdD75FUyfPv3ll18u/wjbuX0tWrSoa9euSUlJFQdWv1Vjm1eZ0LadO3c+88wz1tfp6em+vr4VP01NTTUajTExMUajccCAAdanUkAt3XathoaGJiQkCCESEhJCQ0PVqNGZHT9+/P777y8rKysfgu3cvnbs2PHNN99U2rtWv1Vjm1eX1o8IY2NjP/roIyKSZXnGjBmDBw+u+On169c7d+788ccfX7p0ycvLa8KECSqV6VRuu1avXr0aHBxMRNbfyGrU6LTMZvMzzzyzePHiijd0Yzu3r169evXr16/SwOq3amzz6sIN9UREP/3005QpUx588ME5c+ZU1eNDampq27Zts7KyFK7Nudlcq56enpmZmW5ubkVFRSaTqbCwUK3ynM+8efOuXbu2ZMmSqkbAdm4vjN20d61+q8Y2ry6td/MjhHj11Vf37t27fv360NDQSp8ePny4pKSka9euRKTX69EBmF3cdq02adLk8uXLbdq0uXr1atOmTdWo0TnJsrxs2bLt27dXGo7tXAHVb9XY5tWl9VOj+/bt27hx45YtW5o0aVLxIek///wzERUWFg4cODAhIcFsNsfHxz/yyCNq1uosqlqr1nVORHFxccuXLxdCLF++fMCAAaoV6nR27NjRvHnz1q1blw/Bdq4Ym1s1tnlHoeYFSgcwZ84cmyvE+oJz/sEHH7Rq1crPz2/48OG5ubmqFuskqlqr5Ss/Ozu7b9++TZs2jYuLy8nJUa9SZ/Pkk0++/vrrFYdgO687lfauNrdqbPMOAtcIAQBA07R+ahQAADQOQQgAAJqGIAQAAE1DEAIAgKYhCAEAQNMQhAAAoGkIQgAA0DQEIQAAaBqCEAAANA1BCAAAmoYgBAAATUMQAlTGGEtMTFS7CrANfx2wOwQhQH2ydevW8oeFaYc2lxoUgyAEqE/i4uKuXLmidhVKq7jU+fn5tz5DG6A2tP6EegCoX4xGo9olgLPBESGorOIln8TERMZYxY8OHDgwePBgHx+fkJCQzz//vPwjIcSyZcsiIiI8PDw6dOiwcuXK8idrMsZ27949YMAAHx+f6OjoX3/99Zdffrn33nuNRmObNm1++OGH8hmdO3fuoYce8vLyat++/dq1a22Wxzl/77332rZtazQaY2JivvzyS+vw1157rWXLluUzTUlJkSRp8+bNNSzgtotgc8GtKyciIqLiWiqf2uLFiyMjIxs2bNijR489e/ZUX79dVtQdLAIRnThxok+fPj4+Pl5eXg8++GD5X//UqVOPPPJI06ZN3d3d27Vrt379eptLXXGDqX7pqioAoDLVHgkMIIQQgogSEhKsrxMSEipuk0R0zz33LFu2bN++fcOHD9fr9QUFBdaPVq1aFRkZuXLlyu+///6VV17R6XRLly4t/1ZUVNSaNWv27t3bu3dvLy+vtm3bfvXVV7/88ssDDzwQHBxcPqPQ0NC5c+du3rz5hRdeYIx98803t5b01ltvGY3GBQsWfPfdd9YZWUc7efIkEf3vf/+zjjZ37tyAgACz2VzDAm67CDYXPDU1lYh2796dmppaaTUuWrTIy8tr8eLFX3755cCBA3U63bFjx6qp3y4r6g4WwWKxNG7c+Kmnnvr888/Xr1//8MMPd+7cWQhRUlJiMpkiIyM//PDDzZs3P/fcczqdzvqs9kpLXZO/TvUbD0AlCEJQWfVBuHDhQuvrvLy8imN27Njx3Llz5WP++9//7t69e/m3vvjiC+vrAwcOENGRI0cqvi2f0YoVK8qn8OKLL/bs2bNSSZxzX19f64GO1SuvvNKjRw/r66ioqIkTJ1pHCw8Pnzx5cs0LuO0iVLXgFV+X45wHBgauXbvW+laW5T59+qxevbr6+mu/ou5gES5fvkxEp0+ftn6UkZGxatUqIUR2dvbMmTMPHTpkHZ6Tk1PVUtfwr1PNOgSoBEEIKqs+CA8ePGhzTE9Pz0rnNgICAspHK9/PWifIOa80feuLGzdulE98165dfn5+lWaUlpZGRJmZmeWj7d6929fX1/p63rx5TZs2lWX5t99+qzjTmhRw20WoasFt7tAzMjIqLY5V9fXXfkXdwSLIsjxy5Eij0di/f/8333wzJSWlfBzO+Z49ez744IPnnnvO2hym+iC87dJVtQ4BKkFjGXAgRUVFlYbcuqu18vDw+OWXX9zd3W1+WukS2q1X1G4lSZLFYqnJaLIsW18/8cQTU6dO3bt374YNG+69996IiIi/VED1i1DVgttkrVyn0/2l+mtY561TKF9Rd7AIkiQtX7583rx5GzZs2LVrV3x8/L/+9a8FCxYQ0bBhww4ePDho0KBHH300Pj7e39//tsXcOvGKS/eX1iFoGRrLgPrK8+/XX3+t4Vfatm179erV8PDw8PDwsLCwJUuWrFy58q/O95tvvil/vWnTpnbt2lUawWQy+fr6fvfdd+VDvv322/LRWrRo0bVr1zVr1qxbt27UqFF/de52WQSrgICARo0abd++3fqWc96pU6e5c+dWX3/NVbWi7mARsrOzR48e7ePjM27cuA0bNmzYsGHp0qVElJGR8dlnn+3atWv27Nl///vfxR+Nbqphr6UDwBEhqCw4OHjKlCmTJ0++fv36J598UsNvTZw48cknn5w9e3aLFi02bty4atWqijvrGpo8eXJaWlpkZOT27dsXL15sbfNZEWPs5ZdfHjt2rHW03bt3v/XWW5s2bSofYciQIRMmTHBzcxs0aNBfnfudLYJOp/vvf/9bVFTUsWPHinVOnjx59OjR169fb9269bp1606cOPHpp5/etv4aqmpF3cEiNGzYcMuWLYWFhYMHDy4pKVm+fHlMTAwRGY1GDw+P+Pj4IUOGXL58ecGCBTqd7rfffgsJCdHr9VUttV2WDgDXCEFl27dvDw8PNxqN999//4ULF+jma4QVr+tUevvJJ5+Eh4d7eHjExMR89dVXNkerdNGx0qWvAwcOdOvWzWg03nXXXeUtTSpNQZbld99913qHQHR09Jdfflmx+OvXr0uS9PTTT1ccWJMCar4Ild5OmjTJ09OzUaNGlVajLMvz589v3bq1h4dH586d//vf/962/tqvqDtbhF9//bVr167WpRg4cGD5ZcKNGze2atXKaDT26tXr8OHD//rXv4xGY3JycqWlruFfp/qNB6AiJmpwCgLAySQmJkZERNR+409OTg4JCdm9e3ePHj3sUpijsdeKAnBkODUKcCcsFktxcXF8fHx4eHj37t3VLgcA7tz/AxQrjNvfdykQAAAAAElFTkSuQmCC" }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "c_vals = range(2, 12, length = 25)\n", "\n", "models = [McCallModel(c = cval) for cval in c_vals]\n", "sols = solve_mccall_model.(models)\n", "w_bar_vals = [sol.w_bar for sol in sols]\n", "\n", "plot(c_vals,\n", " w_bar_vals,\n", " lw = 2,\n", " α = 0.7,\n", " xlabel = \"unemployment compensation\",\n", " ylabel = \"reservation wage\",\n", " label = \"w_bar as a function of c\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Note that we could’ve done the above in one pass (which would be important if, for example, the parameter space was quite large)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "w_bar_vals = [solve_mccall_model(McCallModel(c = cval)).w_bar for cval in c_vals]; # doesn't allocate new arrays for models and solutions" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Exercise 2\n", "\n", "Similar to above, we can plot $ \\bar w $ against $ \\gamma $ as follows" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVwV5f4H8O8zw3bgIAiyiYC5giKgad5LWnrzZlpHtNLyF2o300pJ85b7dq8LinZTs8XbQpL3qr/qZm7lrdyX+ytNEVEWdxEJUZaDbOcw8/z+OF0iPCDLYeYsn/fLP86ZmTPzdZjhw/PMc2YY55wAAAAclaB2AQAAAGpCEAIAgENDEAIAgENDEAIAgENDEAIAgENDEAIAgENDEAIAgENDEAIAgENDEAIAgENDEAIAgEOz4SC8dOlSaWmp2lVYL0mSZFlWuwrbxjmvrq5WuwqbV11djVs5thAOxVZlw0GYkJBw+PBhtauwXnq93mAwqF2Fbauuri4uLla7CptXWlpaVVWldhW2DYdiq7LhIAQAAGg5BCEAADg0BCEAADg0BCEAADg0BCEAADg0BCEAADg0BCEAADg0J7ULAAAA+1RYxc8VNf9eCu3cWLg3s2A99UEQAgCA5RkkvuyUfL6k+UE4IFBAEAIAgK1695x8voS3pFXX3cuyFdULQQgAABa2/Yq87wZ3E2lxH6GjpxKtupbAYBkAALCktEL+SbbMGE2PFK0/BQlBCAAAFnSzgq86LUucRt8nDAi0gRQkBCEAAFiKQaYVqXKJgff2Zc91sZl8sZlCAQDAmnFO69LlC3oe7MFmx4iCbbQGiRCEAABgEf+6LB/KkzUizYsRPWxqICaCEAAAWir1Nt90QWaMpvcSQ7VqV9NECEIAAGiRm5W0Ok2SOY3tLDwYYDtdov+FIAQAgOarlGjJT5LeQP392bOdbTJTbKofFwAsZ9euXQsXLlRgQ7IsM8YYs72GgvXgnMuyLIqi2oWYUWKgSomLjK65sn9Z4ofcvXv3rVu3WmBFjYYgBHBQt2/f7tChw5IlS9QuBOBX2dnZiYmJCm8UQQjguHx8fHr37q12FQC/EgQVOldtsj8XAADAUhCEAADg0BCEAADg0BCEAADg0BCEAADg0BCEAADg0BCEAADg0BCEAADg0BCEAGDDrOrObQsXLgwMDLTU2iz+X7NsefYEQQgAYBnvv/9+enp6C1cyZMgQ04vZs2e3uKLfsEh5dglBCABgGbdv327Xrl0LV7J3717Ti5UrV7a4ot+wSHl2CUEIAL/KucMv6RX6V83rbj02Nvbrr78motmzZwcEBHDOOecBAQFnzpxpoObFixf7+fnFxsZeuHDBNGXHjh0xMTHe3t5BQUFvvvmmaSJjbNOmTXX6Bs0umZycHBQU1K5du7fffrvOtswubzJy5EgiiomJod/2ata8Zoxt3bo1Ojra19d37dq1RHTr1q2nn37az8+ve/fu//rXv8yupLCwcPz48UFBQe3bt58wYUJhYWF9a6th9iO112yi1+snTpwYGBjYvXv3lJSUmjrr23vz5s0LCgpasmTJX/7yl27dunl5edXcHbvhuQ3sNGvBbdawYcN2796tdhXWq7CwsKKiQu0qbJvBYCgoKFC7itaycePG8ePH15n4wgHjE3sU+vdzed2Sli9fPmPGDM553759w8LCzpw5k5aWFhYWJstyff8LIlqwYIHRaJw7d+6wYcNME6OiolavXl1dXX3y5EkXF5eaJadOnZqenl7742aX9Pb2TktLy8rKGjJkSJ3NmV2+djF1XtSZmJSUJMvy3r17NRoN5zw+Pn7mzJlGo/HQoUNt27atrKy8eyXjxo2bOHFieXl5eXn5iy+++Pzzz9e3thoNfKT2YpMnT37hhRfKysoqKiomT55cM7e+vbd58+a0tDQiWrNmTVVV1bfffuvm5taYuQ3vtDpSU1OjoqIaXsbiGOd3/VVmI4YPH56QkDB8+HC1C7FSRUVFGo3Gzc1N7UJsmNFoLCkpsdfepJSUlH379qWkpNSe+NeT0u1KhQr4y/2ij+tvpqSlpT333HNHjhz53e9+FxcXFxQUxDm/ePHi+vXr61sJY+zixYudOnW6detW586dS0pKiEiW5ePHj589e/bgwYOffvqp6bccY+zmzZt+fn61P252ybi4OFmWx48fP2rUKCcnp3suX7uYmm3VzKo9Ua/Xe3p61kxs165dVlaWr68vERUWFnp5eYmiWGcl7dq1O3funL+/PxHl5+dHRUXl5+ebXVtNGQ18pPZipqZ2zWKBgYGmufXtvaqqKmdnZ0EQjEajk5MT51wQhMbMbXin1XH69Onx48efPn26gWUsTtHHMEmS1LNnz8zMzJop27dvnz9//vXr16Oioj766KNu3brVXr6oqGj8+PFHjx4dMGBASkpK27ZtlawWwAEt7qPmo1979eql1+s/++yzhx566JFHHlm/fj3nfNq0aY35rCAIzs7OptdjxoxxcXEZO3bsihUrPv3005pl6qRgfUt+9dVX33777caNGz/44IPvvvvunss3oLS0tPZbU27VkCSp5qlDeXl5bm5u7u7ud6+kdueqJEn1ra0xH6mz9ZrFaj/yt77/o4uLi+mF6Y+DOoNaG5jb1J2mAsXanmvXrn3ggQdqb/Hq1atarfbYsWPl5eWrV6+OjY2t85HZs2dPnTq1srJy6tSpc+bMqTMXXaMNQ9doyzlg16jqXnnllZCQkM2bN5eVlfn5+XXo0KGqqqqB5Ylo4cKFRqNx/vz5Tz/9tGlimzZt0tPTZVn++OOPichoNPK7OgYbWDIsLCwrK+vcuXNt27ZtzPK1izG9cHV13bdvnyzLSUlJVH9/6dNPPz137tzq6upjx45ptdqSkhLTLIPBULNMfHz8pEmTKioqTP2c48aNq29tNRrzEc75hAkTXnzxxfLy8oqKildeeaVmbsN7z+x2G57b8E6rQ5WuUeVahFFRUZ07d9bpdDVTLl269Oyzz/7+978nogkTJtw9RGrbtm3bt293dXVNSEiIi4tbsWJF7bkVFRUpKSlHjhyp86n+/fs/8sgjrfOfsCVlZWWSJFVXV6tdiA0zGo1lZWX22r1cc0XKqjzxxBPvv//+oEGD3N3dIyMjfX19a5oa9XF2dg4KCurZs+emTZtMUxITEwcNGhQQEDBhwoShQ4dOnDixTg9wDbNLzpo16/e//72Tk9Pq1asbs/zdq122bNlTTz0VHBz8xhtvNFD5+vXrJ02a5O/v7+3tvXHjxjZt2hDR8OHDO3fufO3aNdMya9eufe211zp27MgYe/TRR+uMizGrkR9Zs2bNtGnTQkJCAgIC5s6du3nz5ib9HxuvqSvknN+5c6epW3Fzc6vTj914Sl8jrNNJbSJJUkJCgiAI7777bu3pWq22oKBAo9FUVFQEBATo9fracwcPHuzv79+5c+c6a+vfv/8f/vCH1ijethQXF7u5udnrL3FlGI1GvV5vuoRjfzZv3nz06FEr7aqC1rdjx47o6OiwsDAiysrK0ul02dnZahdFp0+fjo+PP3bsWFM/qNFomh2Eil4jNOv777+fNWvWo48+umzZsjqzOOemjmbO+d3d3BqNZsKECRgsU5/q6moMlmkho9Eoy3IDV2Jsmpubm1XdlqVhZku1whatDTly5MiHH3748ccfG43G2bNnP/3002pX9AtBEBQ+6dQMQs75vHnzjh49unXr1jrDZEzat2+fk5PTtWvX3Nzc4OBg5SsEACuBzLO4RYsWvfLKK926ddNoNCNGjJg3b57aFalGzSA8duzYtm3b/u///s/JycnUI6zVaonowIEDgwYNIiKdTpecnJyYmJicnBwXF6diqQAAdkar1dZcVXVwagbhgQMHsrKyan8pwvRH3+DBg00vFi1a9Nxzz4WEhPTp0wc/MAAAaA1KB2Ht/o358+fPnz+/gWW8vb13796tUGUAAOCQcK9RAABwaAhCAABwaAhCAABwaAhCAABwaOp/oR4AVOHq6vrTTz818pbWYH+qZeJEjJGTNd1W4datW8pvFEEI4KCeffbZqqoq03OLWlV5ebmzs3PNoyGgGSRJqqys9PDwsNQKD+bJWSXk6USj7hPc1HziSF1dunSpfUtqZSAIARzXhAkTFNgKHo3ZcpZ9NObOa/J3GXKESKv7ix09rak9qBJcIwQAcCDninhypswYTYtECv4CQQgA4CgKKikxVarmNPo+YWAgUvAX6BoFAHAIBomvSJVLDNTblz3XBa2gX2FfAAA4hPfOyedLeHt3NjtGFNAarAVBCABg/7Zd4XtvcDeR5vcWPdAV+FsIQgAAO3f6Nk/Jlhij13qJoVq1q7E++MMAAMDafX+Dtl9002ik5n08t4xLnMZ2Fh4MQJeoGQhCAACrdq6Ib8jglUbBxcjvvXQ9+vuzsRggUw8EIQCA9bpVyVeeliVijwZWx3XXNHs9HTwYGoP1QRACAFgp0xceiqp4lA/FdzQEtEGWtQq0lAEArNSGDDm7hPu70es9SUQIthoEIQCANdpxVf4ul7sINK+32MZF7WrsGoIQAMDqnCvin2TJpi88dEaPaCtDEAIAWJeCSkpMlXFHUMUgCAEArIhBphWpUomB446gisFeBgCwIuvOSLgjqMIQhAAA1uLLK/Khn3FHUKUhCAEArELqbf5ptow7gioPQQgAoL6blbQ6TZI4PdsJdwRVGoIQAEBllRItPSnpDbgjqDqwxwEA1MQ5rUuXrpTyYHeaEYnxMSrA1VgAADV9flk+8jPXiDS/j5OHs9rVOCS0CAEAVHPqNv/nBZkRzYwWQzzUrsZRIQgBANSRW8aTUiWZU3xXoZ8f+kRVgyAEAFBBhUSJqVJZNf0+gI2+D7+K1YS9DwCgNM5p3Rnp2h0K0bLXIkU8M1ddCEIAAKVtuSgfzedaZ1rYW3THmEW1IQgBABT1w02+9aLMiN6IEoPc1a4GEIQAAEq6XsbXpMuc6E/dxfvboUvUKiAIAQAUUl5NialymZE/HMRGdUQKWgsEIQCAEjjRm2lyzh1+nyd7taeodjnwKwQhAIASNp2XjxfIXi5sQW/BFTloTRCEAACt7sjP8heXZZHRrGjBX4NOUeuCIAQAaF1XSmldusw5TQwXonyQglYHX2ABAGhFpUaemCpXSjS4PdOFou1hjfBTAQBoLTKnv53heeW8cxuWgAEy1gpBCADQWj7Jln8qkL1daEFvwQW/bq0VukYBAFrFgTz5qyuyk8BmRwvt3HBp0HrhTxQAAMu7pOfvnJWJ6KUIIRIDZKwbghAAwMJKjZSYKlVJ9Eh79lgHpKC1QxACAFiSxGlFqpRfQeHebCoGyNgCBCEAgCV9lCmfKeQ+rmxOtOiMX7G2AD8lAACL2X+D77omOws0L0bwdVO7GmgcBCEAgGVc1PN3zkpE9HKE2N0blwZtBoIQAMACig207JRskEkXJjyKATI2Bd8jBABoqWqZJ52Wb1XyCG/2Qnc0MGwMfmAAAC3190yeXsj93Gh+b9EJrUFbgxYhAECLfH1N3pMju4hsbozg5aJ2NdB0aBECADRfRjH/MFMmolcihK5eaAzaJAQhAEAzFVbxlalSNadRHYUhwUhBW4UgBABoDoPEl5+SC6so2pdN6IbfpTYMPzwAgObYkCFnl3B/DZsVJYpoDdoyBCEAQJPtuCp/l8tdBJoXI7TBABkbhyAEAGiac0X8kyyZMXqtl9i5DRqDNg9BCADQBAWVlJgqV3N6+j5hYCBS0B4gCAEAGssg04pUqcTAe/uy+C74/Wkn8IV6AHAgZdWUV86b/fFtV/j5Et7enc2OEQW0Bu0FghAAHEU1p7/+JGUUNz8IiUgj0oI+ogd+d9oRRX+YkiT17NkzMzPznhNNHnzwwWPHjplev/TSSxs2bFCiSgCwUxvOyRnF3NOZBWiav5JnOgshHparCayAckG4bt26zZs3Z2Vl3XOiCec8MzPz+vXrXl5eROTkhD/AAKD5vs7h/74uu4jsr/fjXmjwG2Yu9nLO33///YEDBwYGBmZmZs6dO3fXrl0t31JUVNTChQsbM9EkPz/fYDDExcUFBQXFx8fr9fqW1wAAjimjmH+YIRHRlB5IQajLTDMrJSXl1VdfnT59+pEjR4hIr9frdLovv/xy1KhRLdnS4MGDGznR5Oeff+7Xr99bb70VGho6Y8aM6dOnb9mypfYCJSUlo0ePdnZ2rvPBF154Yfbs2S0p1T6UlJS4ubm5urqqXYgNMxqNpaWlkiSpXYhtU/1QLDayv6S7lxnY8CBDpGjIz1erkObDoXhP3t7ezT7GGOd1rxtHR0c//vjjiYmJjLGMjIzu3btPmDAhMzPzxx9/bHGpxJiZLZqdWFteXl7Pnj0LCwtrTxw2bNjzzz8/ZMiQOgu7ublpNC24AmAvioqKNBqNm5ub2oXYMKPRWFJS0q5dO7ULsW3qHooGic87wbNLeLQv+0sfwUbvhYZD8Z4EofnfZjHTIszOzl61alXNW8bYqFGjxo0b1+xtNM/JkycrKytjY2OJyMXF5e6oZ4x5enr6+voqXJitEP5L7UJsGPahRai7Gz84J5/X8wANmx3t5CyqUoIF4FBsVWZ2a1hYWEFBQe0pV65cCQ0NVaokOnDgABGVlZWNGjUqIyPDYDAsXbp05MiRihUAAPZh+xX5u1zZTaSFfUTcERTqYyYIp0+fvnLlyosXLxJRRUXFrl27Fi9ePHPmTMVqMl04HDBgwOLFi3U6XXBwcFFRUVJSkmIFAIAdOFfEN57njNG0SLGjp212iYIizHSNvvzyy1qtdsyYMR4eHkOHDu3Zs+c//vGPESNGWGR7Zq8F1ploessYmzJlypQpUyyyXQBwKDcreGKqXC3zMZ1wR1C4BzNByBgbN26c8hcFAQAswiDTilTZdEfQ53BHULgXHCIAYG/WnZEu6HFHUGgsM0HI7uLs7BwUFNS3b9/33nvPYDAoXyUAQCP967J86GeuEWl+b9wRFBrFzGFy+PDhsWPHjh8/ftiwYYIgfPPNN1u2bPnggw+uXbv2l7/85dq1aytXrlS+UACAe0q9zT89LzNG03uJoVq1qwEbYSYI33rrrQkTJixbtsz0NjY2lnP+ySefbNq0KSws7IknnkhMTMTXWQDA2tyspNVpkszpf7oIDwagSxQay0ye7d+/f8CAAbWnPPjgg3v27CGifv36lZeX5+bmKlQdAEDjVEq09KSkN1B/f/ZsZ/ylDk1g5nAJDAxMS0urPeXMmTPe3t5ElJOTQ0Q+Pj7KFAcA0Bic07p06UopD3anGZEYHwNNY6ZrNCEhYc6cOYGBgY8//jhjbPfu3UuXLl2yZMnly5enTJny2GOPeXjgYVwAYEU+uywf+Zm7O9H83qKHM3IQmsZMEE6ZMoUxNnfu3AkTJhCRn5/f0qVLp02btnnzZn9//3feeUfxIgEA6nXyFt98QWZEb0QJIVqkIDSZ+S/Um27pcvv27erqan9/f8YYEcXHx8fHxyteIQBAvXLL+KrTksxpfFehnx8uDUJzNHTc+Pr6BgQEMMYMBsOVK1eUKgkAoFEqJEpMlcqq6fcB7On7kILQTOa/blpSUpKXl1fz9scff0xISMAz4gHAenBOa89I1+5QiJa9Fiky9IlCc5kJwi1btsTHx8uyXDNFFMXJkycrWBUAwD1svigfy+eezmxhb8Edd5CBFjDTmbB06dJJkyYVFxd369bt5s2b+fn5UVFRCEIAsB4/3OT/e1EWGL0eJQS5ozEILWImCC9dujR8+HAvL69BgwadPn3a399/5syZr7/+uvLFAQDc7XoZX5Muc6Lnuwn3t0MKQkuZCUIfH5+bN28SUVRU1KFDh4goJCTk+PHjSpcGAHCXsmpaelIuM/KHg9iojhggAxZg5jDq16/fmjVrTp48GRMT89lnn+Xn52/ZsiUgIED54gAAauNEb6ZJN8p55zZsWqSodjlgJ8wE4fLly2/evLljx47Y2Ni+ffsGBgZu3LjxzTffVL44AIDaNp2XTxRwT2eaGyO6oDUIFmJmrFVkZOT169crKioYY//4xz/WrFmj1Wo1Go3yxQEA1PhPPv/isiwymhsjBuAXEliOmSCsrq52dXV1dXU1vfXz81O2JACAuq6U0ltnJM7pxQihlw8GyIAlmelc8Pf3j4+P/9///d/i4mLlCwIAqKPUyBNTpUqJBrdnT4SiSxQszMwhlZCQcPny5fj4eD8/v0ceeWTt2rUXLlxQvjIAACKSOf0tTc4r553bsISeGCADlmcmCJcsWXL06NHCwsKvvvoqJiZm48aNXbt27dGjh/LFAQBszJZ/usW9XWhBbwEDZKA1NHRjIvZfRJSfn69USQAAvziYx7ddkZ0ENjtaaOeGS4PQKsz8fbV48eIBAwb4+Pg89dRTp0+ffvbZZ0+cOFFQUKB8cQDgyC7p+fqzEhG9FM4iMUAGWo2ZFuGSJUucnJwmTpw4b9680NBQ5WsCACg1UmKqVCXRsBDhsRB0iUIrMnN4bd68+YUXXti/f39YWFjPnj0TEhK++OIL003XAAAUIHFamSrlV1C4N5sUjhSE1mWmRTh27NixY8cS0Y0bNw4ePLhr165nnnlGlmXOueLlAYAj+jhLTivkPq5sbozojByEVmZ+sIwkSSdPnvzuu+++/fbbY8eOeXh4DB06VOHKAMAx7b/Bd16VXUQ2v7fg46p2NeAAzATh6NGj9+7dW1RU1KlTJ51Ot2DBgoceesjFxUX54gDA0VzU83fOSkT0coTQzQsDZEAJZoIwPz9/zpw5Op0uPDzc9N0JAAAFFBto2SnZIJMuVPhjMH75gELMBKHpGYQAAEqqlnnSaflWJe/Rlr2AATKgIBxtAGAV/p4hpxdyPzeaFyM6oTUICmrozjIAAMrYe4Pvuc5dRDY3RvDCgARQFlqEAKCyzGL+7lmJiKb0ELpigAwoDkEIAGoqrOIrT0tGmUZ1ZI+0RwqCCsx3jZaVleXk5NSZGB4e3vr1AIADMcqUmCrfrqRoXzahGx6xBOowE4Rbt26dMGGCwWCoMx13lgEAy9qQIWUVc383mhUlimgNgkrMdI3OnTt37NixpaWl/LeULw4A7NjOq/K317mLQPN6i20wQAbUYyYIS0pKXn75Za1Wq3w1AOAgMop5cpbMGL3WS+zcBo1BUJOZIBwyZEhaWprypQCAgyiopOWnpGpOT3UUBgYiBUFlZq4RvvHGG9OmTauoqBgwYICHh0fNdAyWAYCWM8i0/JRUYqD727Hx3TBwHdRnJgj79+9PRD/88EOd6bhMCAAt995Z6aKeB2jo9SiMjwGrYCYIEXgA0Eq2XeF7b3A3kRb2ET2d1a4GgIjwhXoAUMzp2zwlW2KMpkcKYVq0BsFamA/CnTt3Dhw4sF27dj4+PgMHDty9e7fCZQGAnblZSavSJInTs52EAYH4ExysiJnD8fPPP3/yyScHDhz41VdfmRIxLi7uyy+/VL44ALAPBpkST0l6A/X3Z2O7IAXBupi5RpiYmDhr1qzly5eb3j744IOyLC9fvvzJJ59UtjYAsAec0zsZ7KKeB7vTjEgBXaJgbcz8aZadnT1gwIDaUx5++OGsrCylSgIAu7I91+noTaYRaX4fJw9n5CBYHTNBGBYWdvbs2dpT0tPTw8LClCoJAOzHqdv882tOjGhmtBjice/lAZRnpmt08uTJCxcuDAgIePzxx4lo9+7dS5YsWbp0qeK1AYBtu1HOk1IlmdNznXg/P7QFwUqZCcJp06ZVV1fPmDFj/PjxROTr67to0aJp06YpXhsA2LBKiZafksqqqa+P9GQoUhCsl5kgFAThjTfeeP311wsKCojIz8+PMRzEANAEnNPaM9K1OxSiZVO6VTNm/tGnANag3qOTMebv769kKQBgN7ZclI/mc60zLewtulXhZlVg1X4JwszMTCLq1q2bIAim13fDTbcBoDF+uMm3XpQZ0RtRYpA7FVWpXRBAg34JwoiICCIqLS3VarWm13fDPUgB4J5yy2lNusyJ/tRduL8drqqADfglCGuHHAIPAJqnvJqWn6wuM1JsABsVhjvIgG0wc6Qyxn7++efaU06ePOnj46NUSQBgkzjRm2lyThnd58lm9BIxxg5sxa+DZe7cuXP9+nXT6/PnzxcXF9fMOnTokNFoVLo0ALApm87LxwtkT2c2L0ZwE9WuBqDRfg3CAwcO6HQ60+uHHnqo9kKiKM6YMUPRugDApvwnn39xWRYZzYkRAt3RGARb8msQPvHEE6arg4yxvLy8wMBA9aoCAFuSc4evTZc4p4kRQpQPUhBsjJlrhJzzOiloMBiuXLmiUEUAYFNKjXzpKbm8mga3Z7pQDJAB22P+C/UlJSV5eXk1b3/88ceEhAS9Xq9UVQBgG2ROf0uT88p5J09K6IkLg2CTzAThli1b4uPjZVmumSKK4uTJkxWsCgBsw8Zs+adb3NuFFvYRXdAaBNtk5shdunTppEmTiouLu3XrdvPmzfz8/KioKAQhANRxMI9vuyI7CWx2jNjODZcGwVaZCcJLly4NHz7cy8tr0KBBp0+f9vf3nzlz5uuvv658cQBgtS7p+fqzEhFNDmeRbZGCYMPMBKGPj8/NmzeJKCoq6tChQ0QUEhJy/Pjxlm9MkqS7b1hqdqJJUVGRTqfz8fEZMWJEUVFRywsAAIsoMfDlqXKVRI92YMNC0CUKts3MEdyvX781a9acPHkyJibms88+y8/P37JlS0BAQAu3tG7dutjY2KysrHtOrJGUlBQWFpaXlxcaGrpq1aoWFgAAFiFxSjot36zg4d7s5QgMkAGbZ2awzPLlywcPHrxjx47Fixf37ds3MDDQ3d198+bNLdxSVFRU586da76z38DEGtu2bdu+fburq2tCQkJcXNyKFStqzzUYDP/+979v3LhR51MRERG9e/duYbV2oKKignNee9ATNJXRaKyoqCgvL1e7EOuSfEFILWDeLjS9m2ys5Pe86RQOxZbDoXhPrq6uotjMP8uY2VtsV1VVVVRUeHt7E1FBQYFWq9VoNC2qsWZ7zMwWzU4kIq1WW1BQoNFoKioqAgIC6nx/4+GHHyai9u3b1/nUkCFDnn76aYtUa9OKi4vd3Nzc3NzULlgH67gAABuWSURBVMSGGY1GvV7v6+urdiFW5PBN8d1sJ5HRwkhDuFejbtCPQ7HlcCjek7u7u7Ozc/M+a6ZFuH79+vj4+LZt25re+vn5Nb+0luGcM8ZMLyRJqjPXw8MjISFh+PDhapRmA2RZ1mg0+O3TEqZb7Hp5ealdiLW4qOefXJGcnOjVnmL/Do09tHAothwOxVZl5hrhG2+8ERQUNG7cuEOHDqn7SKb27dvn5OQQUW5ubnBwsIqVAECxgZadkg0yPREqPNoBw0TBfpgJwry8vLfffvvatWsPP/xweHj4m2++aRpEqpgDBw6YXuh0uuTkZM55cnJyXFyckjUAQG3VMk86Ld2q5BHebGI4homCXTH/9YnJkycfPHjw6tWrf/rTn1JSUjp06DBmzBjFaho8eLDpxaJFi9LS0kJCQs6ePbtgwQLFCgCAOv6eIacX8nZuNL+36ITWINgX8/caNQkNDR05cqTBYFi/fv3nn39uke2Z7WutM7Hmrbe39+7duy2yXQBotr03+J7r3EVk82IELxe1qwGwNPNdHBcuXEhMTIyJiYmIiNi0aVNCQsL58+cVrgwArEFmMX/3rEREr0QIXb3QGAQ7ZKZF2Ldv359++snX1/eZZ57ZsGFD//79TUM3AcDRFFbxlaclo0wjOwpDgvF7AOyTmSDs1KnTokWLHnvsMRcXdIIAOC6DxJefkm9XUrQve74bBsiA3TJzcF+9erWsrAwpCODgNmTI2SXcX8NmRgkiWoNgv8wEob+//759+5QvBQCsx85r8ne53EWguTGClwtiEOyZma7RRYsWTZkyZdasWTqdrvZtZep7RgQA2JlzRTw5U2aMXusldmmDFAQ7ZyYIH3jgASI6ceLE6tWra09X9y4zAKCMgkpKTJWqOY3uJAwMRAqC/TMThAg8AIdlkPiKVLnEQDG+LL4LBsiAQ8CBDgC/eu+cfL6EB2hoVrQooDUIjsFMEHLO33///YEDBwYGBmZmZs6dO3fXrl3KVwYACtt2he+9wd1EWthH9GzmA20AbI+ZIExJSXn11VcfeOCB/Px8ItLr9Tqdbtu2bYrXBgDKOX2bp2RLjNH0SCFMi8YgOBAzQbhmzZpZs2b97W9/M7195513xo0bV+fp8ABgT25W0qo0SeL0bCdhQCCumIBjMXPEZ2dnmx7+bsIYGzVq1Llz5xSsCgCUUynR0pOS3kD9/dlYDJABx2PmoA8LCysoKKg95cqVK6GhoUqVBADK4ZzeTpeulPJgd5oRifEx4IjMBOH06dNXrlx58eJFIqqoqNi1a9fixYtnzpypeG0A0Oo+vywf/plrRJrfx8nDGTkIjsjM9whffvllrVY7ZswYDw+PoUOH9uzZ8x//+MeIESOULw4AWtWp2/yfF2RGNDNaDPFQuxoAlZgJQsbYuHHjxo0bp3w1AKCY3DKelCrJnMZ1Ffr5oS0IjgsXxgEcUYVEy09JZdUUG8BG34ffA+DQzJwAGRkZ/fr127RpExHNnDlTq9UOHjw4JydH8doAoFVwTuvOSDllFKJl0yNFPHgbHJyZIHz11VcDAgKGDRt27dq1DRs2fPHFFy4uLjNmzFC+OABoDVsuykfzudaZFvYW3c1cHgFwLGZOgh9++OGjjz5q167dqlWrhg0b9thjj5WUlEydOlX54gDA4n64ybdelBnRG1FikLva1QBYATMtQicnJ8YYER09etT0zXqNRmMwGJQuDQAsLbec1qTLnOhP3cX726FLFIDIbBD2799/586dGRkZ+/btGz16tNFo3Lp1a3R0tPLFAYAFlVfT8lNSmZHHBrCRYUhBgF+YCcJVq1Z9//33PXr0eOaZZ/z9/adOnbpv3746D+kFANvCid5Mk3Pu8Ps82Z97YYAMwK/MXCOMiorKycnJz89v3749ESUlJb3//vuiKCpeGwBYzKbz8vEC2dOZ5sUIrjibAWoxP2LMyckpODjY9Lpt27YK1gMAlveffP7FZVlkNCdGDHRHYxDgN/BgXgA7d6WU3jojcU4Tw4UoH6QgQF14MC+APSs18sRUqVKiwe2ZLhR3kAEwAw/mBbBbMqe/pcl55byTJyX0xIVBAPPwYF4Au7UxW/7pFvd2oYV9RBe0BgHqYWawDB7MC2AHDubxbVdkJ4HNjhHaueHSIEC98GBeADt0Sc/Xn5WIaHI4i2yLFARoCB7MC2BvSo2UmCpVSfRIsDAsBF2iAPeAB/MC2BWJ08pUKb+Cwr3Z1B5IQYB7M3+v0S1btihfCgC03MdZclohb+vK5kSLzshBgEYwc6L4+/vv27dP+VIAoIX23+A7r8pOApsTLfi6qV0NgI0w0zW6aNGiKVOmzJo1S6fT+fn51UwPDw9XsDAAaJqLev7OWYmIXokQemCADECjmQnCBx54gIhOnDhR54kTnHOFigKAJio20LJTskGmJ0KFRzsgBQGawEwQIvAAbEu1zJNOy7cqeYQ3mxiOC4MATYNzBsDm/T1DTi/k7dxofm/RCa1BgCYy/xgmALAV3+TIe65zF5HNixG8XNSuBsAGoUUIYMMyi/mHmTIRvRIhdPVCYxCgORCEALaqsIqvSJWMMo3sKAwJRgoCNBOCEMAmGSS+/JRcWEXRvuz5bjiRAZoP5w+ATdqQIWeXcH8NmxkliGgNArQAghDA9uy8Jn+Xy10EmhsjeLkgBgFaBEEIYGPOFfHkTJkxmh4pdGmDFARoKQQhgC0pqKTEVKma01MdhYeCcP4CWAC+RwhgMwwSX5EqlxgoxpeN64oUBLAMnEsANuO9c/L5Eh6goVnRooA+UQALQRAC2IZtV/jeG9xNpIV9RE9ntasBsCMIQgAbcPo2T8mWTANkwrRoDAJYEoIQwNrdrKRVaZLE6dlOwoBAnLMAFoaTCsCqGWRKPCXpDdTfn43tghMWwPJwXgFYL85p7Rnpop4Hu9OMSIyPAWgV+PoEgPX6/LJ8+GeuEWl+HycPDJABaB1oEQJYqVO3+T8vyIxoZrQY4qF2NQD2C0EIYI1ulPOkVEnmFN9V6OeHPlGAVoQgBLA6lRItPyWVVdPv/Nno+3CSArQunGMA1sU0QObaHQrRshm9RIbWIEArQxACWJctF+Wj+VzrTAt7i+4YzQbQ+hCEAFbkh5t860WZEb0RJQa5q10NgGNAEAJYi9xyWpMuc6Lnuwv3t0OXKIBCEIQAVqG8mpafksqMPDaAjQrDiQmgHJxvAOrjRG+myTl3+H2eGCADoDRFg1CSpPDw8NpTioqKdDqdj4/PiBEjioqK6iz/4IMPsv96+eWXFawUQFGbzsvHC2RPZzYvRnAT1a4GwMEoF4Tr1q2LjY3NysqqPTEpKSksLCwvLy80NHTVqlW1Z3HOMzMzr1+/XlpaWlpaunbtWsVKBVDSf/L5F5dlkdGcGCHQHY1BAKUpF4RRUVELFy6sM3Hbtm0JCQmurq4JCQlffvll7Vn5+fkGgyEuLi4oKCg+Pl6v1ytWKoBicu7wtekS5/RCdyHKBykIoALGOVd0e+w3W9RqtQUFBRqNpqKiIiAgoHbapaam/vnPf37rrbdCQ0NnzJhhMBi2bNlSe1UDBw7Mzc319PSss4kxY8a89NJLrfq/sAnFxcUajcbV1VXtQmyY0WgsLS318fFppfXfqaZFZzT5lWyAn/RKl6pW2orqcCi2XGsfinagTZs2Li4uzfusyt/X5ZwzxkwvJEmqPSsmJmbfvn2m1ytXruzZs2edz7q7u7/66quxsbF1pgcEBHh5ebVayTZDlmWNRuPm5qZ2ITbMaDQSUSsdTjKnt09ToUTd2vIZMc6uot3+pHAotlyrHor2QRSbf3Vd5SBs3759Tk5O165dc3Nzg4ODa886efJkZWWlKedcXFzu/nNSFMXu3bv3799fuXJtivN/qV2IbWu9fZicJZ8qlNu60uK+zlq7zggcihaBfdh6VP76hE6nS05O5pwnJyfHxcWZJh44cICIysrKRo0alZGRYTAYli5dOnLkSDULBbCog3l82xXZSWCzY8R2dp2CANZP5SBctGhRWlpaSEjI2bNnFyxYYJo4ePBgIhowYMDixYt1Ol1wcHBRUVFSUpKqlQJYzCU9X39WIqJJ3VlkWwyQAVCZ0l2jdcbmeHt779692+wyjLEpU6ZMmTJFueIAWl+pkRJTpSqJHgkWhofijhYA6sN5CKAcidPKVCm/grp7sak9cPYBWAWcigDK+ThLTivkbV3Z3BjRGScfgHXAuQigkP03+M6rspPAZkcLvhggA2A1EIQASrio5++clYjo5QjWEwNkAKwJghCg1RUbaNkp2SDT46HC0A446QCsi8pfqAewe9UyTzot36rkEd7sxXCkIIDVwWkJ0Lr+niGnF/J2bjS/t+iEPlEA64MWIUAr2nuD77nOXUQ2L0bwauYNgQGgdaFFCNBaMov5u2clInolQujqhcYggJVCEAK0isIqvvK0ZJRpZEdhSDBSEMB6IQgBLM8oU2KqfLuSon3Z891wlgFYNZyiAJa3IUPKKub+GjYzShDRGgSwbghCAAvbeVX+9jp3EWhujODlghgEsHYIQgBLOlfEk7Nkxmh6pNClDVIQwAYgCAEspqCSElOlak5PdRQeCsLJBWAb8D1CAMswSHxFqlxioBhfNq4rUhDAZuB0BbCM987J50t4gIZmRYsC+kQBbAeCEMACtl3he29wN5EW9hE9ndWuBgCaAkEI0FJ55TwlW2KMZvQSwrRoDALYGAQhQEt9m8slToODWGwATigA24PzFqBFqjl9f10mouEhOJsAbBJOXYAW+b98Xmygjp7U3RudogA2CUEI0CLf5MhENAzPnQewWTh7AZovr5yfKeJuIg1qj1MJwFbh7AVovm9yZM7poSDmjltTANgsBCFAM1Vz2neDE9Fj6BcFsGU4gQGa6Vg+LzHQfZ4MT58HsGkIQoBm+uaa6VsTSEEA24YgBGiO3HI6W8zdRMJTJgBsHc5hgOb45prEOQ3CMBkA24cgBGgyo0wH8oiIHsPdZABsH05jgCY78rNcYuBdvVhnPIMewPYhCAGabM91TkTD0BwEsAs4kwGaJucOzyjmGpEGBKI5CGAPEIQATfNNDuecBrdnGlHtUgDAEhCEAE1gkOlAHieioSGIQQA7gSAEaILDeXKpkXf3Yp081S4FACwEQQjQBL8MkwnFiQNgP3A+AzTWlVKeWcw9nNmAAAyTAbAfCEKAxjI1B//Qnrni+iCAHUEQAjSKQaaDeTIRDe2A5iCAXUEQAjTKwTz5jpEivFmYFkEIYFcQhACNsicHd5MBsE84qwHu7UopZZdwDyeKxTAZALuDIAS4t69zZCIa0kHAMBkA+4MgBLiHSumXYTKPBuN8AbBDOLEB7uFgnlxeTZE+LFSrdikA0AoQhAD3YBom81gHnCwA9gnnNkBDrpYJF/Tc05l+7692KQDQOhCEAA357mcnIhoSLLiIGC8KYJ8QhAD1qpToP7dEIvpjMFIQwG4hCAHqdSCPV8osyoeF4G4yAPYLQQhQr+9uMCJ6DHeTAbBrOMMBzMsu4ZdKydOZ/84fzUEAe4YgBDDP9K2JQf6SM84SALuGUxzAjDIjP/yzTMQH+RvVrgUAWheCEMCMA3m8UqJoHxboxtWuBQBaF4IQwIxvczkRPRqsdh0A0PoQhAB1ZRbzS3ru7UIPtENzEMD+IQgB6tpznRPRHzsITgLGiwLYPwQhwG+UVdORn2XG6I946BKAY8CpDvAb+27IVRL19mVB7mqXAgCKQBAC/Ma/8dAlAAeDsx3gV+eK+NU7vK0r6+endikAoBQEIcCvTMNkHg1mGCYD4DgQhAC/KKumY/kyI3oUd9kGcCSKnvCSJIWHh9eeUlRUpNPpfHx8RowYUVRUVGf5hucCWNb3ubxKovv9mL+b2qUAgIKUC8J169bFxsZmZWXVnpiUlBQWFpaXlxcaGrpq1ao6H2l4LoBl/TtHIgyTAXA8jHOF7p2xf//+srIynU5Xe4vdu3ffvn17eHh4ZmZmXFxcnZhseO7QoUMHDRrUv3//OhsKDg7u2LFjnYl7ctkPBZb871i/iooKZ2dnJycntQuxDUaZzhWTjyv/eyyJ/70+aDQa9Xq9r6+vqqXZvKKiIo1G4+aGhnbz4VC8J2dnZ0Fo5l+xygXhL9tjv9miVqstKCjQaDQVFRUBAQF6vb72wg3PHThw4K1bt7y9vets4sknn5w4cWKdiZ9edt6T51iRYDQaRVFs9pHhmEZ1MI4Ora55i98+FlFcXKzRaFxdXdUuxIbhULwnrVbr4uLSvM+qnA2cc8aY6YUkSU2a6+npOXfu3OHDhzdmQ6PdaFBHx7pvpF5vcHNzavaR4Zi6tNFonX99a/pjwsfHR72K7AFjDC3CFsKh2KpUDsL27dvn5OR07do1Nzc3OLjurf4bntskwe4U7O5YA+KLBFmj4W5ujvW/BgBoKpX7zXQ6XXJyMuc8OTk5Li7ONPHAgQMNzAUAALAglYNw0aJFaWlpISEhZ8+eXbBggWni4MGDG5gLjVRUVFRZWal2FbbNaDTevn1b7SpsHg7FlsOh2KqUDsI6Y3O8vb137959/fr1HTt2eHl51VnG7FxopKlTp37//fdqV2Hb0tPTR48erXYVNm/KlCl79+5VuwrblpaW9swzz6hdhd3CkEIAAHBoCEIAAHBoCEIAAHBoNvwdc4PB8N133+Xn56tdiJW6fv36/v37S0tL1S7Ehl2+fLmoqOiTTz5RuxDblpubu3///jo3xIAmuXjxYmFhIQ7FBjzyyCOhoaHN+6zSd5axoBMnTrz99tuiKKpdiJWqqKhwcXHB/mkJWZYrKyvd3fGs+hbBodhyOBTv6dVXX+3Tp0/zPmvDQQgAANByuEYIAAAODUEIAAAODUEIAAAODUEIAAAODUFoJ4qKinQ6nY+Pz4gRI4qKiurM3b59e2RkpLe390MPPZSdna1Khdav4X1okp6e7uHhoXBhtqXh3VhdXT1lyhQ/P78HH3wwNzdXlQqtX8P78ODBgzExMZ6enjExMYcOHVKlQjuDILQTSUlJYWFheXl5oaGhq1atqj3r2rVr8fHxH374YV5e3ogRI/70pz+pVaSVa2AfmpSUlDz//PPl5eXK12ZDGt6Na9eu1ev1V69ejY2NXbx4sSoVWr+G92F8fPz8+fMLCwvnzZsXHx+vSoX2hoNd6NatW0ZGBuc8IyOjW7dutWft37//xRdfNL2+efOmr6+vCvXZggb2IedcluWRI0d+/vnnOGsa1vBu7N27d2pqKudcr9efOHFChfpsQcP7sEePHh9++GFhYeFHH30UERGhRoH2Bt8jtBNarbagoECj0VRUVAQEBJi9i4ckSQkJCYIgvPvuu8pXaP0a3ocrV668devWm2++yRjOmoY0vBt9fX0nTZr0wQcfdOrU6ZNPPunVq5dadVqzhvfhiRMn+vXrZ3p9/Pjxvn37qlGjXUHXqJ3gnDPGTC8kSbp7ge+//75fv35eXl7r1q1TvDrb0MA+3L9//549e1asWKFSabak4UNRr9dzzs+ePfvYY49NmjRJjQJtQMP7cPbs2bNmzbpx48bMmTPnzJmjRoF2R6WWKFhYly5dsrOzOefZ2dldu3atPUuW5Tlz5gwcODArK0ul6mxDA/tw/vz5dU6cw4cPq1SmtWtgN3LOg4KCbty4wTnPy8vz8PBQoT5b0PA+9PDwyMvL45zfunVLq9WqUJ/dQYvQTuh0uuTkZM55cnJyXFycaeKBAweI6NixY9u2bduxY0f79u3v3Llz584dNQu1Yg3sw2XLltWcM0TEOR8wYICKpVqzBnYjEQ0dOnTjxo1VVVUffPAB+vTq0/A+jIqK+vjjj+/cufPpp59GR0erVqU9USd/wdKKioqGDx8eHBys0+mKi4tNE00/32XLluGH3hgN7MPasAMb1vBuzMvLGzJkiJeX10MPPXT+/Hn1yrRqDe/DjIyM2NhYrVYbGxtrGlMDLYTL/gAA4NDQNQoAAA4NQQgAAA4NQQgAAA4NQQgAAA4NQQgAAA4NQQgAAA4NQQgAAA4NQQgAAA4NQQgAAA4NQQgAAA4NQQgAAA4NQQjQUoyxzMzMhpfJzMw0PWGuJbZv3x4UFOTu7l7nNQC0BIIQoKXi4uK0Wq0CG3rttdeGDx9++PDhOq9b265du/D0LrBjePoEgBIyMzMjIiJaeLoxxs6cORMZGVnndWtjjGVkZISHhyuwLQDloUUIYEucnJzMvgaAZkMQArRU7WuEsiyvXbu2Z8+eWq32/vvv/+KLL2ovmZaW9vDDD7dp0yYyMvKTTz4x20Csbw2mS4wRERGMsdqviYhzvmHDhoiICHd39969e6ekpNSsmTF2/PhxnU7Xv3//u8uuPevs2bMjR44MDg7WaDSRkZFbt269e7sNbwvAVqn2SGAAe0FENQ8KX716tVarXbVq1ddffz1nzhxRFHfu3Mk5z8jIIKIOHTqsWLFi+/btCQkJjLHPPvvs7rXVt4a8vDwiOnToUF5eXu3XnPNPP/20R48eKSkp33zzjekj77//fk1tffv2nThx4j//+c+7y66ZVVlZ6efn16NHjw8++GD79u0vvfSSKIqmZ6M3flsANgpBCNBSNUEoy7Kvr6+pkWQyZ86cgQMH8v8GoakVaDJ9+vTf/e53dVbVwBr4bxO39us+ffpcuHCh5iN//vOfBwwYULPYzJkz6yu7ZlZRUdHixYt/+ukn09vi4uJmbAvARuEaA4DFFBQU3L59+4knnqiZMnz48A8//LDmrU6nq3n95JNPbty4salrMCsrK6tLly61pwQEBJjdaB01s7y9vRcvXnzs2LH33nsvLS1t//79zdsWgC1CEAK0IkEQJEmqb5aLi0tL1lDD3d398OHDGo3G7Fw/P7/6Plh71rhx406cODFmzJinnnpq6dKl/v7+zdgWgC3CYBkAi/Hz8/P19f36669rpuzevbv2Nxx27dpV83rbtm29evVq6hrM6tmzZ25ubnh4eHh4ePfu3devX5+SktKkygsKCv75z38ePHhwyZIlf/zjH3n9419avi0Aa4MWIYDFMMZmz549ZcqU/Pz8Hj16HDp0aPXq1V999VXNAn/961/z8/PDw8O///77d999d8+ePU1dg1kzZsz4n//5nyVLlnTs2HHbtm2ffvrpzp07m1S5Vqt1d3dfunTp2LFjc3JyVq1aJYrijz/+2KlTJxcXF1EU//3vf5eXl/fp06fl2wKwOmpfpASweVRrLIkkSW+99Zbp2wUxMTFffPGFabppsEx6evrDDz/s6ekZHR395Zdfml1bfWvg9Q+W4Zx//PHH4eHh7u7u999//7/+9S+zH2mgbM75tm3bOnfurNVqBw8efPLkyVdeeUWr1V6+fJlz/vrrr3t4eLRt27bhbQHYKNxZBqBFOOeCIODGKwC2C9cIAZqPc37kyBEiasywFwCwTv8PM4jzOxAjtXMAAAAASUVORK5CYII=" }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "γ_vals = range(0.05, 0.95, length = 25)\n", "\n", "models = [McCallModel(γ = γval) for γval in γ_vals]\n", "sols = solve_mccall_model.(models)\n", "w_bar_vals = [sol.w_bar for sol in sols]\n", "\n", "plot(γ_vals, w_bar_vals, lw = 2, α = 0.7, xlabel = \"job offer rate\", ylabel = \"reservation wage\", label = \"w_bar as a function of gamma\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "As expected, the reservation wage increases in $ \\gamma $\n", "\n", "This is because higher $ \\gamma $ translates to a more favorable job\n", "search environment\n", "\n", "Hence workers are less willing to accept lower offers" ] } ], "metadata": { "kernelspec": { "display_name": "Julia 1.0", "language": "julia", "name": "julia-1.0" } }, "nbformat": 4, "nbformat_minor": 2 }