{ "cells": [ { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "function cls_solve(A,b,C,d) \n", " m, n = size(A)\n", " p, n = size(C)\n", " G = A'*A # Gram matrix\n", " KKT = [2*G C'; C zeros(p,p)] # KKT matrix \n", " xzhat = KKT \\ [2*A'*b; d]\n", " return xzhat[1:n,:]\n", "end;" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The objective function " ] }, { "attachments": { "image.jpeg": { "image/jpeg": "/9j/4AAQSkZJRgABAQAAkACQAAD/4QCARXhpZgAATU0AKgAAAAgABAEaAAUAAAABAAAAPgEbAAUAAAABAAAARgEoAAMAAAABAAIAAIdpAAQAAAABAAAATgAAAAAAAACQAAAAAQAAAJAAAAABAAOgAQADAAAAAQABAACgAgAEAAAAAQAAAS6gAwAEAAAAAQAAACgAAAAA/+0AOFBob3Rvc2hvcCAzLjAAOEJJTQQEAAAAAAAAOEJJTQQlAAAAAAAQ1B2M2Y8AsgTpgAmY7PhCfv/AABEIACgBLgMBIgACEQEDEQH/xAAfAAABBQEBAQEBAQAAAAAAAAAAAQIDBAUGBwgJCgv/xAC1EAACAQMDAgQDBQUEBAAAAX0BAgMABBEFEiExQQYTUWEHInEUMoGRoQgjQrHBFVLR8CQzYnKCCQoWFxgZGiUmJygpKjQ1Njc4OTpDREVGR0hJSlNUVVZXWFlaY2RlZmdoaWpzdHV2d3h5eoOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4eLj5OXm5+jp6vHy8/T19vf4+fr/xAAfAQADAQEBAQEBAQEBAAAAAAAAAQIDBAUGBwgJCgv/xAC1EQACAQIEBAMEBwUEBAABAncAAQIDEQQFITEGEkFRB2FxEyIygQgUQpGhscEJIzNS8BVictEKFiQ04SXxFxgZGiYnKCkqNTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqCg4SFhoeIiYqSk5SVlpeYmZqio6Slpqeoqaqys7S1tre4ubrCw8TFxsfIycrS09TV1tfY2dri4+Tl5ufo6ery8/T19vf4+fr/2wBDAAICAgICAgMCAgMFAwMDBQYFBQUFBggGBgYGBggKCAgICAgICgoKCgoKCgoMDAwMDAwODg4ODg8PDw8PDw8PDw//2wBDAQICAgQEBAcEBAcQCwkLEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBD/3QAEABP/2gAMAwEAAhEDEQA/AP38ooooAKKrXt0tlZz3jRvKIEaQpEpeRtozhVHJY9gOpr5Wf9sr4VReO7f4XzaV4jj8X3cH2qLSm0W5F01vzmYJj7gwctnHFAH1lRTUbeivgjcAcHgjPrTqACiiigAooooAKK8w+JfxW0T4WWK6pr+l6te2SxSzzT6dYS3kdvFCMu0xjB2DHPPYH0rnvgt+0D4D+PulN4g+HMOpXGjFS0V/c2Mttaz7W2MsMkgAcqwIIHSgD3CiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKK+bfi7+1R8M/gY80vxKtdZ0zToZo7cagNMnlspJpRlESZAVYnpgdwRXsvgnxhaeOtAh8RWFhf6dBOSEi1G1e0nIGMN5UmGCsDwT1oA62iiigAooooAKKKKACiiigD/0P38oorK13W9L8M6Jf8AiLXLhbXTtLglubmZs7Y4YVLuxxzgKCaANWvzD/Zzz8WP2+Pjz8YXAl0/wTb2fhDT5OoDp89ztPTh42z7MK+tPCv7TXwv8f8AwI1f9oLwffSSeFdMtr+dp7mF7ds2KtvGyQA8lcA96+e/+CZfg7VdG/ZsTx94igMOs/ErVr/xHcFhhmW7k2xE55wUQMvswoA/QyuH8Q/Ez4c+Er8aV4p8U6Xo96UEnkXl7DBLsbo2yRgcHBwcV3FfM/xa1vwdp/ilbfXfg7qvjq5+zof7Qs9Itb2IKScRebNIjZXqRjAzQB6J/wAL0+Cn/Q/aD/4M7b/45WhpXxe+FGu6jBpGi+MtGv766bZFBBf28ssjeiorkk+wFfLv/CU/DT/o2nX/APwnNP8A/j9dZ4F8Q+A7rxbplvpnwG1nwzdPKBHqVxodlbRWxwfneWOVmQdsgE0AfXVFFeKfEb9oL4YfCvx14L+G/jDUJIPEHj+6NppVvFC8xkcFVLOUBCJuYDcfX0BIAPGf+ChHxDm+HH7I/j7UbFymoavaro9rt+8ZdSdbchffYzGvZv2cPh3F8JvgP4D+Hca7G0TSLWKUYxmdkDzH8ZGY18X/ALdSN8VfjZ+z/wDs12WZ01bXz4h1aJedun6WucyAdFbLgZ7gV+gvir4heF/B97YaTqtw8mqapvNpY20bXF3OsWC7JFGC2xcjcxwoyATkigDtqK4HQviZ4Q8R2Or3elXMkk2gZGoWbQyLe2zbd4V7cr5mWXlMKQ38JNcRp37R/wAJtV+FuofGXTtRuLjwrpc1xBdXCWdwZIZLVykweAR+aPLYEMSvGDmgD3Wivkj9or9pzRfhh+z1/wALe8KpdakNfs7eXR54LSSWNft3liGebK4jQCRWw+Mn5QM17TD8V/C48E3/AI6uUv4tN0pd1wX0+5SbCqGZlgMfmOADnKqR19DQB6dRXzvpn7VPwR1bwvYeObLXJG8M6gsZXVTZ3K2EZlICiacx7I+Tg7iAp4Yiu18TfGf4feEfGHhjwNrt/JDqnjKTytK228rwXMm1n2LOqmINtUtgtnHPSgD1OivMNP8AjD4E1P4n3/wdtruVfFWnWR1GW1ktpo1NoHWPzY5WQI673A+Vjz9K8J+H37U2l+NPjt4++HZstTt9P8MjSbS2EmmXAL3N0JnmlciPMcbAxhTJgEKWHBoA+xaKK8d+Ifgr4seJNVt7zwF8R28H2UcOyS2Gk2t/5km4nzPMn+ZeMDaOOM0Act8S9a1ex+Pvwc0ezvJYLHUm177VAjkRz+TYho/MUcNtblc9DzX0XX5s/Ej4bfHqD48fCCzvPjE91fXLa79luv7AsU+zbLIGT92Plk8xfl+b7vUc19OaB8Nfj5p2t2N9rnxjfVtPgmR7iz/sCwg8+NTlo/NQ7k3DjcOR2oA+iq5vxf4w8M+AvDt74s8X6hFpek6em+aeU4VcnCgAcszEgKoBLEgAEmukr4C+Nt8/xQ/bK+EnwLd/N0LwvZ3fjXV4M5SaW3b7Pp6SL0ISYmTB7gHtQB9DWXx+8N/8JJ4b8N+JtI1PwxJ4yeSPRZdThSKO8mjQyeT8rs0UrICyxyBSwB7jFe71+b37XEmp/Ej9pv8AZ6+C3hLMl9outHxhqskZ/wCPOw08eWjuf4fMZnVc9Tivrbw1+0F8MPF3xj8SfAjQNRkuvFvhO1ju9RhELiGJJCgC+cRsLjzFyue/scAHxd+3ep+Jvxp/Z5/ZwiJe317xC2u6lGOQbPSV3Yb2bLjnv9K/TyvzE+HCN8Yv+ClPxA8fx5uNF+EegW/h61kHMa6he5knCnpuUNIrfhX1RoeveErn9p/xDpMHjDWJfEEWgQeZ4bnSVNMht45v+P2DcgR3dm2MysehHbgA+k6K8I0n9pD4Wa/qOu6Lol1eXuq+Grk2t/Yx2Fz9qhcIJMmIxhtmwgh8bTkYJyKvav8AtB/CjRvhUnxrn1hrjwayGQ31tbzThEQsrmRI0LoEZSH3KNpBBxQB7TRXiNt+0T8Jbl9IkbV3t7HX5I4dPv7i2nhsbmaUZjjjuZEERZ+i5YBjwpJrf8RfGHwH4U+IHh74Za9eS2uveKS402M20xinaNWd1EwXy9yqpJUtnHagD0+ivLtG+MfgLXviPrPwm0+8m/4SfQLRb66tZLaaLFsz+WsqO6BJFLcAoTXh3wp/aj0r4i/F34geC/sOqW1n4evtO02yWXTbhCZZYTJNJIwj/doxZdvmY+UbuhoA+wqKzNV1vRtCgjudbv7fT4ZZEhR7iVIlaWQ4VAXIBZjwAOTWnQB//9H92PCfhv8A4RPRIdEGqX2sCEsftOpT/aLltxzhpMLkDoOOlaOtaNpniLR77QNagW60/UoJLa4if7skUqlHU+xUkVp0UAfLmq/sseGfEnhXT/hr4p17UNR8BaYY9mgqILW3njhIaOK6kt40knjBAJVmAYjLZr6ZsbGy0uyt9N06BLW0tY1ihijUIkcaDaqqo4AAGABVqigAooooAKKKKACvLPHHwp0jxh4k0XxxbXLaV4n8PRzw2WoRxRTskNzt82MpMrqQSoIOAQRweTn1OigDxHwB8BvCHgfxnq/xOu57nxF421yNYLrWdRZXnFun3LeBVCxwQr/cjUZPJya8V8RWOp/Df9sO4+MXjFZZPBeveEotFtb5I3mj06+t7szvFKqBmjW4Uhg+NpK7SRxX2xXzt+1B/wAkruf+viL+tAHmvw00vWPGv7Wni343+H4prbwQ3hmy0KOeWN4Rqd9FcvO0sccgVmjgQ7BIQASxC5ANeB+CdSk+H37NHxq+FmuWV0PGNxqvjCS20tIJGmuI9Rlmmt5IiF2NG8bht+7AGc8jFfcH7O3/ACSfR/8Atp/6Ea+XPH3/ACXvxP8A9eM//pG1AHMeL9P1H4j/APBMvTNH8F2dzqOo6d4d0KOW0EEsdx5ultatdRiJ1VyyCNuAOccZr6/034hJ8QfB/jK80eAt4fs9N8uzvWR4zcyPas84CyBTiIlVzj724dq81/Z5/wCTYLn/AK4X/wD6AatfAj/k2fUP+vfUv/RZoA+Zfgn8RvClx/wTy8O/D230y61rxFqnhJ9LttFjs5Wlu7m6ieKLb8vl+WzMGMpYKBkkgjFdJ4q+G/jT4afDf9lSPxiZNQb4aaxp48QXcQacWynTprYyNtDMY0kdUL4IA5PFeSfsg/8AIA+CH/YF07/0E1+i/wC0j/ySDWv+2X/owUAfNFn480XVP2/NK1q2juTo2o+Bp9Lsr/7NN9mubw6jHM0SPswNqqfmbCk5AJxU+m61rfws/ae+Nt7JpU1zqHjWHwzL4eiMUhi1AwwtazhZFG0eSxzICQQvzdK4nwV/yVD4Rf8AXC5/9HPX0D8UP+Tlfh1/1xk/9DNAH1zRRRQB4940+Huq+I/iv8OvHdpPFHZ+EW1U3Mb58yT7da+SmzHHDcnPavYaKKACvzL+B3hGD44ftPfHn43G/utOm0DVLTwlol/aMFkij0yHN4oDhkdJJXGQykHHqAa/TSvz5/4J/f8AHh8bf+yka9/6FHQB9a+AfhD4R+H2qav4lsFm1HxH4gZDqOrX0nn3tyIxhIy+AEiT+CNAqL6Z5rndS+CFtH4i8Q+JvBGvXXhK98WlH1V7KG2eSeSOMRCRJJo3eN9igEqcZGcZ5r3OigDzH4T/AAf8B/BXwwfCvgGwNpbzTSXVzNK7TXN3dTHdJPcTPlpJHPVifYYHFfMh8Qae/wDwUCWVVuTbr4H/ALJM/wBln+z/AG8ag1x5HnbPL3eUd/3sY754r7pr4bi/5Oef/r+P/ogUAcz8BPFfhqz/AGj/ANo/xNeJcW1tqd1pVzbXUtjcotxa6fpyQ3DRO0Q3iOXK4XJJ6A14N4Sv0/4dp/GHQriC5tL+xt/Fs0tvdW01s6xXl5dTwMBKi5V42DDGevPNew+C/wDkNeNf+wXqn/odcZ4r/wCTPPj7/wBi/c/+k8lAHV/GzX9P+MX7IVl8FPh5p9zq3jPxPpek2NnZfZpYmsZk8l2uZ3ZQkKW4Qvv3c4ATdmvR/wBoLQtX0Hxp+zv4v1TzbvSvBeszJrF+iNIIftGmyW6TSBQWCPLgFsYBIziuP/Z6/wCSo6T/ANgqH/0lSvpP9pz/AJJFqP8A12t//RgoA+dPCnjfSNT/AG/9U1eGK6XStZ8DWun6detbTC3u7iK/eWRY5NmBtXPLYBwcZq34O1vWfhb+0D8c7OXTJZtZ8aatoF5oULRSmG+t2tIrWaRZVGwCAq5kywKhc9CM8j8P/wDktXww/wCwTJ/6Mnr3Hx5/ydT4L/68W/m9AGr+1/8ACj4a/Fv4T/2P8SLW5vPsN7DPpMVpO8E76q4aG1RCnJ3NJgg8YyeMZH0rotpcWGj2FjdyedPbQRRyP13OiAM34kZrx/48/wDIJ8If9jVov/pQK9zoA//Z" } }, "cell_type": "markdown", "metadata": {}, "source": [ "![image.jpeg](attachment:image.jpeg)" ] }, { "attachments": { "image.jpeg": { "image/jpeg": "/9j/4AAQSkZJRgABAQAAkACQAAD/4QCARXhpZgAATU0AKgAAAAgABAEaAAUAAAABAAAAPgEbAAUAAAABAAAARgEoAAMAAAABAAIAAIdpAAQAAAABAAAATgAAAAAAAACQAAAAAQAAAJAAAAABAAOgAQADAAAAAQABAACgAgAEAAAAAQAAAtigAwAEAAAAAQAAAFwAAAAA/+0AOFBob3Rvc2hvcCAzLjAAOEJJTQQEAAAAAAAAOEJJTQQlAAAAAAAQ1B2M2Y8AsgTpgAmY7PhCfv/AABEIAFwC2AMBIgACEQEDEQH/xAAfAAABBQEBAQEBAQAAAAAAAAAAAQIDBAUGBwgJCgv/xAC1EAACAQMDAgQDBQUEBAAAAX0BAgMABBEFEiExQQYTUWEHInEUMoGRoQgjQrHBFVLR8CQzYnKCCQoWFxgZGiUmJygpKjQ1Njc4OTpDREVGR0hJSlNUVVZXWFlaY2RlZmdoaWpzdHV2d3h5eoOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4eLj5OXm5+jp6vHy8/T19vf4+fr/xAAfAQADAQEBAQEBAQEBAAAAAAAAAQIDBAUGBwgJCgv/xAC1EQACAQIEBAMEBwUEBAABAncAAQIDEQQFITEGEkFRB2FxEyIygQgUQpGhscEJIzNS8BVictEKFiQ04SXxFxgZGiYnKCkqNTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqCg4SFhoeIiYqSk5SVlpeYmZqio6Slpqeoqaqys7S1tre4ubrCw8TFxsfIycrS09TV1tfY2dri4+Tl5ufo6ery8/T19vf4+fr/2wBDAAICAgICAgMCAgMFAwMDBQYFBQUFBggGBgYGBggKCAgICAgICgoKCgoKCgoMDAwMDAwODg4ODg8PDw8PDw8PDw//2wBDAQICAgQEBAcEBAcQCwkLEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBD/3QAEAC7/2gAMAwEAAhEDEQA/AP38ooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKZJJHDG0srBEQFmYnAAHJJPpWXDr+h3A05oNQt5Bq677IrKp+0oE8zdDg/ONnzZXPHPSgDXooooAKKKo/2np39oNpP2qL7ckQuGg3r5ohLFRIUznaWBG7GMjFAF6ivMviD8TtK8E/Dyf4hWMTa7aZtktxZBrgStdzJBG4MKyHywzhnZVbCgnBxXBWPxg8TaH8bdI+CnxC020W58UaZd6npN/pskjRN9gZBcQTxSgMjBZFZHBKsMjCnigD6KooooAKKKKACiiigAooooAKKKKACiiqVxqOn2t1a2N1cxxXN6XWCN3CvKUUuwRScsVUEnHQc0AXaKpWeo6fqPn/2fcx3P2aVoJfLcP5cqfeRsE4Zc8g8irtABRRRQAUVTn1HT7W6tbG5uY4rm9LrBG7hXlMa72CKTliqgk46DmsbVfFekaZ4f1rxFDJ/aMOhRXL3EdpiaXfaoXkiCLkmXjGzrkgUAdLRXybfftB+J/DE/wANdc8W6Ja/8Iv8UL60020e2eZL3T7rUIjNai4hnRd6tt2SFQrI38JGTX1lQAUUUUAFFFUrzUtP09rZL+5itmvJRBCJHCGWVgWCJkjcxCkgDnANAF2iqVtqOn3lxdWtpcxzzWLiOdEcM0TsocK4BypKsGAPYg1doAKKYsiOWVGDFTg4OcH3p+c9KACiiigAooprsqKXchVUZJPAAFADqKyIdf0K5j06W31G3lTVxusmWVSLkbDJmEg/ONgLfLn5eeleLeIPif44vNb8aaN8PdGtZ5vAscDXK6obi2W+knt/tIS2mWMxhdnymT58PkMoAyQD6AorzX4O/E7R/jP8L/DXxT0G3ltLHxJZx3ccM+PMiLZDRsRwSrAjI4OMjrXpVABRRRQAUUUUAFFFFABRRRQAUUUUAFFYyeItAktba+j1K2a3vJ/s0EglQpLPuKeWjZwz7lYbRzkEdq2aACimPJHGAZGCgnHJxyafmgAooooAKKKKACiiqcuo2EF9Bps1zHHd3SyPDCzgSSLFt3sqk5YLuXcR0yM9aALlFeb+P/iXo/gf4faz8QbeN9bttHDK0ViDcO8yyiApiIO2UkOHwpKgNxxivOF+Mfijwv8AGPwr8JfiFplnjx1aX1xpN9pskhVZtORZZ7e4imUMMxtuSRTgkEFVOKAPo+iiigAooooAKKKKAP/Q/fyiiigAqC6toby2ls7ld8U6NG4yRlWGCMjBHHpU9fOlh+098Mtd8B/Eb4gaNczrpPw1mv7TULm6ge3he409C0ggd8CRd3yAjq3HpQB8w/su+EtH8X/tLfG3xpYvdnwt4K1ez8P6FaG+upbWG8s7cNfzJG8pUsZGAOQQPzr9Ka/Pv9jiW2+Bv7J/hjxV8ShLba/8QNSm1ea3EZa7utS1+5aS3gjj4LSPGU4OMAEsQASPrWw+KWl3PiG78GX2nXeneJbewbUotOmEXm3dqjbGe3dJGifa5VGG8FSy7sAgkA9Por5si/ah8DTfDXwz8VU0rVTonibVotFT9xH51rdTXpsF+0J5vyr5425Utjriq2t/G3xLp/7Sdl8HrXwvqFzpf/CP3GqSXEItyJ5DdQQI67pVZY4gzhicEluFIGaAPpyivENV+Oelab8U7n4PQeHdYv8AxDDpi6vH5EMJt5rNpPJLrM0qqu2T5SH2knpmsLUv2nvh5p/gyDxylrqNzYjW4fD19Eluq3OmajNcJbCK8id1KYlkQEru4YMMqQaAPoyivKovi3o8vxQ1b4SDTL8a3pekLrSsY4/IurV5PJHkP5mS/mArhwvPfHNee/s3/GfxH8YtH1zU9e8O3ukC01vV7KJpxAY447C7a2WEtFIxaUbTuO3bnOGPBIB9L0UUUAFFFfIXxZ8L+ANS8cXl5rvwZ17xheOkQfUrFoBBKAgAVd99A2UHB/djn160AfXtFfn1/wAIR8KP+jc/FX/fdt/8tK+t/hFp2i6X4Kt7TQPCt54Ns1llK6dflDOhLcsdk04w3UfOfoKAPTqKKKAOe8W/8irrP/Xlc/8Aotq+OPAgH2H9lH20aX/0w19keLf+RV1n/ryuf/RbV8ceBP8Ajx/ZR/7A0v8A6YaAPrnXvHvg3wvqVjpHiHWLfT7zUpI4reOZ9pd5n8uNc9FMj/Km7G5vlGTxWrc+I9As9csvDN1qNvDq+oxSz21o8qieaKDb5rpGTuZU3LuIGBkZr4o/aO8O2niL41/CH4UQCQQ+N/ETeI9YlZid9t4UtlngtlP8KGcxttHfcepNWvCsUnxP/bU+JPiSyu/Kg+G3hvT/AAxZ3MYSTyb/AFR2vrtkDhkLoiwqwIOOhFAH3NXzfJ/ydJqx/wCpGtv/AE43FZv7MPibxZ4gX4o2XizWrjXJPD/jbVdLtZrkRh0tbeODy4wIlRAAWJ4Uda0pP+TpNW/7Ea2/9ONxQBg/BvxlB8Pv2RvCfjCbTL7WRpmhwSLZaZbPd3lw5+VY4oowWZmJA6YA5OADXnHwL8cWvxB+Mv8AwmOt+DfEyeL9RsJoZb/WdHudM03RtNjIdbKxadQHeWUqZG4eUgscKqoPev2W/wDk3j4f/wDYJt/5V73QAV5l8ZfHN98Nfhlr/jnTLeO7utJhWSOKYkRsWkVPm24P8Xau+1SLUZ9NuoNIuEtL6SJ1gmkj85I5CPlZowy7gDyV3DPTIr4Z/aV8L/H+2+Bni2bWfH2j3lmttH5kUegPC7jzo+A5vX28/wCyaAPvWivCtN8K/tBw6jazap8QNHurJJUaaKPQHieSMEblVzfNtJHAbacdcGvdDnBxwaAFory34QaF8U/D3hJ9P+MPiS18Va8by6kW8s7UWcYtXcmCLyx/Ei8E/hzjJbq2g/FW4+LmheINI8S2tr4BtbC4i1HR3tQ9zc3rn9zMlx1RUGMr7Hg7uAD1SiuF+Jml+Otb8Ba3pPwy1mDw94purcpp9/cwC4ht5iR87xHhhjI6HGc4OMV0Hhu21uy8PaZZ+JbyPUdXgtYUvLmKPyo57hUAlkWPnarPkhewOKANqvNPhr43vvG3/CU/breO3/sHXr7SYvLJ+eK12bXbJPzHdzjisjxR4c+N19rlzdeEvG2laTpT7fJtrjRXu5Y8KA26YXkQbLZI+QYHHvXzh8EvC/x/m/4T7+zfH2j22zxbqqzb9AeTzJh5e51/01dqnjC849TQB9418+fE7/ktnwa/6/db/wDTXNXsXhWz8TWGh29r4v1ODV9VTd5tzbWxtInyxK7YTJKVwuAfnOTzx0rx74nf8ls+DX/X7rf/AKa5qAIvgLNDbW/xKuLh1iii8Za07uxCqqqYySSeAAOpr1zwx448I+M4rybwvqsGorp0giuRE3zQuyCRQ6nBXcjBlJGGUgjI5rwz4ReH9P8AFWh/E/QNW3mzu/GeriVUYoXRZImKEj+FsbWHdSR3r5F8JfEK48M/sp/Hn9qq5Yx658QtS1ltPjY4ZFtWOi6Vbqp/iBjBwOSWNAH6gaB4h0LxVpFvr/hnUINV0y7DGG5tpFlhkCsVJV1JBwwIOD1FbFfGPgX4I23wx+GHh9PGGsX9n4b+HvhO1t0srDULmyieeCEz391ciBoy7lhhMsQAGPVq4f8AZ6+GXjD4hfsy+FvGmteJ9dh8YeJmj1uO5l1i+H2e0ubz7RDb7PNKMi2hWPDKd3c85oA+g/ib/wAlt+DX/X7rf/prlqr8Dr2HTNM+J2ozq7xWvi/XZXWJGkkKoUYhUQFmYgcKASTwBmrfxN/5LZ8Gv+v3W/8A01y0v7P/APq/iL/2Omtf+hR0AfMuifF2z+Lnxd8M+IfFHw88Ym+07UBFoFjqGg3dnpel+cfLm1K7uZF2NceTu2E/LGDsXLMz1+jdFFABRXxT4/8ACXw4vfGmr3eq/AzxF4iu5Zy0uo2rQCC5bA/eR7tQiOD05RfpXIf8IR8KP+jc/FX/AH3bf/LSgD9Ba+efj3/yF/hP/wBjpZf+kd5XqXw6tNNsPBOk2mkaHceGrOKHEenXe0z2w3H5H2ySjOeeHbr1ry349/8AIX+E/wD2Oll/6R3lAC/CHA+J/wAZieP+J9Y5/wDBTZ14V4N8Z+If2xfiF4oGj6pdaL8F/BN8+k5sJWt7nxJqcH/HwWuIyJI7KHIULGwMpPzNt+WvU/CVlqeo+Iv2gdP0V/L1G61CKK2bONsz6JarGc+zEV5d/wAE0bfS9L/ZH8M+G4F8jWNEu9TtNZtn4nt9TS8lM0cynkPgqef4SD0xQB137QHwQi8O/BfxJqX7PFk/hPx1Y2TDTbnSVKzyuxCNFKoyJg6k8uGKthwQyg19SeC/DqeEfB+ieFY5GmXR7K3tPMclncwRqhdickliMknqTWbqni2RPF2meDtDjS8vJCZ9QJJIs7MI213x0eWTasanqNzdFNeHfEvXPGfh39or4QaTZeI7r+w/FN1rUd1puyFbcraac0sYJVBI2JPm+ZyM9uKAPpfV9Y0vQNPm1bWrqOzs4AC8srbVGTgD3JPAA5J4HNef3Hxt+E1pp41W78U2UFqJHid5JNnkvEVVxMpAMW0soYyBQNwz1FX/ABVfeAT4v8J6N4o1WCDW5Zri50iwlnVGupoIW8yRIjzIYY2LA/w53dcV8L/ALQfhh8UvCvxo+LfxJ06DWfD3xO8bS20NvIhlW5s9Imi0+wUKnLF5oi+B14J4oA/SpWV1DoQysMgjkEGszXP+QLqH/XvL/wCgGvmT9rbWfGHgn4Y2/inwV4hudDa01bQbMwWyQhJY7vVLa3kDO0bSAGJ2XCMor6b1z/kC6h/17y/+gGgD4l+HwH/CMfso+1kf/TBNWl+0d8ZjHrdz8HJfDnjA6JcW6/2vqmgaFeXpnhnHzWdrcRLtRnQ4llBJQHanzklM74ff8ix+yl/15H/0wzV920Aea/B2TTpfhf4aOj+HJPCWmrZxpaaTNEYJbO1T5YI5ImAKOIwpZTypJByRmvSqK8v8baF8WtU1KGfwF4s07QbFYtskN3pTXztJk/MJBcw4GMDbtPrmgCXw743vtZ+JvjHwPNbxx2vhu30qaGVc+ZIb9Z2cPnjC+UMY9TmvS6+D/Bvhf4/v8cviTDbePtHjvI7PQDPKdAdkkVo7rywqfbRt285O47sjpjn688EaZ470vT54fH2u2mvXjy7opbOxNgiR4HylDNPuOcnduHpigDtKK8p+KOg/FfXW8MH4W+JrTw2LHV7e41cXVoLr7ZpiZ862jz/q3fjDdvUV6Dr0Gr3Wh6ja6BdJY6pNbTJaXEieakNwyERyMnG5VbBK9wMUAatFcD8LtJ8faF4A0XSfijrcHiPxVbQ7b/ULaAW0NxLuJDLEOFwuB0GSM4GcVk6doPxWh+L+r+ItS8S2lx8P7jTYIbHRltAt1BfqwMs7XHVlZcgL7jgYyQD1SvMvi145vvh54Si8Q6dbx3U0mp6VYlJSQuy/vYbV2+XByqyEr7gZ4rtfEFtrd5ot5beG72LTdTkjIt7ieE3MUT9maIPGXHtvH1r4k/aG8L/H+3+HlvJq3j7R7uA65oChI9AeJhI2qWwjbcb1uFfDEY+YAjIzmgD7zorxTQvDHx5tdZtLnxB470jUNMjkDXFvDoT28ksY6qspvX2E+u1sele10Afn54cA/wCFLfCP/sob/wDp21CvY/2o/j/d/BHwxoul+ELBNb8f+Or+PRvDenSEiOW8mIBmmI5EEAIaQj2GRnI8c8Of8kW+Ef8A2UN//TrqFY/7Q2lrY/t4fs5+MvFR8vw0ttrenWk0nEEeszwkwozHgPMuBHnksuBzQB9GeDv2cfDcOnxal8Xbh/iL4tuUDXuoaoTJCJWGWS0tc+TbQqeEVEBxyxJ5rkvhT8MPGHhf9o3xzqt7qd/ceB9K0nTrLw3ZXMjvBam8Z575YmblwrxxgFixRSEBCgCvqzVdV03Q9OuNX1i5js7K1QvLNKwVEUdyTXjl9oPin4taf4du9Xubzwzooubm+ntrK5msr24hXKWMMzxlXRXRvNmQMpDBUPGaAPc6o6nqenaLYT6rq1zHZ2dspeWWVgqIo7kmvg34C+FW+KXi34z65d67rsnhnTvEL+H9BgXW9QVIBpdukd3NGwn3EyXLNksWxs4r6Rt7Lwv4Z0H4beB/ih4jS51+NreOzF3dZk1TUrS2O84fDTsvzSjIyGCseQKANyf43fCe106bVbrxRZwW1tI8U5kco0DRqrv5qEB4wqsrMzgAKwJOCDXp0E8F1BHc20iywzKHR0IZWVhkEEcEEcg1+cPwh0z4b/FHV/2kvi78RbODV/C3ibxCnhoROpkW4sdBhitQiqnzN51wTgLy5C9eK+ifiT4TsfC+meL/AIoeKtUvrXw74S0Nzp2l2Go3NjaQ22n27TSSSJbvHuldvkGSQqIoHJNAH01Xz74z/wCTj/hl/wBgfxJ/Owrx34FfCHxlq/wB+GniHxF4o1tvG1yulaxqF1Nq96Q8cs6XU9s8RkaNl8hjCFKYOBk969j8Z/8AJx3wy/7A/iT+dhQBy/wf8UW/gn4LeKvFdzY3mpx6Xr/iaY2un273V3ORqtziOGGMFndjgAD6nAya8p+DnxBtfiV8a9P8Z+IfA/imPxbd21xbQ3Gr6LdadpegaZsMjwW8syhZJ55FRZZOGk6DaihT79+zX/yIOrf9jN4k/wDTrc19A0AFFc5F4v8ACs/iifwRDq9o/iC1tkvJdPWZDdJbu21ZWizuCE8BsYro6ACiiigAooooA//R/fyiiigAryPVvgZ8MddkKazpH2ywa7N+2nSTzNpzXZfzDM1nv8hmL/OcoRu+bGea9cooA8U+NHwon+Jdl4WvdGu47HWfBOt2mu6cZgTbyTWodGhmC8hJI5HXcoJQkMAcYOXp3wy8Sa/8adN+M/jY2thL4e0i60nTdPspnuB/p8sUlxPNM8cWTiFFRAmAMkkkgD3+igD4C1b9m/4t23wX0z4OeGp9Ilj8P+KINfgvLmeZPtcEOsHVFgeNYW8l8HYzguMjgc8e1eOfhd4/1L4s6d8VfB15Z2t2/hq78O3KzO4a1a5uIrlLqAiNhKY2RhsYJuyDkcivpOigD4Y1fX7+2/bsFl4am07UL5fh8kEttdXf2eRn/tNpAqsiyEOF+cqUyUORgddDUP2W9f1H4Z+M9MuNYtm8X+NfFVp4uuXVXWwS5sbm2lhtE6v5QitUjMhXcSS+3+Gvsz7BY/aftn2eP7R/z02Df0x97GenFW6APmTRvh58UD+0kfjLrEWlw6JeeGl0KS0iuZZLm3eO6NyJAxhVJQ5OCPk2jpu5qx8GvhT46+GWpano9zqNtJ4bk1vXNYiMLOZ7n+2Lk3McU0bIFT7OzuNyu2/5TgcivpOigAooooAKKKKACiiigAooooAwfFSPL4Y1iONS7vZ3AAAySTG2ABXyL4I0nVYrL9l0S2UyHT9IlW53RsPIY6HsxLkfId3y4bHPHWvtiigD5U+JNl4v1X4pSv4k8HL4h8K6Npsdx4efTlkTVP7ZlLx3AN4skaW0RiwjbiAwbJJwRWh+zf8ABXVPgl8PdbFzFaTeLvFup3uvakkUjrapeXZHl2ySlWcxwRqkfmFSWwXxzivpuigD5p+AXw4+JPw31Tx0/jCPSpLTxh4iv/ECPY3M8kkDXgiUQMksEYbb5ZJcMOv3a0JLG9P7TOq34t5PsreCreIS7D5ZkGoXDFA2MbsEHHXFfQ1FAHhn7M1nd6f8AfAdlfwPbXEOlQK8cqlHU46MrYIP1r3OiigArh/iT4GsfiV4H1bwNqVxJaW2rRrG8sQBdArq+QG46rXcUUAFFFFABRRRQAUUUUAFcT4K8EWXgn+3vsVxJcf2/q11q0nmADZJdbdyLj+EbeM8121FABXg3xIsryf4yfCG7ggkkgtb3WTNIqkpGH0yVVLkcLk8DPU8V7zRQB8teA4/HWj+CPi/deDNPjl8SnxHr8+kw3oaOCecqhgLHgmNmA5BwfWvEr79nq4+LcngPQNU8Lr4f0nSdStfEnii42S2tre6pAfOa20+wkkby0nuT5k0pVflGAWZzj9EqKAPnz9pfwL8TPin8J/Efww+HUunWEniiwlsZr6+nmjMEc2Fk2RxQvuLRllzuXBOea9W8B6Te+H/AAdpHh++s7ewbS7aK0SG1maeFIoFEcYV3jiJ+VRn5Rjpz1rraKAPBfiRZ3k/xl+EN1BBJJBa3msmWRVJSMPpsqqXI4XJ4Gep4pfgPZ3lmnxA+2QSQef4w1iWPzFK742ZNrrnGVPYjg17zRQAUUUUAFFFFABXgXxysry81b4WtaQSTi38YWcspRSwjjFpdgu2OigkDJ45Fe+0UAeCfCiyvLb4lfF64uYJIorrXLJ4ndSqyKNLtFJQkYYAggkdxitvV/gR8NNW8UXfjSLT59J1zUdv2y60u8udOe72jCm4+yyRiUgcBnBOOM17BRQBz/hzwtoHhKxOneHrJLSFmLuQSzyOerySMS7se7MSfevBPij8PPil4q+L/wAPvH3hy30f+zPAc+oTCO6vLiOe7/tC0+zEYS3dY9hJPVt3tX03RQBwXij4aeCPHN9pms+LdHhvNS0pJo7afLLLClyoWeNZEKtskAAdTwwHIqvofwl+HXhnVDq+gaHDYy+c9yscRZbeOeUYeWODPlJIw6sqAnnnk59FooA+cv2l/hx4/wDi14Gi8EeChpsKvf6Zfy3GoTzR7Tpt9DdiNUihk3eYItu4sNuc4Ne0u2sXHha4Os28VvqD20okitpGnjDbWACOyRs2Rjqg5rpKKAPiLwFpGqw+Gv2XUmsp0bT7Mi5DRsDAf7ClTEoI+T5vl+bHPHWvt2iigAooooA4nR/BFlo3jfxH45huJJLnxJDYQyxMBsjFgsqoVI5+bzTnPoMV21FFABRRRQAUUUUAFcR8QfA9j8Q/D8fh7ULiS1ijvrC+DxAFt9hdR3SLzxhmjCn2PFdvRQAUUUUAfB3h/RtXT4O/Cq1exnWa28fvNKhiYNHF/al+29hjIXaQdx4wQe9fZHjPwR4R+Ifh+48LeN9Jt9a0q5Kl7e5QOu5DlXXurqeVZSGB5BBrqaKAPKdB+C/gLw/PbXEUF5qJsiDbLqWoXeoRwEdDGlzLIqkdjjI7EV2vieXxTFpMh8HW1pdak3CC9neCFQQfmLRxSscHHy7RkdxXQ0UAfOH7LHwq8XfBX4S2Pw98ZtY3Wo2k1zdT31jNLIL25vriS5nldZYoyh3yYAy3HcV6z4u+HPgrx3caZeeK9KjvrrRnkkspyWSa3aZPLkMciFXXeh2sAeRwa7aigDzbR/hB8NvD+qnV9D0GCxlMy3Plw7ktxcKgjE4twfKEuxQPM2buBzmvMP2pfhr8S/jN8MNW+FXgeTTbGz8QLBFe3d7cTxyfZ1nWSeFI4oX/ANbGpTcXGNx+U45+mKKAMfQIr6DRbO31G1gsbiGMI0FtI00MYThVR2SMkbQOqD0rx/xfZXkv7Qnw4vooJHtoNJ8RLJKFJRGc2O0M3QFsHAPXBxXvFFAHg/7OtleWPgXVIb6CS3kbxH4icLIpQlH1S4ZWAIHDAgg9xyK94oooA8Zk+GWgTfHOD4p22gW9pqtrpL2VxqwVRcXayuvl25IOTHEELHcOrLjocezUUUAFFFFABRRRQB//0v38ooooAKKKKACvD/jR8f8AwX8BNGk8TfEGz1OPQ4FiM1/a2bXNvE00nlRo5Q7tzMQAAp6ivcK/Pz9t8nx14s+Bv7PkJ3r418XQahfxjndpuhL9rnDD+6W2de4oA+4/C/iBPFOh22ux2F5pqXQLLBfQ+RcKM4BeMklc9RnnFdBXK6X428L6z4p1vwVpl8s+teHFtnv7cI4MC3il4SWKhW3qpPyk474rqqACiiigAooqtbXtpe+b9knSbyJGik2MG2SJ95Wx0Ydx1FAFmiiigAooryHxZ8Wbnwtrk2ix+BfEutrCqN9q06yimtn3jOFdpkJI6H5etAHr1FfPv/C+73/ol/jP/wAFsH/yTXrfhHxJJ4s0WPWZdHv9CaRnX7LqUSw3K7TjJRHcYPUfN0oA6eiiigChqt8NL0u81Nk8wWkMkxUHG7y1LYz74rx/Q/jPDrUHwrmXS2i/4WbZtdoDKD9jC2P23a3y/P8A3MjHr7V6f4t/5FXWf+vK5/8ARbV8ceBP+PH9lH/sDS/+mGgD7nor5l+MPx81P4Z6rpUdhoy6hZ3fiDSPDxV2dLi6u9VddwtABtYW8LebIzcHBUYIJrrdQ+MPkftBaP8AAnTtLN3JdaDda9fXvm7RZwxTpb26bNp3GZ2f+IYC5wewB7bXAfEDx3H4Ei8PyyWZvP7d1mx0gAOE8s3rlBIcg5C4zjv61ZtfiT4Cv/GJ+H9hr1pdeI1glumsYpRJMkMDrHIzBc7druqkEg5PSvMf2hf+PX4e/wDY66D/AOjjQB3el+NrjxN4n8b+CNNg+xXnhcWkUd058xJJL2289G2YGAhIBGTmvkDx/wCFfGPgq6+D2heFdUk1H4xXOr2EmtXVlJOLa80mMn+1ri9ikd1FsynEe/kSlViweB9A/DwTN8afjKtuypKZtECMw3KG/s4YJAIyAeoyPrXl/hn9nv8AaG0vxHJquu/GGzv7TU9Qgu9VWDw7HbXl5BDKJBbLd/aXeKMKPLQAEKpOBkkkA+16KK5248WeFoTJDLrNlHImVZWuYwVYcEEFuCKAOiBBGRzRXzN+zL4w8Oj4C+C/t+uWv2j7D8/m3Me/PmN13Nn86+lYpYp4kmgcSRyAMrKcqwPIII6g0ASUVgQeKvDNz4jufB9vq1rLrtnAlzPYLMhuooJDhJHiB3qjHgEjBpfEfinw14P0w6z4s1W10awEkcX2i8mSCLzJWCou9yBlmOAM8mgDeopAQRkcg1g+H/FXhnxXFdz+GNWtdWjsLiS0uGtJkmEVxCcSROUJ2uvdTyKAN+jIzjvWdqGr6TpQQ6pewWYkzt86RY92OuNxGcV87Xvi/wAPf8NM6ORrlp9l/wCER1DP+lR+X5n2+0x/FjdjOO+M0AfTVcJ8T/G8fw2+HniLx9LaG+TQLKa8MCtsMohUttDEHGfXBrqNP1rRtVZ00u/t7xo8FhDKkhUHpnaTivF/2pf+Tc/iN/2A73/0WaAO4/4T6P8A4WbZfDj7Gd95os2sC438KIriKDytuOp8zOc9sYr0OvnL/m6DQ/8AsSbz/wBOFrVGD4+6jL8afCXw0/sdDY+MrfWbi2fc63dvb6QUX7VPGwCiG4kYpGB833STyQAD6borwDwl8edL8S+OPivoFzZjT9B+FclrBc6q8u5JppLX7XdLsC/KLdCoJySSTwMUyb9p34U29x4fs7mTVobjxYzJpMT6LqKSXzJH5zCFWtwW2x/OeOByaAPoKiuF8GfEbwx49m1a10B7lbnQ5kt72C7tJ7OWGSSMSoClwiEgowIIyPeu6oAKKytZ1zRfDtg+qa/fQadZxkBpriRYowT0G5iBk9h3ql4T8W+G/HXh2x8W+ENQi1XR9SQyW11CSY5UBK7lJxkZBoA8q1/43QaHofxR1o6S8w+GkgR0EoH2vNlDeZB2nZ/rtvOeme+K6Tx/4V8SeOH8LR6VqMVlocN8LjWrSQSbr6z8iQLAjxMhUiVkc5ypC4II4r5a+IH/ACIv7VH/AF8L/wCmSyr6L+L3gX4gfEb4fW/hb4feLYvB09yYvtd1JY/bjLbBfngC+dDsEhwGYHO3IGM5oA8l/Z4Ost8X/ionhOe4m+E9vJYQaR58rzQf2rGjjUhYNIWP2ZW2KcHy/ND7O9fY9eB/Br4ZfEzwLc3tx8RPG1v4pj+zW9np1tY6VHpFrZQxF2YCGOWRWZ8qAeNoXA6175QAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAH/0/38ooooAKKKKAPAPjbq3xg0jWfAc3w4mt7Tw/8A2sG8TXElnJfTDTkTd5UMUQZw0pyvmAHYdpPGa8o8I+BPF3xQ/atl/aG8SaXPovhTwfokmheGre9Tyru7nvJA95fNCfmijKgRRh9rsMsVAxX2tRQB8yeGPEXjTVP2gfin8ObjXHGl2Oj6Le6cy29uJbKW/N2km1tn7wDyVZRLuwcjpXgHgj40fGH/AIUT8Ivi94g8QLqOo+KfEmmaLf2v2WGO3mtL2/ksWf5UDrNgLIGVgoI27cV9oQ/CLwFbeLdc8dW9jNFrviS2Wzv7pLy6VpoEzsTAl2oE3Ns2BSuTtxk1yi/s1fBiPwjpHgKPQZE0DQb1dQsLRb68VLe7R/MSVCJtwKONy84ViWABJNAHA6prnxJ1H9qu7+F9j4um07wzL4Sj1sQR2ls00Vz9tNriKZ42IUqNzbw3PAwK8V8Q/Hr4sw/DfxPcQawkOt/Dv4had4Xnv4LWEwavZXV5ZxsXjdWEcixXWyTyyMSISMA7R7frHwX8S6/+03H8UNStIU8NW/hkaDDNBqdxBqKSi5NyZcRIvyMD5ZUykn7x9K9o/wCFR/Df/hBpfhsug26+HJmMj2q7gGmMnnecXB8wy+aBJ5m7fv8Am3Z5oA8qtvF3jN/2qdc+GDa676BN4Pi1eC38mDfZ3b3jW29HCbmBVdwEhYZz24rzr9jV9Us/DWuf8JJ4uu9Wur3xX4qt4oLxYN0z2mpSq025I1csEUZGdoBwFAAx9FaV8FPh5o3jmL4l2dlOfE8VobE38l5cyyy2xbd5c2+QiRQeVDg7e2K1vD3wt8C+FfEF/wCJtA0tbS+1Ga4uZSskhiE92Va4kjiZjHG8zIpkZFBYjJ70AegUUUUAFeQeLfh7491/XJtT0L4l6r4cs5FQLZWtnp0sSFRgkPcW0kh3Hk5Y+3Fev0UAfPn/AAqP4q/9Fn13/wAF+j//ACHXrnhDRNZ8P6JHpmva/ceJbxGdmvbqKCGVgxyFKW8cceF6DC59a6eigAooooA5/wAW/wDIq6z/ANeVz/6LavjjwJ/x4/so/wDYGl/9MNfaus2LappF9piOI2u4JYQxGQpkUrn8M14j4e+DN/otv8H4ZNSilPwzsXtJyI2H2otp32LcnPyjd83OeOKAOL+PHhHX9S+NnwX8eSaZc6v4T8HXesXF/FaRGeWG+uLPybG4MK/OyoTIuVB2lgSMcjh/h5pHxQb9pr4h/EbVfC11po8Uf8I/YaVcXKK8MOgWEcst0XZWISeSZmHlE7lLKSCAcfdlFAHyJrUMcX7cfhLykCBvAmtE7RjJOpWWSfeu3/aF/wCPX4e/9jroP/o411N58Dvhtf8AjeP4kXVhct4lhR4o70aheq6QyOJHiQCYKsTOoJjC7TgcVpfEnwLc+O4fDcVtdpaHQtd0/V33qW8xLJy7RjHQtng9qAOG+Gv/ACXD4w/9fGh/+m9a+gq838LeBrnw94+8beMZbpJovFclg8cQUhoRZ2wgIY9DuIyMdBXpFABXk9/8CvgpqFxc6hfeAdBuLm5Z5ZZZNNtmd5HJLMzFMkknJJ5Jr1iigD4u/Zp+B3wY1X4EeDdQ1TwHoV1dT2W6SWbTbZ5HPmNyzMhJP1r7HsrKz02zg07ToEtbW1RYooolCRxxoMKqqMAKAMADgCpYIILaJYLaNYokGFVAFUD2A4FS0AcVafDjwLYePL/4n2eh20PivU7SOxudSVMXEttEQUjZu6ggfkPQU7x58O/A/wAUNAPhX4haLba/pBmiuDbXSb4/NgbfG+PVSMj/AArs6KAGqqooRAFVRgAdABXGeCfhx4F+G9tqVp4E0O10OHWL2bUbxLVNgnvLjHmTP6s2Bmu1ooA5LxV4B8DeOUt4/Gvh7T9eWzLGEX9rFciIvgMU8xW25wM464r5avfgd8GV/aV0jS18B6GLJ/CeoTND/Ztv5RlW/tVVyuzG4KSAeuCa+0qiMEBnFyY185VKB8DcFJBK564JAOKAOS8K/DrwB4GkuJvBfhvTdBkuwqzNY2kVsZAmSocxqu4DJxnpXm37Uv8Aybn8Rv8AsB3v/os173Xnfxb8E3HxI+GXifwFaXS2U2v2E9mk7qXWNplKhiowSBQB5t/zdBof/Yk3n/pwta4ebQPEvhX9rbxR8Vtb0O/1nTNR8K6bpOhS2UPniKSK5mlvbd+QImkcxOGfCED73ymveP8AhAbn/hbFh8RvtaeRZ6DPo5t9p3s8tzDOJN2cYAjIx1ya9NoA+Dv2YPh745tvD3iTw18UfDklhc+J9f1/W/Ef2gZjnnurxfsMNvMrATRC2jUsw+XaAp5JFbt8v/Cwf27dNtsl7D4S+Epbpv7q6l4im8pB9Ra27n1w3vX2fPClxBJbyFgkqlTtYq2GGDhlIIPoQcjtXlPhX4G/DTwV4rv/ABz4b064tdd1Yxm9umv7yZ7ryUMcfniWZlk2KxCbgdvagDrtUgXwvbax4k0DRZtW1O+eOWa3t5I0muXREhXDTuka7UUdWA49Tz57/wALQ+In/RJ9b/8AAzSv/kyvb6KAOdsbq51vw7FfavpUmmXE0Zd7O5MUkkLjOAxiaSMnvlWNfOf7Df8Ayad8Nv8AsHN/6Pkr6V1/QNO8TaXLo+q+d9mmxvEE8ts5x28yFkcD1APPeuY+HXwu8EfCfRF8N+ALB9L0qIBY7X7TPPDEoJOI1mkcIMsSQuM96APkf4gf8iL+1R/18L/6ZLKvu3T/APjwtv8Arkn/AKCK8B8R/BLUNc8P/F3RotUiif4lSB4HMbEWuLCCzw4z83zQluMcECvoO3i8i3igzny1Vc+uBigCavKfhj4y8f8Ai268Vw+OfBcvhCLRtWnstMkkuY7j+0rKMDy7xQn+rD5+6eR+derUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQB//U/fyiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigD/9k=" } }, "cell_type": "markdown", "metadata": {}, "source": [ "![image.jpeg](attachment:image.jpeg)" ] }, { "attachments": { "image.jpeg": { "image/jpeg": "/9j/4AAQSkZJRgABAQAAkACQAAD/4QCARXhpZgAATU0AKgAAAAgABAEaAAUAAAABAAAAPgEbAAUAAAABAAAARgEoAAMAAAABAAIAAIdpAAQAAAABAAAATgAAAAAAAACQAAAAAQAAAJAAAAABAAOgAQADAAAAAQABAACgAgAEAAAAAQAAAl6gAwAEAAAAAQAAAGoAAAAA/+0AOFBob3Rvc2hvcCAzLjAAOEJJTQQEAAAAAAAAOEJJTQQlAAAAAAAQ1B2M2Y8AsgTpgAmY7PhCfv/AABEIAGoCXgMBIgACEQEDEQH/xAAfAAABBQEBAQEBAQAAAAAAAAAAAQIDBAUGBwgJCgv/xAC1EAACAQMDAgQDBQUEBAAAAX0BAgMABBEFEiExQQYTUWEHInEUMoGRoQgjQrHBFVLR8CQzYnKCCQoWFxgZGiUmJygpKjQ1Njc4OTpDREVGR0hJSlNUVVZXWFlaY2RlZmdoaWpzdHV2d3h5eoOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4eLj5OXm5+jp6vHy8/T19vf4+fr/xAAfAQADAQEBAQEBAQEBAAAAAAAAAQIDBAUGBwgJCgv/xAC1EQACAQIEBAMEBwUEBAABAncAAQIDEQQFITEGEkFRB2FxEyIygQgUQpGhscEJIzNS8BVictEKFiQ04SXxFxgZGiYnKCkqNTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqCg4SFhoeIiYqSk5SVlpeYmZqio6Slpqeoqaqys7S1tre4ubrCw8TFxsfIycrS09TV1tfY2dri4+Tl5ufo6ery8/T19vf4+fr/2wBDAAICAgICAgMCAgMFAwMDBQYFBQUFBggGBgYGBggKCAgICAgICgoKCgoKCgoMDAwMDAwODg4ODg8PDw8PDw8PDw//2wBDAQICAgQEBAcEBAcQCwkLEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBD/3QAEACb/2gAMAwEAAhEDEQA/AP38ooooAKKKKACiiigAoorzvwl8UPC3jTxZ4t8F6J9q/tLwVcQW2oedbSwxeZcx+ankyOoWUbepQnB+ooA9EoorzvX/AIoeFvDfj/wz8NNT+1f214tju5LHy7WWS322SB5fNnVSkZwRtDEZ6CgD0SiiigAorzzxb8T/AAt4K8U+E/B2t/av7R8aXM1rp/k20s0XmQR+a/myIpWIbehcgH869DoAKKKKACiiigAooooAKK8C+KH7QGkfCjxR4c8Ka14a1e/ufF18um6S9kltIl1dtGZSgDzo6BVUlndQox1r3iCR5YI5ZI2hZ1DFGxuUkZKnBIyOhwSKAJaKKKACiiigAooooAKKCccmobe4t7uFLm1lWaGQbldGDKwPcEcEUATUUUUAFFFfP15+0Fplj8YLH4IzeFtZPiLUbWXUISq2rW50+GYQvdNILjKIHIADKHPZSeKAPoGiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAP/0P38ooooAKKKKAMnXrbU7zRry10W/wD7MvpI2EN0Yln8l+zeWxAbHoSM1+ddh+0H8c7P9juD9pnVdRsrzVdLlxJpq2iR2+oww6mbGQyOCXheRclNhwhA3B+a/RbW49Yl0m6i0CWCDUXjIge5RpIVc9C6IyMw9gw+tfED/so/EiT9luT9mc+MdLEMkzH+0/7Mm3+Q94b4r5X2rG/zTt3bsbO2eaAPUNX+IHxD8E/Gb4YeEdc1C31TTfiKmqw3FutusX2G5sbUXcbW8incyEBo2Em4nhgV5FZ/gPxj8UPEPxd+Nvwu1DxFD5fhZdGl0a7SwjD2o1O3lmZZE3bZthUAFiMgVoeMvg98U/F3xF+GfxBfxFo9tJ8PjdSPbjT7hlu5b63NtPhvtIMaiM5jHzEN94kcU/wv8H/ih4Y+K3xI+J8HiHSJ5PHsFjHHbNYThbR9MiaG3JYXOZAUcmQYXLY2lRxQB89aR+0D8bof2PtK/aZ1vU7KbUrR7dp9PSzRbe+g/tEWUrOwO+J3Qlk8s4QgAh+a+g/if4+8eeE/2gvhB4a0nVIf+EZ8czapa31k9qjSBrKye5jkjuM7lywAK4xgda8ub9lD4jt+yxH+zJ/wmOl+RG6L/af9mzeYbdLoXgHlfasb/NG3duxs7Z5r1Dxt8Hfif4z+Ifwx8fy+ItItZPh5JczvAthOy3ct7bm2nwxuQY1EZygwxDdSRxQBzOjeNvjBq37Unjj4LP4jgTw9puh6drdrcJYxC6tftk88Rt1LFkk/1YO91OAMbcnI89sPjp8adV/Zw+IPxkk1CxstT+HFzrtrHEloHh1Q6DPIjSThm3RCdY9uyJgVOWDEEKPbNL+DfxG0r9oHxT8b4fEGlyReIdIt9HSwaxnBhhsXmlt5DKLj5nLynzPlAIHy4PNed6b+zD8TLL4D+P8A4Jy+LtJkHjq71a4a+GmzqbZdbkllulEf2khihkxESwwB8wNAHSfEj4r+PNK8SfAfWvDl/DaaJ8RNVtrDUtOktklYx3NhNdh45yQyMpjC4AwRzxXO/Fz4l/GzwMj6Ro+q295458aazdaf4Q0COxiKC1ifAu72UuzeRDEDPNINuAVjADGug8T/AAD+KXiK0+EsP/CUaPBJ8Lry3vgf7OuGF7LbW8loqkfah5amKQk4LHfyDjirEHwP+LGnfFLxl8XLXxLo1/4h8QRCw0q4v7C4kGjaWnKWsESXCq26T97M+QZHxnAAFAFLXvjR4hh+MukfAGTWxY3WmeG4tc17WbSxM0k0ssv2eKG1gKTJEsjq8jM4fChVXklq9W+Bvjjxd4w07xHYeMrOVbnw7q89hbag1rJZx6rZhElgu445AMEq/lybfl8xGK4BAHHaf+z9q3hXXvCPxC8I63CPGOhaDH4e1OW6hY2mr2asJSZERg8UiThpI3UtgMysGBGPevDGl67p8Fzc+JNRGoX97L5riJSltAoUKsUCMWIUAZJJJZiSccAAHT0UUUAFFFeMfFP4Y6l458QeC/FGmanNby+DL6W/Wx+1TWtreSPHsT7QYMs4iPIRlZWyQRzQB8/eKd3xD/b68HaACJtP+FXhW+1qZeqpqGtSizgz/tCBJCv1NfcwdC5QMCy4yM8jPrXzz8PfgzrPga7+IHxCvNVi1b4h+PmEs935RitLZbWExWNpBGSzCCAHksSzsWY4zgeB+GtM8Z/8Ij8K9O0TTb+1+K1jqmmv4qvLiCZcwKT/AGu11dMBFPFMN3kqrtktGYwNvygH6Bb1xnIrAfxToaeKY/BZuV/teWze/EA5It0kWIufTLuAM9ecdDX5vfFfw5qugfDj48+E49C1OXUNQ8UWus6FBZ2tzOTbFLAvcWzRKwAV0nLKpzncNuSAfbPFXhzR7H9qDTvimvhGTVotb8LRW9pcx2kmZNRjv43RZXKfuXSBt4Mu3Cqw6gigD7U3qOpHJx+PpS7lzjPNfnZc6b4H/wCFuftC6ZrekX9/cv8A2S+nxWMN1LILq40tXJtGhysNy8uG3goc4YtgEje8IeHvG/hXxz8PvEvxktLi8vdR+H0mj6/cwxS3Udxq6y2reVMIAwMhTzAGIAb5gD2oA+1PEHirQvC6afJrd0tuNUvILC2B5MlzcttjRR3z1PoASeBXRdelfllovgR9e/Zu+B154n8L3lxceCPFsL6pDeWU73llY/aLpZSY2UyNGFeLcUDArggkA1+lfhK4W78NabcJpraRG8K7LRxtaFAMIpXjB244xkdDQB5x4o+OngHw3qGoeH9UTV/tdpujkMGh6ncxZ25+WaG2eNxz1ViPevEv2efj78PLD4NeCNGnTWTcR6fbxHZoOqvHuPHEi2pQjP8AEGx719ozqzwSIn3mUgfUivNfgt4a1jwb8KPCvhbxBGIdR0uwhgnRWDhZEHIDDg/UUAeng5Ga+StV+L/ivxbpvxO8V+CtSt9C8NfDX7fZpdzQLcHUdS02DzbotuICWsL4i+XDu4chlCjP0T498WWHgLwP4g8b6owS00DT7q/lJ6bLaJpD+e3FfEXwW+APjbxZ+y54b+G3jbW4rXQfGMH9sa+sMLJqNwNWkN9dWRl3bEV3lMbyBS5jyoAPzUAfYfwg8dS/E/4VeEPiNPZHTpfE2lWeovbHJ8lrmFZCgJ5IBbAJ6ivkH4b6vc+If2mv2gfjfY6Vc6+vgm10/wAH6ZaWmwz3EllGb29igMjIm5ppUX5mAyOTXvuv/B/xBafEm0+JPgHUkt5dP0A6FZabdzXCaXZjzN4nW1gISVsbV2ttICgBh2818Q/A7SfhL+yb8Q/A2nteeJdY1uz1fUbu6WCSS71HXNQVn85IoAzKxm2CNVzsAXngmgD7Itrgz2sNxNG1u0iK7RyY3ISMlWwSMjocEip9ykAgjDdPevhXxRpdvrnjv9ni9vdL1KaGKyu4tVcWd5sWGTT9kUV6RHhVa5A+SXHzDJGM14j4oEw+C88OvWWqW9vpHxZthpsctveW06aRdavH5S2oZUkMLwsyqqcbTtx0FAH6sBlK7gQR69qXcpwMjmviW38Hal4i8cfEey8D6XPD8ONW8KC0+yzxyW1tc+IS8uJLSGYIVAg2LK6hVZ9vJZWI8+8CNfp46/ZsN5pWqrd6P4av9H124ksrsJDd/YbdI7e6cpjd5qyYZvlByd3zDIB9+6D4n0TxM2ojRLlbr+yruSxuCvIW4hCl0z327gD6HI6it+vjX9mzwnp/w88V+OPC9v4SbS7y58Sa3efaVtmhhTTp5IpbTZIVCSLNuO0IxIKtuAIr7KoAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooA/9H9/KKKKACiiigAoriviN480T4YeBdc+IPiPf8A2boNq91P5YBbYg7ZIA+pIA6k4rlNM+Iniq+1nTNP/wCEbgurHWbG7vLbULDUFurUPbiMxwSP5SgNL5h2uu5MKeaAPYK43U/iH4F0bxXpvgXVtesrTxFrCs1lp8syLc3CqCWMcZO5gApJwOxr580/9p7U9R+Cfjj4uw+CJluvAGoatY6hpRv4d5GjZ+0SRz7djcAlVxz61558bvGX/CQeLP2XPH2maTLcS63rUt3b2isnmr9t0W4cIZGKqAu4F26YGQCcCgD74qhHqulzajNo8V5C9/bosstusimZI3yFZkB3BWIOCRg44ryzRvi5A3xKHwj8YacdC8SXdi+p6eBMtxbX9pE4SYwS7UbzIWZfMjZAQGDDcuSPMtJ8WeB9D/ad8a2es+DIfD/iBPDEOq3XiU3KSC90m1naFVkQAGLyypY57AcnjAB9X0V8vWv7SkR8KeGvijqnhyWw8A+LLy0tLPUXnU3UceoyCKzuri12ARwTuyAESs6h1LKOdutf/H1rvUfG8Xgbw++vaZ8OTJFrN6ZjEpu4YRPLaWcccUz3E0aEbxhFDkICWzgA+i6K+HPiD8aviFffFz4F2Xg3w076B4xg1HVQl1d/YbuZo9MaRbeaB0PleUJdzbiSXAAAxmvtmxlup7KCa+gFrcSIrSRBxII3I+Zd4ADYPGcc0AWqK+NPij8T/ifpH7Unw8+HugeH4r3R7zSta1AH7eIHu5bdIY8OpQqiRebkZyWJyMbefatb+Kz2/jaz+F3hrTE1XxfNpw1W6tnuRDbWFmX8pXuJwkjZkkysSrGxfazcKM0Aew0V5d8Pvihp/jfVfEXhO6tW0nxP4Snig1TT3kWUxi5jEtvNFIuBJDNGco+FOQysoZSK9RoAKKK5DxL4nuNIu7LRtHsDqmr6kJHhgL+TEsUOPMlll2tsRSyjhWYswAHUgA6+ivmH4e/tJSfEifw1caD4J1ZtG8QX2o6ZLfhVeKyutL8xLh5ivy/Z/NjMUcu4GRvupjJr6cDoXMYYFlAJGeQD04oAdRSZGdueT2qC6u7Wxt5Lu9mS3giBZ5JGCIqjkkscAAUAcV4c+GvhXwr4s8Q+NtHjuE1bxS0T6jJJdTzJM0ChIiIpHZE2INq7FGF4rvaht7i3u7eK7tJVmgmVXjkQhkdGGQykcEEcgjrXJeM/HWleBba2u9VstSvUunKKNO0+51B1IGcuttHIUHoWwCeKAOzorwr/AIaE8If9ALxR/wCE3qn/AMj16b4R8W6f4z0o6xplrfWkIkaPZqFnPYzZXGT5VwiPt54bGD2oA6iue8QeLPDXhU6aviPUodPbWLyHT7MTNtNxdz58uGMfxO2DgD0NdCTjk1+Wv7YV9Lr1r8PPjqJydL8K/EHw/BoyKxCvbtcmC6u2HQ+dJhYyR/q0DD/WGgD9P76wsdUs5tO1O3ju7S5QxywzIJI5EbgqysCGBHUEYqyiJEixRKERAAqgYAA6ACh3SNd0jBV4GScDngVwvivxjeaLdjSNC0z+1tSFtJeSI8wtoYrePjdJKVbBdvlRQpJIJOACaAO8or5u+Fv7Q/8Awti58KX2g+DtWh8OeMNOuNQtdVlQeTALYqhS5AysZkYkQ4di4UtgKQa+jw6MzIrAsvUA8jPrQA6uB+IPw08J/E/T7HS/F8VxPbadeQX8CwXU9qVubZg8MhMDoSY3AZckgEZxXdySxRDdK4QHuxx/OvGh4mvP+F+Hw9/aP/Eq/wCEZF19n3jy/tH20p5mP72zj6UAezKoRQoJOBjk5P51ieI/E/h3wfpUuu+KtTt9I06EgPcXUixRKWOBl2IAzW1HLFKN0ThwO6nP8q+Wf22fHN58P/2X/HuqaVk6rqVj/ZFgqjc7XeqOtnEFHdsy5A9qAPTLT4//AAPvkiltPH2hyxzkCOQahBsYngYfftOTwOa9bR0kRZI2DKwBBByCD0INfn78QvEvh74Yfsx6b8D9I8LzgXei2/hm1n1uzew0O3nmhFsst9dXChEXedwGCzvhVGTkfVfwN+H118Ivg94N+Geqaw2t3fh7Tbexe9kJzPJEnJXcSdo6IDyFAoA9Zorwzx58cdM8CeLNA8NXGlzXkWua1Z6D58bqCl7fRvMgSIjMixxrvmYEbFIPzHIHUaB8T9I8SfE3xX8MdOtZmuvB9tp897dHb5HmaiJHjhXndvVIw7ZGMMKAPS6Karo4JRgwBIODnBHUU6gAooooAKK43x/4ovPBXhLUvFVpYLqK6VDJczxvcLbKsEKF5HMjKw+VVzjHNfPH7FLeOLj4HaXrvjmylhvfFMlx4g8+e8+0vKNYnku0QIRmJYonRAmSBjigD65ori/GfjvSvAtvbXOq2WpXq3TsijTtPudQZSozl1to5Cg9C2Aa8+/4aE8If9ALxR/4Teqf/I9AHutRTzw20MlxcyLFFEpZ3chVVVGSSTwAB1Nc94S8V2HjLSRrOm217aQmRo9l/Zz2M2Uxk+VcIj7eeGxg9q+a/wBuXwvoetfsv/EvWNTtzNd6T4c1SS1bzJFWNzATu2KwViCBgsDjtjmgD63iljnjSaFxJHIAyspyGB5BBHUGn1w/wx/5Jt4T/wCwTYf+k6V3FABRRRQAUUUUAVb29stNs59R1G4jtbS1RpZZpXCRxxoMszs2AqgDJJOAK8w/4X18D9kUv/CwfD+ycgRt/alrhyegU+Zzn2rqPHng+18e+Gp/CepsP7Nv5IVvY2XcJ7VJFeWA8j5ZVGxv9kkV8YfFTRNI8Z/tofBX4X6dYQJpPgHS9V8W3kEUSLFG7BbGwBUDAw7OVHtkdKAPs7w18Rfh/wCMry407wj4l03Wru0RZJobK7huJIkY4VnSNiVBPAJGK7KuFtPAWj2XxCu/iFbQQw3l3psOmsY4lSRkimeXLOOWGWAAPTB9a7aSaGHBldUz03ED+dAElFVvttn/AM94/wDvoVYVldQykEHoR0oAgnu7S2eGO5mSJ7h/LiDsFLvgttUHqcAnA7A1Yr4P/aC8LaFpf7Sv7PXiGytimoap4m1IXErSSPuC6Pc4AVmKqOOigDPPWvvCgAooooAKKKKACiiigAooooAKKKKAP//S/fyiiigAooooA83+LqeLJfhvr0PgfSbfXtblgCQWF35f2e6DOokhl835NjxlgS3GOx6V81/C79n/AELwF8YdJ8Z/CDwpffDjRDa3sXiHTWuAum3skir9mENok0sYkiky3nRqi7MrzuwPtyigD887r4b/ABT8OfBz47/CbS/C8+pat431PxNf6TdJJCLKW31mMtHvkaQOkiElChTJYDB2ncNmT4ffFDXPDP7OWtweFJ7K7+Gk6zaxp93cWyXCpDpUliREUkeN2d3yg3AEfeKV950UAfKUngDxD8Rf2jfC/wAa9a0+bw7oPw90rUrSxju2jF1fXeqBFmkaONnEcEMSYG5gzOc7QqgngvGPhWw+I/7Vvi3wbcXDR2XiP4VNYvPCclY7vUZY96kcE4OR64r7Y1nR9N8QaRfaDrEAubDUoJLa4iJIEkMylHUkEHlSRwc15z8N/gr4A+FT+d4TtJFnWxtNLjmuJWmki0+xBFtaxlvuxR7mIA5LEsxJ5oA+fvhX8M9Q8OfDnwz8LPFHw0tLvxD4bhs7B9Tmjgm0idLHaqXykv5u5lQOIjGHWTjgDfW38LPCPjX4CeKviHoieH7vxN4f8Y6/d+I9LvbJ4S8U2ohDcWl2s0kbJslUmOUbkKHBIYYP19RQB8ofFDwj8Qrn4ifBr4rjR/7YfwdNqiaxZ6c6eai6pZGBZIBO8YkWKQAMNwYqdwHBFfRvhJdfHhuwPikg6q0Ya4A2na7EnZlPlO0EKSODjNdFRQB8tfFjwd4vg+Pfw3+Meg6RNrum+H9P1vS723tGiFzG2pLA0EyrK8atGHhKvhsruDYIBxz+m/Drx78Pfjd/wvW6sZfE3/CU+GbDRtdt7IxfabW+sJGkjuIUdo0khcSOkiq25SFZVIJx9i0UAfMXwe+GXiK0+MHxJ+PPiu2bSLnxyumWNlpjukktvYaVEyJJcGMsgmmd3YqrMEXaCS2QPp2iigArxT4j6/8AF2y14aH4F8MQ6jpd/pF4E1JrlUe31Ysq2ySRnGINpaSSQbj8u1VLYB9rooA+XP2ZvCnjX4f/AAw8EfDW/wBHfR4vDOmmHVpbkxyNd6i+WkNsY3bMbSs8hkbGQVAGd2OP+BFx4huvHnxY+J2p6r/bT6lq9xpOnWMcMMM13beHI1tyYC8gUAXDzA8hckFiM19p1yegeBfCHha9vNR8PaVBY3F/JLLM8a4y88hllIByF8yQl3C4DN8xyeaAPmnWfiX4yb41+FZz8NNfR10XWFEBn0vc4M1kS4xfFcLjByQfmGAecelfFtLfxl8APGF74n8OtZSromqTLZaitvPLBLHbShH/AHTzR7u6lWJGex4r1G78KWF54v03xnJJILzS7S6s40BHlmO7eF3LDGcgwrjnuam8UeFPD/jTR5vD/iezW/065BWWB2ZUkUgqVbaRuUgkFTwR1FAHI/BH/kjHgH/sX9K/9JI63vGPhzXvEdtbwaD4ovfC8kLlnksobWZpQRgKwuoZgAOvygH3rQ8L+FPD/gvRoPD3hezWw022AWKBGZkjUAKFQMTtUAABRwB0FdDQB4f/AMKu+Iv/AEVvXf8AwC0j/wCQa9K8J6JrGgaWbDW9fuvElx5jP9qu47eKTacYTbbRRJgdjtzzya6aub1bxf4c0PXdF8Narerb6n4iedLCEqxM7W0fmygEAqNqc/MR7c0AeVfGfWfiJINP8FeC/B+oaxp+tzQw6rqlrc2cC2VhI+Ljy1nnjleUp8o2r8obcCWAU/Nv7Yv7OFrr/wAGx4S+BPwztrrxMb3TrmznslsrNLQWF1FMd8k8kRAaNWVQgbPfHWv0Lrm4/F3hyXxbN4FjvVOu29mmoSWu1ty2skhiWTdjbgupGM59sUAfJfx/vfFHjq8+EngkXE3hP/hJNaXUNRs7uKCaSC20BGv5JJZI5Wj2CWOFSAxUhxkjpXa/GG++OXiHR/FfgT4d6Da3VprvhlodJ103Sov9qXQkjkMyf8s4o4yroy7izMAOASPcNZ8I+B/Hj2Gq6xY22rtpxnW2mzv2CT93OgZTyr7AsiHKttwwOK7JVVFCqAABgAdAKAPCvgzpPiPQfBfhzwda6O/hrRPDeiW2nLHdiN7l7uGNI9yiJ3XyowrZJOZGbIwBk+Y/slNrOrWvjD4ha9rqajN431zUb22hMSRSGxsZRp1tMAGJEbR24IGNuWyDzX2JXJ+FvAvhDwTDJB4U0qHTI5eqwrgBdzOEXJO1AzMQi4UEkgDJoAk8V+CfB/juwj0vxpotnrlnDIJUhvYEnjWQAgMFcEA4JGa+XR+z38GP+GgDZ/8ACvdG/sz/AIRcPt/s6HyfP+2kZ+5t37fxx7V9lUUAct4T8EeDvAljLpngrRLPQrSaQyvDZQJAjSEAbiqAAnAAz7V8m/tP6D8SviJ40+G3hrQ/A2o6x4Q8NeJLXXtauop7BEuEsI2ktoYY5rqOR/8ASGUvuVQAvGa+26KAPiT9oTwh8WP2nvBE3wR0/wALz+CfDOvy2/8AbesatPaSTLZQTJM8NlbWk05aaQoBvkZEQZPzHgS+PNK1G8/aN+HPhGy8RLp2geANHuNdmSeJHVJJNul2YkkZxuaUPPtyOCuRk19rVyeoeBfCGq+IoPFmoaVDPq9vEkK3DL85jicyRq3OGCOxdNwO1iSuCc0AeIfFnwd4lvPjZ8MfiPDpM2u+HvCUWseba2vlGeLUL2GOK2udkroGVYxNHkHKlwemSPMfhjo3xR+HEnxf8d+OrFPDr6zquq+Ibu/ZoruMWNvawR6fFDtf5jFBE/m7goDDAyDmvumqOp6Zp+tabd6Pq9tHeWN9E8E8Eqh45YpVKujqeCrKSCD1FAHzH+yxpfibw58A9Pl1K9j1vxLqltJr09qwWB47zWN9+YZWyxG55cBmHA7YFev/AAh8X+KvHPgHT/EHjnQV8M+IWaeC/wBOjuFu44Li3laJ1SdPlkUlcgj1xXU+H/Cnh/wrph0jw7ZrYWx6iMncSFCAlySxIUBQSTgAAcAU5Y9J8F+G2FtDImn6VAz7IkeeUogLMQq7nkc8k9WYnuTQB0FFcx4L8Y+H/iD4U0rxt4UuDd6PrUCXNrKUaMvE/Kko4DL9CAa6egDwn9pvwl4u8e/APxv4K8DR+brOt6c9pGgcRs8cxCzqjMQodoS4XJA3YyRXb/D6HVrfSxazaedH0ezht7XTLKUL9pjgt49m6YozKC3AVAflVRnkkDv6KAON8Y+G9f8AEdvbQ6B4pvfC7wuWeSyhtZmlBGArC6hmAA6/KAfU1wP/AAq74i/9Fb13/wAAtI/+Qa9wooA5rwpour6BpQsNa1658R3O9m+1XcdvFLtbom22jijwvY7c+prxD9rPRvGfjH4EeM/hz4E8NXXiHV/FmkXunweRNawxQyzR7EaZ7maLC5bPyhjweOmfpSigDwzwvrHxF0v4KD+zfBc8Xi3QdJjgtdK1C6tYkvLu3t1VUFxBLMiIzjbuYj1xjmvQfh54l1Dxj4G0LxTq+n/2TfanaRT3FoJBMIJmX95GJF4cK2QGHDDnvXT39nDqNlPp9wWEVyjRvsYo21xg4ZSCDjuDkU6zs7XTrOCwsYlgtrZFjjjUYVEQYVQPQAYoA4T4aa58RNf0a+uviX4ah8LajDqF3Bb28F4t6s1lG+Le4LoAFaVeSh5Xv6U7Stb+INz8R9d0LVfDcNp4Qs7S1l07V1u1klu7mTd9oia2A3RiLAwxOGzx7ehUUAeeeKdb+IeneMfCel+FvDcOreH9SluV1rUJLxYJNNjji3QPHCQTOZJPlIUjb1NQfF7xP458G/D3VvE3w58PReKdcsFSSPT57tLJJYw6+cTM4KqUj3MAepGO9elVmaxo9hr2nTaTqiGW0uMCRAzJvUEHaSpBKnGGHQjIPBNAFGTXbgeGY/ENpplxfSyW8c62UBjE7GQA7FMrxx7hnncwHFfHPwc0D4op+098S/i34+8Bajpdn4rh0rStFne50+ZbXTtPjkaQzrFdO6GSdy2FV+2e9fcwAUAAYA7UtABXHeLvh74E8fx20Pjjw/Ya8lmzNAt9bx3AjZwAxXeDgkDnFdjRQB4v/wAM5fAL/onWgf8Agut//iK9Y0nSdL0HTLbRtEtIrCws41igggQRxRRrwFRVwAB2ArQooA+LPj7o/wAS/Efxm+EniHwr4E1HV9J+H+tXl/f3Mdzp8Syw3NhLbL5CTXKOxDyjIZV4BxnjPt/xT8f+OvBnhXRvFHhPwj/bXnX9pHqlrcXkVrNZWExxNMvLpLLGSoEasdxPBNeyVmapo+n6ylvFqUfnR208dwqEkKZIjuQsBwwVsMAeMgHqBQBp0UUUAFFFFABRRRQAUUUUAFFFFAH/0/38ooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigArI1vXtG8OWX9oa5dx2dvuCBnP3nboqgZLMewAJrXrwz4jfE34Y+DPF0Ft4kuWbxPpeh6lrFlbbJWU2sRSOVlIBj8122xoM+Y24qo5OQDttP+KPw51bUNO0vTPElhdXWrLutI47hGM+FLlUwcFwqltn3sAnGAa72vg39kTwr4AT4TfCubWtLttW8Y3UeoeKftCRpI+nXess81zI7Z/dMRMIB/EdpUcKcevaT4i8ZWv7VOo+BL/XZb/QJPCSatDaPFDGkFw+ovB8rRortiNQPnY0AfSlN3L/eFOr4b1T4YifU7yf8A4Z4srzzJpG8865bIZcsTvKnpu647ZoA+4gyngEGnV8aeCvhyNN8W6TqH/Ch7TQPs9wj/ANoJrNvO1rj/AJaiNeWK+gr6e8feI/DnhDwTr3ifxhdNY6HpdlPcXs6By8duiEyMvlgvuC5xtGc9OaAMyb4q/Di3ubmzm8RWSz2gLPH5o3kB/LJQdXw/ynZnDcHnivIPH/xF3fEf4KX3g7xQsmg+L9SvraaO3lhktb63GnXE8bA4LErLGuCjDuCDmvD49d+G3i/9qfRLhHgj8IfDfwBLcWfmqURpdblSNkQPgs0drCquvLKZACAxxX2P4R8BfDy18MeGodI8I2elWWjW0R0y1ltIlk09CAwRBhvKYHqFPX3oA+Jda+KPjzw74S+MXi6PxxPcS+APGdtp+mxTNbmMWs7WBkguAEHmKRPIq9CoHB3Amvoi/wBe1i8/aRv/AIfweJLiPStV8GPfJawyQhre6W8WAT252FwxjOfmLDPIFe0v8O/h/LDc20nhnTHivZBNOhsoCssq5w7jZhmGTgnnk81pQeFPC1rqkeuW2j2cOpRRCFLlLeNZ1iUbQgkC7goXgAHGOKAPlL9h+9gX4H+GtKudbvdW1ZrFrm5ju3WXyCbu4jI3BFIdnRt4clsgmvsysnTtD0PRDdS6Rp9vYG8kM1w0ESRGWQ8l32AbmPcnJrl/ht8UPAnxe8Mjxh8OtWj1nSDcT2vnxq6YmtnMcqFXCsCrDuORgjgigDvqKKKAEJCgsxwBySa85l+L3wxhS6eXxNYr9jAaUGYZCNuw4HVkOxvnXK/K3PBp/wAWPFPgnwZ8O9b8QfEa6ez8OQQiO9kjWVn8udhFtVYQZCWLBflGea+P4L74d+MP2nPiV4i8RRW914f8F+D9P8M29sUBF0upyTz3UcEfHmAhUt8Ln5tydcigD740/ULDVrC31TSrmO8s7yNZoZoXDxyxuNyujLkMrA5BBwRVyvm/x5oNn4K8L3/je41K88N+DPBHhx2tNH024ayiT7NE0jmXycE7ERI41Vyo+bgkiua+BHgT4h3fwv8Ahd4u8V+LtXufEb29vqWtLc3sjxXC3du7m3aE5T920iAYCn5M5znIB9a03en94Vj+I7b7Z4f1K0+wDVfOtpU+yM4jFxuUjyi54UP0yemc18TD4WLgf8Y32Q/7j1tQB95BlPQ5pHdI0aSRgqKCSScAAdSTXy/8IvBA8P8Ai37ePg/beCP9HkT+0ItUhvG5I/d+WnPzevbFesfF3xT4E8H+ANS1f4lXL2vh5zDbXLRrK7P9qlSFI9sILkO7hSAOhOeM0ALN8Yvhbb291dT+KLCOOzAaYtMAVRgWEmOpjIVjvA24UnOAceL3fjnVPGPxI8feHr3xA3h7w/4Z0bTr7SZLaSFF1GO+hlklvTLIrrLDGyCMKp2DBL53LjybT7r4deOv2iPjL4o8W21vqmh+HtA0vwbaWRRWN9FciW6u44IzjzVkldYF28FkZc8GvuH/AIQzwhcWGm2VxoNk9vpcKQ2kUtvE4to1UKI4wQQgAAGF44FAH55fB/4g+LrT4S/s4eC9BlEOgeJ9Eunvb6G4S2aa5s4kMFnHcMHWJpSzyEDDsIyqkc16JrPjr42Wh0bQrS7sfFl3CNYe40+w1iLTdUltbe5RLW4S4khFtPLbofKnjzGpkIYk8gfZkXgrwbBoY8Mw6DYR6OrbxZLaxC2DZ3bhEF2ZzznHXmluvBng++gt7a90OxuIbRdkKSW0TLEnXagKkKPYcUAYXwm8Qf8ACVfDLwx4i8+6umv9Pt5WlvolhunYoNzTRplFcnOdpK5+6cYr0KvO/HXxP+H3wpj0CLxrqkWixeIdRg0fTgyNslvJwfKhGxSFztOC2FHrXolABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAf/9T9/KKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKoXWlaXfXFvd3tnDcT2jboZJI1d4m9UYglT7ir9FAGbp2jaPpBnOk2MFkblzJL5ESx+Y56s+0Dceep5rxeP4U+L1+Ob/GU+JLQwyaWuinT/AOz2/wCPJLlrkHz/ALRnztzEbtm3H8Ne9UUAFFFFABUcsUU8TwToJI5AVZWGVYHggg9QakooAwl8LeGVFoF0izA09i9ti3j/AHDMckx/L8hJ5JXFbtFFABRRRQAVyPg7wyPDNnerJ5RudSvJr2fyE8uIPJhVVV9FRVUk8sQWPWuuooAKKKKAILm1tr23ktbyJJ4JRteORQyMD2IPBFZaeGvDkc9pcx6VaJNYLstnEEYaBf7sZxlBz0XFbdFAHhvx/wDhR4h+Nfw/1H4baf4ii8PaVrcQhvnNmbqaSISK5jQ+dGqq4Uo+QSVJwRXr2iWt/Y6TaWepywzXMEYR3t4jBEdvA2Rs7lRjAxuNalFABRRRQAVXu7S0v7eSzvoUuIJRh45FDow9CpyDViigDGj8O+H4rm1vYtMtUuLKPyreQQoHhj5+SNsZVeTwMCtmiigAooooA5Lxl4aPizTrfSJDCLb7VBLOZY/MbyoW3lY+yuxAXd2UnHOK62iigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAP//V/fyiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooA//Z" }, "image.png": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAl4AAABqCAYAAABtVV8UAAAAAXNSR0IArs4c6QAAAHhlWElmTU0AKgAAAAgABAEaAAUAAAABAAAAPgEbAAUAAAABAAAARgEoAAMAAAABAAIAAIdpAAQAAAABAAAATgAAAAAAAACQAAAAAQAAAJAAAAABAAOgAQADAAAAAQABAACgAgAEAAAAAQAAAl6gAwAEAAAAAQAAAGoAAAAATAzH3AAAAAlwSFlzAAAWJQAAFiUBSVIk8AAAO9hJREFUeAHt3QmYdEV1BuCLa1xiNHEnib97UAEVN4TAj4KyiYDsKgIiKrKoEURARYMGUFFBUAQEBREBQTSGRVEQQQxBEQUXFFBjXGL2xeyd+9aTM7ncv3ume2Z6+s70Oc8z0zO3763lq7p1vjrnVNVavVqqlEQgEUgEEoFEIBFIBBKBsSNwl7HnkBkkAolAIpAIJAKJQCKQCBQEknhlR0gEEoFEIBFIBBKBRGCJEEjitURAZzaJQCKQCCQCiUAikAgk8co+kAgkAolAIpAIJAKJwBIhkMRriYDObBKBRCARSAQSgUQgEUjilX0gEUgEEoFEIBFIBBKBJUIgidcSAZ3ZJAKJQCKQCCQCiUAikMQr+0AikAgkAolAIpAIJAJLhEASryUCOrNJBBKBRCARSAQSgUQgiVf2gUQgEUgEEoFEIBFIBJYIgSReSwR0ZpMIJAKJQCKQCCQCicDdEoJEoAsI/PM//3P1H//xH9U973nP6h73uEf1P//zP9Wvf/3r8v+97nWvLhRxKsugHf7pn/6pcqSrtrnb3e5W/ed//mdpq/ve977l/6kEJiudCCQCicA8EZgK4vWrX/2quvjii6t11lmnevaznz1PqKoKOfj0pz9dPfShD61Wr149FqVz3XXXVbfccku17bbbVg9+8IPnXdbl9OB//dd/VRdeeGH17W9/uyj5v/7rv67ucpe7VA9/+MOr7bffvmC9nOqzksrq3fnIRz5S+fTzL//yL9Vv/MZvVI985COrV7ziFdXv/d7vraTqZl0SgUQgERg/AvVMdsXLmWee2auR7D3rWc/q/eu//uu863vFFVf06hl/79GPfnTvL//yL+edzqAHa6XW22abbUpZP/jBDw66bcVdr60qvb/7u7/r/eIXv+h99KMfLfV/4Qtf2PvJT36yoPZacUBNoEK1datXE+Hez372s95+++1X2ua9731v75e//GXPdymJQCKQCCQCoyEwFRavJz/5ydVWW21VbbbZZsWNNV86+9jHPrbacccdq5p4Vfe///3nm8zA51gSnve85xW3ztOf/vSB9620L9Zaa60ZPGtiXKqnvX73d393pVV1pPpwvV555ZXVeuutV6ysIz28SDdzLT7wgQ8sqbFE6vdbbrll9aAHPWiRcshkEoFEIBGYLgTWwtOmocrihRAbSn4h8m//9m+FvHGFjUP++7//u+J6E08zbYJo7LPPPtVnP/vZ6qqrrqoQ5mmWv//7v68OOOCA6lWvelW18cYbTxSKH/3oR9Uf/uEfFtfiJZdcUt3vfvebaHky80QgEUgElisC42EPHURDgPZCSZdqIW/Dkq5BnHbQdenf9a53XRTSJQ8EbjnJ3/7t31Zf/epXqz/4gz+oHvOYxyynoo+trPraYvTbhRbwe9/7XvVXf/VX1VOf+tTqN3/zNxeaXD6fCCQCicDUItBZV+O3vvWt6otf/GJFGT/lKU8p7o1rr722uu2226p/+Id/qOp4rWqjjTYqq6u+9rWvVT/4wQ+qv/mbvyn3rq4D3xEYUsemlMD6n//85yVYe++9967ufve7V9/4xjeqL3/5y+UZM3kWhWb60t5www1nOobn/+zP/qyqY7uK4mGZ+a3f+q2y8u5P//RPS/6Izm677VYhedKS973vfe+qjlcqgfI//elPqz//8z+vfFKmroc7zeqxOoasuvHGG6t//Md/rLbeeuuZ/GHxzW9+s9TpPve5TwnqVz9uOXVRXyvMCIuZckqHm8iPfB73uMfN1KWrf9x8883Vj3/847KwQD1TqkUnXd6Vm266qUD7pCc9aaaPzYW1/oTMe0+6QATnKm9+nwgkAolAVxHorMXL4C6m5MQTT6w++clPVuedd17F9bLJJpuUeJd99923uuCCC6pPfepTVR2UXVYrWmn1yle+svrEJz4xg7d0rEb8wAc+UFZnhRWIJcFz73rXu6qLLrqopGXF1qabblpI0ktf+tJC2CIh6XCFnXHGGaVM0iRhvRKLc+yxxxbSg/goi3gtFhwk7Zprrikr96yI3HzzzcvKRSSNFSFEmSi4o48+uhDDuP7hD3+4es973lO2WODiecADHlBI3ctf/vLKdyHKf/jhhxdMdthhh4LFQx7ykJI/8tZ1QUq1D6Kdyn1xW+vf//3fS//0zmywwQbVuuuuW7397W8v/X6unGwfYaIivuuJT3ziXLfn94lAIpAIJAKzIVATh86KVVP1dgK9msT0ajIzU07XaytPr7YW9c4///yZ67W1p7fddtv1tthii14d0zVz3R9WyT3jGc+40yq5Ol6rV28v0XvCE57Q+8IXvjBzf02wenVwe+9FL3pRz4q7prz61a/uPepRj+rVVqvm5d6ll17aq61bvZe85CVlFVh8Kd16X6pyvSZZcblXbxvRqy1mvXqp/sw1f3z9618v10866aSZ63/0R3/Uq4ndzP+1Eu3VxKqsrqy3npi5XpOz3tprr92rSdbMNX/sueeevXpRQM9zXZWacBW8a8td7y/+4i+6WswlLVdt2S1tV5P2BeUL2z/+4z8u70ZtFZ5J64gjjujVJLdXT2hmrvX74/bbby/vSG3t6ll5S+q9vXo33HBDr7bO9nskryUCiUAikAgMQKCzFi9kkduMFciGmk972tNm+GNdl2L14Z7jcmyKva9YyszwmyKdtkhHGlZttVcR/s7v/E5Vb3EwVDrSlb48xSbFKjDXucxYb6yEfNjDHuZSEdfViwuzLc2ywuC3f/u37+QqfMc73lEJcH7rW99a9ibzPDfrxz72sWr99defuRbpWs3J4sZN21XhUhZHBKdx7Q3FSvmd73yntPk4ceAq/u53vzt0FiypNilt/igrV7L+qV/5v/m9vy0Y0YfnEgsVWEbf+MY3lr4U93NP16SqqrftiEt9P7nxv//975d3hOuc2G+u3vqkuMf7PrSEF4Ue/PCHP1zCHDOrRCARSATmj0BnY7yiShQPd1mTzPjOdfFR/VZXISvDKCTpcG3VFqxCkPwfgizJw88wIj+EqbY4rXG7gPxHPOIRa1z3jLI2pV1uaXJ7xvJ98WTco6997WurPfbYY+ZRikf8GZw+97nPld3FfSnGi3J1nau2q4IQie/adddd12jrxSoz4mmj1iOPPHLoBRLzyRvpuvzyy0s+wzx/2WWXVbWFsxD0potV37SZLuKDfDf7BrJmM2B1QeAHif5b72NX+nhz8uK6/qL/tScp7bTEhSmLuMcQW1zU+3mtMfGJ75fyE97c1AceeOBSZpt5JQKJQCIwLwQ6T7zUCvloKqSoqQDzftfj+2E/2+lQcE0lN2w67oug/nhGOsqIAM1HPBsWIFaH2u1Y4twOPfTQkhdLiK0nkCpWECTVqkBxOSH2H9trr73K4oK41rXPWFRQu4P7kiJWJARBfFsTSwQiiLaFBvDyvx9/u0ZYj5w6wDLoGfhIxz1IhftjCw84+rtJaCJNaUU+npOW/hllEmcnxtD/8pB+u39JoykWccQii7juOXVGsu0dZzWhvEL0KzFXUb+43v60DcTVV19dvexlL7tTfdSxdumWiUL0r/az/lcH1lKW5GZ8F4uwGMXmexIYBR6wcK05OVIHuEXbaFN5xOKQKEOk5f/Z8Nb/zz333Ep5Au9oi0grPxOBRCAR6BIC82MDXarBIpSlqTwWmtygtAZdHzY/SoyriKI5/vjjZ1ytyIRgaZY2y/wpnS6vYLTSEzFhfbFQglDEFiFYJfrMZz7zTpD4jkUIgSCUuufquLzyv9WeFjwgSocddlhRwDDxIy37YCFdJ5xwQllg8fjHP76qTwUoee20005FqVu8UcfiVc9//vPLKlSuN3hzPz/3uc8tJAH5OPvsswvOb3jDG8q2IuqhbM95znPKAgZE6X3ve19ZxKFNPvShD5V86ljBNSyqpfD/94s1sm3R9ZUyuK6uXMjzEYsqpGPlblNY0qxuPOigg2asqdyH6rnLLrvMkH1WMRjrU2G15eJ0jBBLmIUjq+tVtcSiEotd5IUoes4Clt///d8v5FEbIdgscP62SS5LLZe+vrvzzjvP4PCVr3yl+vjHP14mCyYZSHB9qkFxbdZxnIX0cTFadHL66aeXPLWrxSuOmWqS5lK4/JUIJAKJQEcQWDPwqSMFUwyzYj/9JL5rfh/X4tnmc8374npci8+43ny+/V38H5/xTPwfnwu9Hs/H5/vf//6i5I877rgZ0sF6QHEhIxQjy0kdnF+sOvFcfIoJ486bpCCfSCOrXXM1JqLz+c9/vpws0CSNSJf6IgRckKx2VjxajWe7DiKmjpWPCxHBIKxmFP6V9UpTQgl73mpSBNU5mOLeuIBZU7jsbGWBKInD48JFJt7ylrcUd5py23oBOQg3LszFF95xxx2FrMjHNiK77757yR9JEQMlHfnMR1hwwkI0n+c9wwWn3KyeTTnnnHMKblYB67PiIl/zmtdUSGVg637kyspbK3TjsHKkidvRfaywIYgqrKwg1qYI0gte8ILS5p/5zGfKbcoBY1gjYbZN4UrnDrdKN2IerbqUD9c6HFgN4Y3MsdQR5akXs5T2ExsoL22ZFq8CT/5KBBKBjiLQaeLFDUGZ2i1eTEsIhczy47vmdQO06/3uZ/Vw3bMh/vfju+Z1f0tb/s30KcF+1ykb1+XfDHh2PfL1vf9J87rv29flKe+4bif3Y445pjr44IPLTD/SsA/Z9ddfX1w5CIRge8+ecsopd3JLOXyakkNQJinwE0RvXzFWDAJrVguKVf0o2xCK3YIBlhCuJIHdrFgOO7fAQF2RCvFGzee4zij4cPtSxLb34JpjVaOkESPXESjkLQ4+p9xZmZBYh0AjvBFjhQyEaw9Z8ZyfyMd3/uc24xIVO8jaE99HvZbqUx9CWvXxZnwfy5Q94/QpuBD9Vvlf97rXVfWq4HKNdVL97d2F9IbAjDvbIeb+DoGh62R1TTjFpa1atapY+1jYiPZyD1xseYIII2KsuXfUJPbkk08u93m2jTdS7rnIE6HWltLUP9SFyza+Lwnlr0QgEUgEOoZAZ12N4lLOOuusMhNGQCiEeluEYr047bTTyoyaYjnkkEMqm6JSdKeeempxn8R1s2HKwP0GddfN6OstIQpB47pA1ATnSsf98jqzdoWY5fubQmABkQ73CosSpfvmN7+57JMVLiXKjAKpt7coik78C+uIGbuycWdxqYi14dZiAaAwbBKLRLCUWH3mupgY+y0hBvJmQXCv69xqSJxA9HoLi6LYEArCfYeo2E9M0DWrBKLDXaQ8lNkkBQFheaq3ISj1sYqR5YVVRhtR1k1h2UAaYmFBfEe5qierCsWL0LUtjUFa4xmfrsX1+JuSFk/k/zZBovgR41tvvbUQOfc1JZ7rd02ZSHxOggxod6SRtUrf8j8yJfD/T/7kTwoBjrLru9pGH/FOqDPrJBIqiL4dB9av7tIy+UDIwjoGV20TuLvHs7Butpl3ZFVN0ljYQobFG8aRfvSFZtqRXn4mAolAItAFBDpLvChUm6SKQSEUYMx2ETDkyeBKSVAKLB5ICncFJYecUCasJC9+8YvLdfe7bpZNQbBoWB1o0Ja+dNyDHO2///7l77hfOmJ15EFpuI4AIEfcVuKOKBtEQf5IEreHWBfXWWcoBRYXrjBxO9JUDgqGgmNtEcsT6bBA2HaCRUi6kQY8uHUoSjP9piuLRUjMC3IjgFpduXPc1wXh7qPckS11gsN+++1Xytkun4Bu2GibpoilggssQ7SbnxAku/m/60EC/C2+SRtqo5A2OWIl8oy2IaHc437/q0MT/7gv8rbxKMsMEjcfYRFs5ztsOqxMQfZZANWZpU8cmnejKfqa9wHxEV8Fv71qK5d4rXbgezzXxrx5vV3mwCPu8dnEWzvDe1VNvprSfE6ac7VrLAQQy5eSCCQCiUAXEegs8TJr9tNPKOR+MohcDLo+KP1BLrn2Xl9RhkGBz81VYHGvz3DHNK/5u0kkmt+JLxpFWA8QnC4KZcsyx1JImoq1XV51YJVkIUTECeXMxcqCF4QICSNhXUI2xGxFMHj5sv7FMsjySbhfm3hT6uKMEB1kGsll+USogzSxOkYe0vjlL39ZtntoHi3lWemaEBDB7VH2cmGEXwgbnKKeIzxabmUxRHK5PxFwP7MJgmUxgTo2SVG/Z0w+TBqaeLgPjq7DgbgPpk3rlTbXRixwJlOE5fdXv/pV9ba3va387xe8m885aksbmbiEIIzwFmhPEF2u0ZREIBFIBLqKwF2PqqWrhctyrVwEKN/ZSJeaIwKsgIiWoGokhiv2yU9+crFIUroEuaDEuYwFiYc1C/GQBuKKACBe3LPuRaxYAlkdCVcwUkfRc4EK1veduLkg7gi/BQoCyrns5IdgsS5yh7LeqZN8xFAhCkgHAhmut5LZkL/UD/GS9lxYtZNkGeLCRrZYgdV/WJkrLxZHVlXWMfUXe2WyIhbx4osvLgQKcUN8rUx0H+srC/Gq2qIFN655EwQ4WtFqcYVFF4hfiHrD2/5rsOQ2jfZVL/kqK+IlPWmxPkoj+kaklZ+JQCKQCHQFgbXqGercW193pbRZjqlEgDvP6k0B4GLtKPC2IBqUNGsL5Y5AIVisRQhAkAnKmdWL1ZHiJ6w2VsSxpHC3WVnHVSw4vm354Z5EHLw2vlcm9yNlzTgo8WfiBOXTtKy1yz2u/5FD21xYtSgecTEFmVQ3Vi1/q58fK0kRWiQP/sgYq2BYvRBZJNaWE/URRoWUKZc0tFHEKjbLys0cq3HhzbJlLzau0qbF2nYj8kd+B7lGm+nm34lAIpAITAqBJF6TQj7z7QwCSIIAdPFx73znOztTroUUhEvOCkHWrkGu7YWkv5Bn7X9mqw7Wx/m6YReSfz6bCCQCicAkEfj/teCTLEXmnQhMCAEWK3FkLFncmVxWK8EIbGsFu953jXRx/1rty/p1Zh1Dx3WYkggkAonANCGQFq9pau2s6xoIsHZR/lyRXI4C2iOma42b88KCERAfJt4LztyS4rj8nZIIJAKJwLQgkMRrWlo665kIJAKJQCKQCCQCE0cgXY0Tb4IsQCKQCCQCiUAikAhMCwJJvKalpbOeiUAikAgkAolAIjBxBDq7gepiI2PJuu0FxJiI67GppA0aUxKB5YaAfa9sYeGsw2kUW0yIy7NNhe0+vM+xkauFEba5sI1I+0SBacQq6zxZBGyFo29G/9Rf9Vv9VL/1Kb7Uj+v2n2tvYTPZGmTu40BgamK8DNTHHXdc2eTRnk/nnXdeOXZnWFAdcWP/oV122WVm1/Vhn837EoHFQsDZm9ttt10hFhdddNHM/mSLlf5ySAfxdFanvd1MpuwBZiIVe7VZqWqvL8ck2cfM/l8p3ULARNimuTY/to3LShN98LDDDquuv/76smWKfQERK5s067MWlFh5HHvh2e/OHoP6qyPrJiUWvtgI2VFhzX0JJ1WeFZtvzbinRuoZRq8+HLhXb5zZq3fTHqne9X5INprt1Yd1j/TcqDfXs/leveS+V5PDUR8t99cbW/bqzTPn9Ww+1H0E6nM7Sz+sN0ft1ZvFdr/AYyhhbT3oeZdf//rXFyzq3fHL/7UyL5+10uvVh9P3asXRq4/56tUHyo+hFMsjyVqR9vSZ+vSHiRdYu9RHQ/XOPffc3iGHHNKrN93tvfvd7554ucZRgOuuu6633nrrFX1TH2PVMy7Xk4FefU5q6bP1+cG9epPlXk24ymd9LFqvPkptInjUGw/36m10evW2Or2XvOQlvZok9q6++upxwJJp/h8CUxXjxYRr9juf40QctG1Xcwdlj1NY5uqBqRwGPJ98vv71r5cjdubzbD7TbQTs4O74HW40Fh3bMSwnufzyy8tRSgstM8sWF47d+Vm7Nt100+KeCdcjF+M222xT7VUf8s3icMEFFyw0y2X7vPHgyCOPLGeKTroS2sxZmqwqq+rTJXxqs5UorLKOwGL1csSZkxZY9+pJQOmrO+64YzmBw3YqTuJwzJXNjidhnTWW3HjjjcUKpzxx+sRKbJeu1GlqYrwCcH71+chmm21W+Rm3eDEdaTPfAckmoHH+4LjLmukvHQI2eP30pz9d7bnnnmXwdjyOQ6XjkOmlK8n8c7J5qklP8zzG+abmyCIklEuEUusn3D1EXOe0CuLFzeVg+kmLI6FqK2UpRm0RmvcYN+l6zJW/A+AdbdV2GXqH1dvEyfmxTakNIYX4DOrLzXsX++911lmn8kOMKcqSMl4EloXFy0G8Zkchi9ExIoBRrIHz/RYqCJ0BXvzYbKLs6hJKoXmvspx99tmVF3c+5IkFwKHEw1j0lFN55Tkp6WexWYy2mFR9xpkva5E+i/xTYN4JM9X5inaPgN9Io197xHeL8WkyMd8JRTt/B6azDm+44YZ9g5EpEJjBaquttmo/PvD/NiYDb1wmX7BkiHV7ylOe0qkST3LcGTcQzg11ZmicBRv5ieG69dZby+HuzXNGfU8vsOQ683WSMl/DxCTLvBzz7rTFC5FgskUS/DiI1+HDO+20U/X4xz++uuyyy6prrrmmKCAz39122620wZe+9KXKDwLz1Kc+tcw8gmi5wd9WRtV+7PIiICDS5p5od3wvRO1/L+476dUxI8XqEIG80pPXlVdeWQgP5bXBBhuUQ5ebebpPfo6koXwc5Pu4xz2u2mKLLcrqyptvvrmcrecYFebmk046qZSTefqFL3zhrCu05KmMxx57bHEpfOUrXyn3M+2rP+UUAkfBk9/+9rfLJWU0G6acmnWK+8fxCUflZbWAA/ctwnXJJZcUUipQmlVE4GlKVfq3fr7//vtX97rXvcrO+mbP8yFe+nMd71P6if6x+eabV0984hML9iytVvpuueWWa7wHXWsH44IA5X6EgrXhTW96U7Eceyccfj6MUIpXXXVV6YMWMDTPkYS398q70vWTDVhVEC7tq98o7+mnn17GOO95ruYepjfM/x7jl7G9LUIFvLPrrrvuGge5G3uN89k2bdRW5v+d1Wzf/e53y1lzlE2YQZnNrUxcvXp1IV780Rj6hz/84dI6QbwMNPe85z2rU089teJLb5p8dXDWpjPOOKMsxzcQUfrnnHNOtfvuu5f0EaemMA0jCO9973vLqhTuHkKJnXDCCYVMvf3tby8KzGxHmRGbww8/vJAZs+jjjz++uvDCC8s1JBFZqoNeqzrYtKSLjKmXGAiD/Nprr12eRQjbBK5ZtiiHF9aLa8D1jBmVPPwd4qUXcyDtffbZp8zIbrvttur9739/ZdUmZQW3cYoynHXWWdX6669fbbzxxgUXsQ0wrwPGi+KDL3I4rMIcZ3nbabOkcPONQlL1E5YXbTLKc5E3LJDniP/Qv00cRiVe+qF4J/0LsdDXDjzwwII78mUyI45RH37zm98c2Xfuk8VYf/VOejf0YXWzYoxyY+kiiPyzn/3socr/05/+tKoD8qvnP//51THHHFPSZ30OCx1Ctt9++xX8uk68jAXe8Ztuuqm0pfdJ34PXXNZwq+tYEkfpp/o3TB760IcOhfVKvwnO/cT7Rkze2wJvk+yU6UCgs8TLQGcQaJIgilqQuxedsOZY8nr++effqbXq1SQl9oXbrd8AQmlttNFGRfnEg4ceemixAlBEYmnC8uV5s2rkz8DcTE+MCSKIuMSSaAqM8kJikD7PsYYZzBG3bbfdtmRJebD6UJ5HHXVUOczYLAkB9N32228fRZvzE1lC5gy473znOwupabtXYCb/733ve+UzBgcE6C1veUu19dZbl+XNTZI6Z8Yj3oAk28ZDfEPEy8m/XtlUSOOznvWsQk4d7IzsLoXApdmmc+VJqVv8MIo1jlsFcRLfMiqxNQHRZkccccRM0cKFYWuJUUR/9U696lWvKuUP0qIvIFqIir+5psYhcJ5rEjFMviylMBHf9oMf/KBYtrQjy7U4MhZVk6hh45qQNlsbaCMuIriatGi3IF6s1fe5z33KOzJMGSd5jzHIzx21BfMhD3lIwYIVZi6BoTHONiVzEbRmWnBiuTd2pvRHgPWRlda728UJZf9S59VxIdBZ4qWDIl9vfOMbZwYOSmGPPfa4k/IyWMTg2ASJFWvQIM9dY3bfFPfuuuuuxe2FsNXbRjS/LlaxSI8CMVjXy2/LPQhDU5BF+TPzI14nnnhiOQyYWydEDJe6GeBDkXIDSteP66MeHixP0i9+wuz3Qx/6UKlXkK4oy6p6hRFli5ghYCyJ4xAKEXZN1yflYKbM2kUQAG0qILgp8EAaBJ/2a2/3wg2ppiAHkSlpGwRZR1gbPfPSl760mdWsf3OLcsWNKlbajUq61NmkwoRDX2Cp1QcDG3gOEuSEtRMWIdz0JhxBGsWcyON5z3tewUv/5DbvZ9Fh5dNn53KFINcmDj6bbeBvbSM/5dYOIfEOS3tQ28a9Plm1WGUOOOCAat999535KvJDjk3QuFRZlZtllhdsWBeMA4QFEK4mbAgdaxoSEe0VhM6EAZFpC2uZd3g+cZnttBbrf/3axA45bccTDcoDfsbXUSZ9kVaX6h5l6tKnd80YrA/Nd0GMvquv0hPR1+eqo/uM96MQ6bnSzO8XjkBniReT/84771ysSYgLyxbFYIM5A+RChWJoiwGKgkSYXvva187aucWRiMsiFLiBx2Cno1PqrDY6OyXD7WjG2XT7UTDqt1Qi7oN1LSx57XwpFDFz3DbjIl6UExdnEEpKGNZm55FnU0k2y/i5z32u+shHPlLclNx2/QTuyCOlOaieBi0b4VK+LJsI+CjEizIOhdyvDIt5jeK84oorSn8yCTHwIk0GcDLI1YgIsFy+9a1vLZbQKBNXdHMA5pbWZ8Naq136WRrDSsZS5r2cTZA7FlyKJgie+xEbMVTeD5Y3+Yb4W5ty1w9jbeOykR6Ldz8FxPriO3FN4jbVO4QFUSgAQuYegpyadGlXmCtP02JspTGsbMvQFm5ZrjyxZO0JWPvepfzfyQbqilSPQoqMf35SFhcBusLkRV9pT3yHzcnYZXwzoYgxdLZnjRfuE/KSVrbZkFr67zpLvJCUU045pXRUbgCDLQuTOKnPfOYzff3kC4UPGfLDIqLT9hvUIw/3cElQGOJImi+C58waEQ0DMwsAhTfMbD7Sj08vG5mPslcHVhIDr1k7GVQG1+XVb7VleXARfsmjmT/yw3rRtmL2y+rSSy8tVopBpMszLH7SDMz6pQNHhFd73X777aWt+9036BqyOFv6g57T/vKerU81nxVLxr3FjY1EBlFBOJBVJBkhMoFoYioNxAzxac+sm31IegiQPYRMEmYTbcQSNAwpkh63Oqtrs67+PvnkkwvurFT6Zoi/kbR+lra4Jz6t5mSRkk8zDCG+9ylvuKgjy1hTBMgj6E0rEExhwyKnn1FSXOAhxh7f9YvNQci0lf3EuiQWJnkXRlW48+3f3qdm/+oSFl0oS0w2THKa78UoZUOIX/Oa11R77733UGl4r/TtYd6rUcqR9y4cgc4SL6yeFcbmiH4QHfEbArEFxscgOKgTG3TD9TYsTMy4iAfrmg47m3D3rKpddLfcckshBP0CIylonZ41hwKQfntwCgXRtETIN+p17bXXFsW6ySabzFacvt+pv3gkgdQsbtIMAtZ+AIlUh7A8tb8fx/9mgRRp0/XIkkiRWVEGG/ew4FC2Bi2WO1aqcM+2y4WEBHbt7/zvu2gD6c92b7/nWcnE4TWtOf3ua15DBMTgscQ1CXrznvbfF198cYktbBKAuEdbSkefUocgXnfUbltE38REf0KY1K9NwKSDfFAGJg3hujRQs0rp/9xwrFcsoHE0EWVOIsi//NP6BZdBQdYmU8rdJD2tx+f8lyXHIht4DlIo6qDu+rIVZATJhg03pLp5b2EnCD0EUUGytFPzfUR0ubj1SVZG5MSPNGFNIdrewpjTTC/SncSn/fxYVsI7YCxyxJIYttnkk5/8ZFlA0Kz/bPf7Tr1Z1lhEU9ZEgO4ybhmz5sJ/zafvfMUY3U/X3Pmu/K/rCHSWeOmoyE/EYhncxLqYLRt4QwwQ/ZQg5W2Q9H2bRBk0+8VBGay4HcSHzCYUVNxnsDFYtzeFtGM2xShmSnrveMc7ium/vcpKQDPFaICkQJUXYQpCgJg043RmKxelZhBUPkKxIFrSW12vmERwwo3aTIc7VN1f8IIX9F0G3bx3IX9T4oiL7UDUOUzmXI0hXD3KHcQLEVM2CpfC166DFD/8tHUQkUhzMT/FooUyGyVdLtRhlRkriqDxfq4teSJK+gSyDKsQxN45cHC1W7ZJhGuE4hW7CHcue+5v5OSggw6a6Wvax+KHiG9k5WE5u7JeHILsNdOLPEf51Deb5R3l2bg3LEyIV792Vsajjz66WHssNIkJGhwQU8+byHkvKMSmIKNwam4F4Jp32TiE6Om/CB1yiZiaGJkQsGrLuyn6qvdNGRZCNptpDvO39189tVmsCFdOfWouxc+VbFFCjD/D5Oee2SzR/dKIMbtfG7pfO4g3NAlojg/ttLwDLMMI9qBxwTNIt3HaODhbWY01XILuGzRJ0oe5/nkShnEve4dMIJVRfbos0R7RPu2y6lsmbPSneNFR+0k7vWn9v7PESwPbDkIgc+zyizBRwE2SQwE5bNQgaGCneA2oFLg0uGMQILMEnco9ZrxeaqsVY9ZsoGJJo3QQvGFkhx12KC4f8R0GuXA3iAmhrBAZgizKrz4nsrhPYxCm/OQbFh8vuroinV5+s2fWHh18GPFSU6zSJZSGFyQUtRgasWsCtiO+DKaC7lkjrJwblhwMU572PR/84AeLuwn5lI9ympUHMebOovzr88LKo/CgYJEJWIgPixgk7ah+TZJqENb22tvAof39+NtsM/Jpl2uU/+UfZRjluWHupbxZa2z5wcLZb+CnkNSPUOzIBIXDikdBaUcrbcWtNRcB3HDDDSW26Q1veENRKt6PZqwJgs99qS8G0Q/L4rve9a7iOtffJyXqrYz2wdN3kCNERz8g2h12YrcoBu/awQcfPNPm3ivPwUk/ak+ApMGl6wfZJ9I2MdIntYd3idUI0fKuSUsbiaeLVbrlwf/75SBvx8ZwDQmbGOe71cxXueHBSqd8xkbu52GOOwsMmukt1t/eReUy5txR91vvrnfbBNnE2k8ofmOg/oY4iu8ctKM7HWEFuXAFePd7x+GhP9h2SLwwd10/Qcb1DRMSlszmYqjm/Sy/ti5SXhbUfq56dVNP/dZ7hXzpQ/6Hg74TdW2mPYm/EUnt4kd7KDdrrvdfn1XPIFj0r/5urGUJH1ZXTqJeXc6zs8SLhcDMywwzluHq8Aa8pkVKx7AVBIuSHwOyTu4+zxuEvWgCDFkBdCLL+r0Ep512WlFaXnyxH1wMsRfYMI0mLS+ygEcrFJEESl7+SFe8kBS1+DQETWCv4F8DBEWrLkH+5GkrCoNkfYhsmSlTgmH1m6tMyCXFc9RRR5XVgZ5F2kKRIqwGH4MVqyHCClOK5hOf+MSsM8a58h7mezE56k55XVkTUysYta8yGWARCoSwPciaXcKyOVNFOsQSUSowNzAYOMxYETJ1NuBSlkibQPNB7slhyj7uewxk+geCbrKgT8KkqSzVqz5kd+Z9UG9KRxvru/qj+hv0230GBgLN9TUkgBXGp+0DKGcDqmDzNoHgkmMtmo+Vb7Ew8554x9RNv2Vt8j5RjmFBU0b9Rz+hqMPF2CwDy5V+0bRoNb8XN+Y9lq73Ux+Tl+1W9FcTBxYjWBIWFDIo9s1kDKaUmbLFhKs8NMZfFLrtNJTXghT9hiVL/SYpiLN9D72zSLI+6t3W972bFHpYt4SZhCWRtao9JkQ9pIFUx+Q8rjc/jQ3aTdrNTXGb9/gbGaIzTMxnayvjkLIZx9Qjxvlmely2LHHeRe+QcqqHcZ3OEmQ/jLWsmea4/jaxMAbHQjDvBx1hLFpVh9NY0BZjr/dH3KC+T08l8Zpfq6xVs+//j3KdXxpjeYri0cgG1oifQKb89BOzFeZcMwoDnheHRQV712kMll4oL4JB1eAkDy8OCAyes1kyDP724GJRMrOKGUCUxQzOi6jMFCbLQz9BdLyAFKAXtt+KIzN5ljADAQUy6kxZrIoVTYiVn7ZwIcEK4aJYlIMSXwphyTPwKhcMECMvsHYxOLZXNRqs4Y6AIBfaE6GCifKH4tUelBv3EqsOJeM+oq31gfZs2IaY8qecJi3qYeBTRnVTT+3ftErFPcrqO6KvqB/SDYOjatLNvYV8SMf30ce0O3eT9yDIgkEXMWA166eYBffbdNWs3XPS89nu/6Uwc/wyMULyKJ9RRDtSzuqvTjCKfhDpuKaeTbziO5/ecdZe77HVjtKUhrq0xXvM0mwMMZYQ74s6e7ej7kivvhtEzdiD/DbFZMB4IQZqVa3EllKMbd4JZe7CpAPm+rhP44221KbefeLdj37tf+8mC6cxarZFAki3dmy/39IIkae+a3yeTeTpp92O7WfcY7KuD2nbthib4t1UT/X1TPRb+qFf32unsxT/K5Py6tfKajyJsvq77S3QZogXHRPeiaUo50rKY2m07TwQa1qBhpkZcKeFyy6y6zdL17FC5NHMJ67HJ4LAehWrpwysOmEMvHGfT4Oqn7nEbGLQjDue9UK26xLfDfNJYcw2Y6OgIvZlmPQW8x4zWT8hBqRY1h/Xmp/wZ3kwwzR4UmKsMmau7SBTylXdEOjZ2jXSN6hIswuiLLMRf2Xsd0+TqKqLRQjcaAgONwh8Y4802LTfCYRrUGwMPFmYWAE8a9sU70DThTkKdvqcOowqFOpc2MyVJgLE3crS4P3lJmFF6efG17falpGwxEQ+yAIL2uo6FkidEF1jEKt3Uyg1bTEJ4oNcB8FulmlSf2vHYd7LKJ/7EbXZrFnunYtMuUdaw9znffEzl7gHiR9kCECsBk2+50p7qb/XP0dZVKW/G5cH1X2py78c87vLciz0UpTZwMo9yT3DQuUl447hrkxZOgQMmJQ+d4mVfgaJiKVrlwLxQBb8DBL3sARcWbuOxHKYtfk7LJ+DnlsO19UbaacUrABmKe7nchu2LtKDtx8WWCRM/+838RgmTZaBZnzmMM8s5j36kdm6+CGKo02mRs2LVUT/RHb1pX4uJ7hRastFCY+KwTjvZ3lEmPtZYseZ7zBp6z88C22CPsyzy/0eZPiO2sMzFyFe7vUcZ/nvelQt48xguaZNuXDDeLnM+J3bxookYHcuM/RyrXMXy81FKD6C64ebl7Vl0IyUy0HQ+OraCjFI0QXx0q4sMCxpXHDSNIObL6noAnZIgL7KbYA0cdEuxFIEC0QLPiw34hab1sou1HnYMrB0w4brzYxdXQb1kWHSlEZgDSdbtrStOfosd5Q+2xW30jB168I9+pz4Q9baUawxS1F2k3KTQMRjkLV4KcoxiTyMKyzpXPDG5ZT5IdDZGK/5VWdxn/KC8WVzOSIAXDg5gC4uxouZWgTcW0k2KBh3MfPLtBKB2RAwERArsxCCN1v6K/k7RJ+XoU1mu1BnVlMW+C6Wbdz4mLiqe7+Y2XHnvZLST+K1klpzyutiQLQi0JJ4FomURCARSAQSgUSgawgk8epai2R5EoFEIBFIBBKBRGDFIjD38o0VW/WsWCKQCCQCk0FAGIMAbS41cXPN1amTKVHmmggkAkuFQK5qXCqkM59EoEbARoU26RT/M4pYSSR2zUa3TRHsysWasrwQECdjU9jtt9++7Lu2vEqfpU0EEoGFIJDEayHo5bOJwAgIIEmWYTsayX5Yo4j7bVlgoUdTpGX3/5TlhYDVps551SdsEpySCCQC04NAuhqnp62zphNGwLYDdm+30q25ke8wxbJ7tx3S24sGbFdga4zm0ULDpJf3TBYBfcEKaXuL2QYkJRFIBKYHgSRe09PWWdMJI+DYG8dYIV0OaLcfHNLEkmVnc0cm2fTUthg297RvWYijfhy+joC5jwXMcSoOsLZvmQ1THZfjVIRU5IFa9z61NyslwmUndQSsLTYNtSmrPcfsE9XctJZ1zEa2+pKTG2xVYePW5lE77fTy/0QgEegWAjnV6lZ7ZGlWMALOBnVw8UEHHVTOCVVV5ySK9TnssMOqCy+8sARci+dysDsyFYKciQ3zPHE2KQVuo0lpUNR24c94r0CsW59ciuedd15pa/F9iNP5559f+kGTKNt01VmjcdbqZz/72er4448v7YqIaX9p2bxSux999NHpquxWU2dpEoE5EUiL15wQ5Q2JwOIg4Ny8vffeuzr00ENnzoh0+O+LX/ziciKCnfOd9Yc8UdJ27rY7OrHT+rXXXluOraJ47eb9+te/vuzUbwf+/ffff3EKmamMBQFnOSJQ55xzzsxu56xVZ5555swRV6yYRxxxRDmiyQHuxHE5u+66a+W8Wi5q52XuvPPO1ar6XFiWLn2ln9VsLJXIRBOBRGBREEjitSgwZiKJwPAItOO07AZt09c4fgSx4oZqB11TsE3riC0JPOsnpbsIIEznnntucR1HGyst8mwbCe1NLLxwbuguu+xSPrWvZ93ju80337yshkXOHSzvc6eddlrjsPiSWP5KBBKBziKQxKuzTZMFmyYExOg0z6CkjEe1ZHA9IWzixVK6g4C9uriP22cOauMgXUortosFi+WTNSy+O+aYY4rlS9uymH3+858vcV6sZayiH/3oRwceHN8dFLIkiUAiEAgk8Qok8nNJERDD8qMf/ahae+21y2GzrDY33nhjiV1aZ511ltVhzI4puvnmm0vQPEuEwGkWih/+8IflYHUB72Hl8olQxf9Ab1/zP8tW07rlvrgWhMz/oZx9zw217rrrdop4sdp885vfLLFMFgUgFc7gu+mmm0p9HDS8nDcPdSC5uiBXG2ywQTnD7uc//3lZMOEsP31Zf+Ay/vKXv1zBI1a06vMsWkG4uQ8f/ehHF3wsvgj5xS9+UQLtvR/2/+Ku9iPma9999y0u6Ec84hFxe34mAolAxxG461G1dLyMWbwVhIAZ/QUXXFBZpccyY7Z+2223VbfeemtRJK4fd9xx1YYbbrgoh9AiJvIMl9ywn8hNEJzZ4HeIOquDfZmuv/76EjCNEF133XVFob7tbW8rina99dYrrsNLLrmkuvzyy6uNN964WCmQkIsvvri6+uqrq0033bQQNYSNNYMydzD7/e53vxJIbfNUwdebbbZZsYhQ2FZBIn1Wvgmwh1tXDmVGEj7+8Y9Xv/71rwu5FOMEm2uuuab8b1WmWDYutHGsytPWo7a9/qKMw4i+KkDealKrFU899dTS5upnxar6uke7CobXJyywMNmwOELft0WI2D/ETDs/8pGPLJjd4x73KCTOStZPfepT5R4WzWOPPbasapU+0ofgb7HFFstqojIMtnlPIrCSEUiL10pu3Q7W7Qtf+EJRGHvttVchJBTLAQccUB155JHV4YcfXv7+2te+Vlwzi1F8yt3KsKaFaa50WSGe+cxnVso4mxK2man0X/aylxVX0GMe85gSd6NOlOpXv/rVEqtDMRIWPiREusgVkoRwIlAsGFYu2kaCZQO5UmbE7mEPe1j5jGBqZIvlCDG0m/3pp59eNlGl4F3vgiA8SBWSsc0225Qi+V9bn3HGGQWn973vfaWdYYJ0LKZYOWgFIOIzbNsjXe7dc889KwsWZhOkB+my6IHVTh1OPvnk0k7iuaSl/VnzYPHwhz+8EDGTDkTKdVYtm6giY/qGtteGSLxgfP0LLvBD3BBxQfesxcg2YrnbbrsVK+dsZc3vEoFEoFsI5CHZ3WqPFV0aFpATTzyxuEeQCXLSSSeVrRRsBLr++usXa4ig8uc+97kzLpgAxfYJSBEFNax4hpIaxnoVaVJoXDdcR7M9x5pz//vff4ZYyGejjTYqK9MOOeSQsu+WeiFe7fieyGuxPpV5NpK4WPkMmw5CiZRabcm1hohYjYdwXlkHkItXEqvEHWdlJ0E0zjrrrNLu4X7UZ7S5+5AiFkKxUHvsscfMYoR+ZeLSQ3y444bFJYgX0oUozSaXXnpp2RJCOYgYrtX1fmpcvWeffXb53yRDH2oG1Lt32LZC2AaRxqbLUpopiUAisHwQSIvX8mmrZV9S7iTWIZuAEoqF+8VsXmwLQVyaQsly5XGpCCxmCXrPe97TvGXWv6XtZxwS1olI28aW6sg9SGyAGoo57hnX57DkYlz5t9NFNG2TEfFM9q1iqYEZKw5Cu9VWW93pMdYge5fZJsPB0Sw8yJotFT7wgQ8UQsalZ9PYTTbZZA1C00xMvu4Zl3DtcjGGIJpisQ488MByiRtd2fvJsG01iHRJM3Dtl35eSwQSgW4jkMSr2+2zokrniJTmbuzcNciKPYqaSqxZaUpKcLoYmMsuu2zkMw6baS323whBU5BIFjAux2kXODRFfJKgc+7VflZEMU/cg+KiVtVB5kTsn+tbb7116QOusR7Z/X/cFkR5zSYxeYh7uH9Z5sJ6F9fzMxFIBBKBNgJJvNqI5P9LhgALiPiWppWLG8ZGkkgaEWQsMJ0IKB5VrCTj8pnNetBOkyWORWOHHXboSxLa9/ufe5S1hhuUWyxE7A/rRKxci+vT9nnDDTcUHKIt1Z+7DNbaFSnjam6uzkPK4YeshXiGNayJcXzX/GQptRjh9ttvH6ntEf0dd9xxJHc2F6XFEWLwmuVXdu2e1qlmy+TfiUAiMNzyncQpEVgEBChDAeQUEqGMWQnEwYRQtuJnKLO29LvWvqf9v/Q9N+pPPNdOr/k/pf6zn/2sXBJQL37JdhJBssQnnXbaaYu2UKCZ97B/wxq5XWpRd+43RNqPeCtxU1bthYgBs3iAcDULIg9rmPZCZMUCNmOkkJjttttuTouXPKMNh217zyCCfuYSFjxtTsR3qZ9+HIsE5H1mvTIVoVxq4dZlKUxJBBKBbiKQFq9utsuKLJXl9jaDtLpNHBTXoSDqcNs4d9BeVLYXCAW8UCCe85znVH4WWwTSv+hFLypuMO6xL33pSyXYumnxcM0eTg4znoSIkXK+IwJjUUOT4I6zPCyWb3rTmyoLC7SngPgrrriiWBEjaF4gPYtnWLPaFkluR1YkFrK2dat9b7+6wP0Vr3hFv68WfM2WHnaXt2LTykZHOXGLavvot1bmkujbC850yASsjtxnn31KnzvllFMm1veGLG7elghMJQJJvKay2SdTaRuNCj5HsD72sY8V8oUcIGDcTAiCDSdZjbouyASrAlKFXHA17r777kUJ2yaANUywtWvDEIVx1JfV5ZZbbinWmIsuumjJiJd8BcE7BolVyL5WFlUgYn64j8XDIcTu6SeeQ86cUzhsMHq/dMZxjRWRVcnqR5YuMWe2hbBRrB/9nGUXMV9qN6O4SWXwjtmWZByTjnFgmmkmAtOEQG4nMU2tPeG6cslQuJQWa4EVbvZZYg2JWK44FLpfUbfddtsSeG2F26QFuUAUEQSuRVY6BMEKTMpP4L0tMeYTl7aYdYOvjV25RA8++OAlI4Fw4TZmuUJE7VrPgoWksHqxeDZdiO06i8+yFYXNZufaU6v97Lj/57q0Wz1XKrekfoxA2gyVa1HgP8IziaOb9EvuW20O3y233HLccGT6iUAiMCICSbxGBCxvnxwCNqsUI3TCCSdMrhDLMGduMISX9S1cYV2vxute97qyKAJZa6+Q7HrZu1C+q666qpB+GwGnJAKJQLcQyOD6brVHlqYPAtx4LAmCxOPcOq6UlOEQ4Hpi/eg66WKtYS206a0FFqyGLEoWZaQMjwBXKLLd3Lpl+KfzzkQgERg3AnlW47gRzvQXhIDtAwTjc+NYdcatJyhbjNUoO9gvqBDL+OE4y9FmpbHasqvV+f73v19WAooD47YVVC9GDSFznmHK3Ahwg3J3a+u0ds2NV96RCEwCgXQ1TgL1zHNoBCgScUIsH+Kl/G/VHMXi2JmU2RGw+tICgDiiafa7J/uthRYsmRYjWIShzbW1+D+rFFPmRgBeLJzOj0zM5sYr70gEJoFAEq9JoJ55JgKJQCKQCCQCicBUIpAxXlPZ7FnpRCARSAQSgUQgEZgEAkm8JoF65pkIJAKJQCKQCCQCU4lAEq+pbPasdCKQCCQCiUAikAhMAoEkXpNAPfNMBBKBRCARSAQSgalEIInXVDZ7VjoRSAQSgUQgEUgEJoFAEq9JoJ55JgKJQCKQCCQCicBUIpDEayqbPSudCCQCiUAikAgkApNAIInXJFDPPBOBRCARSAQSgURgKhFI4jWVzZ6VTgQSgUQgEUgEEoFJIJDEaxKoZ56JQCKQCCQCiUAiMJUIJPGaymbPSicCiUAikAgkAonAJBBI4jUJ1DPPRCARSAQSgUQgEZhKBJJ4TWWzZ6UTgUQgEUgEEoFEYBIIJPGaBOqZZyKQCCQCiUAikAhMJQJJvKay2bPSiUAikAgkAolAIjAJBJJ4TQL1zDMRSAQSgUQgEUgEphKBJF5T2exZ6UQgEUgEEoFEIBGYBAJJvCaBeuaZCCQCiUAikAgkAlOJQBKvqWz2rHQikAgkAolAIpAITAKBJF6TQD3zTAQSgUQgEUgEEoGpROB/AVJZ3PGQdTnXAAAAAElFTkSuQmCC" } }, "cell_type": "markdown", "metadata": {}, "source": [ "![image.jpeg](attachment:image.jpeg)![image.png](attachment:image.png)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "function lqr(A,B,C,x_init,x_des,T,rho) \n", " n = size(A,1)\n", " m = size(B,2)\n", " p = size(C,1)\n", " q = size(x_init,2)\n", " Atil = [ kron(eye(T), C) zeros(p*T,m*(T-1)) ;\n", " zeros(m*(T-1), n*T) sqrt(rho)*eye(m*(T-1)) ]\n", " btil = zeros(p*T + m*(T-1), q)\n", "\n", " Ctil11 = [ kron(eye(T-1), A) zeros(n*(T-1),n) ] -\n", " [ zeros(n*(T-1), n) eye(n*(T-1)) ]\n", " Ctil12 = kron(eye(T-1), B)\n", " Ctil21 = [eye(n) zeros(n,n*(T-1)); zeros(n,n*(T-1)) eye(n)] \n", " Ctil22 = zeros(2*n,m*(T-1))\n", " Ctil = [Ctil11 Ctil12; Ctil21 Ctil22]\n", " dtil = [zeros(n*(T-1), q); x_init; x_des]\n", " z = cls_solve(Atil,btil,Ctil,dtil)\n", " x = [z[(i-1)*n+1:i*n,:] for i=1:T]\n", " u = [z[n*T+(i-1)*m+1 : n*T+i*m, :] for i=1:T-1]\n", " y = [C*xt for xt in x]\n", " return x, u, y\n", "end;" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "A = [ 0.855 1.161 0.667; \n", " 0.015 1.073 0.053;\n", " -0.084 0.059 1.022 ];\n", "B = [-0.076; -0.139; 0.342 ];\n", "C = [ 0.218 -3.597 -1.683 ];\n", "n = 3; p = 1; m = 1;\n", "x_init = [0.496; -0.745; 1.394];\n", "x_des = zeros(n,1);" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "20\n", "\n", "\n", "40\n", "\n", "\n", "60\n", "\n", "\n", "80\n", "\n", "\n", "100\n", "\n", "\n", "0.0\n", "\n", "\n", "0.1\n", "\n", "\n", "0.2\n", "\n", "\n", "0.3\n", "\n", "\n", "0.4\n", "\n", "\n", "\n" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "T = 100;\n", "yol = zeros(T,1);\n", "Xol = [ x_init zeros(n, T-1) ]; \n", "for k=1:T-1\n", " Xol[:,k+1] = A*Xol[:,k];\n", "end;\n", "yol = C*Xol;\n", "using Plots\n", "plot(1:T, yol', legend = false)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.7738942551160483" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "rho = 0.2;\n", "T = 100;\n", "x, u, y = lqr(A,B,C,x_init,x_des,T,rho) \n", "J_input = norm(u)^2" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "3.782998646332306" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "J_output = norm(y)^2" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "20\n", "\n", "\n", "40\n", "\n", "\n", "60\n", "\n", "\n", "80\n", "\n", "\n", "-0.2\n", "\n", "\n", "-0.1\n", "\n", "\n", "0.0\n", "\n", "\n", "0.1\n", "\n", "\n", "0.2\n", "\n", "\n", "0.3\n", "\n", "\n", "t\n", "\n", "\n", "u_t\n", "\n", "\n", "\n" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ " plot(1:T-1, vcat(u...), legend = false, xlabel=\"t\",\n", "ylabel= \"u_t\")" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "20\n", "\n", "\n", "40\n", "\n", "\n", "60\n", "\n", "\n", "80\n", "\n", "\n", "100\n", "\n", "\n", "0.0\n", "\n", "\n", "0.1\n", "\n", "\n", "0.2\n", "\n", "\n", "0.3\n", "\n", "\n", "0.4\n", "\n", "\n", "t\n", "\n", "\n", "y_t\n", "\n", "\n", "\n" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "plot(1:T, vcat(y...), legend=false, xlabel = \"t\",\n", "ylabel = \"y_t\")" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Julia 0.6.4", "language": "julia", "name": "julia-0.6" }, "language_info": { "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", "version": "0.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }