{
"metadata": {
"language": "Julia",
"name": "",
"signature": "sha256:e7e9080b45a1c05cfe83e17f04ffa6650d374b0e2b727aa16190b4ce207fb617"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"using Stochy, Stochy.GadflySupport, Gadfly, DataStructures"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"@pp function rainprior()\n",
" flip(0.3)\n",
"end"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 2,
"text": [
"rainprior (generic function with 1 method)"
]
}
],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"@pp function sprinklerprior()\n",
" flip(0.2)\n",
"end"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 3,
"text": [
"sprinklerprior (generic function with 1 method)"
]
}
],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"@pp function wetgrassprior(rain, sprinkler)\n",
" flip(\n",
" !rain && !sprinkler ? 0.05 :\n",
" !rain && sprinkler ? 0.9 :\n",
" rain && !sprinkler ? 0.9 :\n",
" rain && sprinkler ? 0.95 : :unreachable)\n",
"end"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 4,
"text": [
"wetgrassprior (generic function with 1 method)"
]
}
],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"@pp function model()\n",
" local\n",
" r = rainprior(),\n",
" s = sprinklerprior(),\n",
" w = wetgrassprior(r,s)\n",
" tuple(r,s,w)\n",
"end "
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 5,
"text": [
"model (generic function with 1 method)"
]
}
],
"prompt_number": 5
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"@pp model() # Generate a sample from the model."
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 6,
"text": [
"(false,false,false)"
]
}
],
"prompt_number": 6
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Joint Distribution"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"joint = @pp enum() do\n",
" model()\n",
"end"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 7,
"text": [
"Discrete({(false,true,false),(true,false,false),(true,true,false),(false,true,true),(false,false,false),(true,false,true),(false,false,true),(true,true,true)}, [0.014,0.024,0.003,0.126,0.532,0.216,0.028,0.057])"
]
}
],
"prompt_number": 7
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"bar(joint)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"\n",
"\n"
],
"metadata": {},
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAAhcAAAF6CAYAAACqW3pRAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOzdfXxcVZ348e/33DuTTB6aNmna0qbJJH1ACOnMpIUiVGixLFaBXcXAwoo8rFS0hfW5uvtTqqzrsqvgLtQVcFlYVhTqw64iLOLaolBbdjL33IlTK2lyZ5JS2rTpU5pMMjP3fH9/JAOiVuo6x9OeOe/Xi9eLtEPymZDv9HbuuecCGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZRYtu2DYaSycwM1R2ypFKD9b29vRWqO2RJJAbmqm6QZcsWsl139yzVHbJs3947LZUaqlHdIUsq5c0hIqa6QxadZ6/UtP0hME4sFCosF0J8SHWHLPl8YcPYmB1T3SELIn1DdYMstbX9c4nsL6nukKWy0r4ulxu7UnWHLPk8PNDd3V+rukMWnWev1MzBRRmyLDYCAK+o7pAFkb2MaI+q7pCHXlJdIEsoxHIAkFbdIQsiHmCMDqnukAWR9YVCtb7qDnn0nb1SQ9UBsny468OhfHX+8vsevm/zb/v9T9x8c+2oH1qDiLnAaOCZezbfk/1jNxqGYRiGjrR85+KDN31wfr564guA9OHf9vtr164NjInQTwGgCwAuzVdNfO+PW6hWItHbyHmmXXWHLD09/Ut27NjToLpDFsfxVqpukCUe31vFef95qjtk6ekZWBCP9zWr7pDFcTIXxuPxgOoOWXSevVLT8uCCCfwKAFt+ot8P5oPXAEJm0yObuu57+L51gDhz/fXrW/+YjSoh2u0A4mrVHbL4Pv5FZWV+keoOedhnVRfIk50JwG5TXSGL7/uXWpa1QnWHLIhiA8CMKtUd8ug8e6Wl5cHFfY9suoIEfeSEDyBoB0L3tY/xF2iLpX+MtlOBZbEBIva86g5ZGMOt+bzYq7pDFkTxmOoGWSyLHQOgJ1V3yIIIDgDtVN0hCxF8u6GBTajukEXn2Ss1W3WAEkSzEeHVP1wR6RUAmKOw6I9qyZKWfgDoV90hy5Il4adVN8gUjbber7pBllis9QgAPK66Q5ZIpHWH6gaZYrHWR1Q3yKTz7JVaWR5cEEISgOYVPxYA9ZYPr/vbEueZDQDUJQR9tLOz9TnO048CwFmBgPUn2ez+Edue+TMiejkWa/1TztMxAHgQAJ6ORsOfdhzvBkS8DYD+MRptfdxx0v+ICKuEwPd2drbs5DzzDABMj0ZblicSA3MZE98DAB6Nht/HuXcZAH4eER+NRFr+yXXTHyOCP0eET0Qi4R+7buZhIjpHiPyaWGzhsOtmXgSAfdFo+PKenv4lvs8eAsBno9GWT3GeeQ8AfYgI7o7Fwo9x7v09AK5mzPq074sCAH0YEWZFo+Flrrt7FpH9FAD2RKMtN7lu5q1EdBcAPBaNhu92HO9DiPgeIvxULNbyLOfpfwWASCBAl7e3t+7j3HsRAIej0fCayfUc9AgRbInFwh93Xe/PifBjAPhP0WjLo46T/jwiXAYAt0SjYYfz9PcB4IyXXmpZ3tLy8vSKisIzAPCLaDR8veN4KxHxi0TweCwW/kfOM7cB0A2MwaeXLAk/7brpB4igE8D602h0/sucp38GQBWFgvizYJBVCYH/AQA/iUbDH3GcTBcibSCi+2Kx1oc5T38OAN4uBLu1s7M5znn6uwAwPxAYvSCXq6pCxB8RwUuxWPg61x14C5G4B4C+FY22/r3jpD+ICDcD4MZotOVJzjP/AkDnFgr+u5YtWzDAefp5AMhFo+FLHKdvEaL1DQB6IRpt/SvO0+8CgL9GhK9GIuGvuW7mDiK6ApHWRSKtOzhPfwsAwoVC8KKammPW+HjlFgDqi0Zbr0kk0m9GpG8i4oPRaPhvOffeD4C3ENGdsVjrf3GevhcA3owIXZFI2HOczHMAALFYy8Wum24lgs0AsD0aDa93HO9PEfHTAPRgNNp6P+fp/wcAf0YEt8di4W2cZ74JQAsrK8dXHT8+zbft3E8AIB2Nht89uS6CfQURvx+JtHzWddPvI4JbAeDvotHwdxzH+zIiriDyr43FFvRynv4xAASj0fCKeLyv2bat7wDg/0ajLR/gPHM5AG0kgocA6DHG4CtEuBiRfTgSaf6p66a/TgRnEtHqYHBsLJ+v3gYAg9Fo+J3JZGapEHQ/ADwVjYY/4zjejYi4ngjvisVaNnOe/hIAXOz7cP3SpeFfcO79EABro9HwmzkfnAfg/xcROLFY+JZkMr1GCLgTAB+JRlvudZz0xxHhGgD4eDQa3sJ5+t8B4OyJCfuyTGbekcWLMzsA4JVoNHyF43hRRPwaETwTi4X/hvPM9QD0V4j0xUik9Zuum/4HIriECO4GoJ2I+AUAaohGW89Lpbw5+Tw+CQDJaDR8s+NkLkWkLxDRf8RirV/mPP0RALhOCPxkZ2fLjzjP/BsAdSAW3h6JLBziPB0ngqFYLPz2ZHKwQwj/3wDof6LR1g2Ok74OET4CgF+ORlv+g/PMFwDoUiHEX3Z2trmcp58EgDmRSMt5jrO7gbHA00SUisVab+A8vQoA/hERvhmJhL/oupm/IqLrAehvotHWZxwn/SAixIRgV3Z2Nu/lPLMDgKrq6ujcQ4ew1bLgUQDYGo2GP8a5dw0AfpyI7o3FWh/hPH0nAKxBhLWRSDjhON5/IeK8QuHgm0Oh2bX5vP9DANgVjYbfk0z2XyQEuxsAN0ejLXdxnl4PADcSwWdisfBTnGe+CkDLLKvwzo6OhYOum95GBGPRaHg1532LAazHEPGnkUjLh103824i+iQAfCUaDT/kON5GRLyciH0gFmv+X87T3wGA5ro6WrF/fyFYWRn4MQD0RqPhax0ncyGAeIwxvD8SCf8d596tAPg+RPxsJNLyfddNbyKC5UT07lisNc259xMA9KPR8KqenoEFvi8eB4CfRaPh2xzH+zNE/H9EeH8s1vIg5+lPA8CfCgG3dXaGf+Y46c2I0IpYeTEAANH4c0TgxWLhLtf1zifC+wDgv6LR8J2Ok7kFkd5PRH8bi7X+J+fpfwaAC3zf+vOlS+fv5jy9BYCsaLT1Isfxwoj4LSJ4MRYLf9B1M1cQ0R0A9LVotPWrf/ifrmVg3XvXvXn9Teu2FT/u6uqyPnDzB9oAANbfsP7C9Teu2wIA8IHrPjBj3Y3rvHXvXaftAsBf5zjeSs49bc8dcu7d5bre+ao7ZCkeMOgoHu9rnjqQ1xLn3q2Ok75OdYcsnHvfi8f76lR3yKLz7JVa2bxz0VDZ0GAJ1gcAuH9s//bZ1Y37b7tx3YsE0IwEf3fvv28aVt34x0JUSCEGD6jukMWy6OtjY8GXVXfII+5QXSBP6CDA+L2qK2SxLOvZiYl8XnWHLETsLoDhMdUd8ug8e0bJrL95/dxP3HyztrvJGYZhGIZh/FEkEv0Rx/FuUN0hi+NkbunuTp+lukMWztN3q26QxXV3z+Lc+6TqDlk4z1zuuulLVHfIwnn6TtfdV626QxadZ6/UtLwU1fjdEHEGIrSp7pAFUSy0bdL2vC8RanvZdD6PlQCo7QZvAKJJCJ2vTKNIPj+q7el2nWev1LT9ITBOLJu1d9TUCPeNH3l6CgTsu4LBcW3vLUKE16pukGVkpG1vfX3fR1V3yDI+XnistrZOqO6QJRCAtZFI24jqDll0nj3DMAzDMAzDOLWYS1FPbzpfDmcuRT29mUtRjSKz5qIMIbICIoyr7pCFCLK+bxVUd8iCSNq+7RwIVAgA1PaUFgBMMIY51RGyEOHxQkGQ6g5ZdJ49wzAMwzAMwzi17Nq1qzaRGJirukMWzgfn6Xw5nOv2n6m6QZZUKhV03bS2dyhOJHobU6nBetUdsnR3Dy584gmyVHfIovPslZo5LVKGstmKpYz571fdIU/hdoBsh+oKWYSwHlDdIEs2WzmHCD6nukMWxuyrcjn/bao7ZLGswt1tbf01qjtk0Xn2Ss1cilqGENkBAEip7pCH9QhR0HY7d0SxVXWDPKExgPyLqitkIWJ9iJBV3SEPewGgQtvtzfWePcMwDMMwDMM4lSQSA3MTiYFlqjtkcV1vuevunqW6QxbXzVyhukGWVGqoZupW31rq7k6f1d09uFB1hyyce5elUqmg6g5ZdJ69UjNrLsoQor+YMV/bISGCdwFY2m5vLgR8THWDLNnsSD0A3Ky6QxbLoosZ889T3SHRumy2MqQ6QhadZ6/UzJqLMuT7ot+2bW33gSCCp4WwBlR3yEP3qy6QpVAQR2zbflx1hyxEsJ1I330uEPGRujqh7R46Os+eYRiGYRiGYZxaOM+0c565WnWHLJxn3uM4fYtUd8ii89bt8fjemZxnblPdIYvjZC7lPLNCdYcsnGc2xON7q1R3yKLz7JWaWXNRhohEI4DQ+bbWHYyxBtUVshCxlaob5MlWAZC2axIQxQIialbdIY+4ECAbUF0hi96zV1pmzUUZsizmIuKg6g5ZiMTXCgUYUt0hCxFqe0vyhgY2dOwYu0t1hyyIgScDgbzO650+MzLSpu29YXSePcMwDMMwDMM4tbiut9xxvA+p7pCF88wGx/Giqjtk4Tz9DdUNsiQSA3M5T9+tukMWx0lf57oDV6rukMVx0g/u2rWrVnWHLDrPXqmZ0yJlSAgIIcIM1R3yiHrGoFJ1hSxEqO1N54TI24xZjao7ZEGkaUII1RnSINLs48cD2q7l03n2DMMwDMMwDOPUsmUL2du2DWq7i962bYOhLVtI23fldH7bmYiY6+6rVt0hS29vb4XO22OnUkM1RISqO2TRefZKTdu3r4wTmz49vaKqqvBJ1R2yVFUVNtbXp7W9d0o2G3pSdYMs3d39TUTjX1XdIcvoqH1TLlf9btUdsuTzo491d/dPU90hi86zV2rm4KIMWRYbAYBXVHfIgsheRrS1vRwOgF5SXSBLKMRyAJBW3SELIh5gjA6p7pAFkfWFQrW+6g559J09wzAMwzAMwzi1JBK9jZxntN2hs6enf8mOHXu03aHTcbyVqhtkicf3VnHer+0OnT09Awvi8T5td+h0nMyF8Xhc2x06dZ69UjOnRcoQot0OILS9t4jv419UVua1vbcIANP4/gbZmQBM23uL+L5/qWVZ2t5bBFFsAJih7b1F9J690tJ2Rb1xYohiL2KgW3WHPLgDAParrpCFMfq+6gZZQqHg8Xze36K6QxbGaKcQpO0tyYnwmVwuqO0t5XWePcMwDMMwDMM4tSSTmTbHyVyqukOWZDK9Rufz2px771fdIIvjeNM5965R3SHL1Nb72m5N7zjeDZ7nabs7rs6zV2pmzUUZ8n3RjCi0Pe8rBK0MBJi22/QSsetUN8ji+2IaAF6uukMWIogB4NmqO2RBhKuGh0WF6g5ZdJ69UjNrLsoQUSGFGDygukMWy6Kvj40FX1bdIY+4Q3WBPKGDAOP3qq6QxbKsZycm8nnVHbIQsbsAhsdUd8ij8+wZhmEYhmEYxqkkkeiPOI53g+oOWRwnc0t3d/os1R2y6HxLctfdPYtzT9ut6TnPXO666UtUd8jCefpOne8No/PslZpZc1GGEHEGIrSp7pAFUSy0bapT3SELES5V3SBLPo+VAKjtBm8AokkImKO6Qh6K5POj2p5u13n2Sk3bHwLjxLJZe0dNjXBVd8gSCNh3BYPj2t5bhAivVd0gy8hI2976+r6Pqu6QZXy88FhtbZ1Q3SFLIABrI5G2EdUdsug8e4ZhGIZhGIZxanEcbyXnnrbb2HLu3eW63vmqO2RxnMxzqhtkicf7mjlPP6q6QxbOvVsdJ63t5Yyce9+Lx/u0PSWp8+yVmllzUYYQWQERNN6CGLK+bxVUd8iCSNq+7RwIVAgA1PaUFgBMMIbabo9NhMcLBUGqO2TRefYMwzAMwzAM49Sya9eu2kRiQNsdLDkfnKfz5XCu23+m6gZZUqlU0HXTrao7ZEkkehtTqcF61R2ydHcPLnziCbJUd8ii8+yVmjktUoay2YqljPka75FfuB0g26G6QhYhrAdUN8iSzVbOIYLPqe6QhTH7qlzOf5vqDlksq3B3W1t/jeoOWXSevVIzl6KWISI6DID9qjtkIWK7fZ+Oqu6QBZG6VTfIEgjQOBGlVHfIw/YwRhpvj41uIFCt83onbWfPMAzDMAzDME4ticTA3ERiYJnqDllc11vuurtnqe6QxXUzV6hukCWVGqrhPL1KdYcs3d3ps7q7Bxeq7pCFc++yVCoVVN0hi86zV2pmzUUZQvQXM+ZrOyRE8C4AS9vtzYWAj6lukCWbHakHgJtVd8hiWXQxY/55qjskWpfNVoZUR8ii8+yVmllzUYZ8X/Tbtq3teVEieFoIa0B1hzx0v+oCWQoFccS27cdVd8hCBNuJ9N3nAhEfqasT2u6ho/PsGYZhGIZhGMaphfO+xa47cKXqDlk4T79L570SXDet7VuzqdRgPedpbU+LJJP9F3Her+1pEcdJr9u2bVDb0yI6z16pmTUXZYiIzSXyNb51MC0HoNmqK2QRArVdL5PN5moAQNsFnULg2URM2wWdiHRZMJjTdkGnzrNXambNRRmyLOYi4qDqDlmIxNcKBRhS3SELEWp7S/KGBjZ07Bi7S3WHLIiBJwOBvM7rnT4zMtKm7b1hdJ49wzAMwzAMwzi1uK633HG8D6nukIXzzAbH8aKqO2ThPP0N1Q2yJBIDczlP3626QxbHSV+n83onx0k/uGvXrlrVHbLoPHulZk6LlCEhIIQIM1R3yCPqGYNK1RWyEKG2N50TIm8zZjWq7pAFkaYJIVRnSINIs48fD2i7lk/n2TOMP9iWLWTrvKJ727bB0JYtpO2Bs85/MyQipvMdbXt7eyt03sEylRqqISJU3SGLzrNXatoeYRqGYRiGoYY5uChD06enV1RVFT6pukOWqqrCxvr6tLb3TslmQ0+qbpClu7u/iWj8q6o7ZBkdtW/K5arfrbpDlnx+9LHu7v5pqjtk0Xn2Ss0cXJQhy2IjAPCK6g5ZENnLiLa2l8MB0EuqC2QJhVgOANKqO2RBxAOM0SHVHbIgsr5QqNZX3SGPvrNnGIZhGIZhGKeWZDIzo6dnYIHqDlkcp29RPN5Xp7pDlkRiQNtTPp7nVbruwDmqO2Rx3T1NqZQ3R3WHLI7jRXVeTK3z7JWaOS1ShnxfRHzff6/qDlkQ2fsCAXaW6g5ZEOlLqhtkGR4Ws4jEBtUdshDlL8/l8BLVHbIgwudqa/u1vdpH59krNW2PMI0TQxR7EQPdqjvkwR0AsF91hSyM0fdVN8gSCgWP5/P+FtUdsjBGO4UgbW9JToTP5HJBbW8pr/PsGYZhGIZhGMapJZnMtDlO5lLVHbIkk+k18Xhfs+oOWTj33q+6QRbH8aZz7l2jukOWqa33td2a3nG8GzzP03Z3XJ1nr9TMmosy5PuiGVGsUN0hixC0MhBg2m7TS8SuU90gi++LaQB4ueoOWYggBoBnq+6QBRGuGh4WFao7ZNF59krNrLkoQ0TWS0RwXHWHLIjwHQDfU90hC2PwRdUNsoRCtYfy+bGHVHfI4vv4HICVV90h0aZQaDyrOkIWnWfPMAzDMAzDME4tiUR/xHG8G1R3yOI4mVu6u9PaXoqq8y3JXXf3LM49bbem5zxzueumtb0UlfP0nTrfeE7n2Ss1s+aiDCHiDERoU90hC6JYaNuk7SZaRLhUdYMs+TxWAmC76g55RJMQoO0mWgAUyedHtT3drvPslZq2PwTGiWWz9o6aGuGq7pAlELDvCgbHtb23CBFeq7pBlpGRtr319X0fVd0hy/h44bHa2jqhukOWQADWRiJtI6o7ZNF59gzDMAzDMAzj1MJ5ZoXrep9S3SGL43gbdb4HAOdpbW/77Lp7mjjPaHvLdc69m1w3o+0t1x0n/dj27b3a3nJd59krNbPmwjAMwzAMwzAMwzAMwziF7Nq1qzaRGNB2B0vOB+fpfDmc6/afqbpBllQqFXTddKvqDlkSid7GVGqwXnWHLN3dgwufeIIs1R2y6Dx7pWZOi5ShbLZiKWO+xnvkF24HyHaorpBFCOsB1Q2yZLOVc4jgc6o7ZGHMviqX89+mukMWyyrc3dbWX6O6QxadZ6/UzKWoZYiIDgNgv+oOWYjYbt+no6o7ZEGkbtUNsgQCNE5EKdUd8rA9jNGY6gp50A0EqguqK2TRefYMwzAMwzAM49SSSAzM1flSTdf1lrvu7lmqO2Rx3cwVqhtkSaWGajhPr1LdIUt3d/qs7u7Bhao7ZOHcuyyVSgVVd8ii8+yVmllzUYYQ/cWM+doOCRG8C8DSdntzIeBjqhtkyWZH6gHgZtUdslgWXcyYf57qDonWZbOVIdURsug8e6Vm1lyUId8X/bZta3telAieFsIaUN0hD92vukCWQkEcsW37cdUdshDBdiLMqe6QBREfqasT46o75NF39gzDMAzDMAzj1MJ532LXHbhSdYcsnKffpfNeCa6b1vat2VRqsJ7ztLanRZLJ/os479f2tIjjpNdt2zao7WkRnWev1MyaizJExOYS+RrfOpiWA9Bs1RWyCIHarpfJZnM1AKDtgk4h8Gwipu2CTkS6LBjMabugU+fZKzWz5qIMWRZzEXFQdYcsROJrhQIMqe6QhQi1vSV5QwMbOnaM3aW6QxbEwJOBQF7n9U6fGRlpG1XdIYvOs1dqKOsTX7PmmgXA/CtI4EWEMBcBJhAhKQB6LBLuN5/69osAQLK+vmEYhmEYapT8tEjXFV2tV6/p+g6h6CXCTwPAHATcCYQuAbQhwUcFsJ90vaOrp+sdXdqe9z+VJZOZpZx7t6rukMV1M7cnk4Pabv/tOOkHVTfIkkp5czhP36m6QxbXTV/lup6223+7buaeVGpI2+2/dZ69UivpaZGr11x9C/j0MWD4IBF9cPNTm/f9tsetXLnSnl0980oC/HjX27tu2fzUZnMe64/I90UtIpyhukMWIjEPQGh74zIAXKy6QJZsVgRt2wqr7pCFiBqJsEJ1hyxEYkE2O6Ltjct0nr1TWtfbulasWbPmDQdn48aNr75j0rWm62K5Vcav27KFbJ1XdG/bNhjasoW0XU+0a9euWtUNshAR0/mOtr29vRU672CZSg3VEJG00+2q6Tx7pVbSF+DN/735+eK/X7Om67MA7OnHn358+68+puuKrtbUi6lHAOAiAIDNT29+rpQN0myhlQAwXXXGHyAOq3CP6gjDMAxDf/L+docQJBDPX7Om6/P7swfu3Lp1a+HqNVffAoLuBoDT70oFhC8BQafqjP8zgr8AgMcAAKZPT69AhFUAcIfaKDmqqgobq6vT3wWA7W/44NNQNht6EgC0fMevu7u/ybatzwPA9apbZBgdtW8iChyDqVnUTT4/+lh398j1AKDlXYl1nr1Sk3Zw8fhTmz91zTve7RLgvbOqGtdc/Y6uA0B0GRD+wwiNfFbW1wUA+MTNN9eO+qE1iJgLjAaeuWfzPdlff8wHb/zg2xiwV9/iygVz//nAAw/kZXadKhiDrBBwWHWHPOyQEPpuQYxIe1U3yMJYoAAgDqjukIUIjzHGjqvukIUI99fU5IXqDll0nr1Sk35u7Oo1V98MSP86+dXozid+8K3PyPx6a9euDQRzgf8lgl5EHAISi+975CuX/vrj1t24zkOA7xc/Pg6jn3j44YdP/AfSVuo+7d+5uAS1/NuSYRiGcWqR9s7FtVdcO1P4hX8ioOsQ4F8JaT8Q/vXVb+86mwDWn+hKkj9UMB+8BpAymx7e1AUAsP7Gdc7669e33vfofV7xMeveu64BELz7/m3T7TIaTnXJZGYGItZ3dDT3qW6RwXH6Fvk+DC1btkDLt2YTiYFlnZ3NcdUdMnieV3nsmLUwEmn+ueoWGVx3T5Nt5wvt7a1SXv9UcxwveuRI+OerVqGWG4XpPHulJm37b9/PbyKASxDoisef2vy+J37wrb9hIC4EgA4E+IWsrwsE7UDovvYx/gJt8fqtri1YiAQz1t+47ofrblz3H+tuWHeRtJ5TkO+LiO/771XdIQsie18gwM5S3SELIn1JdYMsw8NiFpHYoLpDFqL85bkcXqK6QxZE+Fxtbb+2V/voPHulJnFBJ3sqH8x/8Lvf/e5w8Ze++dS3d1xxxRWxKlEpb3tfotmI8OpVK4j0CgDMeV2aQBuRXiCGfwc+nYUIm9euXdv+wAMPHCw+Jpn03lQosCaAnNvZuehANYjq0dP8VizJ5GBHoSBmC5E/bNvBbtcdeIvvU0VnZ8uPent7K0ZGgm9BxKOxWPP/uu7uWb4fWMJYYSAaXfCS4/QtIrJbAgG/p6OjbX8iMbAMgKbPmCGeb21tHXfdzFsBMBeJNP90+/beacFg8LxAwB/q6GhL9vTsnp/PB84E8Hs7O9syrjtwju/THEQRj8Vaj3CeWSEEVMZizf+zc+fOwMREzUWM+cei0bYXE4neRoBgxLL8wUik7Zfd3YMLEUW4okL8vL29dV8ymVlaKMCM8XH2wgUXzM+6bvoSRBghgv27du2qHRsLLQfwD3R2trmuu6fJ9/03IYrdsVhrmvNMuxBwhm1D95IlLYcdJ30BEVbFYs0/7u7uthhrvNiyxEgk0rojHt87k7F81LbFniVLWnf19AwsyOepVQhKLVsWfsV1052+j/VCBLYtWzZ3jPP0KiFIdHa2PpdKDdVMTGTPRxQHY7FWzvngPCHEWZZFfZFI2EskMmcDwNyKCpZob59/KJFIvxkAq3fvbt7S2AhYVzew0rLE8UikdfuOHXsaEP1d3d3ps5YuDf8imcy0FQrQBoA7Ozub93KejgmBDZZV8bNIZM6o43grAQBisdatrruv2vcn3swYDUejYSeRGJgLQGfbNvQvWdLS392dPgsR5+XzlrN8edOw63rn+z6rOXq0eeuBA0ALFw6sAqDRzs7wz1KpwfqJCdEJAHs7O1t2um661Wzb3oUAACAASURBVPdxAWPsF9Ho/Jcdx4sSsZkVFaHt7e2zjicS3sWMIYtGw1vi8b1VjOUvsCw6FImEE/F4+gzGsD0QQI8xPJzL+bsSicxqIQJ82bK5B13XW+77rFaIA88tXbrUd5yBSxBpLBYLb0smMzMKBVjKGLwSjbakHMcLE7GFlmXtikSa9iQS/REAqzGXy714/vmLjiWT/RcRMTsSCf9427bBUGWluBAAj3R2NsdTKW/OxAQ7h4illy6dv9t1+8/0fWu+ZeWTkcjCIc77zxPCmlZRcfwnZ599dt5xBt7KGIxHoy3PO443nYgtsyzcF4k0/zyR6G8BsBYFAvlfdnQsHOzp6V+Sz1uzGKOMEGLYdQfeAkDBSKTlfzzPqzx8mK0ozl5PT//sfN7qQCxkYrEFvZz3LRbCbi7OnuMMnEtEdbW1uZ8uWrRoIpHIrLYsnIhEmn8aj/fVMWafW5y9eLyvmTF7cXH2iq8Bv232du/eHRwZCb6lOHvF14Di7BVfA4qzV3wNeG32Mm8FEH25XDBXfA0ozl7xNaA4e8XXgNdmL3MhEYRiseYf79y5056YqLmoOHvF14Di7BVfA4qzV3wNKM6e66Yv8X3yOztbnyu+BhRnr/gaUJy94mtAcfaKrwG7dzdvaWvrZlOvAccjkdbt8fjemYj5V2ev+BpQnL3ia0Bx9jhPryIiisVatxZfA4qzV3wNKM5e8TWgOHvF14CjR5u3AgDU1Q2sLM7ejh17GgIBP1acveJrQHH2iq8BxdlzHG8lImI0Gt5SfA0ozl7xNSAQQK/U72SX9E/Lriu6mov//sQPnnjkVw8sir7//e+PPf6Dzbf9yn9T0rtXEkKSgOYVPxYA9ehD6lcfc98j971w78Ob1t/30H17Nz2y6X8Q4JmKXOB1u4X6PpzJmFgNEGwEAKhEqCplpwq+7y9hTKwmwgORSPP3APwVk88R4OhRVjn57+JcAADEitmTH9tTm8ZYi6ceOxsAwLL8cxkTq48cgUoAACJYDSDeAgBQXY11jInVvm8tAQAoFALNjInVtg0tUyUdU5+reGnvhVMfo2VZFZP/zs4DAGAsOGuq40wAgEDAX8SYWJ3LTR4wCiGWMSZW19VBaPJjeKsQAJFI2Bsfr57GmFiNiNHJ51+YP/V1wpNfVpwz9fv1v9qxdSuwQKApONVxPgBAMDgxc/I5wZumvpcLGROrA4FXNyPrnPrcxZ+TSxhjU6vKR2smvw7EJr+3haapx7ZOfi/h7MnnlGuYfM5wAWNidXv7TquiYk9g8v8ZvBkAoKJiogFAHLYsOGuyAxZMfe7iz3xs8nMfr578WrBq6sogADhePfXYzqnfmzf1nBZMdZzFmFhdWZmbOfW9PJ8xsbqiYk+gvX2nxZhYzRhcAAAwPg71jInVlgVnT33d1snPXWia+jjKmFidzY7VTj4ndjEAFP/GXsWYWC0ELp38fwpnTHb4C9vb5x8CoNTk58o1Tj6cLZ/8XjcFt24FNvW9u3Dqc82Y+rk9Z+rj8OTnKsyffE5WZPI5VdZOPid2kRDwVgCAujoITT6nyZv45XIwZ/Lr+IsmP5V95uT3PlDsOG/yOVsVAIC/1jF98mO/AwDAtqGFMbG6UAg0T/5/spZMfX9+Ho22vQgg3jI5MwBHjkDlVMeyqc81NXvW1OzZi6c6Zk1+LM5lTKw+epRVTn5vxWoAf8XUf1s3+XXZVIfdPNlsT3VMzh4R1U19rldfA4aHoWLy99i5k79XnL1ih7WIMbF6YmLypoCM+csYE6trakTl1M/LW4XAgxdcMD8bCoWmTX2vIgAAhYI9f+rnuPgacM5Ux+teAzZvBgRoDE51nAcAQBQozt6Zk8+psHDq9WTO5HOanD3fF8W9ey5hDC8CABgbq6qden7Ryd/KN011hKce286YWD0+DvWTz2ly9trauhnA3MDkzymcP/lzOl6cvanXgMKCqc89d+pzxSY7RoqvAasAYOXk92Zy9gAmXwMA8vOmOqZeA+isyVmbKL4GvJkxsbqpabfV1LS7OHtvBgAIBiemZo/Omvz+TM4eQH7qNYCikx+PFndKXTnVAr4/UjXVUVw/OHdqXhZAiZV0Qec1b+/6AgFcQAj/0H5u+9MbN2484arha9921RKfWR8ioLM3P7X5/FI1rL9h/YWA9Lf3Pbxp1Qeu+8AMFmQJELBsKDt0ZGbtzJZ/eehf+tffsP4vEemcex/e9OGpBaDdSOyd9z5y74mP3DRa0BmP9zXbtt0cjbY8/0b/2enIcbyVRNZLnZ3NWq7sdpz0dbFYWMvFudu3906rrAxeFI22PKm6RYapd3RynZ0tO1W3yOC66auqqvJPLlq0aEJ1iww6z16plfS0yONPbf7UNWve/ScAeM/OF1P/evXbu3YC0m4QsBsZZAFwMRG8CQDO9AFqAeBzUA23lLJh/9j+7bOrG/ffduO6FwmgGQn+7t5/3zR86/W3zrIE6wMArBSV3xy3s8+uv+mDT0GORQHgu7/zwEIzlsXaAMQqANDy4AIR1jDmjwOAlgcXAPh+0HSfBNtm0wHoGgDQ8uACEc4HoGMAoOXBBRHdcPQo+xEAaHlwofPslVrJ11w8/vS3frhx48aOndt3vg0YnEcEHQxhJRHMA4D9APQCAjwmxnDz5q1PlPx6782bN/sA8Ofrb14/twrGRv7hoYdGAAC++uhXh2DqnZovPvrFUQC44EPvWXsG2FXZLz/85SOl7jiVEVkvEYG219ojwncAfO+NH3l6Ygy+qLpBllCo9lA+P/aQ6g5ZfB+fA7B03k9nUyg0/hv7CulC59krNW33gC85jU6LGIZhGIZMUm/u1NXVFYRRuIYBnE1IAgX8UiD+UNYeF8bJSST6I4gYjcVaH1HdIoPjZG4Rgp5fujQs75JnhThP3x2Nhj+iukMG1909i8i6ORpt/XvVLTJwnrkckcYikfCPVbfIwHn6TsTKv49E5oyqbpFB59krNWnXVl575bWzcQxSCPA1AngHEF5BiA8jQKrrHV3vkPV1jTeGiDMQoU11hyyIYqFtF1fD64cIl77xo05P+TxWAmC76g55RJMQr780Xi8UyedHtb0jsc6zV2rSfgj8gv8FADha8P3wd575zisAAO/8k3fOCgSsv2WE9wDAD2R9beN3C4UmusfGql5S3SGP/c8AAW3X0TDmr1XdIEsoNL6vUKiWeosAlYQofLuiotJX3SGL79sf6e9v0nY9l86zV2rS1lxc/fauOCDe+8QPnnjdW+/vfOc7GwIT9qDli3O/8cy3Uyf67085Zs2FYRiGYZwUiVtO4g4Qr27s8yp7xC4AAI5W5AfkfW3jd+E8s8J1vU+p7pDFcbyNk7sH6onztJaXaQJM3nuD88xXVXfIwrl3k+tm3q26QxbHST+2fXvvNNUdsug8e6Um7bSIJfz7fcZ+0PX2LsGIbUXCfT7zFyHAHQjw9YAfaLzqHVc1AgB0nNuR/l0bbhmGYRiGcfqQdlrkmrd3/TsBXH8yjw1QsO7rT3/9mKyWkjCnRQzDMAzjpEg7LZL3/Q0k4E0n88+i5Yu0XQB0Ktq2bTCUTGZmqO6QJZUarO/t7a1Q3SHL5M2G9LRlC9muu3uW6g5Ztm/vnZZKDdW88SNPT6mUN4eITu87PP4OOs9eqUk7LTJ1hcgrJ/PYzf+9WVaG8VuEQoXlQsAqALhDdYsM+XxhQ6FgfxcAtqtukQGRvgEAF6vukKG2tn8ukf15OMl3PU83lZX2dbnc2DHQdAvpfB4e6O7uvx4AjqpukUHn2Ss1ba9HNk6MiA4DYL/qDlmI2G7fJy1f3AAAEKlbdYMsgQCNE9HpcxXZ743tYYzGVFfIg24gUF1QXSGLzrNXamb775Nl1lwYhmEYxknR9tyYcWKJxMBcnS/VdF1vuc7n7V03c4XqBllSqaEaztOrVHfI0t2dPqu7e3Ch6g5ZOPcuS6VSQdUdsug8e6VmDi7KEKK/mDFf2yEhgncBWNpuby4EfEx1gyzZ7Eg9ANysukMWy6KLGfPPU90h0bpstjKkOkIWnWev1MyaizJkWWzA9+F51R2yMIZbczl/r+oOWRCFtqe3LIsdAyBtNypCBIeIJlR3yEIE325oYNo+P51nr9TMmouTZdZcGIZhGMZJMadFyhDnfYtdd+BK1R2ycJ5+l+umW1V3yOK6aW3fmk2lBus5T2t7WiSZ7L+I835tT4s4Tnrdtm2D2p4W0Xn2Ss0cXJQhIjaXyNf41sG0HIBmq66QRQjUdr1MNpurAQBtF3QKgWcTMW0XdCLSZcFgTtsFnTrPXqmZNRdlyLKYi4iDqjtkIRJfKxRgSHWHLET4UdUNsjQ0sKFjx9hdqjtkQQw8GQjktd0Hggg+MzLSNqq6QxadZ6/UzJqLk2XWXBiGYRjGSTGnRcpQMplZyrl3q+oOWVw3c3syOdihukMWx0k/qLpBllTKm8N5+k7VHbK4bvoq1/XeprpDFtfN3KPzvVN0nr1SM6dFypDvi1pEOEN1hyxEYh6AqFbdIQ8uVl0gSzYrgrZthVV3yEJEjUSo7U31iMSCbHbEUt0hj76zV2rm4KIMHTkSfr6iYs8O1R2yjI3ZGycmmvKqO2QJhbKXq26QZenStj3J5H5t31Wrri78Wy6XI9UdsgQC1ddFIo3arrnQefZKzZwWMQzDMAyjpMzBRRmaPj29oqqq8EnVHbJUVRU21tentb13SjYb0nYHy+7u/iai8a+q7pBldNS+KZerfrfqDlny+dHHurv7p6nukEXn2Ss1c3BRhhiDLBEcVt0hDzskBIyrrpAFkbTd2pyxQAEADqjukIUIjzHGjqvukIUI99fU5IXqDll0nr1SM5einixzKaphGIZhnBTzzkUZSiYzM3p6Bhao7pDFcfoWxeN9dao7ZEkkBrQ95eN5XqXrDpyjukMW193TlEp5c1R3yOI4XnTLFtL2QgGdZ6/UzMFFGfJ9EfF9/72qO2RBZO8LBNhZqjtkQaQvqW6QZXhYzCISG1R3yEKUvzyXw0tUd8iCCJ+rre3X9jJwnWev1LQ9wjRODFHsRQx0q+6QB3cAwH7VFbIwRt9X3SBLKBQ8ns/7W1R3yMIY7RSCtF0PRITP5HLBnOoOWXSevVIzay5OlllzYRiGYRgnxZwWKUPxeF8z55kVqjtkcRxvZSIxMFd1hyyOk75OdYMs27f3TuM8o+1GRY7jRROJzNmqO2Rx3fRVvb292u5AqvPslZo5uChDlsXaAMSlqjtkQYQ1luU3q+6QB9+vukAW22bTAega1R2yIML5iBRV3SELEd1w9CirVN0hj76zV2pmzUUZIrJeIgJtr7VHhO8A+J7qDlkYgy+qbpAlFKo9lM+PPaS6Qxbfx+cALG23pgeATaHQeFZ1hCw6z16pmTUXJ8usuTAMwzCMk2JOi5QhzjPtnGeuVt0hC+eZ9zhO3yLVHbJw7n1WdYMs8fjemZxnblPdIYvjZC7Veb0T55kN8fjeKtUdsug8e6VmDi7KEJFoBBDtqjvkER2MsQbVFbIQsZWqG+TJVgHQeaorZEEUC4hI4/VA4kKAbEB1hSx6z15pmTUXZSgUmugeG6t6SXWHPPY/AwSOqK6QhTF/reoGWUKh8X2FQvVnVHfIIkTh2xUVlb7qDll83/5If3+Ttuu5dJ69UjNrLk6WWXNhGIZhGCfFnBYpQ5xnVriu9ynVHbI4jrdR53sAcJ7W9rbPrrunifOMtrdc59y7yXUz2t5y3XHSj23f3qvtLdd1nr1SMwcXZYhI2ESg7bXoiBCyLF/bU35EWKu6QZZ8foIBkLb3pgCACiEoqDpCFkSqsW2m7TviOs9eqWn7Q1By5rSIYRiGYZwU885FGdq2bTCUTGZmqO6QJZUarNd5C2KdtzbfsoVs1909S3WHLNu3905LpYZqVHfIkkp5c4hI2z9XdJ69UtP2h8A4sVCosFwI8SHVHbLk84UNY2N2THWHLIj0DdUNstTW9s8lsrW9rXVlpX1dLjd2peoOWfJ5eKC7u1/bUwc6z16paXte2jgxIjoMgP2qO2QhYrt9n46q7pAFkbpVN8gSCNA4EaVUd8jD9jBGY6or5EE3EKguqK6QRefZKzWz5uJkmTUXhmEYhnFSzGmRMpRI9DZyntF2h86env4lO3bs0XaHTsfxVqpukCUe31vFeb+2O3T29AwsiMf7tN2h03EyF8bjcW136NR59krNHFyUIUS7HUBoe28R38e/qKzMa3tvEQCm8f0NsjMBmLb3FvF9/1LLsrS9twii2AAwQ9t7i+g9e6Vl1lyUIctiA74Pz6vukIUx3JrL+XtVd8iCKLQ9vWVZ7BgAabtRESI4RDShukMWIvh2QwPT9vnpPHulZtZcnCyz5sIwDMMwToo5LVKGOO9b7LoD2l4Ox3n6Xa6bblXdIYvrpj+mukGWVGqwnvP0zao7ZEkm+y/SeU2J46TXbds2GFLdIYvOs1dq5uCiDBGxuUT+UtUd8tByAJqtukIWIfAK1Q2yZLO5GgBYpbpDFiHwbCK2UHWHLIh0WTCY03Z7c51nr9TMmosyRFRIIQYPqO6QxbLo62NjwZdVd8gj7lBdIE/oIMD4vaorZLEs69mJiXxedYcsROwugGGN9/HQefZKy6y5OFlmzYVhGKX2Y2oHhL9VnfF/hnAIVuJfqs4wTj3mnYsylExmlgohzo1GW7W8tbXrZm5HZFuWLJnfo7pFBsdJPxiLhW9R3SFDKuXNyedxXTQa/rTqFhlcN30VAI1GIq3/DQAACI0A8Gdqq/4ABPt+9UPXzdxj26FPt7fPOq4qSSadZ6/UzMFFGfJ9UYsIZ6jukIVIzAMQGt+2GxerLpAlmxVB27bCqjtkIaJGItT2pnpEYkE2O2Kp7pBH39krNXNwUYaOHAk/X1GxZ4fqDlnGxuyNExNN2p7XDoWyl6tukGXp0rY9yeT+W1V3yFJdXfi3XC5HqjtkCQSqr4tEGkdVd8ii8+yVmllzcbLMmgvDMEptC60EgC2qM/4A+2AVavsuqPF/Zy5FLUOO463k3NN2G1vOvbtc1ztfdYcsjpN5TnWDLPF4XzPn6UdVd8jCuXer46SvU90hC+fe9+LxvjrVHbLoPHulZg4uyhBjkCWCw6o75GGHhIBx1RWyIJK2W5szFigAgLaXSRPhMcaYlosdAQCIcH9NTV6o7pBF59krNXNa5GSZ0yKGYZSaOS1iaMq8c1GGksnMjJ6egQWqO2RxnL5FOr81m0gMLFPdIIvneZWuO3CO6g5ZXHdPUyrlzVHdIYvjeNEtW0jbCwV0nr1SMwcXZcj3RcT3/feq7pAFkb0vEGBnqe6QBZG+pLpBluFhMYtIbFDdIQtR/vJcDi9R3SELInyutrZf28vAdZ69UtP2CNM4MUR2AABSqjvkYT1CFIZVV8iCKLaqbpAnNAaQf1F1hSxErA8Rsqo75GEvAFRoexm43rNXWmbNxckyay4Mwyg1s+bC0JQ5LVKGJi/3y6xQ3SGL43grE4mBuao7ZNH5Usbt23uncZ7RdqMix/GiiUTmbNUdsrhu+qre3l5tdyDVefZKzRxclCHLYm0A4lLVHbIgwhrL8ptVd8iD71ddIItts+kAdI3qDlkQ4XxEiqrukIWIbjh6lFWq7pBH39krNbPmogwRWS8RgbbX2iPCdwB8T3WHLIzBF1U3yBIK1R7K58ceUt0hi+/jcwCWtmsSAGBTKDSu7ZoSnWev1Myai5Nl1lwYhlFqZs2FoSlzWqQMcZ5p5zxzteoOWTjPvMdx+hap7pBF563b4/G9MznP3Ka6QxbHyVyq83onzjMb4vG9Vao7ZNF59krNHFyUISLRCCDaVXfIIzoYYw2qK2QhYitVN8iTrQKg81RXyIIoFhCRxuuBxIUA2YDqCln0nr3SMmsuylAoNNE9Nlb1kuoOeex/BggcUV0hC2P+2tf9wha6EQBO30WCCD+AlfgsAEAoNL6vUKj+jOokWYQofLuiotJX3SGL79sf6e9v0nY912/MnnFCZs3FyTJrLoxT1Rb6FgBcpTrjD/D/YBV+XnWEEmbNhaEpc1qkDHGeWeG63qdUd8jiON5Gne8BwHn6SdUNsrjunibOM19V3SEL595Nrpt5t+oOWRwn/dj27b3TVHfIovPslZo5uChDRMImAm2vRUeEkGX52p7yI8Ja1Q2y5PMTDIC0vTcFAFQIQUHVEbIgUo1tM23fEdd59kpN2x+CkjOnRYxTlTktcvoyp0UMTZl3LsrQtm2DoWQyM0N1hyyp1GC9zlsQ67y1+ZYtZLvu7lmqO2TZvr13Wio1VKO6Q5ZUyptDRNr+uaLz7JWatj8ExomFQoXlQogPqe6QJZ8vbBgbs2OqO2RBpG+obpCltrZ/LpGt7W2tKyvt63K5sStVd8iSz8MD3d392p460Hn2Sk3b89LGiRHRYQDsV90hCxHb7ft0VHWHLIjUrbpBlkCAxokopbpDHraHMRpTXSEPuoFAdUF1hSw6z16pmTUXJ8usuTBOVWbNxenLrLkwNGVOi5ShRKK3kfOMtjt09vT0L9mxY4+2O3Q6jrdSdYMs8fjeKs77td2hs6dnYEE83qftDp2Ok7kwHo9ru0OnzrNXaubgogwh2u0AQtt7i/g+/kVlZV7be4sAMI3vb5CdCcC0vbeI7/uXWpal7b1FEMUGgBna3ltE79krLbPmogxZFhvwfXhedYcsjOHWXM7fq7pDFkSh7ekty2LHAEjbjYoQwSGiCdUdshDBtxsamLbPT+fZKzWz5uJkmTUXxqnKrLk4fZk1F4amzGmRMpRMZtocJ3Op6g5Zksn0Gp3Pa3PuvV91gyyO403n3LtGdYcsrustdxzv9L3J3BtwHO8Gz/O03f1X59krNXNapAz5vmhGhBUA8KzqFhmEoJWBADsMAAOqW2QgYtcBwP2qO2TwfTHNtq3LAeBx1S0yEEEMAI8BAFfdIgMiXDU8LP4TAMZVt8jwG7O3lb4EBKfvX2QIvgiX4A4Zn9ocXJQhokIKMXhAdYcslkVfHxsLvqy6Qx5xh+oCeUIHAcbvVV0hi2VZz05M5POqO2QhYncBDGu8j8evzR7BZQBwOl95900AMAcXRml0di46AADaHlx0dLQlVTfIFIu1blXdIMuyZXPHAOBF1R2ydHQ096lukCkWa3lBdYNMOs9eqZk1F2Uomcws5dy7VXWHLK6buT2ZHOxQ3SGL46QfVN0gSyrlzeE8fafqDllcN32V63pvU90hi+tm7tH53ik6z16pmXcuypDvi1pE0HaFN5GYByA0vm03LlZdIEs2K4K2bYVVd8hCRI1EqO1N9YjEgmx2xFLdIY++s1dq5uCiDGWz9o6aGuGq7pAlELDvCgbHR1V3yEKE16pukGVkpG1vfX3fR1V3yDI+XnistrZOqO6QJRCAtZFI24jqDll0nr1SMwcXZeiCC+ZnASD76i/8iBrAgufUFZVAAN4MK3AEAKC9ff4h1TkydXY2a7tB2KpVWACAodf9YpwC0H8a78nTBXlAJACA889fdEx1jkzt7a37VDfIpPPslZqWBxefuPnm2lE/tAYRc4HRwDP3bL4n+395jK4cx1uJCKui0dbiymcbTu8VzwB5ePWtWM69uxDhu5FI63aVSbI4Tua5WKzlYtUdMsTjfc22bX0+Gg1f/+ovjsBRaISQwqw/zE+gEwAcAADOvVuJ8FgsFtZyQzvOve8VCuL6ZcsWaHlXYp1nr9S0O7hYu3ZtYCwX+CkA9ALAUL5qYh0AXPr7PkZnjEFWCDisukMedkgI8dp19kQIm0/jxcs7gWAjvvpWOiJp+7cnxgIFAKHtlUxEeIwxdlx1hyxEuL+mJq/taR+dZ6/UtDu4COaD1wBSZtPDm7oAANbfuM5Zf/361vsevc/7fR6js0ikdQdIurb5VBCNttz1ul94Dm6ERnhIUc4f7mL4GQBcUPwwGg1re9536m3nj6jukEXXdyyKYrHwLaobZNJ59krt9P3b3IkQtAPha4sVCX+Btlj6ez9GY7t27apNJAbmqu6QhfPBea67T9urRVy3/0zVDbKkUqmg66ZbVXfIkkj0NqZSg/WqO2Tp7h5c+MQTpO3VIjrPXqlp984FEM1GfO2On4j0CgDM+X0f47retUSwUgj6Smdnm3uOlbMOCevgPmG9TAB0BvObfCJ/P9l7K5AqGlDMGgccOyTYcC2IabWM6o4RO3yc8HgDEzMrgEIHBdufA8zNYYV5DJDtFdagDWTPYuKMHOHEQWJDIaCqGUw0jBKOHCV2ZDqKGVVINYcFDmeBjTUyf3YAILjPt/YKBH8u8+f7BP5+sn6zA6m2FsX0EWJHRgAPcZ5ZCyCWTkzAVsb8N3HuVQDAjOg++JRtQ88sLMzNEeYOkrU/hBSagf7MUWIjR4kdqUMxvRpF7WGyhrOEYzOZmBUEUTFE9t4CgX8GFuYLRH+/sPYGkYIz0Z89Tpg9RNbBGqSaaejPGCF2ZITYyAwUDSEUVQeFNZQDnJjN/DMsIHuvsAcZkDWH+XPzgLkDwtpfCRSqZ/7MMcLjR8g6XMeorhr8aa91+LOCQBU32cfrPg9wxHXTm4gK5wIUbo/H+/baNvubTWNH6cGJOv4rHUdHiB2bwUR9CET1MFlDE4QTs9A/w0ayXxH2HkTAOViYlweWOyDY/kqEynosNI4RO36E2OFpKOpqUEw7SuzQKLHRBhSNFSgqD5D9Sp6gcAYrNBEg7RPWywGAQCMrzJkgHB8m60ANipppKF7tmI6ivgpF9bCwDkwAjs9i/hwbKFDsONceb3jAzdwRibR8lvNMuxDwLOeZDdFoy6Oum7ni4yOji14oVPWOEo42oN9YgVR5QFj78oD5OcyfBwAw2UGBRua/2lGNVF2Hfv1xsI4eE3hsOvr1VUjVh8g6ME443oj+nABSYB/ZLxMBncEKTQXA/JCw9lUAVTYwv3GM2OgRYodqUUyrRTH1M8+mfuZF6ICw9+UBECGNDAAAH0NJREFU8nOYPw+B8BVh7wkg2I1YOCMHbPygYAfOtcarH+Te/UTwVDYrHNuGf+fc22lZtKmjoy05l+UPAmDV6zoIC0NkvTL5M+/PygIbPSx+tcM6fJzweD2KmZUoQgfJ2p8jzM1m/lwGwF4R1p7J2fPPyAFOHBTWUAipagb6DaNgHTsq8Oh09GdUIdUcEtbBccDs5OxRcJ+w9gpAfy4rzPcBC/uF9UoQqGIm82dliY0dJjZci6K2FsX0Y2QdPu7j2NQ9KToBaDCX8/s59yIArC56CB602dTsAZs4KNhrHb8xe2w4S+w3Zm+yA/z9wn5t9oCNHRJsuAZE7TQ29RrwhrMH1hxWeHX23vg1wJ9Vj754dMeehuXLm4Y5z/wLkX9pW1v/0kSifzpj+Ndfzs7Ah8frsNkqnNHECgv2CpZO+8E9Z9sTZ05HatxVCPQcIuvoMntiWRCpcnuucpuNwloWyC3PEhx3CpV8FvPrF1r5s4fJeuWXhUDfAisfns38prRvv7RX2ENLrPGOGgZ1bj7QPQpW9nx7/AJC9HfkK3ZMQ1F9jp2LjRA71FMI7mxihTOarcKCVwTLeH5w8CxrYvEMRrN+WQj8fJisI8sCE0uDQKEdhcptFgC7ueJI+NbQsS3RaMvnXHfgHCHoh46T/ngsFv666w5c+d5jE6HdfrAEs3eC1wDfehkA4AzLn3ei2ZvGaFoN+HVHyTr0+o7XZu9EHe+sOD7vDu7dj2j9IBJp/l4p/yjW7uCCEJIANK/4sQCot3x48vd9TD4vXggG7ZeIRAYA4JGqV64NBLBqeLjZXbkS/J6elzuJcCISaf55IjVU4/vZM4noUCQS9ranvDm+z+YR5QcjkYVDW3sGFhDRdKKKXZHInNFne/qXCMHsSCSccFOpoO/XdBD5xyORtl/GU4P1vi9aiaz9kUjTnp/27J5PFJgFAP1LlrQc/p+k9yYAVp3LHUguW7Ys/9/JzFIAlluyZH5PsQMADi9Z0tK/rad/NpHVhOjv6eho28+SmTYA1l0oIDJGPmPsFwAiCG/DQ24qtcz3azoAxOiSJa27diQzMwCgrdjxE3dPE6I/+9UOt/9MRKvGso73tLe3555JZpYi+vmOjrbk/7r7qhEn3oSIRzo6mvtecHfPQgzML3ZsddOtiFhvWaFftrfPOv5scrADQAQ7OpoT3d0JOxhsXJLP09jSpeFf/MzxplsWW4CYH+roWDj4HB+cx5iYY1nMa2+ff+hHvG8xY3ZtKJTbN3VbzX8lwpeJ/GGA0EHG8g/cXnvsyIMXTO/7SaK30baDzZYlXm5vb933Y8cLWxZrCAb/f3v3Hxx1fe4L/P18dsOPoIhCao2I+INeb5XyQ+bo8dy5raeWU9vS41hpHaBkE7kpZL8rlmtt53amXefoFE+vv9jdwKVCNuH2dqiX22tLRXs8FY/3nGpPJIuKpSCaSEWBIGCYIMl+P8/9A1ZiSqyajZ98P7xfM8yQTXTezzx5lk++z3c3R3dcdtllXf+09bUrRHTk1KmT2jZvhhl/zp5ppRz/r3XXWSNGxC/t6SnunzXrktc2b3mtOh7X844ejbVfddXEA/98IseZZ9ptF1100TubCu0zrLV25syLtz71b7tHn3GG/XRvrzl85ZUXvFzK0dOje2bNmvzGkydydHf37Lz66ilv/6bQcbkxGDV16qS2hx+GXHbZvunF4vFX97zzTs/u0aMrNoah/RcACEN57r+PPVDX0/N2x6xZ1Z1PtO2aEovFx76b48QvyZox46LCjldfHdXVZS4Pw+LbM2ZcsvOp1j0TRozovbBYlDdmzpy057dbXrkwHo9N6Okpvjxr1iWHH9/S8el4HKO3b59UmDsX+sILe2ZYi3emT79w21PP7BxbWTliShjaAzNmXNT+ZGv7eSNGSHWx2PPazJlT9v/mud2XVlTYs44cMS9dc80FRx/d8so0Y4yZPn1yWymHtcWu6dMv2ZF99k/jjTGTi0V5AxjZA/T+kzHy62PHih0A8Ksz9vxdn9mzL7ywZ0Zp9rZu335mT8/oT5VylGavlOOpE7PX01Pxh1mzqrs3bXllWixmYtOmTd7yh507Rx49OuKK/rNnrXlz+vQLXi/NXhjaXTNmXHSo7+xdeeWVxRde2DOzNHtt/Z4DSrNXeg4wz3f0AmhVRTWAg7GYfVEVFfhbea40e6Ucf2n2fnsiR2n2+j4HlGav9Bzwu37PAU8933ExgLP7z95nPnPhc79vba0YMaLqM6XngH/rN3ulHKXZ6/Mc0AUAIvqT479OfnSvtdgfj/eu/q9nvH0o/9fjdv3qxHNAafZ+duI5oDR7a07O3pbNm01s/Nh9787e4ydylGZvw4nngNLstZx8DnhxypQpx1b1mb3ftu6pHDGi9z+WZu+X/Wbvf/V7Dnioz+xt3gwzfszb00qz19vbu7uiIraxtzd8+vjsoXXdWXu/XsrRf/YeK7TPUFWdMeOiQuk5oP/slXKUZq//c8D27ZMKAHDZZXum95+9YjHsnDnz4o7fnpi9np6Kjlmzqjv7z96mtleni4j0nb1Sju+9O3t4Y/D/+r5XdF/eNYCgJvgbiN6VzeeuXTJvydlmhNkCi1n7ju47NOHMCReuXLvylYG+JteSO+A6PxERUdR5d8/F3u69zwC6N5VI/j42wvzBKO7PteQOjB81fnzMml3v9zWus39cWlt3TSoUOv6T6xxDpa3t1c/5fE9JW1v7PNcZhsozz+wcWyh0fMV1jqHS1vbq9C1bOj7tOsdQ2bq1/Ws7d+709h1IfZ69cvNuLfLwww+HAG4O6oLqSnR3/ePatV0AsGrdqn04caVmoK85XcRi5mLAXgucvO/EJyK43pjwHQCevmxMvgXAy1cdxONmHKDfAN67pvSFCK4G9G0AL7nOMhRUtebwYfMEgGOuswwNf2ev3Lw7XJRk12b/4j8sH+RrfBSG9pV4PF50nWOoqGKTtbHXXOcYOvo/XCcYKsWiPRSPx9e7zjFUVPGMqvS4zjFURKT5rLP6vMeMd/ydPSIiIiKi4aVQ6Li8UOj4uuscQ6VQ6FjQ1rZriuscQ6VQePVO1xmGSmvrngmFQkfKdY6h0tbW8QWf73cqFDq+29q6p9J1jqHi8+yVm7drERqYtXIQwCuucwwVEbtTRL383QYAIGJaXWcYKqNHx9/p7e150XWOoWKt/gmI9brOMXRsYfTog96uXH2ePSIiIiIiIiIiIiIiIiIiIiJyz7u3/46SWxc2XGsN/h6KyyESQvUFFXleK/UXjY2NR1znGyyf6/O5NoD1uc43WD7X53NtPuHhwoEgEUwV6P2q+oaIPA/odqty2IhebCFTBLgaKuv3du9dc+LdRCPF5/p8rg1gfaxv+PK5tr58OTzxpaguCKb1FsN5J96SvK9/Kf0lWZv8yrlnnnsuovgW1j7X53NtAOsD6xu2fK4NJw9PVvUNgTwP0X+2isNGcLECl5lueSSoCSJzeOKVC8eCuqAaIa4HsF1UujMtmQIAdZ2rXHyuz+faANYXdT7X52NtQW2woFgs/mbVulX7BqovWZv8iohsicKvruDhwqEgEUwFdCMEBahsheIMFT0rl8/d4jpbOfhcn8+1Aawv6nyuz+faAH/q8+5XrkeJqC6Byt1W9McAkG3OLhPg6sULF5/vOls5+Fyfz7UBrC/qfK7P59oAf+rjPRcOqUhcRfe+50HB4dEVo486ilRWPtfnc20A64s6n+vzuTbAn/q4FnGoIdEw3UA2iKBVrZ4DwRGI7Ms25b7lOls5+Fyfz7UBrC/qfK7P59oAf+rj4cKxO+rqzuy2lXMATBWEj2TyK59FxG9M6svn+nyuDWB9UedzfT7XBvhRHw8XDqVSqZHapfdNmDxhaeerncsguhAi92absk2us5XD0kVLzw2L4dn9H8/ms9td5Ckn9i7a2D8arnzpHe+5cEiP6FwABw50HDgbovUaw+clRHM6nV6XTqcj/2uLbW/vAoi56fhHOhbAJQL8EcA0l7nKgb2LNvYvunw/GPrSOx4uHBLV2dbiLjX6RQM8ll2T6whqk28dan9zIoB21/kGK9PceC+Ae0sfB3VBNayucRipbNi7aGP/osv3g6EvveNLUV1SfUJiWC5AGjDrUzUNswGd/EB+VbvraEMhuza7R4FPpOanxrrOMmjsXbSxf5ElqrPV6v+0ar8owGO5NbkOiJYOht6Jau945cKhTHNjy60LG3ariRUz+czTQW1Qq6He7DpXuQS1wQKoXnfyETkP0O2Zn2bedpeqPNi7aGP/Ikz1CYnJciimASaRqmmYrR4dDH3pHQ8Xjq1oaXyy9HdfdoYlEso2NTj27sdGDoaj7DMuM5UTexdt7F80+X4w9KV3PFw4kKptCNSaawb6fLY5O+/jzDNkpFgdh/67Lz9RAOxd1LF/0dW3dxYAVBHUBEtgARH5AgD2bhjh4cIFiT0N1V2uYww1K2aKqn4OwHdcZykb9i7a2L/oYu8ihe9z4Vh9fX3FyN6R404+YoOYVtz/QP6BQ+5SlUeqpmG+imRF8aJCdpce9+WnQ/Yu2ti/6GLvhj9euXAoWZu8Snp0nYX9pACHAEwE8OS4yePudJ2tLEysDaHOi9Tbyn1A7F20sX/Rxd5FAw8Xbs1TlUVGcKkVhAixwxh8N51OW9fBysHCftbG7eMr1658pfRYkAh+VF9f/9Tq1au7XWYrA/Yu2ti/6GLvIoCHC4eMRaU1OCbQ7UZlbqYl1xwkgjCoC6qza7N7XOf7qFKLUhO1aB+FojpmzbIgkSz9Nj8D6Kjq6urvOw1YBuxdtLF/0cXeRQMPFy7FzHqxdoXEK75qe8PrUonkcgu9vGpS1T7X0QZj/MTxb+7fvX+OCXGbBTYihpdLn6u6oOp1L37CYO+ijf2LLvYuEnhDp2Op+amxXRVdPWPsmEli9EaR2C8zTZmXXOeiv4y9izb2L7rYu+GPhwsHkrXJB2Fwn4Ryg4F5dEV+xU7XmT4OS+qWXByzZpcA98Ca9ZmWTJvrTB8Wexfd3gHsX5T7x95Fq3dcizggVrZBtQCosQi/HSSSx/p+PpvP/QdX2cotqAuqEeJ6ANsBHIvZ2IUPtjz4WtTeJ7/kdOpdX51dnR1VZ1RNamxq3B3V3gHsX5T7x95Fq3e8cuHI7d+8fcyx+NE7oHaTteY9J/BcS+6Aq1zlFCSCqYBuhKAAla1QnKGiZ+XyuVtcZxuM06F3wHsPhqLSnWnJFABE/RVy7F+EsXfRwcMFDZlUTbJRIQUbsy8Za2Zn89kfBInktqINZ69qWfW663w0MF8PhqcL9i+6fOkd1yLDSFR3awNRkbiK7n3Pg4LDoytGHx3gP4ks33onqksUcrc19iWjMjvbnF0WJJLbFi9cfL6PB0P2L7rYu+HJuA5wugvqguqgJrglqAn+JobYOTEbvzCTz30PISL/HvoWtlEU98VUUlD710Gi4RcAXrh/zf1vuc5Wbp1dnR1WdJIvvVORuJrT42AIsH9Rxt4NT1yLOOTL5a/3c0dd3ZndtnIOgKmC8JFMfuWziNjucCA+7EUH0pBomG4gG0TQqlbPgeAIRPZlm3Lfcp2tXNi/6GLvhj8eLhziPQnRxYNhtLF/0cXeRQPXIg75cvlrIKlUamSQCHLpdDoe1AR3BInki0FtUOs6VzmI6hKo3G1FfwwA2ebsMgGuXrxw8fmus5XD0kVLz+22lecD2AKgWRE7FCQCb17q53v/UqnUyG5buXzC5Ak/h+KgwjwU1AYJ17nKwffe+TJ7PFw45Ps9CXpE5wI4cKDjwNkQrdcYvgxFTTqdjvyNxL4fDG1v7wIATcf/6AZACwJd7zpXufjeP85edPkye5H/Rouyxnxj4Y66uundtnIORI5f/mpqfNZ1rnIR1dnW4i41+kUDPJZdk+sIapNvHWp/cyKAdtf5BsPCNhqVDQbSqmrP8e1gmGluvBfAvaWPg7qgGlbXOIxUVr73j7MXXb7MHu+5oCGTqmlYqEZugGKawCSg4Wg1WJ5tapzpOls5+LAX/TCCRPI56TXXZn6aedt1lnLwuX+cPb9EcfZ4uHAolUqN1C69b8LkCUs7X+1cBtGFELk325Rtcp2tXG5d2HCtmlgxk888HdQGtRrqv+Zacjtc5xqspYuWnhsWw7P7P57NZ7e7yFNuQW2wAKrXnXxEzgO0M5vPzXeXqnw4e9HF2YsGrkUcOsVe9PMSojmdTq9Lp9NF1/k+qlRtQ6DWXAMAFgBUEdQES2ABEfkCgHlOA5aB7e1dADE3Hf9IxwK4RIA/ApjmMle5CHQygG0KaQcAMXIwHGWfcRqqjDh70cXZiwYeLhzydi8qsaehGvk3s3k/vuxFB2IVR0T1/Gxz7seuswwFzl50cfaigWsRh3zfiwJAfX19xcjekeNOPmKDmFbc/0D+gUPuUg2NKO5FB5KqaZivIllRvKiQ3aXHs83ZyP/kC3D23KUaGpy94YdXLhzKNDe23LqwYXe/vejNrnOVS7I2eZX06DoL+0kBDgGYCODJcZPH3ek622ANsBfd7sOTGwDAxNoQ6jxf75Dj7EUXZy8aeLhw4HTYi54wT1UWGcGlVhAixA5j8N10Om1dBxssX/aiA7Gwn7Vx+/jKtStfKT0WJIIf1dfXP7V69epul9kGg7PH2RvufJk9Hi5cOA32ogBgLCqtwTGBbjcqczMtueYgEYRBXVCdXZvd4zrfYPiyF+0vtSg1UYv2USiqY9YsCxLJ0hsTGUBHVVdXf99pwMHi7HH2hinfZo+HCwcyazNbAWwFTr0XvS1x2++82IvGzHqxdoXEK75qe8PrUonkcgu9vGpS1T7X0QbLqO5XkR+mapJXR3kv2t/4iePf3L97/xwT4jYLbEQML5c+V3VB1etR/8mXs8fZG658mz3e0OlQsjZ5laiuU8h79qITJld9IWrfSANJzU+N7aro6hljx0wSozeKxH6Zacq85DrXYKVqU5/WUC/s/3i2JbvJRR76cDh70cXZiwZeuXDLy71osjb5IAzuk1BuEMij+Xx+J4AdAJa7zlYuvuxFP6gldUsujlmzS4B7YM36TEumzXWmQeLsRRRnLxqzx8OFQ77uRcXKNqgWADUW4beDRPJY389n87nI/Ya/Et/2oh9UZ1dnR9UZVZMamxp3p+anxrrOM1icvejh7EVr9ni4cMnTvWi2Obv69m/e/tNj8aN3QO0ma81O15nKxbe96PsJ6oJqhLgewHZR6c40ZQoA4MVL/jh7kcPZi9bs8Z4Lx3zdi1K0BYlgKqAbIShAZSsUZ6joWbl87hbX2cqFs0fDkS+zx8OFA333ogbm0RX5Fd78dPF+oro7/CB8qy1Vk2xUSMHG7EvGmtnZfPYHQSK5rWjD2ataVr3uOt9Hxdnz4/uzL99q82X2uBZxwOe96PuJ6u7wg/CtNhWJq+je9zwoODy6YvTRAf6TSODs+fH92Zdvtfkye7xy4cjt37x9zEB70VxL7oCrXOX2Z7vDlkwBQNTf2RaA37U1JBqmG8gGEbSq1XMgOAKRfdmm3LdcZxsszl70+VybL7PHwwUNGV92h6fic20ld9TVndltK+cAmCoIH8nkVz4LT57Afefz96fPtZX4MHvGdQA6aUndkouDRFJTieTy1MLUDNd5BktUl0Dlbiv6YwDINmeXCXD14oWLz3edbbB8rg0AUqnUyG5buXzC5Ak/h+KgwjwU1AYJ17mGCmcvOnyuDQCWLlp6bretPB/AFgDNitihIBFEbl3Hw8Uw0tnV2WFFJ2Xyue8hROR//4GKxNVEf3d4Kj7XBgB6ROcCOHCg48DZEK3XGL4MRU06nfbyPi3OXnT4XBsA2N7eBQCajv/RDYAWBLreda4Pi2sRx7g7jCafawOAVCLZYi3ugsFVBrgqk88FQW3y/8Q1XPZAflW763zlwNmLJp9rO5WgLqgWq2sy+dz1rrN8GDxcOMTdYbT5XFuqpmGhGrkBimkCk4CGo9VgebapcabrbOXA2Ys2n2s7lSCRfE56zbVRehMtrkUc8n136PPe3ufaACDT3NhiQs0ITCKTzzytxpyvodzsOle5+D57vuztT8Xn2gAgqA0WBIlk/uSf4HEA26N0sAD4PhdO+fJ65oGcYm//eQnRnE6n16XT6aLrfIPhc20lK1oanyz9PduUbXKZpdx8nz3b27sAYm46/pGOBXCJAH8EMM1lrnLwuTYAkFC2qcG7778iRg6Go+wzLjN9FDxcOGRhG43KBgNpVbXnBImGXwB44f4197/lOls5iOpsa3GXGv2iAR7Lrsl1BLXJtw61vzkRQLvrfIPha22p2oZArblmoM9nm7PzPs48Q8X32cs0N94L4N7Sx0FdUA2raxxGKhufawMASLE6Dv33qN/bxLWIQ435xkKlOTpdVR6BmFaB3pNtyi12natsVJ+QGJYLkAbM+lRNw2xAJ0d9aAD4W5vEnoZi3YB/POH97PWTXZvdo8AnfHgHy/58q82KmVJUk3SdY7B4Q6dDqVRqpHbpfRMmT1ja+WrnMoguhMi9Pl2CvnVhw7VqYsVMPvN0UBvUaqj/mmvJ7XCdqxx8rg0A6uvrK0b2jhx38hEbxLTi/gfyDxxyl6o8li5aem5YDM/u/3g2n93uIk+5BbXBAqhed/IROQ/Qzmw+N99dqvLwuTYASNU0zFeRrCheVMju0uNRu2rItYhD3NtHm8+1JWuTV0mPrrOwnxTgEICJAJ4cN3ncna6zlQP39tHlc20AABNrQ6jzov7SFx4uHOLePnp8rq2feaqyyAgutYIQIXYYg++m02nrOlg5cG8fYT7XBsDCftbG7eMr1658pfRYkAh+VF9f/9Tq1au7XWb7MHi4cEn1CYnJciimASaRqmmYrd7s7TXy73J4Sj7X1oexqLQGxwS63ajMzbTkmoNEEAZ1QXV2bXaP63zlll2b3RMkkp9IzU+NjdpL/k7Fipmiqp8D8B3XWcrN19pSi1ITtWgfhaI6Zs2yIJEsvXLJADqqurr6+04Dfki858Ix7u2jy+faUnWp69Tau2Px+FfD3vAJEf21BW6omlx1hQ8rO+7to8vX2urr6ysrwooqE+I2q9iIOF4ufa7qgqrX0+l0MZ1OG2xDPP1wusdh1A+EVy4c494+mnyuDQAyazNPpOanfn8Yh3vG6JivQfRGI7EbfThYANzbR5qntY3oGXGHQEeFMc01rml8ue/nkguT41OJhhs72zv/PnZW7BYAewf43wwbPFw4wL29F3t7L2tL1iYfhMF9EsoNAnk0n8/vBLADwHLX2cqKe/vI8rW2bD6bTi1M/ZWB5oJE8jMA9kNwBIpPQXFQ1fzkiBy5Kf9Q/h3XWT8IHi5c4N4+8nt7X2sTK9ugWgDUWITfDhLJY30/n83nvHibZe7to8fn2koyLZnfA/g7AHJbYvGFPWImxmKx9sxDmT+5zvZh8XDhQGZtZiuArcCp9/a3JW77nQ97e8TMerF2hcQrvmp7w+tSieRyC728alLVPtfRBs3T2iZ0T8gfGXXkp8fiR++A2k3Wmp39vyZKe9+BGNX9KvLDVE3yap/29sfMsbcqYhVzTIjbLLARsT/b29uo9s/n2k5BT1xVa3ec4yPj4cIh7u2jy9faOsd0/jeBjgqNPtS4ZuWp9r7/JUp73wFxbx+5/vlcm4/4ahGHkrXJB2Gxof/ePpPP3eA622D03dsbmEdX5Ff82U+/UeVzbSWpham/UqP/AOif7X0B+ckRObIin4/G3ncgydrkEit/vrfvGdHzD1He2wN+98/n2nzDw4VDqZrkT6zBQzFoTFXmZvK5bweJYAMMUlHe2wc1QT1E78Hx311zEIA3e3ufazuF43tfYybGTDT3vv29u7cHqqE4CEGfvT1GTZhc9amo35Tbh3f968Pn2rzAtYhL3NtHbjfqc22nEPm9b3/c23vD59q8wMOFQ9zbR2836nNtpwPu7Yk+HlyLOMC9fbR3oz7Xdjpg/4iGHg8XDnBv7w2fazsdsH9EQ4SHCwfSc9Mjjow6UjHQ3j7Xkjvg0d6XiIhOM7znwgHu7YmIyGe8cuEI975EROQrHi7c496XiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIhoGBLXAYjo9Pb1L3+9Bmq/FItVJH/2q591AsA3vjT3R1D0rN/08A9d5yOiD8+4DkBEp7feEb0bAfnPNiw+CADf+PJNNyvwHcBscp2NiD4aXrkgIue+8aWbvqKQX0HlFojeI4pGXrUgiq6Y6wBERNt2vrTjiimXXwBBGsC2fUf317S3t1vXuYjoo+FahIiGBRXdCwAQ7N+8eXPRcRwiGgReuSAi527+0teuUpi1EL0bKrWXX3rFG9te3rbFdS4i+mh45YKInJozZ06lhWkB8H9//uv//QNA/hFG75s7Z+5FrrMR0UfDwwUROVVpR90D4BPFMEwBQJd23QnFboRoTqfTfI4iIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIi8tb/B8wWqZPxwJhwAAAAAElFTkSuQmCC",
"prompt_number": 8,
"svg": [
"\n",
"\n"
],
"text": [
"Plot(...)"
]
}
],
"prompt_number": 8
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Marginalization"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"prior = @pp enum() do\n",
" local triple = model()\n",
" triple[1] # rain\n",
"end"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 9,
"text": [
"Discrete({false,true}, [0.7,0.3])"
]
}
],
"prompt_number": 9
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"bar(prior)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"\n",
"\n"
],
"metadata": {},
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAAhcAAAF6CAYAAACqW3pRAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deWBcVd3/8c+5szSTpZTuFGgmUNwCzEypgMsPqSIKUvdY4BFoi9TSpsqq6KNYQcQqq21lFQq4UIK4AEVwKeJWMc3cKQaULjMpWJq0pYUsk2Rm7vn9UYsL4KPchntT3q+/mOQ28w4l7ZdzzsyVAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA8BqxapWN5nLrxwfdAcAfJ+gAANitrm7jJGujVwbdAcAfhgsAoZFIOIOSCkF3APDHBB0wVM5tOjdRqimdtHT50paX+vxn5syp660kTjDGDMZ6Yw9e3XJ18dVuBABgb7RXrlzMnz3/wFLNwOUy9tyX+vzcuXNjfV7i15KaJL27VD3wk1e3EMBLaW3dXO26G48MugOAP3vlcOF45luSc9TLfT5eis+UUcey25Y1LV2+dIGMGdt8WnPDq9kI4KUUx0rOwqArAPizVw4XS29bNsN69ryXvcCqUdbk/v7YPGGi3hGvRhuAlxeJOM9L9r6gOwD4Ew06IBDWTjBGv9n90Bj7jKSJARYBkJTJNOyUtCLoDgD+vCaHC2u0VrL7737sSaMjFf3T/y25bsdnJdvkefb8qVMbfuW6hTskvTEWixxfLHZ2R6Njf2+t/Wsm0/AB1y1kJN0k6YF0OvnFbDZ/hjFmoWS/kU43rMhmC98wRtM9z5w+dWr9467b8aCkUel0/VFtbZsmOY73E0luOp38hOvm3yOZy4wxd6RS9dfmcoULrNXJxugzqVTyl7lcx3Jr7aGeVzohk5myPZfreFTSlnQ6edJjj208vFJxbpHMz9Lp+s+5bsfHJXuOtboqk0l+z3XzX5PMcY4TmX344Qc+ls0WVhqj8el0clout368tdGVknksna6fnct1vMtau1jS99Lp5FXZbP4cY8zHrTWfy2Tqf+a6hW9LSsVi9qTGxoYtrpt/VDLb0+nkCa7b0SjZ26zVqkwmeWEulz/ZWnOBZK5Np+vvyGYLlxmj90g6K51OZl23cK+k/Z58sv6o+vq/jhoxovygpCfS6eRp2Wz+WGPMFdZqRSaT/IbrdiyU7BmOoy8efnjygVyucKO1mipFPpBOH/hX1y38XrLd6XTD8WvX5t/geeY7kh5Jp5PnZbMdTcbYz1prl2YyDctdt3CJpBM9z5k3derkVtct/FDSgbFY71sHB6urjTE/t1ZPZjLJU3O5Tf/PWu9qyd6dTjd8LZstzDdGcySzKJ2uv891O66T7JvL5cqHp007eJPrFn4jaTCdTr4zm91wiDGR70v2t+l0w6ddt/BhSZ83RtenUsmbc7mOL1lrZxhjF6RSDX9w3cLdkpLlcvyY2trnI/39VaskuyGdbpjZ1lZ4i+NoiaQfpdPJr7hu/pOSOctae2km0/Bj1y0skfQWY9SUSiXz2WzHryQpk6l/Ry5XaLBWLZJWp9PJ5mw2/wFjzBcle1M63XCD6xa+IOmD1upTmUzyd67bcadkp1RV9U/v6RlZiUYHH5FUSKeTH911LsL5ljHm3lSq/su5XOET1mqepK+m08l7stn8NcaYt1tbOSWTOXid6xZ+KSmeTiff3tq6YXI0GrlHMn9Mp+vPdt2OkyS7yFrdItnvOY6+Za15nTHOuanU5F/ncoXvWqvXW2uPi8f7+kqlmt9JeiqdTn5o7dqOIzzP3iBpZTqdvDibzc8yxjRbaxZnMvUtrlu4UtI7KhWddsQRySdcN/+QZOrS6eRbXPep/aXKj61VNpNJnrV2beEEz9Olkrktna5fks0WLjRGMyVdmE4nV7lu4XZJbxoYiL6no2P/na97XccfJD2TTidnZLP5tDHmZmv1YCaT/F/X7ThNsp82xl6RSjXcmcsVvm6t3mmMMyuVmvwn1y08INkx6XTDke3t+YmlkrlP0tp0Ojknm+14tzH2cmvtdzKZhmtct3CepFM9z1w0dWr9z12341bJHmZM+cRUakqX6xZarVVXJpM8ce3apw7zvMqtkv1FOt3w2Wy2cKoxOk8y16TT9d9x3Y7LJftuz/POnDr1oJzrFu6TNDGVqj8ym10/xnFiD1hr2zOZhjNctzBd0jeM0Z2pVPKKXK7j09ba0yT7v+l0w4PZbOEmY5TxPOf9U6dO3uy6HX+QtDOdrn/PmjWFN0YiukPSw+l08gLXzc+UzIXW2iWZTMNtrlu4VNIJxmhuKpVsy2bzPzbG7F8ub3tLIjGhrlSqPCTpz+l08uNr1248xvOcqyTTkk7XL3bdQrOkWdbq4kwmudJ1O66X7LRIpPyhww6b8lQuV/idtepLp5PHue6G10mR7xljfp1K1Z+by3V81Fp7kaRvpdPJW7LZ/CJjzEnWOmdnMpP/6LqFeyRN3mcf+/bOznK8qir2S0nr0unkKdlsx9uMsdcao3tSqeRXXTc/TzKfMMZ8OZWqvzeXKyyzVkdZaz+ayTQUXDf/iGQq6XRy+mOPbTq4UvFWSPp9Op1cmM3mP2iM+YK15oZMpv4m1y18UdIHPE8Lp05N/j6bLbQYowZjqt4hSdb2/8pa5TOZZFMulz/aWrNU0o/T6eSl2WzHWcbYT1prv5LJNPzIdQvflPTWSiVy8hFHHLjedQurJBtJpxuOyWbzSWPM3dbq0UwmOT+X65hhrf2SZG9Opxuu9/t362vCgtMXvKV59oLf7X7c1NQUOXvO2QdJUvMZzW9rnrVglSSdferZ+y6YtSC/4PQFY4JqBbBLa+uGyX8b5AEMY3vlmYuXMqZqzJiI52yQpM6+ztWS7Vw4a8GjkbjzhGN19bLbl20PuhFAYpvkLQm6AgBeseY5zZM+M2dOXdAdAAAAAIZALrd+vOvmLwq6A4A/r5ltEQDhVyqZKsk0Bt0BwB+GCwCh0d190GZjyucH3QEAAAAAAPBivBQV2DuwLQIgNGKxEZ5keoPuAAAAAAAAeLH29vZ4LlfgDsXAMMe2CIDQKBarJlqrS4LuAOAPwwWAEEn0SebRoCsAAAAAAABerL29q/Zvt/oGMIyxLQIgNIrF7tGS5gTdAcAfhgsAoVEuezslsyLoDgAAAAAAgBdrbd081nU7FgbdAcAftkUAhEixWrJHBl0BwB+GCwChMWaM02WMszjoDgAAAAAAgBdra9s0yXULVwXdAcAftkUAhIbnlaKSxgXdAQAAAAAA8GLWWieX21ITdAcAf9gWARAaa9ZsPMDa/uuD7gDgD8MFgNBIJJxBSYWgOwAAAAAAAF6stXVztetu5B06gWGObREAIVIcKzncWwQY5hguAIRGIhHvkbQq6A4AAAAAAIAXy2bzo1w3PzPoDgD+sC0CIDQqFW+kZE4KugOAPwwXAEIksU3ylgRdAQAAAAAA8GK53Prxrpu/KOgOAP6wLQIgNEolUyWZxqA7APjDcAEgNLq7D9psTPn8oDsAAAAAAABerLV1w2TXLdwRdAcAf9gWARAasdgITzK9QXcAAAAAAAC8WHt7ezyXKzQE3QHAH7ZFAIRGsVg10VpdEnQHAH8YLgCERixm+yXbHnQHAAAAAADAi7W3d9W6bmF60B0A/GFbBEBoFIvdoyXNCboDgD8MFwBCo1z2dkpmRdAdAAAAAAAAL9be/tRo1y2wLQIMc2yLAAiNYnGwVhIHOoFhjuECQGiMGeN0GeMsDroDAAAAAADgxdraNk1y3cJVQXcA8IdtEQCh4XmlqKRxQXcAAIC9hLXWyeW21ATdAcAfVi4AAMAexXABIDTWrNl4gLX91wfdAcAfhgsAoZFIOIOSCkF3AAAAAAAAvFg+n6/K5TYdGnQHAH/YFgEQGtu3e+Ot9T4bdAcAfxguAIRGIhHvkbQq6A4AAAAAAIAXy2bzo1w3PzPoDgD+sC0CIDQqFW+kZE4KugOAPwwXAEIjkah7VtItQXcAAAAAAAC8WC63frzr5i8KugOAP2yLAAiNUslUSaYx6A4A/jBcAAiN7u6DNhtTPj/oDgAAAAAAgBfL5Z4+wHU7uOU6MMyxLQIAAAAAAAAA/6f29vZ4LldoCLoDgD9siwAIjWKxaqK1uiToDgD+MFwACI1YzPZLtj3oDgAAAAAAgBdrb++qdd3C9KA7APgTDTpg2HjY/kpWmaAzgL3QJZpurpCkYrF7dDQamSNpVcBNAHxguPjP1UqqCzoC2AuN2P0P5bK3MxqNrggyBoB/DBcAQuPoow95XtJ9QXcA8IdXiwAIjfb2p0a7bmFO0B0A/GG4ABAaxeJgrSQOdALDHMMFgNAYM8bpMsZZHHQHAH84cwEgNBoaGvol/SnoDgD+sHIBIDTa2/MTXbdwadAdAPxhuAAQGsWiF5eUDLoDgD8MFwBC44gjDnramKp5QXcA8IfhAgAA7FEMFwBCY82ajQdY23990B0A/GG4ABAajhMrS9oadAcAf3gpKoDQmDp18mZJ5wXdAcAfVi4AhEY+n6/K5TYdGnQHAH8YLgCExvbt3nhrvc8G3QHAH4YLAKGRSMR7JK0KugOAP5y5ABAajY0HPivplqA7APjDygWA0Fi9et1I1+04KegOAP4wXAAIjWjUGSXZmUF3APCH4QJAaCQSdWyLAHsBzlwACI3GxvEc6AT2AqxcAAiNXG79eNfNXxR0BwB/GC4AhEapZKok0xh0BwB/GC4AhEYi0b/FGF0cdAcAfzhzASA0GhsbByXlg+4A4A8rFwBCI5d7+gDX7eCW68Awx3ABAAD2KLZFAIRGKnXA05LmBd0BwB9WLgCExqpVNprLrR8fdAcAfxguAIRGXd3GSdZGrwy6A4A/DBcAQiMWs/2SbQ+6A4A/nLkAEBqp1JQuSV8LugOAP6xcAAiN9vauWtctTA+6A4A/DBcAQqNY7B4taU7QHQD8YbgAEBqRiPO8ZO8LugOAP5y5ABAamUzDTkkrgu4A4A8rFwBCo739qdGuW2BbBBjmGC4AhEaxOFgriQOdwDDHcAEgNMaMcbqMcRYH3QHAH85cAAiNhoaGfkl/CroDgD+sXAAIjfb2/ETXLVwadAcAfxguAIRGsejFJSWD7gDgD8MFgNA44oiDnjamiluuA8PckJ25mHnCzIPlVGZYzxxjjSYZacAYrfWkxyLWy9258gePSrJD9fwAACAYe3y4aJrR1GDKutIa74OyZoekvxiZx2XVY4092Fgd78lJNr2vaZ2kz7fc3/KTPd0AYHhas2bjAdFo5DJJpwXdAuCV26PDxcdO+NhZqtgL5JibrLXzW1a2bHmp64499tjohJqx77cyFzad2HRWy8qWGXuyA8Dw5DixsuRtDboDgD97dLiw1j7Ro57DH1j5wMC/u+7hhx8uS7pH0j1NJzS9Y082ABi+pk6dvFnSeUF3APBnjx7obPlpy28eeODfDxaStGjRoheet+WBll/tyQYAw1c+n6/K5TYdGnQHAH+G7NUiM09o+vLME2Ye/a8fb5rR1ND+aPvDQ/W8AIav7du98dZ6nw26A4A/Q/dSVKO4Nd5vZp7Q9OVjjz02Ku06k2E8rTXS2CF7XgDDViIR75G0KugOAP6YofziM9/30ZOtNUsk5WW0VVbvkczXu233l/+T7ZNX6jNz5tT1VhInGGMGY72xB69uubr4r9fMnzX/vY6cut2PB+ODP7rxxhtLL/tFH7ZrZDV1iJKB17IvaLq5LOgIAHvOkN5bZMX9d9/5sRM+Vi1jvy0rydhL77q/5eKhfM65c+fG+gZjv5a0TlJXqXpggaR3/+t1RuY6yd67+/Hg4OC9kl5+uAAw5FavXjeyqip+TDpdf1/QLQBeuSEbLk6ZccpYr1K+1sqeaqRvW2M7Zc3nP3Zi05us1PxyL1P1K16Kz5SxHcuWL2uSpOZZC7LNpzU3LL1jaX73NQtOXzBGRvmlty771FA0AHhlolFnlGRnSmK4AIaxIRsuKpXSMskcY2RnrFh5932SdPKJH/mJJ+d2Iz0had8heWKrRsnk/v7YPGGi3hGSXhguFNEUY7Vv86wFD1mpS1Y3Lrtt2SND0gPgP5ZI1D1bKvXdEnQHAH+GblvEOCtL8dL8H/7wh9t3f+jOlT/4w4wZMzLVXtXiIXteaycYo9+8kGHsM5Im/lOaZ6LG2N9ax3xVFftGY9Qyd+7cxhtvvHHb7mtyufwp1upYz7Pfmjr1oNx4pzKxqxIZsmzgteoNzsDkO938DdZqZWPj+B9ns4X9XDd/QyRilx122EFrXbfwBckeGIvVnF8sbhqIRscslUxXOp38outueJ3knC9pdTrdcGsul3+vtfqQtfp+JtPwcDZbmG+MTZXL3mXTph28yXULVxijRCqVXNDaunlsNDpwmTFmXSqVvMJ1Nx4pmTMl82A6nbwnl8ufbK2mW6vrMpkGN5fLf85aJfv7yxceddSUnlyucJ21zvZMpv7za9Y8NSUSKV9ojPljKpW82XU3Hi+ZjxhjVqRSyV+6bv6TkqYaY76WSiXzrptfLDn7pNP189rbnxpdKpUvt9bZkMnUf72tbdM0x6mcZYzzs1Sq/u5stqPJGO84x3FuPPzw+jWu2/FZyTvIcZyLDj+8fofr5m+QtCOdbrgolys0WGsvktSWTjfc0NbWcZzjeE3WOndnMvU/y2Y7zjLGmxaJRL5+2GGTN7hux+WSNzqVSs5z3cI+xmixMSqkUg2Xu24hI9l5kn6ZTjesyOUKH7HWHm9t5OZMZvIfs9nChcbYKQMDsc8fddQB21234zprbU8mk7ywrW1jveOYzxsjN5VquM51C9Mle7Kke9LphgdzucKZ1tojJe/KdPrgJ7PZwmXG2LE7diQX7LffXxL9/SOusNZsymSSl7W1bUw5jplvjB5OpRq+77qFD0n2vcbo1lSqYbXrFs6T7OuNqXwxlZrS5bqFpZIG0unk+Y89tv7ASiXyBck8lk4nl7a15d/hODrVWvPjTCa5MpvNzzJGb3EcXX344Q1/dt3CJZKdEIv1LSyXx8SsLV4lOX9Np+svyeU2HWptZaG15pFMJvndXG7T+62tvM9ac1smk/xdNps/xxi9MRbTlxobG7a4buGbxphKKlV/blvbpkmOU/mSMU57KlX/TdfteLvknSY596bT9fflch2nW+u9zfOca6dOrX88m80vMkb77bOPPi1Jzz2na63VM5lMw6K2to43OY73aWOc36ZS9be7bsdJkjdDcu5Ip+t/47odCyXvUCl6STp94F9zuY6rrbWRdDr5qfb2/MRSSV+2Vk9kMg3XZLOFtxpjzzAmcn8qNXmPvlv2Hh0ummY0TW65t2WTJN11/123vdQ19957b5+khf/waxpa7m3Jv9S1r4Q1WivZ/Xc/9qTRkco/L7EuvW3pbyX99m8PNy+cteDBEYOx90t64f+YSqURP4vHS38ol2NbJGmHF9kuadKe6gSwy9M21uk4zvJIxNnZ2rp5rDGlyY5jFg8MxP62dRq51XG8Ee3t44pNTeO8xx7btNiYUkmS+vriT9XWeoutLfVIkudpdSTiPNnXN7hNkuJxe0+l4vx0n33Knbs+7y2JRJy/vUpu83OOM25xuWz6JSkWq328UikujkScnZJUKo34eTxeetTaEZ2SZG30dsfxRmzaNKX3qKNkHcdZPDhYLkvS6NHlp7u7nRc6rDWPRiLO+sHB8nZJikTsj6yN/CyRGNy8q8N+KxpVRJK6ug54fsyYTS90jBhR9edKpbi4WIw8J0mDg5FfJhJas7vD88wd0ahTtX375G5J+seO6urS5mIxvtiYSu+ua8t/jEajG3d3OE7px8bEflFbW9ksSeVy+bp4PBo1xti77rLdb3jDpsWe5wxIkjFVTxoz8EJHNBpZVal42UikqmvX9+h8NxKxVX19+z+363r7dc/zPEmqq6ts+ceO/v7SmurqeEel4j276999+V7Hia2qq9PmXV+7fL21sdixx6rS0vL6vn/s8LyqddFoafHgYOz5Xb9PkV9VKl4uEkl07foe9P143EkMDm7fsatDV1YquzocZ6DT2trFkYjXt+vfZTlbXR1/aneHteX7I5H4Iz09zjO7fp/KN1kbi73pTW8qSSo/9timxdZGBnf99xDdGI/bFzqMsb82xvlTPF7cKkmVilkRj5tEsbh1+66vba+ydtfts0aM6N5Wqfy9o6qqmBscTGyWtGPXrx18IBqN/6avb1eH48RuNqYSTyYnD0rSP3b09zv52lotLpcH//b7r99KzuO7O2Ixr6VSce6PRJ7f+rf/5q/d/bNWLG7fHo+PWzw4aIuSlEj0PzY4mFgs2R3aw/boq0Vmnth0uZXeao2+3vjmxgcWLVrkvdy1p7z3I4dXnMg5VvZNLStbXvR+GK9U8xnNb5OxX1m6fNn0s089e18n7rTJ07SuYtfOsXVj66+75bqNzWc0n2mMPXTJ8mXnzp07NxYfjK0x1vnQktuWbHjZL8yrRYCh8sKrRVpbN0yORiOXpdNJ7i0CDGN7dOVixcqWz8084aPHS+bqxx9t//bHTmx6XMaul6f1xlFRMq+zVm+Q9PqKVCfpEtXorD3Z0NnXuXpCzbjOhbMWPGqlycbqq0tuX7Z93mnzxkc8Z4MkU+VV3dkfLf6sefb8lRp00pJ++G8HCwCvikSif0u5XDOkrygDMPSG5H0uFi1a5Dy++vH3ytGRVvYwRzrMSvtL6pTsb430S6/PtLQ83NIzFM8vSc1zmidVq6/767fc0v1y15zz8bn7KVpdvGb5NTv/zy/IygUwVHifC2AvM6RvorVXYbgAhsoLw0Uu9/QB1la+kE7Xzws6CsArN6RvotXU1BRXr2Y60pussZ7x9BfPmIeG6j0uAAxvpdKAE41GaoLuAODPkK1cnPL+UyZUKuXfyGqypL/87cOHStphjU5vub/l/qF67iHBygUwVNgWAfYyQ3bjskq5crmsnitXKsm7VrYcftfKlsNL5fJEGfsDx+rqoXpeAMPXqlU2msutHx90BwB/hu6uqLKHy5gl9zx4zzO7P/LDh37YVYpXPmelA055z0cah+65AQxHdXUbJ1kbvTLoDgD+DOFwYf4gT/v/60ej3dGyJNM7orRp6J4bwHAUi9l+ybYH3QHAnyE70BnxKjdUHOf+phObPMc6DxtrtlScyiFG+pKRvhurxMZ95H0fGSdJh735sMK/e8MtAK8NqdSULklfC7oDgD9DNlx4jnOBpAOMdLk1nqz5++lRK70tYp0zd1+77g/r9pH0/FC1ABgeWls3V0ej/Yem0wc9GnQLgFduyIaLUqXy2YiJ/EcnwA85+pAePTBUJQCGj+JYKbJQEm//DQxjQzZc/O0g5zP/54WSWn7aMlQZAIaRSMR5XrL3/d9XAgizIX0TLQD4b2QyDTslrQi6A4A/Q/hqEQD477S3PzXadQtzgu4A4A/DBYDQKBYHayVND7oDgD8MFwBCJLFN8pYEXQHAH85cAAiNadMm9UniZajAMMfKBYDQaG/PT3TdwqVBdwDwh+ECQGgUi15cUjLoDgD+MFwACI0jjjjoaWOq5gXdAcAfzlwACA1jjCepN+gOAP6wcgEgNFpbN0x23cIdQXcA8IfhAkBoOE6sLGlr0B0A/GFbBEBoTJ06ebOk84LuAOAPKxcAQiOfz1flcpsODboDgD8MFwBCY/t2b7y13meD7gDgD8MFgBBJ9EmGd+gEhjnOXAAIjWnTJm2TxL1FgGGOlQsAobF69bqRrttxUtAdAPxhuAAQGtGoM0qyM4PuAOAPwwWA0Egk6p6VdEvQHQD84cwFgNBobBzfI2lV0B0A/GHlAkBotLZuHuu6HQuD7gDgD8MFgBApVkv2yKArAPjDcAEgNBKJ/i3G6OKgOwD4w5kLAKHR2Ng4KCkfdAcAf1i5ABAaudzTB7hux/VBdwDwh+ECQGiUSgOOZGuC7gDgD9siAEJj2rSDN0k6LegOAP6wcgEgNFatstFcbv34oDsA+MNwASA06uo2TrI2emXQHQD8YbgAEBqxmO2XbHvQHQD84cwFgNBIpaZ0Sfpa0B0A/GHlAkBotLZurnbdjbxDJzDMMVwACJHiWMnh3iLAMMdwASA0IhHnecneF3QHAH84cwEgNDKZhp2SVgTdAcAfVi4AhEY2mx/luvmZQXcA8IfhAkBoVCreSMmcFHQHAH8YLgCESGKb5C0JugKAP5y5ABAa06ZN6pP0aNAdAPxh5QJAaLS35ye6buHSoDsA+MNwASA0ikUvLikZdAcAf9gWARAa3d0HbR49esP5//TBVfazkmLBFAF7MasWvdP8ZSi+NMMFgNCYPt2UJXX9y4e/JCkRQA6wt3tc0pAMF2yLAAiN1tYNk123cEfQHQD8YbgAEBqOEytL2hp0BwB/2BYBEBpTp07eLOm8oDsA+MPKBYDQaG9vj+dyhYagOwD4w3ABIDSKxaqJ1uqSoDsA+MNwASBEEn2S4R06gWGOMxcAQmPatEnbJHFvEWCYY+UCQGisXr1upOt2cFdUYJhjuAAQGtGoM0qyM4PuAOAPwwWA0CiXvZ2SWRF0BwB/OHMBIDSOPvqQ5yXdF3QHAH9YuQAQGq2tm8e6bsfCoDsA+MNwASBEitWSPTLoCgD+MFwACI1Eon+LMbo46A4A/nDmAkBoNDY2DkrKB90BwB9WLgCERlvbpkmuW7gq6A4A/jBcAAgNzytFJY0LugOAP2yLAAiNadMO3iTptKA7APjDygWA0Fi1ykZzufXjg+4A4A/DBYDQqKvbOMna6JVBdwDwh+ECQGgkEs6gpELQHQD84cwFgNBobGzYIumLQXcA8IeVCwCh0dq6udp1N/IOncAwx3ABIESKYyWHe4sAwxzDBYDQiESc5yXLXVGBYY4zFwBCI5Np2ClpRdAdAPxh5QJAaGSz+VGum58ZdAcAfxguAIRGpeKNlMxJQXcA8IfhAkCIJLZJ3pKgKwD4w5kLAKExbdqkPkmPBt0BwB9WLgCERnt7fqLrFi4NugOAPwwXAEKjWPTikpJBdwDwh+ECQLCbHdAAABBASURBVGh0dx+02Zjy+UF3APCHMxcAQmP6dFOW1BV0BwB/WLkAEBqtrRsmu27hjqA7APjDcAEgNBwnVpa0NegOAP6wLQIgNKZOnbxZ0nlBdwDwh5ULAKHR3t4ez+UKDUF3APCH4QJAaBSLVROt1SVBdwDwh+ECQIgk+iTDO3QCwxxnLgCExrRpk7ZJ4t4iwDDHygWA0Fi9et1I1+3grqjAMMdwASA0olFnlGRnBt0BwB+GCwChUS57OyWzIugOAP5w5gJAaBx99CHPS7ov6A4A/rByASA0Wls3j3XdjoVBdwDwh+ECQIgUqyV7ZNAVAPxhuAAQGmPGOF3GOIuD7gDgD2cuAIRGQ0NDv6Q/Bd0BwB9WLgCERlvbpkmuW7gq6A4A/jBcAAgNzytFJY0LugOAP2yLAAiNadMO3iTptKA7APjDygWA0LDWOrnclpqgOwD4w3ABIDTWrNl4gLX91wfdAcAfhgsAoZFIOIOSCkF3APCHMxcAQqOxsWGLpC8G3QHAH1YuAIRGa+vmatfdyDt0AsMcwwWAECmOlRzuLQIMcwwXAEIjkYj3SFoVdAcAfzhzASA0GhsPfFbSLUF3APCHlQsAoZHN5ke5bn5m0B0A/GG4ABAalYo3UjInBd0BwB+GCwAhktgmeUuCrgDgD2cuAITGtGmT+iQ9GnQHAH9YuQAQGrnc+vGum78o6A4A/jBcAAiNUslUSaYx6A4A/jBcAAiN7u6DNhtTPj/oDgD+cOYCQGhMn27KkrqC7gDgDysXAEKjtXXDZNct3BF0BwB/GC4AhEYsNsKTTG/QHQD8YVsEQGikUgc8LWle0B0A/GHlAkBotLe3x3O5QkPQHQD8YbgAEBrFYtVEa3VJ0B0A/GG4ABAiiT7J8A6dwDDHmQsAoTFt2qRtkri3CDDMsXIBIDTa27tqXbcwPegOAP4wXAAIjWKxe7SkOUF3APBnr9wW+cycOXW9lcQJxpjBWG/swatbri6+kmsAvLrKZW9nNBpdEXQHAH/2uuFi7ty5sb7B2K8lrZPUVaoeWCDp3f/tNQBefUcffcjzku4LugOAP3vdtki8FJ8po45lty1rWrp86QIZM7b5tOaG//YaAK++1tbNY123Y2HQHQD82euGC1k1yprc3x+bJ0zUO+K/vgZAAIrVkj0y6AoA/uyFw4WdYGQLux8aY5+RNPG/vWbt2vwb2to6jmtrWzdOkmrk1QxhNfCaNUblura2juMee2zTwWPGOF3GmLva2jqOa23dPFaSnL3xzykgBKaY0ujdP3t7+mvvdWcurNFaye6/+7EnjY5U/nkP9z+5plJxDnQcTbXW6ZC09YbarcvkaL9v9I5a3lEeMbh0n865A56eP7N3wndnx7dPeteI8gd2WGfdwp5xP7+0elumIVI5en0l/ptFffs+9s1E53tHR9WwcrDqB98b2Kfr1potsxxrq87o2+/6d40YrJsd3/HxXk9bPtk74Yfnjdg5ZWp84N2by07uM8Vxv7uqbtvbxqtyeOtA5KFrBsduuLFm64erHW/Csv59b3e9eN9N1Z3zSlY9s3sm3PE/0ecnnpAofmhHRRsW9k14aFHVjtSU2OBb15fiv1vUv29uSXXn8ftGdPCDpap77ujfp/OW6mdOj0ecmo93T7jumHhP7dwRvaf1eU7n3N5x95wT33bwtBGV4zutkzu/Z9zvrqzd+tYJxku5lREPXdE3asP1ic4P10Y14caBfe94ZDDec3vNM2eXpb45vfvdduqI58afGO//yLNl5T9VnPDTi6t2HP662ODbNpQjv/9Scay7tG7ru0fJm/JQf+JHt5dGPnNLTefH447qmosTbphiB6rOqd55Rm9FWz/ZN+HuTye2Nrw56r13SyXy2AV9Y3/ztcTWtxwQ9dJuqernV/Tvs+76ms4P1jra745yzXcfLNY+f3vtM5/0rBmc1Tvx1pOr+8aeFOlueq5iOhb0jV95cfWOQ18XGfx/hYr+8IW+CW3frN36rtHGe93P+2I/WV4Z/ddb67aeGpO3z/ndE24arYHo/9btnF307Pazeife9emaHfVvdgZP3GpN+7k94x9ZXL31qP0j3tQ/Vap++bW+ff5yQ+2W99cYs//3+mu/v7JUs/O2mi1nORGnctrz42/5YHTnmI8mBj72vGc2ze8df/8XElsb3xD1jtnkRR/9fO+YNVdVbX3n+Jj3+l+Xqu+9ob/u6Vtqt5wSN2bU53tH3Rx3ZBYldp45YO2OM3sm3jk/0X3gW6N9J22tmCfO7Rv/8FcT2988OVqe9mfPPPyV3vFPXFfTdVKdYw+8s7d2xX1ezbO31W45U5LO6Jn47ZOc3tEn1/TM7PbMU2f3jr/vCzVdb3yDY4/dVI62fr445o9XV3cdOy5i37i6Un3/0r66TTfXbJlZ5ZjRi4qjvj3oyX61ZucnBq3dOadn4vc/WdV9wP+L9c3oKjl/Oa9/3C+/WrP9iMlO+cg/l51HvlIc135d7bYT60yl/sfl2rtaijXb7xjZNcereJEzeifedGKsd9SpVT2n9Fr710/2TPzJN2q2T3IcHWttebCh4aANrrtxwHEiR1hbKUjadnPtM5+LR5xxl3WPuvVZjShfWdd5VknOc7O7x31vVuTZ/Y+rLr3/Wes8+amecb/4SnXn1GRERz1Zif/6kr59/7S0uuuEURGb/KlXd/d3equ3Lq/ZMtsxNn56z343vCfRM/K0aO//9Hh6Zl7vhB9dUPXcIelY/3FPlx33ouK4319Rve3tEyOVw/7QH3lwSWnsxptqt30kYSrjr+kbddt6M6J/aaLzk4Oeuuf0TvjO6bHn9zu+qvjBnXLWN3eP+9mXE9vSB0crb3myFP/tJf37rl1a0/meUY4OeuFnr+aZM6JS9em9+113THywdu6IHaf1lNU5rzjhnguqdx6cjgy87M/e7j8DbhyoueORwdqe79R1nj1Y8Xrn9O13+2lVz014T6z/wzs9bWzunfDgy/0Z8EAx8cPvlkduubW287SYUe1ZfROuTzuD1QuqdpxetJGus3rG/mBhbNtBR1VV3tOlyNrzusf+9uuJrW+dFPVSbYMjfnbVwKj1N9R0fqjG0cRbB/f9zi8G4t23VT8zzzOmf3bvxOUfr+kb916n+6M7K6bQ3Df+gUXVOw6bEhl8e74SWf3FvrHZJbVbj9vXeIf8YiD641sHx2z+dk3n/4xwNLL5uQk31kcH4hfW7JzV53nb5vbu17KwakfyqNjgCV1e5E/n9Y799eVVW48+MOZl1parfvH14j5P3lC75QM1xky6vVjzvYfKtc/dXtM5V44pnd49/tamRO+YD0R7PtZtIx1n94xd+YXEs41viJaOyZfMo1/sH7/m2pqt7xzjeK9/eDDxk5sHRv715potp1Y5Zp//LU24qbYyEPlc1c45/Z599hO9E1c0V3dPPjrS977Osnn8/OL4X11e/eyRB0ZKRzxeia36at/oP19fs2VGrWMO+E5P7Z0/tTU7bq/b8glrZV/qZ+/zVV1velPMvuMpG/3j53rGtF5T3TV9bMS+4Xfl6vu+Vax76tu1W04eYcy+i4qjvi1JL/VnwLaK+fM5feNXXV67fdqBpvzmx0vmV1/tH//4t2q63jfSsZPvLo6460flUdtvq91ypjEyp3dPvPm9pnffj9f2nNzj2afn9U6898rRz01wKnqX52lA0oY9+Xex2ZNfLAyaz2h+m4z9ytLly6afferZ+zpxp02epnUVu3aOrRtbf90t1218uWuW3b5se9D9wGtZW9umSY7jXZBOJ88LugXAK7fXLTd29nWulmznwlkLHo3EnSccq6uX3b5s+5iqMWMinrPh310TdDvwWud5paikcUF3AMBLap7TPOkzc+bU+b0GAAAAwDBlrXVyuS0cngaGub1uWwTA8LVmzcYDrO2/PugOAP4wXAAIjUTCGZRUCLoDAAAAAAAAAAAAAAAAAAC8jOYzmucumL3g2pf7/MI5C1MLZ89vfjWbgNeKRYsWRefPnn9g0B0YHni1CIaFC067oEbGLpajq17uGlu2k6znvPXV7AJeK7Y+tXV/x5r7g+7A8LDX3bgMe6eBaPEzsnJMxXxQ0rUL5yw8znqVk60xvZK+t+zWZX/411+z8Iz5p8uY46wxD5la07JkyZKB+bPmpx058421ZWvMdUuXL33s1f9ugOHHVOy5kpnUPGv+fGMiD8urZDxjxjnWWycb3WwjtnHprUu/I0kLz5h/vhOLfefam6/t5GfutYmVCwwP1ntA0g5HzsqF/7NwpPW8y6wxN8nTD4y1d8ydOzf2j5d/6vT5061jPuh5+oqsXme77VWfmTOnzpH5gTX2fivjSva+5jnNkwL6joBhxVpzr6x2qOI8YCu23hqz1Fi7f8TYduvYKbL2uBeuNc5HK+XKvvzMvXYxXGBY8DxnnaSBby7/5rruWPdgJBp9f0Q24hhNsTITR5RGjPrH660TKcsqLUdHGWOWjE2O/XSfVz3DGLVaY7d7Ee9xGbmq6ISAviVgeIlqvYyKS+9YmpckY/Wnpbd968Jrll9feLlfws/caxfDBYadGq9mcqVc/oW1pskzqhhp579es2T5kl8bOWc40lHWer/flt92nqTDrGdHO55zvOM5x8uanKQ/v/rfATD8WZmn/s1nR/7tH/iZe41iuMCwYxz7YSPdt2T5snNV0ZOSDvjXa5pnN8+WrSSWLF/WbCN6l4w93ajyYxn1LF2+9OKly5deLKtaY01fAN8CMFy95N8ZxjE7JLOftOtu05IOliR+5l67ONCJYceYyE88693TPKv5ECMb8YweNvKaJbN69zW2Yn+riO5snr1gnipmtIyuXHLrdaubZy/oap694Hey2sfKrl5621I3yO8FGC7GHTjur9sKW6uaZzVfLk+P/OPnKlWV1U6f2dY8a8Ea49kuSX+xkpYs52cOwDCyaNGi6O6DYYsWLXLOmXXOqJe67pxZ85KLFi36pyH63DPPHT3v9Hn7vxqdwN5k0aJFzty5c6tf7vML/2fhyJf6OD9zAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPCKmaADALy2fex9HztD1jsxEokt+P69398mSTNPbLpcVoMrHmj5UtB9AP573HIdQKBK8dJ9kjnGq5SvlaSZ7/voyVa6UHIeCLoNwCvDygWAwM088aMnWZl7Zc2ZMnaxsfoWqxbA8BUJOgAA2tc9/uShhzQeKKNFktq7ilvPKBQKXtBdAF4ZtkUAhII1tlOSZLT14YcfLgecA8AHVi4ABO7kEz9ylJVzi4y9TNbMbpxy6DPt69vbgu4C8MqwcgEgUDNmzKj25Nwu6Ud33X/3xZL5uhx7VdOMpoag2wC8MgwXAAJV7VUtljS+XKkslKRu2/1lWT2lim5btGgRf0YBAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAwF7r/wMiw+Y7SgPEUAAAAABJRU5ErkJggg==",
"prompt_number": 10,
"svg": [
"\n",
"\n"
],
"text": [
"Plot(...)"
]
}
],
"prompt_number": 10
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Conditioning"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"d1 = @pp enum() do\n",
" local triple = model()\n",
" factor(triple[3] ? 0 : -Inf) # Observe the grass is wet.\n",
" triple[1]\n",
"end "
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 11,
"text": [
"Discrete({false,true}, [0.360656,0.639344])"
]
}
],
"prompt_number": 11
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"bar(d1)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"\n",
"\n"
],
"metadata": {},
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAAhcAAAF6CAYAAACqW3pRAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deYBVdf3/8dfn3IW5syg7iMrcQSxr0jkXSG35mZZlmrSKqIUCJSkM5ZbtRloalZoB7gsulUjZomLagmUZ+h2450KjFcu9g8YygCKz3Jm5957P7w/EFrTvN4/jOYPPx1/dmcPc5xADbz+fz71HAgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAwOvE8uU2nsutGxl2B4BgnLADAGCPuroNY6yNXxl2B4BgGC4AREYq5fRJKoTdASAYE3ZAfzl/yvmpUk3p5IWLFy59qc9fPHNmXVcldaIxpi/RlXjo6qVXF1/rRgAA9kX75MrF7BmzDy7V9F4hY89/qc/PmjUr0e2nHpU0RdJ7S9W9v3htCwG8lJaWTdWet+HIsDsABLNPDheOb66VnKNe7vPJUnKqjNoW3b5oysLFC+fImOHN05obXstGAC+lOFxy5oZdASCYfXK4WHj7osnWtxe87AVWjbIm94/H5ikT9ye+Fm0AXl4s5uyS7P1hdwAIJh52QCisHWWM/rDnoTF2s6TRIRYBkJTJNOyUtCTsDgDBvC6HC2u0WrIH7nnsS0NjFf3Lfy15XtvnJTvF9+2FEyY0/M7zCndKelMiEXtfsbi1Ix4f/idr7d8zmYYPeV4hI+kmSQ+6bvqr2Wz+LGPMXMl+x3UblmSzhe8Yo+N835w5YUL9k57X9pCkwa5bf9SqVRvHOI7/C0me66Y/5Xn5EyTzTWPMnU1N9dfkcoWLrNVpxujipqb0b3O5tsXW2rf4funETGb8jlyu7QlJW1w3ffKaNRuOqFScWyXzK9et/6LntX1CsudZq6symfQPPS//Lckc7zixGUcccfCabLawzBiNdN30pFxu3Uhr48sks8Z162fkcm3vsdbOl/RD101flc3mzzPGfMJa88VMpv5Xnle4RVJTImFPbmxs2OJ5+Scks8N10yd6XlujZG+3VsszmfTncrn8adaaiyRzjevW35nNFr5pjE6QdLbrprOeV7hP0gF/+1v9UfX1fx88aFD5IUlPuW56WjabP9YY811rtSSTSX/H89rmSvYsx9FXjzgi/WAuV7jRWk2QYh9y3YP/7nmFP0m2w3Ub3rd6df4w3zd3Sfq966YvyGbbphhjP2+tXZjJNCz2vMKlkk7yfeecCRPGtnhe4aeSDk4kut7e11ddbYz5tbX6WyaTPiOX2/j/rPWvluyPXbfhW9lsYbYxmimZea5bf7/ntV0n2beWy5WPTpp0yEbPK/xBUp/rpt+dza4/1JjYjyT7R9dt+KznFT4q6UvG6PqmpvTNuVzb16y1k42xc5qaGh73vMKPJaXL5eQxtbW7Yj09Vcslu951G6auWlV4m+NogaSfuW76G56X/7RkzrbWXpbJNPzc8woLJL3NGE1pakrns9m230lSJlP/rlyu0GCtlkpa4brp5mw2/yFjzFcle5PrNtzgeYWvSPqwtfpMJpN+zPPa7pbs+KqqnuM6O/erxON9v5dUcN30KbvPRTjXGmPua2qq/3ouV/iUtTpH0uWum743m81/zxjzTmsrp2cyh6z1vMJvJSVdN/3Olpb1Y+Px2L2S+R/XrT/X89pOluw8a3WrZH/oOLrWWvMGY5zzm5rGPprLFX5grd5orT0+mezuLpVqHpP0tOumP7J6ddtE37c3SFrmuulLstn8dGNMs7VmfiZTv9TzCldKelelomkTJ6af8rz8w5Kpc9302zzv6QOlys+tVTaTSZ+9enXhRN/XZZK53XXrF2Szhc8Zo6mSPue66eWeV7hD0pt7e+MntLUduPMNb2h7XNJm101PzmbzrjHmZmv1UCaT/rLntU2T7GeNsd9tamq4O5crfNtavdsYZ3pT09g/e17hQckOc92GI1tb86NLJXO/pNWum56Zzba91xh7hbX2rkym4XueV7hA0hm+b74wYUL9rz2v7TbJHm5M+aSmpvHtnldosVbtmUz6pNWrnz7c9yu3SfY3rtvw+Wy2cIYxukAy33Pd+rs8r+0Kyb7X9/1PTpgwLud5hfsljW5qqj8ym103zHESD1prWzOZhrM8r3CcpO8Yo7ubmtLfzeXaPmutnSbZL7tuw0PZbOEmY5TxfeeDEyaM3eR5bY9L2um69SesXFl4UyymOyU94rrpizwvP1Uyn7PWLshkGm73vMJlkk40RrOamtKrstn8z40xB5bL29+WSo2qK5UqD0v6i+umP7F69YZjfN+5SjJLXbd+vucVmiVNt1aXZDLpZZ7Xdr1kJ8Vi5Y8cfvj4p3O5wmPWqtt108d73vo3SLEfGmMebWqqPz+XazvFWvsFSde6bvrWbDY/zxhzsrXOuZnM2P/xvMK9ksbuv79959at5WRVVeK3kta6bvr0bLbtHcbYa4zRvU1N6cs9L3+OZD5ljPl6U1P9fblcYZG1Ospae0om01DwvPzvJVNx3fRxa9ZsPKRS8ZdI+pPrpudms/kPG2O+Yq25IZOpv8nzCl+V9CHf19wJE9J/ymYLS41RgzFV75Ika3t+Z63ymUx6Si6XP9pas1DSz103fVk223a2MfbT1tpvZDINP/O8wvclvb1SiZ02ceLB6zyvsFyyMddtOCabzaeNMT+2Vk9kMunZuVzbZGvt1yR7s+s2XB/039bXhTlnznlb84w5j+15PGXKlNi5M88dJ0nNZzW/o3n6nOWSdO4Z5w6ZM31Ofs6Zc4aF1Qpgt5aW9WNfGOQBDGD75JmLlzKsatiwmO+sl6St3VtXSHbr3Olznoglnaccq6sX3bFoR9iNAFLbJX9B2BUA8Io1z2wec/HMmXVhdwAAAADoB7ncupGel/9C2B0AgnndbIsAiL5SyVRJpjHsDgDBMFwAiIyOjnGbjClfGHYHAAAAAADA3ngpKrBvYFsEQGQkEoN8yXSF3QEAAAAAALC31tbWZC5X4A7FwADHtgiAyCgWq0Zbq0vD7gAQDMMFgAhJdUvmibArAAAAAAAA9tba2l77wq2+AQxgbIsAiIxisWOopJlhdwAIhuECQGSUy/5OySwJuwMAAAAAAGBvLS2bhnte29ywOwAEw7YIgAgpVkv2yLArAATDcAEgMoYNc9qNceaH3QEAAAAAALC3Vas2jvG8wlVhdwAIhm0RAJHh+6W4pBFhdwAAAAAAAOzNWuvkcltqwu4AEAzbIgAiY+XKDQdZ23N92B0AgmG4ABAZqZTTJ6kQdgcAAAAAAMDeWlo2VXveBt6hExjg2BYBECHF4ZLDvUWAAY7hAkBkpFLJTknLw+4AAAAAAADYWzabH+x5+alhdwAIhm0RAJFRqfj7SebksDsABMNwASBCUtslf0HYFQAAAAAAAHvL5daN9Lz8F8LuABAM2yIAIqNUMlWSaQy7A0AwDBcAIqOjY9wmY8oXht0BAAAAAACwt5aW9WM9r3Bn2B0AgmFbBEBkJBKDfMl0hd0BAAAAAACwt9bW1mQuV2gIuwNAMGyLAIiMYrFqtLW6NOwOAMEwXACIjETC9ki2NewOAAAAAACAvbW2ttd6XuG4sDsABMO2CIDIKBY7hkqaGXYHgGAYLgBERrns75TMkrA7AAAAAAAA9tba+vRQzyuwLQIMcGyLAIiMYrGvVhIHOoEBjuECQGQMG+a0G+PMD7sDAAAAAABgb6tWbRzjeYWrwu4AEAzbIgAiw/dLcUkjwu4AAAD7CGutk8ttqQm7A0AwrFwAAIBXFcMFgMhYuXLDQdb2XB92B4BgGC4AREYq5fRJKoTdAQAAAAAAsLd8Pl+Vy218S9gdAIJhWwRAZOzY4Y+01v982B0AgmG4ABAZqVSyU9LysDsAAAAAAAD2ls3mB3tefmrYHQCCYVsEQGRUKv5+kjk57A4AwTBcAIiMVKruWUm3ht0BAAAAAACwt1xu3UjPy38h7A4AwbAtAiAySiVTJZnGsDsABMNwASAyOjrGbTKmfGHYHQAAAAAAAHvL5Z45yPPauOU6MMCxLQIAAAAAAAAA/6vW1tZkLldoCLsDQDBsiwCIjGKxarS1ujTsDgDBMFwAiIxEwvZItjXsDgAAAAAAgL21trbXel7huLA7AATDtgiAyCgWO4ZKmhl2B4BgGC4AREa57O+UzJKwOwAAAAAAAPbW2vr0UM8rsC0CDHBsiwCIjGKxr1YSBzqBAY7hAkBkDBvmtBvjzA+7AwAAAAAAYG+trfnRnle4LOwOAMGwLQIgMopFPykpHXYHAADYR1hrnVxuS03YHQCCYeUCAAC8qhguAETGypUbDrK25/qwOwAEw3ABIDIcJ1GWtC3sDgAAAAAAgL3l8/mqXG7jW8LuABAM2yIAImPHDn+ktf7nw+4AEEw87AAA2COVSnaWSpXl//LB+2y1kjIhJQH7rufVo1NNpT++ND+wAKJtue2WlAo7A9jnWH1M7zb39seXZlsEQGSsWLF2P89rOznsDgDBMFwAiIx43Bks2alhdwAIhuECQGSkUnXPSro17A4AwXCgE0BkNDaO7JS0/H+9EECksXIBIDJyuXUjPS//hbA7AATDcAEgMkolUyWZxrA7AATDcAEgMlKpni3G6JKwOwAEw5kLAJHR2NjYJykfdgeAYFi5ABAZudwzB3leG7dcBwY4hgsAAPCqYlsEQGQ0NR30jKRzwu4AEAwrFwAiY/lyG8/l1o0MuwNAMAwXACKjrm7DGGvjV4bdASAYhgsAkZFI2B7JtobdASAYzlwAiIympvHtkr4VdgeAYFi5ABAZra3ttZ5XOC7sDgDBMFwAiIxisWOopJlhdwAIhuECQGTEYs4uyd4fdgeAYDhzASAyMpmGnZKWhN0BIBhWLgBERmvr00M9r8C2CDDAMVwAiIxisa9WEgc6gQGO4QJAZAwb5rQb48wPuwNAMJy5ABAZDQ0NPZL+HHYHgGBYuQAQGa2t+dGeV7gs7A4AwTBcAIiMYtFPSkqH3QEgGIYLAJExceK4Z4yp4pbrwADXb2cupp449RA5lcnWN8dYozFG6jVGq31pTcz6ubuX/eQJSba/nh8AAITjVR8upkye0mDKutIa/8Oy5jlJfzUyT8qq0xp7iLF6ny8nPeUDU9ZK+tLSB5b+4tVuADAwrVy54aB4PPZNSdPCbgHwyr2qw8WpJ556tir2IjnmJmvt7KXLlm55qeuOPfbY+Kia4R+0Mp+bctKUs5cuWzr51ewAMDA5TqIs+dvC7gAQzKs6XFhrn+pU5xEPLnuw9z9d98gjj5Ql3Svp3iknTnnXq9kAYOCaMGHsJkkXhN0BIBgTxpPOmzfPmTdvnh/GcwOIrnw+X7VrV2x8U9PYf7zXxXLbLSkVXhWwj7L6mN5t7u2PL91vrxaZeuKUr089cerR//7xKZOnNLQ+0fpIfz0vgIFrxw5/pLX+58PuABBM/70U1Shpjf+HqSdO+fqxxx4bl3afyTC+VhtpeL89L4ABK5VKdkpaHnYHgGD6dVtk6gdOOc1as0BSXkbbZHWCZL7dYTu+/uCD//lcRhAXz5xZ11VJnWiM6Ut0JR66eunVxX+/Zvb02e935NTtedyX7PvZjTfeWOqvJgCvENsiQP/ox22Rfr23yJIHfnz3qSeeWi1jb5GVZOxl9zyw9JL+fM5Zs2YluvsSj0paK6m9VN07R9J7//06I3OdZO/b87ivr+8+SQwXQIhWrFi7X1VV8hjXrb8/7BYAr1y/DRenTz59uF8pX2NlzzDSLdbYrbLmS6eeNOXNVmp+uZepBpUsJafK2LZFixdNkaTm6XOyzdOaGxbeuTC/55o5Z84ZJqP8wtsWfaY/GgC8MvG4M1iyUyUxXAADWL8NF5VKaZFkjjGyk5cs+/H9knTaSR/7hS/nDiM9JWlIvzyxVaNkcv94bJ4ycX+ipBeHC8U03lgNaZ4+52ErtcvqxkW3L/p9v/QA+D9LpeqeLZW6bw27A0Aw/bctYpxlpWRp9k9/+tMdez5097KfPD558uRMtV81v9+e19pRxugPL2YYu1nS6H9J803cGPtH65jLVbFvMkZLZ82a1XjjjTdu33NNLpc/3Vod6/v22gkTxuWy2cKXjbFjq6p6L9q8+Y3FIUMKi6w12zOZ9Jez2fWHGuNcZIx5oqkpfYvn5U+Q9FHJ3O266eW5XP5ca+VaqysymYZCNlv4jjGm1nXrz3388WeGDRpUutxasy6TSX8nm934VmMqnzLGPNzUlP6J5+WnSnq3MeaGpqb0Ks/Lf0FSg7X6vOumn8/lCtdLzrOuW//FNWs2HlKpVC621mnJZOpvymbb3muMf4rvO0snTKj/tee1zZL8iY7jzD/iiPoNnpf/lqQhrtvw6dWr24b4vv8tSXnXbfhWLleYYK39tLXOrzOZ+qW5XNsp1vrv9f3YTRMmjG3JZtsuNsY/JJGIf7Gx8eBnPa/tesl/3nUbPp/N5tPG6IuSsq7bcH0uV3i3tXaqZH/iuuMezuUKn7LWvrVSiX9n4sSD12WzbZcb4w9rakqf+/jj62qrquLfMca0NTWlL89m864xOtcYLW9qarjb8woflewJkr3Fdcc94XmFCyX7hnJ50JcnTRqzPZcrLLJWRddNX9TSsn5sPO582VqTy2TS12az+WON0enG6KdNTQ2/9Lz8DElHG2Ovamoa99fdd+O0I8vlHc2p1NhBpVLXlZJ52nXT31i9+unDfb/cbK35XSaT/mE2m/+QMTrJ983iCRPSf/K8/PmSDovF7CWHHz5uq+cVFkgquW76As97+kCpfInk/Nl16xesXr3hGN83H7fW/CKTST+QzebPMkZvr1TM9yZOTD/lefmvSxpdU1P+zPPP18Ti8d6rjXE2NTXVf93z2hol/zOS8wfXrb8zl2ubbK1/srXOHZlM/R9zubbPWuu/uVw28yZNSm/2vPw1kuS6DZ9taSkcEI/becY4TzY11V+Tzba9wxj/TGOc+5ua6u/zvLZpkv9OY2ILmprG/tnz2i6R/APL5UHn779/V6WrK/59SVtct+FrK1cW3hSL2fOs1WOZTMPt2WzhA8bYDzqO/cERR4z7vecVmiV7uDGJy5qaDnrG8wpXSUq4bnrumjUbRlUq5lLJ/NV101etWlV4m+PY6dZqWWPjyJ9ns4UDPC9/QyxmFx1++LjVcaN4mRsFAK+69yR73nill7/BmNgDTU1jX9V3y35Vh4spk6eMXXrf0o2SdM8D99z+Utfcd9993ZLm/tOvaVh639L8S137Slij1ZI9cM9jXxoaq/zrEuvC2xf+UdIfX3i4ae70OQ8N6kt8UNKL/8VUKg36VTJZerxcTmyRJGNiix3HH7R69Ru7p0yRv2aNM9+YUkmSBg92nu7ocOZbW+qUJGv1eCzmrO3u7tsuSY5j77U29lAq1bd59+f9hY7jOJLU3X3g86nUxvnlsumRpGSy6qlKpTg/FnN2SlJvb+LXqVTlf6wdtHV3WfxOx/EH/eUvYzsyGWNXr26b39dXLktSbW3l7//cUamUn0gm4+vL5fIOSYrFKj+3Nvbrujp/kySVy/61yWQ8Lkk7doztGDZs43zfd144aFv1V8fpnV8sxp6XpJ6e2PJUSqsSiar23d+DuSsWc6ra2w/atft71LfLZVvZ3VHeXCwm5xtT6dr9e1lZmUzGC5WK/+zur13+heMkfrv//uVNL3Re/0KH3bhxfNdhh/2jw3FSa43pnd/Xl9glSYlE7JFKxfdisVS7JPm+86N43FZJm55/4f+67/q+70vS/vv7W4vF+PxYzO/e/XtZXlVdndy4p8OYygPGJH7X2els3v37U77R2kRi4sSJ5aVL5R92WHG+tbE+Serri61PJv0XO5LJ+KOVir+mqqq4bffvpbk7mTSp3t5tL3xtXeX71u5u3rWtUql9saO7u+JVVyefkfTc7t/L8rJYLPlod/eejsrN1iYS48ePL0kqrVmz8cWOcjmRTyZL88vlvo4XnucPxjityeTujlLJ3pNMOqlUqnPHC38Wr97zZzqV6tpRqdTO7+uzxd2Pi6v7+lLz93T4ft8v4/HkH/v64nv+zN9qjJITJx7QI0lr1mx88c98b2+sUFvrv9gRi5nHJPPUnj/ziYT9caXiLIvFnn/hz4v9vuMYs/vXPvtsMjnixY7q6p4/9/Wl5sdizs6Wlk3DjSmNdRwzv7d3989exaoiKSEAr6pcydnoOM4SyT73an/tV/XVIlNPmnKFld5ujb7d+NbGB//TG2Wd/v6PHVFxYudZ2TcvXbZ0r/fDeKWaz2p+h4z9xsLFi44794xzhzhJZ5V8TWovtu8cXje8/rpbr9vQfFbzJ42xb1mweNH5s2bNSiT7EiuNdT6y4PYF61+tDgD/vd2rTbFvum76H/cW4dUiQP8YKK8WWbJs6RennnjK+yRz9ZNPtN5y6klTnpSx6+RrnXFUlMwbrNVhkt5YkeokXaoanf1qNmzt3rpiVM2IrXOnz3nCSmON1eUL7li045xp54yM+c56SabKr7q7J178VfOM2cvU57iSfspgAYQvlerZUi7X9OsrygD0v355n4t58+Y5T6548v1ydKSVPdyRDrfSgZK2SvaPRvqt322WLn1kaWd/PL8kNc9sHlOt7o5v33prx8tdc94nZh2geHXxe4u/t7O/OgAExMoF0D/6ceUilHuLAMBLyeWeOcjayldct/6cFz/IcAH0j4GyLfLvpkyZklSXpjrSm62xvvH1V9+Yh/vrPS4ADGylUq8Tj8dqwu4AEEy/rVyc/sHTR1Uq5T/Iaqykv77w4bdIes4anbn0gaUP9NdzA9iHsHIB9I+BeFfUSrlyhayeL1cq6XuWLT3inmVLjyiVy6Nl7E+cf3p5HADssXy5jedy60aG3QEgmP67K6rsETJmwb0P3bt5z0d++vBP20vJyhetdNDpJ3yssf+eG8BAVFe3YYy18SvD7gAQTD8OF+Zx+Trw3z8a74iXJZmuQaWN/ffcAAaiRML2SLY17A4AwfTfmYv3f+yIiuM8YKVFjnUeMdZsqTiVQ430NSP9pWz8y/dce/hbDy/8pzfcAvA6xpkLoH8MxFeL+I5zkaSDjHSFNb6s+cckY6V3xKzzyT3Xrn187f6SdvVXC4CBoaVlU3U83vMW1x33RNgtAF65fhsuSpXK52Mm9s3/y7WHHn1opx7srxIAA0dxuBSbK2na/3opgMjqt+HihYOcm//XCyUt/eXS/soAMIDEYs4uyd7/v18JIMr69U20AOC/kck07JS0JOwOAMH046tFAOC/09r69FDPK8wMuwNAMAwXACKjWOyrlXRc2B0AgmG4ABAhqe2SvyDsCgDBcOYCQGRMmjSmWxIvQwUGOFYuAERGa2t+tOcVLgu7A0AwDBcAIqNY9JOS0mF3AAiG4QJAZEycOO4ZY6rOCbsDQDCcuQAQGcYYX1JX2B0AgmHlAkBktLSsH+t5hTvD7gAQDMMFgMhwnERZ0rawOwAEw7YIgMiYMGHsJkkXhN0BIBhWLgBERj6fr8rlNr4l7A4AwTBcAIiMHTv8kdb6nw+7A0AwDBcAIiTVLRneoRMY4DhzASAyJk0as10S9xYBBjhWLgBExooVa/fzvLaTw+4AEAzDBYDIiMedwZKdGnYHgGAYLgBERipV96ykW8PuABAMZy4AREZj48hOScvD7gAQDCsXACKjpWXTcM9rmxt2B4BgGC4AREixWrJHhl0BIBiGCwCRkUr1bDFGl4TdASAYzlwAiIzGxsY+SfmwOwAEw3Dxf/Vb+yEZDQ87A9jnWLXo3SYnSbncMwdZW/mK69afE3YWgFeO4eL/ytElspoQdgawzzH6iqScJJVKvU48HqsJuQhAQAwXACJj0qRDNkqaFnYHgGA40AkgMpYvt/Fcbt3IsDsABMNwASAy6uo2jLE2fmXYHQCCYbgAEBmJhO2RbGvYHQCC4cwFgMhoahrfLulbYXcACIaVCwCR0dKyqdrzNvAOncAAx3ABIEKKwyWHe4sAAxzDBYDIiMWcXZK9P+wOAMFw5gJAZGQyDTslLQm7A0AwrFwAiIxsNj/Y8/JTw+4AEAzDBYDIqFT8/SRzctgdAIJhuAAQIantkr8g7AoAwXDmAkBkTJo0plvSE2F3AAiGlQsAkdHamh/teYXLwu4AEAzDBYDIKBb9pKR02B0AgmG4ABAZHR3jNhlTvjDsDgDBcOYCQGQcd5wpS2oPuwNAMKxcAIiMlpb1Yz2vcGfYHQCCYbgAEBmOkyhL2hZ2B4Bg2BYBEBkTJozdJOmCsDsABMPKBYDIaG1tTeZyhYawOwAEw3ABIDKKxarR1urSsDsABMNwASBCUt2S4R06gQGOMxcAImPSpDHbJXFvEWCAY+UCQGSsWLF2P89r466owADHcAEgMuJxZ7Bkp4bdASAYhgsAkVEu+zslsyTsDgDBcOYCQGQcffShuyTdH3YHgGBYuQAQGS0tm4Z7XtvcsDsABMNwASBCitWSPTLsCgDBMFwAiIxUqmeLMbok7A4AwXDmAkBkNDY29knKh90BIBhWLgBExqpVG8d4XuGqsDsABMNwASAyfL8UlzQi7A4AwbAtAiAyJk06ZKOkaWF3AAiGlQsAkbF8uY3ncutGht0BIBiGCwCRUVe3YYy18SvD7gAQDMMFgMhIpZw+SYWwOwAEw5kLAJHR2NiwRdJXw+4AEAwrFwAio6VlU7XnbeAdOoEBjuECQIQUh0sO9xYBBjiGCwCREYs5uyTLXVGBAY4zFwAiI5Np2ClpSdgdAIJh5QJAZGSz+cGel58adgeAYBguAERGpeLvJ5mTw+4AEAzDBYAISW2X/AVhVwAIhjMXACJj0qQx3ZKeCLsDQDCsXACIjNbW/GjPK1wWdgeAYBguAERGsegnJaXD7gAQDMMFgMjo6Bi3yZjyhWF3AAiGMxcAIuO440xZUnvYHQCCYeUCQGS0tKwf63mFO8PuABAMwwWAyHCcRFnStrA7AATDtgiAyJgwYewmSReE3W2LJdkAAA/5SURBVAEgGFYuAERGa2trMpcrNITdASAYhgsAkVEsVo22VpeG3QEgGIYLABGS6pYM79AJDHCcuQAQGZMmjdkuiXuLAAMcKxcAImPFirX7eV4bd0UFBjiGCwCREY87gyU7NewOAMEwXACIjHLZ3ymZJWF3AAiGMxcAIuPoow/dJen+sDsABMPKBYDIaGnZNNzz2uaG3QEgGIYLABFSrJbskWFXAAiG4QJAZAwb5rQb48wPuwNAMJy5ABAZDQ0NPZL+HHYHgGBYuQAQGatWbRzjeYWrwu4AEAzDBYDI8P1SXNKIsDsABMO2CIDImDTpkI2SpoXdASAYVi4ARIa11snlttSE3QEgGIYLAJGxcuWGg6ztuT7sDgDBMFwAiIxUyumTVAi7A0AwnLkAEBmNjQ1bJH017A4AwbByASAyWlo2VXveBt6hExjgGC4AREhxuORwbxFggGO4ABAZqVSyU9LysDsABMOZCwCR0dh48LOSbg27A0AwrFwAiIxsNj/Y8/JTw+4AEAzDBYDIqFT8/SRzctgdAIJhuAAQIantkr8g7AoAwXDmAkBkTJo0plvSE2F3AAiGlQsAkZHLrRvpefkvhN0BIBiGCwCRUSqZKsk0ht0BIBiGCwCR0dExbpMx5QvD7gAQDGcuAETGcceZsqT2sDsABMPKBYDIaGlZP9bzCneG3QEgGIYLAJGRSAzyJdMVdgeAYNgWARAZTU0HPSPpnLA7AATDygWAyGhtbU3mcoWGsDsABMNwASAyisWq0dbq0rA7AATDcAEgQlLdkuEdOoEBjjMXACJj0qQx2yVxbxFggGPlAkBktLa213pe4biwOwAEw3ABIDKKxY6hkmaG3QEgmH1yW+TimTPruiqpE40xfYmuxENXL726+EquAfDaKpf9nfF4fEnYHQCC2eeGi1mzZiW6+xKPSlorqb1U3TtH0nv/22sAvPaOPvrQXZLuD7sDQDD73LZIspScKqO2RbcvmrJw8cI5MmZ487Tmhv/2GgCvvZaWTcM9r21u2B0AgtnnhgtZNcqa3D8em6dM3J/4X18DIATFaskeGXYFgGD2weHCjjKyhT0PjbGbJY3+b69ZvTp/2KpVbcevWrV2hCTVyK/px2rgdWuYynWrVrUdv2bNxkOGDXPajTH3rFrVdnxLy6bhkuTsi39PAREw3pSG7vnZe7W/9j535sIarZbsgXse+9LQWOVf93D/L9dUKs7BjqMJ1jptkrbdULttkRwd8J2uwYvbyoP6Fu6/dVavr12f7Br1gxnJHWPeM6j8oeess3Zu54hfX1a9PdMQqxy9rpL8w7zuIWu+n9r6/qFxNSzrq/rJD3v3b7+tZst0x9qqs7oPuP49g/rqZiSf+0SXry2f7hr10wsG7Rw/Idn73k1lJ3dxccRjV9Vtf8dIVY5o6Y09/L2+4etvrNn20WrHH7WoZ8gdnp/svql66zklq84ZnaPu/Hh81+gTU8WPPFfR+rndox6eV/Vc0/hE39vXlZKPzesZkltQvfV9Q2I65KFS1b139uy/9dbqzWcmY07NJzpGXXdMsrN21qCuad2+s3VW14h7z0tuP2TSoMr7tlond2HniMeurN329lHGb/Iqgx7+bvfg9dentn60Nq5RN/YOufP3fcnOO2o2n1uWumd2HXD7GYOeH3lSsudjz5aV/0xx1C8vqXruiDck+t6xvhz709eKw72FddveO1j++Id7Uj+7o7Tf5ltrtn4i6aiuuTjqhvG2t+q86p1ndVW07dPdo3782dS2hrfG/fdvqcTWXNQ9/A/fSm1720Fx3/VKVb/+bs/+a6+v2frhWkcH3Fmu+cFDxdpdd9Ru/rRvTd/0rtG3nVbdPfzkWMeU5yumbU73yGWXVD/3ljfE+v5foaLHv9I9atX3a7e9Z6jx3/Dr7sQvFleG/v22um1nJOTvf2HHqJuGqjf+5bqdM4q+3XF21+h7PlvzXP1bnb6TtlnTen7nyN/Pr9521IExf8KfK1W//Vb3/n+9oXbLB2uMOfCHPbU/Wlaq2Xl7zZaznZhTmbZr5K0fju8cdkqq99Rdvtk4u2vkA19JbWs8LO4fs9GPP/GlrmErr6ra9u6RCf+Nj5aq77uhp+6ZW2u3nJ40ZvCXugbfnHRk5qV2frLX2uc+2Tn67tmpjoPfHu8+eVvFPHV+98hHLk/teOvYeHnSX3zzyDe6Rj51XU37yXWOPfjurtol9/s1z95eu+WTknRW5+hbTna6hp5W0zm1wzdPn9s18v6v1LS/6TDHHruxHG/5UnHY/1xd3X7siJh904pK9QMLu+s23lyzZWqVY4bOKw6+pc+Xvbxm56f6rN05s3P0jz5d1XHQ/0t0T24vOX+9oGfEby+v2TFxrFM+8i9l5/ffKI5ova52+0l1plL/83LtPUuLNTvu3K99pl/xY2d1jb7ppETX4DOqOk/vsvbvn+4c/Yvv1OwY4zg61tpyX0PDuPWet6HXcWITra0UJG2/uXbzF5MxZ8Q3Owbf9qwGla+s23p2Sc7zMzpG/HB67NkDj68uffBZ6/ztM50jfvON6q0T0jEd9bdK8tFLu4f8eWF1+4mDYzb9S7/ux3d1VW9bXLNlhmNs8szOA244IdW537R418c7fW0+p2vUzy6qev5QN9Fz/DNlx/tCccSfvlu9/Z2jY5XDH++JPbSgNHzDTbXbP5YylZHf6x58+zozqGdhauun+3x1zOwaddeZiV0HvK+q+OGdctY1d4z41ddT291D4pW3/a2U/OOlPUNWL6zZesJgR+Ne/Nmr2XxWXKo+s+uA645J9tXOGvTctM6ytp5THHXvRdU7D3FjvS/7s7fn74Abe2vu/H1fbedddVvP7av4XTO7D7hjWtXzo05I9Hx0p68NzV2jHnq5vwMeLKZ++oPyfltuq906LWFUe3b3qOtdp696TtVzZxZtrP3szuE/mZvYPu6oqsoJ7YqtvqBj+B+/ndr29jFxv2lV36BfXdU7eN0NNVs/UuNo9G19Q+76TW+y4/bqzef4xvTM6Bq9+BM13SPe73ScsrNiCs3dIx+cV/3c4eNjfe/MV2Irvto9PLugdtvxQ4x/6G964z+/rW/Ypltqtn58kKP9mp8fdWN9vDf5uZqd07t9f/usrgOWzq16Ln1Uou/Edj/25wu6hj96RdW2ow9O+JnV5arffLu4/99uqN3yoRpjxtxRrPnhw+Xa5++o2TpLjimd2THytimprmEfinee2mFjbed2Dl/2ldSzjYfFS8fkS+aJr/aMXHlNzbZ3D3P8Nz7Sl/rFzb37/f3mmi1nVDlm/y+XRt1UW+mNfbFq58we3z77qa7RS5qrO8YeHev+wNayefLC4sjfXVH97JEHx0oTn6wkll/ePfQv19dsmVzrmIPu6qy9+5e25rk76rZ8ylrZl/rZ+1JV+5vfnLDvetrG/+eLncNavlfdftzwmD3ssXL1/dcW656+pXbLaYOMGTKvOPgWSXqpvwO2V8xfzuseufyK2h2TDjbltz5ZMr+7vGfkk9fWtH9gP8eO/XFx0D0/Kw/ecXvtlk8aI3Nmx+ib32+6hnyitvO0Tt8+c07X6PuuHPr8KKei9/i+eiWtfzX/LTav5heLguazmt8hY7+xcPGi484949whTtJZJV+T2ovtO4fXDa+/7tbrNrzcNYvuWLQj7H7g9WzVqo1jHMe/yHXTF4TdAuCV2+eWG7d2b10h2a1zp895IpZ0nnKsrl50x6Idw6qGDYv5zvr/dE3Y7cDrne+X4pJGhN0BAC+peWbzmItnzqwLeg0AAACAAcpa6+RyWzg8DQxw+9y2CICBa+XKDQdZ23N92B0AgmG4ABAZqZTTJ6kQdgcAAAAAAAAAAAAAAAAAAHgZzWc1z5ozY841L/f5uTPnNs2dMbv5tWwCXi/mzZsXnz1j9sFhd2Bg4NUiGBAumnZRjYydL0dXvdw1tmzHWN95+2vZBbxebHt624GONQ+E3YGBYZ+7cRn2Tb3x4sWyckzFfFjSNXNnzj3e+pXTrDFdkn646LZFj//7r5l71uwzZczx1piHTa1ZumDBgt7Z02e7jpzZxtqyNea6hYsXrnntvxtg4DEVe75kxjRPnz3bmNgj8isZ35gRjvXXysY32ZhtXHjbwrskae5Zsy90Eom7rrn5mq38zL0+sXKBgcH6D0p6zpGzbO7H5+5nff+b1pib5Osnxto7Z82alfjnyz9z5uzjrGM+7Pv6hqzeYDvsVRfPnFnnyPzEGvuAlfEke3/zzOYxIX1HwIBirblPVs+p4jxoK7beGrPQWHtgzNhW69jxsvb4F681zimVcmUIP3OvXwwXGBB831krqff7i7+/tiPR0ReLxz8Yk405RuOtzOhBpUGD//l668TKsnLl6ChjzILh6eGf7farJxujFmvsDj/mPykjTxWdGNK3BAwsca2TUXHhnQvzkmSs/rzw9ms/973F1xde7pfwM/f6xXCBAafGrxlbKZd/Y62Z4htVjLTz369ZsHjBo0bOWY50lLX+n7bnt18g6XDr26GO77zP8Z33yZqcpL+89t8BMPBZmaf/w2f3e+F/8DP3OsVwgQHHOPajRrp/weJF56uiv0k66N+vaZ7RPEO2klqweFGzjek9MvZMo8rPZdS5cPHCSxYuXniJrGqNNd0hfAvAQPWS/2YYxzwnmQOk3XeblnSIJPEz9/rFgU4MOMbEfuFb/97m6c2HGtmYb/SIkd8smRV7rrEV+0fFdHfzjDnnqGKGyujKBbddt6J5xpz25hlzHpPV/lZ2xcLbF3phfi/AQDHi4BF/317YVtU8vfkK+fr9P3+uUlVZ4XSb7c3T56w0vm2X9FcracFifuYADCDz5s2L7zkYNm/ePOe86ecNfqnrzpt+TnrevHn/MkSf/8nzh55z5jkHvhadwL5k3rx5zqxZs6pf7vNzPz53v5f6OD9zAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPCKmbADALy+nfqBU8+S9U+KxRJzfnTfj7ZL0tSTplwhq74lDy79Wth9AP573HIdQKhKydL9kjnGr5SvkaSpHzjlNCt9TnIeDLsNwCvDygWA0E096ZSTrcx9suaTMna+sbqWVQtg4IqFHQAArWuf/NtbDm08WEbzJLW2F7edVSgU/LC7ALwybIsAiARr7FZJktG2Rx55pBxyDoAAWLkAELrTTvrYUVbOrTL2m7JmRuP4t2xuXde6KuwuAK8MKxcAQjV58uRqX84dkn52zwM/vkQy35Zjr5oyeUpD2G0AXhmGCwChqvar5ksaWa5U5kpSh+34uqyeVkW3z5s3j7+jAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAGCf9f8BSITmztsYLdMAAAAASUVORK5CYII=",
"prompt_number": 12,
"svg": [
"\n",
"\n"
],
"text": [
"Plot(...)"
]
}
],
"prompt_number": 12
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Explaining Away"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"d2 = @pp enum() do\n",
" local triple = model()\n",
" factor(triple[3] && triple[2] ? 0 : -Inf) # Also observe that the sprinkler is on.\n",
" triple[1]\n",
"end"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 13,
"text": [
"Discrete({false,true}, [0.688525,0.311475])"
]
}
],
"prompt_number": 13
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"bar(d2)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"\n",
"\n"
],
"metadata": {},
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAAhcAAAF6CAYAAACqW3pRAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deWBcVd3/8c+5szSTJm3pki5AM4GCS4DMlAq4PEofEQWpG4QCCpQCFdpUZVP0p1hRxCKL2JZdaAGVEsQFKCJqq4hWTDN3WgNKl5kU7JK2tDRpJsnM3PP7oxQXQB+5DfemvF9/Mclt5h1K2i/nnJkrAQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA4E1i6VIbzWbX1ATdAcAfJ+gAANijunrdOGuj1wXdAcAfhgsAoZFIOH2S8kF3APDHBB3QXy5qvChRHFw8af7C+c2v9vnPT59evaucOMEY0xfbFXvshuYbCm90IwAA+6J9cuVi5jkzDywO7r1axl70ap+fMWNGrNtLPCGpUdIHipW9P3tjCwG8mpaWDZWuu+6ooDsA+LNPDheOZ26SnKNf6/PxYnyqjNoXLFrQOH/h/FkyZmTTmU11b2QjgFdTGCk5s4OuAODPPjlczF+0YIr17MWveYFVvazJ/v2xecZEvSPfiDYAry0ScXZK9uGgOwD4Ew06IBDWjjZGv9vz0Bi7UdKYAIsASEqn63ZIWhx0BwB/3pTDhTVaKdn99zz2pOGRsv7p/5Zct/0Lkm30PHvJxIl1v3Hd/D2S3haLRY4vFDZ3RqMj/2Ct/Vs6XfdR182nJd0u6dFUKvmVTCZ3tjFmtmS/nUrVLc5k8t82RpM9z5w1cWLt067b/pikYalU7dGtrevHOY73M0luKpU8z3VzH5TMVcaYexoaam/MZvOXWqvTjNHnGxqSv85m2xdaaw/zvOIJ6fSEbdls+1OSNqVSyZNWrVp3RLns3CmZx1Op2i+6bvunJPs5a3V9Op38gevmviWZ4xwncs4RRxy4KpPJLzFGNalUclI2u6bG2ugSyaxKpWrPyWbb32+tnSvpB6lU8vpMJvc5Y8ynrDVfTKdrH3fd/PckNcRi9qT6+rpNrpt7SjLbUqnkCa7bXi/ZRdZqaTqdvCybzZ1mrblUMjemUrX3ZDL5q4zRByWdn0olM66bf0jS2GefrT26tvZvwwYNKj0m6ZlUKnlmJpM71hhzrbVanE4nv+267bMle7bj6CtHHJF8NJvN32atJkqRj6ZSB/7NdfN/kGxnKlV3/MqVubd6nrlX0m9TqeTFmUx7ozH2C9ba+el03ULXzV8p6UTPcy6YOHF8i+vmfyzpwFhs17v6+iorjTG/tFbPptPJM7LZ9f9jrXeDZB9Ipeq+lcnkZxqj6ZKZk0rVPuy67TdL9h2lUvkTkyYdvN5187+T1JdKJf83k1l7iDGRH0r2yVSq7rOum/+EpC8Zo1saGpJ3ZLPtX7XWTjHGzmpoqPuj6+YfkJQsleLvraraGenpqVgq2bWpVN3U1tb8Ox1H8yT9JJVKfsN1c5+WzPnW2q+n03U/dd38PEnvNEaNDQ3JXCbT/htJSqdr35fN5uusVbOk5alUsimTyX3UGPMVyd6eStXd6rr5L0v6mLX6TDqd/L3rtt8n2QkVFT2Tu7qGlKPRvt9KyqdSyVN2n4twbjLGPNTQUPu1bDZ/nrW6QNI3U6nkg5lM7jvGmPdYWz49nT54tevmfy0pnkol39PSsnZ8NBp5UDJ/SqVqL3Td9pMkO8da3SnZHziObrLWHGqMc1FDw/gnstn8963VW6y1x8Xj3d3F4uDfS3oulUp+fOXK9iM9z94qaUkqlbwik8lNM8Y0WWvmptO1za6bv07S+8plnXnkkclnXDf3C8lUp1LJd7ruc/tL5Z9aq0w6nTx/5cr8CZ6nr0tmUSpVOy+TyV9mjKZKuiyVSi513fzdkt7e2xv9YHv7/jsOPbT9j5I2plLJKZlMLmWMucNaPZZOJ/+f67afKdnPGmOvbWiouy+bzV9jrf7XGGdaQ8P4P7tu/lHJjkil6o5qa8uNKRbNw5JWplLJ6ZlM+weMsVdba+9Np+u+47r5iyWd4Xnm8okTa3/puu13SfZwY0onNjRM6HDdfIu16kinkyeuXPnc4Z5Xvkuyv0ql6r6QyeTPMEYXS+Y7qVTtva7bfrVkP+B53rkTJx6Udd38w5LGNDTUHpXJrBnhOLFHrbVt6XTd2a6bnyzp28bovoaG5LXZbPtnrbVnSvb/pVJ1j2Uy+duNUdrznI9MnDh+g+u2/1HSjlSq9oMrVuTfFonoHknLUqnkpa6bmyqZy6y189LpukWum/+6pBOM0YyGhmRrJpP7qTFm/1Jp6zsTidHVxWL5F5L+kkolP7Vy5br3ep5zvWSaU6naua6bb5I0zVpdkU4nl7hu+y2SnRSJlD5++OETnstm87+3Vt2pVPI41117qBT5gTHmiYaG2ouy2fZTrLWXS7oplUremcnk5hhjTrLWuTCdHv8n180/KGn80KH2PZs3l+IVFbFfS1qdSiVPz2Ta322MvdEYPdjQkPym6+YukMx5xpivNTTUPpTN5hdYq6Ottaek03V51839VjLlVCo5edWq9QeXy95iSX9IpZKzM5ncx4wxX7bW3JpO197uuvmvSPqo52n2xInJP2Qy+WZjVGdMxfskydqe31irXDqdbMxmc8dYa+ZL+mkqlfx6JtN+vjH209bab6TTdT9x3fx3Jb2rXI6cduSRB65x3fxSyUZSqbr3ZjK5pDHmAWv1VDqdnJnNtk+x1n5VsnekUnW3+P279U1h1lmz3tl0zqzf73nc2NgYuXD6hQdJUtPZTe9umjZrqSRdeMaF+82aNis366xZI4JqBbBbS8va8S8N8gAGsH3yzMWrGVExYkTEc9ZK0ubuzcslu3n2tFlPReLOM47VDQvuXrAt6EYAia2SNy/oCgB43ZqmN437/PTp1UF3AAAAAOgH2eyaGtfNXR50BwB/3jTbIgDCr1g0FZKpD7oDgD8MFwBCo7PzoA3GlC4JugMAAAAAAOCVeCkqsG9gWwRAaMRigzzJ7Aq6AwAAAAAA4JXa2tri2WyeOxQDAxzbIgBCo1CoGGOtrgy6A4A/DBcAQiTRLZmngq4AAAAAAAB4pba2jqqXbvUNYABjWwRAaBQKncMlTQ+6A4A/DBcAQqNU8nZIZnHQHQAAAAAAAK/U0rJhpOu2zw66A4A/bIsACJFCpWSPCroCgD8MFwBCY8QIp8MYZ27QHQAAAAAAAK/U2rp+nOvmrw+6A4A/bIsACA3PK0YljQq6AwAAAAAA4JWstU42u2lw0B0A/GFbBEBorFix7gBre24JugOAPwwXAEIjkXD6JOWD7gAAAAAAAHillpYNla67jnfoBAY4tkUAhEhhpORwbxFggGO4ABAaiUS8S9LSoDsAAAAAAABeKZPJDXPd3NSgOwD4w7YIgNAol70hkjkp6A4A/jBcAAiRxFbJmxd0BQAAAAAAwCtls2tqXDd3edAdAPxhWwRAaBSLpkIy9UF3APCH4QJAaHR2HrTBmNIlQXcAAAAAAAC8UkvL2vGum78n6A4A/rAtAiA0YrFBnmR2Bd0BAAAAAADwSm1tbfFsNl8XdAcAf9gWARAahULFGGt1ZdAdAPxhuAAQGrGY7ZFsW9AdAAAAAAAAr9TW1lHluvnJQXcA8IdtEQChUSh0Dpc0PegOAP4wXAAIjVLJ2yGZxUF3AAAAAAAAvFJb23PDXTfPtggwwLEtAiA0CoW+Kkkc6AQGOIYLAKExYoTTYYwzN+gOAAAAAACAV2ptXT/OdfPXB90BwB+2RQCEhucVo5JGBd0BAAD2EdZaJ5vdNDjoDgD+sHIBAAD2KoYLAKGxYsW6A6ztuSXoDgD+MFwACI1EwumTlA+6AwAAAAAA4JVyuVxFNrv+sKA7APjDtgiA0Ni2zaux1vtC0B0A/GG4ABAaiUS8S9LSoDsAAAAAAABeKZPJDXPd3NSgOwD4w7YIgNAol70hkjkp6A4A/jBcAAiNRKL6BUl3Bt0BAAAAAADwStnsmhrXzV0edAcAf9gWARAaxaKpkEx90B0A/GG4ABAanZ0HbTCmdEnQHQAAAAAAAK+UzT5/gOu2c8t1YIBjWwQAAAAAAAAA/qO2trZ4NpuvC7oDgD9siwAIjUKhYoy1ujLoDgD+MFwACI1YzPZIti3oDgAAAAAAgFdqa+uoct385KA7APjDtgiA0CgUOodLmh50BwB/GC4AhEap5O2QzOKgOwAAAACEiAk6YMBYZm+V1VuCzgD2QXdosrlXktranhteLJY/lkol7ww6CsDrFw06YACZJGli0BHAPujxPf9QKPRVRaORyZIYLoABjDMXAEJjxAinwxhnbtAdAPxh5QJAaNTV1fVI+nPQHQD8YeUCQGi0teXGuG7+60F3APCH4QJAaBQKXlxSMugOAP4wXAAIjSOPPOh5YyouCLoDgD8MFwAAYK9iuAAQGitWrDvA2p5bgu4A4A/DBYDQcJxYSdKWoDsA+MNLUQGExsSJ4zdIujjoDgD+sHIBIDRyuVxFNrv+sKA7APjDcAEgNLZt82qs9b4QdAcAfxguAIRGIhHvkrQ06A4A/nDmAkBo1Ncf+IK4aRkw4LFyASA0li9fPcR1208KugOAPwwXAEIjGnWGSXZq0B0A/GG4ABAaiUQ12yLAPoAzFwBCo76+hgOdwD6AlQsAoZHNrqlx3dzlQXcA8IfhAkBoFIumQjL1QXcA8IfhAkBoJBI9m4zRFUF3APCHMxcAQqO+vr5PUi7oDgD+sHIBIDSy2ecPcN12brkODHAMFwAAYK9iWwRAaDQ0HPC8pAuC7gDgDysXAEJj6VIbzWbX1ATdAcAfhgsAoVFdvW6ctdHrgu4A4A/DBYDQiMVsj2Tbgu4A4A9nLgCERkPDhA5J3wq6A4A/rFwACI22to4q181PDroDgD8MFwBCo1DoHC5petAdAPxhuAAQGpGIs1OyDwfdAcAfzlwACI10um6HpMVBdwDwh5ULAKHR1vbccNfNsy0CDHAMFwBCo1Doq5LEgU5ggGO4ABAaI0Y4HcY4c4PuAOAPZy4AhEZdXV2PpD8H3QHAH1YuAIRGW1tujOvmvx50BwB/GC4AhEah4MUlJYPuAOAPwwWA0DjyyIOeN6aCW64DA1y/nbmYesLUg+WUp1jPvNcajTNSrzFa6UmrItbL3rfkR09Jsv31/AAAIBh7fbhonNJYZ0q6zhrvY7Jmu6S/GpmnZdVljT3YWB3vyUk2frhxtaQvNT/S/LO93QBgYFqxYt0B0WjkKklnBt0C4PXbq8PFqSecer7K9lI55nZr7czmJc2bXu26Y489Njp68MiPWJnLGk9sPL95SfOUvdkBYGBynFhJ8rYE3QHAn706XFhrn+lS1xGPLnm0999dt2zZspKkByU92HhC4/v2ZgOAgWvixPEbJF0cdAcAf/bqgc7mnzf/7tFH//1gIUlz5sx5+XmbH23+zd5sADBw5XK5imx2/WFBdwDwp99eLTL1hMavTT1h6jH/+vHGKY11bU+1Leuv5wUwcG3b5tVY630h6A4A/vTfS1GN4tZ4v5t6QuPXjj322Ki0+0yG8bTSSCP77XkBDFiJRLxL0tKgOwD4Y/rzi0/98CmnWWvmScrJaIusPiiZazpt59f+L9snr9fnp0+v3lVOnGCM6Yvtij12Q/MNhX+9Zua0mR9y5FTvedwX7/vJbbfdVnzNL7rMrpDVxH5KBt7MvqzJ5qqgIwDsPf16b5HFjzxw36knnFopY78nK8nYr9//SPMV/fmcM2bMiHX3xZ6QtFpSR7Gyd5akD/zrdUbmZsk+tOdxX1/fQ5Jee7gA0O+WL189pKIi/t5UqvbhoFsAvH79NlycPuX0kV65dKOVPcNI37PGbpY1Xzr1xMa3W6nptV6m6le8GJ8qY9sXLFzQKElN02Zlms5sqpt/z/zcnmtmnTVrhIxy8+9a8Jn+aADw+kSjzjDJTpXEcAEMYP02XJTLxQWSea+RnbJ4yQMPS9JpJ578M0/O3UZ6RtJ+/fLEVvWSyf79sXnGRL0jJb08XCiiCcZqv6Zps35hpQ5Z3bZg0YLf9ksPgP+zRKL6hWKx+86gOwD403/bIsZZUowXZ/74xz/etudD9y350R+nTJmSrvQq5vbb81o72hj97uUMYzdKGvNPaZ6JGmOftI75psr2bcaoecaMGfW33Xbb1j3XZLO5063VsZ5nb5o48aBsjVMe01GO9Fs28Gb1Vqd3/H1u7lZrtaS+vuanmUx+rOvmbo1E7ILDDz9opevmvyzZA2OxwZcUCut7o9ER8yXTkUolv+K6aw+VnEskLU+l6u7KZnMfslYft1Y/TKfrlmUy+ZnG2IZSybtq0qSD17tu/lpjlGhoSM5qadkwMhrtvcoYs7qhIXmt6647SjLnSuaxVCr5YDabO81aTbZWN6fTdW42m/uitUr29JQuO/roCV3ZbP5ma51t6XTtl1aseG5CJFK6zBjzp4aG5B2uu+54yZxsjFnc0JD8tevmPi1pojHmWw0NyZzr5uZKztBUqvaCtrbnhheLpautddam07XXtLaun+Q45fONcR5vaKh9IJNpbzTGO85xnNuOOKJ2heu2f0HyDnIc5/Ijjqjd7rq5WyVtT6XqLs9m83XW2ssltaZSdbe2trYf5zheo7XOA+l07eOZTPv5xniTIpHINYcfPn6t67ZfLXnDGxqSF7hufqgxmmuM8g0NdVe7bj4t2Qsk/TqVqluczeZPttYeb23kjnR6/J8ymfxlxtgJvb2xLx199AHbXLf9ZmttVzqdvKy1dV2t45gvGSO3oaHuZtfNT5bsaZIeTKXqHstm8+daa4+SvOtSqYOfzWTyVxljR27fnpw1duxfEz09g6611qxPp5NXtbaua3AcM9MYLWtoqPuh6+Y/LtkPGaO7Ghrqlrtu/mLJvsWY8lcaGiZ0uG5+vqTeVCp5yapVaw4slyNflsyqVCo5v7U19z7H0RnWmp+m08klmUxumjF6p+PohiOOqPuL6+avlOzoWKx7dqk0ImZt4XrJ+VsqVXtlNrv+MGvLs601v02nk9/PZtd/xNryh601i9Lp5O8zmdznjNHbYjF9tb6+bpPr5r9rjCk3NNRe1Nq6fpzjlL9qjNPW0FD7Xddtf4/knSk5D6VStQ9ns+1nWeu92/OcGydOrH06k8nNMUZjhw7VZyXpxRd1o7XamE7XzWltbX+743ifNcZ5sqGh9m7XbT9J8qZIzj2pVO3vXLd9tuQdJkWvTKUO/Fs2236DtTaSSiU/09aWG1Ms6mvW6pl0uu47mUz+XcbYs42JPNLQMH6vvlv2Xh0uGqc0jm9+qHm9JN3/yP2LXu2ahx56qFvS7H/4NXXNDzXnXu3a18MarZTs/nsee9LwSPmfl1jnL5r/pKQnX3q4Yfa0WY8N6ot9RNLL/8dULA56PB4v/rFUim2SpO1eZJukcXurE8Buz9vYZsdxFkYizo6Wlg0jjSmOdxwzt7c39tLWaeQux/EGtbWNKjQ2jvJWrVo/15hiUZK6u+PPVVV5c60tdkmS52l5JOI8293dt1WS4nH7YLns/Hzo0NLm3Z/35kUizkuvktvwouOMmlsqmR5JisWqni6XC3MjEWeHJBWLg34ZjxefsnbQZkmyNnq343iD1q+fsOvoo2Udx5nb11cqSdLw4aXnOzudlzusNU9FIs6avr7SNkmKROxPrI08nkj0bdjdYW+KRhWRpI6OA3aOGLH+5Y5Bgyr+Ui4X5hYKkRclqa8v8utEQiv2dHieuScadSq2bRvfKUn/2FFZWdxQKMTnGlPetfva0p+i0ei6PR2OU/ypMbFfVVWVN0hSqVS6OR6PRo0x9v77bedb37p+ruc5vZJkTMWzxvS+3BGNRpaWy14mEqno2P09Ot+PRGxFd/f+L+6+3l7jeZ4nSdXV5U3/2NHTU1xRWRlvL5e9F3b/uy895DixpdXV2rD7a5dusTYWO/ZYlZub39L9jx2eV7E6Gi3O7euL7dz9+xT5TbnsZSORRMfu70E/jMedRF/ftu27O3Rduby7w3F6N1tbNTcS8bp3/7ssZSor48/t6bC29EgkEv9tV5ezcffvU+l2a2Oxt7/97UVJpVWr1s+1NtK3+7+H6Lp43L7cYYx9whjnz/F4YYsklctmcTxuEoXClm27v7a93trdt88aNKhza7n8946KikK2ry+xQdL23b+279FoNP677u7dHY4Tu8OYcjyZHN8nSf/Y0dPj5KqqNLdU6nvp919PSs7TezpiMa+5XHYeiUR2bnnpv/kb9/ysFQrbtsXjo+b29dmCJCUSPav6+hJzJbtde9lefbXI1BMbr7bSu6zRNfXvqH90zpw53mtde/qHTj6i7EQ+Z2Xf3ryk+RXvh/F6NZ3d9G4Z+435CxdMvvCMC/dz4k6rPE3qKHTsGFk9svbmO29e13R207nG2MPmLVxw0YwZM2LxvtgKY52Pz1s0b+1rfmFeLQL0l5dfLdLSsnZ8NBq5KpVKcm8RYADbqysXi5c0f3HqCaccL5kbnn6q7Xunntj4tIxdI09rjKOCZA61Vm+V9JayVC3pSg3W+XuzYXP35uWjB4/aPHvarKesNN5YfXPe3Qu2XXDmBTURz1kryVR4Fff1RAuPN50zc4n6nJSkH//bwQLAGyKR6NlUKg3u11eUAeh//fI+F3PmzHGeXv70h+ToKCt7uCMdbqX9JW2W7JNG+rXXbZqblzV39cfzS1LT9KZxleruvObOOztf65rPfWrGWEUrC99Z+J0d//ELsnIB9Bfe5wLYx/Trm2jtUxgugP7y8nCRzT5/gLXlL6dStRcEHQXg9evXN9FqbGyMa5emOtLbrbGe8fRXz5hf9Nd7XAAY2IrFXicajQwOugOAP/22cnH6R04fXS6Xfier8ZL++tKHD5O03Rqd1fxI8yP99dz9gpULoL+wLQLsY/rtxmXlUvlqWb1YKpeT9y9pPuL+Jc1HFEulMTL2R47VDf31vAAGrqVLbTSbXVMTdAcAf/rvrqiyR8iYeQ8+9uDGPR/58S9+3FGMl79opQNO/+DJ9f333AAGourqdeOsjV4XdAcAf/pxuDB/lKf9//Wj0c5oSZLZNai4vv+eG8BAFIvZHsm2Bd0BwJ9+O9AZ8cq3lh3nkcYTGz3HOsuMNZvKTvkQI33VSN+PlWOjTv7wyaMk6fB3HJ7/d2+4BeDNoaFhQoekbwXdAcCffhsuPMe5VNIBRrraGk/W/P30qJXeHbHOuXuuXf3H1UMl7eyvFgADQ0vLhspotOewVOqgp4JuAfD69dtwUSyXvxAxkf/TCfBDjjmkS4/2VwmAgaMwUorMlsTbfwMDWL8NFy8d5Nz4Hy+U1Pzz5v7KADCARCLOTsk+/J+vBBBm/fomWgDw30in63ZIWhx0BwB/+vHVIgDw32lre2646+anB90BwB+GCwChUSj0VUmaHHQHAH8YLgCESGKr5M0LugKAP5y5ABAakyaN65bEy1CBAY6VCwCh0daWG+O6+a8H3QHAH4YLAKFRKHhxScmgOwD4w3ABIDSOPPKg542puCDoDgD+cOYCQGgYYzxJu4LuAOAPKxcAQqOlZe14183fE3QHAH8YLgCEhuPESpK2BN0BwB+2RQCExsSJ4zdIujjoDgD+sHIBIDRyuVxFNrv+sKA7APjDcAEgNLZt82qs9b4QdAcAfxguAIRIolsyvEMnMMBx5gJAaEyaNG6rJO4tAgxwrFwACI3ly1cPcd32k4LuAOAPwwWA0IhGnWGSnRp0BwB/GC4AhEYiUf2CpDuD7gDgD2cuAIRGfX1Nl6SlQXcA8IeVCwCh0dKyYaTrts8OugOAPwwXAEKkUCnZo4KuAOAPwwWA0EgkejYZoyuC7gDgD2cuAIRGfX19n6Rc0B0A/GHlAkBoZLPPH+C67bcE3QHAH4YLAKFRLPY6kh0cdAcAf9gWARAakyYdvF7SmUF3APCHlQsAobF0qY1ms2tqgu4A4A/DBYDQqK5eN87a6HVBdwDwh+ECQGjEYrZHsm1BdwDwhzMXAEKjoWFCh6RvBd0BwB9WLgCERkvLhkrXXcc7dAIDHMMFgBApjJQc7i0CDHAMFwBCIxJxdkr24aA7APjDmQsAoZFO1+2QtDjoDgD+sHIBIDQymdww181NDboDgD8MFwBCo1z2hkjmpKA7APjDcAEgRBJbJW9e0BUA/OHMBYDQmDRpXLekp/7pg0vtfZIGBRIE7NvmarJZ3h9fmOECQGi0teXGFItmViqV/Mo/fPgjkhJBNQH7LKt7+utLsy0CIDQKBS8uKRl0BwB/GC4AhEZn50EbjCldEnQHAH/YFgEQGpMnm5KkjqA7APjDygWA0GhpWTvedfP9tg8M4I3BcAEgNBwnVpK0JegOAP6wLQIgNCZOHL9B0sVBdwDwh5ULAKHR1tYWz2bzdUF3APCH4QJAaBQKFWOs1ZVBdwDwh+ECQIgkuiXz1H++DkCYceYCQGhMmjRuqyTuLQIMcKxcAAiN5ctXD3Hddu6KCgxwDBcAQiMadYZJdmrQHQD8YbgAEBqlkrdDMouD7gDgD2cuAITGMcccslPSw0F3APCHlQsAodHSsmGk67bPDroDgD8MFwBCpFAp2aOCrgDgD8MFgNBIJHo2GaMrgu4A4A9nLgCERn19fZ+kXNAdAPxh5QJAaLS2rh/nuvnrg+4A4A/DBYDQ8LxiVNKooDsA+MO2CIDQmDTp4PWSzgy6A4A/rFwACI2lS200m11TE3QHAH8YLgCERnX1unHWRq8LugOAPwwXAEIjkXD6JOWD7gDgD2cuAIRGfX3dJklfCboDgD+sXAAIjZaWDZWuu4536AQGOIYLACFSGCk53FsEGOAYLgCERiTi7JQsd0UFBjjOXAAIjXS6boekxUF3APCHlQsAoZHJ5Ia5bm5q0B0A/GG4ABAa5bI3RDInBd0BwB+GCwAhktgqefOCrgDgD2cuAITGpEnjuiU9FXQHAH9YuQAQGm1tuTGum/960B0A/GG4ABAahYIXl5QMugOAPwwXAEKjswNVOwgAABBFSURBVPOgDcaULgm6A4A/nLkAEBqTJ5uSpI6gOwD4w8oFgNBoaVk73nXz9wTdAcAfhgsAoeE4sZKkLUF3APCHbREAoTFx4vgNki4OugOAP6xcAAiNtra2eDabrwu6A4A/DBcAQqNQqBhjra4MugOAPwwXAEIk0S0Z3qETGOA4cwEgNCZNGrdVEvcWAQY4Vi4AhMby5auHuG47d0UFBjiGCwChEY06wyQ7NegOAP4wXAAIjVLJ2yGZxUF3APCHMxcAQuOYYw7ZKenhoDsA+MPKBYDQaGnZMNJ122cH3QHAH4YLACFSqJTsUUFXAPCH4QJAaIwY4XQY48wNugOAP5y5ABAadXV1PZL+HHQHAH9YuQAQGq2t68e5bv76oDsA+MNwASA0PK8YlTQq6A4A/rAtAiA0Jk06eL2kM4PuAOAPKxcAQsNa62SzmwYH3QHAH4YLAKGxYsW6A6ztuSXoDgD+MFwACI1EwumTlA+6A4A/nLkAEBr19XWbJH0l6A4A/rByASA0Wlo2VLruOt6hExjgGC4AhEhhpORwbxFggGO4ABAaiUS8S9LSoDsA+MOZCwChUV9/4AuS7gy6A4A/rFwACI1MJjfMdXNTg+4A4A/DBYDQKJe9IZI5KegOAP4wXAAIkcRWyZsXdAUAfzhzASA0Jk0a1y3pqaA7APjDygWA0Mhm19S4bu7yoDsA+MNwASA0ikVTIZn6oDsA+MNwASA0OjsP2mBM6ZKgOwD4w5kLAKExebIpSeoIugOAP6xcAAiNlpa14103f0/QHQD8YbgAEBqx2CBPMruC7gDgD9siAEKjoeGA5yVdEHQHAH9YuQAQGm1tbfFsNl8XdAcAfxguAIRGoVAxxlpdGXQHAH8YLgCESKJbMrxDJzDAceYCQGhMmjRuqyTuLQIMcKxcAAiNtraOKtfNTw66A4A/DBcAQqNQ6BwuaXrQHQD82Se3RT4/fXr1rnLiBGNMX2xX7LEbmm8ovJ5rALyxSiVvRzQaXRx0BwB/9rnhYsaMGbHuvtgTklZL6ihW9s6S9IH/9hoAb7xjjjlkp6SHg+4A4M8+ty0SL8anyqh9waIFjfMXzp8lY0Y2ndlU999eA+CN19KyYaTrts8OugOAP/vccCGrelmT/ftj84yJekf+19cACEChUrJHBV0BwJ99cLiwo41sfs9DY+xGSWP+22tWrsy9tbW1/bjW1tWjJGmwvMH9WA28aY1Qqbq1tf24VavWHzxihNNhjLm/tbX9uJaWDSMlydkX/5wCQmCCKQ7f87O3t7/2PnfmwhqtlOz+ex570vBI+Z/3cP8v15TLzoGOo4nWOu2Sttya2LRAETP22l3D78qXBvXNH7rx00Vp5zmdY+89d9C2cZPjfR/b7kVXz9416vFvDN6aTjrFd661FU98tWu/VTdVb/7QEHkH/bw49IF7eyo7FlZvnhY1qvjUztG3HFfRVz0ttu3MbjkbZ3SO/vGlFTsmpGKF4zd6UfeyXaN+f33V1vfUmOIRK/pij93QO3LtHdWbT66QN/rmnuGLWsuDum8fvPFCL+J0nrVj9D3TK18c87+R7k/ssM6apq7Rv/ha5faGgyM9715brnjyq937ZedXbT5+mPEmPF6q/NGiwtDNi4ZsPitivapPdY696dhBnVXnxbvO6pGz+bzO0T+6aNDWg4+MFz+42YtmL9k16snrBm9592in1OAWY7+4tmfkmturN38iIW/MHX0j7l7WG++6p2rzhWWpe1rX6EWfquiu+VDsxVN2ylk3s3P0z+ckth8xIdrznlw59vuvdI90bxqy5QNDbOmQX/YO/vHCviEbF1ZvPDMqVc/uGXvLIba34jOJF6b1WG05r2ts80WVW+qOjJRO2FR2Vl3aPfqJayq3vHNcpJReWUo8fk1h2OpbqzZ9bLCx475fqr730ULVzu8P3fzpkmf7zu4cc9cZg3eNPNHZeeouE8l/emfNkisqthx+aKz0P+ttbPmXuka2fnfwluOGO6VDf1WI//Su0oi/3VW96ZMx2aGXdo69bYR6o1+sfmF6ryJbz+2suf9zVZ21k0zXh7fJ+fNnO0f/9pqqjqPHmfKRT3uVv/rmrqF/va264yOVKh/ww97qHzzSV7Xj3qGbzrdW3pk7x3zv49EdI05OFKbuLEfWz+yuefiKyo76QyPl960vR576UndNy3cqO/53ZKT81idLg392c2HI83cN2XR6zNr9/t+u/W4f5BhzReKF8/pkXpjeOea+WZWdB74z0jVlu4k8PXtnzbJvD9n2jrG27x3P2kFLr+wa/szNgzumVDvlA+/vHnLfz8qDX1hUteE8STq7a9wdH4nsGn5q5c7TOr3IcxfuqnnoiqoX3nao6Z280cT/dNnOEX+aN6Tj2P1s+e1PeVUPf3dX9frvVXdMHaTyiCsLw+/o9ay9avD284vGbD9n55gfXpjYecC7o7s+srUc+cvnumt+/c3KjknjI+Wjni1HfnNld03bzVUdH6425dqf9ibub+4btvWeIZvONUbOp14cc/uH413DTh/UeUa3Is/P6Kz52bert+3vODrW2lJfXd1Ba113Xa/jRI60tpyXtPV7VRsvjxnVXN05/M5tGlS6tnrjjKLMi+d0jvn+OdFt+78/0ffRF7zos5/ZNeqX36zaOnG8KR7zbDH6xJU9o1bdPGTzCdXWq/uFHdp8d1fllkXVm86JOib+yRdH33pComvIJ6Odn9plzYZPd435yecTOw45Ilr4wIZyNPP57lF/uLZy8/+MiXiHZ7zoz6/bNWrd94ZsPmWQ9Wq+Wxi+cLUZ1DOvYuMFJalzWufYe6bFd449btCuj+800dUzd456/OuVW1N1keK71pQqfjensN/KBVWbPzjUeAfv+dlbWLX57IhUeWbX6JuPHdRXdV5821kFOZvO7xz94KUVWyekYsXjX+tnb8+fAXf0Vd29rLe6697qjTPLxuk6e+fou89OvDj6A9Huk1+0ztpZXaMfe60/A35drnzwzu6hm+4etvlMp+xVn79r7M0TI72VF1a8cHavcTrO3Tn6gUsGbzko7ZQ+1GFjKy/uGvm7bw/e8q6xTinlFhO/uLZn2Jrbqjd/vFLe2IXFEff8sifeee+QzReUrHqmdY5eeFZV96jjzYuNncbJXbhz9KNfq9p++MGm53/yXuwPX941MjNv8JYP7OeUDlnaF//J93pHbLireuOnYtKQphfH3pqM9sYvHfzCOT3W2XJe1+jmzya2J98R7Tlxi42tuqhr5BNzE1uO2T9amvjncuUvv9U99Nnbh2z6aMLa/e/pqf7+Y8WqF78/dNOMsqfSWZ1j7myM7xj50UGFUzttpP3CrppH5gzectgEp/TefCn+xy8XRqy4sXLL+0dESm9Z1jPoZ3cUhz9/Z9WmM+LGDvtycextVV5v5PJBL5zbq8i2cztrFn9mcOf4o5yuk14wkbbP7Kz5zbVDth01xvZN+out/PU3uob+5dbqjimDVT7wB7uG/HCJN3j73dUbzrdW9tV+9r6c6Hj7W6PlY58rR/70xe6aP91Y1TF5hCm/7Q/lqocWdFc/d2f1ptPissOvLAy/Q5Je7c+AbTbyzGe7apZeXdnxjgMj5Xf8pRRZ9o1CzdM3VXacNCRSHv+jQmLxj0vDti2q2nCeMTJndY67/URn135nDN55+i5Fnvt0Z81D1w3dOsZx9H7PU6+ktXvz72KzN79YGDSd3fRuGfuN+QsXTL7wjAv3c+JOqzxN6ih07BhZPbL25jtvXvda1yy4e8G2oPuBN7PW1vXjHMe7NJVKXhx0C4DXb59bbtzcvXm5ZDfPnjbrqUjcecaxumHB3Qu2jagYMSLiOWv/3TVBtwNvdp5XjEoaFXQHALyqpulN4z4/fXq132sAAAAADFDWWieb3cThaWCA2+e2RQAMXCtWrDvA2p5bgu4A4A/DBYDQSCScPkn5oDsAAAAAAAAAAAAAAAAAAMBraDq7acasc2bd+Fqfnz19dsPsc2Y2vZFNwJvFnDlzojPPmXlg0B0YGHi1CAaES8+8dLCMnStH17/WNbZkx1nPedcb2QW8WWx5bsv+jjWPBN2BgWGfu3EZ9k290cLnZeWYsvmYpBtnT599nPXKp1ljdkn6wYK7FvzxX3/N7LNnniVjjrPG/MJUmeZ58+b1zpw2M+XImWmsLVljbp6/cP6qN/67AQYeU7YXSWZc07SZM42JLJNXTnvGjHKst1o2usFGbP38u+bfK0mzz555iROL3XvjHTdu5mfuzYmVCwwM1ntU0nZHzpLZn5w9xHreVdaY2+XpR8bae2bMmBH7x8s/c9bMydYxH/M8fUNWh9pOe/3np0+vdmR+ZI19xMq4kn24aXrTuIC+I2BAsdY8JKvtKjuP2rKttcbMN9buHzG2zTp2gqw97uVrjXNKuVTej5+5Ny+GCwwInuesltT73YXfXd0Z6+yLRKMfichGHKMJVmbMoOKgYf94vXUiJVml5OhoY8y8kcmRn+32KqcYoxZr7DYv4j0tI1dlnRDQtwQMLFGtkVFh/j3zc5JkrP48f9FNl31n4S351/ol/My9eTFcYMAZ7A0eXy6VfmWtafSMykba8a/XzFs47wkj52xHOtpa7w9bc1svlnS49exwx3OOdzzneFmTlfSXN/47AAY+K/Pcv/nskJf+gZ+5NymGCww4xrGfMNLD8xYuuEhlPSvpgH+9pumcpnNky4l5Cxc02YjeL2PPMir/VEZd8xfOv2L+wvlXyKrKWNMdwLcADFSv+neGccx2yYyVdt9tWtLBksTP3JsXBzox4BgT+ZlnvQebpjUdYmQjntEyI69JMsv3XGPL9klFdF/TObMuUNkMl9F18+66eXnTObM6ms6Z9XtZDbWyy+cvmu8G+b0AA8WoA0f9bWt+S0XTtKar5em3//i5ckV5udNttjZNm7XCeLZD0l+tpHkL+ZkDMIDMmTMnuudg2Jw5c5zPTfvcsFe77nPTLkjOmTPnn4boi869aPgFZ12w/xvRCexL5syZ48yYMaPytT4/+5Ozh7zax/mZAwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIDXzQQdAODN7dQPn3q2rHdiJBKb9cOHfrhVkqae2Hi1rPoWP9r81aD7APz3uOU6gEAV48WHJfNer1y6UZKmfviU06x0meQ8GnQbgNeHlQsAgZt64iknWZmHZM25MnausbqJVQtg4IoEHQAAbauffvawQ+oPlNEcSW0dhS1n5/N5L+guAK8P2yIAQsEau1mSZLRl2bJlpYBzAPjAygWAwJ124slHWzl3ytirZM059RMO29i2pq016C4Arw8rFwACNWXKlEpPzt2SfnL/Iw9cIZlr5NjrG6c01gXdBuD1YbgAEKhKr2KupJpSuTxbkjpt59dk9ZzKWjRnzhz+jAIAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACAfdb/BxbI4mMT7C4QAAAAAElFTkSuQmCC",
"prompt_number": 14,
"svg": [
"\n",
"\n"
],
"text": [
"Plot(...)"
]
}
],
"prompt_number": 14
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 15
}
],
"metadata": {}
}
]
}