{ "cells": [ { "cell_type": "markdown", "id": "dbcd8f8a-db6f-41d3-99f4-7f131a7304f5", "metadata": {}, "source": [ "(aggregation)=\n", "# Aggregation of various sources of information, expert elicitation, or evidence" ] }, { "cell_type": "markdown", "id": "7404277c-c81b-47ff-b267-7958e0d9c841", "metadata": {}, "source": [ "## Installation\n", "\n", "Install the [pyuncertainnumber](https://github.com/leslieDLcy/PyUncertainNumber) library from [PyPI](https://pypi.org/project/pyuncertainnumber/).\n", "\n", "```shell\n", "pip install pyuncertainnumber\n", "```\n", "\n", "```{important}\n", "A virtual enviroment is recommended for installation.\n", "\n", "Follow the [instructions](https://pyuncertainnumber.readthedocs.io/en/latest/guides/installation.html) for additional details to install `pyuncertainnumber`.\n", "```" ] }, { "cell_type": "code", "execution_count": 1, "id": "06aef6b3-95b0-483a-b142-dddcf2c59307", "metadata": {}, "outputs": [], "source": [ "from pyuncertainnumber import UN\n", "import pyuncertainnumber as pun\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": null, "id": "c34c48ce-7d69-4efa-85b2-a1788d05eb5e", "metadata": {}, "outputs": [], "source": [ "from pyuncertainnumber import stochastic_mixture, envelope\n", "from pyuncertainnumber.pba.dss import plot_DS_structure_with_labels" ] }, { "cell_type": "markdown", "id": "d7444f84-b01a-4c63-ab2f-370b2b3c893c", "metadata": {}, "source": [ "## expert elicitation is represented as probability distribution" ] }, { "cell_type": "code", "execution_count": 12, "id": "3aaf3b1d-49a9-423c-b7f4-73b9c548cc28", "metadata": {}, "outputs": [], "source": [ "a_d = pun.D('gaussian', (-5, 2))\n", "b_d = pun.D('normal', (1.5, 3.))\n", "c_d = pun.D('uniform', (1, 12))\n", "d_d = pun.D('gaussian', (10, 1.5))\n", "fixed_value=0\n", "p_env = envelope(a_d, b_d, c_d, d_d)" ] }, { "cell_type": "code", "execution_count": 13, "id": "a6b89059-a9fe-4609-bbfa-db60cb515961", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgUAAAE6CAYAAACGWCHoAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAY7pJREFUeJztvQecVNX5//+Zur3DNnaBXXovuxTp1QKiogJGk/zQKJr8v4kdxRSjJkEIRhMTI+j3a0k0CoixEI2AgIBSliK97VIWdhdYttep/9dzhrsO48zuzO7szi3P+/W6rzkzc+bOmbnnnvOc5zxF53Q6nWAYhmEYRvPoQ90AhmEYhmHkAQsFDMMwDMMIWChgGIZhGEbAQgHDMAzDMAIWChiGYRiGEbBQwDAMwzCMgIUChmEYhmEELBQwDMMwDCNgoSDI7NmzB9OnT4dOp0NBQUGom8MolIqKCtx///1YsWIF5syZE+rmyBK+1xgm+LBQEGSGDx+OdevWhboZzBV8TRY02fbo0QNy5b777kNOTg7mzp0rJj8SEtRAMCdvX/daa66tr8+0tb1K7X+MdmGhgFE1q1ev9vp6dna2mFTk3O5p06YhPj4e+fn54lHN1yOYtOba+vpMW9ur1P7HaBcd5z5oH0ilSYM53fxM6CD1shI1N9R/ysvLVSMMtOf1aM97ra3tVWr/Y7QLawoYv/BHfd3eKu5Az0+rtPXr17dbexio+nq0tb1K+70MQxjV9jfQXp3E7t27sWTJErHaon3ZJ554QtyktKqgx8TERLz//vvCoItUtXQTL168WNRdsGABli9fLs5Dn1u6dCluv/12rFq1qtnvaQ46P+0xUj2a4Eh9SN9LUHvoHPRIKwtpL5LOTasN+m5/fqe//w99/+XLl8XnPKHfm5SUJMpUh8q0Ctu1a1dT/ebqUHvpffoN1Db6TfRf0n8nreZaar8/bWjut0j/o3QuYsSIEeJ/9OwL7ivM5q6RP33IH1r6Djqf1G6qs2jRomavra//kl6n/5zaSeeX+hXZKlC76b+gaxVIv2vu++gcvq47vebrerREIH29uWsr/Ubay6f/nf4T+u+9faa5/uPZLjn2P2obvUeUlZWJa8AaC8YvnCpi4cKF4pDYvXu3Mzs7u+l5eXk5bZU4lyxZclWd+Pj4puf5+fmiDj16ntvf7yE8z0Hf6f4ZYsGCBc7ly5d/77vdXyPo3KtWrQro+71BbXA/Dz2fNm3a936n+7lvv/32q/4vf+t4/hb6rdROf9rvz/n9+S3S9fZFa66RP32oOfz5Dm9t84U/fYH+F8//k35HoP3On+9r7rq3dD1a+/u8/Veer9Hv8Owfw4cPb/YzzbVXzv2Pnnueh645w/iDaoQC6WbxHBxoAJEGJYLqrFu37nufcx8kPQdRqi+dN5Dv8fyM+3e4D6DNtYWgQUG66f39fm/Qb3IfHKi+56DlORHQ4OM5ePpTR2qnZ5v8ab8/5/fntwQyKPt7jfztQ94I9DtaEgr87QtUj/oPtZn6kuf3+9Pv/P0+X9fd/T1/ac295us16fdL0Hv0mvtvDkQokHP/o3uH7hf339LS2MAwEqrZPpD27kjFRoc7pI5zt/R1V9d5U0OSKo5cwiR1IJ3PXc3v7/d4ts3zu6R2SOpdX+Tm5gpVIp1f+s5Avl/CXb1JdfPy8r5Xhz7v7kZFalFPAy5/6nj+Rgl//j9/zu/PbwmEQK9RS30oGN/h7/la6gv0faQ+pngH9Ohve937naRS9+f73H9TW2jNveYNyaWTfo97+8iQs7XIuf/RFgVtKdA2Cb1Hn/O2Tcgw3lCNUED7ZoTnHqU/e5ae0GdIKJDcwtxvwGB+T2toy/fT4EV7kvR75s2b53UCIoGIDul8NBht2LAh4Dq+Bix/2u/P+f35Ld5+v5q8QQLpC/S7aVKkySEQwaO13+eP4NHS9QjWvSadp624t1fu/U+y45BsVEhAUKM3CxN8VON9IN2U3izUW2MVLxkarly58qpBqDXfI61oPN+XVj/uKxhv0CqEbmYaUNryO8nATDIs81xluWsgyIiJHiWDMc+BxJ86vvCn/f6c35/f4nne5izB23qN/CHY3xFIXyABlwzNaIIko1l/cO93gX5fc/hzPYL5fdL/6i2QUCBeNe7tlXP/IyNDSQCRjKNpPGNPCEZTQgHdANTxyXvAczBsbqUg3YSedcjim24iz5u0Nd9Dn1m4cOH3PkNChzdLas+bl+q89tprbf6dnipUacCh1yX1J6nq6XUaTHytfvyp4/69nv9FS+1v6fz+/hbpf/VnpRjoNfL2G1v6nrZ8h6/z+dMXSOtC30uQxoXqe5sgmut3gXyfr8k2kOsR6Pc1B30v/X7JC8B98mzuPL7aq4T+57ld4C7cMYymghfRjU43iXQT0Y0ruR3RDUbv02RDAyW57Hi+5j4J0R4s3VzebqZAvkc6p/tnaCXs6eJE7yUkJIgVHZXpJqdVMrXDc3L09f0t/Td0Plrl0G+igwZYaov0HbTCcA+/Suenc9OAKrWhpTqSKxqdW3J7o8Hd3/b70wZ/fgtBr5H6lNrg6xqRACh9d3PXqDV9yNd18Pc76D1pQvfnfO7/JWkE6FyStkX67VOnThX1KISyNHn42++a+z5/rru369Ge97TntaX/W+pb/vQHX+2Vc/+j6y2dg96n6ym5oDKM5oQCJSMJBaHa+6PvpwGNBhka5KQJggYjsrGgQZAGl5bqtGVF4k8beMWjrn7HMIx8UI2hIdN2pIx80qRLE4SkdqRVpLTH31Kdtkza/rSBhQKGYZj2QTU2BUzbIZWnFLHRE9onlaKrtVSnvdvAMAzDtA+8fSATPPdj3fcZOxJ3FyYJ2tukvUpphe5PnfZuA6OufscwjDxgoYBhGIZhGAFvHzAMwzAMI2ChgGEYhmEYAQsFDMMwDMNozyXR4XCgqKgIMTEx0Ol0oW4OwzBMmyGzsOrqaqSnp0Ov53Ue0zY0JRSQQJCZmRnqZjAMwwSdwsJCZGRkhLoZjMLRlFBAGgLp5omNjQ11cxiGYdpMVVWVWOxI4xvDtAVNCQXSlgEJBCwUMAyjJnhLlAkGvAHFMAzDMIyAhQKGYRiGYeSxfUAJdijMKoW09Uyz6g1KB0uhbulz9MjpQBmGYRhGBUIBxV2ntK0U6z4pKanF+hT73j1XuZRNj2O1MwzDMIzCtw9ocqeVvr853CmtriQQENOnTxdaBoZhGIZhVLB9EIhWwZvwQK8zTKiCYdEhWX6z9TfTlgBEUm46CkDEQYiYUKEYoYC2GTxJTExEWVmZz880NjaKw92fl2H8GZjdsdqtOFt5FpfqLuFi7UXxWFpXigs1F1BSWoJfDftVSNrLAOaqCuitFuhtVuitVhiayhbYw8JR0aMf5AZ1sfJyA0pKTLhwwdT0+LOfXYTZ7ETfvn1ZKGBChmKEApr8SQhwhzQH3oQFicWLF+OZZ57pgNYxShYAai212HZqG/656Z+oCq9CWWMZrsu4DrO7zxbvH688jtvW3+b1PAadAU8NfQp6HQ/iQcHpQMKJQzDXVENvaYThygRPE77eZkFjXCLOj53eVD3zq89gbGzweqqGuERZCQWvv94Ja9YkCAHAYvl+f7nttnJkZlpC0jaGUZxQ4CkQECQQNGePsGjRIjzyyCPfi/zFaBsSBrYe2Iq/7vsr9pftx4mqE7A77VfV6R7dvamcGJaICEMEEsISRJke6UgKS0KCOQE6gw7JnZJD8EsUiNMJffllGC8UwXCpBIbKMjiiY1E3eUZTlbj/rIS+od7rx3VR0UhOdvuvY+Nhb6iH0xwGp9l85TEMMJnhiI2/uq4f2BvtOLP+DLJmZLV6O+jMGT1efx3YuDEar79+ComJrr5VV6dHYWFYU71OnaxISbEiNZUebTCZnAgLC2MtARNSFCMU0OTvuVXgTXvgDt1gdDAMCQLl9eWID4+H3W6H0+rE6lOr4YDLJiAlIgWDEwejk7UThvcajuHpw5GanComhlRnKs797JxPuwG2J2gBpxPGHV9BX1QIfck56Ny29AhHfCLCw8ObntuH5IKmUWdEFECTfFi465Em+4iIq+paf3i/z6+lK/JdTf/IW5GHLb/egp6zemLWP2cF9Nl9+3RYtiwMa9ea4HC4+sO2bdGYNasSZrMZP/95DObObQStSzIy9AgLo34TBp1OamUC9yUm5ChGKCCvA8+tAnru7o3AMJ7YHDa8/e3beGn7S4jWRWPFuBXidVrxLxyyED1TemJE2ggkGhPRu3dv1NfX44WiFxAVFRXqpisXUvWXnIcj84q2RaeDoeA49BeKxFOn0QRHajocKV3gTOoEZ0Knqz5uGzMlFK1G3aU67PjjDlHOvj7b78/RWuW3vw3H22+b4HS6JvQxY6pxyy0VuOuuTkhJcWUvpMl+0KB2az7DqF8ooPgFFKRIClBEj+RtIAkC69atw5IlS0LcSkaumoE1R9bgqQ1P4XjZcfEabQFcbriMpHBXTIx7Btwj1Ms0WNfW1gqBgGnDtsC50zAc/haGE0cAuw0NDywkdZ1425o7Brr6OjjSMuDslEIm9pAbX//+a1iqLEgekoz+d/b36zPl5U6MHRuDoiLX77nhhgosWHAJPXs2IjIyEmlp4bzyZxRFyIMX0cRPjzT5SxM/BSQi3n//ffG+JBSsWrUKTzzxhKhLWgIKZsSBixhPTpWfws/W/gyf538untO+/92978ZtWbchKSoJKSkp4nVW1QYBiwWGw/tg/HYX9GWlTS/Tfr6uuhLOMNeevqP3AMiZSwcv4eBbB0V54vMTodPr/BI8GxpKMHasHbt2ReGZZ85j9GibiM4q9S3uX4zS0Dm9+WCpFDI0jIuLQ2VlJWdJVCnbC7dj2tvTUGurhUlvwj297xECQXxkvBAGfA3UpClIT08X5aKiIt4+8AP9mXyY//MBdFeMAp0mM+x9B8HebzAc6Zli20AJ0BD4wawPULilEL1u6YUb37qxxc/YbE5cvFgMm82G+nqdcDNMTAxrEgg6Eh7XGM1sHzBMIAM7BRKKqI5AakQqYswxeDbnWfRO7N2sMMC0HkdSZ8BqgSMuAbbho2HvN6Rpu0BJFPynQAgEhjADxj87vsX6y5aZsXmzDX/6k1383Lg4kxAGJLsBhlEyLBQwiofiDBSdKRKBqgwwYMX4Fegc2RmZXTJZGAiyZkBfeBq2cVNdL0THonHu3XAmp8nSRsBfIjpHoPOgzsi6Ngtx3eKarfvWW0Y895zLW2DjxhjcdFMDevXqxX2MUQ0sFDCKhtwMJ7w+AdPSp+EnfX4iXusS2wWpqS53QiYI1NbAtPm/MB5z7bk7umXDkZklys7ULlA66SPTcefmO+GwutxTfbF9ux6PPhohyg88cJEFAkaVsFDAKF4gOFh2EEW1Rbiz153omdmT1bhBRH/iCMzrPoausQFOnQ72oSPhSHbZXqgJvUEvDl+UlgJ33RUOq1WH6dMr8dBDFSwQMKqEhQJGkdRZ6jD1/6YKgSDeHI//m/h/6N2td6sHaQqI8+mnnzaVNY/VCtPmz2E8sEc8dSSnwTJtFpwpaVALm5/ajLD4MOT+PBfGCN9DocPhxE9/qkdpqRHZ2Q14/vkL6N2bBQJGnbBQwCgOu8OO2/55G/aW7kWsKRZvTHoDk/pPatMgbTAYMH58y0ZmmsDphPmDf8BQXAhyTbKNGAvbNZPpT4JauHTgEvb+fS+cDicyJ2Siy2jf2yCrVxvxxRdRMBqdWLasBMOGsUDAqBcWChjFeRk8+OGD+Lzwcxh1Rvx13F8xqV/bBALGA51OCAL6DWthuf4WOLr6H91PKX1o86LNQiDoPbt3swKBy6ulDOnpJtxySzlmzcrgvsaoGhYKGEWxt3gvXjn4iij/bsTvMGvwrKAkkLFarXjzzTdFef78+TCZTNAcDQ20dyKKjh590NAtGzCq73/IX5vf5II47plxzQoEJSUlGDPGio8+qkJERBgMBk58xagbFgoYxUCDdFR1FH4z/DcoqCrAgtELgpZRzmKx4LHHHhPlO++8U3NCgXHXNhj3bne5GMZfSTKmQoHA1mjDV7/6SpRzfp7TrAsi9TcSFomYGCN69+74wEQM09GwUMAoAhqgT5w4IQbp27NuF5N2MFPMkk3BzTff3FTWVAbDrRtgytsmnhpOHoUtdwzUyr5X96HyVCUiUyIx4uERzfa3hx5yoGvXBMyeXYG+fVkgYLQBCwWM7KEB+t2v30V2ZDZiTDFi0g52HALyOHj77behNYxfb2wSCKzjp6taILA12JD3lzxRHvf0OJijzT7rHjqkw9tvJ8DpTERurhVDh2pIUGQ0jXLDkDGa4VTZKSzYtAC3rLsFxfXFSEtL41VbEDBu3wzTzi2ibJl8g6oFAsIYbsQd6+5Azi9y0P8H/ZsVQp99lvIZuGISzJjBxoWMdmChgJE1lM9g/gfzUWerQ0ZkBoZlDwvqtoFWMeZ9DdM3m0TZOmG6CEqkBeKz4zHhuQnNZkE8eVKHL76IFuVf/KJCW9tJjObh0ZWRLbRie3Xzq9hSvAVmvRm/G/k7mNrJ+I2yJFKmOTqorGpsNhiO7BdF69gpsOWMUX0/unz0st91KdERaQkmTKjC9dd3YS0BoylYKGBkS721Hn/I+4Mo39PnHlzT+xoeoIOB0Si8DCzTZ8E2Uv0Bm/I/zcfbo97Ghoc3tFi3vBz48EOXR8I991SyloDRHCwUMLLlpe0v4XzdeSRHJGPRhEW8bdBWbC73OkFYGOwDh0PtuLsghie2HL76rbfMqK/Xo3fvetx5J9uuMNqDR1lGlpTWluIPW1xagocGPoTY8NhQN0nZNDYi7J0VMO7+WrghagUKZVx5uhJRqVHNuiBKWwdZWZeEceH8+ZdhaCZBEsOoFXZJZGS7dTApbRJOVZ/CLdm38IqtrbkMvvgI+rJS6PZsh400BGHqT/pUe7EWO5ftFOWxT49t1gVREgqGD68SR1hYGPc5RpOwUMDIDhqc6y/W4/mRz8NityAtldW4bcGYtw2Gk0fg1OvReONcTQgExNfPfQ1LtQUpw1LQ/w7fLohSn7tw4ULT86ysLO5zjCZh/RgjK2hwzs/PF2GHiajwKB6c24D+/FkYt30pytZJN8CZlgEtcHH/RRz8x0FRnvj8xGZdEImKCmDZskScPm0WWgI2MGS0CgsFjKw4V3kOj255FGeqz7RL5EJN0dAA0+droHM6Yes7GPbBOdAKNcU1iEiKQJ/b+jSbBVHigw9MeP31znj44a6sJWA0DW8fMLLSEjy37jl8fPZjFNUV4b8/+C8Pzq3F6YTpy7XQV1XCEZcA65QZIiWyVsi+Lht377kbdovdr3739tuu/2b27HLo9Skd0EKGkSesKWBkQ0V9Bf51/F+ifE/fe1iF20YcqelwGk2w3HCrcEHUGmFxYYjsHNlivTNndPj22wjo9U7ccks9C6KMpmFNASMLaLW2ZP0S1NhqkBWThTlD5/Dg3BZ0OtiHXwN7vyFARMsTo1o4+PZBmKJN6D27t1/9h/rdP//ZCCAWOTm1GDmyK/c7RtOwUMDIgkZbI944+oYo/6TvT9otnLHqoRgEdruIWijQkEBALoibn9osPA6MEUb0uKGHX0LB55+78hzMmFELgyGqA1rKMPKFtw8YWfDewfdwsf4iOod3xt25d/NqrZVQToOwd5ZDX1QIreHugkg2Bf5w9qweBw5Eiq2De+6J537HaB7WFDAhh1ZrL257UZR/0OMHiDR3/OqW3NDefPPNprIiqa2BadPn0DU2QH/uNBzpmdAKgbogSv1ux45qxMREoW/fenTpEtEBLWUYecNCARNyLDYLRncejdL6UsztOTckqzWj0YjZs2dDyUgCgSM5VfWZDz0n982LNgNO+O2CKH3ummsqsHlzJaqqyNAwq93byjByh7cPmJBCA3PhmUL8fMDP8dn1n6F/t/6swm0F+jP5MB4/BKdOB8v0mwANeW6c/OQkzm09B0O4AeOeGRfw500mJ0aNSud+xzCsKWBCjcPhQENDgyibTeaQZUK02Wz45JNPRHnWrFlCc6AY7HaYNn7mKg4dCWdyGrSC3WpvyoKY+4tcxGbG+i2MFhRchMnkCt+gqOvNMO2ILO6EpUuXIjs7GwUFBeLx9ttvb7F+fHw8Kioq/KrPyBMamP93y/+KkMbj08YjPTl0q7XGxkbMnz9flIuKihQ1SRj3bIe+/DKckVGwXjMJWsJgMmD6X6Yj7y95yH0wN6C+t3BhCo4ejcDTT1/EgAHa0awwTHOEfOS7//77MWfOHEybNk08pzJN9MOHD/cpECxcuPCq53v27PFZn5G3lmDpnqXIr87HM7nP4MGuD4asLaShGDfOpXoOlbaiVTidTZ4G1vHTNZPsyJ2uk7qKIxAaG4FvvolGXZ0B/fp14q0DhrlCyEe/FStWNAkExPTp07F8+XKf9Xft2nXVc/rs+vXr27WNTPuw/dx2IRCEG8Lxw6E/DOnAHBERgbVr14qDyoqBbAhumofG2XfB3m8wtAS5H7aWb74xCoGgUycrRoxgLQHDyEIooMmctgG8ve4L0gqQICGRl5fHWgIFQurbV7a9IsrTukxDQmRCqJukXHQ6OLr31FRug4vfXsRr/V/Djj/uEH0pEKj+v//tEijGj6+BycRCAcPIQiggmwBPEhMTUVZW5vMzTzzxRNOWAwkIdA53TYPnPnFVVdVVByMPGqwN+PT0p6J8S9YtylLZywGnE4a9O0QmRK0hXBApcmGVBZePXg5Yw0Sf37jRFblwyhTOdcAw7oR0JKbJn4QAdyQDQl8sWLBAHKtXr8bUqVN9CgTE4sWLERcX13RkZmonmIvcWXtiLaqsVUgOT8ZNA28K+cBcW1srbFnooLLc0ecfhXnT5wj7x99dYY01xMmPXS6IFMp43G8Dd0E8dUqPs2fDYDQ6cMcdSSHvewwjJ0IqFHgKBAQJBN62FCRo64DsDvLz85Gbm4ucnByf2w2LFi1CZWVl01FYqL3Qr3KEVmrLt7vsRmZ2nSmbPAeXL18WhyJcELducBUHDNFUTAJbgw1f/drlgpjzixy/XRDd2bTJZV89eHA9EhJCbmvNMLIipHcETf6eWwXetAcS5LK4bt06rFq1SjynMnkf0JbC7t27v1efwtUqNmStirHarSiqKRLl2dmzeaUWIIaDe1wuiBGRsOWMhZbY88oeVJ2pQlRaFEY8OKJVAmmXLpdwxx3R6NWrATqddmI6MIzshQJS/XtuFTRnI0A2BCNGXD0QkHvi+++/367tZIKLQWfA+1PfR0FVAcb1HsdCQSBYGmHavlkUraMnkuQLrVB7oRY7X9gpyrRtYIoytUoo6N+/Why0YOC+xzBXE3LrLgo85K7+p9U/GRK6CwJkP0CQsEDve2oP5s2b14EtZtoCDcqnTp0S5ezYbB6UWxGoSFdXC0d8IuyDcqAlzn99HvZGO1JzUtFvbr82n6979+7c/xjGg5BvqNFWAKn/aXInLQEJBO4uhqQFIMGAhAfabliyZImo36PHd7nS3YMZMfKmqqEK5TXliDBGiKiBPCgHQEMDjHu+EUXbmMmasiUges/ujc6DO8PeYPcrC6I3duwwoKQkUtgTsMcLw3wfnTNQJ18FQy6J5IVARoexsYEbKDFtZ+nWpfjtpt/iJ31+gmenPQuDTCY28jhIT09vCnMcFeVyWZMVdbUwbV0P/cUSNN61QFNxCYIBDXWzZhmwZUs0Fi4sxuLFKaoQDHhcY1SlKWC0NSi/s+8d1NvrEWeOYy1BoFBug2tvdrkgaui/K9xSCHO0GSnDUtp0nsZGJ/LyIkV53DiLpvqf3W6H1WoNdTOYEGE2+59sjoUCpsMoKCvA/sv7oYce13W9TlODclCRiXalo1wQv/jZF6gqrMKsf8xCz1k9W32uvDwD6uv1SEy04dprtZEqmQTxkpKSZmO/MOpHr9cjKytLCActwUIB02GD02vbXhPlnM45GNh9oCYG5aDQUA/Tps9hyx0LZ6dkaM4F8WwVotOj0W1Ktzb1v3XrGgHEIDe3FkZjNLSAJBAkJycjMjKS7zmNJp4rKipCcXExunbt2mIfYKGA6RBoUF57eq0o39D1BlXs5XYUZFxoPLIf+tKLmrIlqCmpabMLonv/27XLlUFy5MhG6PWxmtgykASCpKSkUDeHCSGdO3cWgoHNZoPJ1Px9xCMz0yGcKj+Fg+UHxdbBvMHzeMXiLw31MFKOA4pLMGqCZgQC4uvnvoa1xipcEPvO6dumc5EZxr59LnuCG2/Uhj2LZENAGgJG20jbBiQotgQLBUyHbh3kds5FWgxHkfMX476d0FkscHRKgaNn2yZGJXFh3wUceueQKE9aMqnVLogShw4ZUFtrQHS0HUOGaGvY04IAxASvD/D2AdMhQsHE5IloHNiIHnE9ZLl1QJL0smXLmsqywGoVQgFhGzFWM1oCkQVx0WbACaEhSBuR1ubzJSYW4913jSgpMcFk6hK0tjKM2mChgGl3aFDOjM4UsQlSU1NluXKhfbb77rsPcsJwaC909XVwxMbD3nsANMMVYYA8DlqTBfF7pxOhWKwYNMiKnBy7LIVSJnRQRF3Kt+MeNE/L8N3BdFhYY0IuwYpkj90O4+6vRdGWO4Z8iqAVaKtg8D2Dcc++exCTERPUc1P6dDkKpUzgtMbN0ttnKEou8x2sKWDaXSh4df+rSApPwvTM6bJdpZEBztdfuybhMWPGhF54cTphHzwCOHYA9gFDoaX+Ik3aemNw+sqZM3o891w6Ro2qxYMPslBaa6n1+Z5Bb0C4MdyvunqdHhGmCJ91o8z+RwWl/DbuGXMXLFggXlu8eDGys7NFOPzp06c3hbqnsPhz5swR9ShJHj13D3nv7XykEaDPvPbaa9i1a5fwyCDtANWjcPp5eXminjdWrFgh2kG5d6gdvpL2qQEWCph2pc5Sh1cOvyKiGI7KGiXbVVpDQwNuvPFG+YQ5NhqFHYHQEsj0P2sPF8QPb/0QoxaOQq+bewWlrwhX2LV1+OCDNJw6FYZHHmGhIHqx7xgNM3rNwNo7Xa7DRPKyZNRZ67zWndhtIjbN39T0vPufu6O0rrTpufNp/yLoU26b5cuXNyW7o0l37ty5It8NTdo0kdMKn16XJn2aoGliJoGA6hE5OTmiPm0FeDsf1aeDhIBFixYJQYLq07koqV5z2wdLlixBfn6++HxCQgLKy8uhVlgoYNoNGpDf3f6uEAhSI1IxJHkI5ApNQH37uqz7ZSW4yKktHeCCWHqoFLtf3o1eN/UCdMHpgzt3ula+I0Y0QK9n9zy5Qat00gBI2XKpTKt2moBpwqYkeVOnTsXu3bu/91mqK0ETOwkD9Blf5yOk9wOxIcjPzxeCiaSRoLL7d6sJFgqYdoMG5C/OfCHKU7pMEVkR5Qr5cu/Y4YoHEGpMGz+DvUtXOHr204wtwYW9bi6Iz7fdBdEdKT7BtddGyUvgCxE1i2qa3T5w5+JjF5vdPnDn9IOnW9UemmBpxS9N2p6qeZrECZrkm1Pb0yRNWoKWzkd1fFFQUND0fe5QZl4pg29zn1cD2hhxmJDgcDqwuWSzKN/S7xYekP1Ad7FYuCGaP1sD1PoevNUmPG5atMnldTC37S6I7pSV6XDmTJgoX3MN9z9pr9/X4W5P0FJdd3sCb3X9hSZbSdVPeK7I6diwYYOo52ko6P6c7A5oq6G583mDJnnaUqA6FV4MEck+gV6XhAXJVoFeVyMsFDDtxq6iXShtKEWUMQrjMtruWqYFjHtc2gp7r/5AjPpD8RIn/n0CRd8UwRhhxLing9tPdu92rXy7d29E585sTyBHaPVNE/nSpUvFREsaAZqAqUx2AqQBkA6a9MkGQYIMBqkefZbeI62Ar/PRo2S/4C4kUF0SKKjucC9bCpIxIn2e6pA9AmkO1OrCKF99LqN4Pj72sXgclzoO4aarVyByo66uDpMnTxbljRs3hiY0bE01DMcOiKJt+GhoJQvilt9sEeXch3KD7oKYl+ca4gYNqoNeHxfUczPBQzIW9HzN/XVvNgW+PAG8nY/qkW2AJzS5k6Dgi+wr3g9agYUCpt1UwgfPHRTlSWmTZL91QO09evRoUzkUGPfnQedwwJ6eCWeqNqLu5a/Nd2VB7BKN3F/kBvXcdB2PH6fMiOEYPLgeAAsFDNMSLBQw7QINyItzF+PRgY8iNjxW9kJByLFZhVAgisO0oSUg+tzWB+EJ4bBb7TBFtj4Loq8++Nxz5/Dww+RiauY+qCKkrQBawefm5qrWEyAUsFDAtAvSartTeCd06dKFB+QWMBw9cCWkcZymEh8R3aZ0a9fzJyba0adPy3nkGeXgayuAaTssFDBBR6ht8483PefBuGWcMXGwp2XA0Usbbojl+eUwR5sRldIxQaLkGkmTYeQG3ylM0LHYLBi3Zhzu/epeVNoqWSjwA0e3HrDc8RNNbB2I+BU/+wJvDH8DBZ8XtNt3PP64Hffd1x1bt0ZzH2QYP2GhgAk6Xxd+jbLGMhyvPI4+mX14QA4EDaxoj394HEXbi+C0O9F5UOd2Ewq2bo3E9u3RsFjYnoBhOmT7YN++fcJ/Uwr5SK4b3bt3b8spGYVDg/HKPStFeUzKGBgNvEPVHLqKchhOHoZtwDAgQv0heG31Nmx9eut3LohdguuCKFFbC5w44XKDnTkziYUChvGTgEdsiixFPp2UDpeEgKysLPE6CQYkINDrZA1KwR1YQNAeIjrdOVeSlAnpE3gwbgHD/jyYdn8N/bmzsNzyA6idPX/b024uiO7s22eA3a5DcrIV3bpx0CKGCbpQUFlZKSZ6mvBXrnStBH1BgsGrr76KTp064bHHHvO7MYzyKa4uxtHKo9BBh5sH3sxCQUtuiIf2uoqD1BkdzZ2a4hrs/NNOUR7/zPiguyC6C6abNzcAiMHgwRS0qH20EYx23B/vvxL1UK1RDAMWCkggoBzUNNH7A2kPnn/+eSEcLFu2jAUDDfHf/P+Kx/4J/ZEclRzq5sgaw7FD0DXUu9wQs3pB7Wx7bhustVaR26DP7X3a7XtIKNi715XvYNgwCwumTIuTPuU/8DXhT7sSOlkr+GXVFBcX16qJnYQDFgi0Aw3GH3z7gSiPSxnHg3ELSMGK7INyVW9gSH0jIiECBrMBE5+f2O594+BBV7KeadM4cJYnTkuj78Nm9b+utfm67Y235EWt+cySJUuC1CJ1wFZgTFAH/tGdR6OqsQqTMybzYNwMupIi6EvOw2kwwDZwGNQO9YUJv5+AnAdzEJXcvrEJGhuBbt0sqK/XY/hw7oOe2BY/5fM9Xa9+MN5573d1l/0WsFq81+3WA8b5P/uu7p9/D9TVNj03Pf2C322iRENS9kFiwYIF4rXFixc35R6gPAdk0E6TOCU0ogRIVI/SJEsJjhYuXOjzfKQRoM+Q1psSKSUlJTUlO3r//feRl5cn6vmC6tD3SOdtrq4mhYIHHnhAJKigP5fYu3cvEhIS2LhQ40LB7O6zxZGenq4oocBkMuHJJ59sKrc3xv2u+8beawAQ2TEBfORAewsERFgY8Prrp0FBNePjtRUdUolImQuldMc0+c+dO1ckNaJJmyZyWuHT69KkT4ICqfVJIJCSH1FGRapPWwHezkf16aBJnTId0gRP9elc8+bNa3GLYITbd5GNAQkZ3pIxaVYooAtAF0ti2LBhwjOB/ugpU6YEq32MggQCsiGRUJJAQJjNZjFQdAg0W+l0cOr1sA9pPwt8ubggfvE/XwhPg+QhyR3SDy9cuCDK1AU5kuH3MS76g+83Pf4v42O/9V1X51H3wV+2qj20AicNAE2yBJVp1U4TLk3YNAFPnTrVa5ZE95wHNLGTMECf8XU+Qno/UDuBeLfvovmPtBcsFLhBf6znwE8Xbs2aNQGfi/Je0/lIoKBHb2kv3aF6dPFJciPownBCjNBCg/HnBZ+jT3wfZMZm8mDcHDodrNNvgnXcNCDctfetVnb/dTeOrT6G4p3FuHvv3dAb9e3eD6uqbIiIcAl6ShNOOwKdOSzkdd0hLQCN5dIE6znR0pxAtLQypzmAtAQtnY/q+KLgyhykZfRtUfnQ5G0wGHDdddcJLwMSCKTtBH8hKZAkNjoXqYZIaqRz+0LaS6J9JfoM1V+xYkVrfwYTJMiO4KHtD+Haz65FQ1iD4gZjh8OBI0eOiIPKHQIFK1LY/xSoC+KuF13jwdinx7a7QCBxyy29cP31vVFTw0mQlADNAZKqn6BJXbIRoDIdpIWmep6Ggu7PaeVOc0Nz5/MGCQm0pUB1KpoxXnR/j87vrilXE63WFJBUJhlc0AWjP4lcFiWVjb/QhE6rfgna/6Hn7q+5Q7ES6KJLkMpX65KdHNh8ZjPsTjsyojLQNa4rlEZ9fT1Gj3blHSgqolS77bP3rSu96No6SGqf8L5yYtuzbi6It7WfC6I7ly/rUFRkFuWuXb0byDHyghaFNKZLGmOCFnxkLEjjPc0FNN/QIS0IJdW/tAiVFouSVsDb+aR0y3Q+Oq/0HtWl13r06NFks+CJpF2QDBhpnlLj1kHQvA9o24AOik0QSFwCukje1P7NCRZ0UdxdSLTkPypXSGX772//Lcqjk0crduuArJHbG9PXX8KQfwyWSdfDPmwU1ErJnhIcfvewKE9aMqnDVuwUyZDo3r0RCQnsXKUUvG0Z02vur3uzKfA1OXs7n690yzSH+FqESrT0vppo9V1D1px//OMfxUUZOnRok6FHIKt2b6oaSZXjDWlbgaRCKrfkGtLY2CgOiaqqKr/bxgSYfKbIFc9+bNpYRapsSTPQnIoxKNRUQ1/gSint6Kpe7ZaIKLhosyj3u6MfUnNSO+y7JaGgX7966HSxHfa9DAOtCwUU0Ojxxx//nqAQiFBAk7qn0QdpDnzt67jvM0mSIAkldA5vkiH5uD7zzDN+t4dpHSU1JThZdVKENp7Rb4YihYKOgEIa65xO2NO7qnr7oOCzApEF0RhpxNjfjO1QYWTHDpso9+9fD4CFArUibQWQHQGF3mdD8+ARVP3abbfdFlB9b1agUsbF5nDfMiChgCZ/b0IB2Rs88sgjV2kKMjMzA2oj0zJfnvpSPPaN74vOUeqd7NqE0wnDQVeeA7vK8xxkXZuFKcumwG61t1sWRF9CwaFDLgv4IUPsLJyqGF9bAUzbCemmm7uxYnPaAwlJC+GpjfCl9g0LCxMH0872BPu/sydQ6kBMhoaSYEl2KxHk0xZE9GcLoK+qgDMsHPbe/aFmyMtgyH1DOvx7y8q+MzK8/vpkxfZFhtGsUEDSnjcXE19WnZKGwN2X9PLly+x9EGKh4KH+D2FM5zHoGd9TsQMxuSFu3eqyi2gPl0TjfpeRlL3fIMDY/hETQ0FDeYPYMjCGhWZYaWggA7MylJcbkJgYHZI2MIzSCbmZuOQqIkGuje4uh7RvRCs3CXIZca9P73dYJDrGK/Fh8bg241qM7aVMI8N2x2aDvrjQVRyUA7WyadEmvD3ybRRucf3WjhZODYYSPP10EV56qZD7IcO0kpD77JChCPmMSoEjpGBGElIwI0m1S+6IVJ8OciEj39SWIiAy7TsYSyjVFbHdMRrRcM+D0BeegrNTCtRIye4SHPnXEVE2RZpC0g+tV7L2cSRDhukgoeDLL79sNq+B+/v79u1rclVsieZSV3p7j1Ndymjr4KOH4LQ7cXO3m5Ghywh1k+SL0QhHVi+otR+QloDo94OOdUF05+jRcGRlNaJnT45kyDCtJaClHUUsbM7Xf+XKlZoM9qBVbHYb3jr6Fv52+G8os5bxQOxro9tNm6JGjn9wHMU7ioU9wbjfjAtJG8rLdZgzpydGj+6H2lrWWDGBuzcyLgK6e2jvntz/vEGRDN1tAdzVyow62V28G9XWasSYYjC5z2QWCrxg2vApwt54GfrTJ6HWLIhbnt4iyiMfGYno9NAY+H37rStoUWqqDfHxrjKjbprLUxDIZ1jz3IbtA0qP7Eso2Llz51XhjXmCUDck9H2w5wNRzu2UC7PR5QrGuFFXC0P+Uejsdjgj2yeXQqjJezkP1eeqEZMRg5z/yQldRM2tFKwoGgMGUCTDjouNoEQstb5zQugNehjDjX7V1el1MEWYfNY1R/k/Jkg5BSQoSi29RvMNeZeR7RnFpCE3dprEpVwHVI8yIkpu6VLuAm/nI40Afea1114TORPIJo3s16ge2a5ReuUFPqLjStoEKZuvrxwJmjQ0JG2BZ34Des4eANqCBuIt57coOrRxe2M88q0QCBzJaXAmp0GNfeDy4cuiPP6Z8TBGGEPWjgMHXBPQoEEW7ostsDja+8KO6DWjF+5ce2fT82XJy2CtcxlwetJtYjfM3zS/6fmfu/8ZdaV1Tc+fdj7tV3ukJEVSZkOa/Ck6LhmQ06RNEzmt8Ol1aTKmyZlc10kgkAzNc3JyRH2Kc+PtfFSfDhICaL6iyZ3q07koRL+vPDoFBQXCsF3KvUDfQ9+pVlf4gO9ib9oCTy0Bo37qrfXYe9kVoW9GXw5t/D3IRe7AHlW7IdI1n/nmTAx9YCjSR6WHtC2HD7uCTU2eHMt9UWHQKp00AJKrOZVp1U4TOE28tC1NCfe8JURyj35LEzsJA/QZX+cjpPf9Taa3evVqEUpZwls71ESrRHu6SK+//jruvfde8ehuS8Bog22F22BxWJAckYzeib1D3RzZoT9/Fvryy3CaTLD3HQg102V0l5B+f0WFDufOuTQFOTksELTEoppFzW4fuPPYxcea3T5w58HTD7aqPaQFoBW/NGl7Bq+TVuQ0yTeXrpgmetIStHQ+XxFzPQPjuaOl3AqtMtMlqU1SzXzxxRfiOaMdSF27r2AfooxRGN15NAwGNuzyRNIS2PsMBMxhqrv+O/+0E7UXayEHJCPDjAwLOnXivtgStNfv63C3J2iprrs9gbe6/kKLSmk+IWhSd09+R8eGDRtEPW8RcCXI7oC2Gpo7nzekzLxSrBxPSKhwD5hH9do9o2oIafUmIBlkXHfddWKvxRvsfaBe6NrOSJ+Ba2ddi3onGXbx6uwqGuphOHFYtVsH5IK47Zlt2Ld8H+7Zf0/IwhpLfTE6+gIefbSaQkFAp1NncCg1Q2p8msiXLl3atEqnPXtS29P8QlsCtFKngyZ9MjSUVP9kMEhIhoeSVsDb+SRjQTofnVd6j+rSaz169PBqQEjfRTYI9BnSQJAA4csgUQ3onG2YvZ988kk8//zzXt/74IMPAs6a2N5QjAVK+VxZWYnYWE6r2lrsdjuOHHFFr0tNTRUR5JROY2MjnnrqKVH+wx/+0LZEWk6na/vg9EnYxk6hmQpqgYzO3hrxlvA4GPOrMRj1+KiQtofyVJw7d06UqR/26tVLc0Kqr3GtoaEBp06dQlZWFsLDw6E2JK1Ac1sKTOB9oU0ivi+BgJCbQMAEB5Ihj5081vRcLVsHJAS88MILwTmZTgdHRjdxqI3df93tckHMDJ0Loi8oLbrWBAKGCTYc+osJWCh4duezmPH5DHxR9AXnO9AQNUU12PXirpC7ILpTXQ385z9xOH3arBoBlWkZaSuA7AhaE8SI8U3o72pGcey4uAOFtYXoFNdJNSsz4XN/2eVzT0FNWvu7TJv/CzjssA27Bs74BKiJrc9sha3OhvTR6eh9qzw8TvbuNeKJJzLRpYsFp0+zgKoVaMsgPz8/1M1QJSwUMAFxtvIsTtechkFnwPjM8VALdXV1wtCIKCoqQlRUKyIQNjbCsH83dDYr7L0HqkooEFkQ33PZkUxcPFE2wuC+fS7tQP/+9dDrOZIhw7QVFgqYgNhQsEE8DkgYgPgI7fju+oPh6AEhEDgSO8GRngk1EZ8Vj2E/HQZrrRWpw0OTBdGbdmfHDoq2FyGEAoCFAoZpKywUMAENwh8d/EiURyePls1qMRiQZoCst9uC8aAUwXC4qjwOiPDEcEx6fpKsXI2pLQcPuiyphwyxq6o/Mkyo8HsTjvIbMNqG3L++Lv5alDnfwdXoLhZDf7EYToMB9n5DoBYcNsdVgoCcrjnJcGfPulxHr7++s6zaxjCqFwooPjWjbQ5dOoTLjZcRYYjAtL7TeBB2wyhFMOzRF4iIhFrY+cJOfDDrA1w6eAlyY/9+l6IzPd2Czp3ZyJBhgoHfdxIlgfBHW8AaBfViNpgxL3seZnadiQiTKwGNmoJ7/PjHPxYHlQPCahX2BISdtg5UQvX5aux6aRcKtxSi7Nh3aWjlwtVGhiwUMB3D/fffLzIlqhW/7yTyB6UcB81N+mvWrBHhIhn1QSpkY6URvxr2Kzw9/GnVaQkoSuNHH30kDioH9mEbbENyYU/LhCMzC2pBji6I7v1xypTzeOmlM7jzTpcrKcO0JeaBvyxZsgRqxm+hgCIUUtrk++67T4Q3ljh9+jR++tOfCt9uii+t5kQRWoYGYWkFbTQaVScUtInwCNjGTYPljntUY2BYnFeMo+8flZ0Lont/jItrwNSp1Rg3jo0MtUhrghZ5+4zaJ/l29z6gGNuUHIKSIVFiCJKwKJ4yCQqUJGLFihUBN4KRP6fKT2Fv6V4MTByI9OR0HoRVDE24m5/cLMr97+ovGxdEX3Tt2pX7YwDUNpPckoJCuofGb64u7dhERPiuG0ioD0p+RPOJBM0l9NrixYtF4iLSVE+fPl0kRaJJXEqARPUoSZG0GJUSGnk7H2kE6DOvvfaaSKREC1lKdkT1yGYuLy/PZ6IjSZtAbVH7wtdvoYC0AVI2KZr4aeCgP/TVV18V2gOJxx9/vL3ayoQIutZ/3vxnvHzgZczMnIl/ZP4j1E2SDYbD38IZEQlHtx6uUVIFHFt9DMW7imGKMmHsr8dCjhw6pMe//tUZOTm16NuXwxsHQnS07/dmzADWrv3ueXIyBfbyXnfiRGDTpu+ed+8OlJZ+99xf71Upc6GU7pgm/7lz5wrNM80xNJHTCp9elyZ9mpwpqiEJBFSPoH1+qk+pkL2dj+rTQUIALWxpcqf6dK558+Y1ZV70hOpRhkSyq5Oeq9nw3u9RjP5k+tPpjyZBoLy8XEhWxL59+9qzjYwMhIJtRdtEeWTKSDbqkrDZRFjjsH+/C/0Z9YRcPfyuK+3zyEdGIjqtmRkkhGzcaMLf/paCd95J4v6ocGiCJQ0ArcbpoLI0t9CETYtRsmfzltaY6krQxE7zVHPnk85Jr/kSAjwhrYN7JkYSOtSM35oC+gNJ7UJ2Be6QloDSJJP0dOutt7ZHG5kQU91Yjf1l+0V5Zr+ZrKq9giH/GHQN9XBGx7g0BSrh5pU34/A7h9FvXj/IVUjdvl2KZEh2LhzJMBBqany/55lT6uJF33U9ZbHTp1vXHtIC0Ipfmng9UyHTJE7QBN9cmmSa6GnCbul8zU3qBQUFTd+nVfwWsUla8xQI3I0Qya7gj3/8I+64445gto+RAZvPbIbNaUNGVAay4tVjXd9WDFIEwwFDVbN1QBhMBgyaP0gWWRB9CQWHDrmCFg0ebGUhNUBor9/X4W5P0FJdd3sCb3UDmVskVT9Bk7q0b09lOjZs2CDqeRoKuj8nuwPaamjufN4gIYG2FKhOhRdDRNqecPdO8FVPLfg9krnbDXiDBAY66MIw6oEG4A+//bAptDGral3oKsphOFsA2ja1D1BHbIKCzwtga7RB7lRV0aqUIxmqBdJC00S+dOlSoaonjQCt1qlMW9akAZAOmvTdJ2gyGKR69Fl6j7QCvs4nGQvSFoO7kEB1ad6iusO9bClIWxh0PmlLgj5P9dWIzhnkYOZk0LFy5UrIkaqqKuE9QTHuY2NjQ90cRUChjfv+uS9OVJ3Ai2NexN0j71blIFxbW4v09HS/syQat22AaedW2Lv1gOXWH0LpkGHhe9PeQ3yPePxwyw+FkaFc2bJFhxtvjEFamgVnz+qFi6yW8TWukQvxqVOnhBY33FMFoAIkrUBzWwrtRUVFxVX2DHInkL4Q9GUfWXUy6qGkpkQIBDrocEPfG1QpEASMwwHjIZdxrW2g9y01JUHrgk1PuszIKVCRnAUCautXX1FGRI5kyISOeAUJBIESdBHbl90Bo0w6RXTCO5PfwbGKY0iOTg51c2SBrroKzvBIIRw4KNeBwjm66ihK8kpk7YLoLhQcOeISWgYNYnsCrSJtBZDaPzc3V9WTtCyFAlJLkefBY489FrBai2IaBPo5Rj4D8NkzZzE4cbA4eAB24YyLR+OPHgBqqr9vrq0wrLVWbH16q+xdEN159tnzmD+/FIMGZXCf1Ci0ZZCfrx43YMUJBbRfRR4GDzzwgLAZmDJlSoufISGCJLm///3vLdYlAw4pUhQ9SsEoWoLqk9EIh6lsH7QU2thgMODmm29uKrcI/RcxyrdLyXs5DzVFNYjtGovh/58yDCbJhKB370Z07cpbBwwTsu0DMlCgoEU02T///PNISEgQvqCS2oYMLy5fviyiPtF7ZADSkscCQfUkq1GCyiQY+BNYgj6rdZ/S9iS/LB+/3f1bjEsdhx/l/kjVQgEZ37z99tst1tNdKIYzIQkwm6F0qs9VI+8lV1CX8c+Ol60Loi/8Et4YhgmIgEVtmui/+OILsS1AggKtJim6IWkTaGKn9yiilD/aBILO4249SiEp/cm0KLmZMO3Hf/P/iw9Of4B3899lgy7C4YD54/cQ/toL0JWch9Jx2B3IGJeBLmO6oNctvaAE3nvPhF/+sgu++ipa1UIqw4SKVi8NpC2FtiCFoPT2uj9bBz169OB9pXaChL1PDn3SFJ+AB2BAf+oE9DVVIteBs1MKlE5ctzjMXj0blmqLIq4v9clPP7Vj7doEdO+ujDYzjNJo9fLvyy+/bHPOA29RoaToUi1pF3xls3KnsbFRGDu6H4x/2Ow2fFPyjSiPTRur+gGY4hSQoEsHlb1h3J/3XQRDFfnGm2OUsRVCQsHBgxzJkGFkKRSQXYF7konWQJO/Zxxq0hw0F0KS3vN324DSbkoDPR2ZmZltaq+W2FOyB9XWasSYYjCp9yTND8AUwVB/+qQo2wflQOkuiOsfWo+6Uh/p72RKZSVw5gxHMmQYWQoFZBB47733fu/1ZcuW+X0Ob4kpWooURdES/Y1gRYGUyJ1SOgoLC/1um9ZZX+DawhnReQTMRmWsJNtCZGSk2Iqig8qeGA7sBk1BFMHQGZ+oaBfELb/ZggNvHMChfx6CkjhwwKWdoUiGKSlsZMgEN+YB46LVOlAaPHv16iUmcAoeIan3KHGFv3EJ6LOeWwXetAcSLWXJ8iQsLEwcTGDQdVx72JVUfVTnUZpYkdFv7NSpk/c3bTYYD+11FQe7+rpSyfvLdy6Iwx5QVqCxfftcgsCAAfXQ6eQfT4GRX6hhb58hl3Z2aw+CUEATtLf81oFsKdAE7y3rVXMTv3sSCvJykOIckI0BR7UKDnaHHeUN5aI8Pn28JoSC5tAXFwL1dXBGxcCR3RuKdkH88xUXxOfGwxiuHLsIz3TJOh2nS24tvmxmmoMWV1KOCZvNJuy1yCMpwi1Voud5W8of4jmuuy8QaTyn12gLmLaLKXIheabRGE8TOBmak7aa6pFrvJTgSJqTvJ2P5iz6DLnVUyKlpKQk4fpO9WguoblrgQ9bNW/nUyutHhXowkydOvV7rwfqJkiBitw1AJTy0l1qI7UOXXCqR3XcBQaKi0DveRNOmNZDeQ5WTVuFkroSDOsxTBNCAQ1yTz31lCj/4Q9/uErD5MjMQuM9D0JXWaboFMlbf7sVtnqbywXxZmW4ILoLBVVVeuh0TjYybCPR0YFrWWjbliZU4sMPPxRB7CZOnIhNm1w5M4ju3bujtLS06bm/ufakzIVSumOa/On8NObTpE3fS4tFel0a62meobmABAIp2B1lVKT6pGn2dj5p/qDJnbaWae6g+nSuefPm+YyNs8dH+9S6CA1YKFizZo149LWa9yYoNAdJgE888URTjmoKSOR+cUiCo4viGeVQSosplVkwCA50I1M2LSI1MlUzAWJo9fP666+L8rPPPvu9bScKbUyHUinaWSQMDMkwYuLiiYqcVF955Qzq6vTo27eHItvPeIfGeJpgJVd0KtOqneYYKW0xzSsUGM8T94mZJnaavOkzvs5HSO/7EyCvpfZpWiggQz2SxCQ1DcUIoIsUjBTEze3n+HqPhAAWBNpHKKiuq4ZJb1J9aGO/sDQCZuXbpexYskM8DrhrAFKGKjfGQmSkA5GR2hBU24uampqAP+MuJM+ePVucwzOg2enTp1vVHloM0opfmmQ9J1tJ+9ySTRlN1qQlaOl8vmzWiIIrofYDaZ/a8FsX+uSTTwqJjQwMKYLh448/7lcYY0ZZFJQXYOzHY/GLr3+BxKREbQsFVgvC//fPMH/0L2FToGSuf+16kdtgzK/HQImC6oULF5qec3TNtkF7/YEekj0BQWV6zd2ewNt5/YXmFUk1L03C0uKTynSQATvV82aD5q51pq2G5s7nDSk2jqSt9iTQ82lGU0A3JgkC7oYWav5jtBzauN5ejyprFUxGV4parWI4dgi6hnroLl8Cwq8eAJVGRGIEJv5hIpQIjT0PP5yCkhITHn64DAMGaFhQVSGkxqeJV0qMR9B2MW0P09YybQmQFoAOmvRJeyyp/slgkJAMD6VVvLfzSa6HdD46r/Qe1aXXSPu90Iv22Vf7oHWhwJtRBf2JjHqgwfffB/4tyuNSxmlbS+B0wvjtzu+CFSn0vyjPL0d8dryiryXZq+3aFYXSUhPS03lLS414m2TpNffXvdkUkNGfN3W+t/P5SrdMk35L+XZuV7EQ4InfejhvPtzebs5Aghcx8qLR1oivS74WZa27IuqLCqG/WAKnwQjbQGX580tUFVbhH2P+gTW3rEFjZSOUyvnzOiEQGAxO5ORot08yjKw0BWSBSSmR3d1M3PdZpL0Wqudv8CJGXmw7uw11tjokhiViYh9lWqgHC8O+K1qCfoOAiO9HOFSKC6K9wQ671Q5zrHKjUu7Z4xqmevVqQHS0tre0GBfSVgDZEVDwPLW6B8paKCDVjbtNgS/BgLwUGOVBwt6qPatEeUzKGBj02rXw1tdUw3DisCjbho6CEinaUYRjq4+5XBD/oFwBj/rlV1+RliMKgwbVQ69XvicI03Z8bQUwHSgUkAEGJUHyx0uBUebgu/n8ZlEenzZe0xbe4Uf3Q+d0wp7RDc7OynPfczqc2LxosypcEKlffvutSxAYNozTJTOMbIQCCgwRzHqM/AbfGZkzEGuOxY39b9T04Fs3ZASM8QlwxidBiVCQopLdJTBFmzD2N2OhZOx24NChcFG+7roETfdLhpGVUDBs2LCg1mPkx9197hZHWkwaNI3JDPuQEVBqFkSyJSBGPjoSUSn++4vLkcpKHUaNqsWpU2Ho148FAoZpb5STEYXpkNDGBK/GlEv1+WqYokyI7RaL4T/zL4yrnElMdOLll10pz83mfqFuDsOoHhYKGOGKuPLYSoxLHYfkqGRNCgVkQ/HAzOvxYN/uCD95BBiizBTJib0T8aNvfoSqs1WKyoLoTyRDLfbL9sDhcPidrKgt0PXSsm2SUlH2qMEEhc2nN+OXeb9EcngyDt97WJODL4VsffEHt8Fw8ghsF4tASXqVisFkQEKPBCgdmriKi51ISqLY+2ZN9sv2EAiqqqpgJ2ONdoaSqVFuHDkIBuTCKEUulFvugvVX2kbulf4maWpPWCjQODTwvrv7XVEmTYEcbuBQoKsog540BOSGOGw0lEbxrmLhhjh0wVAYzOpwJ6W8PVOn9kFiog0HDzpYKAjS/U4CQXuv4kn4oO/pCI2EP0hpk+XItGnTZCEMSLBQoHHo5t1QuEGUJ3eZrFmhwLhnO7n0w969J5ydkqE0F8SNCzfiwp4LqC+tx7jfjoMa2L3bCIdDB5OJtAXqEHTkAt3n7X2vB0MbQQHxQhGYqCII3xuqtrcVFgo0zpHSIzhfdx5mvVm7roj1dTAc3CuK8178G1699paAsryFmiMrjwiBwBxjxrCfqcf7Z8cOlyAwdGgd9Pq2p2hn5AslP6JMhe4J9yhiISU5ouRFlMmQAuXl5OSIMmXopTqUHMlT/e7tXN6gz1EiJTofPUrJkNavXy++lz5HKZOlxH/0fiDn9nUOX1A0YKonnd/93N7aSm1ZvHixSNJEv53yQJAQQv+JZ/rnQGChQON8evxT8Tii8wjERcRBixi/3QWd3Ya9xRfw8beH8CoU7IKYrBxhpsVgWpsdojx8OKWtZqFArUiZC6XouDS5zZ07V0zwUuRCmiCpTEIBPafsiFJCPgpzTJ+n+r7O5blip4mV6tFkKrXB3eZg2rRpYjKXEiHR90ZHR+PDDz9s8dyEr3NIv8kb7nWpLSQIUN3m2krnI+GDtBLUnuaEDn/Rpq6YaRp4PzjwgShPTFNuKNw2YbPBeCXPQfx1N2P79u3fyxMvZ3a9tAu1xbWI6x6nKi2BzUaRDF1Bi0aOtGqzb2oEWiHTxEqTIB1UzsvLu2qyJOh1aRVNkyKtiAlp8vTnXBI0oUrnJWhyXbFixVV13Cd7Csr35z//2a9zN3eO5rIxutclAUISApprK2kE6L+YOnVqUAQCgjUFGqayoRL7SveJ8pSMKdoceA0GWG64FYbD+5E6aRpSDQZFZUHM+4trUBr/3HgYw9RzOx8+bEBtrQHR0XZMn56mzb6pEWiVS5OeNLH7YxBIkyFNorRqdlfnt+Zc/hAfHy9sMNpybjoHqf9pZS9tJ5Cg0NbJXNoqcBeO2gJrCjRMjDkG62aswwujXsCwrGHaHHh1Oji69YD1htlCQFAS257ZJrIgZozLQM9ZPaEmtm93XYvBg+tgVok3BeMdWum6J9ajiV2aNFv6HNkWuE+E/p6LJmPagpCgCVVS3bt/VoJW7T//+c8DbqfnOUggoEdKMEiHu0DgXpe+h+q21Fb6DB0bNmwQv939HK1FPUsLplVRDDuFd8K1Gddq0+uA3KWuCEIWiwUvvPCCKD/66KMwm+WfanjEoyPQUNGAsb8eqyqBjvpmt26X8KMfRYl0yTqdOuwk5OZ1JJfzkzqcJrSlS5c2rXpp0iMtAKnnSRNAEz+pzGnSIwM7ep/sDGiCdTeq8+dcZINA9WiylerRpCup6yWkiZgmfsloMDk5+Xvnbg7Pc/hayZMGgZAMGck+QKrrq61Ul4wwaXuBtBB00HfQtkpbXBx1Trk4knYAFLQjLi5OpHemoBpahm7aw4dd6YGNRiPS0jSmonU6Yf7gH3CkdoEtdwxq7Q6kp6eLt4qKihTlfaDGvnnu3DlRJuGsV69e2uqbQRrXGhoahOCflZWF8PBwTQcvCpQ5c+YI4aIt6vhgnCNYeOsLvmBNgUb56OhH+P3m3+O2rNtw/+j7NTfo6s/kw1B4CvriQtiGKytYEXkcUH4DLdCtWzfN9c32hCZomqg5zDHjC75iGoQGhDd2vYHdpbtxtOKo9m5cpxPGHV+Jon1wLhCpHK2ApcaCt0a+hQ2PbEBjVSPUyOHDeuzYEYW6Op1YbTLBhe53+l/b+1DquLJ+/Xqx3UAq+tbu0QfjHKGCNQUapMHagI3nNory9V2v19xKTH/uDAxFhXAaDLDmjIGSyHspD9XnqnFmwxlMXDwRahRYX37Ziffey8IPf1iKt95S5sTCKJdpV2IjhPocoYLvOA2y/tR61NhqRAKkaX2naU4oaNISDBwGRMdAKVDmw7yXXS6IE343QVUuiO5CwY4drj3Pa65p1FzfZJhQw0KBxqBB9+1db4vylPQpMBrUN7E0h76oUNgSOPV62HKVlSNgy9NbXC6I4zPQ40ZXNDe1UVKiw6lT4dDpnLj11k4sFCjA04CRP4HYkGhrRmBgtVuxvnC9KF/X9TrF7vu1FuNOV0hge/8hcMYqJ6zz+e3ncXzNcVDWpkmLJ6l2sty61TUk9e3bgE6dtGFM2V6Q5wbd3+RN07lzZ/Fcrf2GaV4guHTpkrj2JlPL9xQLBRrjy1NfosJSgXhzPG7of4PmBgnr5OvhjI4RbohKyoK46YlNojzwxwPReVBnqBVJKMjNrYVe7/LdZloHCQTkglZcXCwEA0a76HQ6ZGRk+GW4y0KBxiTG6tJqjE8dj4yoDJiN8g/QE2yccQmwTrsRSqLsRBkqT1WKLIgUqEjN/XPTJpeQOnJkLYV0CXWTFA9pB7p27QqbzdYhsQkYeUIaAn89eVgo0BA06PaL6YdXxr6iaJehVmGzAkZlqqOT+iTh7r134+L+i4jsHAm1QvGKzpwhlbcTY8bYNKfFai8ktbE/qmOGkYVQIIVvpFCQ9Nhc6Ejy+ZQyRFF8aIoY1VKoSeb7xiYpKSnaGXSdToStfBPOmDhYJ16nKFsCiYikCHSb3A1qJi3NidWrT+L48TAMHsxBixhGk0IBTeruMaGpTIKBr9jNUqxngmJRJyQkiMQSbYn1rBWB4O+b/44+0X2QFpmmKS2B/uRR6C8UwVl2CU7jTCjJBbH0cCmyrstS/QTpMoa6gD59rOjTpwF6fZdQN4lhNEnIZwZa9bvHhqZEEL5yTpMmgQ4pQhQlgKDPLl68uMPaq1TK68vxyLZHcO1n1+J8/XnVTzJNOBwwff2lKIpwxj6iF9L/0bdvX3HI5b8hF8SP5n2Ebc9ug9ohocBqtYoyW8kzjEaFAgoFSRO7t9d9IWW6kpC2HZjmWXNkDawOK3rG9kRudq5mBl3DkW+hLyuFMywctuG+PQ4iIyOxY8cOcVA51Jz/5jsXxD639oHaOXZMjyeeyMDnn8eie/fumumfDCM3Qrp94C0mNKWQdJ/03SEBoLy8/HsChK8sVI2NjeKQoOxgWoRWYf+36/9EeUbmDO3Ek7dYYNp6RUswcjzQQnYwWbkgPqkNF0SJL74w4T//iUBVlQEPP6yR/skwMiSkmgKa/KU80hKkOfA3gQQlnKBzUP5ob9C2AqUUlY7MzExokZOXT+KbC99ABx1uyrpJM6sw466t0NXVwBGXANvQkVAKh/91GBf3XVS9C6K70PrZZ66oe2PH1mimfzKMHAmpUOApEBAkEHjbUvDGfffdhw0bNvisv2jRIpFjXDoKCwuhNWjAfWnzS6J8TfI1GJ49XBuDrt0Ow4nDomidMB0wNq8Uq6urw6hRo8RB5VBmQdz6jCvq4qjHR6naBVGipsaJvLwIUZ48mfMdMIxmtw9oMvfcKvCmPfDlhfDaa68163UQFhYmDi1js9uwOn+1KN/e43bteB0YDGi8634Yjh+Eo0dfv4Sno0ePNpVDxa4/7ULdhTrEZcVh6ANDoQUoiqHVqkeXLhZMntyFhQKG0apQQLYAnlsF9NyXjYC7x8K8efOaBILm7Aq0zqFLh1BtqRZhjecNnaetAddkgn3AML+qhoeH49NPP20qh4r00emI7xGP8c+OV2UWRE9IAPvkE5fXwZgxNTAa/dMSMgzTPoR81KHAQ+6TOgUkcrcRILsB8i6QAhRJngmkZZDcE6kOCwXeB9zo2mhsnLkR+dX5iDC5VLSq3zY4dhD2voMo+LvfHyPjy/HjxyPUZF2bha6Tu0Jv1IZGx+FwYuNGl5vopEn10OkSQt0khtE0IRcKVq1aJbYCpAmeghm5bwm8//77YtInoYDepzgGnvgyNNQ6lDK1oaEBUaYoDE8eromtA+PeHTBtWQf7sYOwzL4LShLgJC2OwaQd6/uyMh3i4+0oL3fgrrs6a0uTxTAyROcM5QZqB0MuieSFQEaHsbGxUDN0Wbcd2IY4XZwYaNPT02FswdhO6eiqKhD21ivQ2aywXHuT31sHBAXOefPNN0V5/vz5HRonnlwQ18xeIyIXDrlviGaEAuqjJSUl4r+vrDRg9Oje2nGXDSJaGteY9kf9S0eN0mhtxK2f3Yq5G+aiuKFYE4OtaePnQiCwd+kGe//AjPQsFgsee+wxcVC5o10Qz246i28Wf4PGiu/iamgpimHnzhpL0MUwMkXdS0eNQoPtq1+9iksNl+CEEwO7DlS9WpbyGxgKjsGp18M6dSbFLYYSsFRrzwVRgmKJ1dToER3tEOl91d5HGUYJsGiuUqHgzSMuVfgdPe5AhFnlBob1dTBvcHkO2HLGwJmknAiAu150c0G8XxsuiE0Juv5uwfjxffHSSymq39piGKXAQoEK2XV+F74t+xZGnREPjHhA9WpZ06bPoaurhSOxE2yjJ0IpVJ6pxO6/7hblCb+boAkXxKujGEbDZtMjJcWpie0thlEC2hmFNDTYPrfhOVG+LuM6pMemQ+3YRoyFvvwyLFNmthi5UE5s+c0W2BvtyJyQiR4ze0BLnD6tx7ffRkKnc+Lee2N564BhZIJyRlDGL46XHsdnZz4T5fv636d6LQHh7JSCxh/cqxg7AqLiVAVOfnISOr0OExdP1NykuHq1y7tj5MhaZGYqI1EVw2gBFgpUxj/2/wMOODAxdSIm9Jmg3snG4YDu0gU4U9JczxX2O+Oz4vHDrT/Eua3n0HmgcmwgghWw6L33XMLqzJkV0OvVr81iGKXAQoHKtg5+3PXHSL0mFSmRKeoVCKjj7vhKHJTsyD78GiiRTv07iUNrHDigQ35+GMxmB2bO5ARIDCMnWChQWQRDS6MFk9InieA7ah1s9WcLYNy+GeLXRbhC5CoFckGsPl+NpL5J0CorV5rF4/jx1Rg8uJtq+ynDKBEWClRCSXUJis4WwQzXgJuSolJNQU01zJ+tEQKBbeAw2PsNhpIgF8RdL+3CmF+NwchHRkKL2qy77ipEZGQ0+vZtgF4fE+omMQzjhvqt0DQy0C5YswBTP56KDec3qFdLYLPBvHaVy/2wUzKsk2+Akqg87XJBdNqdSOqTpN0kXdGN+NGPLmPcOJs6+ynDKBjWFKiAPUV78MnpT0Q5IyYDqamp6htsnU6Y1n8KQ1EhnGFhsMycAxg7Lj9BMNjy9HcuiNkzsqFFgeDChQtNz7Ozs9XXTxlG4bCmQAV2BA9++qAo35BxA6b0m6LKgVZ/Jh/GI9/CqdMJgcCZGHwDvaSkJHG0B+e2ncOJf5/QrAsicfSoDnfckYFPPolDWFiYJtxlGUZpsKZA4aw5sgbbSrbBpDfhoSEPqXagdXTrAQttF+h0ohxsoqKiRPru9sBhd2Dzk5tFeeD/G6g5F0RJS7BsmRP79sUiIcGOxx9P06RgxDByh4UCBVNrqcWD/3FpCeb3no8RPUaob6ClzN70m3Q62Icq0zDv8LuHcXH/RZhjzcLAUIsUFwMffeRK63vffRUwGNjAkGHkiDqXlRpZeT360aMoqitCWmQafjrgp6rTEujPn4X5w3eAxgYoGafDibC4MIxeOBqRnbSTBdG9r774og1Wqx5DhtTh9ttZS8AwcoU1BQrFbrej3lovyk8NewrZmeoy2tIXFcL873ehszSKmAS2ide16/fV19fj9ttvF+XVq1cjIiJ4mSUH/b9BIrdBWGwYtEhlpRP/+le8KC9YUA6jkSMYMoxcYaFAoSsv2v9+fPDjuLX7rRjfd7y6BIIzBTB/8h50VivsXbrBNmZKhxhsbt26takcbLSoIZD66pIlNlRXG5CV1YCf/CRZVX2VYdQGCwUKHGSPHjsKu80unvdJ7KOqtLP6k0dh/s9q6Ox22Mm4cNZcwNT+rodkDf/mm282lYOVBTFjfAaypmdBq1y+DLzxRoIoP/xwGczmK7kqGIaRJTonzTIaoaqqCnFxcaisrERsrMvoSUnQpXr5y5fx+qHX8afRf0KXmC5IS0tTjS2B4eAeEYtA53TC3rMfLDfcqqhUyO5QoqNVM1cJF8T5e+aLBEhag/prcXEJNm4Mw3//G4eVK6NhNisrtoQSUPq4xsgLdcwmGuHAhQN48psncaD8AP599t+qEghgscC4/SshENj6D4Fl5u2KFQiEC+IilwvioPmDNCkQSEKBzWbF+PE1WLr0EkwmZV5PhtESfJcqhLK6Mtz87s2ot9djVOdR+PXEX6tHICDMZlhuugOGU8dhGzm+w1Mh22w2fPKJKyrkrFmzYGyDQCK5IJLHwTW/VGYGx7Ziszlx8GApEhNdzzl6IcMoAxXNKuqlprEG096YhtPVp5EakYo/jfkTws3hUINBoeHQvqbnzuRU2EZN6HCBgGhsbMT8+fPFQeVWn6eqEdue3SbKoxaO0qSBIWkIXnihATNmZGHNmgRho6EmuxeGUTOsKZA55HY4460Z2Fu6FzGmGCyfsBwDsgYoe9VltcL4zUaYdn8Dp8EAR0o6nJ2SoQZ2/WkX6i7WIT47HkMXDIUWOXBAhxde6IzGRj2cTgN69Oih7P7KMBqChQKZr7juXXUvthRvQYQhAsvHL8eEvhMUPcDqT52A6cv/QF9VIZ7b+w+FM85lna50aopqsOdve0R5wu8nwGDW3uq4qsqJH/84XAgEEyZU4cknk9S1zcUwKoeFAhlD/vI/6fUT5F3Mw6Jhi3DDoBsUO8DqKsph2roehhOHxXNHdCysU2bA0aMP1EJUWhRmvjkTBZ8XIPsG7WVBtFqB+fMjceqUCcnJVmFcaDZr739gGCXDQoEMOXH5BLLjs3Hy5EkkRyRj1dRVyMzIVKxAAJsNYe+ugK6xQWQ5tA0fDdvoScK4UE2QBociF9KhNSje04MPhmPDBhMiIhx46aWzyM3trmitFsNoEYXOMurE7rBjydYl6Pe3flj232Ww0tKLgumYlZdmVldZ7kpmRBiNsA3JFcGIGu9aANuEa1UlEJALYkO5svMztJWPPzbgnXfM0Oud+OMfCzFihFNxfZZhGNYUyIa9xXvxwNoHsPP8TvE8rzQPM7vOhMlkQmpqqjJWXHY79KdPwrg/Tzxabv1hU5pj2zWTAZVOEoffOYyvfv0VJjw3AQN/PBBas3uhY8SIIsyZk4xhw2oxfXojevTopYw+yzDMVbBQEGJOXj6Jpzc9jX8d/BeccCLKGIXHBj+GeT3nieBENLDKenB1OKAvPgfDsYMwHD8EXX1d01v0uiQUqFUgEC6Iz21DY0UjLNUWaImdO/WIji5BdLRLS/Kb3xTBbDajVy8WCBhGqchCKFi6dKkIbkJJfuhRylYXrPpyhFZXS7ctxaINi4QwQFyfcb1IctQltosytAP1dQh/86/QNbiyNRLOyCjY+w6CbXAunAlJUDs7X9gpXBATeiZgyH1DoAVKSoDnngsT2wVTpnTGiy8WitAS4eHh7H7IMAon5ELB/fffjzlz5mDatGniOZVpoh8+fHhQ6suFRlsjvin8BilRKeid1FsINKm2VCEQjEsZh/8Z8D8YmjwUKSkp8tIONNRDf/kS9CXnoSs+J+wDrNfPdr0XEQlnRKSwHbBn9xbCgKNrtmq1Ap5UnKrA3lf2asIFkcxD9u/XYflyM1auNMFqdV3jqCjKKGlGv349hA2BbPotwzDKTIhEg4h7E1asWIHdu3dj+fLlQanf0YlDqG211lqcrzqPvKI87Di/Q9gJ7C3ZC4vdgh/0+AGeGvpUU91zteeQnZAdGmGA/kcyZrRagKjoppeNW9bDcP4MdOWXr9ICiI+YzWj46RNNE7+uqgLO6FjFCwK1tbVIT08X5aKiIkRFRbX4mU9+9AlOfnwSXSd3xa0f3qq6CVGyF/jf/zXjL38Jw9mz3wk9Q4bU4bHHijF6tJO1AyGGEyIxqtEUrF+/HvHx8V5fD0b9YMUKsNltOFVxCpfqLqG0thQX6y6itK5UPL9YexEDOg/AwjELXQlg7DbELolt2hJwJzEsEeGG78ITkxHhqN6jmgZUaRD2m9oa6BrrobNYhLsfLI2iDCpTW3K+i7tv3rJO7PHrLI0U09f1SPXpe40m1P1/TzbV1V0qEXWb/oOoGDhS0uBIzYA9rYv4T5oggcD1R0HJ0P8eERHRVL7qN3rh3LZzQiCgLIjjfzc+8GsXQqidFGTo1CkDioro0Ivj/Hl61InHv/2tFtnZJcID5vLlRJw9mw6z2YFJk6rxox+VCmEgKyuLtQMMozJCKhRUVLii2rmTmJiIsrKyoNSnGPbucexJog4UmhyOHD2CYWuGeZ3oCUpQdFPiTU3P48xxaLA3oFdsLwxOHIxBiYPEY0ZUxlUDKCXhoVVpa8nctBaRly96b7fBiNMpmU3PuxSdQ3TJdxP9VdhtOHf2bNNqPzKzB/SpmbBGx8ISHSuEhqtoQ5vlzM6dLs8P6k+++pTEsS+PCYferDlZaIhtwLlzPv5bmbJqVQKefbaLz/f37atEZqbLJXb69Cp06WLFiBE1SEwMY2GAYVRMSIUCGnhpUneHNAHeJv/W1F+8eDGeeeaZNrfToDMgNTIVOuiQEJaABHMCEsMTkWhOFM+zYrKuqr9uxrqrNALthd0c5jqMJjhMdJjhMH73KLYHrgzcZb0HobJ7b1c9owl2qiPKZpF/wD0JUV2K78mCcdHn7j5IHZOK8GRlJqZKS7MiKcmK1FQrUlJsSEmRyq7Hnj1dwjQZD06YkIWJE139Q1b2LgzDqEso8JzgCZrgvW0RtKb+okWL8Mgjj1ylKcjM/G717A+0Iurbty9O9z7dpB72d1Bs9wG0b99mL2KCl7pMEFHoX0r9uHdvJ+65hwIMUcAh34INCwEMoy1CKhTQZO6ppvWmDWhtfUrZSkdbIKGAI7MxDMMwWiCksx25FXqq/um55G7Y1voMwzAMw/hPyJfAFHjI3Xtg3bp1IhaBxJ49e7B69Wq/6zMMwzAMo9A4BcQTTzwhfJ1p1e8ZoZDeI8GAJn9/6jcH+/MyDKM2eFxjVCcUdBR88zAMozZ4XGNUtX3AMAzDMIw8YKGAYRiGYRgBCwUMwzAMw8gjS2JHIplPtCbcMcMwjByRxjMNmYcx7YimhILq6mrxGGhUQ4ZhGCWMb2RwyDBtQVPeB5TciBIQxcTEtHvoVimkcmFhoeosgtX629T6u9T829T6uwL5bTSEk0BAqb85+irTVjSlKaAbJiMjo0O/k25mtQ1Wav9tav1dav5tav1d/v421hAwwYLFSoZhGIZhBCwUMAzDMAwjYKGgnaDsjE8//XSbszTKEbX+NrX+LjX/NrX+LrX/Nka+aMrQkGEYhmEY37CmgGEYhmEYAQsFDMMwDMNozyWxIykoKMDy5ctFiucFCxZ87/2lS5eKtM9UL5D0z3KD2k8prIcPHy4e8/LykJubi/j4eCgJtVwPtV4fNd9Tzf0uNV0/RhmwpqAdWL9+Pfbs2SMOupE9uf/++8VNToPWwoUL8f7774u6SmT16tXIyckRwaASEhLE71DagKWm66HG66Pme6ql36WW68coBxYK2oFp06aJwcnXzbtixQpRR2L69OlipaBU8vPzsXv3bhFZjQZkpaG266G266Pme6ql36WW68coB94+CMHKwNsAQK8rFVLVKhU1Xg81XR9/UPs1VPv1Y+QFCwUdjDcVYWJiIsrKyqBUSMVJ7Nq1CyNGjFDMXq5ar4earo8/qP0aqv36MfKChYIOhgYqGrDcoVWOt4FNCdA+rmT4RIMVGUvRyoZeVwJqux5quz5av4ZauH6MvGCbgg7Gc/AiaPBSqvEQ7Ym6t50GKyXs5ar1eqjt+mj9Gmrh+jHygjUFfkBGTGTo0xxz5sy5ytDJF3SDe6o1va10lPJbqb67GxWtYpS0lyv369FWlH59tH4NtXD9GHnBQoEfePOJbi00mXqqNem5PwKF3H4r+VCTKxi1XTKGot+iJMMouV+PtqCG66Pla6iV68fIC94+CAG0N+gu7a9bt07c/EqDBqclS5ZcNUjR71q0aBGUhFquh1qvj1avoZauHyMfOCFSOwYkWbx4sbih582bJwYt95v7iSeeEEZDkuSvVItiKRpbUlKS8Kf2dxtFbqjleqj1+qj1nmrpd6nl+jHKgYUChmEYhmEEvH3AMAzDMIyAhQKGYRiGYQQsFDAMwzAMI2ChgGEYhmEYAQsFDMMwDMMIWChgGIZhGEbAQgHDMAzDMAIWChiGYRiGEbBQwDAMwzCMgIUChmEYhmEELBQwDMMwDCPg1MkME0QoGc/KlSubysOHDxdJbXbv3i0S2zAMw8gZFgoYJoisWLECCxcuFOWEhASsWrVKZLxjgYBhGCXAWRIZJohQGlzSDhA6nQ58ezEMoyTYpoBhgogkEKxfv57z3jMMozhYKGCYdmDdunWYPn36VRoEhmEYucNCAcMEidWrV2POnDlNZbIlkLQG8fHxIW4dwzBMy7ChIcMECRIC6Fi6dKnQFEjGhSQQSAICwzCMnGFDQ4ZhGIZhBLx9wDAMwzCMgIUChmEYhmEELBQwDMMwDCNgoYBhGIZhGAELBQzDMAzDCFgoYBiGYRhGwEIBwzAMwzACFgoYhmEYhhGwUMAwDMMwjICFAoZhGIZhBCwUMAzDMAwjYKGAYRiGYRgQ/z/LiLB/zjJCpwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(4, 3))\n", "a_d.plot(ax=ax, label='expert a', color='green', ls='--', zorder=100)\n", "b_d.plot(ax=ax, color='salmon', zorder=50, ls='--', label='expert b')\n", "c_d.plot(ax=ax, label='expert c', color='purple', ls='--', zorder=100)\n", "d_d.plot(ax=ax, label='expert d', color='blue', ls='--', zorder=100)\n", "ax.vlines(fixed_value, ymin=0, ymax=1, color='black', linestyle='dashdot', label='expert e')\n", "p_env.plot(ax=ax, bound_colors=['lightgray', 'lightgray'], fill_color='lightgray', label='envelope p-box')\n", "ax.set_title(\"envelope aggregation of expert elicitations\")\n", "ax.legend(loc=\"center left\", bbox_to_anchor=(1, 0.5))\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "c3b03297-bd53-45af-9904-5f0f75e8844a", "metadata": { "jp-MarkdownHeadingCollapsed": true }, "source": [ "## expert opinion is represented as intervals\n" ] }, { "cell_type": "code", "execution_count": 14, "id": "bae7e75a-6d8e-4003-a139-ffe5bcfab509", "metadata": {}, "outputs": [], "source": [ "a = pun.I(1,5)\n", "b = pun.I(3,6)\n", "c = pun.I(2,3)\n", "d = pun.I(2,9)" ] }, { "cell_type": "code", "execution_count": 15, "id": "2ead5c6f-9ff3-4389-8ca3-501fe377dad6", "metadata": {}, "outputs": [], "source": [ "dss = stochastic_mixture(a,b,c,d)" ] }, { "cell_type": "code", "execution_count": 16, "id": "b572838f-4ca4-4814-a2cc-d7e9f4971f7b", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAn8AAAElCAYAAACVs4T3AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAKU1JREFUeJzt3U1MXGee7/Gn7DjY2LwnHQ2KpTao76h71QGjaGayCiDvO2Bas7pXt43dLd27iBIcZjGtzKIJJFdqqTcBWprdSBiU/QSTZSQ3L8kmUUYdsDV4Io2SACYNNn6rq9+TPPShOFV1qiiqzjnP9yOVqZdTh1NFcfzj/7xlstls1gAAAMALJ2p9AAAAAKgewh8AAIBHCH8AAAAeIfwBAAB4hPAHAADgEcIfAACARwh/AAAAHiH8AQAAeOSZKBs9ffrUfPXVV6ahocFkMpnjPyoA3tF88999951pb283J06k7+9SzqMA4nIejRT+dMI6f/58JY8PAEKtr6+bF1980aQN51EAcTmPRgp/+kvV7ayxsbFyRwcAP9je3rbhyJ1v0obzKIC4nEcjhT/XRKETFictAMcprU2inEcBxOU8mr6ONQAAAMiL8AcAAOCRSM2+AJBGa2trZnJy0nR2dprh4eGi209MTJiOjg77PH0dGBioynECQCUR/gB46ebNm2Zra8usrKyYtra2ottfvXrVDA4Omr6+Pntb1xUAu7q6qnC0AFA5NPsC8JJCnCp3zc3NkbafmpraD37S399vq4YAkDRU/gAgQpUwLCTqfhzfpNjf7H5jr9efqk/tKHAgTDZrzO7u99efa6o3J05U9vNP+AOAItQ8nKu1tdVsbGzkfc7e3p69BOffQnQKfi/8vxdqfRhAzf33//2L+VHL2Yruk2ZfAChCIU9hL0iVwLBQ6IyNjZmmpqb9C6t7ACjHk6dPTKVR+QOAInKDnyj4FeovODo6al5//fVDM+8jGjX1Op/9z8/M2WcrW/kA4mx3x5if/ew5e73u/1Q+qhH+AFRkahOFIQ2KkPn5eTs6Nri99qFtNDpWX5eWlszFixcjD7ioJR1jbhNvWDUwqK6uzl5QnmAfv8YzjabxNKuiwB/P6p9H3//Bk8k8rPj+CX8AKjK1yfXr1/dHv2rOvJaWFrO8vLy//dzcnN3GGR8fPzB6Ns50nLlNvLqdlOMHgCD6/AE48tQmquq5yp6rlOm56vcWtLq6agNhNps1IyMjJs40/58Cq6MqZnB0r6tuAkDSUPkDUJGpTdSMq6ZQ9zxVCXVfkO6LC70WBTx9VXB1Ac8d48zMjH3cNV3Pzs7ayqULuQp+TPAMIIkIfwCOPLWJAtPm5uaB+xSqcptFXSVtcXHR9PT01HR5NB2bLvkqkGqWjnIfACQN4S8CNVFl79+v9WEAFZE5c6bohLnlTG0SpIqZ9hEMS6qSuQEeCn1aT5fl0QCg+gh/ESj4/UdXd60PA6iIv11ZNpn6v06jUampTYKuXLliFhYWDmyfWwVU6FMfQpZIA4DqYsAHgIpMbeKoX9z09PShip6bBsYJ6xMIADh+VP4iNpOpWgKk5fN8XFObKOANDQ3tBz/X70+DJDRAQtfdgArtL04DQADAF4S/CNQ/qlgzGZA2bmoTF/g0tUmwD5/69SnUuUEbbiSwqoZuRKy2cYFPzw2GPW2vEbQAgOoi/AEIVWxqk+BUKHpc8wDmCoZFbaf9tbW12fn+1NePwR4AUH2EPwB5FZraJPiYqn0aFV+Iq/4BAGqLAR8AAAAeIfwBAAB4hPAHAADgEcIfAACARwh/AAAAHiH8AQAAeITwBwAA4BHCHwAAgEcIfwAAAB4h/AEAAHiE8AcAAOARwh8AAIBHCH8AAAAeIfwBAAB4hPAHAADgEcIfAACARwh/AAAAHiH8AQAAeITwBwAA4BHCHwAAgEcIfwAAAB4h/AEAAHiE8AcAAOCRZ2p9AABQSxMTE6ajo8Osra3ZrwMDA0W3b25uNltbW5G2B4C4IfwB8NbVq1fN4OCg6evrs7d1XYGuq6srb/AbGRk5cHtlZSXv9gAQRzT7AvDW1NTUfvCT/v5+Mzk5mXf7xcXFA7f13Js3bx7rMQJApRH+AHhJoU3Nt2H356MqnwKjs7S0RNUPQOIQ/gB4SX32crW2tpqNjY28z7l+/fp+U7GCoPYRrBwG7e3tme3t7QMXAIgDwh8ALynkKewFuYEc+QwPD9vL3Nyc6e3tzRv8ZGxszDQ1Ne1fzp8/X9HjB4ByEf4AeCk3+ImCX1hTsKMmX/ULXF1dNRcvXjTd3d15m4lHR0fNvXv39i/r6+sVPX4AKBejfQF4SSEvt4k3rBroaCqY+fl5Mzs7a2/rukb7qil4eXn50PZ1dXX2AgBxQ+UPgJfUZJvbxFuoD5/6+PX09By4LzjtCwAkBeEPgLc0QXOw2VbVPA3oCAY+9e8ThUI9nlsNHBoaquIRA8DR0ewLwFtqwlWzrUKcqn4KfsGpW2ZmZmwAVEhUM/H4+LjdvrOzc38bqn8AkobwB8BrCnRRH1MwZF4/AElHsy8AAIBHCH8AAAAeodkXABAqm82a3Ue7Vftewa9/efiXqnxf5KcfxW51fvypkv3hM3yU7XZ3M+Y4Ef4AAKH/Mb3yr6+Yj9c/rvWhmEzmeP8jxGHKJZcu1Ztbt4gJafz80+wLADhEFb84BL+X/+Zlc+7Zc7U+DO+o4kfwq72///usaWo6VfH98pMFABT0xf/+wpyrO3eoClGNilz9M/VU/mrsP/9zzzQ1hceFQj+bfI/x84yuvj5jjuPtIvwBAApS8Guoa6j1YaBGGhpOmMbGk7U+DFQQzb4AAAAeIfwBAAB4hPAHAADgEcIfAACARwh/AAAAHiH8AQAAeITwBwAA4BHCHwAAgEcIfwAAAB4h/AEAAHiE5d1SYGJiwnR0dJi1tTX7dWBgIO+2W1tbZmpqyl6fn583V69ePbC99qFturq67NelpSVz8eJF09zcXJXXAgAAjhfhL+EU3gYHB01fX5+9resKgApvYa5fv24mJyft9eHhYdPS0mKWl5f3t5+bm7PbOOPj4/v7BgAAyUezb8KpihcMZ/39/fvhLpeqeq6yJ6rm6bljY2MHtltdXbWBMJvNmpGRkWN+BQAAoJqo/CXYzZs3Q5tjdX8+asbd2NjYf56qhLovSPcBAIB0IvwlmKvgBbW2ttpwF0ahbnNz81BQzG3WVdOvLC4ump6enoJ9CAEAQLIQ/iJQ8+f9x/er+j3PPHPGZDKZgtso5CnsBamiFxYKw6ysrNh9qF+fo75/boCHQl9nZ2fBPoRArXz66af286vPuz6v+pz++Mc/rvVhAUDsEf4iUPB7+d9erur3vPWPt0z9qfqC2+QGP3H/EUZx5coVs7CwcGD73CqgQp/6EObrRwhUkz6v+izevn3bhr0LFy7sf+4VBHW//njRoCWCIACEI/wlmEJbbhNvWDUwjP5znJ6ePlTR0wASjQJ29B9soT6EQDXcu3fPfmYV7G7cuFFwWwXA999/3zz33HPmjTfeqNoxAkBSEP4iNsGqElft71mMqnS5Tby6XWxqFgW8oaGh/eDn+v1pJLCmjtF1N+hD+2MACGod/PSHigJdFKoGvvPOOzYEvvfeewRAAMhB+ItAfe+KNcHWivrlBQdtaOLmYB8+9etTqHODNlwVT1VDN+2LtnGBT88Nhj1tPzs7W/XXBThNTU1lBTiFQIIfABxG+Es4BTM1h7kgp8pdsCl3ZmbGhjuFPz2ueQBzBcOittP+2tra7Hx/6l/FYA8AANKD8JcCwfBW6DFV+zRyuRBX/QPi7Nq1a3Yick1HJJ988oldrYZBHgBQHCt8AEic7u5u26/Peemll2yl+qOPPqrpcQFAElD5A5A4qlDnzoPZ29trPvjgg5L3NTExYfenrhP6WmxSc22n7hCaAF3UXzbq9EoAEAeEPwCJo36sWpNaI4EVvtSXVcFNzcC/+MUvIu9HfWQHBwf3B0zpeqFJzRX8tI2anN32uo81sAEkCc2+ABI7x+WTJ09s8Prmm2/sV4W5Umjao+DUSAqRhSY012Co4PcYHR09MC8mACQB4Q9Aoqm5V/3/vvzyy/11qaPQNEZhzbWFJjXX/oNhURVCmnwBJA3NvgAS5/Lly+bdd9+1lbqf//zn9j5NXF7KhORha2BrdZzcVXOCTc2iZl63Lrbkq/zt7e3Zi7O9vR352ADgOBH+ACRy4uc333zzUCAsJfyFLYWoKl5YKHShT/S4GxSi8Kl9hA0SUZ/Et99+O/LxAEC10OwLIBVee+01O+VLVGFrYCvYFWvGDQ4GUfhTyAuj/oAakOIu6+vrkY8NAI4TlT8AXg8aKVYNdFxVMbe66CqCuerq6uwFAOKGyh8AL2ngRm4Tr24HB3SEVfyCYe/bb78tqakZAOKA8AfAW+qrFxzdOz8/f2AqFw3sCI4g1nQywe31uJp3ASBJaPYF4K3Z2Vk7d5+qear6KfgF+/TNzMzYgOcGdGjda22vS1tbm53kudiKIAAQN4Q/AImgdXtfffXVSI9/+umn+1PAFKNAV8pjhbYHgCSg2RdAIrz//vsF58q7cePG/vVCq3QAgO8IfwASQX3r8k2r8t577x3oq5fNZqt4ZACQLIQ/AImgOfxWV1dDH/vTn/50YI6/TCZTxSMDgGQh/AFIVPVPVb4g3WbELQBER/gDkBiq7qnKV6jqBwAojPAHIFHUt++Pf/yjva6vwb5+AIDiCH8AEqW3t9dOxiwffvihvQ0AiI7wByBxhoeHzaVLl8y1a9dCH2e0LwDkR/gDkDiq9qmfX75Jn/v7+6t+TACQFIQ/AIn0zjvv5H3stddeq+qxAECSEP4AAAA8QvgDAADwCOEPAADAI8/U+gAAAMBhGrS+u3vUfWSLXg97bHeXJRLTjPAHAEDMKItdulRvbt2q/X/TJ0+erPUhoMJo9gUQe7nr+QJpp4pfHILfP/yDMefOERXSpvafLAAoYmZmxrzxxhu1PgygJv7rvx6bpqba/HddX29Mhhbg1CHOA4i95eXlSNU/KoRIo7Nna3ch+KUT4Q9A7M3OztpVPQqFuw8++MBMTk5W9bgAIIlo9gUQe27Fjo6ODvPWW2/tr+5x584dMz4+bm7cuGE2NzdNhjIFABRF5Q9AYjQ1NZnR0VFz6dIl09PTYzo7O838/LwNhAp/hZZ8AwB8j/AHIPZ+/etfm08//dR+bW1ttYFPc5G9//775ssvvzRvvvmmDYb6CgAojPAHIPbUl6+7u9uGPgU+VfmWlpbsYwqFAIDoCH8AYq+rq8uGPVX5rly5Yqt8ouurq6t2sAcAIBrCH4DYu3r1qnnppZfyDga5cOGCeffdd80vf/nLqh8bACQN4Q9A7KnCV4iCoS6aEgYAUBjhD0Aq9PX17U8JAwDIj/AHIDU0DQwAoDDCH4DUyNcvEADwV4Q/ALF27969stbs3d7eZq1fAAhB+AMQa5rWRX35rl27Zj766KNIz5menjbXr183b7zxRtFtJyYmzNzc3P7XqNbW1uz3AICkYW1fALGnqVw0ubNCnZZwa2lpscu7NTc328e3trbMt99+a5aXl+1jmhqm2Ahh0XaDg4N2sIjoutYP1ryCUZ6rbQEgaaj8AUgMBboPP/zQTE1N2UCoJd602oeqgwpwemxmZsa8+uqrkfan/bjgJ/39/XY1kWJu3rxJ8AOQWFT+ACS2KfgoFOBc5TD3/ihNvp2dnXZ1EQBIGip/ABJHff+Ouqavmopztba2mo2NjaLVwuHh4aL739vbs4NOghcAiAPCH4DEUb8/rfV7FAp5CntBqgSGhUJHj0Vt7h0bG7MVSnc5f/78kY4XACqF8AcgcTQw41e/+tWh+0uZ2iU3+LlwF9YU7Ny4ceNAH8FiE05rmhp3WV9fj3xsAHCc6PMHIHHU1+4nP/mJDWoXL16092nwx8LCQqTpXUTPzW3iDasGBvsCRg1+UldXZy8AEDeEPwCJoyA2MjJy6P5SmoIV5HKbeHW7UMALzgOoUcWi+QHVB7BQxRAA4oTwB+/oP2v129KITX0dGBjIu63CgDr4y/z8vJ3bLbi99qFtNC+cvip8qBJFEDhe4+Pjpre399D9pU6/op9lsKKnn7H27aysrNifsbbTNsFgqHkF9VhYCAWAOCP8wSulTuqrFRzcvG+q7mgCYU0k7LZXJSi4yoOCQylNgyjNBx98YL/me4/DAmEhs7Oz9ufnQrw+H8HPgqp7CoC5fyAEVwPRdQIggCTJZNVRpghNUaDRauq03NjYWJ0jA45BJpOxfcMcVfUU5sIm9lUgUBhQQHCVPE0CrOu6z/3Hr2Dgqn8oX6HzjO7r7u62PxPRHHv6uSXpfJS08+jOwx1zbuycvX73N3dNQ11DrQ/JKzs7xrS3f/852dp6bJqaqNWgcucZRvvCG+VM6qtm3OCgANdcHBR1OTCU76233rJBXAM9tKLHm2++GWn5NgDAYfwpAW+UOqmvQp2CRrERn675b3Fx0a43W6gPIcqjaq0Cn6Mm+NwQDgCIhvAXJ2qOfLRb66NInlP1as89lkl9g9T3S/sIDghQxc8N8FDoU3MklcDKC6vY6r0GAJSO8BcnCn6/a6/1USTPP31lzLNnj2VS3yA1M2oeueD2uVVAhT71HwzrQ4jyPffcc6H9N8MmeY46zx8A+IrwB2+UOqlvkEaETk9PH6ro5a7zqqpfoT6EKI9G3WqkdXCwjqZlyQ3y2o7wBwCFEf7i1nypKhZKf9+OaVJfF/CGhob2g5/r9+dGA+u6m1+ulLVfEZ1G9gb7/OULgBrhBgAojPAXJ2rGitB8ifKFTer7zjvvmAcPHtjbn332mbl79665dOmSvf3xxx/bJboaGhrMn//8ZzuM/osvvjCvvPKKaW9vN3/4wx/sV/f8zz//3Pz+97/fv42DTp06ZU6ePFny8zSPnn5OUUYFAwAKY54/eEdNuBosoCrdT3/6U7tGrPs10Ojehw8fmhdeeME8ffrUrK+vhzYf6/dBHj9+bL777jtz4sQJe/3s2bPm9OnTVX9NSaF+ei+++KI5d+77+eOinmc++eQT89JLLxXdf9TtaiFp51Hm+ast5vnDcZ5n+DTBO2607pMnT2w1r76+3jz//PM2mFy4cOHAtoworRwF7K+//tpWVhW4S6kARg10cQ1+ABAnhD9469GjRzaQKPidOXOm1ofjBb3Xd+7cse99Oc2/AICjY4UPeC9syhAcD95rAKg9wh8AAIBHCH8ACtLoaK1uAgBIB8IfkFJRl60r9pzgcnYAgORjwAfwAw3+uP/oSdW+35lTJyP3gZubmzuwOolWFdF9Y2NjdlLp2dlZ09/fb6ehUVjTBNSDg4N2u56eHnvbzZeXb3+q8Ok5WslkcXHRtLW12YmttZ1WzlhaWjqwmknuRNg6Ds2bqOMoNnE2AKB2CH/ADxT8fvbP/1617/f5v1wy9c8W/xVUk6vWCnarWShcXb582U5YrXCmwKaKne534U5BTAFMwU/bSXd3t91ey9mF7U/b66KwNzo6agOjtte+giuchFHgXF1dtc/XMmyaLxEAEE+EPyDmVHVTRc+tGazrqsK5ZeW0xFxvb69dAi2XtnUU4BT69Jx8+xP3eKGwl0vBTwHUVRh1Pfi9AQDxQfgDAs2wqsZV8/tFoSClCp4LZ7lNqm4t4eCydWEUxlT1K7Y/bZOPwl3Y2sVaNUUh1FUWAQDxxYAP4Afqf6dm2Gpdovb3U6hyTbSSW2HTZWFhwW6XO2AjeFv9AtVEXGh/YRTm1BSsbcIGhKj/oO53odD1JdT9AID4IfwBMadqmgLbxMSEDVSq8Clo6br68ami5y4Kd8FpWTRwQ9vpuXpMVb58+3NTuqhpOBgGta2Co7YNawp2g0L0fG2j/oKqBJbSbAwAqJ5M1q1on6IFyYEoHjx4YG7fvm3X8z19+rRJG1fli9PI20LvedrPM0l7fTsPd8y5sXP2+t3f3DUNdQ21PiSv7OwY097+/edka+uxaWqilxYqd56h8gcAAOARwh+QQq4JV8215Uz2DABIL+rIQAqpqVfTryDZ1Ctn99FuTb73zqMdU2vqlLRbm5dfc7u70QaEAeUg/AFATIPfK//6ivl4/WPjIwW/S5fqza1b/DcFVBrNvgAQQ6r4xSH4vfw3L5uzp85W/fuq4kfwM+bv/u6paWiINicoEBW/WQAQc5//r8/N2WfPHpobMni7lMdKUf9M/ZGeXwl37tw3DQ0HaxVhx1ToPaj0NtVSX3/C1PjtRwoR/gAg5hrPNHo91Upz8ymmOgEqiGZfAAAAjxD+ABx5WpnOzs4DK4sAAOKLOjqAouFO6/vmW67NLRmXVFrmTsvbaUk7fR0YGMi7reZMnJqaste1PrJWUCm0PQDEEeEPSCkFFa33e9TnjI+P20saKby5NY9F1xUA84VZrVmstY9leHjYtLS0mOXl5USHXwD+IfwBwYnFqjmh7ql6DSOMtOnc3JzZ2NjYv63gofvGxsZsWNFKHv39/Ta4KaipiqUgo+16enrsbRkZGcm7P1X49Jzp6WmzuLho2trabKjRdjMzM2Zpaclul4+20fdx+y20bVyoiufCnOg91O3gfY5emy4uIOui0Kifgd5/AEgKwh/gKPj9rr163++fvjLm2eLzp6kvncKImhldQLl8+bJtblQ4U2BTINH9LtwpECqYKPi5Zsnu7m67vZpww/an7XVReBsdHbVBR9trX0NDQ0WrW8HvpYqawqSrqMWRji+sMqr781EA1vvjnqf3RvcBQJIQ/oCYU0VNYcOFEl1X4FCwUvhQ0Ort7bXNj7mC4UYBTqFPz8m3P3GPl9qUGfxeCpqqhsU5/IWteaxgHKyIBul92dzcPHBfoYC7t7dnL8729vaRjxkAKoHwBwSbYVWNq+b3ixhSVFVzISM3bCiUSLFKm8KZwk2x/WmbfNygiDRQyMt9rXqPwkJhvoqs9pGvP6Sag99+++2KHCsAVBJTvQCO+t+pGbZal4j9/VTZc020onDi+vDpui4LCwt2u9zgErytSpyaiAvtL4yrhrn+bvkEH9P+9b3iLCzkljJI5sqVK/Z9z7e9ms7v3bu3f1lfXz/yMQNAJVD5A2JOza8KbG5KElHfOg3acKNP3QAEBS5VolyTrQZuiBsA4qp8YftT5dD1L9R+3WPaVvdpLj/XpzBfkHIDSdSPMM5NvqL3K7eJN6waGEbvjwbGFGoar6ursxcAiBvCH5AAYXPJ6b7g/WF9/vKFsLD9abvV1dVD9yvghI1+DSr2eBzp9YZVSouFVo0QDg6AifvAFgDIRbMvAG+5iqfjJm52VAlVNdMJDpJRNVWPs7IJgKSh8gekkGvCVT+/ixcvljzZsy/0/qgJ1/VnVPALNuVqpLXeR4VEN51OrrROgA0gvQh/QArla8KFKSm8BR9TgM5qInAASDiafQEAADxC+AMAAPAI4Q8AAMAjhD8AAACPEP4AAAA8QvgDAADwCOEPQKQ5AwEA6UD4A1Iqd+mycp/DJMYAkC5M8gz8QBP43n98v2rf78wzZ0wmk4m0rZYY29jY2L89PDxs7xsbGzMdHR12pQqtPqGJiBXWtGLF4OCg3a6np8felpGRkbz7U4VPz5menjaLi4umra3Nrnah7bTSxdLSkt2uUHVQx6Lv5b4PACB+CH/ADxT8Xv63l6v2/W794y1Tf6q+6HYKVZOTk3bdWVHIu3z5sl1yTOFMgc0tPeZCl0KYVvlQ8NN20t3dbbdvbW0N3Z+210Vhb3R01IY4ba99DQ0NHVj2LEjbaYm05eXl/e+j76nnAQDih/AHxJyqbqroqbomuq4qnIKaApbWo+3t7d0PX0HBNX0V4BT69Jx8+xP3eL6wl0tVRK0f7IQdBwAgPgh/QKAZVtW4an6/KFTVUwXPhTP31XEVNoW53MeCFOhU9Su2P22Tj6p8YRW9YMgEAMQbAz6AH6j/nZphq3WJ2t9PlT3XRCsKb64Pn67rsrCwYLfLHbARvK1+gWoiLrS/MAqDagrWNmEDQhQeXRVRtF2h/QEAaovwB8Scml8V2CYmJmwTq4KWqm+6rv51qrq5i8JdcFoWDdzQdnquHlNQy7c/N2hDTcPB8KZtFRy1bVhTsO5TH0H1+wvuDwAQT5mshjgWsb29bZqamsy9e/dMY2NjdY4MOGYPHjwwt2/fNhcuXDCnT582aeOqfIWaguP0nqf9PFPq69t5uGPOjZ2z1+/+5q5pqGswPtnZMaa9/fv3aWvrsWlqopcSUKnzDJU/AAAAjxD+gBRyTbhqri1nsmcAQHpRRwdSSE29q6urtT4MAEAMUfkDAADwCOEPAADAIzT7IlY0/YhbH1Zf3dJkYdSXbWpqyl7XvHUa2Rrc3s1Lp6lI9FWrWGglCiYkBgD4jPCH2FB4c3PRia4rAOZbZkzzymlOOhkeHjYtLS12aTG3veac0zbO+Ph4rKY9AQCgFmj2RWyoihcMZ/39/fvhLpdbRcKNZFU1T88dGxs7sJ0GPSgQajrLkZGRY34FAADEH+EPsZmaJKw5NrhsWC4142rZMcc1FwcVqhzi+Cq4WnkEABBPiW/2VUXn/qMntT4M5HHm1MlIa9iGzUXn1pQNo1C3ubl5KCjmNuuq6dctc9bT01OwDyHC6X3VzyJqiFbzem9v77EfFwDA0/Cn4Pezf/73Wh8G8vj8Xy6Z+meLf8wU8hQwglQJjDpBsSY01j4UPByFFTfAQ6Gvs7OzYCVQf0hk79831ZI5cyZSMC6X3rtSB7eEPUfvafB9BQAkW+LDH9IhN/iVGl6uXLliFhYWDmyfWwVU6FMfwnz9CBX8/qOres2Vf7uybDL19ZG2VQUzWAXVABfdpz6OCrRayUN9JPX6FdTU/K0BM9pOFU/XHO76PYbtTxU+PWd6etpWStva2ux7pu1mZmZsM7u2K7SiSFjTOwAgXp5JQ7OiqkuI788nCoWW3CbesGpgGI3oVWDJrehpAEkwrCiYFOpDGFcKVQqsms5GFPIuX75sq5l6zQpsCsq634U7vVaF32BTt/rhaXu9p2H70/a66H0fHR21IU7ba19DQ0N5K6baTj8DDaxxtxUWAQDxlPjwp2azKM2KiDeFjtwmXt0uNjWLAl4wmLh+fwogGnig6wovbn/uer5mWFXjqkXfLwoFKYVjF1x1XVU499r0OtXHzoWvoGAlVO+TQp+ek29/4h6P2sdPVcTgzylKYAcA1A6pCbGhClVw0IYqU8G+ZqqAKdS5SlYwvLhpX7SNC0V6bjDsaXs1jxb6QyJqM2w16XWpgufel9xA7F5j2ICXIL1PCmbF9lcovLnJtwEAyUX4Q2womKn50AU5VbSC1SdVwBTuFP5cM2euYFjUdtqf+q5pvj9VvZI47YveB70O16Sr166mWYUwXddF/R3VrKvqX7DaF6ymuvdX4S7f/sK4UdfBeRWD9D7rGJ1826VhVZlytgeAuCH8IVYKjSoNPqaAo9G5hbjqX9IpsCpcudAhChxuBROFWr0fuqj/n16zC7kauCFuAIir8oXtzw3a0P60X/eYttV9Gi0dNlG2a3rW/vR9XYVWxxf3YFTqqjKlbg8AcZTJFvsf1Bizvb1tmpqazL1790xjY2N1jgw4Zg8ePDC3b982Fy5cMKdPnzZpo2Di+j1WW76R2oXe81qcZ9TUHzwFqg+pqqf5RoSXuv1RXt/Owx1zbuycvX73N3dNQ12D8cnOjjHt7d+/T1tbj01TE7UKoFLnGX6bAFRcqfMLJmFVmXJWoTkKmzEfft8HdXcnY048Nl7Z3f3r9eIlCgClIPwBKeSacNXPz010jaOtKlPq9nt7e/YS/Iu85PDzux17/X/8znhN7wUfYaByCH9ACqmpV4NcULlVZUrdXhNwv/322xU8Yn8d50o4gI8IfwC8VOqqMqVur4myX3/99QOVv/Pnz0c+vuea6s1/b+yYp0+fmvpT9d4FIDX1uqbfH/0o2mTxAKIh/MF7EcY8IYXvdamrypS6fV1dnb2U68SJjPlRy1njM8YXAseD8AdvnTp1ylZTvv76a/P88897V1mpRfDTe633We990laVKXcVGgCIG8IfvHXy5Enz4osvmrt375o7d+7U+nC8oOCn91zvfRJXlSm2PQAkAfP8wXtPnjwxjx49qvVheEEVv3zBr1bnGU1orQms3drPwYmp9ZgCoEJelO0L4TwK4LhFPc8Q/gDEQtrPM2l/fQCSc545UdWjAgAAQE0R/gAAADxC+AMAAPBIpNG+rltgqcsTAUBU7vwSp7kAK4nzKIC4nEcjhb/vvvvOfi1ldnoAKIfON+qwnDacRwHE5TwaabSvlhf66quvTENDQ8kT4boljdbX1xM9wi0tryNNr4XXka7XolORTljt7e3mxIn09Ug5ynm0VtL02SwHr5/Xfz5hrz/qeTRS5U870MSsR6E3Lilvng+vI02vhdcRP+W+ljRW/Cp5Hq2VNH02y8Hr5/U3Juj1RzmPpu/PawAAAORF+AMAAPDIsYe/uro689vf/tZ+TbK0vI40vRZeR/yk6bWAnyevn9f/25S+/kgDPgAAAJAONPsCAAB4hPAHAADgkUhTvZRrYmLCdHR0mLW1Nft1YGDAJNHW1paZmpqy1+fn583Vq1cT+1oc/UwmJyfN+Pi4SfLx9/T02Nt9fX2mubnZJJF+T3Ts+pwl5ffEvf+dnZ1meHg4tb/7vkrjOc/Xc2W50nSO9eW8XJLsMRkeHs7Oz8/v3x4YGMguLy9nk0ivxdnc3FQfycS+Fqevr+/A60qS1dXVbFdX14HP1vj4eDaJco9bt+P+2dLv9ezsrP0Mhb3vafrd91Uaz3k+nivLlaZzrC/n5VIdW7Ov/mrUXwpOf3+//SsiiX/96KL0L/pLQK9rbGzMJNXNmzftXzJJdf36dVuJcEZHR0OrT0mwuLh44LY+W/r5xJmOUX8F56sCpOV331dpPOf5eq4sV5rOsb6cl0t1LOFPb1LYfwxJffOWlpbMxsbG/m3XnJVUOnY11yXV3NzcgXDR1dWV2OaIlZWV/eY191nT60mqtP3u+ypt5zxfz5XlStM5thwrKTsvVy38ub8Yg1pbWw+cTJJCJ73Nzc0Df/3pP7KLFy+aJNIHOsl/wemX0p2UdYLS6wn+kib1L+zBwUH72vS7EzzpJk2afvd9lbZznq/nynKl7RxbjuspOy9XbcCHTvQ64Qe5jpNJpw+CXl8SO/+6jqtJ5qoPei2uA66aFfV5S2KHXP3nsry8bE+u+g92YWHBJFmaf/d9leRzns/nynKl7RxbjuGUnZerVvnLPfm7D1IaysZXrlyxH4QkvpYbN26k5q+XYAleJ6ak9kfSyUXHv7q6aisr3d3diW4iTfPvvq+SfM4rV5rOlb6fY8sxlbLzctXCn04Suc08YRWBJJaCp6enE9n2rw9uGk5m7q/x3L/Kk9gfScesaTT017Rej66ruqLPWVKl9XffV0k+5/l+rixXms6x5VhL4Xm5auFPvzi5zTxJbzPXXwJDQ0P7J8Ek/hWg/huau0iXmZkZ24lV15PUJOfe/+CJ6Ntvv01kE42a09wcWs7IyIhJsjT+7vsqDec8n8+V5UrTObYcKyk8L4c5tqlelJqDJws3UWgSudehqoZ+IfThcJ1ik0L/+eoD7C66rV9mXU9ac46OOfjZ0s9CUxEkjX4G+r0I0udL/+EmWZp+932VhnNeudJ0rvT9HFuOvpSel3NlNNnfce1cZVINk0/yDNk69paWlkP3qwyc1L8G9Besm3dN/ykn8XW4EnxbW5s9ISd1VJ5OqqosBKeTiPtr0X8KOm71AdLvtU6KrokkTb/7vkrjOc/nc6Xv51hfzsuxCn8AAADwpNkXAAAA8UP4AwAA8AjhDwAAwCOEPwAAAI8Q/gAAADxC+AMAAPAI4Q8AAMAjhD8AAACPEP4AAAA8QvgDAADwCOEPAADAI8/U+gAAAMBfbW1tmRs3buxf7+rqMmtra2Z5edlMTk7W+vCQAoQ/AABiZGpqyoyMjNjrLS0tZnZ21nR0dBD8UDGZbDabrdzuAADAUaysrNhqn2QyGcN/06g0+vwBABAjLvjdvHnT9PX11fpwkEKEPwAAYmh+ft709/cfqAgClUD4AwAgJubm5szg4OD+dfX1c1XA5ubmGh8d0oIBHwAAxITCni4TExO28ucGeSj4uSAIHBUDPgAAADxCsy8AAIBHCH8AAAAeIfwBAAB4hPAHAADgEcIfAACARwh/AAAAHiH8AQAAeITwBwAA4BHCHwAAgEcIfwAAAB4h/AEAAHiE8AcAAGD88f8BEbuGvpRspbEAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, (ax0, ax1) = plt.subplots(ncols=2, figsize=(8, 3)) \n", "ax0 = plot_DS_structure_with_labels([a,b,c,d], masses=[0.25, 0.25, 0.25, 0.25], offset=1.5, ax=ax0)\n", "dss.plot(ax=ax1)\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.10" } }, "nbformat": 4, "nbformat_minor": 5 }