{ "cells": [ { "cell_type": "code", "execution_count": 25, "id": "b4f91fdb-98df-49a0-a20c-c99fe4fd3be5", "metadata": {}, "outputs": [], "source": [ "load(\"www_utc_fr.sod\");" ] }, { "cell_type": "code", "execution_count": 26, "id": "d84e5294-c29c-4e61-8866-b7b9918e26cc", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiQAAAFQCAYAAACPqGU6AABQkElEQVR4Xu3dCZgU1dk2YL98cYsrUaImrji4EBYjKqCIMqOCiJKIG0g0v6BGcYUgfiYKKhpBUDQBE4G4AuJulCUSccFEEB1BFkUUFxQVAVGQVXz/efpwpqrOVHVXdVefqul+7uuqa2ZO1cx015yqeucs79lKiIiIiBK2lVlAREREZBsDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEsrb5s2bpbq6WpYuXWruIiIiioQBCeUFwUhlZaX07NlT2rdvL6NGjTIPISIiCo0BCeVl6tSp0r1798znGzZskCZNmsjGjRuNo4iIiMJhQEJ5GTJkiAwfPrz26w4dOsiiRYtcRxAREYXHgITyMnDgQBk9enTt1127dpUZM2a4jlC22morny2onBs3bmnaiGxijaO8jBgxQoYOHVr7dVVVlSxcuNB1RDDc59QNT38ucthhIv/zP2pzl7s3c9+PfyzyzjvmTyeiODAgIdtY4ygvc+fOlU6dOmU+X7VqlVRUVGQGuobhF5Bsu61I8+ZO0IGPKPvZz0T23ltku+1E9tlHpFEjkQMPFNlvv7oBy+67O/v1VvOyasvcn2PD540b+3/PL38pcscd5isnKh8MSMg21jjKW79+/TIzbFq0aCETJkwwdwfyC0iw/ehHdYOMbBsClgMOUAFL1O8NuyEwQjD01VfmuyAqbQxIyDbWOCrIunXrZNOmTWZxVuphXzcgcQcB5qb3HX6483OOOUbko4+cr+P0v//rfS368z32MI8kKk0MSMg21jiyTgckP/mJ86A3A5LzzxdZuVJtDz8s0qIFxq2o7+/VS+S88yTz/TXxUFHo1/TDDyJt29YNjMzgCS00eI1EpYIBCdnGGkfW6YDE/VA3H/I77eQcj94gBAZa796qpQTHjh+vyt57T+Sbb5xjCqVfByxZIvLllyJoCHK/znbtVCCCTQcm06aJfPKJ+h581Bu+9ivL9vWKFd7XRGQTAxKyjTWOrDMDEr9NP+D1htaHZcucn4GuGowhueUW9fWpp6rAJS477KBex7XXOoNqDzqo7uss9obXcd115qsjKj4GJGQbaxxZpx623oAEM1vAbC3R25lninz2mffn/Pa3SNDmLYvDm2+KHHqoSOfOIsOGicyaJdKqlUjHjiowQlfR2LEi331nfme8MGBXv3+cl4suMo8gKh4GJGQbaxxZ5xeQzJmj9qFbxG/GzKWXen8GzJwp8t//mqXFg8S0775rlhYPuq0QiLRs6Q3U0GpCVGwMSMg21jiyzi8gwRgNWLWqbjCCDQ/lefNETjtNBPnXMKbkpz8VOfFE78+O0447iuy6KxYSVF8/8IAa4GoLWkgw2+ess0S23lp1UW2zjXNOjjzS/A6i+DAgIdtY48g6v4AEgQhgIKcZjGDDg/nf/xbp1k3k+utFPv9c5M47RYq5nh/GjaBlwj2gFvB7zbJiQKI3tBZNmSIyYIDImDGqm6hZM+e8oLWkWDONqLwxICHbWOPIOr+A5Ljj1L733687oFUfgyyqPXuqQZ7meJJieOUVkUGD1OevvSYyf75IZaVKWY9AYfly7/FxQqsMuogwe8jdTYSgDIssY4zLwQc75wbjWp5+2jmOqFAMSMg21jiyzi8g0VNsAYFH69YiL76optG6gxK0WqDb5vHHneOL4fvvnUBg3Dj1O2+7TeTGG0XWrnWm+hYrMZuGn+9awzDjr39VgZm2227OOSpmkETlhQEJ2cYaR9blCkj8+M2+KaYvvlAzf9Aa0rWryP33q+6hK65Q+xGwYDwHXtfs2d7vtQWB0mOPqc8xvgZBkp6tRFQoBiRkG2scWWcGJNtvbx7hz2ZAAhjUmitQ0i0lSQQlmJ7s/r26i4soDgxIyDbWOLLODEjcD1E9kBStDybbAYnuCskFA25xnJ66bMvgwSpPioaZSBhLQhQHBiRkG2scWWcGJDoocecf8QsE0HWChz82LKyHcR2jRplHxaNJE5Fzz1VjWODee0WeeUa9RnTj/N//eQfW6hTyScLrweJ/6E4iKhQDErKNNY6sQ04N3Oz0LBEMYsVHvcIuukrQHZENxnM0bOgdW1JdbR6VvyefFHnnHefrCy8U6d9fZYe96iqR448X+c9/nP2AfCHFnIacy/r1Ig0aiKxebe4hio4BCdnGGkfWzZihbnbIyoqHPB78Xbqo2SRIPPbBB+Z3+EPuEuTk2HNPkUMOEdmwwTzCLptZXImKjQEJ2cYaR9ZhRghudkjyNXGiuTc/337rZFRNClomEGQliUnSKC4MSMg21jiyTnWxbFXb3eIeL6K7baKOx0Badaz4mySsPYPXkaSdd1ZjbYgKxYCEbGONI+t0QKLHfrjve35BShhomSj26ru5YAouXjuCqqQgGEEXFlGhGJCQbaxxZF0xApK4PfxwfrNVsPpwkkEJZtpgkDBRoRiQkG2scWRdfQhIDj88/zEpatBu8YOSV18VufVWb9nrr6tzN2mSt5woKgYkZBtrHFmX9oAEK/lidd1CJNVSgsG9mK20bJm5hygaBiRkG2scWZf2gAT5PPSaNW6YwRJ2SjIgKMH7sB2UDB0q8sILZilRNAxIyDbWOLIun4Dk6quLv7KuhhYSdLto772nPr7/vsjf/+6Uh4HkaXg/fqnwiwXn6rrrzFKiaBiQkG2scWRdPgEJWi2Scs01hSVde/RRu0HJrFnRAyciEwMSso01jqzLJyCp78aPV+/LRp6SIUO4yB4VjgEJ2cYaRxmLFy+WRYsWZbYvvviitnzz5s1SXV0tS5cudR0dXB5GfQ1IsE7Nxx+bpeHZbClZu9YsIYqGAQnZxhpH8sMPP0jDhg3l7LPPzmx33XVXphxBR2VlpfTs2VPat28vo7YsrRtUHlZ9DEimTxd5+mmRu+8290QzYULxgxLkT0HGVoyFIcoXAxKyjTWOZOXKlXI8lq81TJ06Vbp37575fMOGDdKkSRPZuHFjYHlYaQxIcj28MbD1ww/N0vzYaClB11ChU5epvDEgIdtY40gWLlwobdq0kXHjxsm0adMyLSAwZMgQGT58eO1xHTp0yHTpBJWHlbaA5JFHRMaONUuLq9gtJaefLnLBBWYpUXgMSMg21jiSVatWyR133CFjxoyRCy+8UDp27JgpHzhwoIxGlq0tunbtKjNmzAgs96MCD3NLV0DyzTdmiR1oKcH7tJ2nhCgMBiRkG2sc1dGsWbPMYNURI0bIUGTZ2qKqqirTmhJUHlbaApJLLxX5+mvna/Q+3X+/83UxJZXRlSgXBiRkG2scecZ/rFmzRho1aiTr16+XuXPnSqdOnTLlaEWpqKjIdOcElYeVtoDEhEGhzzxjlhbPa6+p9/vjH5t76tpy2omKjgEJ2cYaRzJx4kRp3bq1dOnSRVq2bCkTMMBhi379+mVm0rRo0SJUeRhpD0iSgKAkzpYSNFhNm2aWEoXHgIRsY42jjO+//16+dvdbuKxbt042bdpkFgeW5xJXQDJsmFlSv8XZffPppyJvvWWWEoXHgIRsY40j6+IKSEqxBSCuBfmQN+Xkk81SovAYkJBtrHFknV9AgrQmyJsRNiDBkJViLrb3yiv+eUd69BDp08csjZduKdl+e3OPF1pBsHifhtWIt+S0k/nzRX70I5EDD3T2E0XBgIRsY40j63RAggemOyjRG4KRQw4xv0vklltEvvpKfX7vvcVdr8UdGOEjftfy5XVf6x57mN8ZDwQlO+2kBtgGQUD24ovO1wjS0DKibb21eo0IXIiiYkBCtrHGkXWYKYKbHf6j32EHkWOOUf/N7767CgCuu878DuX999VD14Y33xTZZx+Rxx4TOfRQkQsvVNlcL7tMZNw4kYYN1Wt1zX5OJaS6X7LELCXKjQEJ2cYaR9Z9+61zs/vgA5HvvlPjQbCm37PP2gs6csFrg9WrMb1Zfa4XrVu2TGTMmNwp55OGlPcRsvoT1WJAQraxxpF1qqtmK5k7V3UpIH26u3tkxx3N7/A3dapZEp8zz1Sv5aWXVH4QfI6xI3pch36te+7p/b60OeIIkV69zFKi3BiQkG2scWRdgwbqZocWCDzYf/EL9XDXY0oaNza/w19Fhb207xiPgdeqB4nitV51lfeYNDrlFJG99zZLiXJjQEK2scaRdViJFjc7JO/CQ37XXdVHvW23nchf/uIcX10tMmeOGruxdKkqwwyYysrsgz4LgbEteJ3o7mjbVuTgg0V22UUFIu7XisG1xbJihVlS18svqzE4QVauFPn5z7Eqs7mHKDsGJGQbaxxZh/scbnb7768+190f7g1lWufOaowJ1u9bs0aVHX+8Ou5Pf3KOi5N+Hcig2revSJs2amyL+/UhV8iDD5rfGY/XX1fdLea5cMM52XdfkRNP9Ja/8473c/yMXFOIiUwMSMg21jiyTgck7dt7AxL3R/f4EDxM0TLhTgq7fr3IlVd6y+KkW0AwEwhTjU86SeTyy52ABBtykhRrAC4G0WK6LlqGMI7FzxtvqNeIcSxPPqlaQzDI9g9/cI5ZsECkWzf1enEOXYs0E2XFgIRsY40j63RAcsUVdQMR/dE9ewWzXFCGB6obHsRhujXyoXN4gA5O0CLifp36IW/bzJkiN97oLdOv0RxTc9ZZIkceKXLttc7r5qwbCoMBCdnGGkfW6YDk97+v+4DXX5uQS0MHBTYgBwley847O6/p44+d12a+5qTv3S+8oMa9PPqouUdBaw7gdad9ZhClAwMSso01jqzLJyABdF/YDEow/TjotYD5upOGdPGLF5ulXnomE1EuDEjINtY4sk4HJJdcEhyQYF0bPwhKcAxyg3TpYu6ND8at4HeYAQke6Lqbxnzd9QGmANeX10rJYkBCtrHGkXU6IEGisaCAJNvYDD1d2H2/HDIk3sX2ggIjd/CBsRk4Rm8PPZS7haKYkOn2uOOcmUhEhWBAQraxxpF1OiBp2jQ4IMmVrRWrA7tbL372M/W1e/XbQuggw+QOSEytWomMHWuW2oMZQWjZyZaXhCgsBiRkG2scWacDkosvDg5I/IIBNyQqM++XehZMHEGJfh2mbAFJGvTrJ5n8LkSFYkBCtrHGkXU6ILn00vwDkkMO8Q8MdFCCJGqFSDogwSJ+N91kltaFVhE3dNtg1eRi5Ueh8sGAhGxjjSPrdECSbVCrXzDg5tdCoumg5L//NfeEl3RA8vnnKodINnh/w4Z5y5AoDoNxkSjND3KQnHeeWUpUFwMSso01jqzTAUncXTZueoVepH7PR9IByfvv101+Fha6rJYtM0sd335rlhDVxYCEbGONI+t0QJItD4lfMOB26KG5AwO0lODnTJ9u7snNLyBZvtx5jUhdX0yzZ4sMHmyWhnP77SLPPGOWEkXDgIRsY40j63RActpp+QckuVpINJ0CPiiDaRB3QDJhglpVGN0o+jVedZXIvHne74kTulbyTYv/5z+LDB9ulhJFw4CEbGONI+t0QHLCCfkHJEGDWv1ss030oMQdkIwfrwISXY6tf3+Rr792jk+TBx5QiwESFYIBCdnGGkfW6YCkqir/gCRMl42bDkqee87c488dkJjlUX5vEp54wjv1d+lSNSaFKAoGJGQbaxxZpwOSxo2DA5JsmVqhSZPogQGShuF70AWTi19A8tVX9SMggccecz6fNUtk6lTna6IwGJCQbaxxZeS7776T5RiZ6WPz5s1SXV1d89/00oLKw9ABSSGzbKJ02bjplpJcQYkZkOAB/8gj9SMgQffSrrvWzVFCFAUDErKNNa5MzJ8/X5o3by59+/Y1d2WCi8rKSunZs6e0b99eRo0alVd5WDogKWSWTdhBrX7MoAQtH+Y0WTMgwVRZJCurDwEJbLutyO9+Z5YShceAhGxjjSsD39f8y9ytWzcZOXKkb0AydepU6Y7FYWps2LBBmjRpIhs3boxcHlbSAQnooAR5StCd8a9/efebAYm7vJDfawum/mLBQaJ8MSAh21jjyshTTz3lG5AMqXlyDXfNE+3QoYMsWrQocnlYOiApJFNrvl02bpgSjN/jt/aN7YBk7lyzhChZDEjINta4MhIUkAwcOFBGjx5d+3XXrl1lxowZkcv94KZWd1PlSbaQaEFp5vXrMOnXGLfrrvN+vWCBmr5bDMzUSmEwICHbWOPKSFBAMmLECBk6dGjt11VVVbJw4cLI5WHpgKQYi+vlQwcl7jTzuQKSCG83L2vWiHz4oVkaj3POMUuI6mJAQraxxpWRoIBk7ty50qlTp8znq1atkoqKiszA1ajlYemApJAum7haSDSzpSRXQILU7mn2xhsi48aZpUThMSAh21jjykhQQAL9+vXLzJhp0aKFTHDNiY1aHoYOSNLQZePmXpAvV0Dip3dvsyQ56Jb59FOzlCg8BiRkG2sc1Vq3bp1swvr1BZbnogOSQlpIomZqDUsvyJdPQFKsLpZ8Bb2eX/xCrctDlA0DErKNNY6s0wHJ1VcHByTIo6Hde6+zlozWrVtwYFAo3VISNSCx7U9/Ern8crNUQVI0zCJCdlqzN61Pn/wX7qPywYCEbGONI+saNFA3u88+Uw/9/fZzWkWwoTvGbfJkkR9+8JYVW9u2knl9aYbssVi3Jkjr1iq4Wr/eWz54sEiEWdpUphiQkG2scZQI3uySgwAvj4z/VGZ4jZJtrHFkneqicRbXa9iwbpeMGxbaw3EbNph7knHLLWrBujRA69KwYWYpUeEYkJBtrHFk3eLF6maHrgRkSV292jzCC8djGmtapCUwAo4FoWJhQEK2scZRInizI0o3XqNkG2scJYI3O6J04zVKtrHGUSJ4syNKN16jZBtrHCUi280OYzTMqaq5YFDs2rVmKfn57ru6uUmITNmuUaJiYI0j63Cfw81um23U51g/5qWXnFwjgwaJjB/v/R7YZRf9vWrmjduLL4rcc4+3rFDV1Sr5WNKQ5Oy3vxWZObPwoOu550R23VWdQ8x2+vOfzSOIFAYkZBtrHFnXvbu62Z18sshuu4ksXy4yfLjIxo1q/+23R0+E9vzzInfdZZaWDrRqPPCAyBdfmHuiw6wmTBfWCwoefrh5BBEDErKPNY6s0y0k7pTx2R60+pi33zb3UFhYcgjne8QIb/ncuSKNGol06OAtJ2JAQraxxpF1OjHaXnupz/fcM3uLCLpncNy0aSLHHy/Svr3IffeZR8UP41KWLTNL0wstTbqVyc9VV4m8+qq320evCcRnD5kYkJBtrHFknW4h0Q/Co482j/BnLsJXbIcdpro1wvrqK7UlBWsAYf2abI46yrtwoVZVVb+CLyo+BiRkG2scWadbSA44QAUWr7xiHuEPY0t23llt551n7o0fZvt88IFZGgzjWKZMMUvjcc45qtsFvvlG5IIL1IBb/E4NrUzZUvAD3hPG7mi9e3uDvG+/dfZReWNAQraxxpF1GLOAm92ECSL77KNWrc0GM0xatDBLvTADBT/Phk8/Ffn6a7O0uMzZNfgaW48eIv/+t3efn9deE3n0UfW5+2dh1d9DDhE54QRnzSB8RCvKiSc6x1H5YUBCtrHGkXVbb61udscco1pLmjfP/p+9nh6cLTcJWg1srWD7+OMqAEqDSy4R+d3vzNK6Vq3yPz8ICHXryIIFatYT/ia6e6xYLT6UfgxIyDbWOEoEb3bph4UPMVCWyhOvUbKNNY4SwZsdUbrxGiXbWOPIOiRBw83uD38Q6ddPld15Z2GJ0Sh+yE1SWWmWUrlgQEK2scaRdTvsoG52f/ubGkCJcQwYs6Dzjey9t/kd9mEmSsOGZqlkZgbpMRZI6Z4GSP/unjkTl1tuEbnxRrOUygUDErKNNY6sU4Mot5KnnlKf6wGU7o9JQ9CE12EuQufOhfLeeyKffKK2JKfLXnmlypdivlY/bHmisBiQkG2scWSdDkiwNos7ENGbuXBeUpo2zf5azNddTOjm0jORsK4NFhK87DKRkSOdY37xCzVL5tlnnTLTsceKVFSYpSKXX64WE3TD9OZHHvGWUflgQEK2scaRdToxGqbx7rSTyBVXqAc/Ep5h38CB5nckAw/kbDlSkIodrxvbT35i7o0X0uZrSJD28ssqK+uPf+yUf/aZyP33q/VpgiAxmu4aw/ciuAGk4jfXE3rjDSfYwvHorrKdf4WSw4CEbGONI+t0CwnSrN90U7iuhnKHAb8XXugtQ8K4fJ4Z11/vBBpduph7vZBSXo35yb7tuKNKrmYGNVR/MSAh21jjysh3Nf8OLw9ILLF48WJZtGhRZvvC9VTZXBMtVFdXy1Ijq1ZQeRi6hQT/leNhF/CSyLBihfdrrH6MZwZacqJAplbd3XTggeZef1gAcbvtRH76UxUcIdHaqFEiDz4o0qSJE5hk6+Ki+oUBCdnGGlcm5s+fL82bN5e+ffuau+SHH36Qhg0bytlnn53Z7rrrrkw5go7Kykrp2bOntG/fvuYBNCpreRS82cXrmmvU4FrAANtJk0Reesl7jB+dLt4MbDD4FV1CUQbBIluu2XKiZyRhw8BbbOgqQmr63XcX6dVLtdgguEIa+yQXJyQvXqNkG2tcGfj++++lW7duMnLkSN+AZOXKlXL88cebxTJ16lTp3r175vMNGzbU/CfcpOa/442B5VHwZhevmthQliwR+egjkf791ZgQ5HbJxT0wF60lq1c7+zBIFoFCFHp8CjaMD0IA4i5zb2bwku82eLD5KigOvEbJNta4MvLUU0/5BiQLFy6UNm3ayLhx42TatGmZFhAYMmSIDMf0ji06dOiQ6dIJKo+CN7viwKDTbLNsTI0bOwFCoTAL6OqrVQCCVhAdo9ZUK5kzxzkO+Vvmz1efYxbPsmXOPj8YiPvHP0ombw1mFWEw8amninTqJLLXXiro6dPH/C4qFK9Rso01rowEBSSrVq2SO+64Q8aMGSMXXnihdOzYMVM+cOBAGT16dO1xXbt2lRkzZgSWR8GbXXFgbAe6a8Jq1coJSKJ0z/hBQHLttc6sIz1YGYnbglYkPv98kbfeMkuzw3RxtAbBrbeq144p2hQvXqNkG2tcGQkKSEzNmjXLDFYdMWKEDB06tLa8qqoq05oSVO4HN7WgjYjSi9co2cYaV0aCAhL3+I81a9ZIo0aNZP369TJ37lzphHZxUa0oFRUVme6coPIoeLMjSjdeo2Qba1wZCQpIJk6cKK1bt5YuXbpIy5YtZcKECbX7+vXrl5lJ06JFi1DlYeFm99e/irz+urmHwkKj1DPPmKXKo4+KPPGE87VOgGY6+ODgMSS33Ra+G+cf/1CZXo84wvl52ZLKmZDsDWvnaJhNni3Bm9vzz6tZOxhrQvFhQEK2scZRBmbifB2QhnPdunU1D4xNZnFgeS5qhsVWmTVYMPjRHISJH9mjh/ocU0mnT/fuL7Z33nHStCfpv//NnjQO5wj5Qfw8/bSIjhMXL/YGHBgEqr/GTBg8d7CZwYeeCYNU8xhUGmTIEHUcfq57TIo7IPKDTL16LPTKler36OqEmeQXXeQcm82bbzrvAStIZztnFB4DErKNNY6sUw8PVfW+/FLkySe9+/Hf/NFHq8+xP9eDLW4YHpOGfBhIB4NzgUChc2dzr2oVwIP/l7/MPT0Xg101HWjkiiXRwqEf9Fj5OO4g7eOPvWvxFALBjX5f2FxDnChPDEjINtY4sk63kFBhWrRQD1937hAkR8uVEkY/uNGiYQZeaKnApj3+uPOQx/dtvbVqVdllFzWbBh/RyjV1qvM9ScK5wPvC69XrDOlpyHjt2BDIjR9vfieZeI2SbaxxZB0DknggN515Gu+4Q+SDD7xlpvfec4KMESO8+zBl2Jw27G550IGJTnaGhzsCgEsu8X5PUvA6dJdRtk1PTcbqyGj90duhh6pAr1kzkebNVbbagJ7MksdrlGxjjSPr1IMt3qqHcQNIuFVOMF4in9OIgEUHF0jd7obWFb8Wlg8/dL4HMIBVD6jFAn1+35MEtBBhMC5e56BBquyf/1RZaHWX0/vvO4GJO9DKtgUNCi5lcV+jRLmwxpF1yFKf62aHVOiff26WBsN/uvhPvZxgsb2TTzZL/WHGijs5GdLCY6E8cyArWjyw+UG3B4Kg+gDJ1jCuJA46eCk3ua5RorixxlEi4r7ZbbNNeT40wsLCeQsWOF9juvVllzlfa3rMRX22Zo2TyTUODEiI7GCNo0TEfbPTgyzLCbqpZs82S8NB99Z115mlKrDDVp/hnDz4oFmaPz1Net48c09pi/saJcqFNY4SEffNDt0J5fZfLHK03H67WVqYKC0k7qnEpUy3kMQZ5NQHcV+jRLmwxlEi4r7Z7bpr+AcpBYvSQoIuH4xjKXUIdrFNmWLuKW1xX6NEubDGkXW/+Y262aHbAAsL6zwamKkxdqz32LAwCwJJ1MoJ8oUcc0zuBGcmdPXgWYOuCMw4ccMU11Kf5oo6FqWeIejq318FvLvvbu4tXQxIyDbWOLJujz3UzQ6zILp2FVm7VpXjwWpmbSXHmDEiX3yhPsc6NkgJv/32IsuXe48LCw/YoNTzUaH7Bj9v773NPaVh2DAn9wryl5izk0oRAxKyjTWOrMN9jje76M45R62zA2gZOewwdS6fe857XFj4XneysEIgqFR/V/Wz0PpipprHAnzffut8jWnIYcahoIUCeVDCOOAAlZK+GDBmR79HbBhIXcp4jZJtrHFk3X338WYXF5zLfLVt600O5k4sh9ar//s/5+swRo/2PrC7dfOOMUGQcvPNztdYADDMGJQ33ghe1Rj+8x+RU081S4sLuVrMgE4HYth0inq0QGEV45deUknb6hNeo2Qbaxwlgje74kLXTpixIDogQXfEr3/t3de+vciGDd6ybDDo0x2QPPyweUT8kCofY5IaNKjbIlNMP/uZyA47iOy8swpAdHeO3+Y+J/lu5t/GBl6jZBtrHCWCN7viufdekV/9SuSss8w9dbkfmD//udOlgjESyOYaBQYlYw0Y98MY3SdobcE4F7/WkIMOEmnZUmTp0miZef1gDZ4ZM8zS4sAgarw/zL7BIOF84XsnThR59lmRAQNEzjtP5Pzz1cdWrUQOPlikSRORd981v7P4eI2SbaxxlAje7NLBHZDkOxbFD7pjdFBy+ukip5zivwDf0Uer7YYb/BO1uaHrA+bMUdlYTRhnYq5eXExvvqneX9hp0vnC3+Xii83S4uM1SraxxlEieLNLB3dAMmRIYf/tmxCU6J+NFgA9Q0jDtOV99hHp29db7mfWLJF77lGfYxwKZhilAdbMwTncd19zT3z+8hf1O9A1ZROvUbKNNY4SwZtdOrgDEnQVuFe1RYuEe1ZMPh55RHVr6N/hhm4c7Ae0oGDMShAM3sUgXMAD2mZLSC5nn63eW7t25p74NG1qf/FIXqNkG2scJYI3u3RwByTmANaLLhJZtsxbFhVyy2Dgqf4d7nEk+H3bbadm5wwa5F2N2A2DVRctErnxRnOPI2pyuLiddpp6f8irUwzz56vWJJsDd3mNkm2scZQIfbMbONB/sCPZ4Q5IMCg1bsgzMny4yJFHiuy3n/o9mJ2CGTiYCotZKmhhCOqGwUBPZEm97Tb10PeDVpw0JGTDAFS8v8pKc0/hMLX5ppvM0uJiQEK2scZRIvTN7oMP4h23UE4wBgPdLIXI1kISN3Sz6N+lH66YrostCHKQYAO0kqQdum3w/vbf39xTuH/+025LEAMSso01jhLBm108Cm3CtxmQgM7ZgY/HHac+Nm8uMnOmmj2TDcaZ1AeYTVSMlhIsIDlypFlaPLxGyTbWOEoEb3bp4A5IzP++kQsjTHI1t9tvdz5Hi8aCBepz5CLR0DWkf6euBtXVqvsuKCh59dW6CwGmGVqu8N4w++aJJ8y9+fntb6Nnzy0Er1GyjTWOEsGbXTq4AxKM89ArL+cLYzn0wnOYEvv662r9nS5dRK65xjlu8GA1+wa/H2NKLrhArc3jToWPVOt33y0yapRqQfnoI1WOlOymQluKigErCiOjK97j9debe6P73e/Uz7I15orXKNnGGkeJ4M0uHdwByeTJ5t74vP22yL/+5S3Dg1X/bqx7gxkq7oXx0IX03nt1X5e5IB+CHtsDPqNAJlp0U2GAbiEweFenqHe3OBULr1GyjTWujGzatEnmzp0rS5YsMXdlbN68Waqrq2Up8njnKPcri4I3u3hhwCNWo43KHZBMmODdh5TmUSF5FwIJZDE1E6EBgg79MEXOEyxCh9+N1+G3XktQl9Hf/66SrdUXSMuP9+n3HqNA6xNSypvda8XAa5RsY40rE/Pmzau5kbWS3r17S7t27WQI0nK6IMCorKyUnj17Svv27WUU2skDyv3KouLNLl4IHgoNSNAVooMIPPj+8If8u0LQGoIZVIB1WHTa988+Uzk19OfnnqtaQDC4Fa/FTP6F/aUCXVJ4j+ieygfO6UMPOX8zdHUVE69Rso01rkwgiFi3ZX351atXS+PGjT37p06dKt27d898vqHmX9wmTZrIxo0bfcsn1zxBzDIcG4W+2SG9OCUn26BWjPFAy0uhMK0bgQ66WvBQff55VY7qqFtL/vY353V8+qnzvRrqCY71W8OmPjnmGPUezzzT3JMb3r/OpHv11dEXP4yKAQnZxhpXhmbNmiWdO3f2lKHFZDgyWG3RoUMHWbRokW95//7965Th2Cj0zQ5ZOvUgSLIv27RfzBDJd5VZxKe6ewKLw7VurWbQuLuFevdWC9PtsYf6GtOAd9lF/eePsRLbb+8cO2aMSFWVs55NfYbzivMdZg0fbdw479c4vzhfxcSAhGxjjSszaB1p27atvI1Rhi4DBw6sCQ5G137dtWtXmTFjhm95r1696pTh2Cj0zQ7N1wxI4rd2bbh1aLIFJMiMWmjq+FzQReN+7iGDq/s15TlEKfV++lP1/u66y9zjDzN23BCQIKFcMa8dBiRkG2tcGVlb85Tq2LGjTJkyxdwlI0aMkKFDh9Z+XVXz7+jChQt9ywcMGFCnDMf6wU0taANMDaX8IIAwZ65oN9ygEnTl4n74m71ue+0VvL5MXPTUXw2tBvgaG/ZhyrCmp/2WCqTNx3l/5hlzTzhYbPDLL83S+DAgIdtY48oExo+gm2bSpEnmrgzMvunUqVPm81WrVklFRUVm3Ilf+Zw5c+qU4dgo9M0OzfbF/C+vlGE8AQY5BgnzJ3EHJAcc4M1DgpV4tww7KhpM2d1pJzVW5dZbVYvBbrupdW9+8hP1+ho2VIN2L71UfU9NTJ0JlNBrGJRIrb7Qs4zw3qM6+ujiZm5lQEK2scaViUGDBkmDBg2kZcuWtdsCnUZzi379+mVmzbRo0UImuDr7/cr9yqLgzS4d3AGJOZj0scdU148NNXGvb0ZT/dpOPNEpw7gWtOZgohimF9d3NfF85j1GnRKMFiwMki0WXqNkG2sceaAlBflKTH7lfmVh8WaXDu6ARM/g0F54wft13PD7dJdd0O/Srw8fzYGdbsjqWp9hSjDe5+WXm3uCIRcLBgTneQnmxGuUbGONo0TwZpcO7oAEYzZ0IjJ0o/3qV95jC4XuHz2rCmnQkYekTRu17+CDMfXcezwG1CIXCmDWDTYkVkPm0+nTRfSwJQzeRStDfYeZSPh7hJ0KP3u2Oh4Dg4sx+JjXKNnGGkeJ4M0uHdwBCRa4K+Z4HsyY0fk39HosaBkAdBf5LZ7nXrcFAZN+rQ8+KDJ+vLOvVOgpwVhoMAwMiMXfEGvmxI3XKNnGGkeJ4M0uHdwBiTnttxjQCpJv0IMWEv1a//Mfc2/pQMIzvMdTTjH3+MNx+DvGPcCX1yjZxhpHieDNLh3cAQm6TcLMzMkGU43dC+RpL72kUskjhTwGpC5frrpskHfE5B5Ii+AFuWpeflm1kGBDYIKEaliQ75xznGNLSbt26m/So4e5p67PP1d/R51gLi68Rsk21jhKBG926eAOSDD91p3XApOwCglQ/vQnZ9Dq3Xerqb2YLY605wcdpFbyxUyR229Xq/Xuvrv62nwIY3zErFnemTZICqbHT5QqnCP8XTDeJhek1I97YC+vUbKNNY4SwZtdOrgDEnNJe6Ruf+UVb1kur76qHqTHHmvuUTAAFYsA6oX3kPukSxf1OVo9jjrKORZOP119xPgS90J/yO2nX3exc6XY8tVXIitXesv0mJInn/SW28BrlGxjjaNE8GaXDu6ABK0N7u4SzHCJsoIwgo0dd3R+nm5t2X9/1fpx333qa3TVNG2qWl+wBJIeU9K/v7PwnqaDJMyi0ZlkUYbgBD9XrxJ8/vm132IdBt5icK6Z6VbDOQyTUfXZZ72ZcXVw0qiRam1CcOYe5OuGBbdxzrHmT1x4jZJtrHFkHW60+maHDJzuQY5LljifU3zmzTNLFHdAgm3mTGcfMqRiamlYGBfiHnjaqpVqycCCedtuqwZf/vnPzvEnnaS6iZARFh5/PFy3AxbfQxcPYByJ/n3otjCTu8UFLTRo9dG/F/A1upvAWBrK47bbRLYsjp0JuPR0ZXRFjRjhHOeG84AgRMPvP+88FYRtSZLsgZWU8bfEFhcGJGQbaxxZpx4gW2VuohiIhwcZ4L/zq67yHkuFw3/dQc8W5BrBPjzIzP/woz7gEBBg0Tj9fbp1Ay0hQTN48HDG70VQ+stfqum8uZiJwJBgTQdWhxyiBnnGDd1K+PnuMSuTJ4t8+KHzdRjXXacCL0C3VZj1htwQxAS1tuC14DW2aGHuyQ8DErKNNY6sO/JI3uxsQmCAcx4VxpBgPZn64Prr1RiUcnfuuSq1fhx4jZJtrHGUCN7s4pHvLBizlSEbjNdAq4PZgpJWaG3J97yQg9co2cYaR9ZhPAFvdoVDqvF8TiPGWUT5Pvc4E4wRMSHVeb7JzooB3Sp4rRddpL7GlOOdd3aCMKwwrd+PuenuJveGcRv1BcbtROlmy4bXKNnGGkfWqdVN4696Uf7rLwWTJqlEYVFh7E7Q6UfLgtm6oAMSfPTrwsHsGfeU3KQ995x6vVjzBtBl9dpr3mMwVgnjZ/bZx9l+/nM1SwiDZvE+saEbaNdd0/X+ssGMHgwUXr3a3BNdMa5RomxY4ygRcd/srrzSWaiN8veb3zi5P0jBrBbM5ik3cV+jRLmwxlEi4r7Z4SGKJFtUmCOOUBs5MGulbVuztPTFfY0S5cIaR4mI+2aHboPLLjNLKaqbb1YbOTCDR+cbKSdxX6NEubDGUSLivtmhSR2Lw1Fh0BIQlPa9XGGg6AknmKXphYAS+UoKFfc1SpQLaxwlIu6bXVWVyJ57mqUUFRKLYSPHgQeKNG9ulqbX4sV1BybnI+5rlCgX1jhKRNw3u7POYgtJHI45pjzHS2SDMTWVlWZp6Yv7GiXKhTWOrFPrnWyVmf6LqaTIS+KG1g4zF4Q5bdOE9OFBKbXz9dZbag2StE35/OgjNUV14ECRHj3MvbkhwRnOKfJ1mFk9v/5abW4IUHC8uRKvNnas2koVFrSLsshgqWBAQraxxpF17dqpmx0eqOhmwSJrbuPGqURWOhcE1rsJWuW0mPA7k1j2PZd169QA3jffNPeEN3iwyPTpZimRgwEJ2cYaR4ngzY4o3XiNkm2scZQI3uyI0o3XKNnGGkeJ4M2OKN14jZJtrHGUCN7siNKN1yjZxhpHieDNjijdeI2SbaxxZWTTpk0yd+5cWbJkibkrY/HixbJo0aLM9sUXX9SWb968Waqrq2Xp0qVZy6LgzY4o3XiNkm2scWVi3rx50qpVK+ndu7e0a9dOhgwZ4tn/ww8/SMOGDeXss8/ObHfddVemHIFHZWWl9OzZU9q3by+jRo3yLYuKNzuidOM1SraxxpUJBBHrkMCixurVq6Vx48ae/StXrpTjjz/eUwZTp06V7sgOVmPDhg3SpEkTmTx5cp2yjci2FQFvdkTpxmuUbGONK0OzZs2Szp07e8oWLlwobdq0kXHjxsm0adMyAQygJWX48OG1x3Xo0EH69+9fpwzdPFHwZkeUbrxGyTbWuDKD1pG2bdvK22+/7SlftWqV3HHHHTJmzBi58MILpWPHjpnygQMHyujRo2uP69q1q/Tq1atO2YwZM2q/dsNNre62MvabHYbFGG+p5H3wgUizZt4yZJedOdNbZlq/3knJb6aOR1wZFFsi1Ty+p5ysWaOyBb/6qrknvfA3wuVV6N8q7muUKBfWuDKydu3aTKAxZcoUc1cdzWqedBiwOmLECBk6dGhteVVVlQwYMKBOGVpYwvrd7+K/2T36qMgNN5ilpe2zz0Tat/eWITD75z+9ZYBjFyxQn//wg3pYbb21yPLl3uP+8heRv/7VW6bhwbzNNmZpadu0SWT77UUaNTL3pNf/+38iy5aJbLcdrnlzb3hxX6NEubDGlQmMH0E3zaRJk8xdGe4xIGtq/i1sVHMHXl/zrzRm5XTq1ClTjlaUiooKmTNnTp0y3cUTxsMP173Zoez990VOO03kgQc8u0LB27rjDrO09OG84aGZC1qPpk41S+vCuX/wQbNUwT4EM37uv98sKR1ohdptN7M0/bD69bXXmqXhmdcoUbGxxpWJQYMGSYMGDaRly5a12wL9L3ONiRMnSuvWraVLly6ZfRMmTKjd169fv8xsmhYtWtSW+5WF9fjjdW92KMMqtt26qdYOCgfnLQ2rEUesAmTBkUeK3HyzWRqeeY0SFRtrHNX6vubJ9rW59vwWaGFBHpNcZWHxZkeUbrxGyTbWOEoEb3bp8MknarxBXEaMENlrL5HDDjP3lCcMIC5kHEeSeI2SbaxxlAje7IoHD8GwDVeYORM0UBXdQVFMn46/q7PlMxYoqu++UzOFQozTjhUCrh/9SOTcc51xNRj0izK3554TGT/eWxbk1FORL8gsTQ6vUbKNNY4SwZtd8fzjHyKTJ5ul/jDbxnyIwrffivjkycvqySe9AQl+9ocfmkfFa+RI5MFRY4/GjRN54QXziOL45huRH/9Yvc933lHv87//VZ+b/v1vZwru4MHBA4OhXz+R+fPN0vzg9yA4zRevUbKNNY4SwZtdOgQFJPl4+mlvQIJpxfj5rhniRbVqld3ukYkTnXwuPXuK3HSTGphtwrJQe+6pPsf0YdcyUXUMGyby3ntmaX4Q/Oy3n1kaHq9Rso01jhLBm106xBmQPPOMNyCZNUtk993V57/9rXl0PNAygRaIpLz1lhOUNG0q8tJL5hEimCGvp6RfemlwqwXeC7qCMAMfaX2QE6YQ6MrCuQ/bWmbiNUq2scZRInizSwc8SDGOJA4YL+EOSN58U5VXVKiv8TFuSPiGJMFb0uIkRmdHRdI/E9ap1MnqmjcXCVogG10sSGo2d64KSgIW5Y4E5+bTT83ScHiNkm2scZQI3uzSIc6A5NlnvQHJ7NnOPqzFiLITTnDKSgUy2zZooN7fttuq1PtoOUF+mCeeUK1EF10k0rChSlSG8Sdu+HrxYpVpuG/f+MaQALqKWrY0S8PhNUq2scZRInizKxz+o8aDrBBxdtmYY0jeeMO7/9e/VuXnnectr+82bFAfMUMG5xLvER8//1wNuO3RI3vQh66uXr3UGkTZBrzmAzOfnnrKLA2H1yjZxhpHich2s0Mf+7p1ZimZkMPu+uvN0miKGZBgDInpoIPUvpNOMvfkB10bNgey5oJWELw/nFcdLCJbqp5lk2s69r/+FZwX5sILsTimWZrda6+JbFknM7Js1yhRMbDGUSKy3ezwH+Mjj5il2T3/vMjw4WYp5WI7IIFWrdR+PGALhfV5MM05LcsNYOpxmzYiO+2k3uOvfqVm1ugxLj//ucjFF3u/xw0BZlCQifcZdVozBsqiGygf2a5RomJgjaNE8GaXDnEGJOYsGz2o1c/eezutBnG47z6zpLj8ZsqsXOl8jtejW0pwfjt3VuVXXFE3ELv1VpXDBGNOsEgkxuL4QUAStC8IphgHtbjkwmuUbGONo0TwZpcOxRzUmi0gAfxudHEUAjk7br9dZMUKc09x7b+/yC23OF+jW+aYY7zdR+iewbnF+8QUXuQoQcsJ0uu7IRj58kvVxXPbbaoMg2Avu8x7HN5r1CnOCHBOPtksDYfXKNnGGkeJ4M0uHfDfOxKYxeHVV70BCfJgZHPvvephfcop5p7wkDL+jDP8E5LZhJk1eB2uBbQzvvpKZPvt1flAq9Dll6sEau7Bq/j6zju934suSDOnCaYLY/ZOWK+8os5xwHqZOfEaJdtY4ygRvNmlQ5wtJGYekupq84i6dLp5PKzrE4xdMYMPQMDw4IPesp/+VL1HnOuuXcMPTF2zRk0FzheyvuL35YvXKNnGGkeJ4M0uHeIcQ2IOan39dfMIfxhTgeMHDTL35Ia1XzCI9MorzT3FheRlmMZs5hRBsGEmNKuqErnmGqf7BtNwMV5Ew9gbTBHGANQ+fZxytKL4tTKh9cQ9XiUIfseQIdHXJNJ4jZJtrHGUCN7s0gHBCBaJi8PLL3sDEr8WhCCYeYLviZrNFa0yf/pT/WhhwVR2nacEA1k15ArRAYbfYFkTxuoggMH4kFwQIDVrll9+E16jZBtrHCVC3+ww+yCfmyXFo5gtJEHTfoMg7Tq+D1Njo9puO7Ok+DB7BWM04IILVECGNXvcrSbotmrdWuSPf1QrKB91lHqPGFeCBHFo3UFCNcDMnIED1crFJrMlBt05/ft7y+LGgIRsY42jROib3eOPGzvIqjgDEnPar5mpNYxLLlHfi3EXURx3nFlSfGjx0DOJ0FqBAaxm4rPvvnPS5run+6L7BmU4/9OmOeVIG49kZib8DHfgjtk8SFNfzGCeAQnZxhpHidA3O/xXWMybainDefv4Y7M0mmwBySGHhB+ACRgbUUgLiYZEXvj+bbYx9wTDw1mncE8bjI3B+8HUYA3nHOceW74DV7FGDVpdioUBCdnGGkeJYJdN4TCd87rrzNJosgUkUcaAABaSiyMgAd3agq6NbNCagGP1arppgtaSF19U40PQIoLuHKzmi3EiCLZw7hFI4SMW3UNrIQa/zpunumgmT1YDU5HczM9DD0ULGKNiQEK2scZRInizSwcEI3FN+8UD1B2QuFf7zQfSsONhfeSR5h4HHvaffVY3+2kaoLtGjzEBBN5HHy3y/vtOGeiBrghONASbCLTQAuZuWbGJ1yjZxhpHieDNLh2ytZBEZY4hKaSFRMPCcPhZhx1m7ikdOiDBZgYr2SBRmnv8Sdx4jZJtrHGUCN7s0gEBSVzTfpEfwx2QzJljHpGfY49VP88vzTwytfbunb0VxSbMfgHMuFm4UH2OlO8YhNqrl8j55zuzagDdN3jtGA+y664iu+yipjKjtWniROc4PwhI0CWkf6efu+4S+cUvwk0nNvEaJdtY4ygRvNmlA7progwezQapyt0ByTvvmEfk729/U8GTOSUYrTAYGI0ZL2mALhkEFJh9owccN22qul02blSJytzjPtCNg5wkehzVwQerFhN03+y7r3NcEHzfWWeZpQ78rh49gsehZMNrlGxjjaNE8GaXDnjIxzWGxEwdn2txvaiQAA0/t5S7b1q1Uu8Rf5ewawyhpQQ5YIIgsAmTRM3Ea5RsY40rI2vXrpXZs2fX/Oe25V83w+bNm6W6urrmBrc0Z7lfWRS82aVDMQOSMGvZRIX/9vGzkWAMMKBVd43YhNTxM2d6y9BN88gjIoMHe8uRnyQXrPSLgawNGjjTgdFK8umnIh98kH0mGgb2YsXgIOiS22EHszQ3XqNkG2tcmZg+fboce+yx0qdPHznxxBPlj0gd6YIAo7KyUnr27Cnt27eXUaNGBZb7lUXFm106xDmotdBMrWGdeqr6+b/8pciZZ4ocdJB5RPFhdV+MkUHVRxp3wNo9NZdFBqYh61lGmOqLgAOw2B3GkSB7q9tjj6lxI+3aqfVsbrtN5Gc/U38fjCvJFnDkgsDwnHPM0tx4jZJtrHFlAkGEtn79etnbWPxj6tSp0h3pIAUJpjZIkyZNZOPGjb7lk2vunGYZjo2CN7t0qI8BCVxxhfodeFgjhTuSuCUBQQcGrJowkHX5cvU5ukv0wOGTT1bn+/77nWPHjnUW5MMqwng/Gv4+eJ8YJ2Mbr1GyjTWujCAQ+bzm37nBgwfL73//e8++IUOGyPDhw2u/7tChQ81/gYt8y/v371+nDMdGwZtdOugHHjb37I984OHsDkjQM4ifny0NPFoUrr7aLA3njjvUz8eDf7/9zL3pgQRpuoXjpJPqng+MtdEDXbF6MdLHazhW/42itnIgiEHAhq6gnXfO3u3jh9co2cYaV0aervkX9uSaf9GOOuooeQp5vl0GDhwoo0ePrv26a9euMmPGDN/yXr161SnDsX5wUwvaKHnugASNXMgQai7kFhamlh56qPPz8CDebTf1uR4T8dZbKpDA2A9Apl6kVg8j7CDPJGBBvL32Evn3v809hcOlpc9plCFbmA6M865bk9q3N4/Ijtco2cYaV4YwuLVRo0ayYsWK2rIRI0bI0KFDa7+uqqqShQsX+pYPGDCgThmOjYI3u3RwByQ77aTWo8Hy9vnA1FvM6HAHODBhgho427y5yreB7gu0BGC1Wp1SfdgwNTYkG6yoa3rhBZVCPWl4r6tWmaXKAQeolYyjtlC44Rzq84opvBhfUnMJ54RzhmndO+4YPd8Mr1GyjTWuTGxyLUOKrpt9a54c7oBk7ty50qlTp8znq2rurBUVFZlxJ37lc+bMqVPmHqMSBm926eAOSNq0KeyhidkeeOjpn4fU6W4Y/OlO9oXcHPpYTOlFEq+333b2B3nwQedztMpgPR889NMKuUcQFCDfCAIoPbYkKowv0ecLPa6uy9cXZv2MHCly1VUqMMH3YRZQWLxGyTbWuDJx9913Z2bZnHHGGXL44Ydnvjb1q/m3FbNmWrRoUfMf2YSs5X5lUfBmlw7ugKTQ1XLRQrL//s7Pu+gi8whHVZVqGcHYD7SS6NeB/+QxwwTw/e4Bo+eeqz6OH++UaRHj4djlCuQQlCC1PhKf5TEpLQMr+6LLS5/f3Xc3j/DCwnwICjEbCd+Lc4wtLF6jZBtrXBnBjJjlNf+efY87VYB1NU8Vd2tKtnK/srB4s0sHd0CS55+yFh56P/mJ8/Pc2VPxQEaacxNmoGC9GgQl+J4DDxT58ku1Dy0uyPWB6bAIlrJlLr3ySrPELrTShGndQcCGsSZ6GrCG96rTu6MLym+GUk3sL3vuqaYM41yFCS5w3tC1U1Ghpirj+5DCPgxeo2Qbaxwlgje73PBgRjN9MbkDkiefNPdGg5kiGHjqF5DAv/7lTRL26KMqrToexoCGNnOyln59+HjDDd59gHOEBenM70srBH0vveQte+01FTCcfbZatwcDV/3Wp8H/ETpowTRhnUTNNeEtJ8y2wffk6u4BXqNkG2scJYI3u9xsByT43L3OSlSYnYOxEvrnobvlP/9RrSCY3vvEEyJ/+IOzxguyneYaaIlAA4NF9es84wz1c3T30htvqBaUUoBEa5jmi+Rov/61GmCcDcaR63ONLLlhvPqqOj5Mdl5eo2Qbaxwlgje7dHAHJFl68kLBf/VYsVb/PPw3jxYRpD9//HEVRKA1xD3eQrd6oOXAbxwI9iPJGPa5W0uQTr0UYbZRs2Yq0du115p760IOE31eMEbFD7p50C2mHXecOt7vfLvxGiXbWOMoEbzZpYM7IMn1gMoFQYL75+mxIGE88EDdFgEEIpiKjLVhkHQN3UH77KN+tl7hFqv/4ncGTbm1Aasao6UHeT50DhfkWkEqeQQMaMnINehVw0wjDHrF+8PMIw3jU4JaQdyrLG+Z/OaB9X9wnvB6AH9nHJtrNWZeo2QbaxxZp1KMe6seyj75ROS88wofy0BeeAAFrQaLh3lcqeOxfot+MGKLe7VfP0i2ht+F91HIei+FwIMdv//oo1WiNwzuxfo1mG6LAaxHHqlW8QXkeEHrR64EZ5dfrn6mhnTy777rfG3CjCV93v1ccol31hPG7mDMSjbmNUpUbKxxZB3GAZg3O6ztgRsu1voYMyb8f5QaElOZeS9KHfJZ4Fzmgu4QtEBoOLfbb68Slf35z951UtDNogdOapgZgqAFD1d0IzRsqDKw3n679zh0+eBYPEixRf0b5gPdQBgUikRhuhsH7xevN471X/DzMf7FPCembCv6YgAqBvTCRx+pHCJI554NgkgE6FHgukE+FuSTQUuNmTUWg3+1ww/PnW7fvEaJio01jqw7/vj4b3aYuXDvvWZpaZs3L7+EYAjcdNeKmaQLYz0waNQNS9fjWAyExCwN/Z84ApQ0wsBcvL++fc090WFgLgKwoPEZaYT0/cjm6l7Az4R09Ph7Zgsa475GiXJhjSPr1DiAeKsellivTw+NOKBFIEz6cBP+29etGAhq3F58se60VEDrB1od3N0Ibo89ZpYkL9sDOSy8b4xjwewUBGv1BbqNsKC3XjPIhFYlBCR4b7r1xhT3NUqUC2scWaeSZ8Vb9fCfvrtJuhzoRGTu5erDwH/FCC7QbYP//hs3dvZhRozfQwzfg9+FzQ+6cAoZ1FoMyHMSh9/8RqW9x3s0oQsGg2179xb58ENzryNXF002N90kst12qksmLHQ16QCyTx9zr7peMMUYfyfkM/ET9zVKlAtrHCWCN7t0wEMLuUM05AzBlg+0rOhgBBvyjJQ6BGoI6hDcYYAtVmRAtll34IZU9zqY2WUXtTJwtuAZY3Ouv94JDDE+Bd+PvxWyrGJhwmuuUa0cgG4lv8yuzz+v/g4ISjCdGINuETz94x/OMQhqcZwfXqNkG2scJYI3u3TAQ86dnAx5QnTm1KimTSu/gASwgjGCjDPPVANYsZDdaac5+//yF5WBFRC4nH66Os4Pcobo84cAB9AChoGqSJiGAbG33KI+18tRqVlr/oO6kYROd8+hCwcfkeZeQ8tN0Kw2XqNkG2scJYI3u3TAAyquab/II+IOSPz+a4/biSfmN7A3DmgdwYY1ZvQCgSNGeI9BN44OCDCrpW1bNeA0CAYQI+hAMIeF9AAJ59CyorvDEKhgvJS7lQV5WBC4uFs/ECDec4/KTaL/JhhzZI4bCsJrlGxjjaNExH2zQ3KoOAYxlps4AxL8p+0OSJDWvdgwQBf/5Zvr5tgwcKBI9+7qcwx+xcBXDa0WaDVBbpaTTlIBC1pScF7QRYaB3fpnuM8TkphNn65S7mNhQTdMBcZ+/Az83W680bsf58AdBOL1YOYZjkeLjf67oLUljLivUaJcWOMoEfpmh9wjhWYIBdyMk8zWWQwYMPnww2ZpvOIMSHTXgd5stJAMGqQe0l27mnvsQnfJzTerFgqdrVVDMKLrOM4LFtLDQFkNLUtoAcHAVfwtzPwhfvAz9QrJ6MIJolto0C2H5Gz6b4P8M7kwICHbWOPIOty89c0ON9N8xyyUOgw4DLOkfRh+4wsgzoAE3Qi2A5JTTlGZT4Peny1YcRevBQv99e+vyjDTZfRoNcC1Z09VhhwvZu4PBFO4DlR+HrVhOm6urigMdtXH77ijGgzrdx4wywbBLTz4oPM9OpFcUAI2BiRkG2scWaf6wln1CoVxBmGSk2Eqr3vgKmZlIAjBeikY05Ar2yu+F38z/Eeu15ABdNHolXsBS9rr/9qxYeyDLXgg33qr+hwtZfj9Qe9LrxQcFtaXyQemPWMAqxsy5pqLGKLlAmnvkQ8EaeYxGwYDT5H1NgwENAh6zj9fnX+ci4kT1Ucwu7MQgOD84G96wgnqc9QHLOrnTr/Pa5RsY40j6xiQxGPxYnUu3UFBEGQv1fBA1kHDokVOeZBzz1UPNzMxGjKh6u4Fd8CTBPf7A8xmCZrlg6AK6eXD2mMP1eWiAzN0r+CBn0Y4D0OGqLwpeK16rIoJdQYzf3AMzpUeMIsNrTrAa5RsY40j69SNj1UvDjiNmE2RDZJg4SHqB+MXcn1/GHgdGGB63HHOwy3KQ79QCJiiJA7Lh35gIzDDRzMICgNrLuVKaf/QQ/l3d+G1NWmiPtev1Q9aUhBg6fekVxbWX6vPeY2SXaxxZF0xAhIsT4+VVMsJ3vOuu5qldWEWh/t0o8sGX6OZH4OK3YOBX35ZzVhy0yvq4uGGtWw0NPffd5/6XI9RQJM/HohY9dYmBEEdO6rP1RglbwZaN6woHSVYQrcYZqzo91ioXK1SWAk43+4u/XfFFGBcD0HdUzpoxMw0fMSA2ueeU4sU6iAm7muUKBfWOLIOg+nivtlhgN9775mlpQ2DI+fPN0vDuewykXHjzFI1vsBM2oWysWPVgEhsWtDDLg0w2yVo1pU5CyYXTMH1S6efVmEHiWNasRmYIUDVY07ivkaJcmGNI+sw2JA3u8Ih/wb+248Kgyr1VFAEcmFhzRN8jzvfRjk49VSzpDTceafICy+YpSJz5qiPvEbJNtY4sg7/2fNmFw9zCmlYGO+BwATTTf/+d3OvPxyPMRBBME0ZA0AR7GAzZ3cUE8auFEu+5zhfTZuKHHWUapWKA6Yk+7Xw4H1hWvk55zhlOE53t/EaJdtY48i6ww+P/2a3cKHKcFlOXn893BgSEwZjImDAgm1Tp3pTiVdXq81Np0XHhvwYSEXuBw845ODAGASMSbDJnUAOXTKoXnElS8OCdqNGZV/NN05//avadEtFVBjzgqneGq4Ns2tGQ11ACxDGnQwYINKhgzPoNu5rlCgX1jhKBG523LhxS/dGZBNrHCUiLTc7vg4vvo660vJa0vI6iIqFNZwSkZabK1+HF19HXWl5LWl5HUTFwhpOiUjLzZWvw4uvo660vJa0vA6iYmENp0Tw5urF8+HF81EXzwmVOtZwSgRvrl48H148H3XxnFCpYw2nRPDm6sXz4cXzURfPCZU61nAiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDErJm8+bNUl1dLUuXLjV3lZTvvvtOli9fbhZnBJ2DqOX1zaZNm2Tu3LmyZMkSc1dG0Pv0K/crq2/Wrl0rs2fPlo8//tjclRH0Hv3K/cqI6iMGJGQFbpqVlZXSs2dPad++vYzC8qklaP78+dK8eXPpq5dMdQk6B1HL65t58+ZJq1atpHfv3tKuXTsZMmSIZ3/Q+/Qr9yurb6ZPny7HHnus9OnTR0488UT54x//6Nkf9B79yv3KiOorBiRkxdSpU6V79+6Zzzds2CBNmjSRjRs3GkfVb99//71069ZNRo4c6RuQBJ2DqOX1DR6a69aty3y+evVqady4sWd/0Pv0K588eXKdsvp2TnA+tPXr18vee+/t2lt+54NIY0BCVuC/4uHDh9d+3aFDB1m0aJHriNLx1FNP+QYkQecganl9NmvWLOncubOnLOh9+pX379+/Tll9PCcIRD7//HMZPHiw/P73v/fs83vfpX4+iIABCVkxcOBAGT16dO3XXbt2lRkzZriOKB1BAUnQOYhaXl+hdaRt27by9ttve8qD3qdfea9eveqU1cdz8vTTT8vJJ58sRx11VKa+uPm971I/H0TAgISsGDFihAwdOrT266qqKlm4cKHriNIRFJAEnYOo5fURBnF27NhRpkyZYu4KfJ9+5QMGDKhTVl/PCeC8NGrUSFasWFFb5ve+y+V8UHljQEJWYIZFp06dMp+vWrVKKioqPH3ppSQoIAk6B1HL6xuMH0E3zaRJk8xdGUHv0698zpw5dcrq2znBjCMNXTf77ruvJyDxe9+lfD6INAYkZE2/fv0yMwFatGghEyZMMHeXjKCABILOQdTy+mTQoEHSoEEDadmyZe22YMECzzFB79Ov3K+sPrn77rszs2zOOOMMOfzwwzNfm4Leo1+5XxlRfcSAhKzCf8vu/xDLUdA5iFpeaoLep1+5X1l9ghkxyFWDmVlBgt6jX7lfGVF9w4DEB5b55laaGxERpRPv0D744CpN/LsSEaUX79A++OAqTfy7EhGlF+/QPvjgKk38uxIRpRfv0D744CpN/LsSEaUX79A++OAqTfy7EhGlF+/QPvjgKk38uxIRpRfv0D744Coc8ivMnz9fPvjgA085ciUg4+SSJUs85bBmzRp54403ahcH8zv2q6++yux3b3ol2Vz4dyUiSi/eoX3wwVWYb7/9NrOIGjJIYmn0Sy65JFM+b948adWqlfTu3VvatWuXWb1Ue+KJJzJZK6+++mq57777Ao994IEH5Oyzz67ddtttN5k9e3btz8mGf1ciovTiHdpHqT24vvnmG5k5c2am9eGee+6R6urqrOXwzjvvyEMPPVS7LVu2rHZfLi+//LL06NGj9uuDDjoo8xFrbOjWDKz82rhx48znWGCsadOmmdejBR3r9sknn2SClrBK7e9KRFRKeIf2UWoPLgQXe+21l9x4443y6KOPZgKEd999N7Ac3n///cyaLOPHj88c8+mnn9b+PJR36NDBsy1durR2P7pejj766My6GjfccINcc801tfu0WbNmZRZcg2eeeUauvPJKefbZZzOB0Zdffhl4rNtVV10lY8eONYsDldrflYiolPAO7cN8cF1wAR6ynqKCvP22yM03m6WF6ddP5OOPzVIFgQdaILTrrrtOHn744cBytyuuuELGjBnjKQsDC6q1adMmE+S89tprnn1o8UCXzts4EaKWW2/evHkmuEAAhNeEVhO/Y7WVK1fKoYceGmn9DvPvSkRE6cE7tA/zwVXzTIxd3D8z289D4FFVVVX79YABAzLjNILKteeff15OP/302q81tGiccsopns3dQvLkk0/KxRdfnPkcg1ARbGBpdECg0bFjR5kyZUrt8WgVwe/WzjzzzEyriN+x2s01Ed3QoUPN4qzMvysREaUH79A+Su3BFRR4BJXDihUr5IgjjsisSBrVsGHD5M4776z9Gkuj69kw6HqZNGmS62jJDEpFt4925JFHZrqM/I4F/JyKigrPmJMwSu3vSkRUSniH9lFqD66gwCOoHPr27ZvpbunSpUtmmzFjRu1xuaBVpLKyUrp16yYnn3yy9OnTJ1OObpwGDRpIy5Yta7cFCxZk9l1//fWZ14LumVtuuSXrsSNHjqz9mVGU2t+ViKiU8A7tgw+ueGBw68aNG83iQOii2bBhg1kcG/5diYjSi3doH3xwlSb+XYmI0ot3aB98cJUm/l2JiNKLd2gffHCVJv5diYjSi3doH3xwlSb+XYmI0ot3aB98cJUm/l2JiNKLd2gffHCVJv5diYjSi3doH3xwlSb+XYmI0ot3aB94cHErzY2IiNKJd2giIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJKHAMSIiIiShwDEiIiIkocAxIiIiJK3P8HOogJvGLj4NcAAAAASUVORK5CYII=", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%plot --format=png\n", "PlotSparse(C)\n", "isoview on" ] }, { "cell_type": "code", "execution_count": 27, "id": "d5cfac5d-8048-46a5-b61a-bfe493ad9278", "metadata": {}, "outputs": [], "source": [ "n = size(C,1);\n", "d = full(sum(C,2))';\n", "d(d>0) = 1./d(d>0);\n", "e = ones(n,1);\n", "f = ones(n,1);\n", "f(d==0) = 0;" ] }, { "cell_type": "code", "execution_count": 28, "id": "abc77963-f38c-4e25-a434-bcd84549eba8", "metadata": {}, "outputs": [], "source": [ "ITMAX=1000;\n", "TOL=1e-6;\n", "pi = ones(1,n)/n; // distribution initiale\n", "q = 0.85; // parametre magique des auteurs de l'algo\n", "for k=1:ITMAX\n", " // faire le produit\n", " newpi = q*(pi.*d)*C + 1/n*pi*(e-q*f)*e';\n", " if norm(newpi-pi) < TOL // condition de convergence\n", " break\n", " end\n", " pi = newpi;\n", "end" ] }, { "cell_type": "code", "execution_count": 29, "id": "e615b4d6-7d2d-4646-bb5b-1882af9a3485", "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "02 0001 0003 0005001 5002 5003 50000.010.0020.0040.0060.0080.0120.014" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%plot --format=svg\n", "plot(pi)" ] }, { "cell_type": "code", "execution_count": 30, "id": "4b9c9003-7abb-4334-8a72-0d4cc77ded3f", "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "02 0001 0003 0005001 5002 5003 50000.010.0020.0040.0060.0080.0120.014" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "[pis,k] = gsort(pi,'g','d');\n", "plot(pis)" ] }, { "cell_type": "code", "execution_count": 45, "id": "1f57d6df-f8db-4bfa-a5f0-18a8948a6a7b", "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[4lSearch result for \"classement\":\n", "\n", " Rank | URL \n", "---------|---------------------------\n", "0.001221 | https://www.utc.fr/utc/lutc-en-bref/classements-de-lutc.html\n", "0.000266 | https://www.utc.fr/utc/actualites-de-lutc/breves/lutc-au-classement-qs-2020.html\n", "0.000266 | https://www.utc.fr/utc/actualites-de-lutc/breves/lutc-au-nouveau-classement-des-ecoles-dingenieurs-du-figaro.html\n", "0.000203 | https://interactions.utc.fr/thematiques/prix-et-concours/lutc-4eme-du-classement-de-lusine-nouvelle.html\n", "0.000190 | https://interactions.utc.fr/thematiques/prix-et-concours/lutc-2eme-ecole-post-bac-du-classement-lusine-nouvelle.html\n", "0.000189 | https://interactions.utc.fr/thematiques/prix-et-concours/le-classement-des-ecoles-dingenieurs-de-letudiant.html\n", "0.000181 | https://interactions.utc.fr/en/thematiques/prizes-and-competitions/lutc-4eme-du-classement-de-lusine-nouvelle.html\n", "0.000175 | https://www.utc.fr/en/utc/news-utc/news/lutc-au-classement-qs-2020.html\n", "0.000175 | https://www.utc.fr/en/utc/news-utc/news/lutc-au-nouveau-classement-des-ecoles-dingenieurs-du-figaro.html\n", "0.000170 | https://interactions.utc.fr/en/thematiques/prizes-and-competitions/lutc-2eme-ecole-post-bac-du-classement-lusine-nouvelle.html\n", "\n", "\n" ] } ], "source": [ "search(U,pi,\"classement\")" ] }, { "cell_type": "code", "execution_count": null, "id": "87778648-ed59-42e2-b3ca-be96962146d9", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "c7cc2e90-3855-4812-a985-bb464ac6ea73", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Scilab", "language": "scilab", "name": "scilab" }, "language_info": { "file_extension": ".sci", "help_links": [ { "text": "MetaKernel Magics", "url": "https://metakernel.readthedocs.io/en/latest/source/README.html" } ], "mimetype": "text/x-scilab", "name": "scilab", "version": "0.10.2" } }, "nbformat": 4, "nbformat_minor": 5 }