From sources-request at octave dot org Tue Jan 24 22:22:35 2006 Subject: automatic differentiation From: "Thomas Kasper" To: sources at octave dot org Date: Wed, 25 Jan 2006 00:39:26 +0100 (MET) This is a MIME encapsulated multipart message - please use a MIME-compliant e-mail program to open it. Dies ist eine mehrteilige Nachricht im MIME-Format - bitte verwenden Sie zum Lesen ein MIME-konformes Mailprogramm. --========GMXBoundary327611138145966 Content-Type: text/plain; charset="us-ascii" Content-Transfer-Encoding: 7bit Hi, am I right to assume that currently there is no support for automatic differentiation in Octave? Well then here is what could be a start... I wrote it - you guess - 'in the hope that it will be useful'. If you tell me it might in fact be, I'd be glad to contribute. If not, never mind. Included you find a brief example of how to use it. Regards, Thomas -- 10 GB Mailbox, 100 FreeSMS/Monat http://www.gmx.net/de/go/topmail +++ GMX - die erste Adresse für Mail, Message, More +++ --========GMXBoundary327611138145966 Content-Type: application/gzip; name="ad.tar.gz" Content-Transfer-Encoding: base64 Content-Disposition: attachment; filename="ad.tar.gz" H4sICLDE1kMAA085QkdGYi50YXIA7F15d9s4ks+/5qfAuI9IGUm25GvX6cxGbTuJt33kRU7367c9 q0eRkMQOD4WkbHky+e5bBYCneEmWZGcbeN0xhaNQKBTwqwJBQNV3nq077EI4OjjAv+2jg132u72/ z/6K8Ky9u3d4cLDf6XQgX3v38HD3GTlYO2cQpp6vuoQ888eOpXpF+ahbkPytBlXfsdatAQv0f/vw qA393253ZP9vJLD+H7mqrmqO17LWUgd28CHv76z+hyG/H/X/4R70/97u4cEzsrsWblLhL97/331H TpzJvWuMxj6pndRJB/qK3DBhkF9Ub0LdBvmJC+cT+/l6ZM1aNvX/oXz3HfwHeQ2PTFwHlMgi8Dh0 KSWeM/TvVJe+JPfOlGiqTVyqG57vGoOpT4nhE9XWdxyXWI5uDO+RDsRNbZ26xB9T4lPX8ogzZD/e Xn0kb6lNXdUk76cD09DIhaFR26MEmJxgjDemOhkwOljiDfLQEzyQNw4QVn3DsV8SakC6S26hM+E3 6QR1CIIN4rhIpKb6yLlLnAmWqwO798RU/ahoK6f5USt1YtiM9tiZQIvGQBLaeGeYJhlQAvo0nJoN JAGZyW/nN++uP96Q7tXv5Lfuhw/dq5vfX0JmkDyk0lvKSRnWxDSAMrTLVW3/HthHCpdnH07eQZHu z+cX5ze/QyPIm/Obq7Nej7y5/kC65H33w835yceL7gfy/uOH99e9sxYhPYpsUSRQIOIh6yUQo059 1TC9oOG/Q8d6wJ2pk7F6S6GDNWrcAm8q0UCjyjsPiaimY49YMyFzJMiXxBgS2/Eb5M41QF98Z75b sXjUsw1ybmutBjloQy7V/mSC7Hs+5Acab4wh0H9jOg7o8s+O52P2yy4S2O2027tNmG/ahHzsdRWM a75ogv7NDHvo4DNGvdbp0L+f0KFNvlyqk4kBTL+Z2hrW/JV8+UqCGZTUXt+q7pfZ1zoWc0BXTEfV Pehiw/SbwJMFpUGJXmuOTr9gia9MwCqjYFDbJ4IAq5baOgmrZjEeparpObyoogwFEyB9/xYk/Crk pDaqKwpBOdZs1R1B1X97Rdp1iCKge+qIktp2xPWovl1/CWlQoTHEcr41AWK1Zpu0doj32YXJoU2a kLE1qx3XSet/Sader7eeYyEXcnrGvyim6rMG1IKxkBPiBfnWTMTpGMn+tF5gJTXHppCh3SBuvUGO Wa6oLS8VYChoo/LYc+X/xxDhv2fYj4b/e3P4D2agxP8NBIn/Ev8l/j8E/0P4hwl0QfiHEsvCPxbN hX/GSVX8F5mz8f/B8M+pS/h/qiHCf199NPzfP5L+/yMFif8S/yX+r8b/hxl0QQMASixrAGDRXAOA cVLVABCZCwwA+PfvS2A/Jyyx/ymHEP81x/7zkfC/fXQ47//vS/zfRJD4L/Ff4v9K8B9n0MXwH0ss if+saB7+c04q4n+QOYX/HMNFosBwpGV4mgNKSGcc8uucoED3ILs+SxKTsP5UQwz/vfFj+f+dDPzf k/i/iSDxX+K/xP8V4b83XhT/vfHS+A9F8/EfOamM/zxzCv+5s44ef7mnL0hke/qeYUPirC7d/acY QvwfGrb+SPi/t5+x/29X4v8mgsR/if8S/1eC/ziDluJ/CNjcAsAyS1oArChjeGKqSNlxQSw6sBZZ Bf8DQH3fIP/6pzAMOIsVDYMgc8owCMqhXv5EOrwcAr7IL6wAano0nvXVqyAvY+qfVQrsxQqIVsyX YVl4dIxRaU1UCyH+G5Y6eqz1/735/X8Hcv/fRoLEf4n/Ev9Xgv84gy7m/2OJJdGfFc3z/zknFWE+ yJy5/i8Sk559EClW+aUT/w0Hhv+godSfTtYE/2X4f7A39/1f5+hArv9vJEj8l/gv8X8R/A/3/b8x TMrAX0yf4H9i3h7MJ6DexKMeUzGQgzr1Yfz40HRQ9iF1AeYNxnWgQw42/NbQKemeNgMAVU0DOljs FORVfI2ZFCqxp9YAFJKrgOZYlhOtL3hc4ZgO4xgwjU/UvEdBguJhF5ITFRQxQZ3UnmsmVd3n9a+Y 0aUgKSg6cQ1QAZgmfZpplCTsEmZR2AZI3fDwOdNGCevTLB0NEnD2yXZgAG2jTYGqoevEpndgkcB4 oh6BYUuT9syrV2SX2zMw6Caqr43BpMGPFbYbnJx49pmZkshFZ5MgE3/MyGM6oyAPf8zI4xl2WBl7 zMijOV6Qhz9m5PHVkA5/zKTjR3Ry2qXGGFJzOVJjLKm5PKkxptRcrvD1TkwE41wZjGNCyMkFlYxj YsjJpcarVPPrVOOVqvm1qvFq1fx6cXNL1Ab2nJHLhXEQ5BLPGbnQhA5yiecszgZRNw1yesmbWqH4 2WNGHlwdCzKJ5yDXLHAZYMzWdusvCQCxDWMd5gZALVeFGcBDrCbglRjqwKQeKxYO5W3n03Z89c7w 6OcppODQhhrZfLJdzxuluSMzdzjmjsHcgZc72nKHWP64yh9L+eMnf8zkj5P8sVEwHgqGQIHW5yt6 vnLnK3SuEudqbr66BimRrjFtonqocIS6LkAT1G5DDkMX1lrck5UOaXEI13+xa5/Q/u/DQ+n/bSJI /0/6f9L/W8n6L86gi63/Yokl139Z0bz1X85JxfXfIHPm+q9IrLr/O8gu939/KyHEfzQin87+7/au 3P+1kSDxX+K/xP+V4D/f6LwI/mOJJfGfFc3Df85JRfwPMj9g/7cgkb3/m28Ol/u/n2YI8R+Xgx4L //cyvv/qSPzfRJD4L/Ff4v9K8B9n0MXwH0ssif+saB7+c04q4n+QOdP/F4kC2fmZMHMHwRUbCPLY l6ccIv//s+s/Gv7Pff/V7hxJ/N9EkPgv8V/i/2r8f3ZM6kL+P5RY1v/Horn+P+Okqv8vMmfiv0hM 4H9tt3WAp8JxE0Ci/zcdovNfB+u6/qUU//HOn/T+77Zc/99IkPgv8V/i/2rOfx0sev3LYOnbXwYF l78MvOpnv/K8JOMTb8MLNgbwYoae+Mab/ER2xf48cdrrIHbPCyfxN8Oj1gSUtGaEGwW4TVA7biDB OpRrxn/zspyP8Otu0BLPm1ps1jA8pggT1fUN4A6mC+MW+v8WbApopG6AEjLTAUsJBWTNUGeG2BF/ R9l+eMOGsUw9PkpZrdy4STVDnWLbatExeGj4QM56+F4kXCmxJmK7YqYllLFFAi0ioJ+2iGJdkdj9 uB7jKH7+22PhP4R5/Jf+/0aCxH+J/xL/V3X+28LHvy1/+lvR4W+LnP328KPfct/815rs9heAfun1 P9UQw/91Lf+Xn//WSd//0jnqyP3/GwkS/yX+S/xfEf7nLv9bokzq9DcosbQB4POz37A0DEr4Y6m2 Roljm/cgNcf95JE/YW7DEQJOO+UyZpN8vuFQ/aWByJv5zoBdGZO4+4W71s3UpgHSJG356uCRQ4j/ dPZY5790YAzO47+8/20jQeK/xH+J/yvBf5hBF/P/ocCS8I8l82CcsVERxkXeTBjnaQLBqy4JJDB7 bmlbIvaTCyH+m866jn+tsP6/O7/+L/f/byRI/Jf4L/F/JfgPM+hi+A8FlsR/LJmH/4yNivgv8mas //OLXxNb/LhZwIsIs0CaAd96iH///3j+//z6/57c/7eRIPFf4r/E/1V9/7/w5//Lf/1f9PH/It/+ P/zT/6Iv/+WH/088xL//fzT838vY/yff/28kSPyX+C/xf1Xf/y/8+f/yX/8Xffy/yLf/RZ/+z335 l3iDT/4OJOXXf992iPz/qfVo+//a6fN/Okf78v6XjQSJ/xL/Jf4/5P4Xhtj3X78G7vfUSloBUdnZ QoUb3BbQDavqZbJQetn1BCiZu56API0aoNhWPQXX6dtYxAWvGBlYAzNeEFPbwReAQGLq2slP/cRf Tp2vLLB6W/E75rg9IRLQzujUpc3wsBB9/4/nNK7HAij1/w/T7//3cElA4v8GgsR/if8S/x/i/4ef /7Nzbhf6/h+n3GVPAGBlc88A4IfuVj0FIMidWgcQTj/bCcAPAkr6/00gVuEVQUA++yWBXBV41BDh //ouACjH//T+/73Ovnz/v5Eg8V/iv8T/leD/4uf/qw+4AEAtvAFAXegKADXnDoB5A6BN/p48+7eS AVB4P4A0AB4zRPi/vgsASvH/4GDe/9+X+L+JIPFf4r/E/9Wc/7f4BQDqA24AUAuvAFAXugNAzbkE IG4A1ObO/a+E/ckbBCT2P7EA+H+pfqJDw6Rrq6MM//fj9/8ctfH9/95eW+L/JsKvZx9659dXMNY6 rf9s7SuXv1yf3Lw5vziDmJ2p5+6YjqaaOwPD3rE+OZqPiqK8u+7d3Pz+HvMYh/9x2JxoTZhvp7Pm yJ7CFBrcIM8mZjJRYfLFGa4FxZtY3oMsFigd8aYAK2ggMABhU31T1TTqeXimiWV4CLTK+fsuYOPJ TZIj0xjQGdV2gCgU3vm+JlpS3/GADEZDXMBofaeALYssy9RlkiVvDLCUwZC1A6MMGVB1nTRPr677 PcDj3hmrXpgTHgdB24mQwpuoLuAvTv3K6dmbHh6hcoqtOT/tv+v2+icX3V5PUW66H96e3Yhp37AN H+WsXJxf/YIlDA+jWdT1z/8NEc5tk8ew/E1n4rUc5d1Z9xS4jSWPFeX0+uayFztRHgQVHhQSfGkm ntjHQ8Hpc+KJbSgmf4h3viQ4YkakqlERNSqj8kLRfaQBzfEcKe6uRM5ujNg4Ri2iYf8pHvlN538o W0xglho0gp2AHrzyFk943i08QhT1pxM8O+cHkNwx+aGlaUTZ+r4Wjpc6+b6GvVQnTY18/xNpOuT7 14rSmowd+/4YDB1Tgf+PIRfvMMzPOgkQOorDdOin+hxpjAxoil+YV3Qco8GpkVgVypZpk6ZHmkPk KM6OGAqK+JvJ1lbT+qQbLkRwZYcYbUJeoC4FccAgj22Bkr8AQcXyvqba2EmOOX5zsE8jPqZ2wEn4 dKxsuRZjOapjjrvjH17Fkn+oR2WY1kbJlyyx6VrJloQMaCZVbYX9G1YMTeTN/AsYHDAzfQAdujxb Yx0l+N/e39tL4X/7YF+e/7ORcM3QipzNfHDJcIwiJnWnPgjDB1ftjeOCr6WTU2M4pC44JwYfybXu aV1pLh8UBWFOP2aQqyQSmoTqeNL3HFgLbzowV9Ep9KYG99MBGsbg2CFW4soC4SDcFOCKRNPYC6UZ qALVkGIcj4FiISIrSpxD3pBwWi1qeHMBM4NxhC46ScqDH3TuuAaYLUAQybB1BZuC//1cN4BRXxs/ J0PXsQJZAKkRPd6xVDTlDE+jQMymzhRbWgMPEYnoDnrdnCVcozEBK/V7mC4V5i0eK1EjuM86ABaG 4KBaIDHf5czBf2PnDgUMY4Z0TxsEXFaQf5/OVJz/AUUJrlB4uJhxz6RJtsfUnIQWzDZnHOUvOp0d 1z4Bv5C6CurOyBO8KANKIB46huqtVgs508ZU+0TUkYoiI8/tqRV41s+xtzPEgYs8yNTdmNrsHHln OvgrTP7P2Pwf2oGatp46iuf/o6ODvaO0/3ck939tJuy8UJZe/1VWsvS7knXfhy/6VlnxnWvwoou9 D1zpfdgy70JrvMrDl3cfsrb7sIXd1KouUV6g+23YmjnVAWtCH3c7ivyJ+yazf8SihCcP2UVia5yR PJs4d5kJzFjAFOU7nbIzWE+6vZt+9wrd96vTi7P+9ft2DZphqxatM/eWZfj5/Or6fR/8nR67eAQw jBPsE3y7QX5skERsAG4/1l8yGiL6VjWnlJEGF/623dKgN2itXsRMZxFmwmpT7Agmi5jpzDMD6gtj 4bR/c12LGOD7xJFGf2DYzoTRxuTabbsBZFLl+z9fX19AM7AVDqcAFh5rAKk5WLY/4uw1wrX2BvF5 zi9ieSFDKgQJ3okGRVxxdhwCzLRmtXqU5Wtm1VwujIlGnIOy6ttl1d+hKJCDoPp5qXYvLphwzrpX QLFXE5XGhcZqMXMTaE4C/ZyTMMorMcqrw2YlQt67vd75W9Q9XPEKOAahBvGkVrEz4xq8kO4SHDYq 2OEjGy8wYhqX7oWYetfm5X9yfXXSvem/7f/2rntz9uvZh1g7ICHRiBRrjQRPq22WYN6md+laa7rT 11S/70xCrU40MRq6kOKqfbyzqYFi0hP6t/OCUNdFTAF3BvwGj8AEDHYWepScM20MRtcLgIG+NQtr 98grZWsbgY9a8FO4Og4UR6fEm04mjuuD9wSQO2syhvSw773tl0pmFeDqWOgMGfYImsWrAIgGKNER v9GtRegGL2EAT84wfqdUjKaSmMl0Qx0h59bUTM6KLPlHYjWUra0tEHVW4m1d+SJeyt22DK/vaaqp un30g2rilizRQ7VbbAB7JRYn0Qe7x+cXcgVLonPJtU6y646PLXVkaH3Ngf7r+0BTYXd61XbxEq6r 0y5YGfc13bD6t8yzxBdyVst17jzW1+0Wv+1rZwd488Yq2Djoa3pgDdiO3WTNBCsFEtkygcZgnYnZ xzPJ00jA3yvevlTi8xk1mTxBduI6rpaoKsXWbcuk9sgf1/AFJCisAdbjDFWWDVLfnVK8Q4ur4dRW 3fvozSoqANoCMP4+sjlkahn21IvNxqxr2AD7KMZXHu7VY9xnDSVoEacOvTxjQozHsAHD2QyZAQfa 9iYOGqTr4SesIGRJMJLkY4zrD76h2svzwZUbaaOG44p7TL+/lM1CwGrIwzyrOMVkjSym9z8K3RpS iKtt47L/dlA2eDPN7qZjE1QfRzeNeFqCLXyFvRuxhXfW5aBDQsbMAH5gLwM/QEZwEmh8qOuDe8Jn llDrQ2g9Pe2/7fV7HFnXYCx5VaAmw4LyIjAhMWNKiTNMarxVDVzRq+ckedNBOkmY8JnFwjRWLpTT 5ceLtcmpspiypZSI1UU0iRA4znokFphe6zlJunGbTgqlklEsTGPlEprVe1sqsRXbwiXaFRnI7CGh XVy3ApYzdGs+KdStKGletzLS5nUrklRaTCDx5YVUWbVgVhCSgQL8OSWcOJeBBOqpWNHAWNtOz3+N 2ubNNQ5UhjXOC1vnLaYDXqKBWSYWxwCP7Ij2BXnnZn4RH1QngKszb5bF3j4XSBTqbZBaE6sPqq5H gg3ZINw8rm13T4WRBLbyzjE8d2CwQWsGlHifp7iMhS8z0MLywH6dAj/bIZGvlVyReF9EPZiMne/B 99e/RTPffA9OnDsx63mJac+rgg9Vew8GMNQTm/XKOrFCz6CWp+niZ2xg3bJ+aif7qVi0cSFFok3G Zou2YHCgaMPBscaxIaTgpRYPqgyRMin/oWzhQ9rwI/+FKQlZYi4xzIRwoRNMZ1TzWsHsIOLjw4gc ZxISol6I0oKDKd53yR4vGkxnFyUTIjg+G5sTuZ8F1a2r62OuMZtJuXESVdtkGUXNC/dAUpZRH6Tj s3qhZFID/jY1r/GqVjq1zUk+NAv5lJeqEuSfGhfhNFhx9kvKM94TZTMgz1E8HhLT4PrHwxqnwuzx gFMTLson5yycmvBy8YUHRfbElI6PugL8xIHj4Kan0P7gLmLOinFEMzdDgjhaLDqdUFtHT5SvoeKS HlaRWNiNCKeiE+Sg+zSQv81XBJHI3MpqRCgjKSJ26xh6sJUC+spjS4I1jOU++PlV74Y1Tbgrk75Z sDYbuASmD9OGqP9lNhVagQoto0I/l1Ohn8uojCrwMirlZVRBLqNSudgVeLETvGSSAb+rnA5kKmMH fLRyOpCpjA4M9nI66OSV0AHALKeDtkMJHZjfyulw9CtWQHNlTePmwEpaxyfwCkMjQPhiZTLnS85R hhEvDNdR+YjPp0LLqMRGfC4VGPElVEYVeBmV8jKqIJdRqVzsCrzYCV5KR3wuHRzxJezER3wuHRzx JXTiwyKXDg6LEjrxMZFLB8dECZ34gMilg6OhTAHNlTUtNeLz9Zm7Q+XUKrVRDPpsfYpeJqNKeXYl hI3mDxIjJd7nVi+M6/VJc0R7oDkSLjUyeyQwehY3SCI6tJROkcBCOswmKaZTZJSEdEbl/BSZJRGd cvkUGSYhHTvJz+KmSXzhuJSlQuMkvsxcSqkQw+PL/qWUCiE8/pKglFIhgoeUOH6XqGSJkbJQA8vM lIXaWGaoRAMmsFRKlMvMKJtlqwR+ZgVjpYAOLaUTmw3y6cBsUEZnVIWfUTk/oyryGZXLx67Cj53k p3Q2yKeEs0EZS/HZIJ8SzgZllOKDJZ8SDpYySvGRkk8JR0oZpfgwyaeEY6RUJc0VNjA1GxToOLdh KtCr1lIxIeTo10JmTISCmbNViR2TWXp+5wHfpJW98+Cyf8kW/qyyV+pWYhHWqvBSHUtYifW723bu u7zkSz9cm+PbnfhZxeTf/2brdZpjTi1bROOBh9XeBOLS39b8u3urzhYFt4JdDi+ytl5BreEOJ771 isy9OkztnMtbMSzaUZBgKr3tQmy6EH1FarxD47sn5pKiTRdhUqCvBaXTWZKvyTHv5dtyheHqEqwV xxWn/HV5htLkvwBOK007Q2naK1IaK7ZxgSlNp0hp2ptRGqt0N0XQXxlKM58UKk2UlKs0+Vky9u1E s0zW5oolJxjyf+1dfVvbxrLvv/hTbDlPEznBxjYQWrvhhEMgh5yE5AJp06ft8RW2bKvYkmvJYDeX +9nvzOyLVtLqxUDS3ufgpw22tDua3Z2Znd35zSohLNlxBa5bCN9cEXuQba/0okpOY9IHQlIck4hM z0QTEEUwJrqr04yhDCfV2BPyhTASwziKbqWGlAuY6AIiRbGauBoXsYR0aSYpC7pzS2tUXsImCr1z j+iWlIQ1byNhrRwJa95Swlp/hoQVPrSMvDEucEnbV01czRQ4Cm3nSNxYggc+u8jhShFj+BNzqNTS ZI/9D2aueYi6F7L46JE2UAnJ0qXCMNWWHafN5Dil4QUG4NUvvyDyqlkeeZVvXGLDFQ124nLmaJ9i uZzZi+8zfObZi6/lNAfxDuOsz1H6/HKbcW7KcU5b/HtA1OWPa2xgonFNXM4c1/en715+ODjPG1q+ VL2v0c3xJtaidLu4O1QaMIGJGhnScSfAhBnDYlLkAqRWorej8UrfKRyyHNsrhuwerW/OFK0NW9zN wP2AFQYuw3zf48C17jxwaQOavpM5cP/14d358eFJkbJ9ZlNqzxcxwB2syVZTllWUEB6m+Tm3UzOk sep4pbo6GjDDrcwRK3BxvqSHk+PglPStY6hXk2dtev4AD9o01CfU7JpVA+nZtDhTacys8XpVIGUT ldM4WfONqsl3TywG8M5KzruFq5aBkfS9LgvST1lxafBFAfppY1fSWVQqlj9F3bP6FK4itXkq+eiV IMqrq+FaCUQmy5is7gujnBoUg0ksMbYCSJs7g0WA5c+7GEigh+/T48sELsf8FQMPErB8Cwy53rHR 4CSvF4xMvspJAPOXmLFEz9xm0srMouXOi8iiHfvDdeCG8mdvP7Ct1MBSOXiSdfthTOtY8npqGFcH P8uqReDn2DM8PLwgjoHOgEDLWmYIdIxoKSS0rJCNhI5oypMoPlpDFNKssJ8S3E7ycAuTGogYYiT7 goo0VEmiaW3IPQNnUnwIzgrhmaZhIyJBVahalHWkJD06WCfVqnZbnoqQVhooiudsIEKhO5SejDpR I7W1ETkkWf5LvqokhiolISVmIyhp9tDxBo2RvpYqsHVELJ3kLsc2MYaGISwewRzrdLXKkMW3lwq6 3/S02CiYXL3E5dQIJFIWJnfDCMqIGUEExTNWRwgqKk4RlTwkgqRC8MBcKnnoQEllWMhLHjZQUSns lzxkoKTixXgxksnFBWqhzyJ2clGBWny0iE4uZE7S4XuFuXRy0XKSDt95yRecAhSfEh65aVpErRRj aleoiFoujk+jxj30fFEYp2uaQHz8XhkMXzYVp4iKpq+ZVEBfC6gMS/AyLORlWKJfhoX94pXgxYvx UqivmXRQXwvY0fU1kw7qawEdXTky6aBmFNDR1SKTDupEkfiVIjQuQSmh+dliyDW/mFopzoTmF1PT NT+PGmq+WahWguupiclkkgrAeivUzRahnCZw36+IQnoJEdH58fjl4UkZNoSDqdqRTp64o2OUhBTp ORS3d5RSVJ2yVEtAOOMOVCmqJdIs4g5VOaor9OuwdL+WSMGIO1wpqqs7XgbsWVl2y6RnJByyUnTL 5DIkHLRSdMukNSQctnKCWy79Iu3AlaW+EuPKoStLvUyCRtrBKyd66XSNXIcvXmaF5I1cB7Ak1XQq R65DWJJqOrEj10EsS3WFfh2W7td00keuA2mgWmh/CunqmSArOJaFdPW8kBUczUK6ehLFCo5nIV09 mWIVR7RYIVagbE4lKXRQy1NfifNEislqjmsZ6nrCyT04silHIM8kl8tCuQ2twpaXaXLM8S2kuIoD XII9zzidpXNw1Hft1NsYBH6ohwyj7droMNLMXXi13Z69DV9RgX18ucOVwv+JX9GZqyUw9ALmbI5p 8Y3fSgxLYD5VGYutGXDzBUADPMoZt+jxrxWBzuXxfk9ldo920B8TzMrd+9s9IQZpUNFo0/PU4270 fjCE6BP9EA+PVAqwHHE2FfKCc8MZa9XVtbt3hBnbIR5X9Dy9J1JxoEePTLkOUe8UQQkqis0MgShg PZPLJAhGMKozz7lcy4pIiKjF5qYh+lshNHosoJQ4Hr2qHy+sqWR0fHESq2o2IvdkQyp04njPDvAQ cfh/ii9dEaca1/bwtSX0drLQFwufYSXTBGDgyKgRBjCssPl3tV9plJ6Mz5ewZAXhdRA9M2wikkJz 9F1qq5SUinZu9d1Egq+zPqs83NucsrnqnJJnS6WBw4J0bP1Kc0vSkNR080GnwSluV7emucRjtome ZHqUNFF3hsXdj7zFUolf0dtiwEuaJiE49K4Znkqc4eaUxxyUkMw45t4snSvge4WUgr2LZ3jlnRBp SN1K2UCYTr7ORe6t5c16BCZdk+CegiRmekNPKiWVt0g7ZHctzvTXGVzfP9MRJ9n8p1Jq00yXTVlJ NZxP05lAj5RSlAJ6cKVgFu0xPK2q37Q3UKsmJvC/kMFWszxYkqwpXrAbKVq2TdandFGNT+krIPVX mKtXBuoXrVbK+Z6xCd5aXRtxOjA7CTgZqHwAzf6XEfsM/qIZqOiR9++Z8D3av4Sglz7vvED88iWE Xp+zoa2PjMTTJ5+zp1g35zhtTiCdKhA7U/v+Bk7tr3+ZsbvdmkFwyA3MKjBqbYizX0mTidCF/rvC t1nFMLqSdFMIToa1a2JFzTaJl0BksSgFJZdmK06TBBB5tFTFYqORnZjUzF7R4LOjF/Bki1Um7vdv +P7RQeyNdKu/WPALepmZbyBMvgcwE6/MwcSmE5rlRqH+mr1VlO8zmEzDW/juamGMSGr9pXxloNDU AeL1cDThihd/IoaXLMUef0dbxk18Q1uDG42/wvSTfvOftAoSlVzCEUKhAqoBoyXg+NpeBlpHO302 8mfuH76HkI3lhvImYKSG4OSN8eXn/oDB+nmTb4QxfMW5fBMgjpHwm0Dxm/jeOovzhi+x2xN8Wk1K motutKMb0pUpTI/U+0Jtq/Av2BX35aEkoL2chSSMJXFaOd+zSJwZmrjIQcLJkupiJQo/fJDRG/ki vJTU8xuxmIWqpL2wLlVP3TNX1d4xl6qq7pmr0qvT0hFYPywRds/UZgH4KRFjzyPh5JDIi43FQunZ JPLQO7G4eQ6JUn0xzOuLPFxOLCI+LBcEzyHC4TdlIt45RDjWpkx4O4dItDG8MpYmHsXOJpKLPIku RO7mbdA3qeh0LqmS7Yr2Em6Jq0nFnVOis1KoWRPkBEsFceVkRXIVKp1Oh73xYQXAfrBnrn0B81Sb PXnCb0z8vtNmB0+fqiuHXp/fhgnrq7/6x+5vOgsbV3Obn+0ZDfjs7uzg3+buToN+N7e36a/4fNVs bD3b2dluPtttftVotraau1+xnc/GkfaZw/Q4Y+yrcORP7CCvnDPLuf3/9aONv93viq/1yb0+Awf4 GR9v0/i3tne35Pg3dhtbMP7bO7s7X7HGvXKR8fkPH/9mpwJrSDHu6HyP/GvcbIXmsv2XNWcBnnsA 3ne7Qm9jmtnucBTKTdkT5zoEx9wNcTWNLjpUDPzxlYOhnzEsn6Fnnd/ndK9SGcy9HpVagud+ZC1w 6YBf/3BmPni8uxvcQV9y175R33lCXv5TtrBa+Idf2cLdiIX1rLq5sHZFhRZWoFtYeBv/tKDoDhZt aeW2eDmdqCjV1EptY6la69tvt3Y5B3B767tG47snVKfGdr9tNbeeiOc1v/2uheWgVr1ej5Y0T9m3 23Rjh0ptbX/HOUFqjWffNZ4k2rBj4k1vkeSUF3+GxbcbDeKQHv/vLdqm2UW2mw7vKno8r7Ar6cs+ bNWfNXa+pTZR/U7F8fpyjLTR+hnhyqEzhcXBwHbHTv9XIOTxkbcGPXDgFw1w9/0xjicvAQWITVz8 Lho0j4PIcCp4D9ZkzUaD+mops8wFLZqBPTe01Cqr74xDG0qxmrWsv2a/sGUdU2TZ3xj+dAP22u75 F7SSgDt4gZZYVJcYwI0jpCECAH3Pn02wCfjHojva2pbf/f65bBB+VKMacgF5MXPsS/4DOs0dVOgv tDHeh6Ay4DWMwY0NYA04hWU5KgjYm5AFsNbF3gWNs+chmB7cReCvowf/w+Uqg2+AHUGDrl0g4l85 s7Fv9/GV6+pl9RN7OgW7BP0/8ufDEeodaXE4cibYFXYvnONaG0mBSuNqE+9JhW9XFGsdYnd3scsC IDXuYyzXHs9Q3Wnw+mAHHGDen9lDp86sD2+P3u8f/AtWt7/P3RlQphGe2jNY8VXgUV3+vfubGJ2A 9x+1yWcD1+vDOt8PoXGO6BJhTeywsmhA4Z9B3ztQhf6TX/Fv61eighHmoRNWlHzqcvn46DGXy6ZT +65aqeCodv1BF7pr5i/cCW5GdMFUzUnEpTAcMRQ7oj73+r7qcljNq4563BuDWXsM/fVnm+57+cD8 71/VUOvqvd5nekb+/A+u3tYzOf+3mts4/7d2tnYf5v8v8YHlzYE/XXI9tw6qrAVjxc6pM9i/7ACs ywb7nnfOJf18MZws6p4T7lUq52icQKFAfMjaDGYOKLM/CME7cDps6c/B/nlgIvouuA7uxTx0wFlA zd0EYwErJ3ewrMAF3H+fkV0CT2ISCAPGXp18YK8cD5yLMXs/vxiDhXzj9sgOAXNTvBKMwPJcLCtY /Aiffiaezo58oEpmtMMcF+7PGGgz2dyWfICgBsu/WcWyQ2QYrOoUK1WByyUbg51Q9erpBkftUpZ1 5E+hGSOgBg0jww2GFIRnAKvyCpRkPx6f//Pdh3O2f/IT+3H/9HT/5PynDpRECx4y58rhdFywzy6Q hcbMbC9cAs+Vt4enB/+E8vv/OH5zfP4TsM2Ojs9PDs/OME7B9tn7/dPz44MPb/ZP2fsPp+/fnR3W GTtzkCGnktOdAxoO6LW+E8JkF0BTf4LBEzPBCJbGMIg9x71ycPrpgcAUj1HFHvvekBrGZzHRbx2c Z2kP73oGniPOBqnRq0Sjt8GOvV59g+00oYjtXYJbyc5CKAwEjtwBED8a+z4I6T98mJyg+Nv9SqPV bDZqza1Gk7EPZ/usQlEe1+uN532HrUuDN1qPLort+c0+DJRX64/h7p5290ou/PEyRgqOTw677w7O 9384xOjKO9xXOE3tK9Ck6g4wuISHOxzjCT9n3YM3+2dnisYPXXHr6MPJwfnxu5MzOgro5f75vmGf Yl1+Xde/w3P+RuH5uxIlQuTQYIhEBjd6YztAr0p0AW7EYyAtuWnSbhumfZhbw9nc6VQS4Y90ZXBq umAG3Ct45pUTcB8S2aCddTYBhwk8D3y4ZItiGb0RmG88bSeGg8VY5DqqkDNR8YO+D1RB8Fgwn059 cDg4dr5G7mJfBbaC9Y4MEIFyz3sUJnwEzyefS7nFwEeqCclgginacrXB0L/lnRBieGKB8dsrsdmP EEIM3qjYXBS7SXfQlUJ0CFcZ/ejnKHLsax6/7Fsn77pnYBbODquVNUsMjvDiurh3RWcMUu34K7XP qOhb6iLLRU7ccCnSDawEKxtJ3ujsQSlKcbqrk9Kjfaozbt/UnGameIk6eEMoQ7dnB2E6ltfIaXEZ 4g0ZXr9RcaGUdE1nLshUEPZB0EAyHXvyiPlAh3RxOuvadHIkuLVLkPhFVegHiRViNqh6F81xF8GR FlZdJw7XN0hH6fmbm4wKMgs6BtYh4yqLmLaC6gYvM/DHY/8a/GLKR8ERxX2AGljjmX+N/jRYCtta dN3qOVSQxiCLi4xnRXxBn6QQEOkO6vvQKwNmUDsBmcDYHe+QNDWEKtjjDoV0oYdxXofWebhS5IHE Or/E6zBLBV2xyHP26nT/Zfdj9+z8lGD2TNDDACQuu6k0xQkTVV5TlVSNqMeyKr78aK7Z1x/GTael DDzjFm0OM/3IRmMozdp/fxM8hu72JvUeqAxUqeqhS9k3iWFIj0AwvwhgCMQIaB35iKGeyAA23cBe /z41DntZo6QPEIhPbwSGEkj+3PiVOlyEtxFD+Nh63I4CvWQnMKILK05ngTFdoF8fzHyU8aom+Yl1 PSdVf9xm4h48YS0tNiERR7mL6PLOW1MPt0Kl6SKt4+8s5KHpWAutMAW61hgCK2fPx6FsmhzabwL0 sXFaAz+T2ohOIzpd3/RQgaCPuKZV+dBSm0XHSctqDuor/vkX8PkXYVeOMB9OAmQ0by0qdjD0KOeH ZQlM8mae0OiF43PubBRUVxGoDHn6VFkD03ftoCd7AQZuCSPvPQ6Z50CXo4sLthhc+OuRzy49/zrA yQYqwPLhAkzahPvNICETyccFkOfMRP3BM7OQEyE0LJFKwMcEZhMwlb0uR8YwfTywtZ3KmjIabtB1 JlNYQRBy7ROSSfSEAKhpl9ptqtJF6y6JC7ox0UCkoOTcxMQN50M9NmVblXySRV3TJRscfSnPGj4J 1h3fQEeDSxriQ3gRlX1gwaRBVOKWbO2mEmmVplMwOrAwdHAph2sQDkEeAE1Y9jlB2AYX8ZKWrFMf OADnF1crdq+HQBVuOQMxoPMwAM9GEgUBaaOg4OqXr3JodcUxEfaMIFdenzl24I7xZWnwQPgfvWtW n7CBi1427iYKcmKXj7AcMMMjiy/fvFS+aCDLwWi12djnKJopzqS9ZZ39CPKIDfTZcG7jWhIWWqgl oDiO16OVHGm9pHI9giWobB22AHcCaQBwJGYXLjRjtuS7rMHf48YyhgbSBzvg86g0kMn51FKWDyU2 SMyokQVmUSFdR7iKtHi+HGlx81e8AqqrMOF6fZZrTXAXkRBDHPYE4iQu1GGZGziW/HXhDF2Pyxcj L76O4tKde+7vvEXiMnoAeFeZvHUL51pORaqKZOxGfSMFFvWhTFQkB770BM2QMbswMHoOf/FO7Zs7 tS8m9dt2av/+O/W13qfqbSRqFRxN0UqTKKNh5vTC8XJdEdeTrAqdtm8CaHEgp3UNlKbZOA0SW3pq 1p2lPF8aV7LoOAfuH07Xv+RTLBiQjx8/sqPjj28P6duFM7KvXLDorqesTsg9FL4E5CYcjQRZHrnA oBAjXlUPIFuN9mzJAh8m10oenpE0LuH2YQRrLBBp2YjrKG2C9cUUB5cCaaKSGxneDEXRamTd7tHt psw/IS3rRyr2PdviGoZ01AmnWCtK9pU6lnY/8TG/EfA6ecuKrUHjc/vEHuJK1h/7XjesSs44b6iD irk9ygpP6HmyhWEnbkKgIXsydxk3MwRcs85HQiSkVEGtKQ91jzU0a8J4i+o44Xl964BQoD/wKt4s vgjv+3OYMPfgEU+xp5K5kXGeegaeGvfAU+9WPCVIMs6fN8OVgQC3tgWLGYmYsfq8ZDSOBjcXZSQb 72rQle5vG5pqR7TFYmBiLy+c7mSOqhNN3glLE3Gdd4RxpRC4v8k9Mg+aIcwnuoIZ2yS8RFfWwPVF rhEz3n9EMwnakLz72G/phcV9a2bKogQqFp0wVbjNI01VJVuduSrjnr+VpH35vNVhl9/rteD306fV uIWwLs2qQ5qTMhhiPhNzIBeq6HQSmAe9fg2jHDQPisSmSOij0Dxqw9M8zZXHWNCHd1OzE8mfRqkp KWXYJWiIpCTbkXfegpReHRdeKfCH8F8sj3/r03kwguVg71K4RxWzQ5nnTyKdDSW2RXOckThLmxXo RWnFOKP+lPMZu6hz3/1NGgvQ2/jCTaR14ohjajj/de3PLgMEBDh8+QwyImvjmmzIg0vuH3wT3/Jq TYJReDWQdI5Ngh7xL36D8SOZ9vaaG+gfXIsBxzUeBrAwzoSruR5GByl6ZSNM3ffr7B2GB6/dQC2B EG+ikCg/Nzdav4Ia2FazTSgdG8WHQnsoz8h/JXvql52YFDOhXn29h2m8pGEx5kHjRPH1c2V89Lfp ZXoKwiGSdTqGIpHJury7yYomuMvMOVOkk4my2f48F2pywRLT06V05qrJeaYkUWXA1bZ7vlfM7ZZw iCPD+wieIPbaVQSnj0lvCCLSzG1+8hvps3gC0KPsWf1SlMkrUnpw0sZUFy3fUjy5+SWe3Ig9Ods0 Zk/qhZ09sqfFvZ3j/4OOSSpQJ98cRoP2hPcrheCu6rjLNrYSusU9eCodOeyqVlPzHAvzjKIM30IX yeBJaTMqNjNzSi3yr7KCaLEhepIepHC27HrwFAoE0VYhB0pY3CcSD7SwlYnXmpEpo+sUWxYb1Nh9 6zxut05FYK3PS4lgnjiSTfBe97w/eKWGimHSo4TnBw+cWtXannjbA3jUhhCeaNwerd9uVPRbhJhd QrEycjSHMxe3Kv3ubD52RBPVcSg1vMjmHu77YtZaFCfGAZI4Mqe/znucLxPSAWR+XTQQ43q6rCcY 0NSLV7PyacdeDndL2uJkgTRxkrWVacuDCtIWIf0ELl5hN/Rxu0O594wex4OgGwQGKPfojGUG7q2E I78fEMLB0In40w56rsvSkdjKGlaCVYU38Gyvi+/FdfoiPoubndMuXra9PnA/iOtIRLfbt0ObwqK0 +tC2QTcShMFSp6iKSXED2KlPZ07P5footYmZntOPnvP6rs/hPWnuO9SAWN+5su9cEZ6hnhqO/Qvu pcjOsSisrXGNkRZtcxnnJF5LdhosQtbW1hJh2UeP6HCDFYj1DcT4Rl9eO6mPL1wPN8kzwvWP5EB0 B9ArYWAQBk4gTxriTCOPYkwStEXDjVTTY1+WauFIx3pAG2qhJiy4tqeUeAVmvTdCofLbbSLeBTce TD0bTEKjXOhxhb7fS4uK3koaXnrUAN9KlCs2QAzHOiU3LBKc1Un3TaR1KcqaftUchCCnYQLJQv3C bTPrdyqpV1YBJ04Yr0MTpHUxd8eh6/HZWM0JHNplrZuA2+Ax98WMTZNtPZpuDXvHhCdsG9jh0HB9 ki1Kf+4gnu3DSRdjXZZcjm0wezZEmBD88efhRmW99qTGQtxAHvgMvv/irVfWX0B/oqsx8NinNyCN FFo7EnPyDfv0Apr8aXjDnt+odR6z6OLipkoUDqCDQjzHjOBaLsw0yr2ivWvuOIgqFM2LsCusSSQI BSoCaAS+x+VpJtqfqtA/TxC0+OqwTT/OYOkSuuE8RGkXj6vV9mQDXA9ZtGEmXAJN4NJZgD3muQVH RIAjybA2XJL80uo1xO326czvz3vQUHBTPPh7pLoBGyXa4ShACnu5iErUCVsqk4cunLF/3Y6a8ULc 4D+GM38+pa/64prB6hoR/dUO3bpA3L9dfwzepg2X8c+v8gb9ZaLfnzPxk7Ed9a2lvm2LbxnAIkUM 6qiy+L2hvjcEtRcOphko3umn3i76520iWEoBYgImMzFO9QWJiPz1Gn32KaxogJ/x8hdPDPvJu/PD NvuHwwIM9oAXT9kjL3qgvZ9E/sENjiLPl0BpwoN4Zs4I90GuKGFL+p322A2XEctKHfBSlsFZfxE4 YBUC/5PJDNysS8WNnXxBtpmi+NzDxXwenDD0rUsDzAqXY6BKqMUugkRRqWNoRr5tL26rs1CF1kcL LYHTMB1JMxsKulZDD77AFVyNSACFWpDIVRsVoAa5fbjx9XODLePrCFUouUDEhs34w8VhkB39Vk/c Ups3sQBMXiAoGa79Wut4c0jWsLIeYghObip1tLImU2wJDLMZZUkH4fywf8osk8RssOZGaurip35A TZPdjnhRFvxqxj69oqmUHY3t4U2iGGHJPUwhjI6mEQlTsDzl+VYIvRCGAGvjs3UcrQ3/icwlEZFM sqJp0NVMA0njOVsxyedQZ1kxijrrAD3TboCQUm00DJEU41GcZbcKXA/owHSL/EZ7AlSCQzLniK4W QWcwzOsZoC5tRnYDLHw/8zFMZEQtMQmf8uc3Udjl1IVZbdp5RA3my91bs7nLCtBEXSusQpBA6ekb aGr05FBpFuJ5GQuR3eO84Vn9/bmPH9Dyv0afK8coP/+r2dra2Unm/283Wg/5X1/i85D/9ZD/9R+S /8V0LwYtXndUVef88d+syZ0dUYpZ3S408FW3S96nuvrh7LD7/nT/1dv97vHJ+eHp0f7BYff47fs3 h28PYXQw4QooT8HnmdgcujCwe45a96aSzsB9HbjDeMaZuDVZ1BAlbLxHU7vxDvyp+Re/me9NA/P1 q+wnwT2cynBKz7ifVY3vMGTdnTk18LiMd2HBOzYzOg5q2D7awssrwLdqKGEvy4HN4Ik7g1lM9xZ6 ATmmQo6ijSO2vliPX37NL79OXH4pivehfIxTBc/ATTLs/wAlFr9gRiEuJRIR5yh7j6fsveO3J/a0 I66k3CB5AzWqe22PL51ZpxIvrPY92tzm9tLwdUzX4PcQAezA2kpLctOWaDadgiVT29Al0peJyh3m h9+tidQ4zWHm0oQHGEsYqH5XZJQZjj9O5V415RHIVfA9TSwUp+4xPUnPsAOXcYgjG0YtG9YJqsH3 COEXhaMFQ/9r7JQ4tzIyRxEusWkm075wkfWpBKqVmQgK8H+KbDFVSZCAY/gaohBGDfE+OdxRCFqV tOik18QlpLoGtWjLUCSolc2EkzsS4uEd42qgVO5Ysiq7t3wnI2VCLNyVdH5WTOq5GswivzOS4N8U pWK8hKlKUdSfvxPmiQ4sogSM6CxbdCGZ9VzlYGOmLotFcsVpsvmty9BeO1oWJoGGVM0Wp7zXVWz/ JoNjhEAluCZOQ1A+9Du/50b4fO+2vJ8XcxuLgtp1ESmvJm5oyakyGOMGMhhPxjCm3HG9jpfkmhwR 0QBOBURiUKg4EQUsKCChARDiBNR8a6xOmzqx8vreYelKcgIsXQGm7aKiOv4TcZ/ZzRewUD6EQqBS eDUE+GSwxzFBfC+To404DZNhQfykgYfCen1Dxb6qGLMUUX5ssgJMCfprbTjCizdbhsuSIAvQR8qh IAiLZptE+RRwIqu8cDxYEgwRK18wG7FcnAOb2hJQAEu2nsMRDzQhRPYxy9kVjpm+GSvRoLQIRFAx 350Vm6HPYTkkXbQ4IIYHcE0Z9mhvEEH5+9yZLWEudjDE1ZcPpe0qWm3KAdQRaoY95A4T+2YGQE4q wBlrGUGzBTNQsiwr5iMsdOAQtn3P4mIlmEt4T9m4qMhBJ5ipBjcRap+LL2GrwEs6kmQu7KITe/CK oIX4IzLj/WyFcL985Vt0Btl05odg4ail85Adi5xbQvGGPmE5GUbJMGxs4xYUzNPUn2U8+sXVhtkS XamDORZX0j/Hi+h+U2Jl6LSlO5pOOCjINyhIN8jNNpD+TzhyKFyxpICjtwzxyPG2RFjjdo54gah4 VQdewvlvQyq3PNJNnMIeQFc64eMAyV3ZQ5Bj+wL3o9ywzd2SmIwvxM553HWUgT2hq6npJb4I6yRi iCbtR2ovDw/e7J+mT9vRb+Wcd1NB+yW07s/ebH34PHwePg+fh8/D5+Hz8Hn4PHwePg+fh8/D50/6 /B8g/D5RALgBAA== --========GMXBoundary327611138145966--