minimisation_avec_conditions.ipynb

{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"name":"minimisation_avec_conditions.ipynb","provenance":[],"authorship_tag":"ABX9TyMdtQNFXx8d224ro4xF9eod"},"kernelspec":{"name":"python3","display_name":"Python 3"}},"cells":[{"cell_type":"code","metadata":{"id":"UTvNHvzYadqv","colab_type":"code","colab":{"base_uri":"https://localhost:8080/","height":1000},"executionInfo":{"status":"ok","timestamp":1593025993374,"user_tz":-120,"elapsed":19572,"user":{"displayName":"Clem Lo","photoUrl":"https://lh4.googleusercontent.com/-FzriA16UtQA/AAAAAAAAAAI/AAAAAAAAA2c/zi_BpyERii8/s64/photo.jpg","userId":"12213861745495664979"}},"outputId":"acce04b3-5847-4bbb-b701-1c07ea281978"},"source":["import numpy as np\n","import matplotlib.pyplot as plt\n","import math\n","\n","import keras\n","from keras import layers, models, optimizers\n","\n","from scipy import integrate\n","from scipy import interpolate\n","from scipy.stats import norm\n","\n","import time\n","\n","\n","\n","def calcul_energie (abscisse,fonction):\n","  hbar = 1\n","  omega = 1\n","  m = 1\n","  a = -5.\n","  b = 5.\n","\n","  #Calcul des interpolations\n","  tck_true = interpolate.splrep(abscisse, fonction, k=3, s=0)                        #F.O. fonction d'onde\n","  tck_true_carre = interpolate.splrep(abscisse, fonction*fonction, k=3, s=0)             #F.O. module carré\n","  tck_true_x = interpolate.splrep(abscisse, abscisse*abscisse*fonction*fonction, k=3, s=0)       #F.O. <x^2>\n","  der_true = interpolate.splev(abscisse, tck_true, der=1)                   #F.O. dérivée\n","  tck_true_der = interpolate.splrep(abscisse,der_true*der_true, k=3,s=0)         #F.O. dérivée spline 100\n","  int_true_carre = interpolate.splint(a,b,tck_true_carre)               #F.O. module carré\n","  int_true_x = interpolate.splint(a,b,tck_true_x)                       #F.O. <x^2>\n","  int_true_der = interpolate.splint(a,b,tck_true_der)                   #F.O. derivée carré\n","  #Calcul de l'énergie\n","  Energie = ((-pow(hbar,2)/(2*m))*(fonction[-1]*der_true[-1]-fonction[0]*der_true[0] \n","                             - int_true_der) + 0.5*m*omega*int_true_x ) / int_true_carre\n","  return Energie\n","\n","\n","def normalisation (abscisse,fonction):\n","  tck_true_carre = interpolate.splrep(abscisse, fonction*fonction, s=0)             #F.O. module carré\n","  int_true_carre = interpolate.splint(a,b,tck_true_carre)                           #F.O. module carré\n","  fonction = fonction*pow(1/int_true_carre,1/2)\n","\n","  return fonction\n","\n","\n","\n","\n","\n","hbar = 1\n","omega = 1\n","m = 1\n","pts = 100\n","a = -5.\n","b = 5.\n","x = a\n","h = 10/pts\n","linx = np.linspace (a,b,pts)\n","\n","\n","\n","\n","\n","#vraie fonction d'onde\n","vraie_onde = np.zeros_like(linx, dtype=float)\n","for i in range(0,pts):\n","  vraie_onde[i] = pow(m*omega/(math.pi*hbar),0.25)*math.exp(-m*omega*(pow(x,2))/(2*hbar))\n","  x+=h\n","#symétrisation\n","for j in range(0,pts):\n","    vraie_onde[j] = (vraie_onde[j]+vraie_onde[pts-1-j])/2\n","    vraie_onde[pts-1-j] = vraie_onde[j]\n","#normalisation\n","vraie_onde = normalisation(linx,vraie_onde)\n","#calcul de son énergie\n","energie_cible = calcul_energie(linx,vraie_onde)\n","\n","\n","\n","\n","\n","\n","\n","\n","#fonction d'onde constante\n","onde = np.ones_like(linx)\n","#normalisation\n","onde = normalisation (linx,onde)\n","#calcul de son énergie\n","energie_onde = calcul_energie(linx,onde)\n","#copie\n","premier_target = onde\n","\n","\n","\n","temps = time.clock()\n","\n","\n","#INITIALISATION DU MODEL\n","fits = 2000 #nombre de fits\n","\n","model = models.Sequential([\n","  layers.Dense(200, input_shape=(1,), activation='relu'),\n","  layers.Dense(200, input_shape=(1,), activation='relu'),\n","  layers.Dense(200, input_shape=(1,), activation='relu'),\n","  layers.Dense(1), # no activation -> linear function of the input\n","])\n","#model.summary()\n","opt = optimizers.Adam(learning_rate=0.001)\n","model.compile(loss='mse',optimizer=opt)\n","\n","\n","\n","\n","for i in range(0,fits):\n","\n","  #fit de l'onde\n","  model.fit(linx,onde,epochs=1,batch_size=50,verbose=0)\n","  predictions = model.predict(linx)\n","\n","  #symétrisation et valeur absolue de la prédiction + normalisation\n","  preds = np.abs(predictions.reshape(-1))\n","  for j in range(0,pts):\n","    preds[j] = (preds[j]+preds[pts-1-j])/2\n","    preds[pts-1-j] = preds[j]\n","  preds = normalisation(linx,preds)\n","\n","  energie_preds = calcul_energie(linx,preds)\n","\n","  #sélection de l'onde avec l'énergie la plus faible\n","  if (energie_preds < energie_onde):\n","    onde = preds\n","    energie_onde = energie_preds\n","    print('fit n°',i+1)\n","    print('Energie onde = ',energie_preds)\n","\n","\n","\n","keras.backend.clear_session()\n","\n","temps2 = time.clock()\n","print('Temps de calcul = ',temps2-temps)\n","print('Energie de référence = ',energie_cible)\n","print('Energie trouvée = ',energie_onde)\n","print('Energie 1er target = ',calcul_energie(linx,premier_target))\n","\n","\n","\n","\n","# Ecriture des fonctions d'onde\n","file1 = open(\"premier_target.txt\", \"w\")\n","file2 = open(\"fct_theorique.txt\", \"w\")\n","file3 = open(\"prediction.txt\", \"w\")\n","file4 = open(\"abscisses.txt\", \"w\")\n","\n","np.savetxt(file1, premier_target, fmt='%.18e', delimiter='', newline='\\n')\n","np.savetxt(file2, vraie_onde, fmt='%.18e', delimiter='', newline='\\n')\n","np.savetxt(file3, onde, fmt='%.18e', delimiter='', newline='\\n')\n","np.savetxt(file4, linx, fmt='%.18e', delimiter='', newline='\\n')\n","\n","file1.close()\n","file2.close()\n","file3.close()\n","file4.close()"],"execution_count":1,"outputs":[{"output_type":"stream","text":["Using TensorFlow backend.\n"],"name":"stderr"},{"output_type":"stream","text":["fit n° 16\n","Energie onde =  3.9015874440829883\n","fit n° 17\n","Energie onde =  3.7122208863514006\n","fit n° 18\n","Energie onde =  3.677658055774597\n","fit n° 19\n","Energie onde =  3.5311703694678553\n","fit n° 20\n","Energie onde =  3.421630483481135\n","fit n° 21\n","Energie onde =  3.3426752294156357\n","fit n° 22\n","Energie onde =  3.3302279623711257\n","fit n° 23\n","Energie onde =  3.316275754822528\n","fit n° 24\n","Energie onde =  3.2493029650232126\n","fit n° 25\n","Energie onde =  3.234876818311581\n","fit n° 27\n","Energie onde =  3.209267509415611\n","fit n° 28\n","Energie onde =  3.1903849944536655\n","fit n° 30\n","Energie onde =  3.180898135270433\n","fit n° 31\n","Energie onde =  3.1732017210500083\n","fit n° 33\n","Energie onde =  3.158191427371107\n","fit n° 34\n","Energie onde =  3.1575374949553034\n","fit n° 36\n","Energie onde =  3.14707368606539\n","fit n° 37\n","Energie onde =  3.145884317075071\n","fit n° 39\n","Energie onde =  3.1350433788938217\n","fit n° 42\n","Energie onde =  3.1307086384344527\n","fit n° 44\n","Energie onde =  3.127077340216017\n","fit n° 45\n","Energie onde =  3.1260912931640794\n","fit n° 46\n","Energie onde =  3.1253822394139688\n","fit n° 47\n","Energie onde =  3.1211813276330127\n","fit n° 50\n","Energie onde =  3.1196249698339997\n","fit n° 52\n","Energie onde =  3.1177426012257845\n","fit n° 55\n","Energie onde =  3.1163528026123\n","fit n° 57\n","Energie onde =  3.1154264370104277\n","fit n° 59\n","Energie onde =  3.114242002775449\n","fit n° 61\n","Energie onde =  3.1136424835684595\n","fit n° 62\n","Energie onde =  3.1128291939135404\n","fit n° 64\n","Energie onde =  3.1124718654288874\n","fit n° 65\n","Energie onde =  3.1118877019783984\n","fit n° 66\n","Energie onde =  3.1112786445513327\n","fit n° 68\n","Energie onde =  3.1096537931508332\n","fit n° 70\n","Energie onde =  3.107466653620848\n","fit n° 73\n","Energie onde =  3.106249316942368\n","fit n° 74\n","Energie onde =  3.1062405652835254\n","fit n° 75\n","Energie onde =  3.105080601582878\n","fit n° 77\n","Energie onde =  3.1034113503600245\n","fit n° 79\n","Energie onde =  3.102419391855016\n","fit n° 81\n","Energie onde =  3.102344819106891\n","fit n° 82\n","Energie onde =  3.1020636883233847\n","fit n° 83\n","Energie onde =  3.1009018166835416\n","fit n° 84\n","Energie onde =  3.1002451100415587\n","fit n° 87\n","Energie onde =  3.1000636378894164\n","fit n° 88\n","Energie onde =  3.099095072074958\n","fit n° 90\n","Energie onde =  3.0977707936421397\n","fit n° 92\n","Energie onde =  3.0967549010773157\n","fit n° 94\n","Energie onde =  3.0953833008869878\n","fit n° 96\n","Energie onde =  3.0932374054935496\n","fit n° 98\n","Energie onde =  3.0926066338748472\n","fit n° 99\n","Energie onde =  3.0921867530224367\n","fit n° 101\n","Energie onde =  3.0909257074485814\n","fit n° 103\n","Energie onde =  3.0898306224958523\n","fit n° 105\n","Energie onde =  3.0886419699806904\n","fit n° 107\n","Energie onde =  3.0834165720342623\n","fit n° 109\n","Energie onde =  3.0797841742640086\n","fit n° 111\n","Energie onde =  3.0774393094586565\n","fit n° 112\n","Energie onde =  3.0751090973452326\n","fit n° 113\n","Energie onde =  3.0750093350313934\n","fit n° 114\n","Energie onde =  3.073411963430556\n","fit n° 115\n","Energie onde =  3.073408542573314\n","fit n° 116\n","Energie onde =  3.073079797832711\n","fit n° 117\n","Energie onde =  3.070681384381006\n","fit n° 119\n","Energie onde =  3.0653057232182155\n","fit n° 121\n","Energie onde =  3.0597529381214694\n","fit n° 123\n","Energie onde =  3.0539356505619546\n","fit n° 124\n","Energie onde =  3.0538556414100193\n","fit n° 125\n","Energie onde =  3.051785495162744\n","fit n° 126\n","Energie onde =  3.050082364852266\n","fit n° 127\n","Energie onde =  3.0476951647501496\n","fit n° 129\n","Energie onde =  3.047189810061482\n","fit n° 131\n","Energie onde =  3.0461611196061322\n","fit n° 132\n","Energie onde =  3.0460302384515954\n","fit n° 134\n","Energie onde =  3.044586921525012\n","fit n° 136\n","Energie onde =  3.0434137088407955\n","fit n° 138\n","Energie onde =  3.0420401644253094\n","fit n° 140\n","Energie onde =  3.038923411908107\n","fit n° 142\n","Energie onde =  3.0373492184572486\n","fit n° 143\n","Energie onde =  3.036342792441633\n","fit n° 145\n","Energie onde =  3.0344472367105713\n","fit n° 146\n","Energie onde =  3.033468744097409\n","fit n° 148\n","Energie onde =  3.0308136983268645\n","fit n° 150\n","Energie onde =  3.0276722181523477\n","fit n° 152\n","Energie onde =  3.027511099498543\n","fit n° 153\n","Energie onde =  3.0228862783181274\n","fit n° 155\n","Energie onde =  3.021258306488593\n","fit n° 156\n","Energie onde =  3.018730049952688\n","fit n° 158\n","Energie onde =  3.017272454420846\n","fit n° 160\n","Energie onde =  3.016825307756345\n","fit n° 161\n","Energie onde =  3.015503394582041\n","fit n° 163\n","Energie onde =  3.0145065983743624\n","fit n° 164\n","Energie onde =  3.013305813617229\n","fit n° 166\n","Energie onde =  3.0110200538109635\n","fit n° 167\n","Energie onde =  3.009850050961394\n","fit n° 169\n","Energie onde =  3.007936920046121\n","fit n° 170\n","Energie onde =  3.006910688204242\n","fit n° 172\n","Energie onde =  3.006270933114825\n","fit n° 173\n","Energie onde =  3.002289704722312\n","fit n° 176\n","Energie onde =  2.998009244832981\n","fit n° 178\n","Energie onde =  2.9971212332644606\n","fit n° 179\n","Energie onde =  2.992780492684758\n","fit n° 181\n","Energie onde =  2.9921647021529565\n","fit n° 182\n","Energie onde =  2.9843759988720833\n","fit n° 184\n","Energie onde =  2.9797512030456694\n","fit n° 185\n","Energie onde =  2.978493508442766\n","fit n° 187\n","Energie onde =  2.9743713588240355\n","fit n° 190\n","Energie onde =  2.973531054409819\n","fit n° 194\n","Energie onde =  2.9720888013597264\n","fit n° 197\n","Energie onde =  2.971669553682301\n","fit n° 198\n","Energie onde =  2.970723871491838\n","fit n° 200\n","Energie onde =  2.969644355900347\n","fit n° 204\n","Energie onde =  2.969335574569853\n","fit n° 207\n","Energie onde =  2.967896034427703\n","fit n° 208\n","Energie onde =  2.967695772205718\n","fit n° 209\n","Energie onde =  2.966275937443244\n","fit n° 211\n","Energie onde =  2.9649766964894972\n","fit n° 212\n","Energie onde =  2.9629808533400945\n","fit n° 215\n","Energie onde =  2.961002600354939\n","fit n° 216\n","Energie onde =  2.960239309803819\n","fit n° 219\n","Energie onde =  2.9573234949689153\n","fit n° 222\n","Energie onde =  2.953564583867389\n","fit n° 223\n","Energie onde =  2.9503953508312666\n","fit n° 226\n","Energie onde =  2.9485594384759772\n","fit n° 227\n","Energie onde =  2.9476387729584013\n","fit n° 230\n","Energie onde =  2.946533806418372\n","fit n° 231\n","Energie onde =  2.9435021596316244\n","fit n° 232\n","Energie onde =  2.9421768142603995\n","fit n° 233\n","Energie onde =  2.9421206398904185\n","fit n° 234\n","Energie onde =  2.9418441159957878\n","fit n° 237\n","Energie onde =  2.9389651033329383\n","fit n° 242\n","Energie onde =  2.9339127802328755\n","fit n° 243\n","Energie onde =  2.927459842371439\n","fit n° 245\n","Energie onde =  2.9164763202337824\n","fit n° 249\n","Energie onde =  2.9113111327151757\n","fit n° 256\n","Energie onde =  2.9083090625637236\n","fit n° 257\n","Energie onde =  2.9021753662328575\n","fit n° 258\n","Energie onde =  2.8896589083668447\n","fit n° 261\n","Energie onde =  2.8792214273099193\n","fit n° 262\n","Energie onde =  2.8530249598096793\n","fit n° 267\n","Energie onde =  2.8489859293201762\n","fit n° 268\n","Energie onde =  2.841338567457019\n","fit n° 269\n","Energie onde =  2.841134042257722\n","fit n° 270\n","Energie onde =  2.8124588245052435\n","fit n° 273\n","Energie onde =  2.810954415656583\n","fit n° 274\n","Energie onde =  2.6462202683007083\n","fit n° 276\n","Energie onde =  2.623137282941326\n","fit n° 277\n","Energie onde =  2.444895584412523\n","fit n° 279\n","Energie onde =  2.3120410620446346\n","fit n° 281\n","Energie onde =  1.9578244681154386\n","fit n° 283\n","Energie onde =  1.8450443310750484\n","fit n° 284\n","Energie onde =  1.5689520216314283\n","fit n° 286\n","Energie onde =  1.5038625962909051\n","fit n° 287\n","Energie onde =  1.3224430745876747\n","fit n° 289\n","Energie onde =  1.2317433388480368\n","fit n° 292\n","Energie onde =  1.1884909160966532\n","fit n° 293\n","Energie onde =  1.1778187727305474\n","fit n° 294\n","Energie onde =  1.166410791487421\n","fit n° 295\n","Energie onde =  1.1205539225809913\n","fit n° 296\n","Energie onde =  1.115580467214097\n","fit n° 298\n","Energie onde =  1.098614798225955\n","fit n° 300\n","Energie onde =  1.0391703158718146\n","fit n° 302\n","Energie onde =  0.995655580670785\n","fit n° 303\n","Energie onde =  0.9405377171701804\n","fit n° 305\n","Energie onde =  0.8340278834076544\n","fit n° 307\n","Energie onde =  0.7759568976866213\n","fit n° 308\n","Energie onde =  0.7685886814093216\n","fit n° 310\n","Energie onde =  0.7518316771472968\n","fit n° 312\n","Energie onde =  0.7342829909489583\n","fit n° 314\n","Energie onde =  0.7325628211023646\n","fit n° 317\n","Energie onde =  0.7256947822338746\n","fit n° 322\n","Energie onde =  0.7209932309329855\n","fit n° 328\n","Energie onde =  0.7145578054256216\n","fit n° 335\n","Energie onde =  0.7136503081805985\n","fit n° 337\n","Energie onde =  0.7086123355557531\n","fit n° 339\n","Energie onde =  0.707068278796699\n","fit n° 341\n","Energie onde =  0.706697442305586\n","fit n° 344\n","Energie onde =  0.6976307224926723\n","fit n° 346\n","Energie onde =  0.6965241279861139\n","fit n° 350\n","Energie onde =  0.6964144427099095\n","fit n° 355\n","Energie onde =  0.6906898874159961\n","fit n° 360\n","Energie onde =  0.6895450631774404\n","fit n° 362\n","Energie onde =  0.6868916928386123\n","fit n° 367\n","Energie onde =  0.6857552278330586\n","fit n° 371\n","Energie onde =  0.6845285212907404\n","fit n° 373\n","Energie onde =  0.6821771610678049\n","fit n° 375\n","Energie onde =  0.6797938721938246\n","fit n° 392\n","Energie onde =  0.6780368685406877\n","fit n° 394\n","Energie onde =  0.6757013833348381\n","fit n° 397\n","Energie onde =  0.6742632876006855\n","fit n° 400\n","Energie onde =  0.6725486240005866\n","fit n° 402\n","Energie onde =  0.671932939379853\n","fit n° 404\n","Energie onde =  0.6708647348506047\n","fit n° 411\n","Energie onde =  0.6684553316655083\n","fit n° 416\n","Energie onde =  0.6641671439721983\n","fit n° 420\n","Energie onde =  0.6634624758305991\n","fit n° 426\n","Energie onde =  0.6598332547270627\n","fit n° 428\n","Energie onde =  0.6579450151434344\n","fit n° 430\n","Energie onde =  0.6558956741717166\n","fit n° 432\n","Energie onde =  0.6556328027845474\n","fit n° 435\n","Energie onde =  0.6548620915804853\n","fit n° 437\n","Energie onde =  0.6539054852217699\n","fit n° 444\n","Energie onde =  0.6535936195120534\n","fit n° 449\n","Energie onde =  0.6522361598742122\n","fit n° 453\n","Energie onde =  0.6508115598536088\n","fit n° 458\n","Energie onde =  0.6491396647019072\n","fit n° 460\n","Energie onde =  0.6490648694967571\n","fit n° 467\n","Energie onde =  0.6489306041466023\n","fit n° 469\n","Energie onde =  0.6489291273461338\n","fit n° 471\n","Energie onde =  0.6480808764777508\n","fit n° 473\n","Energie onde =  0.6480159536086206\n","fit n° 475\n","Energie onde =  0.6443962545552769\n","fit n° 480\n","Energie onde =  0.6434122582568893\n","fit n° 487\n","Energie onde =  0.6420476413933498\n","fit n° 492\n","Energie onde =  0.6379412313946887\n","fit n° 494\n","Energie onde =  0.6374332068885815\n","fit n° 507\n","Energie onde =  0.6342432353457952\n","fit n° 509\n","Energie onde =  0.6326348060659402\n","fit n° 517\n","Energie onde =  0.632100001697859\n","fit n° 519\n","Energie onde =  0.6320303858097441\n","fit n° 525\n","Energie onde =  0.631878393528258\n","fit n° 527\n","Energie onde =  0.6304767947585596\n","fit n° 529\n","Energie onde =  0.6262520410484339\n","fit n° 531\n","Energie onde =  0.6257634345510688\n","fit n° 533\n","Energie onde =  0.6238801678191912\n","fit n° 537\n","Energie onde =  0.6238019631470341\n","fit n° 543\n","Energie onde =  0.6197798682028896\n","fit n° 545\n","Energie onde =  0.6182204866959236\n","fit n° 547\n","Energie onde =  0.6171840597236022\n","fit n° 551\n","Energie onde =  0.6165266894230247\n","fit n° 553\n","Energie onde =  0.6121093293421384\n","fit n° 555\n","Energie onde =  0.6107919399032126\n","fit n° 558\n","Energie onde =  0.6098852842320669\n","fit n° 560\n","Energie onde =  0.6092544171274686\n","fit n° 565\n","Energie onde =  0.6079977397945543\n","fit n° 567\n","Energie onde =  0.6078480177086304\n","fit n° 569\n","Energie onde =  0.6034601210166376\n","fit n° 571\n","Energie onde =  0.6025974602625739\n","fit n° 574\n","Energie onde =  0.6024295515402149\n","fit n° 577\n","Energie onde =  0.601502253401076\n","fit n° 582\n","Energie onde =  0.6011861722203329\n","fit n° 584\n","Energie onde =  0.5996755673701935\n","fit n° 586\n","Energie onde =  0.5937927053383704\n","fit n° 588\n","Energie onde =  0.5850332891941353\n","fit n° 593\n","Energie onde =  0.5845049488319096\n","fit n° 597\n","Energie onde =  0.5830058601789319\n","fit n° 599\n","Energie onde =  0.5785311061489258\n","fit n° 609\n","Energie onde =  0.5752095968593153\n","fit n° 611\n","Energie onde =  0.5751254711007606\n","fit n° 613\n","Energie onde =  0.5707691254816988\n","fit n° 615\n","Energie onde =  0.5591153857888967\n","fit n° 617\n","Energie onde =  0.5556379528473553\n","fit n° 619\n","Energie onde =  0.5549086668795242\n","fit n° 622\n","Energie onde =  0.5537470270218507\n","fit n° 624\n","Energie onde =  0.5518077271231087\n","fit n° 631\n","Energie onde =  0.5461665596379341\n","fit n° 637\n","Energie onde =  0.5445507367962437\n","fit n° 639\n","Energie onde =  0.5326833421455747\n","fit n° 645\n","Energie onde =  0.5269677935707757\n","fit n° 651\n","Energie onde =  0.5267176101664827\n","fit n° 653\n","Energie onde =  0.5263548119277445\n","fit n° 654\n","Energie onde =  0.5263091379910991\n","fit n° 656\n","Energie onde =  0.5232832115088469\n","fit n° 658\n","Energie onde =  0.522975747743379\n","fit n° 660\n","Energie onde =  0.5204136267428232\n","fit n° 662\n","Energie onde =  0.5203852584463146\n","fit n° 666\n","Energie onde =  0.520177884087944\n","fit n° 672\n","Energie onde =  0.5200428609524442\n","fit n° 675\n","Energie onde =  0.5193426224069398\n","fit n° 690\n","Energie onde =  0.5172344491313445\n","fit n° 699\n","Energie onde =  0.5170696496299695\n","fit n° 704\n","Energie onde =  0.5170305936266218\n","fit n° 710\n","Energie onde =  0.5160997787196238\n","fit n° 712\n","Energie onde =  0.5159272018394314\n","fit n° 718\n","Energie onde =  0.5150536700476905\n","fit n° 720\n","Energie onde =  0.5137360661042598\n","fit n° 723\n","Energie onde =  0.5119050949235471\n","fit n° 729\n","Energie onde =  0.5118923549154093\n","fit n° 735\n","Energie onde =  0.5108098478316881\n","fit n° 741\n","Energie onde =  0.5108054241753466\n","fit n° 746\n","Energie onde =  0.5095059549282642\n","fit n° 748\n","Energie onde =  0.5084839214151424\n","fit n° 753\n","Energie onde =  0.5072443312175378\n","fit n° 755\n","Energie onde =  0.5070554764101992\n","fit n° 756\n","Energie onde =  0.5053749338655549\n","fit n° 758\n","Energie onde =  0.5048530037648685\n","fit n° 761\n","Energie onde =  0.504782211695959\n","fit n° 767\n","Energie onde =  0.5044980820737459\n","fit n° 778\n","Energie onde =  0.5041896514573672\n","fit n° 780\n","Energie onde =  0.5041050807017018\n","fit n° 793\n","Energie onde =  0.5040889178652069\n","fit n° 795\n","Energie onde =  0.5040875350796046\n","fit n° 796\n","Energie onde =  0.5040661168927804\n","fit n° 798\n","Energie onde =  0.5039143392061544\n","fit n° 813\n","Energie onde =  0.5037131777985288\n","fit n° 815\n","Energie onde =  0.5036869896883428\n","fit n° 830\n","Energie onde =  0.503540906940554\n","fit n° 837\n","Energie onde =  0.5034827249180369\n","fit n° 840\n","Energie onde =  0.5034496013741614\n","fit n° 842\n","Energie onde =  0.5034460130171069\n","fit n° 843\n","Energie onde =  0.5034112711806347\n","fit n° 845\n","Energie onde =  0.5033990335097674\n","fit n° 859\n","Energie onde =  0.5033798134598291\n","fit n° 861\n","Energie onde =  0.5033270938372741\n","fit n° 863\n","Energie onde =  0.5032993277171783\n","fit n° 873\n","Energie onde =  0.5032526483368034\n","fit n° 875\n","Energie onde =  0.5032123935186775\n","fit n° 883\n","Energie onde =  0.5030578295762986\n","fit n° 885\n","Energie onde =  0.502968886057015\n","fit n° 888\n","Energie onde =  0.5028938535625397\n","fit n° 890\n","Energie onde =  0.502805873118249\n","fit n° 897\n","Energie onde =  0.5027753689297936\n","fit n° 910\n","Energie onde =  0.5027646120510351\n","fit n° 912\n","Energie onde =  0.5027058592894771\n","fit n° 921\n","Energie onde =  0.5026729507612547\n","fit n° 928\n","Energie onde =  0.5026545007952814\n","fit n° 929\n","Energie onde =  0.5025963470390962\n","fit n° 966\n","Energie onde =  0.502449935453501\n","fit n° 967\n","Energie onde =  0.502392445398509\n","fit n° 975\n","Energie onde =  0.5023801799854897\n","fit n° 976\n","Energie onde =  0.5023353466018142\n","fit n° 986\n","Energie onde =  0.5021385042060145\n","fit n° 987\n","Energie onde =  0.5020998326416027\n","fit n° 989\n","Energie onde =  0.502053042354485\n","fit n° 991\n","Energie onde =  0.502001276240914\n","fit n° 992\n","Energie onde =  0.5019983941739676\n","fit n° 994\n","Energie onde =  0.5019947718413104\n","fit n° 1007\n","Energie onde =  0.5019743491397528\n","fit n° 1009\n","Energie onde =  0.5018466879849782\n","fit n° 1083\n","Energie onde =  0.5017821121334632\n","fit n° 1084\n","Energie onde =  0.5017468841845167\n","fit n° 1085\n","Energie onde =  0.5017226812415866\n","fit n° 1094\n","Energie onde =  0.5016952135212778\n","fit n° 1096\n","Energie onde =  0.5016662444865554\n","fit n° 1102\n","Energie onde =  0.5016457407937387\n","fit n° 1107\n","Energie onde =  0.5016226845252073\n","fit n° 1110\n","Energie onde =  0.5016193486749332\n","fit n° 1120\n","Energie onde =  0.5016165796067144\n","fit n° 1121\n","Energie onde =  0.5016142371886771\n","fit n° 1122\n","Energie onde =  0.5016125591407955\n","fit n° 1123\n","Energie onde =  0.5016111408288844\n","fit n° 1125\n","Energie onde =  0.5016106769182634\n","fit n° 1126\n","Energie onde =  0.5016092017272963\n","fit n° 1127\n","Energie onde =  0.5016087671194248\n","fit n° 1128\n","Energie onde =  0.5016070873495964\n","fit n° 1129\n","Energie onde =  0.5016059497120744\n","fit n° 1130\n","Energie onde =  0.5016047218686012\n","fit n° 1131\n","Energie onde =  0.5016034810648652\n","fit n° 1132\n","Energie onde =  0.5016020945905729\n","fit n° 1134\n","Energie onde =  0.5016020200242514\n","fit n° 1135\n","Energie onde =  0.5016013999002974\n","fit n° 1136\n","Energie onde =  0.5015998644219011\n","fit n° 1137\n","Energie onde =  0.5015993383577643\n","fit n° 1144\n","Energie onde =  0.5015993019395592\n","fit n° 1145\n","Energie onde =  0.5015992035045277\n","fit n° 1146\n","Energie onde =  0.5015990456488352\n","fit n° 1147\n","Energie onde =  0.501599023602696\n","fit n° 1178\n","Energie onde =  0.501598870515221\n","fit n° 1179\n","Energie onde =  0.501597727436315\n","fit n° 1180\n","Energie onde =  0.5015967317110149\n","fit n° 1181\n","Energie onde =  0.5015957526652263\n","fit n° 1182\n","Energie onde =  0.5015943945213942\n","fit n° 1183\n","Energie onde =  0.5015931638872166\n","fit n° 1184\n","Energie onde =  0.5015924152794522\n","fit n° 1185\n","Energie onde =  0.5015915863638253\n","fit n° 1186\n","Energie onde =  0.5015910837750401\n","fit n° 1187\n","Energie onde =  0.5015907898235651\n","fit n° 1188\n","Energie onde =  0.5015905004413921\n","fit n° 1189\n","Energie onde =  0.5015902434081472\n","fit n° 1190\n","Energie onde =  0.501590144618041\n","fit n° 1191\n","Energie onde =  0.5015898119582338\n","fit n° 1192\n","Energie onde =  0.5015895279914868\n","fit n° 1194\n","Energie onde =  0.5015893477901899\n","fit n° 1203\n","Energie onde =  0.5015892675190089\n","fit n° 1204\n","Energie onde =  0.5015890872432722\n","fit n° 1205\n","Energie onde =  0.5015890312767723\n","fit n° 1206\n","Energie onde =  0.501588995213116\n","fit n° 1209\n","Energie onde =  0.5015888730019151\n","fit n° 1210\n","Energie onde =  0.501588843244535\n","fit n° 1218\n","Energie onde =  0.5015886341826936\n","fit n° 1219\n","Energie onde =  0.5015885240062604\n","fit n° 1220\n","Energie onde =  0.5015884528494011\n","fit n° 1222\n","Energie onde =  0.501588414683313\n","fit n° 1223\n","Energie onde =  0.5015883733175797\n","fit n° 1224\n","Energie onde =  0.5015882744117707\n","fit n° 1226\n","Energie onde =  0.5015882392079286\n","fit n° 1227\n","Energie onde =  0.5015880330394944\n","fit n° 1229\n","Energie onde =  0.5015879449096806\n","fit n° 1230\n","Energie onde =  0.5015877108030622\n","fit n° 1231\n","Energie onde =  0.5015876473720333\n","fit n° 1232\n","Energie onde =  0.5015875635013598\n","fit n° 1234\n","Energie onde =  0.5015875321295628\n","fit n° 1240\n","Energie onde =  0.5015874172053386\n","fit n° 1241\n","Energie onde =  0.5015870971698404\n","fit n° 1242\n","Energie onde =  0.5015868703925747\n","fit n° 1243\n","Energie onde =  0.501586740647896\n","fit n° 1244\n","Energie onde =  0.5015865736704509\n","fit n° 1245\n","Energie onde =  0.5015862814584064\n","fit n° 1246\n","Energie onde =  0.5015859496281977\n","fit n° 1247\n","Energie onde =  0.5015858667973555\n","fit n° 1248\n","Energie onde =  0.5015858212114667\n","fit n° 1249\n","Energie onde =  0.5015856500969614\n","fit n° 1250\n","Energie onde =  0.5015855947335038\n","fit n° 1251\n","Energie onde =  0.5015855639292032\n","fit n° 1252\n","Energie onde =  0.5015855587692434\n","fit n° 1258\n","Energie onde =  0.5015852926987362\n","fit n° 1260\n","Energie onde =  0.5015852375888057\n","fit n° 1261\n","Energie onde =  0.5015848814616538\n","fit n° 1262\n","Energie onde =  0.5015848756329331\n","fit n° 1267\n","Energie onde =  0.5015847918690693\n","fit n° 1268\n","Energie onde =  0.5015845905271956\n","fit n° 1269\n","Energie onde =  0.5015843537936032\n","fit n° 1270\n","Energie onde =  0.5015843174565005\n","fit n° 1272\n","Energie onde =  0.5015841840081344\n","fit n° 1274\n","Energie onde =  0.5015841243758675\n","fit n° 1276\n","Energie onde =  0.5015840957766853\n","fit n° 1280\n","Energie onde =  0.5015840197518472\n","fit n° 1281\n","Energie onde =  0.5015839472434354\n","fit n° 1283\n","Energie onde =  0.5015839189904048\n","fit n° 1284\n","Energie onde =  0.5015838731979132\n","fit n° 1285\n","Energie onde =  0.5015837886274255\n","fit n° 1286\n","Energie onde =  0.5015837242677398\n","fit n° 1287\n","Energie onde =  0.5015836647180706\n","fit n° 1289\n","Energie onde =  0.5015836040863618\n","fit n° 1290\n","Energie onde =  0.5015835992161378\n","fit n° 1291\n","Energie onde =  0.5015835943297363\n","fit n° 1292\n","Energie onde =  0.5015834918866665\n","fit n° 1293\n","Energie onde =  0.5015834432430388\n","fit n° 1294\n","Energie onde =  0.5015833765877427\n","fit n° 1296\n","Energie onde =  0.5015833291392933\n","fit n° 1297\n","Energie onde =  0.5015832512467958\n","fit n° 1298\n","Energie onde =  0.5015831904260111\n","fit n° 1299\n","Energie onde =  0.5015830577356452\n","fit n° 1300\n","Energie onde =  0.5015828369678103\n","fit n° 1301\n","Energie onde =  0.501582658935986\n","fit n° 1302\n","Energie onde =  0.5015825742181229\n","fit n° 1303\n","Energie onde =  0.5015825494529896\n","fit n° 1314\n","Energie onde =  0.5015824655355409\n","fit n° 1315\n","Energie onde =  0.5015823835201232\n","fit n° 1316\n","Energie onde =  0.5015822790745765\n","fit n° 1317\n","Energie onde =  0.501582226919549\n","fit n° 1318\n","Energie onde =  0.501582090139958\n","fit n° 1319\n","Energie onde =  0.5015819485497294\n","fit n° 1320\n","Energie onde =  0.5015818596311056\n","fit n° 1321\n","Energie onde =  0.5015817411986746\n","fit n° 1322\n","Energie onde =  0.5015816604039216\n","fit n° 1323\n","Energie onde =  0.5015815473843825\n","fit n° 1324\n","Energie onde =  0.5015813989390897\n","fit n° 1325\n","Energie onde =  0.5015812119059765\n","fit n° 1326\n","Energie onde =  0.501580911404264\n","fit n° 1327\n","Energie onde =  0.5015807783284376\n","fit n° 1328\n","Energie onde =  0.5015805866000178\n","fit n° 1329\n","Energie onde =  0.5015804530928003\n","fit n° 1330\n","Energie onde =  0.5015803695427778\n","fit n° 1331\n","Energie onde =  0.5015802764874738\n","fit n° 1332\n","Energie onde =  0.5015802035773613\n","fit n° 1333\n","Energie onde =  0.5015801590147964\n","fit n° 1334\n","Energie onde =  0.5015801136247804\n","fit n° 1335\n","Energie onde =  0.5015799867493016\n","fit n° 1336\n","Energie onde =  0.5015798332999459\n","fit n° 1337\n","Energie onde =  0.5015797156056548\n","fit n° 1338\n","Energie onde =  0.5015796336690181\n","fit n° 1339\n","Energie onde =  0.5015795845764818\n","fit n° 1342\n","Energie onde =  0.5015794526330339\n","fit n° 1343\n","Energie onde =  0.5015792463768095\n","fit n° 1344\n","Energie onde =  0.5015790951116362\n","fit n° 1345\n","Energie onde =  0.5015790044569539\n","fit n° 1346\n","Energie onde =  0.501578796509979\n","fit n° 1347\n","Energie onde =  0.5015785807914231\n","fit n° 1348\n","Energie onde =  0.5015784708163322\n","fit n° 1349\n","Energie onde =  0.501578326893479\n","fit n° 1350\n","Energie onde =  0.5015783096103533\n","fit n° 1351\n","Energie onde =  0.5015782846077239\n","fit n° 1365\n","Energie onde =  0.5015782364232718\n","fit n° 1367\n","Energie onde =  0.5015781681714885\n","fit n° 1368\n","Energie onde =  0.5015780739802039\n","fit n° 1369\n","Energie onde =  0.5015780133371569\n","fit n° 1370\n","Energie onde =  0.5015779753126237\n","fit n° 1371\n","Energie onde =  0.501577929302186\n","fit n° 1372\n","Energie onde =  0.5015779232144151\n","fit n° 1380\n","Energie onde =  0.5015778528647528\n","fit n° 1381\n","Energie onde =  0.5015776475781498\n","fit n° 1382\n","Energie onde =  0.5015774684315467\n","fit n° 1383\n","Energie onde =  0.5015773056578148\n","fit n° 1384\n","Energie onde =  0.5015772117159036\n","fit n° 1385\n","Energie onde =  0.501577157516533\n","fit n° 1386\n","Energie onde =  0.50157704230758\n","fit n° 1387\n","Energie onde =  0.5015768707138718\n","fit n° 1388\n","Energie onde =  0.501576781108589\n","fit n° 1389\n","Energie onde =  0.5015767159234595\n","fit n° 1390\n","Energie onde =  0.501576667006648\n","fit n° 1391\n","Energie onde =  0.5015766087549335\n","fit n° 1392\n","Energie onde =  0.5015765981474284\n","fit n° 1393\n","Energie onde =  0.5015765062890697\n","fit n° 1394\n","Energie onde =  0.5015764225342401\n","fit n° 1395\n","Energie onde =  0.5015763254624136\n","fit n° 1396\n","Energie onde =  0.5015761947876018\n","fit n° 1397\n","Energie onde =  0.5015761416716026\n","fit n° 1398\n","Energie onde =  0.5015760320081688\n","fit n° 1399\n","Energie onde =  0.5015757794518346\n","fit n° 1400\n","Energie onde =  0.5015755834858244\n","fit n° 1401\n","Energie onde =  0.5015755545740076\n","fit n° 1402\n","Energie onde =  0.5015755358785335\n","fit n° 1403\n","Energie onde =  0.5015754789857062\n","fit n° 1408\n","Energie onde =  0.5015753880053905\n","fit n° 1409\n","Energie onde =  0.501575289193966\n","fit n° 1410\n","Energie onde =  0.5015752342503585\n","fit n° 1411\n","Energie onde =  0.5015750518777248\n","fit n° 1412\n","Energie onde =  0.5015749521376737\n","fit n° 1413\n","Energie onde =  0.5015748624429883\n","fit n° 1414\n","Energie onde =  0.5015748478013957\n","fit n° 1422\n","Energie onde =  0.501574785042758\n","fit n° 1424\n","Energie onde =  0.5015747756403427\n","fit n° 1425\n","Energie onde =  0.5015747579937069\n","fit n° 1426\n","Energie onde =  0.5015746825163705\n","fit n° 1427\n","Energie onde =  0.5015746764363489\n","fit n° 1428\n","Energie onde =  0.5015746314690986\n","fit n° 1434\n","Energie onde =  0.5015745713267923\n","fit n° 1435\n","Energie onde =  0.501574507403761\n","fit n° 1436\n","Energie onde =  0.5015744096843412\n","fit n° 1437\n","Energie onde =  0.5015742790620584\n","fit n° 1438\n","Energie onde =  0.5015741432116899\n","fit n° 1439\n","Energie onde =  0.5015740701851964\n","fit n° 1440\n","Energie onde =  0.5015740167385694\n","fit n° 1441\n","Energie onde =  0.5015740051194274\n","fit n° 1444\n","Energie onde =  0.5015739094794214\n","fit n° 1445\n","Energie onde =  0.501573893836312\n","fit n° 1454\n","Energie onde =  0.5015738181070328\n","fit n° 1455\n","Energie onde =  0.5015737773895084\n","fit n° 1463\n","Energie onde =  0.50157375661369\n","fit n° 1464\n","Energie onde =  0.5015735809121971\n","fit n° 1465\n","Energie onde =  0.5015735038654412\n","fit n° 1466\n","Energie onde =  0.5015734435334303\n","fit n° 1467\n","Energie onde =  0.5015734401117203\n","fit n° 1469\n","Energie onde =  0.5015734119940553\n","fit n° 1470\n","Energie onde =  0.5015733472493281\n","fit n° 1471\n","Energie onde =  0.5015732020985847\n","fit n° 1472\n","Energie onde =  0.5015731426039404\n","fit n° 1473\n","Energie onde =  0.5015729661829221\n","fit n° 1474\n","Energie onde =  0.5015728407647807\n","fit n° 1475\n","Energie onde =  0.5015728328538257\n","fit n° 1485\n","Energie onde =  0.501572756057095\n","fit n° 1486\n","Energie onde =  0.5015726545760217\n","fit n° 1487\n","Energie onde =  0.5015725647605417\n","fit n° 1488\n","Energie onde =  0.5015724677285278\n","fit n° 1489\n","Energie onde =  0.5015723983344624\n","fit n° 1490\n","Energie onde =  0.5015722984614804\n","fit n° 1491\n","Energie onde =  0.5015722416500716\n","fit n° 1492\n","Energie onde =  0.5015721757705518\n","fit n° 1493\n","Energie onde =  0.5015721418971094\n","fit n° 1494\n","Energie onde =  0.5015721250018085\n","fit n° 1495\n","Energie onde =  0.5015720803419744\n","fit n° 1496\n","Energie onde =  0.501572059647292\n","fit n° 1497\n","Energie onde =  0.501572002478405\n","fit n° 1498\n","Energie onde =  0.5015719706868902\n","fit n° 1499\n","Energie onde =  0.5015719287867552\n","fit n° 1500\n","Energie onde =  0.5015718974723353\n","fit n° 1501\n","Energie onde =  0.5015718593973599\n","fit n° 1502\n","Energie onde =  0.5015718131986043\n","fit n° 1503\n","Energie onde =  0.5015717926434681\n","fit n° 1504\n","Energie onde =  0.501571790958217\n","fit n° 1505\n","Energie onde =  0.5015717527027111\n","fit n° 1506\n","Energie onde =  0.5015716921412433\n","fit n° 1507\n","Energie onde =  0.5015716835909049\n","fit n° 1508\n","Energie onde =  0.5015715796231496\n","fit n° 1509\n","Energie onde =  0.5015715008742451\n","fit n° 1510\n","Energie onde =  0.50157130454354\n","fit n° 1511\n","Energie onde =  0.5015712153252553\n","fit n° 1512\n","Energie onde =  0.5015710178052794\n","fit n° 1513\n","Energie onde =  0.5015707369008043\n","fit n° 1514\n","Energie onde =  0.5015704057099276\n","fit n° 1515\n","Energie onde =  0.501570071233535\n","fit n° 1516\n","Energie onde =  0.5015696597622502\n","fit n° 1517\n","Energie onde =  0.5015690689384098\n","fit n° 1518\n","Energie onde =  0.5015685304315063\n","fit n° 1519\n","Energie onde =  0.5015678626908917\n","fit n° 1520\n","Energie onde =  0.5015669865347226\n","fit n° 1521\n","Energie onde =  0.5015661932014228\n","fit n° 1522\n","Energie onde =  0.5015658528229646\n","fit n° 1523\n","Energie onde =  0.5015653734260208\n","fit n° 1524\n","Energie onde =  0.5015647783544281\n","fit n° 1526\n","Energie onde =  0.5015643280545549\n","fit n° 1531\n","Energie onde =  0.5015636598592527\n","fit n° 1532\n","Energie onde =  0.5015624512313467\n","fit n° 1533\n","Energie onde =  0.5015618726546155\n","fit n° 1534\n","Energie onde =  0.5015614997349744\n","fit n° 1535\n","Energie onde =  0.5015612367460028\n","fit n° 1540\n","Energie onde =  0.5015611573821608\n","fit n° 1541\n","Energie onde =  0.5015609199175557\n","fit n° 1542\n","Energie onde =  0.5015606908226857\n","fit n° 1543\n","Energie onde =  0.5015604912502191\n","fit n° 1544\n","Energie onde =  0.5015604827342149\n","fit n° 1548\n","Energie onde =  0.5015603318771468\n","fit n° 1553\n","Energie onde =  0.5015599867839575\n","fit n° 1554\n","Energie onde =  0.5015595891733979\n","fit n° 1555\n","Energie onde =  0.5015593305691328\n","fit n° 1556\n","Energie onde =  0.5015592111538156\n","fit n° 1557\n","Energie onde =  0.5015591112531421\n","fit n° 1562\n","Energie onde =  0.5015588936312896\n","fit n° 1563\n","Energie onde =  0.5015587446426759\n","fit n° 1573\n","Energie onde =  0.5015584035253906\n","fit n° 1574\n","Energie onde =  0.5015577811939668\n","fit n° 1575\n","Energie onde =  0.5015573666520823\n","fit n° 1576\n","Energie onde =  0.5015570002281564\n","fit n° 1577\n","Energie onde =  0.5015566252887056\n","fit n° 1578\n","Energie onde =  0.5015563145554378\n","fit n° 1579\n","Energie onde =  0.5015561052457586\n","fit n° 1580\n","Energie onde =  0.5015559815868825\n","fit n° 1581\n","Energie onde =  0.5015559276125361\n","fit n° 1583\n","Energie onde =  0.5015558102623027\n","fit n° 1585\n","Energie onde =  0.5015556355467777\n","fit n° 1586\n","Energie onde =  0.5015556167278422\n","fit n° 1587\n","Energie onde =  0.5015556027089286\n","fit n° 1588\n","Energie onde =  0.5015555665349183\n","fit n° 1591\n","Energie onde =  0.5015555293412941\n","fit n° 1597\n","Energie onde =  0.5015554397311617\n","fit n° 1598\n","Energie onde =  0.5015553771299742\n","fit n° 1599\n","Energie onde =  0.5015552711135295\n","fit n° 1600\n","Energie onde =  0.5015551985079272\n","fit n° 1601\n","Energie onde =  0.5015551579011737\n","fit n° 1602\n","Energie onde =  0.501555099320512\n","fit n° 1603\n","Energie onde =  0.501554977802771\n","fit n° 1604\n","Energie onde =  0.501554904489553\n","fit n° 1605\n","Energie onde =  0.5015548334739528\n","fit n° 1606\n","Energie onde =  0.5015547936243826\n","fit n° 1608\n","Energie onde =  0.501554755649791\n","fit n° 1622\n","Energie onde =  0.5015544611687127\n","fit n° 1623\n","Energie onde =  0.5015538906879609\n","fit n° 1624\n","Energie onde =  0.5015535817335356\n","fit n° 1625\n","Energie onde =  0.5015532774244817\n","fit n° 1626\n","Energie onde =  0.5015532097804739\n","fit n° 1627\n","Energie onde =  0.5015531945426842\n","fit n° 1628\n","Energie onde =  0.5015530976806251\n","fit n° 1629\n","Energie onde =  0.5015530628854628\n","fit n° 1630\n","Energie onde =  0.5015528862763393\n","fit n° 1631\n","Energie onde =  0.5015527941078509\n","fit n° 1632\n","Energie onde =  0.5015526451991863\n","fit n° 1644\n","Energie onde =  0.5015520894871077\n","fit n° 1645\n","Energie onde =  0.501551347003493\n","fit n° 1646\n","Energie onde =  0.501550542940609\n","fit n° 1647\n","Energie onde =  0.5015500243245057\n","fit n° 1648\n","Energie onde =  0.5015496463295587\n","fit n° 1649\n","Energie onde =  0.5015494630461576\n","fit n° 1650\n","Energie onde =  0.501549180114197\n","fit n° 1652\n","Energie onde =  0.5015491655890941\n","fit n° 1653\n","Energie onde =  0.5015490957736323\n","fit n° 1654\n","Energie onde =  0.5015489454750192\n","fit n° 1658\n","Energie onde =  0.5015487810042605\n","fit n° 1660\n","Energie onde =  0.5015485704363327\n","fit n° 1661\n","Energie onde =  0.5015482253318004\n","fit n° 1662\n","Energie onde =  0.5015480282137703\n","fit n° 1663\n","Energie onde =  0.5015476843664992\n","fit n° 1665\n","Energie onde =  0.50154746281642\n","fit n° 1666\n","Energie onde =  0.5015473553033208\n","fit n° 1667\n","Energie onde =  0.5015472866937151\n","fit n° 1675\n","Energie onde =  0.5015472654860273\n","fit n° 1676\n","Energie onde =  0.5015471837430764\n","fit n° 1677\n","Energie onde =  0.5015470024148222\n","fit n° 1678\n","Energie onde =  0.5015467913372001\n","fit n° 1679\n","Energie onde =  0.501546701363234\n","fit n° 1681\n","Energie onde =  0.5015466779297437\n","fit n° 1682\n","Energie onde =  0.501546590078182\n","fit n° 1683\n","Energie onde =  0.501546475647443\n","fit n° 1684\n","Energie onde =  0.5015463364958963\n","fit n° 1685\n","Energie onde =  0.5015462086229727\n","fit n° 1694\n","Energie onde =  0.5015458317566798\n","fit n° 1695\n","Energie onde =  0.5015456389218044\n","fit n° 1696\n","Energie onde =  0.501545547665295\n","fit n° 1698\n","Energie onde =  0.5015452533749195\n","fit n° 1702\n","Energie onde =  0.5015452404306754\n","fit n° 1703\n","Energie onde =  0.5015447389505023\n","fit n° 1704\n","Energie onde =  0.5015444305316386\n","fit n° 1705\n","Energie onde =  0.5015441736414946\n","fit n° 1706\n","Energie onde =  0.5015441629685602\n","fit n° 1707\n","Energie onde =  0.501544111088001\n","fit n° 1708\n","Energie onde =  0.5015437321607574\n","fit n° 1709\n","Energie onde =  0.5015437075626918\n","fit n° 1710\n","Energie onde =  0.5015436934816145\n","fit n° 1711\n","Energie onde =  0.5015436295756801\n","fit n° 1712\n","Energie onde =  0.5015433522927033\n","fit n° 1713\n","Energie onde =  0.5015432383707694\n","fit n° 1714\n","Energie onde =  0.5015432007914221\n","fit n° 1715\n","Energie onde =  0.5015431182632258\n","fit n° 1716\n","Energie onde =  0.5015430187389411\n","fit n° 1717\n","Energie onde =  0.5015429373256511\n","fit n° 1718\n","Energie onde =  0.5015428735464902\n","fit n° 1719\n","Energie onde =  0.5015426673430281\n","fit n° 1720\n","Energie onde =  0.5015424105876978\n","fit n° 1721\n","Energie onde =  0.5015418948783552\n","fit n° 1722\n","Energie onde =  0.5015417183863238\n","fit n° 1723\n","Energie onde =  0.5015416299488835\n","fit n° 1724\n","Energie onde =  0.5015414452626604\n","fit n° 1725\n","Energie onde =  0.5015411659233718\n","fit n° 1731\n","Energie onde =  0.5015406201676686\n","fit n° 1732\n","Energie onde =  0.5015401950370804\n","fit n° 1733\n","Energie onde =  0.5015394876439443\n","fit n° 1736\n","Energie onde =  0.5015392477360249\n","fit n° 1737\n","Energie onde =  0.5015378638341752\n","fit n° 1742\n","Energie onde =  0.5015366632088825\n","fit n° 1743\n","Energie onde =  0.5015344279117715\n","fit n° 1744\n","Energie onde =  0.5015341803541111\n","fit n° 1745\n","Energie onde =  0.5015335375977136\n","fit n° 1754\n","Energie onde =  0.5015308389475881\n","fit n° 1755\n","Energie onde =  0.5015285842452295\n","fit n° 1756\n","Energie onde =  0.5015251934642307\n","fit n° 1757\n","Energie onde =  0.5015215696362345\n","fit n° 1758\n","Energie onde =  0.5015208200746218\n","fit n° 1759\n","Energie onde =  0.5015199316860646\n","fit n° 1767\n","Energie onde =  0.5015128353768856\n","fit n° 1768\n","Energie onde =  0.5015013988300234\n","fit n° 1773\n","Energie onde =  0.5014861662613884\n","fit n° 1774\n","Energie onde =  0.5014677093954958\n","fit n° 1775\n","Energie onde =  0.5014397204499867\n","fit n° 1776\n","Energie onde =  0.5014322002997773\n","fit n° 1777\n","Energie onde =  0.501418522454549\n","fit n° 1785\n","Energie onde =  0.5013603698598275\n","fit n° 1786\n","Energie onde =  0.5013286967345365\n","fit n° 1787\n","Energie onde =  0.5013087060539186\n","fit n° 1788\n","Energie onde =  0.5012866295749829\n","fit n° 1864\n","Energie onde =  0.5012686232522274\n","fit n° 1865\n","Energie onde =  0.5012236026177214\n","fit n° 1866\n","Energie onde =  0.5012062991082284\n","fit n° 1868\n","Energie onde =  0.5012026873646899\n","fit n° 1869\n","Energie onde =  0.5011916584018402\n","fit n° 1870\n","Energie onde =  0.5011779871019484\n","fit n° 1871\n","Energie onde =  0.5011754371962526\n","fit n° 1872\n","Energie onde =  0.5011725681263466\n","fit n° 1873\n","Energie onde =  0.5011716923279212\n","fit n° 1877\n","Energie onde =  0.5011679490563004\n","fit n° 1880\n","Energie onde =  0.5011639288213916\n","fit n° 1881\n","Energie onde =  0.5011616986519747\n","fit n° 1882\n","Energie onde =  0.5011602072428949\n","fit n° 1883\n","Energie onde =  0.5011594149107833\n","fit n° 1884\n","Energie onde =  0.5011593818072753\n","fit n° 1887\n","Energie onde =  0.5011575075896636\n","fit n° 1888\n","Energie onde =  0.5011553244919288\n","fit n° 1889\n","Energie onde =  0.5011545885025596\n","fit n° 1890\n","Energie onde =  0.5011540367682668\n","fit n° 1891\n","Energie onde =  0.50115375292032\n","fit n° 1893\n","Energie onde =  0.5011527110884663\n","fit n° 1894\n","Energie onde =  0.5011521346575111\n","fit n° 1895\n","Energie onde =  0.5011506583752329\n","fit n° 1896\n","Energie onde =  0.5011502840565969\n","fit n° 1897\n","Energie onde =  0.5011487745958437\n","fit n° 1898\n","Energie onde =  0.5011483371513381\n","fit n° 1899\n","Energie onde =  0.5011480791014351\n","fit n° 1900\n","Energie onde =  0.501147233390049\n","fit n° 1901\n","Energie onde =  0.5011470492562229\n","fit n° 1902\n","Energie onde =  0.5011467311528938\n","fit n° 1903\n","Energie onde =  0.501146459731648\n","fit n° 1906\n","Energie onde =  0.5011454978439709\n","fit n° 1908\n","Energie onde =  0.5011444855524094\n","fit n° 1909\n","Energie onde =  0.5011438559479092\n","fit n° 1910\n","Energie onde =  0.5011429458324247\n","fit n° 1912\n","Energie onde =  0.501141165709599\n","fit n° 1913\n","Energie onde =  0.5011405612885835\n","fit n° 1914\n","Energie onde =  0.5011402340742617\n","fit n° 1915\n","Energie onde =  0.5011400144686673\n","fit n° 1917\n","Energie onde =  0.501139591226944\n","fit n° 1918\n","Energie onde =  0.5011393777633452\n","fit n° 1919\n","Energie onde =  0.5011385578427918\n","fit n° 1920\n","Energie onde =  0.5011380431240998\n","fit n° 1921\n","Energie onde =  0.5011377388251826\n","fit n° 1923\n","Energie onde =  0.501137724008777\n","fit n° 1924\n","Energie onde =  0.501137532276514\n","fit n° 1925\n","Energie onde =  0.5011370009199303\n","fit n° 1930\n","Energie onde =  0.5011367465070158\n","fit n° 1931\n","Energie onde =  0.5011356979590093\n","fit n° 1932\n","Energie onde =  0.5011351191932998\n","fit n° 1933\n","Energie onde =  0.5011346043998247\n","fit n° 1934\n","Energie onde =  0.5011342223921068\n","fit n° 1936\n","Energie onde =  0.5011341836931145\n","fit n° 1937\n","Energie onde =  0.5011341631896051\n","fit n° 1938\n","Energie onde =  0.5011338259513161\n","fit n° 1939\n","Energie onde =  0.5011336412797537\n","fit n° 1940\n","Energie onde =  0.5011330322631689\n","fit n° 1941\n","Energie onde =  0.5011319511117989\n","fit n° 1942\n","Energie onde =  0.5011319450265238\n","fit n° 1943\n","Energie onde =  0.5011308938937271\n","fit n° 1946\n","Energie onde =  0.5011307693383626\n","fit n° 1947\n","Energie onde =  0.5011305808995419\n","fit n° 1952\n","Energie onde =  0.5011303820046856\n","fit n° 1958\n","Energie onde =  0.5011302958387501\n","fit n° 1967\n","Energie onde =  0.5011302876800996\n","fit n° 1968\n","Energie onde =  0.5011301667797778\n","fit n° 1970\n","Energie onde =  0.5011301514470478\n","fit n° 1978\n","Energie onde =  0.5011297768128378\n","fit n° 1981\n","Energie onde =  0.5011297461733822\n","fit n° 1982\n","Energie onde =  0.5011288708046489\n","fit n° 1983\n","Energie onde =  0.5011284156415855\n","fit n° 1984\n","Energie onde =  0.5011272785304777\n","fit n° 1985\n","Energie onde =  0.501126790775864\n","fit n° 1986\n","Energie onde =  0.5011266253948583\n","fit n° 1995\n","Energie onde =  0.5011254548167212\n","fit n° 1996\n","Energie onde =  0.5011249861761952\n","fit n° 1997\n","Energie onde =  0.501123733719408\n","fit n° 1999\n","Energie onde =  0.5011235063413906\n","fit n° 2000\n","Energie onde =  0.5011226213393575\n","Temps de calcul =  21.868633\n","Energie de référence =  0.5001266854789138\n","Energie trouvée =  0.5011226213393575\n","Energie 1er target =  4.166666666666669\n"],"name":"stdout"}]},{"cell_type":"code","metadata":{"id":"aIJDuTYTbH4b","colab_type":"code","colab":{"base_uri":"https://localhost:8080/","height":278},"executionInfo":{"status":"ok","timestamp":1593026195639,"user_tz":-120,"elapsed":1012,"user":{"displayName":"Clem Lo","photoUrl":"https://lh4.googleusercontent.com/-FzriA16UtQA/AAAAAAAAAAI/AAAAAAAAA2c/zi_BpyERii8/s64/photo.jpg","userId":"12213861745495664979"}},"outputId":"c5f21630-ac2e-4001-c58e-779d863ec7da"},"source":["\"\"\"\n","import numpy as np\n","import matplotlib.pyplot as plt\n","\n","premier_target = np.loadtxt(\"premier_target.txt\", dtype=float, delimiter=' ')\n","vraie_onde = np.loadtxt(\"fct_theorique.txt\", dtype=float, delimiter=' ')\n","onde = np.loadtxt(\"prediction.txt\", dtype=float, delimiter=' ')\n","linx = np.loadtxt(\"abscisses.txt\", dtype=float, delimiter=' ')\n","\"\"\"\n","plt.xlabel('x')\n","plt.ylabel('y')\n","#plt.yscale('log')\n","plt.plot(linx,premier_target,c='darkgrey',label = \"première cible (constante)\",alpha = 0.5)\n","plt.plot(linx,vraie_onde,c='darkorange',label = \"fonction d'onde théorique\")\n","plt.plot(linx,onde,c='deepskyblue',marker = 'x',label = 'prédiction après 2000 epochs',linestyle='None')\n","plt.legend(title = \"Minimisation énergie (avec conditions)\" ,loc = 'lower right', prop={'size': 10})\n","plt.tick_params(axis='both',labelsize=8)\n","plt.savefig('minimisation_avec_conditions.pdf')"],"execution_count":6,"outputs":[{"output_type":"display_data","data":{"image/png":"iVBORw0KGgoAAAANSUhEUgAAAX4AAAEFCAYAAADkP4z+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeVxV1f7/8ddiEgdAVDRTURxSkEkUw8zh5FRppqU0qGk20i3t5zetbjlcrCzrZtqAdW+Z2oiW1r11y+k4lRVqVmo5MIgKKsosINPn98eBLbOoHA7Dej4ePNxnnz2svcF19tl7rfdSIoKmaZrWeNjZugCapmla7dIVv6ZpWiOjK35N07RGRlf8mqZpjYyu+DVN0xoZB1sXoDratGkjXbp0sXUxNE3T6pU9e/acFRGPsvPrRcXfpUsXdu/ebetiaJqm1StKqWMVzde3ejRN0xoZXfFrmqY1Mrri1zRNa2R0xa9pmtbI6Ipf0zStkdEVv6ZdhsXxYE4pPW1OsUxD6WlNq6t0xa9plyHYBUIPWir4YBcYvx/G7bdMm1Ms7wW72LqUmla1etGOX9NsaXG8pTI3uVt+In0sFX4/FxBAAS8eg6gMWO9rWQYsHwRRGTDH05al17Ty9BW/pl1Cyav8YrkCm1NhZkeY0dEynV9iaAt99a/VZfqKX9MuofgqP/QghF0LS0+Ak4KnOlumFTC3aHr8fssHQUSCZZ3iq39Nq0v0Fb+mVYPJ3VLpLzwG+SKsK3wX05HFqPwsJD8LU87PrPe1fBNYeMyyrK70tbpKX/FrWjWYUyAiQRjmkEjUheYQs4aoljexLud+yM8h6mRPglt9j5PL89zQ0o6IBDC11JW/Vjfpil/TKlDyga45BUIPFPJsynLy0+J4rtkFQnt9R6Svg6ViL8iFnyIJzR7FuqO3ExX4Mrd69ib04MXbPfpBr1aX6Fs9mlaBkg90o9ILeTb9AxY5TyTYexCm8UuI9HUgKqNoYXsnojpNJrLnBUyOpwje838sisvjWU9LZa8f9Gp1jRKRSy9lY/369RMdy6zVtuIKO6xgKxEXehPZ8gdMweOqXun8afioH+ZmAwn1+oSwDnb6Qa9mM0qpPSLSr+x8fcWvaZUwuUOY858sLBxKmPoJU7/bL71S83Ywbj2m5K8IS12tH/RqdZKu+DWtEubjCUSkeDA34wMinMdgTlXVW7FdX8xDviKiya3MZRMRCaX7AGiaremKX9MqYE6B0CMtiIybQvjQkUT2VuU6cVW57vmRRMoKwn+9lUjPk9VeV9Nqg674NY3S4WsAUSeP8eyJ+UT1nAUuHY1OXMYD3SpEZRTd079hCovbzYE//l1qXR3kptmarvg1jTKxDCIEH3qNRdc+T3CvG41lTO7Va445x7Ponn6L9gR7difU8TE4d5A5nrqFj1Y3WLXiV0otUUrtUEotLTP/M6XUVqXULqXUPmuWQdOqo2Qsw7y9fxHqOo9I9x8xeTS9uu32n0Bk4qOEHmvPvBgp1bZf02zFahW/UioIaCEigwAnpVRw8XsicreIDAUWA/+1Vhk07XKY3CHsmgIWZngTdn4tpj63XP1GnVpgCryZsNNvsjBe6RY+Wp1gzSv+EGBj0fQmYEAFy4wHvrRiGTSt2swpEHEij7kJ4US0mo45rWY6tps7Tiei3ePMTV1ORILoh7yazVmz4m8JpBdNpxW9NiilHAE/Edlb0cpKqYeVUruVUruTkpKsWExNK773LkSenkF4biSRvk410hLHnAKhf9oT6f4z4dFhRHr8plv4aDZnzYo/DXAtmnYFUsu8PxTYWtnKIvKeiPQTkX4eHh5WKaCmFYvKgEiP/ZiO/wv6PIGplap2K55LbtcHTH5Dwbk1pj/Da2S7mnY1rFnx7wKGFU0PB34q8/54YJ0V969p1TbHE0yHFkKTluAzGah+K55LbtcdcGwK/g9D9FeY7I/psDbNpqxW8RfdwslRSu0ACoB4pdRzAEopheWe/05r7V/TLqVU2/2ME3DkS8y+i1ic2Nw6+2s3B3MLE+x7x5in2/RrtmDV5pwiMlNEBonIEyJySkReLJovItJHRAqtuX9Nq0qptvu/RWBuPphQh4es1sY+2KMlod2/xBz9J+Rl6Tb9ms3oPH6t0TLa7h8QwhLdieixjsje9lZrbmlyh8iOcYQWvE/Ynr+IKAjSbfo1m9A9d7VGzeQOYU5/sNDjKcJaJVu9Ejb18CMsax0Lc4J0m37NZnTFrzVq5hSIyPRk7rllRJzvYvVmluZURUTLKZa+AifydbNOzSZ0xa81WpYhFQuIPHon4e3OEOlT/QTOK97fQYjsmUP4qXAiZYVu06/ZhK74tUYrKgMiHb/ElLEFvCddVgLnFe/PB0zXuEOXUZj+Wkikd6Fu06/VOj30ota4fVQ0Kt3kWv77+vMT+HYS3LUNOg6u3X1rjYYeelHTyjr3J5zeA96Ta3/f3W8Hx+Zw8KPa37fW6OmKX2u8/vwYlB30urv29+3YHHrcAYcjIT+n9vevNWq64tcaFaO3roil4vccjjn3mlrvPbs4HsxdnoALaRD7LaB78Wq1R1f8WqNi9NY9egDS4zB3/z+b9J4NdoHQ5H6Y29wBBz/SvXi1WqUrfq1RMXrrnujCvA4vEpo53Ca9Zy3lUIR6rmRedj9CD4juxavVGl3xa42OqaUQlvw+C6/5O2Ed7GxW2ZrcIazVOUs5nP/Ulb5Wa3TFrzU65tjDRLjey9ym+4hIsF0HKqPXcNJrRGR01B25tFqjK36tUTGnQGh8RyJj7ybcr6MxwHptV7pGL14fRbjLQSLjJhF6QA/LqNUOXfFrjUpUBkSeehxTy0Jo1sbqvXWrLEfxPf0e4zEl/5fI1rt1L16tVuhYZq1RmdPiMJz8EExLjXkm99p/qFpqBK7OI8CxOaaTH2AaHly7BdEaJX3FrzUuR9db/u1+u23LUZKDM3jdAke/Aj02kVYLdMWvNS5H10HbIHDtbOuSlNZ9PJxPhMSfbV0SrRHQFb/WeGQmQOJP0GO8rUtSXtfRYOcIR9bZuiRaI2DVil8ptUQptUMptbTM/FZKqUil1JbiAdg1zVqMmIbory0zuo+vU/EIi+PBnOUGnUyWbyQidap8WsNjtYpfKRUEtBCRQYCTUqrkU6v5wDwRual4AHZNsxYjpiH+GLh5YbbzqVPxCEb5OodB6lHMx0/UqfJpDY81r/hDgI1F05uAASXe8wX+rpQyK6UGlFsTUEo9rJTarZTanZSUZMViag2dyR0ie+YS2nQW87q+Teifqk7FIxgxEhduY177fxAa27pOlU9reKxZ8bcE0oum04peF7sBWATcDbxa0coi8p6I9BORfh4eHlYsptYYmDLNhJ2JYKHdLXVykHOTO4R1sGfhtfMIO7+uzpVPa1isWfGnAa5F065Aaon3DovInyJyGtDt1zSrM8ccJqJtGHM75tk0pqEy5hSISIC5bCTCaSTmpGxbF0lrwKxZ8e8ChhVNDwd+KvHeYaVUe6VUc3QnMs3KzCkQKpOIzHqd8O6ONotpqMzF+AYI72pPZHQooX/a1ZnyaQ2P1Sp+EdkL5CildgAFQHyJFjzzgU+BLcAL1iqDpgFEnTlL5NE7MXXqCGCzmIbKlIpv6HAjpgu7icx7p86UT2t49GDrWsO3dymYn4QHoqFlV1uX5tK+usMyFvBDcaCUrUuj1WN6sHWt8Yr5Flr1qh+VPoDXrZARD+cO2rokWgOlK36tYcs7Dye2WbJw6guvmy3/xv7PtuXQGixd8WsNW7wZCi5YrqLrC5eO4OFvDMKuaTVNV/xag2TENMT+DxybQ4dB9SIGwSi3161wcgfkZtSLcmv1i674tQbJiEE4nQadTJgzmtSLGASj3O3ugsJ8zEf+qBfl1uoX3YZea5BM7hDpmUBo1hLCmscScZB6EYNgxDccCCCs4yIizvgRGVj3y63VL/qKX2uwTMn/ISwpgoW5/etkTENlLPENioXtniEsdVW9KbdWf+iKX2uwzCcSiGj7N+Z6Sp2MaaiMEd/g9DMRLUIxnzxl6yJpDYyu+LUGyZxcQGiTJ4gs/IDwrqrOxTRUplR8Q08XImNCCT3qVufLrdUvuuLXGqSohBNERk/E1NkyxGJdi2moTKn4hlbemDhC5PnFdb7cWv2iH+5qDdKcjJWQuQ081xrzTO51/z7/HM8SL5SCziMxRS/FdPPzgL2tiqU1MPqKX2uYjm2Edn2haWtbl+TqdB4BOSlwZq+tS6I1ILri1xqeC+mQsAu6jLR1Sa5e5+GWf+M22LYcWoOiK36t4TluBimwXC3Xd808oG0fOKYrfq3m6IpfazCMuINjGy0xDe0H1Ou4A+N4Oo+EhB91fINWY3TFrzUYF2MaMqDj4HoT01AZ43jahlriG44eqNfHo9UdulWP1mCY3CHSK4nQrNcIaxZdb2IaKnMxvqEPYR1eJOKML5EB9fd4tLpDX/FrDYop5VtLTENeSL2KaaiMEd9wzd8JS/u43h+PVjdYteJXSi1RSu1QSi0tM/9DpdTPSqmtSql7rVkGrXExHz9JRNvH6l1MQ2WM+AaHnUQ0uwPzqXp+QFqdYLWKXykVBLQQkUGAk1IquMwik0RkqIh8Yq0yaI2LOVkIdXiEyLx36lVMQ2VKxTd0c7LENxxqWm+PR6s7rHnFHwJsLJreBAwo8Z4Aq5RS/1FKda5oZaXUw0qp3Uqp3UlJSVYsptZQRJ05Q2T0BEwdrgXqT0xDZUrFN7TriynvVyJz36q3x6PVHdas+FsC6UXTaUWvi/2fiNwAvAL8s6KVReQ9EeknIv08PDysWEytoZiTFYkpcyt0HmbMM7mXiUGoR+Z4lnhGYWcPnUyY4t6ut8ej1R3WrPjTANeiaVcgtfgNEUku+ncncI0Vy6A1JvGbwc3L8tMQeQ6D9DhIjbF1SbR6zpoV/y6g+NJrOPBT8RtKKdeif3tS4gNB065YYT6c2Aqew21dEuspPrb4zbYth1bvWa3iF5G9QI5SagdQAMQrpZ4revtjpdRO4N/AM9Yqg9aInN4DF9IsV8UNVaue0OJaOLbJ1iXR6jmrduASkZllZr1YNP82a+5Xa4SKr4I9b7JtOaxJKctVf+y3IIWgdDcc7crovxytXjPybOI3g0cANPNokHk2xnF6DoPss5D0e4M8Tq126Ipfq9eCXSD0gGBOcwDPYUbb94aWZ2Pk9rjfAoA55kiDPE6tduisHq1eM7lDZJtfCe38EWGuqfU+n6cyRm7PQQ/Cur5FRNZIIvs0vOPUaoe+4tfqPdPZLwg7+y4L03s0iHyeypjcIexaWOj+N8KSIjC55tm6SFo9pSt+rd4zn0ohou3fmNuZBpHPUxkjt6fFn0S0egBzzCFbF0mrp3TFr9Vr5jOZhLr9g0in9YR7Ue/zeSpTKrend1siY+4i9KRXgztOrXboil+r16JOxhMZE4qpi6W3bn3P56lMqdyepq0xNUslMuX5BnecWu3QD3e1em1O6nuQ8zO0v5gBaHJvePf5y+XzeA7D9OsyTGNeBJrZokhaPaav+LX6LX4zdLgRHJrYuiS1y3MYFOTCyR9sXRKtHtIVv1Z/nT8NZ/dDpwbcW7cyHW4EOwc4vsXWJdHqIV3xa/XXcbPl384NOJ+nMk4toH2IDmzTroiu+LV6p1RMQxM3aBvU6OILFseDueMDlnC6HEvAbWM7B9qV0xW/Vu8Y8QVJ56HjUMxp9o0uviDYBUIL78XcfDAc39pgoyo069CterR6x+QOkZ0TCc1aSliL2AYb01AVkztE9laE5kYSdnw3EWca3znQrpy+4tfqJVPyN4QlRbDwQv8GHdNQFVNrR8JyN7LQ7pZGew60K6Mrfq1eMp9MsMQ0eEqDjmmoijkFIpqNY25COBEnCxrlOdCujK74tXrHnFxIqONjROa/S3hX1WBjGqpiRDh4niQ8cT6RLTY3unOgXTld8Wv1TtTp00RGT8TUqQPQcGMaqmJEOHh1A+fWmE590ujOgXblrPpwVym1BOgH7C07DKNSqikQC0wWET2IqFZtc85/CplbwXO1Ma8hxjRU5WKEg51luMljmzCNEkzuypbF0uoJq13xK6WCgBYiMghwUkoFl1nkQeAPa+1fa8DiN4P7deDS0dYlqRs8h0HmSUg5bOuSaPWENW/1hAAbi6Y3AUaKllLKqej9SoNGlFIPK6V2K6V2JyUlWbGYWr1SkAcntlsGHdcsOhedC92LV6sma1b8LYH0oum0otfFpgEfVbWyiLwnIv1EpJ+Hh4d1SqjVP6d+gbzMxhnTUBm3ruDaWVf8WrVZs+JPA1yLpl2BVACllAMwSkT+Z8V9aw2MEdNwbBOgLD12dUQBAIuPK8ydH4P4LVBYAOj4Bq1q1qz4dwHFl2XDgZ+KptsBnkqp74DJwCKlVCN6LKddCSOm4VQytAvCnNNKRxQUCXaB0CYzMDsGwplfdXyDdklWa9UjInuVUjlKqR3APiBeKfWciLwIBAMopRYAO0VEtz7WqmRyh8jrsgjNep6wVr83ypiGyhjnJjeSsKOHiFD63GhVu2TFr5R6AvjoSirnsk04gRfLvL/gcrepNV6m9C2EJe1l4bXzmKsjCkoxtW9F2C/vstDxEeZ21udGq1p1bvW0A6KUUpFKqZuVUrqhsGYT5mOxRLQNY26n/EYb01AZcwpEuE1i7qkXiTgp+txoVbpkxS8izwM9gPextMY5opR6SSnVzcpl0zSDOQVCZRKR518lvJtDo4xpqIwR39DuT8JPPk9k6yh9brQqVevhrogIcKroJx9wB9YqpRZbsWyaZog6m0rk0TsxdWgHNM6YhsoY8Q3dfcG+CabTn+tzo1WpOvf4ZwL3AWeBfwOzRSRPKWUHHAHmWLeImgZzLqy3xDR0XmrMa2wxDZW5GN/Q1DIW77ENmIbqc6NVrjpX/K2AO0RklIisEZE8ABEpBMZYtXSaVuzYRmjWDtr42bokdVvnEZYB6DMTbV0SrQ6rzj3++SJyrJL3/qz5ImlaGVII8Zss0QS6bUHVOo+0/Buvcw+1yulYZq3uS/odss5Yrma1qrUNgKYelm9ImlYJXfFrddbFmIaiSsxzuI4iuITFx+0wez1hOWcigI5v0MrTFb9WZxkxDQlnoLUP5vwOOorgEoJdILTZU5hVLzi7X8c3aBWy6kAsmnY1LFEEFwjNnkOY63Ad01ANJneI7J5hiW84fJQI9DnTytNX/FqdZso0E3YmgoVqFGE6pqFaTB3aEnZ+LQvzBuhzplVIV/xanWaOO2qJaeiYp2MaqsmcAhEtp+j4Bq1SuuLX6ixzCoQW3kvk+dcI7+6oYxqqwYhvuOaQJb6h1c/6nGnl6Ipfq7OizpyzxDR06gDomIbqKBXf4NAMU+JH+pxp5eiHu1qdNSc70hLT0OVdY56OaajaxfiGJuB5E8R9h2mYPmdaaQ264k9MTCQnJ8fWxdCuUNsDX+DUrBMnUh0gLdbWxal3XFz60Sbmvxz/w0x+iy62Lo52BZydnWnfvn2Nb1ff6tHqpoILNE3aRfY1Q3RMwxXKbjcEgGant9m4JFpd06Cv+K3xSanVkmOboSALV/9QXL28bF2aesoLfulO6/QoWnstsHVhtDrEqlf8SqklSqkdSqmlZeYvVUptU0r9rJQaaM0yaPWLEdMQ9x3YO0Enk44cuEKL48HsNQOOmyHfcstTn0sNrFjxK6WCgBYiMghwUkoFl3j7KREZAoQCf7dWGbT6x4hpSDwHHQZhPt9CRw5coWAXCHV4GLPz9XBih45v0AzWvNUTAhRHBG4CBgBRAMWZ/kAL4DcrlkGrZ0zuENnlNKFZrxDW7KiOabgKJneI9BZLfEPsXiKc9LnULKx5q6clkF40nVb02qCUWgdswPKhUI5S6mGl1G6l1O6kpCQrFlOra0xn1hKWFKEjB2qAycOZsLxNlsiL9qLPpQZYt+JPA1yLpl2B1JJvish44HrgpYpWFpH3RKSfiPTz8PCwYjG1usYcH09Eu8eZ6yk6puEqmVMgotk45iaEE3GyQJ9LDbBuxb8LGFY0PRz4qfgNpVSToslM4LwVy6DVM+akbEKbP0Wk/WeEd1U6puEqGPEN3TMIT5xPpMNafS41wIoVv4jsBXKUUjuAAiBeKfVc0dufK6W2Av8B5lurDFr9E3U8msjoUExdvQEd03A1jPiGDm3BIxDTsbf1udQAUFI0Sk9d1q9fP9m9e7eti6HVhu+mw9EvISwJ7B1tXZqG44d58POLEHYGmra2dWm0WqKU2iMi/crO1z13tbpDCiH2G+hyi670a1q324rO77e2LolWB+iKX6s7En+xDKre7TZbl6ThadcXml8D0f+xdUm0OkBX/JrNGb11Y/4Dyh663Kx7mNawxcftMHedZekRXZAL6F68jZmu+DWbM3rrnjwNHQZizmmle5jWsGAXCHWegdmpr+7FqzXskDatfjC5Q2TnU4RmLdK9da3E0ou30NKLN2aP7sXbyOkrfq1OMJ36RPfWtTKTR1PC8rewUN2se/E2crri1+oE87FYS2/dzujeulZiToGIprcX9eLN1+e4EdMVv2Zz5oSzhLrOI9L5v4R7oXvrWoHRi7dXHuGnFxJZ8IE+x42Yrvg1m4uKP0xkTCimXpZ+Jrq3bs0zevG2cwHPYZiOLibSW/Q5bqT0w13N5uYcnwtOp6C1jzFPD6pesy4Owg70uBM2Powp/zdMnoE2K5NmO/qKX7OtrLNwfJulMtJqR/dxoOzgyBe2LolmI7ri12zC6LQV/RVIAfS4U3coqiWLz3pYhmQ8fLHi1+e+cdEVv2YTRqet2CPg5oXZMVB3KKolwS4Q2uolzLnt4NyfujNXI6TTOTWbMZ/JJPS3bMIc9xHhOKLSDkV5eXmcOHGCnJyc2i9kA5VTUEBSLrioPDKUMx6O4KwvA+stZ2dnOnbsiKNj6XDDytI59cNdzWYsQyzGs/DaecytotPWiRMncHFxoUuXLiilareQDdjJ1LMk2rehtxN0aHLp5bW6SUQ4d+4cJ06cwMvLq1rr6M94zWbMR/dXa4jFnJwcWrdurSv9GpSeD0n2rWifm0BSrpCeb+sSaVdKKUXr1q0v6xuxrvg1mzAnnCO0+dNEOn1VrSEWdaVfc9LzISYHujoX0iEvga6cJSYHXfnXY5f7/0NX/JpNRMUdtHTa6n0DoDtt1YRvvvmG33///ZLLZRVAV2dwdXSAJm645iTS1VnIKqiFQmp1glUrfqXUEqXUDqXU0jLz31VK/aCU2qmU8rdmGbS6aU7MLEzNM6BVT2Oeyb1MR6MqKKWYPHmy8To/Px8PDw/GjBkDwNdff83LL79c5TYSEhKYMGHCZZV73rx5bNq06bLWAdi6dSs//vij8Xr58uWsWrXqsrdTme+++45t27bh5+d3yWWvaQKuxU/3nFtBQS6uhZlcU837/Fdy3rKzsxkyZAgFBQ3j06VFixZA6XOxb98+vv324ghn1fkbvFzDhw8nJaUGcjZExCo/QBDwr6LpCCC4xHteRf/2AL641Lb69u0rWgNy7pDIa4hE/bNaix88eLDcvObNm0tAQIBkZWWJiMi3334rAQEBMnr06Botak2ZP3++vPrqq7YuRnkF+SKn9oikxZWanZeXV6O7eeutt+SNN96o0W3aUvPmzcvNW7Fihfztb3+z6n4//PBDeeGFFyp8r6L/J8BuqaBOteYVfwiwsWh6EzCgxIdNbNFkHtAwLgG0KhkdtgD++gRQmK+976o6Dd1666188803AHz66afcc889xnsffvghjz/+OADTpk1jxowZ3HDDDXTt2pW1a9cCEBcXh6+vr7H8uHHjGDFiBF26dOGtt97i9ddfp0+fPoSEhJCcnGxsq3j9Z555Bh8fH/z9/XnqqacA+M9//sP1119Pnz59GD58OKdPnyYuLo7ly5ezZMkSAgMD2bFjBwsWLOC1114DLFeKISEh+Pv7M378eOOKbujQoTz99NP079+f6667jh07dlR4Hl599VWCg4Px9/dn/vz5xrF5e3vz0EMP0bt3b0aOHEl2djYA0dHR3HzzzfTt25frBw1ld3wm5CQzbdpUHn30Ufr1v56//d8coqOjCQkJwc/Pj+eff964yi153goKCpg9e7ax/3fffbfCMn788cfcfvvtAGRmZjJs2DCCgoLw8/Pjq6++Ms7n22+/baxT8hxVdIwAq1atwt/fn4CAAKZMmVJuv5mZmdx///34+fnh7+/PF198Yfy9+Pn54evry9NPP20s36JFC5577jkCAgIICQnh9OnTAMTGxjJgwADjXBQrPhe5ubnMmzePzz//nMDAQD7//PNSf4NxcXHcdNNN+Pv7M2zYMOLj442/p4r+NhMTExk8eDCBgYH4+voav/uxY8fy6aefVniOL0tFnwY18QP8Hbi5aHo4MK+CZT4Fbqxk/YeB3cBuT0/PK/+I1OqELckibXaKbDlXKPLv7rJl3SzL6+RLr1vZFf9vv/0md955p2RnZ0tAQICYzWbjir/k1dfUqVNlwoQJUlBQIAcOHJBu3bqJiEhsbKz07t3bWL5bt26Snp4uZ86cEVdXV4mIiBARkSeffFKWLFlibGvNmjVy9uxZue6666SwsFBERFJSUkREJDk52Zj3r3/9S2bNmiUi5a/4S7728/OTrVu3iojI3LlzZebMmSIiMmTIEGP9b775RoYNG1buPHz//ffy0EMPSWFhoRQUFMjo0aNl27ZtEhsbK/b29vLrr7+KiMjEiRNl9erVIiJy0003yeHDh0VEZPPOnyR4sEnSzvwlUyffI6NuHS27U/MlLU9k9OjR8sknn4iISEREhHGVW/K8vfvuu7Jw4UIREcnJyZG+fftKTExMqTJeuHBB2rVrZ7zOy8uTtLQ0ERFJSkqSbt26SWFhoezdu1cGDx5sLOft7S3x8fGVHuP+/fulR48ekpSUJCIi586dK3d+5syZY5zP4t/PyZMnpVOnTnLmzBnJy8sTk8kk69atExERQL7++msREZk9e7ZxbLfddpusXLlSRCzfXio6F2Wv+Eu+HjNmjHz44YciIvL+++/L7bffLiKV/22+9tprxpV9fn6+pKenG9vt3r27nD17ttyxXs4VvzXb8acBrrxPMmsAACAASURBVEXTrkBqyTeVUk8CB0VkZ0Uri8h7wHtg6cBlxXJqtaD44W3o/nzCmk4hwn32VY8A5e/vT1xcHJ9++im33nprlcuOGzcOOzs7fHx8jKu4cmU0mXBxccHFxQU3Nzduu80y6Lufn1+5h6Zubm44OzvzwAMPMGbMGOPZwokTJ7jrrrtITEwkNzf3ku2q09LSSE1NZciQIQBMnTqViRMnGu/fcccdAPTt25e4uLhy62/YsIENGzbQp08fwHKFe+TIETw9PfHy8iIwMLDU+pmZmfz444+l9pF94QIxTbqSVSDcMHYiPZrb4+oAu3btYv369QDce++9xreasvv//fffjSvVtLQ0jhw5Uuq4z549S8uWLY3XIsLf//53tm/fjp2dHSdPnuT06dP06dOHM2fOkJCQQFJSEu7u7nTq1ImlS5dWeIy//fYbEydOpE2bNgC0atWqXPk2bdrEZ599Zrx2d3dn+/btDB06FA8PDwAmTZrE9u3bGTduHE5OTsbvsm/fvmzcaLlp8cMPPxjfFqZMmVLqW0J17Nq1iy+//NJYf86cOcZ7Ff1tBgcHM336dPLy8hg3bpzxewRo27YtCQkJtG7d+rLKUJI1K/5dwCNAJJYr/g+L31BKjQRuAO6y4v61OsbkDmG53xd12LpQI+mbY8eO5amnnmLr1q2cO3eu0uWaNLn45FIq6a1echk7OzvjtZ2dHfn5pds6Ojg48Msvv7B582bWrl3LW2+9xZYtW3jiiSeYNWsWY8eOZevWrSxYsOAqju5imezt7cuVofhYnn32WR555JFS8+Pi4kodj729PdnZ2RQWFtKyZUv27dtXavmT6WlkqSa0c2l68cFvNYgIb775JqNGjap0maZNm5ZqY/7xxx+TlJTEnj17cHR0pEuXLsb7EydOZO3atZw6dYq77rqrymN88803q1/QanJ0dDSaRpY959ZqUlzR3+bgwYPZvn0733zzDdOmTWPWrFncd999gKVfS9OmTa9qn1a7xy8ie4EcpdQOLPfx45VSzxW9/SbgBZiVUhXfFNQaHHNSNhESwtz8r4k406RGBgGZPn068+fPr1ZrlpqUmZlJWloat956K0uWLOG3334DLFe8HTp0AGDlypXG8i4uLmRklG+r6ubmhru7u3EPd/Xq1cbVf3WMGjWKDz74gMzMTABOnjzJmTNnKl3e1dUVLy8v1qxZA1gqmh/2/EaScqVZYRbpBRfb84eEhBhXuSWvmsvuPyIigry8PAAOHz7M+fPnSy3j7u5OQUGBUbmnpaXRtm1bHB0dMZvNHDt2zFj2rrvu4rPPPmPt2rXGt5LKjvGmm25izZo1xgd+8XOYkkaMGFHquUFKSgr9+/dn27ZtnD17loKCAj799NNLnvOBAwca5+Djjz+ucJnKfscAN9xwQ6n1Bw0aVOX+jh07Rrt27XjooYd48MEH2bt3L2D5fZ06dYouXbpUuf6lWLU5p4jMFJFBIvKEiJwSkReL5vcUkWARGSoij1xqO1r9Z4wAFT2RcB+PGhtlq2PHjsyYMaNmCnkZMjIyGDNmDP7+/tx44428/vrrgOWB5MSJE+nbt69xCwLgtttuY926dcbD3ZJWrlzJ7Nmz8ff3Z9++fcybN6/a5Rg5ciT33nuv8eBxwoQJlVY+xT7++GPef/99AgIC8O7dm4/WfUXXpopm5NEuP4mYHEtP3jfeeIPXX38df39/jh49ipubW7ltPfjgg/j4+BAUFISvry+PPPJIhd9MRo4cyc6dlru6kyZNYvfu3fj5+bFq1Sp69eplLNe7d28yMjLo0KED7du3r/IYe/fuzXPPPceQIUMICAhg1qxZ5fb7/PPPk5KSgq+vLwEBAZjNZtq3b8/LL7+MyWQiICCAvn37Gg+eK7N06VLefvtt/Pz8OHnyZIXLmEwmDh48aDzcLenNN99kxYoV+Pv7s3r1apYuXVrhNopt3bqVgIAA+vTpw+eff87MmTMB2LNnDyEhITg4XOXNmopu/Ne1H92cs/575ZjIljVhIit8RIoefm5Jtsy/lIoeWmk1IzFHJK245eb5JJHEKEnLzpLEHJHz588bD6o//fRTGTt27BXvZ8+ePTJ58uQaKHHjNmPGDNm0aVOF79WVh7uaZpjT9A84FgFDl0DRvVI9ypbtleq05ewOGfG4XjiFq5sXO3bs4fHHH0dEaNmyJR988MEV7ycoKAiTyURBQQH29vZXX/BGytfXl2HDhl31dnTFr1nN4nhLxrvJHfjjX2DvhLnD/UTFV7+HrlZ7TuXZ06zptbhmJ4BLJwYNGsSOPb+RVUC1e/VWZfr06Ve/kUbuoYceqpHt6KwezWqMwVaSLsDB1Zh7Pk9otJse8KOOamYPMaot6XbNISfZCHNrpi/QGxx9xa9ZjdF2/3chrNX/I6L501fddl+zHlcH6NpUESPd8MhNIalQ6OqsLqt5p1Y/6Ct+zapMLYWw1FUsvHYeYR0ddaVfx7k6KDzsckm098DDLldX+g2Urvg1qzIf+Y2IZuOZ67yXiERVI233NetJz4ckaUr7vFMk5dvrjP4GSlf8mtWYUyD0hBeRJx8ivK93jbXd1y7asGGD0XnsahkDtDRVdHAqoOuFaGKy9ehcDZGu+DWriTpzlsgj4zH18APHpnqwlRpQNgv/xhtv5J///CdHjx696m0bA7Q4AE3b4lqYifnfr/LhSsu4AUOHDmX37t3l1iuZQlldv/76Kw888MBVl/lSymbkX4mXXnqpWsvVWFZ+LdAVv1ajSsYvzznxD0xZOzF7PWnEL1/OYCsN0dUORHLttdcagWgAzZo1Y9WqVXTv3r3C5SvqRVuZkgO0nMp3JN35Wv7fvcOZcf+9ABQInM298rKX9NJLL9VKj+varPinTJnCO++8c1X7qi264tdqlNGE83QmHFiBuXc4obGta7UJZ2JiIrGxsTX6k5iYWOU+4+Li6NWrF5MmTcLb25sJEyaQlZUFQJcuXXj66acJCgpizZo1bNiwgQEDBhAUFMTEiRONDJouXbrw7LPPEhgYSL9+/di7dy+jRo2iW7duLF++3NhPRVn4fn5+Rhb+1q1bGTRoEGPHjsXHx6famfkls+2fnD6FGLtrePa1D3jt5XDS8yG7EL76dLWREf/LL7+U20ZSUhJ33nknwcHBBAcH88MPP5RbJiMjg99//52AgACg5jLz16xZY0QzDB48uMKM/F9++YUBAwbQp08fbrjhBg4dOgRYvrXccccd3HzzzfTo0cNIz3zmmWfIzs4mMDCQSZMmAfDRRx/Rv39/AgMDeeSRR4wP8xrLyq8NFXXnrWs/OrKhftmSLNLGfF7mfvwPabM9r1qZ+1W53MiGhIQEiYmJqdGfhISEKvcZGxsrgOzcuVNERO6//34jb79z587yyiuviIglf37QoEGSmZkpIiIvv/yy/OMf/zCWe+edd0TEMgaAn5+fMT5A27Ztjf2UzMIPDw8XEZHs7Gzp06ePREdHi9lslmbNmhm5+NXJzK8o2z4tT+SRZ+bK8/+YJ7+mF8qNg4fIgw8+KCIi27ZtqzCH/p577pEdO3aIiMixY8ekV69e5c7Vli1b5I477jBe11Rmvq+vr5w4cUJELo6PUDYjPy0tzRhdbOPGjUY5VqxYIV5eXpKamirZ2dni6ekp8fHxIlJ6tK2DBw/KmDFjJDc3V0REwsLCjJx+kcqz8muDjmzQbMrULJ2wM+9a4pc71n67/eJwr9rWqVMnBg4cCMDkyZNZtmyZkWFfHDH8008/cfDgQWO53NxcBgwwBqdj7NixgGUMgMzMTGN8gCZNmpCaWmpICzZs2EBsbCybN282thUTE4ODgwP9+/c3MvGrk5m/ZcuWCrPtm9kVkm7fEg+Vhb3CGOVs8ODBpKenlyvTpk2bOHjwoPE6PT2dzMxMY/QusHwjK87CL16nJjLzBw4cyLRp0wgNDTXGMSgrLS2NqVOncuTIEZRSRqoowLBhw4wgOh8fH44dO0anTp1Krb9582b27NlDcHAwYBlLuG3btsb7NZGVXxt0xa/VOPPur4hwn8bc1glEJFyLqWXj6LRVNq+95OvmzZsDlm/YI0aMqPSWQMkxAMqOD1D2fr2I8OKLL3LzzTeXmr9161Zjf8XLXSozvyLp+ZAl9ngUniepsAkFUvUxAhQWFvLTTz/h7Oxc6XbL5vNfrsoy85cvX87PP//MN998Q9++fdmzZ0+5defOnYvJZGLdunXExcUxdOhQ472y4xdUNv7B1KlTWbRoUYVlq4ms/Nqg7/FrNcp8JpPQ7FuIzHmDcL9rG1UTzvj4eHbt2gXAJ598wo033lhumZCQEH744QejFc758+c5fPjwFe1v1KhRLF++3LhqPXToULks/OLlLpWZXzbbPu5MMjE54O6gcHN0oOuFaLILhdWfWuKGd+7ciZubW7mo5pEjR5YaIKXsgC8A3t7epVoh1VRmfnR0NNdffz3h4eF4eHhw/Pjxchn5JcdL+PDDD6vcXjFHR0fj3A0bNoy1a9caYx4kJycb4wlIDWXl1wZd8WtXrWRLnqhDPxEZPRF6T2VxPI2qCWfPnj15++238fb2JiUlhbCwsHLLeHh48OGHH3LPPffg7+/PgAED+Ouvv65of2Wz8B999NEKr1Krk5lfNtv+6f+bRVdnaGIHGTQBh2Y0LczBvkkT+vTpw8OPPsriiPfL7WvZsmXs3r0bf39/fHx8jIfSJfXq1Yu0tDSjQq6pzPzZs2cbD4NvuOEGAgICymXkz5kzh2effZY+ffpUu8XTww8/jL+/P5MmTcLHx4cXXniBkSNH4u/vz4gRI4wH/zWWlV8LlFQyDF1d0q9fP6mo/bBWNxiDrHTPxLSmE+YufyPU/YUay+X5888/8fb2vvoNWVFcXBxjxoxh//79ti5KjUvPh5jsQrrmHMG1WUvSm7SzdPRy5oojHZYsWYKLiwsPPvhgzRbWhmbOnMnYsWNrJDb5SlT0/0QptUdE+pVdVl/xa1fNCGP7UzGv1f8j1G2+DmNrQCzhbXbENOnGyQtCTLZcVaUPEBYWVuqeekNQU1n5tcGqFb9SaolSaodSammZ+c8ppRKUUi9Yc/9a7THZxRJ2ammjDWPr0qVLg7zaL+bqAB4OhSQ6XoMHGVcd3ubs7MyUKVNqpnB1RE1l5dcGq1X8SqkgoIWIDAKclFLBJd7+NzDJWvvWap9554dEtHmIue3TiUhoHA9zG5P0fEgqdKK9pJIkTUm/UENdeDWbsOYVfwiwsWh6E2A0VhaR00Ddf7igVYv50K+EOj1OZNPvCO/p2qha8jQGRnibM3Ro3pyuuXHEXLDT4W31mDUr/pZAetF0WtHralNKPayU2q2U2p2UlFTjhdOujtGSpzCfqINbiTw9A64LbXQteRqDkuFtp/Idwbk1XS9Ek5WbDVg+GE5dsHEhtctizYo/DXAtmnYFUqtYthwReU9E+olIv5K9/LS6wcjk2f01c2JmQUAYoYebGJk8jT2MrSEpGd7WzB5ixB3sHLkmK5r0vEKbDs+YkJDA6tWrbbPzesyaFf8uoPgR93DgJyvuS6tlJneI9DxJaPog5nl/RGj6oEbfkmfZsmV4e3sbYV41Yf369aUiEObNm8emTZtqbPvFSkYuV9UByRie0akLJ1UrYnIqb+FTMlCuOlJTU0ulW27dutWIZqhs+VmzZl1RS5rly5ezatWqy16vobBaTwMR2auUylFK7QD2AfFKqedE5EWl1APAY0ArpZS7iPzNWuXQrKSwANMP9xBmN4aFbecw99rGXekDvPPOO2zatImOHTvW2DbXr1/PmDFj8PHxASA8PLzGtn2lXB3Aw8mORK6lfW4Crk1cgRaXXO9Siiv+xx57rFrLt2zZslTGT3Xl5+fz6KOPXvZ6DYlVm3OKyEwRGSQiT4jIKRF5sWj++yLSV0S8dKVff5TsocvepZjT7FnW7kmGtaTRt+R59NFHiYmJ4ZZbbmHJkiUkJyczbtw4/P39CQkJ4ffffwdgwYIFTJ8+naFDh9K1a1eWLVtmbKNkLPKUKVP48ccf+frrr5k9ezaBgYFER0czbdo0I2xt8+bN9OnTBz8/P6ZPn86FC5Yb7V26dGH+/PkEBQXh5+dXYc/g7Oxs7r77bry9vRk/fjzZ2dnGeyVvrb7++uv4+vri6+vLG2+8AcAfR+MY5O/N4hkPMvjGEdw06hayz1uipbf9vIfe/gEEBASUimGoTjT0M888Q3R0NIGBgcyePRuwRDZPmDDBiLwu7nC6Z88ehgwZQt++fRk1apTRe3bfvn2EhITg7+/P+PHjjYFRhg4dypNPPkm/fv1YunQpCxYs4LXXXjO2FRBgKfPs2bONbyllB5gZM2YMW7duBag0WrveqCiys6796FjmumFLskibnSJb4uJky79GituWTHHdXihbkku8d5URzBW53Fhm2TJT5LMhNfuzZaZcSufOnY1Y48cff1wWLFggIiKbN2+WgIAAERGZP3++DBgwQHJyciQpKUlatWolubm5FcYii4hMnTpV1qxZY+yj+HV2drZ07NhRDh06JCIiU6ZMkSVLlhjlWLZsmYiIvP322/LAAw+UK+s///lPuf/++0VE5LfffhN7e3uJiooqtczu3bvF19dXMjMzJSMjQ3x8fGT7L3vlm/2xYm9vLzuifpW96QUyYtwEee+9tyQtT6RHbz/5dss2ERF56qmnSkVIXyoaumTktIiI2WwWV1dXOX78uBQUFEhISIjs2LFDcnNzZcCAAXLmzBkREfnss8+MY/Hz85OtW7eKiMjcuXONuOchQ4ZIWFiYse358+cbsdl+fn6ybVv5MpeNdB49erSYzeYqo7Vt6XJimXXPXa3aTO4QeV02oUddePGavyP2TVnvqzC565Y8Ze3cudPooHTTTTdx7tw50tMtjdxGjx5NkyZNaNOmDW3btuX06dOVxiJX5tChQ3h5eXHdddcBMHXqVLZv3268XxxL3LdvX+Li4sqtv337diZPngyAv78//v7+FR7D+PHjad68OS1atOCOO+5g244ddGoCXl5e3NgvkG5N7fAO7MuBY6f4LTGZ7PRUbjENBijVQWvDhg2sWrWKwMBArr/+es6dO8eRI0cueR779+9Px44dsbOzIzAwkLi4OA4dOsT+/fsZMWIEgYGBvPDCC5w4cYK0tDRSU1ONMLey56Q4Gruk1NRUUlNTGTy4fJkrUzJaOzAwkJUrVxpBbfVF3U8T0uoOEUw/TSMs05eF7ecyt1Pp+/rFHwA2Z3rD1iWoUnXif2tqHzW9fVcHcHG4uH1XB3B1sud0jj2t88+gKumeI1cYDV3RuRIRevfubSShFktLS6tyWyWjqqvDwcGBwsJC43VxlLRcIlq7PtBX/FqVSt3Xj1qMOeEMy9rP1vf1L2HQoEF8/PHHgKV1Sps2bXB1da10+bKxyMnJyQDlYoWL9ezZk7i4OCPeePXq1ZeMLS5p8ODBfPLJJwDs37/feAZR9hjWr19PVlYW58+fZ926dQwaNKjUMun5cL5Q4aSEC6274+rSnJ3btwLwweqPyS/6HKhONHRlx1rRsSclJRkVf15eHgcOHMDNzQ13d3d27NhR7XPSsmVLWrZsyc6dOwGM3xlYnpXs27ePwsJCjh8/bgw1WZPR2rair/i1KhW3149s9Qv8+j3jr/svYt+E5zpb3g89SKNvxlmR4oe4/v7+NGvWjJUrV1a5fMlYZHt7e/r06cOHH37I3XffzUMPPcSyZctKDbLu7OzMihUrmDhxIvn5+QQHB19WS5WwsDDuv/9+vL298fb2pm/fvuWWCQoKYtq0afTv3x+wxDv36dPHuHVU3KO3lSOcd1Cg7Hn2nZWEPXY/yt6JINNI7IrGaXnwwQeJi4sjKCgIEcHDw4P169eX2l/r1q0ZOHAgvr6+3HLLLYwePbrCsjs5ObF27VpmzJhBWloa+fn5PPnkk/Tu3ZuVK1fy6KOPkpWVRdeuXVmxYsUlz8WKFSuYPn06SilGjhxpzB84cCBeXl74+Pjg7e1NUFAQUDpau/iB+gsvvGDcdqsPdCyzdknmI/sJPXYNAXlHiWp+Pev9lFHRm1Ms9/Wt2VmrPsQyN0anLlg6bhW34U/Ph+jsQpoVZJJt15yuzexwdVBVb6SOqc/x2jqWWas5p3Zj+m4gYRmfs7lpCDM7qXL39XUP3capZI9esEy3dbIjw94Vj7zTuGbFQz24sGyMdMWvlWPc10/6Hb4Yhdn9Vpa1CdP39bUqpedDUh642AtnnK4h/UIOZJ4AkXqT59PQ47WL6Xv8WjnBLhD6Ry6R0c+DyzDGe36EYKfv62uVKpngCYrobIhu0oNuOUdAThOj2tG1af267dOQ6YpfK8eUuJrIQ6sI7fo5AW7OSLY9630vVvTF7fV1xa8VK5ngCdCtqaXyT3TuTLbY07XgBK527dFVTt2gb/VoQNHtneQC2Pk8/O8+TG4F3NKuOZszmzGzY/n2+vq+vlZSxff7FRk401QJ5GVB8iHIt7SFry+3fhoqXfFrAASrk4T+mo75wA/g+wCvX7+Bj842YUo7fV9fu3zF9/vbO0EWTkQ79yAdJzh3kPTzqcTkiM2inDVd8WtSCL+/h2ldTyLj7yO01/+4z/PfPBXrwGvdYJU3ekQtKzlw4ABff/21rYtR40qN2NUEdn7ybzJS04lu0p2TTToTU9CcrvkncVV6+EZb0RV/I7b49xjMax6FjY/ANf0x3fkOt7R1ZvVpmNwOZnWyLNfQcnhK9UYuYk6xzK8t8fHxvPjiiwwdOrTC90tm0X/99de8/PLLlW6rbI59QkICEyZMqNHyXo6y9/tvHW7i7acfxy43h0T71ripfFwvnIaz+yEzgfS8wnK3fY4fP47JZMLHx4fevXuzdOlS473k5GRGjBhBjx49GDFihJHAKSLMmDGD7t274+/vz969e411Vq5cSY8ePejRo8clO9NZS4sWVx9dXWMqSm6raz86nbPmvHJMZMvRgyLrbpcty4dKm41n5Z9Ru+WVuAL5Z7yIMotMOWi9pM0rcdnpnJdQNknUmsmiZeXn51drObPZLKNHj67WsmVTLeuKvLw8YzotT+TXDJGYLJGodJHE7DxJTD0riUnx8mtanqRlJosU5EtankhijkhCQoLs2bNHRETS09OlR48ecuDAARERmT17tixatEhERBYtWiRz5swREZFvvvlGbr75ZiksLJRdu3ZJ//79RcSSdOrl5SXnzp2T5ORk8fLykuTk2v/jbt68uVW3r9M5tfIKcuHIlwRHPUFoTBvMKQWYAobxbHcXnsrsy6ZUO56KplHc3in+BhN6EObF1kzz1Li4OCMz3tvbmwkTJpCVlQVY2oY//fTTBAUFsWbNmkqz3L/77jt69epFUFAQX375pbHtkrnwp0+fZvz48UZ+/I8//lgux77kyFc5OTncf//9+Pn50adPH8xms7HNO+64g5tvvpkePXowZ86cCo8rPDyc4OBgfH19efjhh408/KFDhzJz5kwCAwPx9fU1cmwWLFjAlClTGDhwIFOmTCEpKYmxd9xJSP9+TB8aTMLeH+jYBL4y/4Bp0DAG33QbU27sg0r8nfTkWGKyC2im8mjfvr0RkeDi4oK3tzcnT54E4KuvvmLq1KmAJYGzOPrhq6++4r777kMpRUhICKmpqSQmJvL9998zYsQIWrVqhbu7OyNGjOC7774rd6yVZfxXdqyVjbmQmZlpnHN/f3+++OILYx/PPfccAQEBhISEcPr0aQDWrFmDr68vAQEBRkqo1VX0aVDXfvQV/+V75ZjIlnP5Ise3ySubvpAtH4yVLcuHyitrFsmWn78Qt+2FMuxXy5XulIMimC3/lrQl2bIdW6vpK/5ic2Msxz035pKLXlJsbKwAsnPnThERuf/++428986dO8srr7wiIlJplntxvv7hw4elsLBQJk6caFzxl8yFDw0NNXL38/PzJTU1tdwVf8nXr732mpFV/+eff0qnTp0kOztbVqxYIV5eXpKamirZ2dni6ekp8fHx5Y6reFwAEZHJkyfL119/LSKWfPsHH3xQRES2bdtm7G/+/PkSFBQkWVlZIiJyzz33yLpN2yUtz1Kunj17iojITbeMkX9v2CkxWSI/JGbI/pRs2ZuWL2ln/pLEs8clLfmYpGWmSGJ2vsTGxkrHTp3k8Jk0ERFxc3MzylRYWGi8Hj16tOzYscN476abbpKoqCh59dVXjbEARETCw8ON302xqjL+KzvWysZcmDNnjjEOgIgY3y4A4/zNnj3bKJOvr6+cOHFCRERSUlLK/Q6q63Ku+HWj2lqyON7SMcrkfnEa4NXjMLvT1U9HZQBSiEPWafLT4wk++y2hTZ7g2YR1RDfz5aUunyB2jqz3swdlT+4fsDkVprSD/yXD3M4XW+8UX/nWmZhlKzCnWI63+LhNLa/+WDt16sTAgQMBmDx5MsuWLeOpp54CLmbBl8xyB8jNzWXAgAH89ddfeHl50aNHD2P99957r9w+tmzZYowVa29vj5ubm3GPuyI7d+7kiSeeAKBXr1507tzZSJIcNmwYbm5uAPj4+HDs2DE6depUan2z2czixYvJysoiOTmZ3r17c9tttwFwzz33AJakz/T0dFJTUwEYO3YsTZs2BWDTpk0kJCTwxkLL9vLz80lIzaT39QN567lZHJo4iXHj7iC7eUfL7Qc3L5rlZhBd2BIKhCbn4rl73O08ueh1nJo5k55XAMCRLGjnBKAQpUjPh3yB4zmWh8tgedaQmQ85hXAu7+L8c3nQponldZZlc/x18BB/FGX8A+TlF9DmmvbGeajoWHfu3GlczZccc2HTpk2lhoR0d7f8YTk5ORnPbfr27cvGjRsBSxjctGnTCA0NNcZRsDarVvxKqSVAP2CviMwsMd8XWA4oIExEymfCXgVrVLLFDzYdFEbU7OVsNzrb8vpZT8v0S8dAgPldYPz+y50WRAqZ7x7LRvUepgAADKdJREFU+N86IVLI+tS/82tWE5669iVeO74YU/Zanu3qzVOerzPZIx9JdkQB5jRYegKcFExoBx+dttzemdXJUvk1hl655pTSx1lTx62UqvR1cRa8VJLlvm/fvivf8RW61LgAOTk5PPbYY+zevZtOnTqxYMECI5MeKj/ekrn3hYWF/O9//zM+CIpb/Lzw3DMcu2M0n//nW+4cNpDVX39P6269iL7gRFun1pAr5OXlMfP+hxgWOpVJN19Ps6xYogu60dKjHVnHDxB9rTdJpxLx8PAgOlto3q4DqYnHiS4aRTIp4QS5rTtg36YDZw9vNeafSzhBuxuGEp0N3ZpCViGcuiD09OnNf7ftAgWncsHdwVLeAoGTF5TxwZEvkFGgKBCIzYY2JeYXLx+XDe2K5p/OtXxIOTo6cjRb0c4JLmBPck4+6fmw+K3lbPjhZ6I2fkNQ376s2bGHXu1bA5YPpmtKfEhdc/FXdlWsdo9fKRUEtBCRQYCTUiq4xNsLgXuA0KLpGlUcJWxOsUyP3w/j9sNw9yufDnaxVPpPRVv+vdzt3t3WUuk/FS1k5xcgUoiikNSsdEQKUFJAasoJpDAPVZhLasI+pCAHVZBN6pFvkLxMVF4Gqb8sQy6koXIzSP3jEyQ/G1WYi9nBl0XXPsdrLfawyOtV5o2IZ1Gzu5ncTrE6yZGZHWFGR1h4zPIHus4XfJtbKv1F8Rev9BtS653KRGWUruRr6rjj4+ONjPhPPvmEG2+8sdwylWW59+rVi7i4OKKjowEqHeRj2LBhREREAJZxbNPS0qrMsS85LsDhw4eJj4+nZ8+e1Tqe4kq+TZs2ZGZmloqFBvj8888By7cKNzc349tDSSNHjuTNN980Xu+K2kNXZ0g6Fk1nbz8mzHiavv2COX7kL7o1tVzcJOaChyMsfuJBOvT0Zcb/+z9OOXUko3kXUHYMvmUMn33yKUgh//14JSNGDQcpYMgto4lc+W+kMJ/ff/mBFi4t8LimHQOGjWDbpu9JTzlHekoy2zZ9z4BhIwEhI7+AUxeEG3r34ExSEj/s2sWJC9Ba5XHq0AGisyG7EDZ++TnR2bB2i+VYzzq50TtkEN+v+ZjobPh8g2XMhSRHV/oMHUHkv962xFZkQ0FGCtHZUIhlEJvobMuHgaOdZXrrgWhuHHA9dz8dTotWHqSfOm6s28z+4odlTfZ7sOYVfwiwsWh6EzAAiCp67S4ixwGUUi1resfGw7t9GYSlrkZcJ6OA1N0fVWN6NeI6pdy0eesqIlzv5bXkCBYVhBGW9jHidh8KReovK5CW0yzTP7+PuE+3rL/rPaT1w5b1v3/7/7d397FV3WUAx79P7217Cy0UCk0Ld8wqSsWWUF5kTEBkuMlQRIaLI4ITpLVUibo3xnbrsu5FM7LROmhEIRiUTDOdLMYJMww33IyabDFO+McwsbasHe9O+nLbxz/OufQWClzKvffQnueT3OS5L22fkzTP+d3n/M7vR2PBWr58+iV2sYpI86MA1HXX9sY9cfG4WCzUjYsQOdUImbnUFawnIvsheyR1wVoixWcgM4+6o2uI3AjfKZnJqSNOgY9v49Q3OV+vbsnvLXCxu28rcnuXYBjK7Z2Y/u46TsZxT5o0iS1btrB69WomT55MdXX1RZ+53Fru27ZtY/HixQwbNoy5c+f2W8zr6+uprKxk+/btBAIBGhsbmT17dp917Gtqas5/ft26dVRXV1NeXk4wGGTnzp19RvqXk5+fz9q1aykrK6OoqIiZM2f2eT8UClFRUUFXVxc7duzo93c0NDRQU1PDlClTiEajzJs3j9tmT+ehzZt5ef8rBDMyKC/7GCuWLKID5380LwD7Xv0jv9y9i9Kycj4zayoBoKr2CZZ/7nbu37CRr664k+d+9lHCN9zAIzufp1DaWXbrPF7bt5dlUycSyslhY+MuCjuPUZgLd99fy1fmO6sT3/3Ad/lwbgd0ttDCOIo7mxnf1cyPt29lw0MPcO70CTq7e6iqWs2nSyYBSlagm7s+MZVoNMrmhqdAo1RtiPC9mlW8MLucUM4wnn72adBuvnbvg2y+p5IvziojIxDgvvu+xYylqwDofr8VAgWgPWi0HbSbzQ/fS/M/DxMlwMc/uYDSDxXTpk4P6uzZ47QFx/LBHOlzZ/S1Stl6/CKyEafF8zsRWQjcrKqPuu+9qqrzLowv+PlKoBJgwoQJ0weyp2XtwYPURecQ6XBGKnXZyxOLO52+XV3WHRfEv+LRruepzVxOXdYyIl0vAEJd5lIi0T1OHFxCpPs3IEJdxmIiuteJuZWVvMlLfITq4JvUR2cgCOtDf6e+vRwRWD/iKPVnPoCIsH7sGerbRjpxWM4X7vVhLhk3NjvfKp48CotG97ZxKnJ7W0W/diZ7DKqWzmBYj38wr+M+EPPnz2fTpk3MmHHRUu8D0neRN863ZWKtmKYOKAjCKbd9UpgFrZ2JxE59KwxGaY0GAaUwo4PWHufkVyjnaNMQRfI+x3Q4Y/kfbQxjJO0cZzjFeprlixdT/dhmbqmYCECL5FOspwClRUZRrCdBoSVjFMU9zvWWy8cn3Hh0YnFmN+NDVx7uX816/Kkc8Z8GYnvNjQBOxb0Xf7bpoR+qug3YBs5GLFf7x185CY0yxx3tLkeASDjR+I5+48bmZeRPWEbjUYiMg/qmL8R97vNx8Wfj4tsQYOUY+Om7FW4hnkODW4jzx09H3nHjMROR/7px7ijkPTcOOsX9cvGn8p3nsSmZUe1t4ywb47R2oHfTFFtozVxP4m/6OtbhFHyAk1HnEc6GjrhKEZBEY+eFQCATos7zQGYOuDeM5YWGE+iBpo4RhLOhKCuPQCc0dQynIAit0ZF0EWRMJrS6zYniLGjtjI9HgcTFXCkefVVxW1eAvChJHfGnsvC/AVQBvwAWAjvj3jshImGcon8m2X84/uIdQENT4gU00cJakXt1vzcno/9C/NS/kxPH2jexot9fGwf8MWPHC35Zxz3mwIEDSf198Rct4+NLnRDe7UxOHJvVE84G3Iuzxzr7nmh+9NsDFGdDi3uySPykc+1xXsB5xL4NJav4p3TrRRGpB6YBbwGPA2tU9XERmQI0uh+rUdXLTmm42q0Xr8dZPbGRdjq2KhxqDh06RGlp6UWzSIxJtvjtJGMx9M7MSXec6KweVeXw4cMJt3psz11z3Tty5Ah5eXkUFBRY8TfmAqrK8ePHOXv2LCUlJX3e86LHb0xShMNhmpqaaGtr8zoVY65LoVCIcDic8Oet8JvrXmZm5kUjGWPMwNkibcYY4zNW+I0xxmes8BtjjM8Milk9ItIGXP2tu94aA7zndRJpZsfsD3bMg8eNqjr2whcHReEfjETkr/1NoxrK7Jj9wY558LNWjzHG+IwVfmOM8Rkr/Klz8fZJQ58dsz/YMQ9y1uM3xhifsRG/Mcb4jBV+Y4zxGSv8aSAie0TkMa/zSAcRqRSRP7mPFV7nk0oi8oyIvOYuPz7kicgsEXldRA6KyDNe55NOIvJtETnodR7JYoU/xdy9B3K8ziON9qnqTcBc4B6vk0kVEZkG5KrqXCBLRGZe6WeGgH8BC1R1DlAoIuVeJ5QOIpINTPU6j2Sywp9664GtXieRLqr6jhtG3cdQdRPwshv/HpjtYS5poarHVLXdfdoFdHuZTxqtAX7idRLJZIU/hUSkFGij737DfvF1YI/XSaRQPr3bhp52n/uC+y12rKr+w+tcUk1EMoH5qrrf61ySydbjTwIRKQKeu+DlYziFoRYoTXtSKXapY1bVL4nILOB2YGn6M0ub08AINx6BT07uIjIaeBa40+tc0mQlsNvrJJLN5vGnkIjsxdlrfTRQAKxW1T94m1Vqich44OfAElU94XU+qeL2+KtUtUpEtgI7VfXPXueVSiISBF4EHhnqxxojIt/H6e8rMAuoVdUfeJvVtbPCnwYiMh9YqKoPe51LqonID4EFwH/clxap6jkPU0oZdzbPNOAtVf2m1/mkmojcBTQAb7svPaiqb3iYUlqJyEH3wvagZ4XfGGN8xi7uGmOMz1jhN8YYn7HCb4wxPmOF3xhjfMYKvzHG+IwVfmOM8Rkr/MYY4zNW+I0ZABGZKSJ/E5GQiAwXkbdFpMzrvIxJhN3AZcwAuXsshHCW3W5S1Sc9TsmYhFjhN2aARCQL+AvQDtysqn5ZptgMctbqMWbgCoBcIA9n5G/MoGAjfmMGSERexFmaugQoVtVveJySMQmx9fiNGQARWQV0qepuEQkAr4vIgqG2YYcZmmzEb4wxPmM9fmOM8Rkr/MYY4zNW+I0xxmes8BtjjM9Y4TfGGJ+xwm+MMT5jhd8YY3zm/ylq9LDes5wcAAAAAElFTkSuQmCC\n","text/plain":["<Figure size 432x288 with 1 Axes>"]},"metadata":{"tags":[],"needs_background":"light"}}]}]}