{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "\n", "Bad key \"text.kerning_factor\" on line 4 in\n", "/home/cmb-06/as/skchoudh/software_frozen/anaconda37/lib/python3.7/site-packages/matplotlib/mpl-data/stylelib/_classic_test_patch.mplstyle.\n", "You probably need to get an updated matplotlibrc file from\n", "http://github.com/matplotlib/matplotlib/blob/master/matplotlibrc.template\n", "or from the matplotlib source distribution\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "source": [ "%pylab inline\n", "import matplotlib\n", "import seaborn as sns\n", "PAPER_PRESET = {\"style\": \"ticks\", \"font\": \"Helvetica\", \"context\": \"paper\", \n", " \"rc\": {\"font.size\":20,\"axes.titlesize\":22,\n", " \"axes.labelsize\":22, 'axes.linewidth': 2,\n", " \"legend.fontsize\":20, \"xtick.labelsize\": 20,\n", " \"ytick.labelsize\":20, \"xtick.major.size\": 8.0,\n", " \"ytick.major.size\": 8.0, \"axes.edgecolor\": \"black\",\n", " \"xtick.major.pad\": 3.0, \"ytick.major.pad\": 3.0}}\n", "PAPER_FONTSIZE = 20\n", "sns.set(**PAPER_PRESET)\n", "fontsize = PAPER_FONTSIZE\n", "matplotlib.rcParams['pdf.fonttype'] = 42\n", "matplotlib.rcParams['ps.fonttype'] = 42" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "alpha = 0.05\n", "\n", "def probability_of_error(alpha, n):\n", " return 1-(1-alpha)**n\n", "\n", "number_of_tests = np.arange(1, 1000)\n", "p_errors = list(map(lambda n: probability_of_error(alpha, n), number_of_tests))\n", "n_errors = list(map(lambda n,error: n*error, number_of_tests, p_errors))\n" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[0.050000000000000044,\n", " 0.19500000000000006,\n", " 0.42787500000000034,\n", " 0.7419750000000005,\n", " 1.1310953125000012,\n", " 1.5894486562500012,\n", " 2.111638927343752,\n", " 2.692636549687502,\n", " 3.3277553124785184,\n", " 4.012630607616213,\n", " 4.743198984958943,\n", " 5.51567894804836,\n", " 6.326552917366437,\n", " 7.172550292382587,\n", " 8.050631547603704,\n", " 8.957973301571755,\n", " 9.891954301273989,\n", " 10.850142267752073,\n", " 11.830281551829163,\n", " 12.830281551829161,\n", " 13.84820584794959,\n", " 14.882262010578401,\n", " 15.930792042324457,\n", " 16.992263415869463,\n", " 18.065260671954157,\n", " 19.148477543890706,\n", " 20.240709576953716,\n", " 21.3408472128507,\n", " 22.4478693112156,\n", " 23.56083708171188,\n", " 24.6788884018805,\n", " 25.801232497327973,\n", " 26.927144962226,\n", " 28.0559630993909,\n", " 29.18708156043375,\n", " 30.319948267623836,\n", " 31.45406060019383,\n", " 32.58896182883776,\n", " 33.72423778311682,\n", " 34.85951373739588,\n", " 35.994451501789236,\n", " 37.12874670539976,\n", " 38.26212625989477,\n", " 39.39434599217678,\n", " 40.52518843558085,\n", " 41.65446076966407,\n", " 42.781992899250014,\n", " 43.9076356639532,\n", " 45.031259169937954,\n", " 46.15275123616434,\n", " 47.27201594784325,\n", " 48.38897231026385,\n", " 49.503552996572786,\n", " 50.61570318347517,\n", " 51.725379469195886,\n", " 52.832548868385835,\n", " 53.937187878983806,\n", " 55.039281616351005,\n", " 56.138823010284035,\n", " 57.235812060782884,\n", " 58.33025514870614,\n", " 59.42216439768511,\n", " 60.51155708389441,\n", " 61.59845509048858,\n", " 62.68288440371359,\n", " 63.764874647889876,\n", " 64.84445865663925,\n", " 65.92167207789396,\n", " 66.9965530103816,\n", " 68.06914166942575,\n", " 69.13948008003953,\n", " 70.20761179541836,\n", " 71.27358163905922,\n", " 72.33743546884743,\n", " 73.39921996155918,\n", " 74.45898241632764,\n", " 75.51677057571536,\n", " 76.57263246312348,\n", " 77.6266162353515,\n", " 78.67877004919892,\n", " 79.72914194107321,\n", " 80.77777971863706,\n", " 81.82473086359187,\n", " 82.87004244475459,\n", " 83.91376104064206,\n", " 84.9559326708289,\n", " 85.99660273539544,\n", " 87.03581596182828,\n", " 88.07361635877932,\n", " 89.11004717613069,\n", " 90.14515087084999,\n", " 91.17896907815702,\n", " 92.21154258755622,\n", " 93.24291132332012,\n", " 94.27311432903873,\n", " 95.30218975587718,\n", " 96.33017485420919,\n", " 97.3571059683183,\n", " 98.38301853388099,\n", " 99.4079470779666,\n", " 100.43192522130896,\n", " 101.45498568262215,\n", " 102.47716028475077,\n", " 103.49847996245995,\n", " 104.51897477168635,\n", " 105.53867390008395,\n", " 106.55760567871259,\n", " 107.57579759472814,\n", " 108.59327630494536,\n", " 109.61006765015404,\n", " 110.6261966700795,\n", " 111.641687618887,\n", " 112.65656398113858,\n", " 113.67084848811777,\n", " 114.6845631344462,\n", " 115.69772919492148,\n", " 116.71036724151313,\n", " 117.7224971604583,\n", " 118.73413816940518,\n", " 119.74530883455621,\n", " 120.75602708776864,\n", " 121.76631024357343,\n", " 122.77617501607833,\n", " 123.7856375357238,\n", " 124.79471336586452,\n", " 125.80341751915188,\n", " 126.81176447369582,\n", " 127.8197681889875,\n", " 128.82744212156575,\n", " 129.83479924041373,\n", " 130.841852042073,\n", " 131.8486125654653,\n", " 132.85509240641318,\n", " 133.86130273185262,\n", " 134.86725429373206,\n", " 135.87295744259396,\n", " 136.8784221408353,\n", " 137.88365797564606,\n", " 138.88867417162365,\n", " 139.89347960306435,\n", " 140.89808280593192,\n", " 141.90249198950508,\n", " 142.9067150477061,\n", " 143.91075957011324,\n", " 144.91463285266042,\n", " 145.91834190802757,\n", " 146.9218934757264,\n", " 147.92529403188524,\n", " 148.92854979873889,\n", " 149.93166675382744,\n", " 150.9346506389103,\n", " 151.93750696860033,\n", " 152.94024103872405,\n", " 153.94285793441392,\n", " 154.945362537938,\n", " 155.94775953627362,\n", " 156.95005342843083,\n", " 157.95224853253166,\n", " 158.9543489926513,\n", " 159.95635878542765,\n", " 160.95828172644477,\n", " 161.96012147639658,\n", " 162.96188154703708,\n", " 163.96356530692256,\n", " 164.9651759869519,\n", " 165.966716685711,\n", " 166.9681903746268,\n", " 167.96959990293675,\n", " 168.97094800247916,\n", " 169.97223729230996,\n", " 170.97347028315147,\n", " 171.9746493816781,\n", " 172.97577689464416,\n", " 173.97685503285942,\n", " 174.97788591501651,\n", " 175.9788715713758,\n", " 176.97981394731158,\n", " 177.98071490672535,\n", " 178.9815762353295,\n", " 179.98239964380636,\n", " 180.98318677084725,\n", " 181.98393918607454,\n", " 182.98465839285197,\n", " 183.9853458309865,\n", " 184.98600287932538,\n", " 185.98663085825297,\n", " 186.98723103208945,\n", " 187.98780461139665,\n", " 188.98835275519292,\n", " 189.9888765730811,\n", " 190.98937712729244,\n", " 191.98985543464994,\n", " 192.99031246845348,\n", " 193.99074916029002,\n", " 194.99116640177178,\n", " 195.99156504620467,\n", " 196.9919459101898,\n", " 197.99230977516092,\n", " 198.99265738885944,\n", " 199.99298946675023,\n", " 200.9933066933798,\n", " 201.9936097236795,\n", " 202.99389918421582,\n", " 203.99417567438928,\n", " 204.99443976758485,\n", " 205.9946920122749,\n", " 206.994932933077,\n", " 207.9951630317682,\n", " 208.99538278825756,\n", " 209.99559266151857,\n", " 210.99579309048286,\n", " 211.9959844948969,\n", " 212.99616727614335,\n", " 213.9963418180279,\n", " 214.99650848753365,\n", " 215.99666763554373,\n", " 216.996819597534,\n", " 217.99696469423634,\n", " 218.99710323227464,\n", " 219.99723550477353,\n", " 220.9973617919418,\n", " 221.99748236163134,\n", " 222.9975974698721,\n", " 223.99770736138467,\n", " 224.9978122700713,\n", " 225.99791241948583,\n", " 226.9980080232837,\n", " 227.9980992856531,\n", " 228.9981864017273,\n", " 229.99826955798,\n", " 230.99834893260308,\n", " 231.99842469586892,\n", " 232.99849701047665,\n", " 233.99856603188394,\n", " 234.99863190862433,\n", " 235.99869478261095,\n", " 236.99875478942738,\n", " 237.99881205860558,\n", " 238.99886671389245,\n", " 239.99891887350407,\n", " 240.99896865036982,\n", " 241.99901615236524,\n", " 242.9990614825352,\n", " 243.99910473930723,\n", " 244.99914601669573,\n", " 245.9991854044971,\n", " 246.99922298847662,\n", " 247.99925885054694,\n", " 248.999293068939,\n", " 249.99932571836553,\n", " 250.99935687017702,\n", " 251.99938659251148,\n", " 252.99941495043706,\n", " 253.99944200608877,\n", " 254.99946781879922,\n", " 255.9994924452234,\n", " 256.9995159394582,\n", " 257.9995383531564,\n", " 258.9995597356362,\n", " 259.9995801339851,\n", " 260.99959959316004,\n", " 261.9996181560825,\n", " 262.99963586372985,\n", " 263.9996527552222,\n", " 264.999668867906,\n", " 265.99968423743337,\n", " 266.99969889783824,\n", " 267.99971288160907,\n", " 268.9997262197582,\n", " 269.99973894188844,\n", " 270.9997510762562,\n", " 271.99976264983246,\n", " 272.99977368836045,\n", " 273.9997842164111,\n", " 274.9997942574358,\n", " 275.99980383381694,\n", " 276.9998129669164,\n", " 277.99982167712136,\n", " 278.99982998388856,\n", " 279.99983790578625,\n", " 280.9998454605344,\n", " 281.9998526650433,\n", " 282.99985953544996,\n", " 283.99986608715335,\n", " 284.9998723348478,\n", " 285.9998782925549,\n", " 286.99988397365416,\n", " 287.99988939091213,\n", " 288.9998945565102,\n", " 289.9998994820711,\n", " 290.9999041786847,\n", " 291.9999086569317,\n", " 292.99991292690737,\n", " 293.99991699824307,\n", " 294.99992088012795,\n", " 295.99992458132874,\n", " 296.99992811020917,\n", " 297.9999314747482,\n", " 298.9999346825578,\n", " 299.99993774089955,\n", " 300.99994065670074,\n", " 301.9999434365696,\n", " 302.9999460868098,\n", " 303.9999486134346,\n", " 304.9999510221798,\n", " 305.99995331851693,\n", " 306.99995550766494,\n", " 307.9999575946018,\n", " 308.99995958407584,\n", " 309.999961480616,\n", " 310.9999632885419,\n", " 311.99996501197376,\n", " 312.999966654841,\n", " 313.99996822089156,\n", " 314.9999697137,\n", " 315.9999711366754,\n", " 316.9999724930689,\n", " 317.99997378598147,\n", " 318.9999750183701,\n", " 319.99997619305486,\n", " 320.99997731272526,\n", " 321.99997837994596,\n", " 322.9999793971628,\n", " 323.9999803667081,\n", " 324.999981290806,\n", " 325.9999821715772,\n", " 326.9999830110444,\n", " 327.99998381113585,\n", " 328.9999845736906,\n", " 329.999985300462,\n", " 330.99998599312204,\n", " 331.99998665326495,\n", " 332.99998728241076,\n", " 333.9999878820088,\n", " 334.999988453441,\n", " 335.99998899802495,\n", " 336.99998951701696,\n", " 337.99999001161467,\n", " 338.9999904829601,\n", " 339.9999909321419,\n", " 340.9999913601981,\n", " 341.9999917681184,\n", " 342.9999921568462,\n", " 343.9999925272808,\n", " 344.99999288027993,\n", " 345.99999321666087,\n", " 346.99999353720307,\n", " 347.99999384264936,\n", " 348.99999413370807,\n", " 349.9999944110542,\n", " 350.9999946753315,\n", " 351.99999492715347,\n", " 352.9999951671049,\n", " 353.99999539574327,\n", " 354.9999956136,\n", " 355.99999582118176,\n", " 356.99999601897133,\n", " 357.999996207429,\n", " 358.9999963869935,\n", " 359.99999655808296,\n", " 360.9999967210959,\n", " 361.9999968764124,\n", " 362.99999702439453,\n", " 363.9999971653874,\n", " 364.99999729972,\n", " 365.9999974277059,\n", " 366.9999975496439,\n", " 367.9999976658188,\n", " 368.99999777650214,\n", " 369.9999978819526,\n", " 370.99999798241674,\n", " 371.9999980781296,\n", " 372.99999816931506,\n", " 373.99999825618676,\n", " 374.99999833894793,\n", " 375.99999841779254,\n", " 376.9999984929053,\n", " 377.9999985644623,\n", " 378.9999986326314,\n", " 379.9999986975723,\n", " 380.9999987594377,\n", " 381.99999881837255,\n", " 382.9999988745153,\n", " 383.99999892799786,\n", " 384.9999989789459,\n", " 385.99999902747913,\n", " 386.9999990737116,\n", " 387.9999991177522,\n", " 388.9999991597045,\n", " 389.99999919966706,\n", " 390.9999992377342,\n", " 391.99999927399546,\n", " 392.99999930853625,\n", " 393.99999934143796,\n", " 394.9999993727781,\n", " 395.9999994026307,\n", " 396.9999994310661,\n", " 397.99999945815136,\n", " 398.99999948395043,\n", " 399.99999950852424,\n", " 400.99999953193077,\n", " 401.9999995542253,\n", " 402.99999957546066,\n", " 403.9999995956868,\n", " 404.9999996149517,\n", " 405.99999963330094,\n", " 406.99999965077785,\n", " 407.9999996674238,\n", " 408.99999968327825,\n", " 409.9999996983787,\n", " 410.9999997127609,\n", " 411.9999997264589,\n", " 412.9999997395052,\n", " 413.9999997519307,\n", " 414.999999763765,\n", " 415.9999997750359,\n", " 416.9999997857704,\n", " 417.9999997959938,\n", " 418.9999998057305,\n", " 419.9999998150035,\n", " 420.99999982383486,\n", " 421.99999983224563,\n", " 422.99999984025567,\n", " 423.99999984788417,\n", " 424.9999998551491,\n", " 425.99999986206785,\n", " 426.99999986865686,\n", " 427.9999998749318,\n", " 428.99999988090764,\n", " 429.9999998865985,\n", " 430.9999998920181,\n", " 431.9999998971791,\n", " 432.9999999020941,\n", " 433.9999999067746,\n", " 434.9999999112318,\n", " 435.9999999154763,\n", " 436.99999991951836,\n", " 437.99999992336745,\n", " 438.9999999270329,\n", " 439.9999999305233,\n", " 440.99999993384716,\n", " 441.9999999370123,\n", " 442.99999994002627,\n", " 443.99999994289635,\n", " 444.9999999456294,\n", " 445.99999994823185,\n", " 446.99999995070993,\n", " 447.99999995306973,\n", " 448.9999999553167,\n", " 449.99999995745634,\n", " 450.9999999594937,\n", " 451.99999996143373,\n", " 452.99999996328097,\n", " 453.9999999650399,\n", " 454.99999996671477,\n", " 455.9999999683095,\n", " 456.999999969828,\n", " 457.9999999712739,\n", " 458.99999997265064,\n", " 459.99999997396145,\n", " 460.9999999752096,\n", " 461.99999997639804,\n", " 462.9999999775296,\n", " 463.99999997860704,\n", " 464.9999999796329,\n", " 465.99999998060963,\n", " 466.99999998153964,\n", " 467.9999999824251,\n", " 468.99999998326814,\n", " 469.9999999840709,\n", " 470.9999999848351,\n", " 471.9999999855628,\n", " 472.9999999862556,\n", " 473.9999999869152,\n", " 474.9999999875432,\n", " 475.9999999881411,\n", " 476.99999998871044,\n", " 477.99999998925244,\n", " 478.99999998976847,\n", " 479.99999999025977,\n", " 480.9999999907274,\n", " 481.99999999117273,\n", " 482.9999999915967,\n", " 483.9999999920003,\n", " 484.99999999238463,\n", " 485.9999999927505,\n", " 486.99999999309875,\n", " 487.9999999934304,\n", " 488.9999999937461,\n", " 489.9999999940466,\n", " 490.99999999433277,\n", " 491.99999999460516,\n", " 492.9999999948645,\n", " 493.99999999511135,\n", " 494.9999999953464,\n", " 495.99999999557014,\n", " 496.9999999957831,\n", " 497.99999999598595,\n", " 498.999999996179,\n", " 499.99999999636276,\n", " 500.99999999653767,\n", " 501.9999999967043,\n", " 502.9999999968628,\n", " 503.99999999701373,\n", " 504.99999999715743,\n", " 505.9999999972942,\n", " 506.9999999974244,\n", " 507.9999999975484,\n", " 508.99999999766635,\n", " 509.9999999977787,\n", " 510.99999999788565,\n", " 511.9999999979874,\n", " 512.9999999980843,\n", " 513.9999999981766,\n", " 514.9999999982643,\n", " 515.9999999983479,\n", " 516.9999999984275,\n", " 517.9999999985032,\n", " 518.9999999985753,\n", " 519.999999998644,\n", " 520.9999999987092,\n", " 521.9999999987714,\n", " 522.9999999988306,\n", " 523.999999998887,\n", " 524.9999999989407,\n", " 525.9999999989916,\n", " 526.9999999990403,\n", " 527.9999999990865,\n", " 528.9999999991305,\n", " 529.9999999991725,\n", " 530.9999999992124,\n", " 531.9999999992503,\n", " 532.9999999992865,\n", " 533.9999999993208,\n", " 534.9999999993537,\n", " 535.9999999993848,\n", " 536.9999999994145,\n", " 537.9999999994427,\n", " 538.9999999994695,\n", " 539.9999999994952,\n", " 540.9999999995196,\n", " 541.9999999995428,\n", " 542.9999999995647,\n", " 543.9999999995857,\n", " 544.9999999996057,\n", " 545.9999999996247,\n", " 546.9999999996429,\n", " 547.9999999996602,\n", " 548.9999999996766,\n", " 549.9999999996921,\n", " 550.999999999707,\n", " 551.9999999997211,\n", " 552.9999999997345,\n", " 553.9999999997475,\n", " 554.9999999997597,\n", " 555.9999999997713,\n", " 556.9999999997823,\n", " 557.9999999997927,\n", " 558.9999999998029,\n", " 559.9999999998123,\n", " 560.9999999998214,\n", " 561.99999999983,\n", " 562.9999999998382,\n", " 563.9999999998461,\n", " 564.9999999998535,\n", " 565.9999999998605,\n", " 566.9999999998673,\n", " 567.9999999998737,\n", " 568.9999999998798,\n", " 569.9999999998856,\n", " 570.9999999998912,\n", " 571.9999999998964,\n", " 572.9999999999014,\n", " 573.9999999999062,\n", " 574.9999999999106,\n", " 575.999999999915,\n", " 576.9999999999192,\n", " 577.999999999923,\n", " 578.9999999999268,\n", " 579.9999999999303,\n", " 580.9999999999337,\n", " 581.9999999999369,\n", " 582.99999999994,\n", " 583.9999999999429,\n", " 584.9999999999457,\n", " 585.9999999999483,\n", " 586.9999999999508,\n", " 587.9999999999532,\n", " 588.9999999999554,\n", " 589.9999999999576,\n", " 590.9999999999595,\n", " 591.9999999999616,\n", " 592.9999999999634,\n", " 593.9999999999652,\n", " 594.9999999999669,\n", " 595.9999999999685,\n", " 596.99999999997,\n", " 597.9999999999715,\n", " 598.9999999999728,\n", " 599.9999999999742,\n", " 600.9999999999753,\n", " 601.9999999999766,\n", " 602.9999999999777,\n", " 603.9999999999789,\n", " 604.9999999999799,\n", " 605.9999999999808,\n", " 606.9999999999817,\n", " 607.9999999999826,\n", " 608.9999999999835,\n", " 609.9999999999843,\n", " 610.999999999985,\n", " 611.9999999999857,\n", " 612.9999999999865,\n", " 613.9999999999872,\n", " 614.9999999999877,\n", " 615.9999999999883,\n", " 616.9999999999889,\n", " 617.9999999999894,\n", " 618.99999999999,\n", " 619.9999999999905,\n", " 620.9999999999909,\n", " 621.9999999999912,\n", " 622.9999999999918,\n", " 623.9999999999922,\n", " 624.9999999999925,\n", " 625.999999999993,\n", " 626.9999999999933,\n", " 627.9999999999936,\n", " 628.9999999999939,\n", " 629.9999999999942,\n", " 630.9999999999944,\n", " 631.9999999999948,\n", " 632.999999999995,\n", " 633.9999999999952,\n", " 634.9999999999955,\n", " 635.9999999999957,\n", " 636.9999999999959,\n", " 637.9999999999961,\n", " 638.9999999999964,\n", " 639.9999999999965,\n", " 640.9999999999967,\n", " 641.9999999999968,\n", " 642.9999999999969,\n", " 643.999999999997,\n", " 644.9999999999972,\n", " 645.9999999999974,\n", " 646.9999999999975,\n", " 647.9999999999976,\n", " 648.9999999999977,\n", " 649.9999999999978,\n", " 650.999999999998,\n", " 651.9999999999981,\n", " 652.9999999999981,\n", " 653.9999999999983,\n", " 654.9999999999983,\n", " 655.9999999999984,\n", " 656.9999999999985,\n", " 657.9999999999985,\n", " 658.9999999999986,\n", " 659.9999999999986,\n", " 660.9999999999987,\n", " 661.9999999999989,\n", " 662.9999999999989,\n", " 663.9999999999989,\n", " 664.999999999999,\n", " 665.9999999999991,\n", " 666.9999999999991,\n", " 667.9999999999991,\n", " 668.9999999999992,\n", " 669.9999999999992,\n", " 670.9999999999992,\n", " 671.9999999999992,\n", " 672.9999999999993,\n", " 673.9999999999993,\n", " 674.9999999999994,\n", " 675.9999999999994,\n", " 676.9999999999994,\n", " 677.9999999999994,\n", " 678.9999999999994,\n", " 679.9999999999995,\n", " 680.9999999999995,\n", " 681.9999999999995,\n", " 682.9999999999997,\n", " 683.9999999999997,\n", " 684.9999999999997,\n", " 685.9999999999997,\n", " 686.9999999999997,\n", " 687.9999999999997,\n", " 688.9999999999997,\n", " 689.9999999999997,\n", " 690.9999999999997,\n", " 691.9999999999998,\n", " 692.9999999999998,\n", " 693.9999999999998,\n", " 694.9999999999998,\n", " 695.9999999999998,\n", " 696.9999999999998,\n", " 697.9999999999998,\n", " 698.9999999999999,\n", " 699.9999999999999,\n", " 700.9999999999999,\n", " 701.9999999999999,\n", " 702.9999999999999,\n", " 703.9999999999999,\n", " 704.9999999999999,\n", " 705.9999999999999,\n", " 706.9999999999999,\n", " 707.9999999999999,\n", " 708.9999999999999,\n", " 709.9999999999999,\n", " 710.9999999999999,\n", " 711.9999999999999,\n", " 712.9999999999999,\n", " 713.9999999999999,\n", " 714.9999999999999,\n", " 715.9999999999999,\n", " 716.9999999999999,\n", " 717.9999999999999,\n", " 718.9999999999999,\n", " 719.9999999999999,\n", " 720.9999999999999,\n", " 721.9999999999999,\n", " 722.9999999999999,\n", " 723.9999999999999,\n", " 724.9999999999999,\n", " 725.9999999999999,\n", " 726.9999999999999,\n", " 727.9999999999999,\n", " 728.9999999999999,\n", " 730.0,\n", " 731.0,\n", " 732.0,\n", " 733.0,\n", " 734.0,\n", " 735.0,\n", " 736.0,\n", " 737.0,\n", " 738.0,\n", " 739.0,\n", " 740.0,\n", " 741.0,\n", " 742.0,\n", " 743.0,\n", " 744.0,\n", " 745.0,\n", " 746.0,\n", " 747.0,\n", " 748.0,\n", " 749.0,\n", " 750.0,\n", " 751.0,\n", " 752.0,\n", " 753.0,\n", " 754.0,\n", " 755.0,\n", " 756.0,\n", " 757.0,\n", " 758.0,\n", " 759.0,\n", " 760.0,\n", " 761.0,\n", " 762.0,\n", " 763.0,\n", " 764.0,\n", " 765.0,\n", " 766.0,\n", " 767.0,\n", " 768.0,\n", " 769.0,\n", " 770.0,\n", " 771.0,\n", " 772.0,\n", " 773.0,\n", " 774.0,\n", " 775.0,\n", " 776.0,\n", " 777.0,\n", " 778.0,\n", " 779.0,\n", " 780.0,\n", " 781.0,\n", " 782.0,\n", " 783.0,\n", " 784.0,\n", " 785.0,\n", " 786.0,\n", " 787.0,\n", " 788.0,\n", " 789.0,\n", " 790.0,\n", " 791.0,\n", " 792.0,\n", " 793.0,\n", " 794.0,\n", " 795.0,\n", " 796.0,\n", " 797.0,\n", " 798.0,\n", " 799.0,\n", " 800.0,\n", " 801.0,\n", " 802.0,\n", " 803.0,\n", " 804.0,\n", " 805.0,\n", " 806.0,\n", " 807.0,\n", " 808.0,\n", " 809.0,\n", " 810.0,\n", " 811.0,\n", " 812.0,\n", " 813.0,\n", " 814.0,\n", " 815.0,\n", " 816.0,\n", " 817.0,\n", " 818.0,\n", " 819.0,\n", " 820.0,\n", " 821.0,\n", " 822.0,\n", " 823.0,\n", " 824.0,\n", " 825.0,\n", " 826.0,\n", " 827.0,\n", " 828.0,\n", " 829.0,\n", " 830.0,\n", " 831.0,\n", " 832.0,\n", " 833.0,\n", " 834.0,\n", " 835.0,\n", " 836.0,\n", " 837.0,\n", " 838.0,\n", " 839.0,\n", " 840.0,\n", " 841.0,\n", " 842.0,\n", " 843.0,\n", " 844.0,\n", " 845.0,\n", " 846.0,\n", " 847.0,\n", " 848.0,\n", " 849.0,\n", " 850.0,\n", " 851.0,\n", " 852.0,\n", " 853.0,\n", " 854.0,\n", " 855.0,\n", " 856.0,\n", " 857.0,\n", " 858.0,\n", " 859.0,\n", " 860.0,\n", " 861.0,\n", " 862.0,\n", " 863.0,\n", " 864.0,\n", " 865.0,\n", " 866.0,\n", " 867.0,\n", " 868.0,\n", " 869.0,\n", " 870.0,\n", " 871.0,\n", " 872.0,\n", " 873.0,\n", " 874.0,\n", " 875.0,\n", " 876.0,\n", " 877.0,\n", " 878.0,\n", " 879.0,\n", " 880.0,\n", " 881.0,\n", " 882.0,\n", " 883.0,\n", " 884.0,\n", " 885.0,\n", " 886.0,\n", " 887.0,\n", " 888.0,\n", " 889.0,\n", " 890.0,\n", " 891.0,\n", " 892.0,\n", " 893.0,\n", " 894.0,\n", " 895.0,\n", " 896.0,\n", " 897.0,\n", " 898.0,\n", " 899.0,\n", " 900.0,\n", " 901.0,\n", " 902.0,\n", " 903.0,\n", " 904.0,\n", " 905.0,\n", " 906.0,\n", " 907.0,\n", " 908.0,\n", " 909.0,\n", " 910.0,\n", " 911.0,\n", " 912.0,\n", " 913.0,\n", " 914.0,\n", " 915.0,\n", " 916.0,\n", " 917.0,\n", " 918.0,\n", " 919.0,\n", " 920.0,\n", " 921.0,\n", " 922.0,\n", " 923.0,\n", " 924.0,\n", " 925.0,\n", " 926.0,\n", " 927.0,\n", " 928.0,\n", " 929.0,\n", " 930.0,\n", " 931.0,\n", " 932.0,\n", " 933.0,\n", " 934.0,\n", " 935.0,\n", " 936.0,\n", " 937.0,\n", " 938.0,\n", " 939.0,\n", " 940.0,\n", " 941.0,\n", " 942.0,\n", " 943.0,\n", " 944.0,\n", " 945.0,\n", " 946.0,\n", " 947.0,\n", " 948.0,\n", " 949.0,\n", " 950.0,\n", " 951.0,\n", " 952.0,\n", " 953.0,\n", " 954.0,\n", " 955.0,\n", " 956.0,\n", " 957.0,\n", " 958.0,\n", " 959.0,\n", " 960.0,\n", " 961.0,\n", " 962.0,\n", " 963.0,\n", " 964.0,\n", " 965.0,\n", " 966.0,\n", " 967.0,\n", " 968.0,\n", " 969.0,\n", " 970.0,\n", " 971.0,\n", " 972.0,\n", " 973.0,\n", " 974.0,\n", " 975.0,\n", " 976.0,\n", " 977.0,\n", " 978.0,\n", " 979.0,\n", " 980.0,\n", " 981.0,\n", " 982.0,\n", " 983.0,\n", " 984.0,\n", " 985.0,\n", " 986.0,\n", " 987.0,\n", " 988.0,\n", " 989.0,\n", " 990.0,\n", " 991.0,\n", " 992.0,\n", " 993.0,\n", " 994.0,\n", " 995.0,\n", " 996.0,\n", " 997.0,\n", " 998.0,\n", " 999.0]" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "n_errors" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAI4CAYAAACWfsh7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzde5hdZ1nw/++eUzJJM5NkkjSnhqShuQnFIlAolAItCFbkVVDAVwQBFUXAAoLoT0RavFB5VQ4tIIggIPxQQcADCCiUFptCsehbCunT0hzaZtI0SSczOUwyO3vm/WPt3Uymc9hrZu/ZM3u+n+uaa81e69nP3Czai5vndBdGRkaQJElqFi2NDkCSJKmWTG4kSVJTMbmRJElNxeRGkiQ1lbZGBzAfFAoFV11LktRgIyMjhWraOXIjSZKaiiM3OdR62/zOnTvZvn17TfvUGb7f+vL91pfvt758v/VV6/dbKFQ1YPMQR24kSVJTMbmRJElNxeRGkiQ1FZMbSZLUVExuJElSUzG5kSRJTcXkRpIkNRWTG0mS1FRMbiRJUlMxuZEkSU3F5EaSJDUVkxtJktRUTG4kSVJTMbmRJElNxeRGkiQ1FZMbSZLUVExuJElSUzG5kSRJTaWt0QEARMQLgWcAPw48FlgGfDql9NJp9LUReAdwJdAD7Ae+CFyTUuqrWdCSJGlOmisjN38AvI4sudk33U4iYitwK/BK4BbgPcAu4PXAzRHRM/NQJUnSXDYnRm6ANwL3AT8iG8G5fpr9fBBYA1yVUrqucjMi3l3+G+8EXj2zUCVJ0lw2J0ZuUkrXp5TuSimNTLePiDgfeA6wB/jAmMdvB44DL4uIpdMOVPNKV1dXo0NoaitWrGh0CE2tu7u70SE0Nd9vfbW1NXbsZK6M3NTCM8vXr6WUhkc/SCkdjYibyJKfJwNfn+3gFrLjg0XuPdDPsROnGR4ZBgp0tLdy+vQwpZFhFrW3PfR7YcyzsZ8na/vwZ+3s/+H9NeinVvE0Xz/3HN4/7X4aHfvc7wceOPbAHIqn2fpp4d6+3nka+9ztp+/oKR548ASb153D+sEiSzvbG/K/O82U3ET5eucEz+8iS262YXJTE8cHi+zuPcLJkyWGGRn3X5YCcOzkEENDJY4NlmAE2loLLOpo4fhgibbWAh0dRU4MloCzn439PFlb+5lf/czn2Ge3n+Ici2fh9jOfY5+tfgaHTvPpr+ykva2FNSs6WdTRwZMuXEsjNFNyUxlj7J/geeX+8vEeRsTVZNNXE9q5c+e0ApvIyZMna95nPXR0dNDZ2Qm003cSBk+WaG2FEydLDBwb4vTwSPYPefvD/2UZocCJwSJDpWGKxezZ0sXtnDh2mmKxdNbvY5/laWs/86uf+Ry7/SzMfuZz7LPVT++h4wyPQGl4hJNDJe6+9zBrlp5kcHCQ2dZMyc1UCuXruOt6UkpXA1eP+8VCYQRg+/btNQ1o586dNe+z1gaOnWJXbz8DR0cYPHWKU+URmK6lHfQdPcnJU6eB8j/kpx7+D/3hIyc4b00XSzrbzsr2l3a2PfT/BCZ6lqet/cyvfuZz7PazMPuZz7HPVj8b17TQUoDWlgKLO1rZel4Pmzc7cjNTlZGZiVaJdY1ppwkMFUvc98AAx0+V6B8Y5MixIksWt3H8xJkRmKMnTrGqu5PWQrYmfaJ/WUqlEW7bdYjLH7+B7iWLGCZbc7OovZWupcMMjwzT0d5Gd/n3sc/ytB37rKNG/dQqHvtp/N+0H/vxn9n69vPaFz6WAw+eYMu6c7jw/MadvtJMyU0qX7dN8PyC8nWiNTkiW0dzx57D9A2cZOmSdg70naRYLNF78OhZIzCnS8MMniqyomsxHW2ttLRAR3vbOP+yLKL7nHbuPXCMx25dxXlrl7F40ewsMNu3bx8bNmyYlb+1EN1///2sXduY/1e2EPT29rJ+/ZpGh9G0sve7vtFhNK277rqrYYuJobmSm8rZOM+JiJbRO6YiYhnwVGAQ+HYjgpvrhooleg8NcLDvFIf7Bzly9BT3PnCUTWu7ODFYoqWl5WEjMKeKw5w8dZr21hbWrVrGuSuX0tHe2uj/KA8ZGBgwuamjvr4+k5s66u/v939868j3W1+nT59u6N+fd8lNRLQDW4FiSunuyv2U0t0R8TWyHVGvBa4b9bVrgKXAh1NKx2cz3vng+GCRH+46xKniMAf7jtPW2srSxR20tLSw5/6jXLhlJQWg+5wO7jtwjIu2rua8tefM2giMJEl5zInkJiKeDzy//LHyfwWfEhEfL/9+KKX05vLvG4CdwF5g85iuXgPsAK6NiGeV210CXEE2HfXWesQ/nw0VS9yx+zD7Dx/nwOHjrO1ZRnt7gVNDpzlvzTkMHD/F/YdPcNH5PVx0wWo6TWgkSXPcnEhuyGpKvXzMvfPLP5AlMm9mCuXRm4s5UzjzuWSFM68lK5z5YM0ibgJDxRK79h1h74GjjIyMsHxZJz/ce5inPGYtrS0FHhw4xQXnLee8teeY1EiS5o05kdxMtg17nLZ7OLOte7zn95IVztQkhool0t4H2bv/KJ2L2jjUf4KuJR08evNK9h08zkUXrOL89d0mNZKkeWdO1JbS7KokNnff18/xk0MMnT7N6uVLGB6GjrYWHh+ruWDjChMbSdK8NCdGbjS77rn/KHfv6+dU8TRDp0sMnS6xdDGcs6SNjauXsWX98jm160mSpDwcuVlghooldvX2cWooS2wAOtraaGkpsH7VOTzyvBUmNpKkec3kZgEZKpbYdV8/J08NM1TKEpv2tjY62lt45MZutm1aaWIjSZr3TG4WkHsPHOX2XYdY1NECI9De2kZbC2xeZ2IjSWoeJjcLxFCxxJ7efoqlEj+4+xDLlizidGmY89Z2sW2TU1GSpOZhcrNA3HdggOMnT8MItLW1sfu+fs5duYRzVyxuaP0PSZJqzd1SC8Rd9x6hWF5AfO7KJbQUoKd7MetXd03xTUmS5hdHbhaAoWKJ4RG49Y4DHBss0nfsJBvXLuMxW1c5HSVJajomNwvAfQcGaG0tsGHNMg4fOcmKcxbT3tpiYiNJakpOSzW5oWKJH+zq4+RQkUVtLZSWZOtrVixb1ODIJEmqD0dumtx9BwYolkr8950PMAK0txRYt2op553rWhtJUnMyuWlilVGbUmmYjWuW0Xf0FGtXncN5a5Y6JSVJalomN01s9KhNe1sLyzrb6elezIY1jtpIkpqXyU2TctRGkrRQmdw0qX0PHGOEYX6w6zCL29tYtsRRG0nSwmBy06T2HzpGR3srG9cso//4KTasXuaojSRpQTC5aVItrQVuu+sgbW0tLO5oZZWjNpKkBcLkpkmdPFXi3J6lHD0xxOqVSygNjzhqI0laEExumtBQscTQ6WFuu+sgD/af5Ps/OsRwo4OSJGmWmNw0od37+hkeHmb75h4Wd7TymK2r2Lh6aaPDkiRpVlh+oQn9aF8/P7qnj9UrOxlhMT3diz2RWJK0YDhy02SGiiXa21q4894+/u+d2bRUR3ur620kSQuGIzdNZu/+AdrbClxy4VoOHhnkEeu6WNezpNFhSZI0a0xumsz9h4+z7+AxlnZ20NpS4NwVS5ySkiQtKCY3TWSoWALg5tv2c05nOyu7FvOozT1OSUmSFhSTmyZy74FjHO4f5EkXruVQeUpKkqSFxgXFTeTA4WMUSyPcdtdBjp0ocs/+oyxb2tHosCRJmlUmN02ktbWV/Q8c46ILVnPOknaeeOG5bHb0RpK0wJjcNJFCYYT1557DyaHTbN+ykiWL2lxvI0lacExumsjISIEf3H2YI0eH+O4PDlCy5oIkaQFyQXGTGCqWGDxVJB6xkiNHT/LoLStpbSk0OixJkmadIzdNYu/+AY4NFkl7H6T/+BC9h465mFiStCA5ctMkDh4ZBOAJ28/l0JFBtmzodjGxJGlBMrlpEiuWLeK/dh6gNDzC8nMWsXbFEhcTS5IWJKelmkR7Wwvnb+hmxbJFrF7RSVub620kSQuTyU0TGCqWSPcc4a57jzAyAsdODHHgwcFGhyVJUkOY3DSBvfsHGBkZ4eCDJ7jznj4OHhlk1fLORoclSVJDuOamCVQWE1/62PUuJpYkLXgmN03AxcSSJJ3htFQTcDGxJElnmNzMcy4mliTpbCY385yLiSVJOptrbuY5FxNLknQ2k5t5zsXEkiSdreppqYj4XkR8rp7BKD8XE0uSdLY8a24eBRTrFYim54G+QQZPnWbFssUMnjrtYmJJ0oKXZ1pqL3BOvQJRfkPFEq0tLeza10/nojZWLFvkYmJJ0oKXZ+TmH4GnR8SqegWjfHb3DnDwyAkuPL+HrqUdxCNWuphYkrTg5Ulu3gncAXw1Ii6pUzzK4XD/IAPHhx6alhoZGXExsSRpwcszLfUloAQ8EdgREQfIpqrGW+QxklJ6Vg3i0ySynVLZmpvORW1s27S80SFJktRweZKby0f9XgDWln/GMzLdgFS9yk6pQ+WD+zraPZNRkqQ8yc0VdYtC0zLeTqmtG1c0OixJkhqq6uQmpXRDPQNRfpUD/CrTUo85v6fRIUmS1HCeUDyPOS0lSdLDTSu5iYgNwNOBDeVb+4AbU0r7ahWYpua0lCRJD5cruYmI5cAHgBfz8G3kwxHx98DrUkpHahSfJuABfpIkjS9PbalO4BvA/ybbLfUd4P8HPgN8u3zvF4Gvl9uqjjzAT5Kk8eUZuXkD8OPADuBVKaWdox9GxHbgw8BTgauAd9UqSD1c5QC/9rYWD/CTJGmUPCtQXwz0AT89NrEBKN/7GeAI2eiO6qinOxscK54e5vjJIiu6Fjc4IkmS5oY8yc0FwPUppf6JGpTX2lxfbqs62rjmHLZuXM6i9ha2blzOhtXWNJUkCfIlNyNk62o0B9z3wDHuuX+AttZW7rl/gH0HjzU6JEmS5oQ8yc3dwDMiYtlEDSKii6xMw49mGJemcLh/8KEpqeLpYQ4dGa/ElyRJC0+e5OazwErgnyPiYdNOEfFI4AvACuAfahOeJlJZc1PhNnBJkjJ5dku9B/gF4BnADyPi28Busumq84EnA63A94H31jhOjZKdcQMb1pzDsRNFzl/f5TZwSZLKqh65SSmdICue+bny954KvBR4Wfn3lvKzZ5bbqk529w5wyw8P8MCDJzh9epj2tha3gUuSVJbrhOKU0mHgxRGxCXgaWfmFAnAf8K2U0j21D1FjHe7P1tcUTw9TPD3MA32DPPI8yy5IkgQ5kpuIuAgYTindXk5iPl2/sDSZnu5O9uwfeOiz620kSTojz4Li/wHeX69AVD3PuJEkaWJ5kpsjZNNPajDPuJEkaWJ5R2621isQVc8zbiRJmlie5OZa4JKIuLJewag6nnEjSdLE8uyW+h7Zmpt/ioiPkR3YtxcYd9jAnVP14Rk3kiRNLk9ys7t8LQC/Xv6ZyEjOvlWl3b0D3HrHgexsm7ZWz7iRJGmMPAnIvWRJixrIM24kSZpc1clNSmlzHeNQlTzjRpKkyeU5xK8LGEkpHa1jPJrCxjXncOTYKe47MMDGc7s840aSpDHynnPzH/UKRNXxjBtJkiaXJ7k5CtxVr0BUHc+4kSRpcnmSm53AxnoFoup4xo0kSZPLs1vqI8CHI+IJKaVb6xFMRGwE3gFcCfQA+4EvAteklPpy9PPTwOuBR4/q51bg3Smlm2sd92zasr6LQgEOHRlk1fJOz7iRJGmMqkduUkofBT4I/HtE/G5EbIuIRbUKJCK2kiUgrwRuAd4D7CJLUm6OiJ4q+3kX8K/A44GvAO8jO4DwZ4GbIuKltYp5tg0VS+zuHaBv4CTnrlzK5nVdnnEjSdIYeXZLlUZ9/OPyDxExXvORlFLeQ/w+CKwBrkopXTfq774beCPwTuDVU8S4FngzcAC4KKX0wKhnVwDfIBsZ+lTO2OaEygF+AHfv66c0fC7bNnnGjSRJo+VZc1PI8ZOnXyLifOA5wB7gA2Mevx04DrwsIpZO0dUjyn/7O6MTG4CU0vVki6JX54ltLqkc4FfhYmJJkh4uzyF+uRKWnJ5Zvn4tpTQ85u8ejYibyJKfJwNfn6Sfu4Ah4EkRsSqldKjyICKeDiwjW8MzL3mAnyRJU5sr9Z8qc1t3TvD8LrLkZhuTJDcppQcj4neBdwM/jIgvAoeBrcDPAP8O/Eatgp5tHuAnSdLU5kpy012+9k/wvHJ/+VQdpZTeGxF7gI8Brxr16EfAx8dOV1VExNVkU2AT2rlz51R/PpeTJ09W3WdLSwt9g+38d7qfRR1tPHDwEBRPsLyzyPDw8NQdLEB53q/y8/3Wl++3vny/9dXo95s7uYmIR5KNfjyFbP3KP6WU3lJ+9mTgIuAfUkpHahhnoXydsnBnRLyFbLHztcD7gfuBRwF/Anw6In68Eu9oKaWrgavH/eOFwgjA9u3bpxH6xHbu3Jmrzx239dK9fCUAi5dAqWUx2XIljSfv+1U+vt/68v3Wl++3vhr9fnMlNxHxq2QLfjvKt0aAVaOarAb+EigCf5Oj68rITPcEz7vGtJsovsuBdwFfSCn99qhH34uIF5BNe70pIj6UUtqVI745wTU3kiRNrepFwhHxVODDwEngd4BLODOiUvEVYIBsfUseqXzdNsHzC8rXidbkVDyvfL3+YX8gpRNk5+e0AI/LGd+csGV9FxdvP5fN685cJUnS2fLsgHoL2UjNT6WU/iKl9N2xDVJKRbJEJe9YVCUZeU5EnBVTRCwDngoMAt+eop/KoYITbfeu3B/KGd+cUQAWtbc+LKuUJEmZPMnNU4BbqihfcC+wLk8QKaW7ga8Bm4HXjnl8DbAU+GRK6ThARLRHxKPKpxqP9q3y9dcjYsPoBxHxU2RJ0klgR5745ordvQN8d+cB0j19fHfngbOmqCRJUibPmptu4L4q2nXk7LfiNWRJx7UR8SyyQp2XAFeQTUe9dVTbDeXne8kSoorPAf8B/ASwMyK+QLageDvZlFUB+L2U0uFpxNdw4x3i5wnFkiSdLc/IzQPAliraBbAvbyDl0ZuLgY+TJTVvIjuf5lrgKdUkJOUDAJ9LVq7hh8ALyv08Gfgy8JMppffljW2usCK4JElTyzPCchPwwoi4OKX0X+M1iIhnky0K/uvpBJNSupescOZU7fbw8MXMlWdF4L3ln6ZiRXBJkqaWJ7l5D/Ai4PMR8Wtk0z8PKZc3+BhwGrju4V/XTHW0t7Jt0wqnoiRJmkTV01Ippe+Q7ZjaCPwbWVmDEeD5EXGAbMfTBuAtKaXv1yHWBW2oWCLt7WPHbb2kvX0MFUtTf0mSpAUoVzHMlNJfkK1p+S+yg/UKZCURVgO3A89PKTXddNBcsLt3gFvvyHZIVa6SJOnhcu9qSil9BfhKRPSQLTBuBe5NKfXWOjid4U4pSZKqM+3CmeXdS/NyS/V8ZOkFSZKqM1eqgmsK7pSSJKk6udbcqLEsvSBJ0tQcuZknKguKK0bANTeSJI3DkZt5YrwFxZIk6eFMbuYJSy9IklQdp6XmCRcUS5JUHZObecLSC5IkVWdayU1EdANPJDuZeG9KaUdNo9JZhooldvcOcLh/kJ7uTras76KjvbXRYUmSNCflSm7KSc17gF8a9d1PADvKz18D/AHwcymlb9cwzgVt9E6pPfsHKBTcKSVJ0kSqXlAcEUuBbwKvAPrIimeOPXLlK8Ba4Pm1CU/gTilJkvLIs1vqzcBjgU8B56eUnje2QUppF3An8MzahCdwp5QkSXnkmZZ6EdALvCqldGqSdvcAF84oKp3FnVKSJFUvT3JzPvDVKRIbgENAz/RD0ljulJIkqXp5pqWKwOIq2m0Ejk0vHEmSpJnJk9wk4HERMWGCExEryNblfH+mgemMoWKJtLePHbf1kvb2MVQsNTokSZLmrDzJzeeANcCfTtLmj4FzgH+YSVA6W2Ur+J79Z66SJGl8edbcvB94OfBbEXEx8Pny/c0R8ZtkC46fQTZq89GaRrnAjbcV3PU3kiSNr+qRm5TSCeA5wHeAS4E/Kz96BlnicznwPeCnU0pDtQ1zYXMruCRJ1ct1QnFKaR9waURcCTyXbAdVK3Av2aF+X0wpjdQ8ygXOreCSJFVvWrWlUkpfITuNWLPAreCSJFXPquBznEUzJUnKZ7pVwVvJDuqbcFt4Sume6QalMyyaKUlSPnmrgl8CvAN4GrBokqYjefvW+NwpJUlSPlUnIBHxVOA/OJPU9AEeuFJnPd2dZ51r404pSZIml2d05RqyxOYjwB+klA7WJySN5k4pSZLyyZPcPAnYmVL6jXoFo4dzp5QkSfnkKb9QAG6rVyCSJEm1kGfk5vvA2noFovG5FVySpHzyjNy8D3haRPx4vYLRw1k0U5KkfPLUlvp74J3Av0fEb0bEpvqFpYrxtoJLkqSJ5dkKXhr18f3A+yNiouYjKSXPuakBt4JLkpRPngSkUKe2moRbwSVJyqfq5CallGd9jmrEreCSJOVjwiJJkpqK62LmOLeCS5KUT+7kJiLagRcClwMbyrf3Ad8EPpdSKtYqOFkVXJKkvHJNS0XEE4AEfAp4FfDc8s+ryvfuiIjH1zrIhcyt4JIk5ZNnK/hG4KvASuBe4NPA3WQ7o7YALylfvxoRP55S2lf7cBcet4JLkpRPnmmp3yNLbK4Ffmfs9FNEvB34M+D15ba/VasgFzK3gkuSlE+eaakrgV3AG8dbV5NSOg28qdzmubUJT5Wt4JdetJ5tm1a4mFiSpCnkSW42ALeklEYmapBSGgZuAdbPNDBJkqTpyJPcDJJNS01lRbmtJEnSrMuT3NwGXB4Rj5qoQWTFpq4ot9UMDRVLpL197Litl7S3j6FiaeovSZK0wOVJbj4KdADfiIhfiYiOyoOIaI+IVwJfB9qBj9Q2zIWpcsbNnv1nrpIkaXJVJzcppb8FPgOsJUtejkfEPRGxFzgB/DXZWpvPpJQ+XY9gFxrPuJEkKb9ch/illH4JeB2wB2gFNgLnlX/fDbwupfTSGse4YPV0n32mjWfcSJI0tdzlF1JKHwQ+GBEbyHZQFYD7PLSv9jzjRpKk/KZdOLOczJjQ1FHljBtrSUmSVL1c01KSJElzncmNJElqKtOellL9DRVL7O4d4HD/ID3dnWxZ32X5BUmSpmByM4dVzrkB2LN/gEIB199IkjQFp6XmMM+5kSQpP5ObOcxzbiRJys9pqTnMc24kScqv6pGbiChFxEeraPeRiDg9s7AEZ865ufSi9WzbtMLFxJIkVSHPtFSh/FNtW0mSpFlXjzU35wDFOvQrSZI0pZqtuYmIFmA78Ezgvlr1K0mSlMekyU1ElMbcenlEvLyKfj85/ZAEHuAnSdJ0TTVyM3rtzAiTr6UpkhXS/ALwthnGteB5gJ8kSdMzaXKTUnpoTU5EDAMfTyn9St2j0rgH+JncSJI0tTwLiq8BvlivQHQ2D/CTJGl6ql5QnFK6pp6B6Gwe4CdJ0vRUndxERAewHBhIKZ0cdf8c4PeAxwJ7gP+TUrq3xnEuOJUD/JyKkiQpnzzTUm8D9gOPq9wob/++Efj/gJ8GXgvcHBE9tQxSkiSpWnmSm2cB+1JKN4+69wLgx4HbgV8j2ym1Hnh1zSKUJEnKIU9ysxlIY+79LNkW8ZemlD4GvIhsdOcFNYluARsqlkh7+9hxWy9pbx9DxbFHDkmSpPHkSW5WAgfG3LsU2JtS+j5ASmkY+A6wqTbhLVyVc2727D9zlSRJU8uT3BSB7sqHiFgDnA/855h2J8jqS2kGxjvnRpIkTS1PcnMn8NSIWFz+/PNkU1Jjk5t1wAM1iG1B85wbSZKmJ0/hzM8CfwzcGBH/SbaAeIhRB/tFRCvweODWWga5EHnOjSRJ05MnuXkP8GzgCuBioAS8IaU0epTmOWRTVzfWLMIFynNuJEmanjwnFJ+KiJ8ALgPOBb6XUto1ptlJ4I3AP9cuREmSpOrlGbkhpTQCfGuS59cD1083mIjYCLwDuBLoIdtW/kXgmpRSX86+nga8gWxH10rgQeD7wHtTSl+eboySJGluy7OguK4iYivZWp1XAreQTYPtAl5PzlOPI+IPyKbGng58BfgL4F+AFcDlNQ1ckiTNKblGbioiYjuwDegCCuO1SSl9Mme3HwTWAFellK4b9bfeTTbV9U6qOPk4Il4E/BHwH8DPpZSOjnnenjMuSZI0j+RKbiLiUuCvgO2TNCuQbRGvOrmJiPPJFiPvAT4w5vHbgV8HXhYRb0opHZ+knxbgXWRn7bxkbGIDkFIqVhtXowwVS+zuHeBw/yA93Z1sWd9FR3tro8OSJGleyFMV/FHA14AlwA5gLbAF+DvgkWQFNVvJ1sj054zjmeXr18qnHD8kpXQ0Im4iS36eDHx9kn4uLcf0OaAvIn4aeAzZQudbxtTFmrMqpxMD7Nk/QKGAu6YkSapSnjU3v0eW2PxGSukyyguLU0q/lFK6BHgs2ZqZbcBVOeOI8vXOCZ7fVb5um6KfJ5avB4DvAf8K/CnwXmBHRNwQEatzxjbrPJ1YkqTpyzMtdTlwV0rpI+M9TCntjIjnAT8C3ga8JUfflbIOE434VO4vn6KfNeXrq4HdwE+Q1bp6BNmi4p8kO4zw8rFfjIiryabAJrRz584p/nw+J0+efFifLS0ttAy3c/DgwYfubT23g5QSw8PDY7vQJMZ7v6od3299+X7ry/dbX41+v3mSm7XAl0Z9LgFExKKU0imAlNIDEXEDWVXwPMnNVCqLlkemaFdZmFIAXphS+r/lzz+IiBeQjQw9IyKeMnaKKqV0NXD1uH+8UBgB2L59sqVG+e3cuXPcPoeKJVauXHHW6cSuuclvover2vD91pfvt758v/XV6PebZ1rqGGfvjKqUqV43pt0gsCFnHJWRme4JnneNaTeRylk4u0YlNgCklAaBr5Y/PilnfLOqcjrxpRetZ9umFSY2kiTlkCe5uQ84b9TnO8rXKyo3ytusLwEOkk8qXydaU3NB+TrRmpyx/RyZ4Hkl+bEKpSRJTSpPcnMTcGFEVEZRvkQ2NfWeiPjNiPhfwD8CG3l4pfCpVE41fk55O/dDImIZ8FSyEaFvT9HPjcBp4IKI6Bjn+WPK1z0545MkSfNEnuTm88A+yotxU0r7gD8hmzJ6P9kW8OeRTR29NU8QKaW7yb37xDwAACAASURBVLaZbwZeO+bxNcBS4JOVM24ioj0iHlU+1Xh0P4eAvyeb3vrD0c8i4tlkC4r7yU4tliRJTShP4cyvc2Z6qHLv7RFxG/BCsvpNd5DVbtozjVheQ3Z+zrUR8SxgJ9kU1xVk01GjE6YN5ed7yRKi0X67/L23RsTTyUo5PIJskXMJeFVKaaJpK0mSNM9Nq/zCaCmlfySbjpppP3dHxMWcKZz5XLLCmdeSFc58sMp+HoiIS4A/IEtongwcJZtG+5OU0lRTWw3nCcWSJE3fjJObWkop3UtWOHOqdnuYoKZV+fmDZCM4v12z4GaRJxRLkjR9uZObiOgGXgo8BVgNfD2l9H/Kz4JsCuhb5a3XmobxTig2uZEkqTp5C2deCXya7KTgSoHMfaOaPK78/CVkC3s1DT3dnezZP/DQ51XL3bkuSVK18hTOfAzZjqk24INk267HJjD/TLZl+2fHeaYqbVnfRaHAWScUS5Kk6uQZufl9YBHwgpTSPwNExFkJTErpRETcQVZEU9NUOaHYqShJkvLLc87N5cB/VxKbSdzLw0sySJIkzYo8yU0PWcXvqQxjeQNJktQgeZKbPrLSClPZChyYXjiSJEkzkye5uQV4YkRcMFGDiHgicBFZHSpJkqRZlye5+QDQDnyufJ7NWSLifOBjZNvD/7I24S08Q8USaW8fO27rJe3tY6hYanRIkiTNK1UnNymlrwLXAT8G/LBcU2oE+ImI+A5ZXakLgfeklPJWBVdZ5XTiPfvPXCVJUvXyjNyQUno9WYHLA8BjyA7y2wg8kaza9htSSm+udZALyXinE0uSpOrlSm4AUkofIktoLgZeDPwicBmwLqV0bW3DW3h6us/eaObpxJIk5TOtwpkppWHge+Uf1ZCnE0uSNDNzqiq4PJ1YkqSZmjC5iYinz6TjlNKNM/m+JEnSdEw2cvNNst1Q0zEyRd+SJEl1MVkCciPTT24kSZIaYsLkJqV0+SzGIUmSVBO5t4JLkiTNZa6LmWOGiiV29w5wuH+Qnu5OtqzvoqO9tdFhSZI0b5jczDGV8gsAe/YPUCjgtnBJknJwWmqOsfyCJEkzY3Izx1h+QZKkmXFaao6x/IIkSTNjcjPHWH5BkqSZcVpKkiQ1laqTm4j4WET8ShXtXhERH5tZWJIkSdOTZ+TmFcBlVbR7KvDyaUUjSZI0Q/WYlmoHhuvQryRJ0pTqkdxcCBypQ7+SJElTmnS31DhrZy6bZD1NG7AdeDzwpRrEJkmSlNtUW8FfMer3EeCR5Z/J3A+8dQYxLWjWlpIkaWamSm5eWb4WgI8B/wl8dIK2Q8A+4NsppaHahLfwWFtKkqSZmTS5SSl9ovJ7RFxNlrh8YuJvaKbGqy1lciNJUvWqPqE4pbS5jnGorKe7kz37Bx76bG0pSZLysfzCHGNtKUmSZiZ3chMRi4ErgG1AF9l6nLFGUkp/NMPYFiRrS0mSNDO5kpuI+HngQ8DKSZoVyHZWmdxIkqRZV3VyExGXAH9HdvrwZ4DHAD8G/CnZ9vBnA91ku6nuq3mkkiRJVchzQvGby+1/LqX0UuC/AVJKb00p/QLZNNWXgeeSje5IkiTNujzJzaXA7SmlcU8fTikdAl4CLAKuqUFskiRJueVJblYBadTn0wAR8dBe5ZTSUeBG4KdqEp0kSVJOeZKbPrJRmYpKccyNY9qNAGtmEpQkSdJ05Ulu7gU2jfp8O9nOqOdVbkTEUuAysjIMymmoWCLt7WPHbb2kvX0MFUuNDkmSpHknz1bwbwKvj4jVKaWDwL8CJ4A/iYi1ZDukfpls+urztQ50IbCulCRJM5dn5OazwA3A4wBSSoeBNwHtZDup3gs8gSzJeVttw1wYxqsrJUmS8slTW+oWsrNsRt/7cET8F/BCsoP97gD+JqV0ZJwuNAXrSkmSNHMzri2VUroVuLUGsSx41pWSJGnmLJw5h1hXSpKkmcuz5kaSJGnOM7mRJElNxeRGkiQ1FZMbSZLUVExuJElSU5kwuYmIqyLiJ2YzGEmSpJmabOTmvcBLKh8iYldEvKv+IUmSJE3fZOfcDAOtoz5vBlbXNZoFbqhYYnfvAIf7B+np7mTL+i462lun/qIkSXrIZMnNg0DMViCycKYkSbUwWXJzM/C8iLgR+FH53mUR8bEq+h1JKf3qjKNbYMYrnGlyI0lSPpMlN78PPBa4rPwD8Mjyz1RGAJObnCycKUnSzE2Y3KSUfhARjwaeBGwCPg78J/DR2Qlt4bFwpiRJMzdp4cyU0gngmwAR8XHgRymlT9Q/rIXJwpmSJM1cnqrgVwD31ysQSZKkWqg6uUkp3VDPQCRJkmohz8gNAOV1OK8HLgc2lG/vA64Hrksp/aBm0UmSJOWUq7ZURPwa8D3g14ALgCXlnwuAXwdujQh3SUmSpIapOrmJiEuAD5GN9nwWuJIsqdkG/CTw9+VnHyq3lSRJmnV5pqXeDBSAX0wp/cOYZz8C/j0iPk+W5LwJeHFtQpQkSapenmmpy4DvjpPYPCSl9FngFuBpMw1MkiRpOvKM3KwEvl5Fux8Bj5teOAuXRTMlSaqNPMnNg1RXemFrua1ysGimJEm1kWdaagfwxIj4uYkaRMTzgUuAm2Ya2EIzXtFMSZKUX56Rm78Afhb4+4j4DPAJYDdZkczzgV8GfhEYLrdVDhbNlCSpNvKcULwjIn4LeB/wS+Wf0QrAaeC3Uko31y7EhcGimZIk1UauE4pTSn8ZETeRnVD8dM4+ofgG4NqU0m21DXFhsGimJEm1kbv8Qjl58RRiSZI0J+UqvyBJkjTXmdxIkqSmYnIjSZKaismNJElqKrkXFNdTRGwE3kFWcbwH2A98EbgmpdQ3zT5fBnyy/PFVKaW/rkWskiRpbpozIzcRsRW4FXglWfHN9wC7yLad3xwRPdPo8zzgOuBYDUOVJElz2FwaufkgsAa4KqV0XeVmRLwbeCPwTuDV1XYWEQXgb4DDwOeBN9c0WkmSNCdVPXITEXUb5YmI84HnAHuAD4x5/HbgOPCyiFiao9urgGeSjQQdr0GYdTVULJH29rHjtl7S3j6GiqVGhyRJ0ryUJ2HZGxF/EBFr6xDHM8vXr6WUhkc/SCkdJSvEuQR4cjWdRcR24E+B96WUbqxloPVSqQq+Z/+ZqyRJyi9PcrMBuIYsyflMRFxWwziifL1zgud3la/bpuwoog34W+Ae4PdnHtrssCq4JEm1kWfNzUXA64CXAL8AvDgibidbK/OplNJMpn66y9f+CZ5X7i+voq8/BB4HXJZSqjpDiIiryabAJrRz585qu6vKyZMn2blzJy0tLbQMt3Pw4MGHnm09t4OUEsPDw5P0oMlU3q/qw/dbX77f+vL91lej32+equC3A6+OiN8BXgH8JvBjZMnNuyLiE8BfppTuqEOchfJ1ZLJGEfEkstGav8hbmTyldDVw9bh/vFAYAdi+fXueLqe0c+fOh/ocKpZYuXLFWVXBO9pba/r3FprR71e15/utL99vffl+66vR73c6hTOPkm2vvi4irgBeC/wM8FvA6yLim8D7gX8au35mEpWRme4JnneNafcwo6aj7gTeVuXfnTOsCi5JUm3MaCt4Sul64PqIWAf8HlmCc3n5pzci3gd8oIrpoVS+TrSm5oLydaI1OQDnjPr+yYgYr81HIuIjZAuN3zBFTJIkaR6a8Tk3EfF44DXA/y7fGgbuAC4E3kU2lfVTKaW7JugC4Pry9TkR0TJ6xCcilgFPBQaBb0/SxyngoxM8ezzZOpz/JEukck1ZSZKk+WNayU1EdJAlM68Bnki2JuYQcC3wwZTSfRHxWLJSCv+L7LTh503UX0rp7oj4GtlZN68lm/aquAZYCny4smg5ItqBrUAxpXR3uY9B4NcmiPdqsuTmE5ZfkCSpueVKbsqH7b2a7GC8lWRJzX+TJSOfSSmdqrRNKf1f4Gcj4jtANdvGXwPsAK6NiGcBO4FLgCvIpqPeOqrthvLzvcDmPP8ZJElSc6s6uYmIL5ONrLQAp4HPAtemlHZM8dXbgYun6r88enMxZwpnPpescOa1ZIUzH6w2VkmStHDlGbm5EngA+CuyLd/7q/zeF8kO1JtSSuleslGhqdrt4cz28Gr6vZoJtnlLkqTmkie5eTnwdymlYp4/kFL6F+BfckUlSZI0TblqSwFbpmoUERdExNOnH5IkSdL05Ulurgd+t4p2b+HM1m5JkqRZlWdaqkCOdS6q3lCxxO7eAQ73D9LT3cmW9ZZekCRpumZ8iN841pAduKcq7e4d4NY7DgCwZ/8AhQKWYZAkaZomTW7GWTuzdpL1NG3AdrLt4pZazeFw/9m54KEjgyY3kiRN01QjN9/k7ErcP1n+mUwB+PAMYlpwero72bN/4KHPq5Z3NjAaSZLmt6mSmxs5k9w8g+ycmzsmaDsE7AO+UN7+rSptWd9FoZCN2Kxa3snmdV1Tf0mSJI1r0uQmpXR55feIGAb+LaX0K/UOaqHpaG9l26YVTkVJklQDeRYUXwHcX69AJEmSaqHq5CaldEM9A5EkSaqFPIf4SZIkzXkTjtxERIlsMfGjU0p3lj9XaySlVI8zdCRJkiY1WQIy9kTiPKcTe5KxJElqiAmTm5RSy2SfJUmS5iITFkmS1FRMbiRJUlMxuZEkSU1lst1S35hBvyMppWfN4PuSJEnTMtluqctn0O/I1E1UMVQssbt3gMP9g/R0d7JlfRcd7a2NDkuSpHlpsuTmilmLYoHb3TvArXccAGDP/gEKBawzJUnSNE22FdxyC7PkcP/gWZ8PHRk0uZEkaZpcUDwH9HR3nvV51fLOCVpKkqSpWCJhDtiyvotCIRuxWbW8k83ruhodkiRJ89Zku6X+sPzr+1NKD476XI2RlNIfzSy0haOjvZVtm1Y4FSVJUg1MNnJzNdmup78DHhz1ebK6UZXnI4DJjSRJmnWTJTfvIEtSDo35LEmSNGdNtlvq6sk+S5IkzUXulpIkSU1l2rulIqIF6CGbqnowpTRcs6gkSZKmKffITURcGRFfBY4C9wMHgKMR8dWIeG6tA5QkScojV3ITEe8FvgQ8G+gkG7UZKf/+bOBfIuLaWgcpSZJUraqTm4h4BXAVcIxs59QFZElNZ/n3a8hGc14bEa+seaSSJElVyLPm5nXAaeAnUkrfHfPsbuCaiPgycBPwGuBvahOiJElS9fJMS20HbhgnsXlI+dkNwKNnGpgkSdJ05Bm5OQE8UEW7g8DglK0EwFCxxO7eAQ73D9LT3cmW9V10tLc2OixJkuatPCM3NwFPjIgJyy+Un11cbqsptLS0sLt3gFvvOMCe/WeukiRp+vIkN28HNgJ/ERHtYx9GRBvw5+U2b69NeM2ttbWVw/1nD3IdOuKglyRJMzFZVfBfHuf2x4HXAy+KiM8Cu8v3NwMvAjYAHwIuAv6nloE2o1KpRE9351mjNauWdzYwIkmS5r/J1tx8nPELZRbIkpjXj3Mf4NXln0/ONLhmNzw8zNb1XRQK2YjNquWdbF7X1eiwJEma1yZLbj6JVcDrrqO9lW2bVrBt04pGhyJJUlOYrCr4K2YxDkmSpJqwKrgkSWoqJjeSJKmp5DnE7yERsR3YBnRxZiHxWVJKLiiWJEmzLldyExGXAn9FVophIgWyhcgmN5IkadZVndxExKOArwFLgB3AWmAL8HfAI4HHAa3AF4H+mkcqSZJUhTxrbn6PLLH5jZTSZcC3AFJKv5RSugR4LHAr2XTVVbUOVJIkqRp5kpvLgbtSSh8Z72FKaSfwPGAT8LaZhyZJkpRfnuRmLXD7qM8lgIhYVLmRUnoAuAF4QU2ikyRJyilPcnOMs3dGVQoirRvTbpCsPIMkSdKsy5Pc3AecN+rzHeXrFZUb5WrhlwAHZx6aJElSfnm2gt8EvDIiulJKA8CXyKam3hMRi8mSn1cBG8l2UEmSJM26PCM3nwf2kS0sJqW0D/gTsoP83k+2Bfx5ZNvA31rTKJvYULFE2tvHjtt6SXv7GCqWGh2SJEnzWtUjNymlrwMXjLn39oi4DXghsJJsquq9KaU9tQyyWbW0tLC7d4Bb7zgAwJ79AxQKWCFckqQZmFb5hdFSSv8I/GMNYllwWltbeaB/8Kx7h44MmtxIkjQDFs5soFKpRE9351n3Vi3vnKC1JEmqxnQLZz6FbO3NRrI6Ur3AN1NKO2oXWvMbHh5m6/ouCoVsxGbV8k42r+tqdFiSJM1reQtnXgD8LfDE8q3KuTcj5ef/BfxySinVLMIm19HeyrZNK5yKkiSpRvIUzjwPuBE4l2xH1L8Ce8qPN5PtlHoicENEPCmldE9NI5UkSapCnpGbPyJLbP4W+K3yWTcPiYgu4Frgl4F3AK+oUYySJElVy7Og+ErgXuBXxyY2AOV7vwbcU24rSZI06/IkN93ATSml0xM1KD+7mexgP0mSpFmXJ7nZDVSz6rUb2Du9cCRJkmYmT3LzSeDyiIiJGkTEo4Bnkq3LkSRJmnV5kps/IyuW+c2I+M3yAmIAImJZRLwa+AbZLqo/rW2YkiRJ1Zlwt1RE7Jrg0blkhTLfHxFHyveWj3r+BOAuYGtNIpQkScphsq3gmyd5Vjm8b7w1OI+gfKifJEnSbJssudkya1FIkiTVyITJTUrJHU+SJGnesSq4JElqKjOtCr6hfGsfWVXwm2sU14IxVCyxu3eAw/2D9HR3smV9Fx3trY0OS5KkeStvVfDNwKeBJ5dvja0KfjPw0pTSnhrF19RaWlrY3TvArXccAGDP/gEKBawQLknSDFQ9LRURK4DrgacAJ4C/B94J/DHwd8Bx4FLgG+W2mkJrayuH+wfPunfoyOAErSVJUjXyjNy8hWyb9z8Cr04pHR79MCJWAh8CXgj8DvD7tQqyWZVKJXq6O9mz/0wd0lXLOxsYkSRJ81+eBcU/C+wnm3Y6PPZhSulB4GXlNs+vTXjNbXh4mC3ru7h4+7lsXnfmKkmSpi9PcrMZ+FZK6dREDcrPvsXkBwBqlI72VrZtWsGlF61n26YVLiaWJGmG8iQ3RWBJFe06y20lSZJmXZ7kZidwRUSsnahB+dkzy20lSZJmXZ7k5lPAUuA/IuKZYx9GxBXA18hGd/62NuFJkiTlk2e31IeAnweeAfx7RPQCu8nOuNlCdqBfgWy7+IemE0xEbATeAVwJ9JAtTv4icE1Kqa+K7/cALwB+GvixckxDwPeBvwH+JqU0PJ3YJEnS/FD1yE1K6TRZ0vHnZGfabAAuA54GbCzf+3PguSmlUt5AImIrcCvwSuAW4D3ALuD1wM3lxGUqLwI+AlwCfAd4L9nW9ccAfw38Q0QUJv66JEma73KdUFzeDfWWiPhD4AmcGa25D7g1pXRyBrF8EFgDXJVSuq5yMyLeDbyR7MDAV0/Rx53AzwBfGj1CExG/T5Yw/Tzwc2QJjyRJakJVJzcR8TNAMaX0b+Uk5qZaBRER5wPPAfYAHxjz+O3ArwMvi4g3pZSOT9RPSukbE9y/PyI+RJYgXY7JjSRJTSvPguIvAG+oUxyVBcpfG7smJqV0lCyRWsKZmlbTUdmefnoGfUiSpDkuz7TUg8ChOsUR5eudEzy/i2xkZxvw9dydR7QBv1z++JUJ2lxNNko0oZ07a7vD/eTJkzXvU2f4fuvL91tfvt/68v3WV6Pfb57k5hayhbn10F2+9k/wvHJ/+TT7/1Oy2L+cUvrqeA1SSlcDV4/3rFAojABs3759mn9+fDt37qx5nzrD91tfvt/68v3Wl++3vhr9fvNMS70LuDAifrVewUyissNpJO8XI+Iq4E3AHWS1ryRJUhPLtVuK7Pyav4qIF5KtwdkLDI7XMKV0Y45+KyMz3RM87xrTrioR8VrgfcAPgWeVi3tKkqQmlie5+SbZyEkB+EmyNTATGcnZdypft03w/ILydaI1OQ8TEW8gOyvndrLE5oEc8UiSpHkqTwJyI9OYFqrS9eXrcyKiZcwZNcuAp5KNEH27ms4i4nfJ1tn8D/DslFK9FkLPSGtbB2lvH4f7B+np7mTL+i6rgkuSNENVJzcppcvrFURK6e6I+BrZaNBrgetGPb6GrKbVhytn3EREO7CV7Nydu0f3FRFvIyvhcCvwnLk8FfXg8QK7dh0AYM/+AQoF2LZpRYOjkiRpfsu75qaeXgPsAK6NiGeRVRa/BLiCbDrqraPabig/3wtsrtyMiJeTJTYl4FvAVRHBGHtSSh+vy3+CnA4fOcHo/woOHRk0uZEkaYamTG4iYjnZGptHAKeA/0kp3VDrQMqjNxdzpnDmc8kKZ15LVjizmhGYLeVrKxMfOHgD8PGZRVsbPcuX0H9g6KHPq5Z3NjAaSZKaw6TJTUS8GPgrYNmY+/8NvCCldG8tgyn398oq2u3hzPbw0fevZoKzauailUtHWLn9XA4dGWTV8k42r+ua+kuSJGlSEyY3EfFY4FPlNsfJpoa6yEZHHk9Wn+lJsxBj0yqdHmL7BSucipIkqYYmO8Tvt8kSm08B61JKT0gpXQBcDOwGnhARl9c/REmSpOpNltw8jWzNy6tSSscqN1NK/wO8kWxa6Gn1DU+SJCmfyZKbdcB3U0qnxnlWOX14fe1DkiRJmr7JkptFZJXAHyaldGRUG0mSpDkjT+FMSZKkOW+qc27WRsTTp/M8Z+FMSZKkmpgqufnJ8s94RiZ5nrdwpiRJUk1MloDcQ/0KZUqSJNXFhMlNSmnzLMYhSZJUEy4oliRJTcXkRpIkNRWTG0mS1FRMbiRJUlMxuZEkSU3Fs2gaqLWtg7S3j8P9g/R0d7JlfRcd7a2NDkuSpHnN5KaBHjxeYNeuAwDs2T9AoQDbNq1ocFSSJM1vTks10OEjJ876fOjIYIMikSSpeZjcNFDP8iVnfV61vLNBkUiS1DyclmqglUtHWLn9XA4dGWTV8k42r+tqdEiSJM17JjcNVDo9xPYLVrjORpKkGnJaSpIkNRWTG0mS1FRMbiRJUlMxuZEkSU3F5EaSJDUVkxtJktRUTG4kSVJTMbmRJElNxeRGkiQ1FZMbSZLUVExuJElSUzG5kSRJTcXkRpIkNRWTG0mS1FRMbiRJUlMxuWmglhZfvyRJteb/ujbIULFE32A7O27rJe3tY6hYanRIkiQ1hbZGB7BQ7e4d4Du372P16tXs2T9AoQDbNq1odFiSJM17jtw0yOH+wbM+HzoyOEFLSZKUh8lNg/R0d571edXyzglaSpKkPJyWapAt67t48mM2UGpZzKrlnWxe19XokCRJagomNw3S0d7K8s4iEec3OhRJkpqK01INNDw83OgQJElqOiY3kiSpqZjcSJKkpmJyI0mSmorJjSRJaiomN5IkqamY3EiSpKZiciNJkpqKyY0kSWoqJjeSJKmpmNxIkqSmYnIjSZKaismNJElqKiY3kiSpqZjcSJKkpmJy00AtLb5+SZJqzf91bZChYom+wXZ23NZL2tvHULHU6JAkSWoKbY0OYKHa3TvAd27fx+rVq9mzf4BCAbZtWtHosCRJmvccuWmQw/2DZ30+dGRwgpaSJCkPk5sG6enuPOvzquWdE7SUJEl5OC3VIFvWd/Hkx2yg1LKYVcs72byuq9EhSZLUFExuGqSjvZXlnUUizm90KJIkNRWnpRpoeHi40SFIktR0TG4kSVJTMbmRJElNxeRGkiQ1FZMbSZLUVExuJElSUzG5kSRJTcXkRpIkNRWTG0mS1FRMbiRJUlMxuZEkSU3F5EaSJDUVkxtJktRUTG4kSVJTaWt0AKNFxEbgHcCVQA+wH/gicE1KqW+2+5EkSfPPnBm5iYitwK3AK4FbgPcAu4DXAzdHRM9s9iNJkuanuTRy80FgDXBVSum6ys2IeDfwRuCdwKtnsR9JkjQPzYmRm4j/1965x1tVVXv8i/lGxUdWaA/UcphpD01STEEtyjRFE/OmFj0sX5lJReUL7OYzU0RT01uo+bk+03wDCojig9S6txQHCFJmqWnGBQQtPfePMZdnsdhrn70PZ++z2ef3/XzOZ5615lxzzjX22nuONeaYY9qWwHBgAXBRIftUYAlwuJn1b0Y9QgghhFh1aQnlBtgzpZPd/Y18hrsvAmYC6wI7N6keIYQQQqyitIpyYymdU5I/N6VbN6keIYQQQqyitIrPzYCULizJz85v2Kh6zGwsMXVVyuzZs7tovj6WLVvW43WKTiTfxiL5NhbJt7FIvo2lt+XbKspNV/RLaUej6nH3scDYihf169cBsO22265k80IIIYRoNK0yLZVZVAaU5G9QKNfoeoQQQgixitIqlhtPaZkvzPtSWuZL09P1LEdHR0e/rkvVj5l1uHtD6haSb6ORfBuL5NtYJN/G0tvybRXLzbSUDjez5fpkZusDuwJLgYeaVI8QQgghVlFaQrlx93nAZGAQcEwhexzQH7jS3ZcAmNkaZrZNikbc7XqEEEII0X60yrQUwNHAA8AFZrYXMBv4GLAHMY10Yq7s5in/T4Qi0916hBBCCNFmtITlBt60unwUmEgoI6OBrYALgF3c/aVm1iOEEEKIVZNWstzg7s8QG152VW4Bncu6u12PEEIIIdqPlrHcCCGEEEL0BFJuepdxvd2BNkfybSySb2ORfBuL5NtYelW+/To6VjborxBCCCFE6yDLjRBCCCHaCik3QgghhGgrpNwIIYQQoq1oqaXgfQEzeydwGvBpYBPgb8DNwDh3f7k3+9ZKmNkmwAHAPsD2RODG14A/AL8Efunub1S4bghwErAzsDbwFPALYIK7v17S1peIiNbbAq8DvwN+4u639fBttTRmdjhwZTo8wt0vr1BG8q0DM9sNOB4YAmwM/IN4hs939zsKZSXbOjCzfYBvEfee/ZY+CvzU3R+sUF7yzWFmBwFDgQ8DHwLWB65298OqXNNwGZrZOsD3gUOA9wD/B0wHTnX32bXenyw3TSRtF/EoEYNnFnAeMJ/4gj6YBnQRjAQuIwIxPgycD9wIbAdcDlxnZsvFOjKz/YEZwO7A3gf1PwAAEXdJREFUTcBFwJqEnK+p1IiZ/YQI+DgwtfcrQpm61cyO7embalXM7F3ABGBxlTKSbx2Y2Ul0yusu4FzgVmAjYFihrGRbB2Z2FnAbsAMh2/HAY8D+wEwzO6xQXvJdkZOAYwnl5tmuCjdDhma2FjAFOIVQasYDdxMvuo+Y2cdqvTmtlmoiZjYJGA4c5+4Tcud/CnwbuNTdj+yt/rUSZrYnsRfY7XkLjZm9g1AM3wUc5O43pvMbEG8RA4Bd3f2RdH5tYCqwC/Af7n5Nrq4hwExgHrBTZjkzs0GEEtof2CYFjWxbkpI4BdgC+DXwHQqWG8m3PsxsJHAd8cN8oLsvKuSv4e7/Sv9LtnWQfgOeBf4OfNDdX8jl7UHI7Gl33zKdk3wrkGT1F0I2Q4mNpytabpolQzP7AXA6cAPw+ey3PylWNwNPANtXstoXkeWmSZjZloRis4DQePOcCiwBDjez/k3uWkvi7lPd/dbiQ+zuzwGXpMNhuayDgE2Ba7IvXiq/jHhDATiq0EymSP44PyWYvmwXAWvRNyJdHwfsSdxr2aaykm+NmNlqwFnAK8AXiooNQKbYJCTb+ngPMXY9nFdsANx9GrCIkGeG5FsBd5/m7nPdvRYLR8NlmF6ysmu+l//td/ffAPcRU1tDa+ivlJsmsmdKJ1cYsBcRGu66xFymqE42MPw7dy6T710Vys8gBpohyexZyzV3Fsq0JWb2fuBMYLy7z6hSVPKtnSGEFewO4GUz28fMxpjZt8xslwrlJdv6mEv43w02s7fmM8xsd8J35O7cacl35WmGDLcC3g3Mcfena7ymFCk3zcNSOqckf25Kt25CX1ZZzGx14IvpMP+lKZWvu/8beJpwoM9M1f0JJ+XF7v63Ck21/eeRZHkV8Gfgh10VT6nk2zU7pfR5wg/kNkKBPB94wMzuNbO8ZUGyrQN3/wcwBng78ISZ/dzMzjCz64DJxBTrN3KXSL4rTzNk2KNjpJSb5jEgpQtL8rPzGzahL6syZxJOxXe4+6Tc+Xrlq88jnPY+Aoxy96VdlJV8a+dtKT0SWAf4BGFN2A6YRDhkXp8rL9nWibufDxxIDKhHEKtrRgLPABML01WS78rTDBn2qNyl3LQO2cofeXiXYGbHAaOBJ4HD67y8u/Jty8/DzAYT1ppzKy2b7QaSbydvSWk/wun9Hndf7O6PE6s+/gIMLZmiqoRkW8DMvkc4nU4kpjP6AzsSq0+vNrOz66hO8l15miHDutqQctM8Mq1zQEn+BoVyIoeZHUMsC3wC2COZpvPUK9+uynf1FrHKkpuOmgOcXONlkm/tZM6T8939f/IZyUKWWRwHp1SyrQMzG0Y4bN/i7ie4+3x3f8XdHyOUx2eB0WkRB0i+PUEzZNijY6SUm+bhKS2bL3xfSsvmG/ssZnY8cCHwR0Kxea5CsVL5psF8C8IBeT6Auy8hfgTXM7OBFepr589jPUJO7weWmVlH9kes3AO4LJ07Px1LvrWTyeqfJfmZ8rNOobxkWxv7pnRaMcPdXyFCRaxGTLmC5NsTNEOGPTpGSrlpHtkXcXhaKvomZrY+sCuwFHio2R1rZcxsDBEk6veEYvNCSdGpKf10hbzdiZVoD7j7qzVes3ehTDvxKvBfJX+/S2XuT8fZlJXkWzsziB/695nZmhXyt0vpgpRKtvWRrcjZtCQ/O/9aSiXflacZMpxHLG7Y2sy2qPGaUqTcNAl3n0d48g8iQlHnGUfMGV+ZNF4BmNnJhAPxo8Be7v5ileI3AC8Ch5jZR3N1rA38Zzq8uHBNFi/nRDPbKHfNIOIzepXY6qGtcPel7v61Sn/ALanYFenctelY8q2R9JxeS5jXT8nnmdkngU8RpvVstZ9kWx/3pfTrZrZ5PsPM9iZeFJcBD6TTku/K03AZpng72TVn540AKYjfboRbwr21dFgRipuIxfYLDxCrKX4DzCa2F9iDMLUNcfeXeq+HrYPFfiQTib1IJlB5nnWBu0/MXTOC+BIuI8KB/wPYj1hieANwcDFglZmdC5xAOHneQIQT/zyxV8033f3CnryvVsfMxhJTUyvsLSX51o6ZvY2IXfVeYjCeRQSfO4BwiPyCu1+fKy/Z1kga9CYRq9AWEVsBPEdMs+5LOJ4e7+7jc9dIvgWSTEakw3cQSvd8OpXHF939O4XyDZVhipMzlYgV9QhwDxH7ZiRhidvT3R+u5f6k3DQZiz18yjbOLDrJ9llyg2w17nX3YYXrdgVOJMKB5zd2u8Crb+x2LBH98g0iNsk53gab49VLNeUm5Uu+NWJmGxPRWw8gYn4sIqb7znD3FaafJdvaMbM1iLf/Q4h7X5cYbGcR8ppc4RrJN0cNv7F/cvdBhWsaLkPr3DjzC4Rik98484na7k7KjRBCCCHaDPncCCGEEKKtkHIjhBBCiLZCyo0QQggh2gopN0IIIYRoK6TcCCGEEKKtkHIjhBBCiLZCyo0QQggh2orVe7sDQvQmZraAiBwLsK+7315S7o/AB4j9raY3pXN1knZLnkaF4Ibthpl9DTgK2IYI4AawkbuXbVaJmY0iwr1f4e6jGt3HVqIvPRtCgJQbIfKcYWZ3uvsbvd0RUY6Z7QtcRoSBn0JEpoXOjRL7HGY2HRhKCyvfQjQTKTdCBK8A2wOHAlf1cl9EdUam9Dh3v6xXeyKEaEnkcyNEcEFKx5nZmr3aE9EV70rp3F7thRCiZZHlRojgRmBPYDBwJJ3KTlW6mg4ws4nAl4AvF3Ywf/M88FtiM9XdgXWA3wMnu/u0VHZf4LvAh4kdj6cDo929dHA3s/7AycDBxKaNLxAbtI4t23k+ber6HWJT13cD/wL+QEwBXVFhx983753YEG8Mscv9xsCB7n5zWf9ydawBfAM4nNjVeQ1gAfAb4Oz8ZrI5mWVMM7Ps/3HuPrar9nJ1rQ+cAhwEbAa8CNwCnFho8xRgHHCpux9ZUtd+qb+/dffB6dwokn8PIdMfAZ8F3go8C/w3cLq7v1JS5z7EZoM7ARsAzxO7JZ/p7rNz5YYRvjQZeZlAhecyyfx7hMwHERsTTgZ+6O5/LulPvc/GhqmN/YEtiOf2JWKjxUnufkZvtCH6DrLcCNHJ91N6opmt16Q2P0rsZLw1cA/gwK7AJDPbzcy+SQyc/YBJhH/JZ4EZZrZJSZ1rprqOBf4I3Ers4Hss8KCZvb14gZntQQwkxxG/C3cBDwMfpHOQLmMkMcC+m/CBuYcYmKpiZmsTg+oEYDtgRurrhoSi9JiZbZm75P7Uj+fT8aR0fAWhENbKAGAm8JV03WTCKflIYEoa/DMuJXx5DjWzDUrqOzqlF1XI24iQ48EpnQRsSuysfI+ZrVu8wMzOAG4DhgOPAzcAC4EvEjLZJ1f8OcplckXKz7MGcCfxrD+V/n+DmI69PykMxf7U9Wyke5oJ/IBQ5u4mFOuniJ2hV9iJuhltiL6FLDdCJNx9mplNJgaV0cQbe6M5hrDC/DQ7YWZnEW+klwPvAIa5+30pL1MIdiMG1R9VqHMXYA5g7v5sum594CZgL0KZODjX3kDCcrUeMAq4MntLTm/TtwCHm9nUvPUpx9HAN9z953Xe+2nAMOBJ4BO5vq4D/Ao4ELg63Q/ufjlwebIYvZ2wYkyvs02AEcAdwBB3X5za3Ax4CNiBkM3Vqc3nzew64DBCubgwX5GZbUU8Ly8B11Zoaz9iEN4xW8mVlMspwM7AWOKzzur7DKF4LAE+4+4zcnnfBc4Grjazrd39BXd/EhhVh0yGAI8AW7n7C6neAYRVaAfiefxxrs3uPBsHEQrG7cAId/93rr63ENa+vAwb3oboe8hyI8Ty/ADoAEab2aZNaO/BvGKTODOlWwMXZYoNgLsvA85Lh3tUqXd0piyk6xYRlonXgc+lQSPjeMLCcK67L2f+d/dngCPS4TdL2ppSr2KTFJij0uFxhb4uJaaqlgA7m9mu9dRdA4uBr2aKTWrzr3QqLnsVyk9I6VGsyNGEVe0X6bMp0gEclV+i7u7PA99Kh0cmhTVjdErH5xWbdN05hDVjAJ2fSb10AF/JFJtU70LgrHRYvPfuPBuZZfDuvNKRrnnd3af2QhuijyHlRogc7v4YcB2wPjF10GjuqtCHlwlLQMV8Oh1pNyup85/ufluFep8irBOrEf49GZ9J6fUl9T1KKAQfLgzEGb8uua4aOxJv6n919ykV+voiMUUFYd3pSR519+J0DYQFCQpydfdZxNThtsnHBXhTQRtFTOtcUtLW/7r7H4onkz/Vs8RztmOqb3ViShJgYkl9v0zpsJL8rvhzpf5Qcu9079mYldIxZnZYpamuXmhD9DE0LSXEipwEfI54qz7P3f/UwLb+UnJ+MbBJSX5mcaikaEA45JaxgBhA35k7l/m1/LbgjFqJTYhBOU935LN5Sp+uUmZeoWxPUdFplnCshcpyvYCYKjuacOgGOIRwnr7d3eeX1Fnt/hYQ95Z9FpsAaxHKUplMV1Ym9d573c+Gu99rZmcTzsFXAR1m9iThM3Wju0/qhTZEH0PKjRAF3P0pM7ucmMY5jeVX6NRLV9bRrgIGNiqgYH7lyVtSei0RGK8ar1Y4t7Qb7fer0I+yMj1Nd2R6PXAuMMLMBrr73+h0JP7ZSvYnk0G/CueKrKxM6r33bj0b7j7GzC4hVjJ9nFCojwCOSH5t++Smk5rRhuhjSLkRojKnEQ6kh5nZOVXKZVFxy1ZXvafkfCMZVEPeX3PnngHeC/zI3R9vUJ+KZBapLaqUyfKKlqKm4+6vmdmlxPLxI8zsLmKl23wqTx1mDKohL/ssXiQG77VSXqWl/s2WSbefDXd/Gjg//WFmHyeWwA8nVqplflrNaEP0MeRzI0QF0pv5eOI7cnqVotkgs00xI62K2aHne9clG6ZVN8X+bEms0Okgll1n3JnSkcVrGkjmR7G5mRWdWEnL3D+bDqc3sV/VuIRY4v51YskywMVdbNfxITPbrnjSzIYSU0uLCVmQrAwzU5EvltQ3KqXTC+czJbunX1h77Nlw9/vp9CX6UJPbEH0MKTdClHMWnXFlyiwM96T0mLSkFQAz25iIzdGseDlFzi30Zz3gYmIK4KZCsLZzCJ+LH5rZMcmxdTnMbGcz6zHlJ62Iypxwxxf6unbq63rAQ+4+s0IVTScpvDcSSsmhxBTKL7q4rB/ws7TcGoC0Cm98Ovx5kkVGtnLu+OIqMTM7gVgWv5AIE5AnU7LfX9vd1Ezdz4aZHWBmu5vZaoVy6wCfSId5n6JmtCH6GJqWEqIEd19oZmcSsUVWCLaWuA44AfgI8LiZzSSC6O1ETDfcTMRVaSYPEkrMHDObSrzVDyWCx80jYpm8ibs/Y2YjiGBxFxJBDB8nVmxtBmyV0mspX9HSHU4mpnaGAXNTX5cSMXwGEs6vh/Zgez3BBYQjMcA1+WjGJdxCBCicl2LRrE4s4d+AiEx9Sr6wu9+e4hyNIQI13kc8R9unepYBh6Xl5HluIqw655jZJ4mI1ADnuLvXe5O5/nTn2RhKLHX/u5n9Dvg7sXx9COGA/SQRHLFpbYi+hyw3QlRnAuUrmnD314g3xYuJgflTxBTVFcQP7cIm9LHIa8RWEpcSEV73S+cuAnautAw6LU3+ADEF9wIxfTWCiDo8l4j/06NL41NcmOHEFM8TxKC/P/EWfzawQ5VVSL3FLDpXFtXiSPwyIcubCKvL3sSgfTqxNcKS4gXu/n3CWjiFUGoOIuLAXEUEA6y0zP8WwsH5SeJ5/Gr6G1gsWy/deDYmElbPOYRCNpLY1uQp4NvA4BRbp6ltiL5Fv46OaosVhBBCZJjZ/oQ1bpa7f6xKuVGkbQPcfVRzeieEyJDlRgghaiD5gpycDotRpYUQLYR8boQQogpm9mUiovNgYj+jhwlfKyFEiyLLjRBCVGco4aw7kFBqDsjvfySEaD3kcyOEEEKItkKWGyGEEEK0FVJuhBBCCNFWSLkRQgghRFsh5UYIIYQQbYWUGyGEEEK0Ff8Pnz/1SFsOz9IAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(8,8))\n", "sns.scatterplot(number_of_tests, p_errors, alpha=0.5, ax=ax)\n", "ax.set_xlabel('Number of hypotheses')\n", "ax.set_ylabel('Probability of at least one error')\n", "ax.grid()\n", "fig.tight_layout()\n", "fig.savefig('p_errors_FDR.pdf')" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAI4CAYAAABndZP2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzde3Sk913n+XepqlSqUukutSW1Wi2p1fo5kp1NHAcSB5zEAySZ2VwAQwwTYICFCUsCbBjCLmTAZ4AZGBwuIRxyyCxkhjBrwFlgkyGBQC42SQgTu20nffnpfu+bpJJUKlXV89Qj7R9PlSzLUrdKrUuV9Hmd0+dxP8/vefSTz8nxJ7/bN7C+vo6IiIjIcVJx1B0QERER2W8KOCIiInLsKOCIiIjIsaOAIyIiIsdO6Kg7cNIEAgGt6hYREdmD9fX1wG7bagRHREREjh2N4BwRbc8XERHZnUBg1wM3GzSCIyIiIseOAo6IiIgcOwo4IiIicuwo4IiIiMixo4AjIiIix44CjoiIiBw7CjgiIiJy7CjgiIiIyLGjgCMiIiLHjgKOiIiIHDsKOCIiInLsKOCIiIjIsaOAIyIiIseOAo6IiIgcOwo4IiIicuyEjroDAMaYh4HXA68A/hegBvhTa+27bvHOA8AHgNcAVcAw8EfA71lrvR3e+SHgJ4F+wAMuAI9Zaz+1Q/so8H8CjwBngWXgC8AvW2svF/2LioiIyKEolRGcDwDvwQ84M7drbIx5O/Ak8CDwl8DvA5XAbwOP7/DOY8DHgDbgo8DHgXuBTxpj3rNN+wjwWeCX8IPN7wJ/D3wn8DVjzDcX8wuKiIjI4Qmsr68fdR8wxrwRmMYfhXk98Hl2GMExxtTm29UBr7PWfi1/vwr4HPBa4PustY9veucB4EvACPBqa20if78LeBqoBu621o5veuf/Av4j8ATwTmvtWv7+24G/Ai4B9xbu71YgEFgHKIV/7yIiIuUgEAgAsL6+HtjtOyUxgmOt/by1dshau5v/6j8MtACPF8JN/hsZ/JEggJ/Y8s6789dfK4Sb/Dvj+KM/EeCHC/eNMYFN77x/c4ix1v418BT+NNfrd9FfEREROWQlEXCK9FD++pltnj0JrAIP5KeYdvPOp7e0ATgHdAKD1tqxXb4jIiJyojmuh51I8OXnZ7ETCRx32yWxh6IcA47JXwe3PrDW5oAx/MXTPQDGmGrgNLBirb26zfeG8te+3fyMW7wjIiJyYjmux7ODN/nE54f4u69O8JmvjDEyvXhk/SmJXVRFqstfl3Z4Xrhfv8f2e31ngzHmUeCXd3hXRETkWEmlXQYnFrg4NsfzQzcJhypoqK1idHaZl3U3HUmfyjHg3E5hAVKxq3iLaX/Ln2GtfRR4dNsX84uMRUREjoNU2uXi6BwT15YJVVSQdnK4XgWxbI54LHxk/SrHgFMYPanb4Xntlna3a7/daE2xP0NEROTEcVyPoakEQ5OLBIMBYJ03vuoMieUMd3c10tFSfWR9K8c1ODZ/fcn6F2NMCOgGcsAogLU2hX+2TtwY07bN987nr5vX2+z4M27xjoiIyInhuB6jM4tcHl8gFg3jev6C4ramau7tbcZ0NnDmrtrbfOXglGPA+Vz++uZtnj0IxIAvW2uzu3znLVvagH9eziTQZ4zp3uU7IiIiJ4K/W2qB4akl4tEwQ1MJWupjxKOV1MUrecX5U9xzrpnKcPDI+liOAecJYA54xBhzf+Fm/qC/X83/9Q+2vPOR/PUXjTENm97pwi/dkAX+uHA/fx5P4Z3/bIyp2PTO24FvxT/o74v78PuIiIiUDcf1GJxcYGR6kZWMw9r6Gt3tdUxfX6a9uZp7ehvpPVN/pOEGSuck43cA78j/tRV4E/4U01P5e3PW2n+3pf0TQAa/NMMC8Db87d1PAN+79dBAY8wHgffhn5j8BH5ph3cCTcB7rbUf3tI+gj9C8wDwNeAf8M/G+R7AAR6y1n612N9VJxmLiEi5clyPockE41eXSa46uDl/Wqo6WklVJEh7c5z+7qZ9Dzd7Ocm4VALOo9x6W/WEtbZryzuvA34RvzTD5mKbH7pNsc334J9CvAY8A/zmLoptfj9+uNlcbPPS7n67F1PAERGRcjU8tcjobILUao5k2iEAhENBIpVBzp+pp6+z8UBGbso24JwkCjgiIlKOVjMuT12YJpVxiVWFuZFYpTIUJBSq4NzpegZ69n/kpmAvAacct4mLiIjIIUqlXYanEqxmPVIZl2AwQEt9jHTWpau1lr7OhiNfc7NVOS4yFhERkUOyvJLFji/w9JXrRCN+iEln1sg4OUxnI/f0tlAdPboD/XaigCMiIiLb8kduFnnGXqOmOsJzw3PURCPkvDVOn4rT1V5XciM3BQo4IiIi8hKptMv47CIXR+eoq4kyOrtET3sdC8tpBnoaGehuLsmRmwIFHBEREXmRQvHMZ+wNmhpijOXDzVIyzSv6TtHX2VjS4Qa0i+rQaReViIiUMsf1uDQ2z4Ur14lXR0iuOrQ1VTOXWGXgXDO9HfXUxiOH2idtEy8DCjgiIlKqCtNSz4/MEQoGGZtdoru9jkzW5T5ziq72+iMZudlLwNEUlYiIiOC4HlfG5hmaXiQerWQ0H26WkmkGepqPLNzslQKOiIjICVeoDD5xPUkm67G2tkZPWx2LyTT93c001leWVbgBBRwREZETbXNl8GgkhOt5LKccKirgdEsNsWiY9qbao+5m0XSSsYiIyAk2eS3JyMwSWSdHVSTon1CcyVEZDtLZWnNg9aUOmgKOiIjICeS4HjcSKwxNJ8g6OZych5PzqK6CeCxER0sNvWdKrwTDbingiIiInDCO63FxdJ6skyOTyeG4HgDhUIhAANqb42UdbkBrcERERE6cyWtJrs6v8I2Rm0SrQhCAcChIMAhdbXVlOy21mQKOiIjICbKacRmZTpDJ5qirifLc8BzxaATPW+dsW2lWBt8LHfR3yHTQn4iIHBW/eGaC0dllVtJZaqsrqQgEWUln6e9uorejoSS3g+ugPxEREdnW8koWO77A01euE434IzTpzBo5z+Pl51p4WVdTSYabvVLAEREROeb8kZtFnrHXqKmO8NzwHDXRCDlvjdOn4nS11x2LaanNFHBERESOsUJ9qYujc9TVRBnNVwZfWE4z0NPIQHfzsRq5KVDAEREROaYc12NoKsEFe4OmhthGuFlKpnlF3yn6OhuPZbgBLTI+dFpkLCIih6FQX+rC4E1CwQqSqw6tTdXMJ1YZONdMb0c9tfHIUXdzV/ayyFgB55Ap4IiIyEFLpV0GJxaYvrkCwOXxBbrb61hNZ3nVy1rpLrPK4NpFJSIicsI5rseVsXlGZhdJpV3W1tbozk9LnT/TSLw6VFbhZq8UcERERI6JwrTUxPUkudw6bq5QGTxA+6nyrQy+F6pFJSIicgw4roedWGDiapJoJEQq4wD5EgyBAJ13lW9l8L1QwBERETkGJq8lGZlZIuvkqIoEaamPkc7kqAwHOdt2ssINKOCIiIiUNcf1uJFYYWg6QdbJ4eQ8nJxHdRXEYyE6WmrKvjL4XijgiIiIlLGxmWUWk2kymRxZ1yMAhEMhAgFob46fyHADWmQsIiJStlYzLjM3k3x95CbRKn/MIhwKEgxCV1vdiZuW2kwBR0REpAyl0i5DkwlSmRx1NVG/vlQsQs5b52xbLX2dJ3PkpkAH/R0yHfQnIiJ3qnCQ34XBG7S3xMm6OSoCQZZXMvSfa6bvTMOxOutGB/2JiIgcc6m0y+CkH25qqiM8PzzH+nqA1YzLQO/xCzd7pYAjIiJSJlJpl6GpBJdG56mpjjA6u0R3ex0LS2n6Ous536FwU6CAIyIiUgYK01J2IkFtPHKiKoPvhQKOiIhIiXNcj6GpBM8O3qAyHGRtfW0j3PT3NNNYX6lws4UCjoiISAkr1Je6PL5AvDrC0FSCSDhENBKip6OBptrYiakvVQztojpk2kUlIiLFGJ5a5MrEAuvr61weX9ioDH5/fxt9nQ3Eqo7/yM1edlHpJGMREZEStLkEQyrtUBUJboSb/u5m6mvDJyLc7JUCjoiISIlxXI+Lo/NknRyZTA7H9XBcj1hVJW0tcWLRsKalbkNrcERERErM5LUkV+dX+EahBEPAL8Gwtr5Ge3MN/d1NJ/qU4t1QwBERESkhqxmXkekEmewLJRji0QieSjAURYuMD5kWGYuIyE5SaZfhqQSjs8uspLPUVldSEQiyks7S391E7wk9yE+lGkRERMrU8koWO77A01euE434IzTpzBo5z+Pl51p4WVfTiQw3e6WAIyIicsT8kZtFnrHXqKmO+JXBoxFy3hqnT8Xpaq/TtFSRFHBERESOUCrtMj67yMXROepqohslGBaW0wz0NDLQ3ayRmz1QwBERETkihRIMF+wNmhpiqi+1j7TI+JBpkbGIiMALJRguDN4kFKwguerQ2lTNfGKVgXPN9HbUUxuPHHU3S8JeFhkr4BwyBRwRESlUBp++uQKwUYJhNZ3lVS9rpbu9XiM3m2gXlYiISIlzXI8rY/OMzC6SSrusra1tlGA4f6aReHVI4WYfKOCIiIgcksK01MT1JLncOm7OYznlUFERoP1UjUow7CPVohIRETkEjuthJxaYuJokGgmRyjiAX4IhGAjQeVcNfZ2N2g6+TxRwREREDsHktSQjM0tknRxVkSAt9THSmRyV4SBn2xRu9psCjoiIyAFyXI8biRWGphNknRxOzsPJeVRXQTwWoqOlht4zqi+13xRwREREDtDYzDKLyTSZTI6s6xEAwqEQgQC0N8cVbg6IFhmLiIgckNWMy8zNJF8fuUm0yh9TCIeCBIPQ1VanaakDpIAjIiJyAFJpl6HJBKlMjrqaqF9fKhYh561ztq2Wvk6N3BwkHfR3yHTQn4jI8Vc4yO/C4A3aW+Jk3RwVgSDLKxn6zzXTd6ZBZ90UQQf9iYiIHLFU2mVw0g83NdURnh+eY309wGrGZaBX4eawKOCIiIjsk1TaZWgqwaXReWqqI4zOLtHdXsfCUpq+znrOdyjcHBYFHBERkX1QmJayEwlq4xFVBj9iCjgiIiJ3yHE9hqYSPDt4g8pwkLX1tY1w09/TTGN9pcLNIVPAERERuQOF+lKXxxeIV0cYmkoQCYeIRkL0dDTQVBtTfakjoF1Uh0y7qEREjo/CtNT0zRUALo8vbFQGv7+/jb7OBmJVGrm5U3vZRaWTjEVERPbAcT2ujM0zdm2JXG6daCS4EW76u5uprw0r3BwhTVGJiIgUqTAtNXE9SS63jpvzWE45VFQEaD9VQywa1rTUEVPAERERKdLktSTD00tEIyFczwPyJRgCATrvqqG/u0mnFB8xTVGJiIgUYTXjMjKdIJV2qIoEaamPkc7kqAwHOdtWo/pSJUIBR0REZBcc12N2bpnEsstq1sPJ+X9iEYjHQpw5VaPK4CVEU1QiIiK7MDazzOyNVZ6+fJVoxA8x4WCItfU1TjVUK9yUGAUcERGR21jNuMzcTHJpbG6jMng8GiHnrakyeIlSwBEREbmFVNplaDJBKpOjria6UYIhsZxmoKeRge5mnVJcgnTQ3yHTQX8iIuWjcJDfhcEbtLfEybo5KgJBllcy9J9TZfDDspeD/jSCIyIiso1U2mVw0g83NdURnh+eY309wGrGZaBX4abUKeCIiIhskUq7DE0luDQ6T021Xxm8u72OhaU0fZ31nO9QuCl1CjgiIiKbFCqDD04kqI1HNtbcLCXTvKLvFH2djQo3ZUABR0REJG9zZfBwOMja+ho9bfn6Uj3NNNZXKtyUibI+6M8Y86+Anwb6gSbgKvA08FvW2q9s0/4B4APAa4AqYBj4I+D3rLXeDj/jh4CfzP8MD7gAPGat/dS+/0IiInJkHNfDTiwwcTVJPBrm8vgCL+9tpoI1Gk830FQbU32pMlK2u6iMMb8BvB+YB/4KmAN6gbfhB7cftNZ+fFP7twOfADLAnwELwFsBAzxhrf2ebX7GY8DPAtPAE0Al8AjQCLzXWvvhYvutXVQiIqXHcT0GJxcYnlok43pEI0Ecd52lZJp7z7VwurWajpa6o+7mibWXXVRlGXCMMa3ADHATeLm19samZ28EPgeMWWt78vdq8Udr6oDXWWu/lr9flW/7WuD7rLWPb/rOA8CXgBHg1dbaRP5+F/4oUTVwt7V2vJi+K+CIiJQWx/UYmkwwfnWZ5KqDm/MH9KujlVRFgrQ3x1U884idpG3iZ/H7/tXN4QbAWvt5IAm0bLr9cP7vjxfCTb5tBn/KCuAntvyMd+evv1YIN/l3xoHfByLAD9/xbyIiIkemMC01czOJ43ob4SYcChIIoMrgZaxcA84Q4ADfZIxp3vzAGPMgUAP8/abbD+Wvn9nmW08Cq8ADxpjILt/59JY2IiJShiavJRmZWWIl7RKtCkHADzehUAVdbXWqDF7GynKRsbV2wRjz88BvAZeMMX+FvxbnHP4anM8C/3bTKyZ/HdzmWzljzBgwAPQAl40x1cBpYMVae3WbLgzlr3378fuIiMjhW824jEwnyDo5nJxHMBigpS7GajZHV2uN6kuVubIMOADW2t8xxozj74L6sU2PhoGPbZm6KqwMW9rhc4X79Xts/yLGmEeBX97hXREROWKptMvwVILVrIeTn5ZKZ9aoiK4x0NNId3u9wk2ZK9uAY4x5P/AfgQ8BHwauAXcD/wn4U2PMK6y179/l5wqLlopd+btte2vto8Cj2/6g/CJjERE5GhslGKxfX2olDeFgiGAQOk7VKNwcE2W5BscY8wbgN4D/z1r7PmvtqLV21Vr7DPCd+DusftYY05N/pTDistMev9ot7W7X/nYjPCIiUoK2lmB4bniOeDRCzlvjbFutpqWOkbIMOMD/mr9+fusDa+0q8M/4v9srC7fz15esmTHGhIBuIAeM5r+Rwg9JcWNM2zY//3z++pI1PSIiUpoKlcHtlhIM84sp+rsbGehu1inFx0i5BpzCbqeWHZ4X7jv56+fy1zdv0/ZBIAZ82Vqb3XT/Vu+8ZUsbEREpYYX6Us8O3qCyUIKhXSUYjrNyDThP5a8/bow5vfmBMeYtwOvwTyz+cv72E/gnHT9ijLl/U9sq4Ffzf/2DLT/jI/nrLxpjGja904VfuiEL/PEd/yYiInKgNteXildHGJpKEAmHiEZC9HSoBMNxVa4nGVcAfwt8G/6hfn+Jv8j4ZfjTVwHgZ6y1v7vpnXfgB50M8Dh+qYa3kS/VAHyvtfZF/zKMMR8E3seLSzW8E7/ulUo1iIiUuMK01PTNFQAujy/QnR+5ub+/jb7OBmJVGrkpdSemVAOAMSaMP5LyCH4hzBh+aPln4EPW2r/b5p3XAb+IX5phc7HND92m2OZ78j9jDXgG+M29FttUwBERORyO6/H1oZuMXVsil1v360vl/PpS/d3NtN8Vo6u14fYfkiN3ogJOuVLAERE5eIVpqUtjCy8qwRCLVhJVfamyc5JqUYmIiOxo8lqS4eklopEQrvdCfalgIKD6UidE2R70JyIisp1CCYZU2qEqEqSlPkY6k6MyHORsW43qS50QCjgiInIsOK7H7NwyiWV3owSDk/OIRSAeC3HmVA29Z3SQ30mhKSoRETkWxmaWmb2xytOXrxKN+CEmHAyxtr7GqYZqhZsTRgFHRETK3mrGZeZmkktjc9TVRFWCQRRwRESkvKXSLkOTCVKZHHU10Y0SDInlNAM9KsFwUmmb+CHTNnERkf1TOMjvwqBfGTzr5qgIBFleydB/rpm+Mw0KN8eAtomLiMiJkUq7DE764aamOsLzw3OsrwdYzbgM9CrcnHQKOCIiUnZSaZehqQSXRuepqfYrg3e317GwlKavs57zHQo3J50CjoiIlJVCZfDBiQS18cjGmpulZJpX9J2ir7NR4UYUcEREpHxsrgweDgdZW1+jp80PN/09zTTWVyrcCKCAIyIiZcJxPezEAsNTS8SjYYamEkTCIaoqg3SfbqCpNkZ7U+1Rd1NKhHZRHTLtohIRKZ7jegxOLjA8tUjG9fzK4K5fGfzecy2cbq2mo6XuqLspB0S7qERE5NhxXI+hyQQTV5NkHA/X9VhecQhWBGg/VUMkEuJUffyouyklRrWoRESkZBWmpa7OreC4Hm7uhcrggQB03qXimbI9jeCIiEjJmryWZGRmiZW0S7QqBAE/3IRCFXS11SncyI40giMiIiVpNeMyMp0g6+Rwch7BYICWuhir2RxdrTWqLyW3pIAjIiIlZ3kly+jMEqtZDyc/LZXOrFERXWOgp5Hu9nqFG7klTVGJiEhJSaVdhqcW+drlq0QjQdbXIRwMEQxCx6kahRvZFQUcEREpGam0y/jsIhdH56irifLc8Bw1sQg5b42zbbWalpJdU8AREZGSUKgM/oy9QVNDjLF8CYb5xRT93Y0MdDfrlGLZNR30d8h00J+IyEs5rselsXkuXLlOvDpCctWhramaucQqA+ea6e2opzYeOepuyhHZy0F/CjiHTAFHROTFCtNSz4/MEQoGGctXBs9kXe4zp+hqr9fIzQmnk4xFRKSspNIuF0fnGJxaJB6tZDQfbpaSaQZ6mhVuZM+0TVxERI6E43oMTSUYmlwkGAxQEYGefLjp71ZlcLkzGsEREZFD57geozOLXB5fIBYN43oey6kX6kvFomFVBpc7ooAjIiKHqlBfanhqiXg0zNBUgpb6GPFoJbGqEH1nGunvbtJ2cLkjWmR8yLTIWEROMsf1GJxcYHhqkYzrEY0Ecdx1lpJp7j3XwunWajpa6o66m1JitMhYRERKluN6DE0mmLiaJON4uK7H8soL01KRSIhT9fGj7qYcE1pkLCIih2LyWpKZuSSO69eXCuBXBg8EoPOuGlUGl32lERwRETlwhcrgK6su0Sr//1uHQ0FCoQq62uoUbmTfaQRHREQOlF88M8Fq1iOVcQkGA7TUx0hnXbpaVV9KDoZGcERE5MAsr2Sx4ws8feU60YgfYtKZNTJODtPZyD29LTrrRg6EAo6IiBwIf+RmkWfsNWqqI35l8KhfGfz0qThd7XUauZEDo4AjIiL7rlBf6uLoHHU1UUbzlcEXltMM9KgyuBy8fVmDY4w5D7wcmLDWfm0/vikiIuUplXYZnFjgG6NzNDXE+MbI3EYJhvvubqWvs1HhRg7crkdwjDHfZYz5G2PMN2+5/wHgMvDnwFeNMR/f5z6KiEiZKNSXenbwBlWRMNfmU9x7rplMxuW+u1vp7VDxTDkcxUxRvQt4EPh64YYx5h7gPwBrwJeAReD7jDHftZ+dFBGR0pdKuwxNLnB5fJ54dYSx2SVqYpXMLa5y392nMGcbqY1HjrqbckIUE3BeCTxnrV3ddO9dwDrwv1lrHwReDbjAj+1fF0VEpNQ5rseVsXmGpheJRysZnV2iOz8tNdDTTFe7Rm7kcBUTcJqAmS33Xg+sAP8dwFo7Cvwj8LJ96Z2IiJS8QmXwietJMlmPtbU1etrqWEym6e9uprG+UuFGDl0xAScCbBS5MsZUAq8AvmKtzW1qdw24a3+6JyIipWxzZfBoJITreSynHCoq4HRLDbFomPam2qPuppxAxeyiugr0b/r7g/ih50tb2sWB5Tvsl4iIlIHJa0lGZpbIOjmqIkH/hOJMjspwkM5W1ZeSo1NMwPki8C5jzPuBzwC/gr/+5jNb2t0DTO9P90REpFQV6ktlnRxOzi+gWV0F8ViIjpYaes+oBIMcnWKmqH4Nf73NfwIuAN8M/IO19n8WGhhj+oAe4Kv72UkRESkdjusxfjWBnfDrSzk5D4BwMEQgAO3NcYUbOXK7DjjW2kHgW4D/CnwaeBR4+5Zm/wJ4DvjUPvVPRERKzNjMMrM3Vnn68tWN+lLhYIhgEFUGl5IRWF9fP+o+nCiBQGAdQP/eRaQcrWZc/unrVxmdSWyUYOg700Aq7dJ3tl4lGORABAL+Hqf19fXAbZpuKOYk4wVjzJN76JeIiBwD/kF+CVKZ3IvqSyVUX0pKUDFrcCqBqYPqiIiIlK5CCYZnrlwnHApQGQ7Q39WE4+R45d13qb6UlJxiAs4w0HxQHRERkdK0tQTD88NzrK0FWM24DPQ203emQeFGSk4xAefjwIPGmO6D6oyIiJSWVNrl4ugcg1MvLsEwv5ji/Jl6znco3EhpKibg/Dbwt8DnjDHvNMaoYpqIyDFWmJYamlwk6+RLMOTrS6kEg5S6Xe+iMsaM4pdqOIt/wB/ADSC9TfN1a+25fenhMaNdVCJSDgr1pS4M3qSqMsRKOgvrUB2tpCoSpL05Tn93k7aDy6HYyy6qYk4y7tr8s/LXnWpO6b/eIiJlqlBfauJqkng0zOXxBV7e20w6kyNaFaK3o4HO1hqFGylpxQQcrb0RETnmHNdjcHKBkelFMq5HNBKku72O6evL3HuuhdOt1XS01B11N0VuSwf9HTJNUYlIqXJcj6HJBONXl0muOrj5EgyalpKjdtBTVCIicoxNXksyM5fEcf36UgEgHAoSCEDnXaoMLuVlTwHHGPNa4A3A6fytGeAL1tqv7FO/RETkEBUqg6cyLrGqMMm0QzgUJBSqUH0pKUtFTVEZY7qAPwVeU3g/fy185CvAu6y14/vUv2NHU1QiUmpSaZfhqQSjs8uspLPUVldSGQqRzrp0tdZiupq0HVyO1F6mqIrZJt4APIO/TTwFfBIYwQ853cBbgTgwDrzKWpvYfddPDgUcESklyytZhqcWeXboBu0tcW4urhIOhohUVtB3poHznQ0auZEjd9BrcN6PH24+AbzbWju/+aExphH4CPAw8HPALxTxbREROWT+yM0iz9hr1NVEeW54DnOmgZW0S8+pWrra6xRupGwVc5Lx24Gr+FNQ81sfWmsXgB/It3nH/nRPREQOQirtMj67yMXRuRdVBl9QZXA5JooJOF3AU9ba7E4N8s+e4sWHAoqISAlJpV0GJxZ4xt6gqSHGWD7cLCXTvKLvlCqDy7FQTMBxgdgu2kXzbUVEpMQU6ks9O3iDqkiYa/Mp7j3XTCbjct/drfR21CvcyLFQTMC5DLzRGNO6U4P8s4fybUVEpISk0i5DkwtcHp8nXh1hbHaJmlglc4ur3Hf3KczZRmrjqqMsx0MxAefjQDXw98aYh7Y+NMa8EWKCfXYAACAASURBVPg7/FGeP9mf7omIyH5wXI8rY/MMTS8Sj1YyOrtEd35aaqCnma52jdzI8VLMNvEQ8Fng9fjn3swCY/l/7sY/9C8AfB74DmutdxAdLnfaJi4ih61QGfzS2AJOvr6U466zuJJmoLuZ9rtidLU2HHU3RXZ0oOfgABhjIsCvAO/GP/NmsxX8beL//lYLkU86BRwROUybK4MHgwHmllZhHWJVYaKREO0tqi8lpe/AA06BMaYKeBUvjNpMA09bazNFf+yEUcARkcM0PLXIN0bnyDo5qiJBIuEQ6UyOynCQs22qLyXl4aBPMn4b4FprP72n3gmggCMih8NxPW4kVvj6yALLK1mcQmXwqkrisRAdLTX0ntEpxVIeDvok478E/h5QwBERKWGO63FxdJ6skyOTyeG4frgJh0IEAtDeHFe4kWOvmF1UC8DcQXVERET2x+S1JFfnV/jGyE2iVSEIQDgUJBhElcHlxCgm4PwzcM9BdURERO7casZlZDpBJpvbqC8Vj0bwvHXOttXSp+KZckIUE3B+AxgwxvzoQXVGRET2zj/IL8Fq1iOVcakMB+jvaiLjuNzT26T6UnKiFLMGB/xt4H9ojHkYf03OBJDerqG19sk77JuIiOzS8kqW4alFnh26QXtLnJU0pDNrRCrh5edaOK+RGzlhigk4X8A/1C8AvAn4jlu0XS/y2yIiskeptMvw1CLP2Gsb01LmTAMraZeeU7V0tdcp3MiJU0wIeRI/uIiISIlIpV3GZxe5ODpHXU2U0Xxl8IXlNK80qgwuJ9eeDvorJcaYbwV+BngAaMTf7fV14HestX+zpe0DwAeA1wBVwDDwR8Dv7VRawhjzQ8BPAv2AB1wAHrPWfmov/dU5OCKyXxzX49LYPN8YvkljfYyvj8zRk68vdd/drZizCjdyPOzlHJxdLzI2xrzcGFNSu6iMMR/AH1l6EPgM8EHgk0AD8IYtbd++qe1fAr8PVAK/DTy+w/cfAz4GtAEfxS84ei/wSWPMe/b79xER2a1CfanL4wtEImGuzae491wzmYzLfXe30tuh4plyshVzkvEa8JS19vUH26XdMcZ8D/Dn+IcPfpe1Nrnledha6+b/uRZ/tKYOeJ219mv5+1XA54DXAt9nrX180/sPAF8CRoBXW2sT+ftdwNP4ldXvttaOF9NvjeCIyJ1KpV0GJxaYvrkCwOXxBbrb61hNZ3nVy1rpVmVwOWYOdAQHWASmiuzTgTDGVOBvW18Fvn9ruAEohJu8h4EW4PFCuMm3yeBPWQH8xJZPvDt//bVCuMm/M44/+hMBfvjOfhMRkeI4rseVsXlGZhdJpV3W1tbozk9LnT/TSLw6pHAjQnEB51ng3EF1pEgPAN3A3wAJY8y/Msb8vDHmp40xr92m/UP562e2efYkflB6IF8tfTfvfHpLGxGRA1eYlpq4niSXW8fNeSynHCoqArSfqiEWDdPeVHvU3RQpCcUEnA8B32yMefNBdaYIr85frwPPAJ8Cfh34HeDLxpgvGmNaNrU3+evg1g9Za3PAGP6Osh4AY0w1fqX0FWvt1W1+/lD+2neHv4eIyK44roedWGB4aoloJITrFepLBQkGAnTeVUN/d5O2g4vkFbNN/Bngw8BfG2P+iNsf9Dd5593b0an89d344eTbgK8CZ/EXGr8J+AteWGhcl78u7fC9wv36PbZ/EWPMo8Av79R5EZFiTV5LMjKzRNbJURUJ0lIfI53JURkOcratRvWlRLYoJuCM5a8B4Mfzf3Zy0Af9Ff5XHAAettY+l//7RWPMd+KP1LzeGPNaa+1XdvG9wqKlYlf+btveWvso8Oi2Pyi/yFhEZDcc1+NGYoWh6QRZJ4eT83ByHtVVEI+F6GipUWVwkW0UE0KmKJ2D/gqLfkc3hRsArLVpY8zfAj8KfBPwFV4Ycalje4VJ66Ut153a326ER0RkX4zNLLOYTJPJ5Mi6HgEgHAoRCEB7c1zhRmQHuw441tquA+xHsWz+urjD80IAim5qfz/+mpmnNzc0xoTwFyzngFEAa23KGDMDnDbGtG2zDud8/vqSNT0iIvtlNeMyczPJ6EyCjrtqSaYdf81NELra6jQtJXILxSwyLiVP4geS88aYym2eFw4kHM9fP5e/brdA+kEgBnzZWpvddP9W77xlSxsRkX1VqAyeyuQ26kvVxCLkvHXOttXSp+KZIrdUlgHHWjsH/Bn+VNEvbX5mjPl2/EXGS7ywxfsJYA54xBhz/6a2VcCv5v/6B1t+zEfy1180xjRseqcLv3RDFvjjffh1RERepHCQ39NXrhMOBagMB+jvamI17TBwromB7maddSNyG0XXojLGvB54D/7pvy3Ax621P5p/9ibg9cCHrLXX9rmvW/txCv+k4V7gKeCf8XdRfSf+WqHvt9b+xab278APOhn80gwLwNvwt5A/AXyvtXZ9y8/4IPA+YDrfphJ4J9AEvNda++Fi+62TjEXkVlJpl8HJBS7YG9RURxibXaKvs4Gs42G6Gjjf0aBwIyfOQZ9kXNj+/Dngu4F2IMwLO5AAXODn888PlLX2BvDN+LWkzgA/hX/w3v8AvnVzuMm3/yv88PVkvn/vzff3fcAjW8NN/p2fBf4NcA1/19gPAheBt+4l3IiI3Eoq7TI0leDS6Dw11RFGZ5fobq9jYSlNX2e9wo1IEYqpRfVW4K/xd1O9Dz8oXAc+Zq39kXybQP7e09bat+z0rZNMIzgisp3CtNTg1CKxqhCXxhdUGVwk76BHcH4Kf93Jm621n7DW3tzaID8KMog/bSQiIrvguB5DUwmeHbxBZTjI2vraRrjp72mmsb5S4UakSMUEnFcB/2StvXybdtNA6967JCJychTqS10eXyBeHWFoKkEkHCIaCdHT0UBTbUz1pUT2oJiD/qLAS0ZttlGzx76IiJw4k9eSDE8vEY+GuTy+QHd7HdPXl7m/v42+zgZiVRq5EdmLYgLOVeDuXbTrx69RJSIiO9hcgiGVdqiKBOkuTEt1N1NfG1a4EbkDxUxRfR4YMMZ8x04NjDHvxN+q/dk77ZiIyHHluB4XR+eZvp4ik8nhuB7LKw4VgQBtLXFi0bCmpUTuUDEjOL8J/GvgL4wxPwd8ovDAGBMDHgY+BKzmryIiso3Ja0muzq8weyNJx121rGQcwkF/cXF7cx393U06pVjkDhV10J8x5hHgY/jn36zjn4Hj8UJ17xzwA9baP9/fbh4f2iYucrKtZlyeujBNKuPircHo7BJ9ZxpYTbucP1uvU4pFtnHgB/1Zax8HXo1/qu8KfsAJ4Z8O/CngAYUbEZHtFepLrWY9Uhl3owRDxnG5p1clGET2U9GlGgryh/o14Y/ezFlrvf3s2HGlERyRk2l5Jcvw1CLPDt2gvSXOzcVVwsEQkcoK+s40cF7FM0V2tJcRnD0HHNkbBRyRkyeVdrHjCzxjr1FXE2V0dglzpoGVtEufpqVEbuvAp6hERKQ4qbTL+OwiF0fnNsJNT3sdC8tpBnoaFW5EDogCjojIASmUYLhgb9DUENsIN0vJNK/oO0Vfp+pLiRwUTVEdMk1RiZwMhRIMFwZvEgpWkFx1aG2qZj6xysC5Zno76qmNR466myJlQWtwyoACjsjxV6gMPn1zBWCjBMNqOsurXtZKd3u9Rm5EiqA1OCIiR8xxPa6MzTMyu0gq7bK2trZRguH8mUbi1SGFG5FDsGPAMcbUGmOqDrMzIiLlrDAtNXE9SS63jpvzWE45VFQEaD9VoxIMIofoViM4CeD3C38xxvySMeZtB98lEZHy47gedmKB4aklopEQrucfDRYOBQkGAnTeVaMSDCKH6FYBJ5D/U/Ao8I4D7Y2ISJmavJZkZGaJVMbByeVoqY8Rj1ZSE6uk90wdfZ2NCjcih+hWxTZXgcbD6oiISDlyXI8biRWGphNknRxOzsPJeVRXQTwWoqOlht4zOqVY5LDdKuBcAb7dGPMjwHD+Xqsx5sHdfNha++Sddk5EpNSNzSyzmEyTyeTIuh4BIBwKEQhAe3Nc4UbkiNwq4HwE+EPgo5vuvSn/53bWb/NtEZGyt5pxmbmZZHQmQcddtSTTjr/mJghdbZqWEjlKO4YQa+1/McZcAx4GzgBvBG7gj+yIiJxoqbTL8FSCVCZHXU2U54bn6DvTQCrt0tNRR5+KZ4ocqV0f9GeMWQM+Zq39kYPt0vGmg/5Eyl/hIL8Lg35l8KyboyIQZHklQ/+5ZvrONOisG5F9dNAH/f1X4B+L7JOIyLGSSrsMTvrhpqY6wvPDc6yvB1jNuAz0KtyIlAqVajhkGsERKV+ptMvQVIKLI3NUVoY2imcmU1leaVQ8U+Sg7GUEZ08LgY0xrwXeAJzO35oBvmCt/cpeviciUuoK01KDU4vUxiNcGl/YqAx+392tCjciJaaoERxjTBfwp8BrCu/nr4WPfAV4l7V2fJ/6d+xoBEek/Diux6WxeS5cuU59bZSKinUcd52lZJr+nmbaT8Xoam046m6KHFsHWk3cGNMAPAOcBVLAJ4ER/JDTDbwViAPjwKustYndd/3kUMARKS+F+lIXBm8SClYwNrvEy3ubyXnrxKpCnG6uoft0rXZMiRygg56iej9+uPkE8G5r7fzmh8aYRvyzcx4Gfg74hSK+LSJScgrTUtM3V4hHw1weX6C7vY7p68vc399GX2cDsSpNS4mUomJ2Ub0duIo/BTW/9aG1dgH4gXwb1awSkbLmuB5XxuYZmV0klXZZW1ujO7/mpr+7mfrasMKNSAkrJuB0AU9Za7M7Ncg/eyrfVkSkLBWmpSauJ8nl1nFzHssph4qKAO2naohFw7Q31R51N0XkFooJOC4Q20W7aL6tiEhZmryWZHh6iWgkhOt5AH4JhkCAzrtq6O9u0pobkRJXTMC5DLzRGNO6U4P8s4fybUVEys5qxmVkOkEq7eDkcrTUx4hHK6mJVdJ7RvWlRMpFMQHn40A18PfGmIe2PjTGvBH4O/xRnj/Zn+6JiBwOx/UYv5rATiRYzXo4+WmpVDpHLBqiu71W4UakjBSzi+ojwHcDrwc+a4yZBcbwz8Dpxj/0LwB8Pt9WRKRsjM0sM7+0yqWxOTruqmUlDeFgiLX1NU41VNN7RsUzRcrJrkdwrLU54M3AY/jn4JwGvgX4VqAjf+8x4F9aa73976qIyMFYzbjM3ExyaWxuozJ4PBoh561xtq1WlcFFytCealEZY6qAV/HCqM008LS1NrO/3Tt+dNCfSGlJpV2GpxJMXl8h4+RUX0qkBB3oScayPxRwREpH4SC/C4M3aG+Jk3VzVASCLK9k6D+nyuAipWIvAaeYRcYiIsdGKu0yOOmHm5rqCM8Pz7G+HmA14zLQq3AjUu4UcETkxEmlXYamElwanaemOsLo7BLd7XUsLKXp66znfIfCjUi5U8ARkRPFcT2GphIMTiSojUc21twsJdO8ok9rbkSOCwUcETkxCiUYLo8vEA4HWVtfo6ctX1+qp5nG+kqFG5FjQgFHRE4Ex/WwEwsMTy0Rj4YZmkoQCYeoqgzSfbqBptqY6kuJHCPaRXXItItK5PA5rsfg5ALDU4tkXI9oJIjjrrOUTHPvuRZOt1bT0VJ31N0UkR0c6C4qY8zbjDFv2UO/RESOjON6DE0mmLiaJON4uK7H8opDMF8ZPBIJcao+ftTdFJF9VswU1V8CP3NQHRER2W+FaamZm0kc18PNvVAZPBBAlcFFjrFiAs4CMHdQHRER2W+T15KMzCyxknaJVoUg4IebUKiCrjZVBhc5zooptvnPwD0H1RERkf20mnEZmU6QdXI4OY9gMEBLXYzVbI6u1hrVlxI55ooZwfkNYMAY86MH1RkRkf2QSrsMTSZYzXo4+WmpdGYN11tjoKeRe3pbtB1c5JgrZgQH4CPAHxpjHsZfkzMBpLdraK198g77JiJStI0SDNavL7WShnAwRDAIHadq6G6v18iNyAlQTMD5ArCOXz38TcB33KLtepHfFhG5Y1tLMDw3PEffmQZSaZeejnpNS4mcIMWEkCfxg4uISMkpVAYfnFqkNh7h0vgCPe11zC+muO/uVsxZlWAQOUl00N8h00F/IvvPcT0ujc1z4cp16mujVFSsbxzk19/TTPupGF2tDUfdTRHZowM96E9EpBRtri8Vr45slGCIRkL0dKgEg8hJtecRHGNMJdAEZK21C/vaq2NMIzgi+6cwLTV9cwWAy+MLdOcrg9/f30ZfZwOxKk1LiZS7QxnBMcb8oDHmfwIpYBp4bNOzh40x/90Y013sd0VEiuG4HlfG5hmZXSSVdllbW9sIN/3dzdTXhhVuRE6wogKOMeZjwB8Dr8LfHr41SU0BjwAP70fnRES2U5iWmrieJJdbx815LKccKvL1pWLRsKalRE64Yopt/hDwg8BzwP3AS0rvWmu/CswCKsopIgdm8lqS4eklopEQrvdCfalgIKD6UiICFLdN/MeAJPBWa+0MgDFmu3YjQNcd90xEZBuFEgyptENVJEhLfYx0JkdlOMjZthrVlxIRoLiAcy/wT4VwcwuzwKv33iURkZdyXI/ZuWUSy+5GCQYn5xGLQDwW4sypGnrP6CA/EfEVswYnDKzsol0T4O6tOyIi2xubWWb2xipPX75KNOKHmHAwxNr6GqcaqhVuRORFigk4k9ymmrgxJggM4E9TiYjsi9WMy8zNJJfG5qirifLc8BzxaISct8bZtlqVYBCRlygm4Pwt0GuMedct2vxboA34H3fUKxGRvEJl8FQmR11NlNHZJXra60gspxnoaWSgu1klGETkJXZ90J8xpgP4BhDDP/vmCeBrwOPArwDfA/wC/kLkAWvtjYPocLnTQX8iu1c4yO/CoF8ZPOvmqAgEWV7J0H+umb4zDQo3IifAXg76K+okY2PMG4FPsM0WcfwzcZaBt1trv7jrj54wCjgiu5NKuwxOLnDB3qCmOsLY7BJ9nQ1kHQ/T1cD5DoUbkZPiwAMOgDGmFfg/8M+66QGC+Af8fRr4TWvtdFEfPGEUcERuL5V2GZpKcHFkjsrK0Ma0VDKV5ZXmFH2dqgwucpIcSsCRO6OAI3JrhcrggxMJolUhLo0v0JMvwXDf3a2Yswo3IieNqomLSFnbXBk8HA6ytr5GT1u+vlRPM431lQo3IrIrxRz0t8EYcxp4EOgA1vEP93tS01MisleO62EnFpi4miQeDXN5fIGX9zZTwRqNpxtoqo2pvpSI7Fqxi4xbgN8DvpuXjv6sA/8v8B7toNqZpqhEXspxPQYnFxieWiTjekQjQRx3naVkmnvPtXC6tZqOlu32NojISXCga3CMMY3Al4HzwBrwT8B4/nEX8Br8BcdDwGuttQu77cRJooAj8mKO6zE0mWD86jLJVQc35xfPrI5WUhUJ0t4cV/FMkRNuLwGnmCmqR4E+4B+Ad1trX3RasTGmG/gI8G3ALwM/XcS3ReQEKkxLXZ1bwXG9jXATDgUJBKDzLhXPFJG9KWaR8TuAm8A7toYbAGvtGPBdwBzwnfvTPRE5ziavJRmZWWIl7RKtCkHADzehUAVdbXUKNyKyZ8UEnFPAF621qZ0a5J99EWi5046JyPG2mnEZmU6QdXKkMi5ZN0dLXYzKcJAe1ZcSkTtUzBTVDFC5i3aV+LuqRERewnE9ZueWSSy7rGY9nPy0VDqzRkV0jYGeRrrb6xVuROSOFDOC8xfAQ/mTjLeVf/YQfjkHEZGXGJtZZvbGKk9fvko04oeYcDBEMAgdp2oUbkRkXxSziyqGv8C4HniftfbTW56/Gfggfj2qh6y16X3u6276+APAf8v/9cestf9lmzYPAB/A3/VVBQwDfwT8nrXW2+G7PwT8JNAPeMAF4DFr7aeK7aN2UclJtppx+aevX2V0JrFRGbzvTAOptEvf2XpVBheRbe3rLipjzOe2ue0BBviUMWaRF28Tr8//81eATwH/Yred2A/GmDP4Z/SsAPEd2rwdf3QpA/wZsAC8Ffht4HX4FdG3vvMY8LPANPBR/Cm4R4BPGmPea6398L7/MiLHUCrtMjyVIJXJbYSbnvY6Estp1ZcSkX234wiOMWbtDr67bq09tDFmY0wA+CzQjX/Y4L9jywiOMaYWf7SmDnidtfZr+ftVwOeA1wLfZ619fNM7DwBfAkaAV1trE/n7XcDTQDVwt7V2fLd91QiOnESF+lIXrlynrSVO1s1REQiyvJKh/1wzfWdUGVxEdrbf5+C88U47dIh+Cn/tzxvy1+08jL+7678Vwg2AtTZjjPkA/vTbTwCPb3rn3fnrrxXCTf6dcWPM7wP/Hvhh/HN/RGQbqbTL+Owil8fniVdHeH54jvNnGsi4LgO9zZzvULgRkf23Y8Cx1n7xMDuyV8aYlwG/DvyutfZJY8xOAadw/zPbPHsSWAUeMMZErLXZXbzzafyA8xAKOCLbSqVdLo7OMXNzhXi0cqMy+PxiivvublW4EZEDU9bVxI0xIeBPgEngF27XPH8d3PrAWpsDxvADX0/+29XAaWDFWnt1m+8N5a99xfdc5PhzXI+hqQRDk4tkHY+1tTV62vOVwbtVGVxEDtaeqomXkF8CXgl8yy52bRUq9S3t8Lxwv7BYutj2G4wxj6JRHTnBHNdjdGaRy+MLxKJhVtJZ3JRHdbSS9lM1xKJhVQYXkQNVVMAxxjQA/zv++px2/G3W21m31p67w77dri/fhD9q80Fr7Vf24ZOFhUvFrv59SXtr7aP4tbte+kPyi4xFjqtCfamJq0ni0TCXxxd4eW8z6UyOaFWI3o4GOltrdNaNiByoXQccY0wvfhmGVl4IAzs50P+Ib5qaGsRfB7MbhRGXuh2e125pd7v2txvhETlxHNdjcHKBkelFMq5HNBKku72O6evL3HuuhdOt1XS07PQ/KRGR/VPMCM4HgTbgKfxzY4bwz5w5CnFeWPuSMcZs1+ajxpiP4i8+/hnAAvfn33t6c8N8YOoGcsAo+HW1jDEzwGljTNs263DO568vWdMjchI5rsfQZIKJq0kyjl8Z3HVfmJaKREKcqt/2iCoRkX1XTMB5A/7Bft9urXUOpDe7lwX+7x2e3Ye/Lucf8UNNYfrqc8C/Bt4M/D9b3nkQiAFPbtpBVXjnB/Lv/PGWd96yqY3IiTd5LcnMXBLH9etLBfArgwcC0HlXjSqDi8ihKqZUwyLwGWvtIwfbpTuzaYHvdgf9jeBPRemgP5F9tJpxeerCNKmMS6wqzI3EKpWhIKFQBedO1zPQ06RwIyJ7tt8H/W31LP76m7JkrV02xvwY8ATwBWPM4/ilGt6Gv4X8CfzyDZvf+bIx5reA9wHPG2OewC/V8E6gEXhvMeFG5DgqlGBYzXqkMi7BYICW+hjprEtXay19nQ0KNyJy6Io5B+cx4FvyoxplyVr7V8Dr8Q/2+27gvYCLH2AesdZutyPqZ4F/A1wDfhz4QeAi8FbVoZKTbnklix1f4Okr1zcqg6cza2ScHKazkXt6W3TWjYgciV1PUQEYY94L/Afgw8Df4heg3LZmlbV2cj86eNxoikqOi1TaxY4v8Iy9tlE805xpYEWVwUVknx30FBXABeA6/vkztzo5eH0P3xaRMlGoL3VxdO5FlcEXVBlcREpEMefgvAG/JlNl/tY8R7dNXESOSCrtMjixwDdG52hqiPGNkbmNEgz33d2qcCMiJaGYUZZfwQ83/xn4dWvt4sF0SURKVaG+1LODN4hXR7g2n+Lec83MJVa57+5WejvqFW5EpCQUs008CVyx1r76YLt0vGkNjpSrwrTU8yNzhIJBxmaX6G6vI5N1uc+coqtd4UZEDsZe1uAUs4sqzQsVtEXkBHFcjytj8wxNLxKPVjKaDzdLyTQDPc0KNyJScooJOE8BAwfVEREpTYXK4BPXk2SyHmtra/S01bGYTNPf3UxjfaXCjYiUnGICzr8HzhljfvqgOiMipaVQGXx4aoloJITreSyn/v/27jw8zqu8+/h3NJJGy2ixZTmWvEp2fIKdUBrWJi0JtKWlUJaWUFqWpgttIJRS0pa2LE1oy1qWEChQeEug5SqkUCh7QklCaKCkJA4QO74tWZYt79a+jTSL9P5xnlHGE43W0Ugz+n2uy9fjedYzxyPPrbPdccrKYGtzHTXVFbQ21c9/IxGRAlvMGJxXAU8FXgt8n/nXwflMnspYUjQGR4pJZ88gj3T1MhlPUhUJE6koJzaRpLIizM4W5ZcSkcJY6XVwbsevbxMCrgbmW9FYAY5IERufSHD05ACT8STxpE+gWVsF0ZpytjXXsWe7UjCIyNq1mADnM/gAR0RKWDyR4nTvMAPDCcYnfWADUBEuJxSC1k1RBTcisuYtKlWDLJ+6qGSts+4B+obGOXSsl22X1HNhcJyKcDnhMOzZtkGZwUWk4FZ6mriIlLjxiQSnLoxw6JhPwfDjzl6i1RGSqSl2tigzuIgUDwU4IgL4hfw6TgwwNpG8KL/UwHCM/e0blTxTRIrKYmdRLZhmUc1OXVSyFsUTKQ4d6+PA4XO0NEeZTCQpC4UZHp1g3+5N7N2+QcGNiKyaQs2imrccwXkKcESKQDoFw6PdfURrI/yks5dLt29gIpFg/55NXLpNwY2IFJ98zKIqA3YCVwK1wJeBoeUXTURW2lgswcGuXk5dGCVaXcmh7n7aWxvoGxzjysu2KLgRkaKVt1lUzrnN+CCoFbjKzEbzcuMSoy4qWSvS3VIHj/YRDoeojoSJJ6cZClIwtF5Sw64tG1a7mCIiqzuLyszOA78DbAVuztd9RST/0vmlHu3up6a6YiYFQ7gsROtmpWAQkeKX11lUZtYP/B/wm/m8r4jkT2Z+qWh1BR09AzQ31hCtrqSmqpy92zeyr01r3YhIcVvMGJyFigMtK3BfEVmmeCLFkRP9HD05yEQiRXUkTFtrAyfPDXPF7ma2bqllW3PDahdTRGTZ8tqC45zbgs9TdSGf9xWR5YsnUnScGOD4mREm4ikSiRTDo491S0Ui5WxuHTtE6gAAIABJREFUjK52MUVE8mLBLTjOuWfOcTgKXAbcCDQC/77McolInp04O8Kp3hHiCZ9fKgRUlIcJhWDHJcoMLiKlZTFdVPcy/zo4IeAA8JalFkhE8i+dGXxsIkFNVQUjsTgV5WHKy8vY1dKg4EZESs5iApz7yB3gxIFTwHeAO8wssdyCiUh+jMUSdPYMMD6ZYmwiQTgcormxhthkgl1blF9KREqTsokXmNbBkUIaHp2ks2eQhzvO09ocnckMHqksY+/2DVyq4EZEioCyiYvIDN9yM8hDdpa62gg/7uylLsgMvnVzlF2tDQpuRKRkKcARKUHp/FIHu3ovygzer8zgIrJO5ByDM8+sqXmZ2X3LuV5ElmYsluDI8X4e6eqlaUMNjxztpb21gaGRGFdetoW9OzYquBGRkjfXION7WVj28NlMz3NvEVkB8USKjp4BHj5ynmhthLN9Y1yxexO9A+NcedkW9mxrVHAjIuvCXEHIIRYf4LQBNUsvjogsVbpb6tHuPqK1EY6dHqKttYHewXGuvGwzu1oV3IjI+pEzwDGzyxd6E+fcfuAdwL5g18lllktEFmEsluBgVy+nLowSra7kUHf/Rd1SCm5EZL1ZVjeSc2478HbgFfgBywPAu4Dbll80EVmIdLdUx4lBwuEQZRFmgpt9bZvY2Fip4EZE1p0lBTjOuSbgzcANQBUwDtwKvNvMhvNXPBGZSzyRouvUII9291NTXcFobJLEWIra6kpaN9dRU11Ba1P9ahdTRKTgFhXgOOdqgJuCP3VACvgY8HYzO5v/4onIXE6cHaHz5BDR6goe7e7niXs2EZtIUl1Vzp5tG9ixpU5r3YjIurSglYydc+X41po3A5uD3XcAbzGzoytXvNKjlYwlH+KJFOcHRvnp0X6GRyepioSJJ6YZGolxxe5mtm6pZVtzw2oXU0QkL5aykvG8LTjOuZcDt+BnSIWAu4C/NrMDSyumiCxHPJHiYFcfk/EkExNJnx08kaKmqpKW5iiRSDmbG6OrXUwRkVU110J/zwXeCVyBD2weAP7KzO4tTNFEZDYnzo5wpm+U0+dH2HZJPaMTcSrCYaamp2jd1MC+tiZ1S4nIupezi8o5N4VfB2cc+BDwxcXc2MweWnbpSpC6qGQ5xicSfO/AScYmEqSmoOv0EHu3b2A8luDSnY1KwSAiJWkpXVQLCXCWYtrMtJLxLBTgyFL55JkDdJ0eZjQ2SX1tJWWhMKOxSfa1NbFn2wYFNyJSkvI9BucESw9wRCSPhkcn6ewZ5OGO87Q2RxmNQWxiikglPHF3M5fu2KBuKRGRDAuaRSX5oxYcWayxWALr7uchOzuTGdxt38BoLMFedUuJyDqwlBacshUrjYgsWzq/1MGu3pngpr21gf7hGPvbNyq4ERHJQQGOyBqVTsFwwM7TtKFmJrgZGonxpL2b2btjo4IbEZEc1EVVYOqikoVIp2A4cOQC5eEyRsbjbGmqpW9gnP27N7FnWyP10chqF1NEpCDyOotKVoYCHJnPWCzBkeP9nLwwCsCj3f20tTYwHpvkyU/YQpsyg4vIOqMxOCJFLp5IcfhYH0dPDzIWSzA1NUVb0C116faNRGvLFdyIiCyAAhyRNSLdLXX83AjJ5DSJZIrhsThlZSFlBhcRWSQtxieyBsQTKex4P8fPjFAdKWdsIg5ARXmYcCjEjkvq2Ltjo9a6ERFZIAU4ImvAibMjHD01xGQ8SVUkTHNjDbGJJJUVYXa2KLgREVksBTgiqyieSHF+YJSOkwNMxpPEkyniyRS1VRCtKWdbcx17tmuVYhGRxVKAI7KKjp0aZnAkxsREkslEihBQUV5OKAStm6IKbkRElkiDjEVWyfhEglMXRvjp0QtUV/nfNSrKw4TDsKulQd1SIiLLoABHZBWMxRJ0nBhgbCJJQ101P+7spa4mQjI1zc6WevYqeaaIyLJoob8C00J/kl7I78ARnxl8MpGkLBRmeHSCfbs3sXf7Bq11IyKSQQv9iaxxY7EER0744KauNsJPOnuZng4xPpFg/x4FNyIi+aIAR6RAxmIJOnoGONTVR11thK7TQ7S1NtA/FGPvjkYu3abgRkQkXxTgiBRAulvKjg9QH40oM7iIyApTgCOywuKJFB09Azx85DyVFWGmpqdmgpt97ZvY2Fip4EZEJM8U4IisoHR+qUe7+4nWRujoGSBSUU51pJz2bRtoqq9RfikRkRWgWVQFpllU60tnzyCHj/czPT3No939M5nBn7Kvhb07NlBTpZYbEZH5LGUWlVYyFlkBmSkYxmJxqiLhmeBmX9smGusrFNyIiKwgBTgieRZPpDjY1cdkPMnERJJ4IkU8kaKmqpKW5ig11RXqlhIRWWEagyOSZyfOjnCmb5RH0ikYQj4Fw9T0FK2b6tjX1qRVikVEVpgCHJE8Gp9IcPTkABOTj6VgiFZHSCkFg4hIQWmQcYFpkHHpGosl6OwZoOv0MKOxSeprKykLhRmNTbKvrYk9WshPRGRJlKpBZJUMj05i3f08ePgc1RHfQhObmCKZSvHE3c08YVeTghsRkQJSgCOyTL7lZpCH7Cx1tRGfGbw6QjI1xdbNUXa1NqhbSkSkwBTgiCzDWCxB9+lBDnb10lBXPZOCoX84xv72jexv26SWGxGRVaAAR2SJ0ikYDth5mjbUKL+UiMgaokHGBaZBxqUhnYLhwJELlIfLGBmPs6Wplr6Bcfbv3sSebY3URyOrXUwRkZKwlEHGCnAKTAFO8UtnBj95YRRgJgXDeGySJz9hC22tjWq5ERHJI82iEllh8USKw8f6OHp6kLFYgqmpqZkUDJdu30i0tlzBjYjIGqAAR2SB0t1Sx8+NkExOk0imGB6LU1YWonVznVIwiIisIcpFJbIA8UQKO97P8TMjVEfKGZuIAz4FQzgUYscldezdsVHTwUVE1ggFOCILcOLsCEdPDTEZT1IVCdPcWENsIkllRZidLQpuRETWGgU4InOIJ1KcHxil4+QAk/Ek8WSKeDJFbRVEa8rZ1lzHnu3KLyUistYUZYDjnGsCXgw8D7gC2ArEgZ8CnwI+ZWZTs1x3FfAW4BlAFdAJ/Atwm5mlcjzrd4EbgX1ACjgA/KOZfS3Pb0vWmHgixcGuPibjSSYmksQT/iNSUV5OKAStm6IKbkRE1qhiHWR8HfAJ4OnAD4EPAl8ELgc+CdzhnLtoKplz7oXAfcAzgS8BHwEqgQ8An5vtIc65fwRuB1qC5/0bPqD6qnPudfl+U7K2nDg7wpm+UR45eoHqqnIIBWNuwrCrpUHdUiIia1hRroPjnHs2UAt8PbOlxjm3BXgA2A68xMy+GOyvx7fWNABXm9mPgv1VwN3AzwG/bWafy7jXVcD9wFHgqWY2EOzfBTwYPP8yM+teTNm1Dk5xGJ9I8L0DJxmbSJCagq7TQ+zdvoHxWIJLdzYqBYOISAGtm3VwzOxuM/tqdjeUmZ0FPha8vDbj0EuAZuBz6eAmOH8C32UF8Jqsx9wQbP8hHdwE13TjW38iwO8t753IWjQWS9BxYoDxyRRjEwkqK0Ls29XERDzB5XuaFNyIiBSBogxw5pEItsmMfc8Ott+a5fz7gHHgKudc5tr6c13zzaxzpESMxRIcOdHPg4fPUR3x3U+xiSmSqRRP3N3ME3Y1KbgRESkCJRXgOOfKgVcFLzMDExdsj2RfY2ZJ4Bh+wHV7cJ9a/MDlUTM7M8ujOoLt3jwUW9aIsViCjp4BDnX1UVcb4cedvUSrIyRTU2zdHGVXa4PG3IiIFIminEU1h3fhBxp/w8zuzNjfEGyHclyX3t+4xPMv4py7Gfjb+Qora0c6v9SRnkHqoxEOdffT3tpA3+AYV162BbdTmcFFRIpJyQQ4zrnXAzcBh4FXLvLy9KClxY78nfV8M7sZuHnWBwWDjGXtiCdSdPQM8PCR8zTWVzM1PUV7kF9qX/smNjZWKrgRESkyJdFF5Zy7EbgVOAQ8y8z6s05Jt7g0MLv6rPPmO3++Fh4pEun8Uo929xOtjdDRM0CkopzqSDnt2zbQVF+j/FIiIkWo6AMc59wbgA8Dj+CDm7OznGbB9nFjZoJxO234QcldAGY2BpwCos65llnud2mwfdyYHikeY7EEB4/20tEzSLS6gmOnh2hrbeDkuWG2bY7yjMtbcLu0kJ+ISDEq6gDHOfcm/EJ9D+ODm/M5Tr072P7qLMeeCdQA3zezyQVe89ysc6TIxBMpDh/r4+jpQcZiCaampmhLd0u1baKxvoKaKnVLiYgUq6INcJxzb8UPKn4Q+EUz653j9C8AvcDLnHNPybhHFfD3wcuPZl2TXk/nzc65DRnX7MKnbpjEp4WQIpPuljp+boRkcppEMsXwWJyyshCtm+uoqa5Qt5SISJEr1pWMfxefQiEF3MbsY2G6zez2jGtehA90JvCpGfqBF+CnkH8BeKmZXVQZzrn3AW8ETgbnVAK/BTQBf2JmH15s2bWS8err7Bnk8PF+wmUheofGYdqnYKiqLGfPdqVgEBFZa5ayknGxBjg3M/807O+a2bVZ110NvBmfmiEz2eaH5km2+Tp8ss0p4CHgvUtNtqkAZ3WlUzAMjk5SFQkTqSgnNpGksiLMzpY6BTciImvQuglwipkCnNURT6Q43TvMwHCC7jPDjMb8cKuaSCXRmnK2b65TZnARkTVq3eSiElmsY6eGOX1+nAcfPTOTgqEiXM7U9BSbN9QquBERKTEKcKTkjU8kOHVhhEPHemmoq74oBcPOlnr27lBwIyJSahTgSElLZwYfm0jSUFdN1+kh2lsbGBiOsb99ozKDi4iUKI3BKTCNwSmcdH6pA0fO09ocZTKRpCwUZnh0gn27N7F3+wYFNyIiRUBjcEQCY7EER0744KauNsJPOnuZng4xPpFg/x4FNyIipU4BjpScsViCjp4BDnX1UVcboStIwdA/FGPvjkYu3abgRkSk1CnAkZKSzgx+5PgA9dHIzJiboZEYT9q7mb07Niq4ERFZBxTgSMnIzAxeURFmanqK9pYgv1T7JjY2Viq4ERFZJxTgSEmIJ1LY8X46e4aIVlfQ0TNApKKcqsowbVs30FRfo/xSIiLriGZRFZhmUa2Mzp5BHunqZTKepCoSJp6YZmgkxhW7m9m6pZZtzQ2rXUQREVmipcyiKl+x0ogUQDyR4vzAKB0nB5iMJ4knU8STKWqrKmndXEckUs7mxuhqF1NERApMAY4UrXgixcGuPibjSSYmksQTPl9qRXk5oRDsuETJM0VE1iuNwZGideLsCGf6Rnnk6AWqq8ohBBXlYcJh2NXSoOBGRGQdU4AjRWl8IsHRkwNMTCYvyi+VSk0rv5SIiGiQcaFpkPHyjcUSdPYM0HV6mNHYJPW1lZSFwozGJtnX1sQeLeQnIlJSlKpBSl46BcODh89RHfEtNLGJKZKpFE/c3cwTdjUpuBEREQU4UjyyUzCku6WSqSm2bo6yq7VB3VIiIgIowJEikc4MblkpGPoGx9jXtpH9bZvUciMiIjMU4Mial84v9fCR81SmUzC0KgWDiIjkpgBH1rTM/FLR2shMCobqSDnt25SCQUREZqdZVAWmWVQLl+6WOnlhFIBHu/tpC1punrKvhb07NlBTpZYbEZFSp1QNUjLiiRSHj/Vx7OwQyeQ01ZHwTHCzr20TjfUVCm5ERCQndVHJmpPuljp+boRkcppEMsXwWJyyshCtm+uoqa5Qt5SIiMxJAY6sOSfOjtB5cojqSDmJVDq/VJhwKMSOS+rY19ak6eAiIjIndVHJmpJOwTAWi1MVCdPcWENsIkllRZidLUqeKSIiC6MAR9aMdAqG8ckU8aT/UxOBaE052zfXsWe78kuJiMjCqItK1oTh0Ums++IUDBXhcqamp9i8oVbBjYiILIoCHFl1vuVmkIfs7EwKhrogBYMyg4uIyFIowJFVNRZL0H16kINdvTTUVc+kYOgfjrG/XSkYRERkaRTgyKpJL+T3kJ2naUMNx4LgZmgkxpP2bmbvjo0KbkREZEm0knGBaSVjL55IcehYHwcOnyNaG2FkPE5LUy29A+Ps372JPdsaqY9GVruYIiKyBixlJWO14EjBZeeXOnZ6iLqaSnoHx7nyss24nRsV3IiIyLIowJGCiidSHDnRT2fPINHqCo6dHppJwbC/fRO7WhvVLSUiIsumAEcKJp5I0XFigONnRhiNJZiamqKt5bH8UhsbKxXciIhIXmihPymIeCKFHe/nTO8o8USKRNL/qamqoKU5qvxSIiKSV2rBkYI4cXaEo6eGGI0lqK4qh5DPLzUNtG5SfikREckvteDIikvnl5qMJ4knU4TDIZobahifTLJrS50W8hMRkbxTgCMrKju/FEBsYoqy6in2t2+krbVRwY2IiOSduqhkxaQX8svOLxUOw7bNdQpuRERkxSjAkRUxFktw5EQ/B46cp642wk86e6mrUX4pEREpDK1kXGDrYSXjsViCjp4BDh7tpbKyfCa/1MjYJD/rlIJBREQWRysZy6pLd0vZ8QHqo5GZ4Eb5pUREpJAU4EjexBMpOnoGePjIeSorwkxNT80EN/vatZCfiIgUjgIcyYvs/FIdPQNEKsqpjpTTvm0DTfU1WshPREQKRmNwCqxUx+B09gxy+Hg/09PTPNrdP5Nf6in7Wti7YwM1VWq5ERGRpdEYHFkV6YX8xmJxpqanZoKbfW2baKyvUHAjIiIFp4X+ZFmyF/KLJ1PURCqVX0pERFaVWnBkyXIt5Dc1PaX8UiIisqoU4MiSZM6Y0kJ+IiKy1ijAkUXLnjF17PQQba0N9A/F2N++kf1tmzQdXEREVpUCHFmUeCKFHe+ns2eIaHXFTHCjhfxERGQtUYAji3Li7AhHTw0xNvH4GVNayE9ERNYKBTiyYOnp4JPxJPFkiuGxOOGyEK2b6zRjSkRE1hRNE5cFyZ4ODn7GVCgEOy6pY++OjRpULCIia4ZacGRe6RlTD2VNBw+HYVdLg4IbERFZcxTgyJyyZ0z9pLOXaLWmg4uIyNqmAEdyyjVjqm9wjH1tmg4uIiJrlwIcyUkzpkREpFgpwJFZacaUiIgUM82iksfRjCkRESl2asGRi+RKoKkZUyIiUkwU4MgMJdAUEZFSoQBHACXQFBGR0qIARwA/Y6rzpBJoiohIaVCAIzMzpsZimg4uIiKlQbOo1rnsGVPxZIqaSCUtzVFNBxcRkaKlFpx1LNeMqanpKVo31bGvrUmDikVEpCgpwFmnNGNKRERKmQKcdUgzpkREpNQpwFlnciXQ1IwpEREpJQpw1hkl0BQRkfVAAc46ogSaIiKyXmia+DqhBJoiIrKeqAVnHUjPmHpICTRFRGSdUIBT4rJnTP2ks5dotaaDi4hIaVOAU8JyzZjqGxxjX5umg4uISOnSGJwFcM5tA94O/CrQBJwBvgzcYmYDq1m2XOKJFEdO9HP05CATiRTVkbBmTImIyLoRmp6eXu0yrGnOud3A94HNwH8Bh4GnAc8CDLjazPoWer9QKDQNsJL1Hk+k6DgxQPeZYUbG4ySCQcW11ZVURcK0booqDYOIiBSNUCgEwPT0dGih16gFZ37/hA9uXm9mt6V3OufeD/wZ8A/ADatUtsdJ55c6PzhOPJGaCW4qysOaMSUiIuuGWnDm4JxrB44C3cBuM5vKOFaH76oKAZvNbGwh91ypFpx4IsXp3mHO9U/Qc3aYUChEbVUFF4bGqQiHKS8vY/fWRva3q+VGRESKy1JacDTIeG7PDrZ3ZQY3AGY2AtwP1ADPKHTB0uKJFMdOD/Cwnedsb4wTZ4dJJKcZm0gwmUiyqaGGyoow7ZoxJSIi64gCnLm5YHskx/GOYLu3AGWZVc/ZYc71j9M/MkHX6SFqqypJpHy3VGwyRSo1zf72jVy+p1mDikVEZN3QGJy5NQTboRzH0/sbM3c6524G/naFynSR8ckkA8NxkqlpqiMVJJJJNjfWEJtMUlkRZmdLHW2tjWq5ERGRdUUBzvKk+wIvGlBjZjcDN896QTAGJ19qIuXUVIWZTEwxPZkiWhMhHk+xsbGKbc31bL8kquBGRETWHQU4c0u30DTkOF6fdV7Bbd9STygMA8NxhkbLmJxM0ba1gR1b6qipUpeUiIisTwpw5mbBNtcYm0uDba4xOiuusiJMe+sGaF2tEoiIiKw9GmQ8t3uC7XOccxfVVTBN/GogBvxvoQsmIiIiuSnAmYOZHQXuAnYBN2YdvgWoBT6z0DVwREREpDC00N88ZknV8CjwdHyqhiPAVWstVYOIiEgpWcpCfwpwFsA5t53cyTb7F3MvBTgiIiKLowCnCCjAERERWRylahARERFBAY6IiIiUIAU4IiIiUnIU4IiIiEjJUYAjIiIiJUcBjoiIiJQcBTgiIiJSchTgiIiISMlRgCMiIiIlRwGOiIiIlBwFOCIiIlJyFOCIiIhIyVGAIyIiIiWnfLULsF6lM6OKiIhI/qkFR0REREpOaHp6erXLIHngnJs2MzUL5YnqM79Un/mnOs0v1Wd+rYX6VAuOiIiIlBwFOCIiIlJyFOCIiIhIyVGAIyIiIiVHAY6IiIiUHAU4IiIiUnIU4JSOW1a7ACVG9Zlfqs/8U53ml+ozv1a9PrUOjoiIiJQcteCIiIhIyVGAIyIiIiVHAY6IiIiUHGUTL2LOuW3A24FfBZqAM8CXgVvMbGA1y7aanHNNwIuB5wFXAFuBOPBT4FPAp8xsapbrrgLeAjwDqAI6gX8BbjOzVI5n/S5wI7APSAEHgH80s6/l+W2tKc65VwKfCV6+2sw+Ocs5qs95OOd+AXgDcBWwEejHf04/aGbfyDpX9TkH59zzgD/Fv9f0/4cPAu83sx/Mcv66r0/n3EuAa4AnAT8D1AGfNbNXzHHNitebc64a+CvgZcBOYBi4F/hbM3t0oe9PLThFyjm3G//D+3vAA8AHgC78D/gPgi/59eo64BPA04EfAh8EvghcDnwSuMM5d1ESOOfcC4H7gGcCXwI+AlTi6/Vzsz3EOfePwO1AS/C8f8MHVF91zr0u329qrXDObQduA0bnOEf1OQ/n3Ft4rI6+BbwP+CqwAbg261zV5xycc+8GvgZcia/LW4GHgBcC9zvnXpF1vurTewvwOnyAc2q+kwtRb865CPBt4G34wOZW4L/xv7T+yDn39IW+Oc2iKlLOuTuB5wCvN7PbMva/H/gz4ONmdsNqlW81OeeeDdQCX89sqXHObcEHg9uBl5jZF4P99fjfQhqAq83sR8H+KuBu4OeA3zazz2Xc6yrgfuAo8NR0i5lzbhc+8KwFLjOz7hV9swUWBIbfBtqA/wT+nKwWHNXn/Jxz1wF34P/j/g0zG8k6XmFmieDvqs85BD/Xp4ALwBPN7HzGsWfh6+iYmbUH+1SfgaB+TuLr4xrgHnK04BSq3pxzfw28A/gC8Fvp/8OD4OrLwCHgitla4bOpBacIOefa8cFNNz6CzvS3wBjwSudcbYGLtiaY2d1m9tXsHwAzOwt8LHh5bcahlwDNwOfSP7TB+RP433AAXpP1mHTw+A+Z3YHBD+pHgAi+da3UvB54Nv69jeU4R/U5B+dcGfBuYBz4nezgBiAd3ARUn3Pbif8u+2FmcANgZvcAI/j6S1N9BszsHjPrMLOFtHSseL0Fv0Clr/nLzP/Dzey/gO/hu7muWUB5FeAUqWcH27tm+RIfwUfMNfg+UrlY+osjmbEvXZ/fmuX8+/BfRFcFTacLueabWeeUBOfcE4B3Abea2X1znKr6nNtV+BawbwADzrnnOefe5Jz7U+fcz81yvupzbh34MXZPc85tyjzgnHsmflzJf2fsVn0uTSHqbTewAzhiZscWeE1OCnCKkwu2R3Ic7wi2ewtQlqLhnCsHXhW8zPyBy1mfZpYEjuEH5KebuGvxA5dHzezMLI8qufoP6u5fgRPA38x3erBVfc7uqcH2HH6cyNfwgeMHge87577rnMtscVB9zsHM+oE3AZcAh5xz/+yce6dz7g7gLnyX6h9nXKL6XJpC1Ftev9sU4BSnhmA7lON4en9jAcpSTN6FH2j8DTO7M2P/YutzPdb/24CfBa43s9g856o+57Y52N4AVAO/hG9luBy4Ez+A8z8yzld9zsPMPgj8Bv4L9tX4GTjXAT3A7VldV6rPpSlEveW1rhXglKb0DCGNIA84514P3AQcBl65yMuXWp8lUf/OuafhW23eN9t02yVY1/UJhINtCD/Y/TtmNmpmB/EzRU4C1+TorprNeq9PnHN/iR+Ueju+m6MWeDJ+ZulnnXPvWcTt1n19LlEh6m1Rz1CAU5zSUWxDjuP1Weeta865G/FTDQ8BzwqatDMttj7nO3++30KKRkbX1BHgrQu8TPU5t/Rgyy4z+3HmgaB1LN26+LRgq/qcg3PuWvyg7a+Y2RvNrMvMxs3sIXzAeAq4KZicAarPpSpEveX1u00BTnGyYJurH/LSYJurH3PdcM69Afgw8Ag+uDk7y2k56zP4gm/DD0ruAjCzMfx/mlHnXMss9yul+o/i6+UJwIRzbjr9Bz9jD+ATwb4PBq9Vn3NL189gjuPpAKg663zV5+yeH2zvyT5gZuP4pSHK8F2soPpcqkLUW16/2xTgFKf0D/JzgimnM5xzdcDVQAz430IXbC1xzr0JvwDVw/jg5nyOU+8Otr86y7Fn4mekfd/MJhd4zXOzzilmk8D/y/HnQHDO/wSv091Xqs+53Yf/IrjUOVc5y/HLg213sFV9zi09a6c5x/H0/niwVX0uTSHq7Sh+IsNe51zbAq/JSQFOETKzo/jZAbvwy19nugXf//yZIIJel5xzb8UPKn4Q+EUz653j9C8AvcDLnHNPybhHFfD3wcuPZl2TXk/nzc65DRnX7ML/m0zi00IUNTOLmdkfzvYH+Epw2qeDfZ8PXqs+5xB8Fj+Pb4Z/W+Yx59wvA7+Cb4JPz/RTfc7te8H2j5xzWzMPOOeei/+FbwL4frBb9bk0K15vwXo86Wvek/kxxx5+AAAMZUlEQVQLfLDQ3y/ghxp8dyEF1krGRcr5VA3fx8/I+C/gUXxqgmfhm++uMrO+1Svh6nE+78nt+JwntzF7f223md2ecc2L8D/AE/glx/uBF+CnLX4BeGn2YljOufcBb8QPCv0Cfsny38LnwfkTM/twPt/XWuOcuxnfTfW4XFSqz7k55zbj16vag/+CfgC/YN2L8QMof8fM/iPjfNVnDsGX4J342Wgj+BQCZ/Hdqs/HD0x9g5ndmnGN6pOZenhR8HILPrju4rGgsdfM/jzr/BWtt2Adnbvx60X9CPgOfm2c6/CtcM82sx8u5P0pwClizucEypVsM3sg7bqR8cU7l++a2bVZ110NvBm/5HhmErkP2dxJ5F6HX11zCr+uyXutSJPvLcZcAU5wXPU5B+fcRvwKsC/Grxcygu/ue6eZPa57WfWZm3OuAt8i8DL8e63Bf/k+gK+fu2a5Zt3X5wL+rzxuZruyrlnxenOPJdv8HXxwk5ls89DC3p0CHBERESlBGoMjIiIiJUcBjoiIiJQcBTgiIiJSchTgiIiISMlRgCMiIiIlRwGOiIiIlBwFOCIiIlJyyle7ACJrnXOuG7/KLMDzzezrOc57BNiPz3t1b0EKt0hB5uV7mGWhw1LjnPtD4DXAZfiF3wA2mFmuJJc4567HLx3/aTO7fqXLuJasp8+GrA8KcEQW553OuW+a2dRqF0Ryc849H/gEfkn5b+NXtYXHEi6uO865e4FrWMMBuEg+KcARWbhx4Arg5cC/rnJZZG7XBdvXm9knVrUkIrIqNAZHZOE+FGxvcc5VrmpJZD7bg23HqpZCRFaNWnBEFu6LwLOBpwE38FjAM6f5ugacc7cDvwv8XlaG85n9wP/hE6s+E6gGHgbeamb3BOc+H/gL4En47Mn3AjeZWc4veOdcLfBW4KX4ZI/n8clab86ViT5I8Prn+ASvO4AE8FN8d9CnZ8kkPPPe8Un23oTPer8R+A0z+3Ku8mXcowL4Y+CV+AzRFUA38F/AezITy2bUWdo9zrn0328xs5vne17GveqAtwEvAVqBXuArwJuznvk24Bbg42Z2Q457vSAo7/+Z2dOCfdcTjPfB1+nfAb8ObAJOAf8OvMPMxnPc83n4BIZPBeqBc/gszO8ys0czzrsWP7YmLbNOYJbPZVDnf4mv8134ZId3AX9jZidylGexn43G4BkvBNrwn9s+fPLGO83snavxDCkdasERWZy/CrZvds5FC/TMp+CzIu8FvgMYcDVwp3PuF5xzf4L/8gwBd+LHm/w6cJ9zrinHPSuDe70OeAT4Kj4z8OuAHzjnLsm+wDn3LPyXyevx/3d8C/gh8EQe+6LO5Tr8l+wO/JiY7+C/nObknKvCf7HeBlwO3BeUtREfLD3knGvPuOR/gnKcC17fGbz+ND4oXKgG4H7g94Pr7sIPVL4B+HYQAKR9HD+25+XOufoc93ttsP3ILMc24OvxpcH2TqAZn7H5O865muwLnHPvBL4GPAc4CHwBGAJeha+T52WcfpbcdfLp4HimCuCb+M96Z/D3KXzX7P8EQUN2eRb12Qje0/3AX+MDuv/GB9ed+IzTj8twXYhnSGlRC47IIpjZPc65u/BfLDfhf3NfaTfiW2Pen97hnHs3/jfTTwJbgGvN7HvBsXRQ8Av4L9a/m+WePwccAZyZnQquqwO+BPwiPqB4acbzWvAtWFHgeuAz6d+Wg9+qvwK80jl3d2YrVIbXAn9sZv+8yPf+duBa4DDwSxllrQb+DfgN4LPB+8HMPgl8Mmg5ugTfmnHvIp8J8CLgG8BVZjYaPLMV+F/gSnzdfDZ45jnn3B3AK/ABxoczb+Sc243/vPQBn5/lWS/AfxE/OT3DKwgwvw08A7gZ/2+dvt+v4YOPMeDXzOy+jGN/AbwH+Kxzbq+ZnTezw8D1i6iTq4AfAbvN7Hxw3wZ869CV+M/jP2Q8cymfjZfgg4yvAy8ys2TG/cL4Vr/MOlzxZ0jpUQuOyOL9NTAN3OScay7A836QGdwE3hVs9wIfSQc3AGY2AXwgePmsOe57UzpgCK4bwbdQpIDfDL440t6Ab2l4n5ld1BVgZj3Aq4OXf5LjWd9ebHATBDGvCV6+PqusMXy31RjwDOfc1Yu59wKMAn+QDm6CZ57mseDlF7POvy3YvobHey2+de1fgn+bbNPAazKnr5vZOeBPg5c3BEFr2k3B9tbM4Ca47r34Vo0GHvs3Waxp4PfTwU1w3yHg3cHL7Pe+lM9GuoXwvzMDj+CalJndvQrPkBKjAEdkkczsIeAOoA7fjbDSvjVLGQbwLQKzHuexwbWtOe45aGZfm+W+nfhWijL8eJ+0Xwu2/5Hjfg/ig4InZX0Zp/1njuvm8mT8b+ynzezbs5S1F99dBb6VJ58eNLPsrhvwLUmQVa9m9gC+G3FfMOYFmAnSrsd38Xwsx7N+YmY/zd4ZjK86hf+cPTm4Xzm+exLg9hz3+1SwvTbH8fmcmK085HjvLO2z8UCwfZNz7hWzdXutwjOkxKiLSmRp3gL8Jv636w+Y2fEVfNbJHPtHgaYcx9MtD7MFG+AH6ebSjf8S3ZaxLz3O5f+yBqjOpgn/xZxpKfWzNdgem+Oco1nn5susA2nxg21h9nr9EL7b7LX4Qd4AL8MPqP66mXXluOdc768b/97S/xZNQAQfMOWq0+XWyWLf+6I/G2b2Xefce/ADhv8VmHbOHcaPofqimd25Cs+QEqMAR2QJzKzTOfdJfJfO27l45s5izdeSOt+igiu16GDmjJRwsP08fvG8uUzOsi+2hOeHZilHrnPybSl1+h/A+4AXOedazOwMjw0u/qdlliddB6FZ9mVbbp0s9r0v6bNhZm9yzn0MP8Pp5/FB9auBVwfj3J6X0bVUiGdIiVGAI7J0b8cPKn2Fc+69c5yXXj0316yrnTn2r6RdCzh2OmNfD7AH+DszO7hCZcqWbplqm+Oc9LHsFqOCM7O4c+7j+Knlr3bOfQs/A66L2bsR03Yt4Fj636IX/wUeCY7NtgxAoetkyZ8NMzsGfDD4g3Pu5/HT45+Dn8GWHrdViGdIidEYHJElCn5DvxX/c/SOOU5Nf9Fcln0gmC1zZf5LN6/GYDZOdnna8TN3pvFTstO+GWyvy75mBaXHVWx1zmUPbCWYAv/rwct7C1iuuXwMP/39j/DTmQE+Ok9qj59xzl2evdM5dw2+m2kUXxcErQ33B6e8Ksf9rg+292btTwfa+f7FNm+fDTP7Hx4bW/QzBX6GlBgFOCLL824eW3cmV0vDd4LtjcF0VwCccxvxa3cUaj2dbO/LKk8U+Ci+O+BLWQu6vRc/BuNvnHM3BoNdL+Kce4ZzLm8BUDBTKj0w99asslYFZY0C/2tm989yi4ILgt4v4gOTl+O7U/5lnstCwD8FU7EBCGbn3Rq8/OegLtLSM+rekD17zDn3RvyU+SH8EgKZ0oH2Exb2bhZs0Z8N59yLnXPPdM6VZZ1XDfxS8DJzjFEhniElRl1UIstgZkPOuXfh1x553IJsgTuANwI/Cxx0zt2PX2jvqfiuhy/j110ppB/gA5kjzrm78b/dX4NfYO4ofq2TGWbW45x7EX5BuQ/jFzo8iJ/J1QrsDrafJ/dMl6V4K76b51qgIyhrDL/GTwt+QOzL8/i8fPgQfnAxwOcyVz3O4Sv4RQyPBmvVlOOn99fjV7B+W+bJZvb1YB2kN+EXc/we/nN0RXCfCeAVwVTzTF/Ct+681zn3y/iVqwHea2a22DeZUZ6lfDauwU+Dv+CcOwBcwE9tvwo/KPswfgHFgj1DSo9acESW7zZyz3TCzOL43xg/iv9y/hV8d9Wn8f/ZDhWgjNni+LQTH8evBPuCYN9HgGfMNkU6mLa8H98ddx7flfUi/OrEHfj1gfI6bT5YN+Y5+O6eQ/gv/hfif5t/D3DlHLOTVssDPDbjaCGDiwfwdfklfOvLc/Ff3O/Ap1EYy77AzP4K32r4bXxg8xL8OjH/il8wcLYlAL6CH/R8GP95/IPgT0v2uYu1hM/G7fjWzyP4oOw6fAqUTuDPgKcFa+8U9BlSWkLT03NNUBARkcVwzr0Q3yr3gJk9fY7zridIMWBm1xemdCLrh1pwRETyJBgb8tbgZfbq0yJSQBqDIyKyTM6538Ov/Pw0fP6jH+LHXonIKlELjojI8l2DH8Dbgg9sXpyZL0lECk9jcERERKTkqAVHRERESo4CHBERESk5CnBERESk5CjAERERkZKjAEdERERKzv8HDWPB99YjEqYAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(8,8))\n", "sns.scatterplot(number_of_tests, n_errors, alpha=0.5, ax=ax)\n", "ax.set_xlabel('Number of hypotheses')\n", "ax.set_ylabel('Number of errors')\n", "fig.tight_layout()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.7.3" } }, "nbformat": 4, "nbformat_minor": 4 }