restart;simpson_n := proc(f,a,b,e) local df4, res; df4 := eval(diff(f(x),x$4),x=a); res:=evalf(((b-a)^b*df4/(16*180*e))^(1/4)); return round(res); end proc:simpson := proc(f,a,b,e) local h, res, df4, n, sum1, sum2, j; #n:=simpson_n(f,a,b,e); n:=4; #print(n); h:=(b-a)/(n*2); sum1:=0; for j from 1 by 2 to 2*n-1 do sum1:=sum1+f(a+j*h); end do; sum2:=0; for j from 2 by 2 to 2*n-2 do sum2:=sum2+f(a+j*h); end do; res:=h/3*(f(a)+f(b)+4*sum1+2*sum2); return evalf(res); end proc:f:=(x)->sqrt(1+x^5);

NiM+SSJmRzYiZio2I0kieEdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUklc3FydEc2JEkqcHJvdGVjdGVkR0YvSShfc3lzbGliR0YlNiMsJiIiIkYzKiQ5JCIiJkYzRiVGJUYl

plot(f(x),x=1..Pi);

LSUlUExPVEc2Ji0lJ0NVUlZFU0c2JDdTNyQkIjMzKysrLCsrKzUhIzwkIjM+aTM5a05AOTlGLDckJCIzKD1FVWhpIW9ZNUYsJCIzLTxfaEtwMi06Riw3JCQiM2UoZnklKT0oSCgzIkYsJCIzKFFJK1QxUXRlIkYsNyQkIjMjSFh6Mnd1SDgiRiwkIjNfLWNUO3E7JHAiRiw3JCQiM0I4KnpPNmIqeTZGLCQiMyopbz0iZT9GLyI9Riw3JCQiMyVmRyopKkhwckM3RiwkIjM4VURpWWIoeSQ+Riw3JCQiM3I1JEc7KlI5bjdGLCQiMzAzLGl4bWpsP0YsNyQkIjMhPT8lW2NaMjY4RiwkIjM4LydwMmlvdz8jRiw3JCQiMzR0PE07ITNsTiJGLCQiM0w5OU05bylcTyNGLDckJCIzU3koZTZkJnosOUYsJCIzQy9VV0JVTktERiw3JCQiM11aaCR5UXokWzlGLCQiM08lWyRbcyo0YnIjRiw3JCQiM0R1OyMqUS9UKlsiRiwkIjM3aCYpW14uMycpR0YsNyQkIjN5ZWZFJW8sY2AiRiwkIjNaW2Y3SUFbKTMkRiw3JCQiM3AhZixgZiMpPmUiRiwkIjNWNCgpeSYpUSFHSSRGLDckJCIzUlkjZXM+em1pIkYsJCIzemdJbkcvKik+TkYsNyQkIjMnPSllJVEqekVuO0YsJCIzYzoiKlxXKUhncyRGLDckJCIzTlA1IykpKj5gOjxGLCQiMyV5O1kpcGFOIylSRiw3JCQiM0VcXCQ9JXpUYzxGLCQiMyEpUWByQyhwITRVRiw3JCQiMylIP1BMLHlSIT1GLCQiM0AvW1VyWyFSWyVGLDckJCIzdCJvZ2B0d2clPUYsJCIzZTU4YyJmN3N0JUYsNyQkIjNnVCozZHFsQSo9RiwkIjNYbkcnSHdVZy0mRiw3JCQiM2EqWz8jRyZbaSQ+RiwkIjNfaXIneWppPEomRiw3JCQiMyd5QDUlKikqUkApPkYsJCIzYSspKio0ZXY1aSZGLDckJCIzcUhoYXdFR0M/RiwkIjNiXFQsJj1sXyJmRiw3JCQiMyV6Wyx6RFIocD9GLCQiMyY9KFJzeCJmTkMnRiw3JCQiMyJlRCFmcWUmcDYjRiwkIjNHVSNbSypccidmJ0YsNyQkIjMlKnl2dHEhZSFlQEYsJCIzeDtXJSpcIylIOXBGLDckJCIzXGMjR0hkXEM/I0YsJCIzaCYqR2cmUS4hb3NGLDckJCIzKjMlR1NRLkpbQUYsJCIzS3cvRDRWPVh3Riw3JCQiM3M/PEZdaDwkSCNGLCQiM1czQEd4MyZlLSlGLDckJCIzPUU+QEloZU9CRiwkIjMqKWZpLVEmUl9TKUYsNyQkIjNSYWx0NmJ5JVEjRiwkIjNMJXosaGZeJFIpKUYsNyQkIjNuRzB6KnolNEdDRiwkIjNPIVtYPz1pNUMqRiw3JCQiM24oZXEqM25MdUNGLCQiM0lxI2UhMydRQW8qRiw3JCQiM0wjeidlXihRaV4jRiwkIjNHVzJGJVsuJDQ1ISM7NyQkIjNxPGglPnZbP2MjRiwkIjNiQzNWKEhCYTAiRmN2NyQkIjM2THVBI1JfXmcjRiwkIjNlSy5hSi8pKio0IkZjdjckJCIzVDkzMjsiMy1sI0YsJCIzRUxVJz1AcHg5IkZjdjckJCIzO19VY0AkZVVwI0YsJCIzdSRcKy4xK2Q+IkZjdjckJCIzYzRgTi9bUFNGRiwkIjMqPTtPUVxwckMiRmN2NyQkIjMwIjRkJGUveiV5I0YsJCIzS1BLY196KnpIIkZjdjckJCIzJD4nKm8nej1ASUdGLCQiM1d1REVNSUVeOEZjdjckJCIzNyFmeicpPWRfKEdGLCQiMy4lUVpLaiFRMDlGY3Y3JCQiMyoqUmJWXCJcbSJIRiwkIjNuJXpJMTFgaVgiRmN2NyQkIjMvTEFhMiEqM2tIRiwkIjNZL2hAJSllImZeIkZjdjckJCIzIlFebkB5PmwrJEYsJCIzJmZNQkYpZV5xOkZjdjckJCIzclwhKXpQKWY8MCRGLCQiMyU+JT5BM20tSTtGY3Y3JCQiM0UsdmBWLjEmNCRGLCQiMzptK1sxckMpbyJGY3Y3JCQiMyEpKioqXC9sI2ZUSkYsJCIzLTU3Ljp4Pl88RmN2LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiIUZgW2xGYVtsLSUrQVhFU0xBQkVMU0c2JFEieDYiUSFGZ1tsLSUlRk9OVEc2JCUqSEVMVkVUSUNBR0ZfW2wtJSVWSUVXRzYkO0ZeW2wkIithRWZUSiEiKjskIjJ1VjxDSGU/NCJGY3YkIjFgTj9VS1QleSIhIzk=

a:=1;

NiM+SSJhRzYiIiIi

b:=Pi;

NiM+SSJiRzYiSSNQaUdJKnByb3RlY3RlZEdGJw==

e := 10^(-4);

NiM+SSJlRzYiIyIiIiImKysi

Int(f(x),x=a..b);

NiMtSSRJbnRHNiRJKnByb3RlY3RlZEdGJkkoX3N5c2xpYkc2IjYkKiQsJiIiIkYsKiRJInhHRigiIiZGLCNGLCIiIy9GLjtGLEkjUGlHRiY=

approx_s:=simpson(f,a,b,e);

NiM+SSlhcHByb3hfc0c2IiQiK0tRV246ISIp

genau:=evalf(int(f(x),x=a..b));

NiM+SSZnZW5hdUc2IiQiKz9oV246ISIp

fehler_s:=abs(approx_s-genau);

NiM+SSlmZWhsZXJfc0c2IiQiJSlHIyEiKQ==

restart;simpson_mit_n := proc(f,n,a,b) local h, res, df4, sum1, sum2, j; h:=(b-a)/(n*2); df4 := eval(diff(f(x),x$4),x=a); sum1:=0; for j from 1 by 2 to 2*n-1 do sum1:=sum1+f(a+j*h); end do; sum2:=0; for j from 2 by 2 to 2*n-2 do sum2:=sum2+f(a+j*h); end do; res:=h/3*(f(a)+f(b)+4*sum1+2*sum2); #res:=h/3*(f(a)+f(b)+4*sum(f(a+j*h),j=1..2*n-1)+2*sum(f(a+j*h),j=2..2*n-2))-((b-a)*h^4)/180*df4; return evalf(res); end proc:f:=(x)->exp(2*x)/x^(2/5);

NiM+SSJmRzYiZio2I0kieEdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYtSSRleHBHNiRJKnByb3RlY3RlZEdGMEkoX3N5c2xpYkdGJTYjLCQ5JCIiIyIiIkY0IyEiIyIiJkYlRiVGJQ==

j:=(x)-> piecewise( x=0, 1, x<>0, exp(2*x)-p4(x)/x^(2/5) );

NiM+SSJqRzYiZio2I0kieEdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkqcGllY2V3aXNlR0kqcHJvdGVjdGVkR0YuNiYvOSQiIiEiIiIwRjFGMiwmLUkkZXhwRzYkRi5JKF9zeXNsaWJHRiU2IywkRjEiIiNGMyomLUkjcDRHRiU2I0YxRjNGMSMhIiMiIiYhIiJGJUYlRiU=

plot(f(x),x=0..1);

LSUlUExPVEc2Ji0lJ0NVUlZFU0c2JDdkbzckJCIzYCoqKioqXG41OyJvISNAJCIzIz1eKUg8N15dPSEjOzckJCIzIioqKioqKlw4QUJPIiEjPyQiM0lyKFI3NlBWUyJGLzckJCIzMysrKy1LW1Y/RjMkIjNNNThEMz1yJj4iRi83JCQiMyIpKioqKioqcFVrQ0ZGMyQiM2pFWSpRRiM+bjVGLzckJCIzKSoqKioqKnpMMGVTJEYzJCIzIT0oWydIZGFSeCohIzw3JCQiM3MqKioqKlxTbXAzJUYzJCIzW3Yob0AvNyopNCpGRTckJCIzTCsrK3N1N29aRjMkIjMsa0xsJSopemtjKUZFNyQkIjNpKioqKioqUiYpR1xhRjMkIjNZdGtSSSIpKj44KUZFNyQkIjN5KioqKioqNCdcLzgnRjMkIjM0eDkqPk9NJG94RkU3JCQiM1MrKyt2MWg2b0YzJCIzPzssKj43K3pYKEZFNyQkIjM5KysrUzx4I1woRjMkIjNDJzNubShHcik9KEZFNyQkIjNXKioqKioqNEckUjwpRjMkIjN4Y1MlSHMiR19wRkU3JCQiM00qKioqKlwlXERPJipGMyQiMyMpNHQiKVtcUWFsRkU3JCQiMyQqKioqKip6cWQpKjMiISM+JCIzXnFiXF5qVklpRkU3JCQiMycqKioqKio+IyoqM0U3RltwJCIzMjRFbENXIypmZkZFNyQkIjMyKysrTkBLaThGW3AkIjNxJj5yUjhkJkhkRkU3JCQiMy0rKytbVmIpXCJGW3AkIjNNZChHSywjR0liRkU3JCQiMyopKioqKio+Yyd5TTtGW3AkIjNRMCoqXGVsb2JgRkU3JCQiMygpKioqKioqKSk0RDI+RltwJCIzNFNDMUpIJUgxJkZFNyQkIjMmKSoqKioqZlQ6KHpARltwJCIzR2w6bEUjUmUjW0ZFNyQkIjMpKSoqKioqejRiUWwjRltwJCIzNShILSUqKT4pSF0lRkU3JCQiMyMqKioqKioqeloqejckRltwJCIzeGQ6UGI7Z2NVRkU3JCQiMyYqKioqKio+WU1AZyRGW3AkIjNaXSUpW3N5R2hTRkU3JCQiMzIrKytYVEZ3U0ZbcCQiM3NuXEFMMzctUkZFNyQkIjM9KysrcU1yVV5GW3AkIjNnOmU9IilwUUtPRkU3JCQiMyYqKioqKio0el8iNGlGW3AkIjMsWSZlWWdzN1ckRkU3JCQiM3gqKioqKlw7aEVHKEZbcCQiM0daKlE2QiJvKUgkRkU3JCQiM24qKioqKipSJnBoTilGW3AkIjNjMEolNCllISoqPSRGRTckJCIzKioqKioqKiopPSlIXDUhIz0kIjM0dio9KSo9YSRSSUZFNyQkIjMlKioqKioqei8zdUMiRmd0JCIzPjQzUndmJDMmSEZFNyQkIjM1KysrSiRSRFgiRmd0JCIzN1hobyIpZngjKkdGRTckJCIzNysrKylSJ29rO0ZndCQiM3RGW2YhUSUpeiZHRkU3JCQiMywrKysxSjp3PUZndCQiMylIdUpjNiI+VUdGRTckJCIzMysrKzNFbiQ0I0ZndCQiM1pPKSpHJGVFNiVHRkU3JCQiMy0rKysvUkUmRyNGZ3QkIjN1diNIelQqXF1HRkU3JCQiMyIpKioqKipcS100XSNGZ3QkIjNIYkxpbScqcHFHRkU3JCQiMyQqKioqKipcUEF2ciNGZ3QkIjMhPUlnem9DKSoqR0ZFNyQkIjMoKioqKioqXG5IaSNIRmd0JCIzYFliImU9YV8kSEZFNyQkIjMqKSoqKioqcCpldjpKRmd0JCIzekVpOyxVMnRIRkU3JCQiMyQpKioqKip6ITQ3VExGZ3QkIjMhUl0kKTNVcFctJEZFNyQkIjM/KysrTFkuS05GZ3QkIjMnPkhYRCVRNnRJRkU3JCQiM3QqKioqKjRvN1R2JEZndCQiM207akxlOkZOSkZFNyQkIjMwKysrJFEqb11SRmd0JCIzc0osJ1IpPjQmPiRGRTckJCIzPysrKyI9bGo7JUZndCQiM1RYXTUnNF1kRSRGRTckJCIzKysrK1YmUjxQJUZndCQiM10kcG9zIyk9eUwkRkU3JCQiMzArKytYaC0nZSVGZ3QkIjN1UCgqKSpSXyN6VCRGRTckJCIzIioqKioqKipRIjNHeSVGZ3QkIjNzT0ApZih5KGVcJEZFNyQkIjMhKSoqKioqSDVrXSpcRmd0JCIzSGMiSCk9N24lZSRGRTckJCIzQCsrKyhSUWJAJkZndCQiM0omKSk0OnAwQG8kRkU3JCQiM2ArKys9PlkyYUZndCQiMyIqUlIiXDIhR3JQRkU3JCQiM1wrKyt5WHU5Y0ZndCQiMzkoZiVIYWw/c1FGRTckJCIzZCoqKioqKlt5KSlHZUZndCQiM0QkKSo0QWAqZSIpUkZFNyQkIjMuKysraV9RUWdGZ3QkIjM7by10I2ZIUDQlRkU3JCQiM0ErKysheSUzVGlGZ3QkIjM4JiozSE5sOzJVRkU3JCQiMzUrKytPIVtoWSdGZ3QkIjMnZlhBIyl6WCpRVkZFNyQkIjNJKysrI1F4JG9tRmd0JCIzNXR4cCMpPnNpV0ZFNyQkIjNzKioqKipSUCtWKW9GZ3QkIjM1TilSVCFmaytZRkU3JCQiM1kqKioqKnBwZSp6cUZndCQiMyZbVyoqemc5NHQlRkU3JCQiM3cqKioqKlIjXCdRSChGZ3QkIjNaQTo7OixFeltGRTckJCIzIyoqKioqKkgsTV5cKEZndCQiMyhIaHk4RW5ZLSZGRTckJCIzUysrKzAjPWJxKEZndCQiM3QnPiNmKiozKkc9JkZFNyQkIjNZKioqKipwPzI3InpGZ3QkIjNLQGNUODElUk0mRkU3JCQiM2ErKytJWGFFIilGZ3QkIjNDV2VsPnVdPmJGRTckJCIzNysrK2wqUlJMKUZndCQiMyFbInpBc0tjJnAmRkU3JCQiM2kqKioqKkh2SmdhKUZndCQiM1FaUyp5JD4jSCllRkU3JCQiM3IqKioqKkhKbmp2KUZndCQiM102JDQ/WiVId2dGRTckJCIzbCoqKioqKltRa1wqKUZndCQiM10mKj5wciVlM0UnRkU3JCQiM3cqKioqKipvMDtyIipGZ3QkIjNQJCpIUEQkeTJbJ0ZFNyQkIjNcKioqKipcdyhHcCQqRmd0JCIzI2U0SG0iekwmbydGRTckJCIzOCsrKyFvSzBlKkZndCQiM2ReQ18jKSk9PiJwRkU3JCQiMzMrKys8NXMjeSpGZ3QkIjNvXC8nWyhlRlByRkU3JCQiIiIiIiEkIjNTXTEkKik0YyEqUShGRS0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiRGX2JsRmhibEZpYmwtJStBWEVTTEFCRUxTRzYkUSJ4NiJRIUZeY2wtJSVGT05URzYkJSpIRUxWRVRJQ0FHRmdibC0lJVZJRVdHNiQ7RmlibEZmYmw7JCIyLyRmXHJvJXlfI0YvJCIyQyF6Wik9Uj0pPSEjOg==

taylor(exp(2*x), x=0, 5);

NiMrL0kieEc2IiIiIiIiISIiI0YmRihGKCMiIiUiIiRGKyNGKEYrRiotSSJPR0kqcHJvdGVjdGVkR0YvNiNGJiIiJg==

p4 := x -> 1+2*x+2*x^2+4/3*x^3+2/3*x^4;

NiM+SSNwNEc2ImYqNiNJInhHRiVGJTYkSSlvcGVyYXRvckdGJUkmYXJyb3dHRiVGJSwsIiIiRi05JCIiIyokRi5GL0YvKiRGLiIiJCMiIiVGMiokRi5GNCNGL0YyRiVGJUYl

a:=0;

NiM+SSJhRzYiIiIh

b:=1;

NiM+SSJiRzYiIiIi

n:=6;

NiM+SSJuRzYiIiIn

Int(f(x),x=a..b);

NiMtSSRJbnRHNiRJKnByb3RlY3RlZEdGJkkoX3N5c2xpYkc2IjYkKiYtSSRleHBHRiU2IywkSSJ4R0YoIiIjIiIiRi8jISIjIiImL0YvOyIiIUYx

genau := evalf( int(f(x),x=a..b) );

NiM+SSZnZW5hdUc2IiQiK3k2YW1VISIq

approx := simpson_mit_n(j,6,0,1);

NiM+SSdhcHByb3hHNiIkIStCVyopSHohIzU=

approx := simpson_mit_n(j,6,0,1) + evalf(int(p4j(x)/x^(2/5), x=0..1));

NiM+SSdhcHByb3hHNiIkIisrZj8zTSEiKg==

abs(approx-genau);

NiMkIip5X0xlKSEiKg==

restart;with(plots,listplot):with(plots,display):t:=(i) -> a + i*h;

NiM+SSJ0RzYiZio2I0kiaUdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCZJImFHRiUiIiIqJjkkRi5JImhHRiVGLkYuRiVGJUYl

eul := proc(f, a, b, h, alpha) local i, n, w, res, t, tmp; n := (b-a)/h; w := [alpha]; t:=(i) -> a + i*h; for i from 0 to n-1 do tmp:=evalf(op(-1,w)+h*f(t(i),op(-1,w))); w:=[op(w),tmp]; #print(w); end do; return w; end proc:f:=(t,y) -> 1+y/t;

NiM+SSJmRzYiZio2JEkidEdGJUkieUdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYiIiJGLiomOSVGLjkkISIiRi5GJUYlRiU=

a:=1;

NiM+SSJhRzYiIiIi

b:=2;

NiM+SSJiRzYiIiIj

alpha:=2;

NiM+SSZhbHBoYUc2IiIiIw==

n:=20; h:=(b-a)/n;

NiM+SSJuRzYiIiM/

NiM+SSJoRzYiIyIiIiIjPw==

m:= unapply(rhs(dsolve({diff(y(k),k)=1+y(k)/k,y(1)=2},y(k))),k);

NiM+SSJtRzYiZio2I0kia0dGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYsJi1JI2xuRzYkSSpwcm90ZWN0ZWRHRjFJKF9zeXNsaWJHRiU2IzkkIiIiIiIjRjVGNUY0RjVGJUYlRiU=

m(1.65);

NiMkIitEI3ppNyUhIio=

appr:=eul(f,a,b,h,alpha);

NiM+SSVhcHByRzYiNzciIiMkIisrKytdQCEiKiQiK180US1CRiokIitLWS5kQ0YqJCIrREEnUWgjRiokIiteSnhzRkYqJCIreFNvTEhGKiQiKyEzPWw0JEYqJCIrTlI/aEtGKiQiKy9ibkZNRiokIis2NyhlZiRGKiQiKyZlTGR3JEYqJCIraSkzcyRSRiokIitSbUM1VEYqJCIrbiYqeiVHJUYqJCIrJTNCM1klRiokIitzXkZRWUYqJCIrXGg2PFtGKiQiK3glM3QqXEYqJCIrJmY7KXleRiokIitrbmdoYEYq

exact:=[seq(evalf(m(t(j))),j=0..n)];

NiM+SSZleGFjdEc2Ijc3JCIiIyIiISQiK3MnSDc6IyEiKiQiKyk+VFtJI0YsJCIrTWlzZ0NGLCQiK29leT1FRiwkIitSJUgqeUZGLCQiK1dOMlRIRiwkIisrNzkwSkYsJCIrSjYxcktGLCQiK2Ryd1FNRiwkIitpdz4zT0YsJCIrVl5IelBGLCQiKzJlK19SRiwkIitEI3ppNyVGLCQiK0Yhbz9JJUYsJCIrSHdLeldGLCQiKygqZixlWUYsJCIrS000UVtGLCQiKyVRQSY+XUYsJCIreHNFLV9GLCQiK2hWSCdRJkYs

i:=(1.65-a)/h+1;

NiM+SSJpRzYiJCIlKzkhIiM=

op(14,appr);

NiMkIitSbUM1VCEiKg==

op(14,exact);

NiMkIitEI3ppNyUhIio=

pa := listplot(appr): pe := listplot(exact): display([pa,pe]);

LSUlUExPVEc2JC0lJ0NVUlZFU0c2JDc3NyQkIiIiIiIhJCIiI0YsNyRGLSQiKysrK11AISIqNyQkIiIkRiwkIitfNFEtQkYyNyQkIiIlRiwkIitLWS5kQ0YyNyQkIiImRiwkIitEQSdRaCNGMjckJCIiJ0YsJCIrXkp4c0ZGMjckJCIiKEYsJCIreFNvTEhGMjckJCIiKUYsJCIrITM9bDQkRjI3JCQiIipGLCQiK05SP2hLRjI3JCQiIzVGLCQiKy9ibkZNRjI3JCQiIzZGLCQiKzY3KGVmJEYyNyQkIiM3RiwkIismZUxkdyRGMjckJCIjOEYsJCIraSkzcyRSRjI3JCQiIzlGLCQiK1JtQzVURjI3JCQiIzpGLCQiK24mKnolRyVGMjckJCIjO0YsJCIrJTNCM1klRjI3JCQiIzxGLCQiK3NeRlFZRjI3JCQiIz1GLCQiK1xoNjxbRjI3JCQiIz5GLCQiK3glM3QqXEYyNyQkIiM/RiwkIismZjspeV5GMjckJCIjQEYsJCIra25naGBGMjc3Rik3JEYtJCIrcydINzojRjI3JEY0JCIrKT5UW0kjRjI3JEY5JCIrTWlzZ0NGMjckRj4kIitvZXk9RUYyNyRGQyQiK1IlSCp5RkYyNyRGSCQiK1dOMlRIRjI3JEZNJCIrKzc5MEpGMjckRlIkIitKNjFyS0YyNyRGVyQiK2Ryd1FNRjI3JEZmbiQiK2l3PjNPRjI3JEZbbyQiK1ZeSHpQRjI3JEZgbyQiKzJlK19SRjI3JEZlbyQiK0Qjemk3JUYyNyRGam8kIitGIW8/SSVGMjckRl9wJCIrSHdLeldGMjckRmRwJCIrKCpmLGVZRjI3JEZpcCQiK0tNNFFbRjI3JEZecSQiKyVRQSY+XUYyNyRGY3EkIit4c0UtX0YyNyRGaHEkIitoVkgnUSZGMi0lJUZPTlRHNiQlKkhFTFZFVElDQUdGWA==

fehler := map(abs,appr-exact);

NiM+SSdmZWhsZXJHNiI3NyQiIiFGKCQiKHMnSDchIiokIihZLVkjRiskIigtO3AkRiskIihWTyNcRiskIigpR2NoRiskIihuJSpRKEYrJCIoP0ppKUYrJCIoJz5kKSpGKyQiKWA7NDZGKyQiKV5rSzdGKyQiKWU6YzhGKyQiKVhwejlGKyQiKSdlS2ciRiskIilnJW9zIkYrJCIpWFhdPUYrJCIpRDN1PkYrJCIpJEd4NCNGKyQiKTJSQEFGKyQiKSNvXU0jRiskIikoZihvQ0Yr

listplot(fehler);

LSUlUExPVEc2JC0lJ0NVUlZFU0c2Izc3NyQkIiIiIiIhJEYsRiw3JCQiIiNGLCQiKHMnSDchIio3JCQiIiRGLCQiKFktWSNGMzckJCIiJUYsJCIoLTtwJEYzNyQkIiImRiwkIihWTyNcRjM3JCQiIidGLCQiKClHY2hGMzckJCIiKEYsJCIobiUqUShGMzckJCIiKUYsJCIoP0ppKUYzNyQkIiIqRiwkIignPmQpKkYzNyQkIiM1RiwkIilgOzQ2RjM3JCQiIzZGLCQiKV5rSzdGMzckJCIjN0YsJCIpZTpjOEYzNyQkIiM4RiwkIilYcHo5RjM3JCQiIzlGLCQiKSdlS2ciRjM3JCQiIzpGLCQiKWclb3MiRjM3JCQiIztGLCQiKVhYXT1GMzckJCIjPEYsJCIpRDN1PkYzNyQkIiM9RiwkIikkR3g0I0YzNyQkIiM+RiwkIikyUkBBRjM3JCQiIz9GLCQiKSNvXU0jRjM3JCQiI0BGLCQiKShmKG9DRjMtJSVGT05URzYkJSpIRUxWRVRJQ0FHRlk=

max(op(fehler));

NiMkIikoZihvQyEiKg==

i:=(1.88-a)/h+1;

NiM+SSJpRzYiJCIlZz0hIiM=

op(18,appr);

NiMkIitcaDY8WyEiKg==

op(19,appr);

NiMkIit4JTN0KlwhIio=

evalf((18-1)*h+a);

NiMkIisrKytdPSEiKg==

evalf((19-1)*h+a);

NiMkIisrKysrPiEiKg==

exakt:=m(1.88);

NiM+SSZleGFrdEc2IiQiK1Q0ellcISIq

approx:=op(18,appr)+((op(19,appr)-op(18,appr))*(1.88-1.85))/h;

NiM+SSdhcHByb3hHNiIkIitZOkJEXCEiKg==

approx:=op(18,appr)+(op(19,appr)-op(18,appr))*(3/5);

NiM+SSdhcHByb3hHNiIkIitZOkJEXCEiKg==

abs(approx-exakt);

NiMkIikmUmY6IyEiKg==

restart;with(plots,listplot):with(plots,display):runge4 := proc(f1, f2, a, b, h, alpha1, alpha2) local i, n, au1, au2, k11, k12, k21, k22, k31, k32, k41, k42, tmp1, tmp2; n := (b-a)/h; au1 := [alpha1]; au2 := [alpha2]; for i from 0 to n-1 do k11 := evalf( h*f1(a+i*h, op(-1,au1), op(-1,au2))); k12 := evalf( h*f2(a+i*h, op(-1,au1), op(-1,au2))); k21 := evalf( h*f1( a+(i+1/2)*h,op(-1,au1)+k11/2, op(-1,au2)+k12/2)); k22 := evalf( h*f2( a+(i+1/2)*h,op(-1,au1)+k11/2, op(-1,au2)+k12/2)); k31 := evalf( h*f1( a+(i+1/2)*h,op(-1,au1)+k21/2, op(-1,au2)+k22/2)); k32 := evalf( h*f2( a+(i+1/2)*h,op(-1,au1)+k21/2, op(-1,au2)+k22/2)); k41 := evalf( h*f1( a+(i+1)*h,op(-1,au1)+k31, op(-1,au2)+k32)); k42 := evalf( h*f2( a+(i+1)*h,op(-1,au1)+k31, op(-1,au2)+k32)); tmp1:=op(-1,au1)+(1/6)*(k11+2*k21+2*k31+k41); tmp2:=op(-1,au2)+(1/6)*(k12+2*k22+2*k32+k42); au1:=[op(au1),tmp1]; au2:=[op(au2),tmp2]; end do; return [au1,au2]; end proc:################################################################ # y' = u2; # y = u1; # # Gleichung umformen nach u3 und ersetzten ################################################################f1 := (t,u1,u2) -> u2; f2 := (t,u1,u2) -> -sin(u1)-0.2*u2;

NiM+SSNmMUc2ImYqNiVJInRHRiVJI3UxR0YlSSN1MkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlOSZGJUYlRiU=

NiM+SSNmMkc2ImYqNiVJInRHRiVJI3UxR0YlSSN1MkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYtSSRzaW5HRiU2IzklISIiKiYkIiIjRjMiIiI5JkY3RjNGJUYlRiU=

a:=0;

NiM+SSJhRzYiIiIh

b:=50;

NiM+SSJiRzYiIiNd

h:=(b-a)/500;

NiM+SSJoRzYiIyIiIiIjNQ==

n:=(b-a)/h;

NiM+SSJuRzYiIiQrJg==

au := runge4(f1, f2,a,b,h,2,0);

NiM+SSNhdUc2IjckN2FqbCIiIyQiKyk0I1smKj4hIiokIitoJUg/KT5GKyQiKy1fd2Y+RiskIitzJHkoRz5GKyQiK3EuOCopPUYrJCIrKkdqMyU9RiskIiteLSwleSJGKyQiK0R5Zz08RiskIis8MHJXO0YrJCIrb3NTaTpGKyQiK2knUj1aIkYrJCIrZDhBdDhGKyQiK3d5Jm9FIkYrJCIrbGI7YDZGKyQiKysjKW9LNUYrJCIrXFM0aCEqISM1JCIrRCoqZVV4RkokIiszWzYiUSdGSiQiK0hgdigpXEZKJCIrLClRWmQkRkokIisnZVZeOiNGSiQiKyg0JFFEdSEjNiQhKyJvQVVcJ0ZXJCErSDBUMj9GSiQhK3VHIik9TEZKJCErLiplP2QlRkokISshPnFvdiZGSiQhK19yVGtvRkokIStUUVEoKXlGSiQhK2grJio+KSlGSiQhKyR5MXhsKkZKJCErLEN2UjVGKyQhK2IhUlA1IkYrJCErKUhDdzoiRiskIStaKXo4PyJGKyQhK3EsL043RiskIStrJCpvZTdGKyQhK0AxWHM3RiskISslW3lrRiJGKyQhK2FSJjRGIkYrJCErNUIzYzdGKyQhKypbJDRLN0YrJCErSlpDKj4iRiskISshPkR5OiJGKyQhKyYzaSIzNkYrJCErKip5aV01RiskIStnMVtjKSpGSiQhKz1ZM1AiKkZKJCErbiY+T04pRkokISsjUUtCXihGSiQhKz44QD9tRkokISsmR3hcbyZGSiQhK3YhUV1yJUZKJCErZEpTPlBGSiQhK18qZnZxI0ZKJCErKUc7JCpvIkZKJCErNUc8WW5GVyQiKy1GX21LRlckIitecydbSSJGSiQiKz0wJjNEI0ZKJCIrbCdRZzokRkokIisjWztFLCVGSiQiKzhbaDhbRkokIitCLShIYiZGSiQiKyJIW2JBJ0ZKJCIrR144Rm9GSiQiK3FNU2F0RkokIitpMSdbIXlGSiQiK01leXciKUZKJCIrK0k7cCUpRkokIispeTo7bylGSiQiK1UjW1YiKSlGSiQiK2RbNG8pKUZKJCIrIW8iM1cpKUZKJCIreCIqKlJ1KUZKJCIrOWwpKnAmKUZKJCIra2ZpQyQpRkokIit5YCU0LClGSiQiK0lnVUt3RkokIittRCwkPihGSiQiKy06NyhwJ0ZKJCIrPFhrXGhGSiQiKy5IJWZiJkZKJCIrQipIPSNcRkokIismXFFORCVGSiQiKzVFbmRORkokIitBQDlUR0ZKJCIrQD4zNkBGSiQiK3F5dnU4RkokIitjOGolUidGVyQhKlhgIWYoKUZXJCErZF8iUSp6RlckISsqKT49KlsiRkokISsmel4xOiNGSiQhKyQ0Pnp4I0ZKJCErISlca2xMRkokIStlNTI0UkZKJCEreSpcU1MlRkokISspPWdxJVtGSiQhK3hwP05fRkokIStnQ0BtYkZKJCErOnVSUWVGSiQhK1h5bF1nRkokISsiKilRQz8nRkokISshSDNQSCdGSiQhK2I1JFxLJ0ZKJCErcmsvKEgnRkokIStDelc2aUZKJCErSGInKnBnRkokIStNNSZbKGVGSiQhK2FTd0djRkokISsuJWVaTCZGSiQhK2pxRScqXEZKJCErbFUzPFlGSiQhK1NLTSxVRkokIStYJClbYFBGSiQhK1QxQ3lLRkokISsza2IheSNGSiQhKydmeWJFI0ZKJCErIXoiZVE8RkokIStlPCJcPyJGSiQhK0k9QCpwJ0ZXJCErVVsuKlEiRlckIisqPkooSFFGVyQiK205SDIqKUZXJCIrZEhqejhGSiQiKzFQRFg9RkokIitTb10kRyNGSiQiKzA4cCFwI0ZKJCIrcyI9TjEkRkokIitIMzgqUiRGSiQiK0l0NiZwJEZKJCIrZmReXFJGSiQiKzNQIjM7JUZKJCIrMUklekslRkokIitQM0ZdV0ZKJCIreiMpZUZYRkokIismNCg0Z1hGSiQiKzhhUltYRkokIisib2lNXCVGSiQiK1ZTaydSJUZKJCIrITNNJ2ZVRkokIitYJ2ZXMyVGSiQiK0E5WXRRRkokIitHUEZIT0ZKJCIrdT8hW04kRkokIitGISk+YElGSiQiK1I1JHlzI0ZKJCIrTHFEI1EjRkokIithUD0/P0ZKJCIrJzNMYWsiRkokIisneSs+RSJGSiQiK0lEOU4oKUZXJCIrJWUiKj0lW0ZXJCIraW5leSgqISM3JCErTk5MPkdGVyQhK0ZOWDhsRlckISs1VysyNUZKJCEra3FuWDhGSiQhK01HUGs7RkokISsoUnEuJz5GSiQhK0QhUjdCI0ZKJCErWV8mW1ojRkokISsqKltUKm8jRkokISt2M1d0R0ZKJCErSXB5RElGSiQhK0whUWM5JEZKJCErJ3AyREIkRkokIStIRkInRyRGSiQhKyp5bHBJJEZKJCErVWc7JkgkRkokIStQIillXktGSiQhK14ncHM8JEZKJCErPXFedElGSiQhK3IqKik9JUhGSiQhKy1XPSV5I0ZKJCErc05ULUVGSiQhKypmKHkpUiNGSiQhK242cHZARkokIStiJWZjJD5GSiQhKyNRXzhvIkZKJCErJUhGYlQiRkokIStvLixUNkZKJCErLi5uMScpRlckIStzIz1QeCZGVyQhKzQ9eVJIRlckIStfWiYpRzhGXF9sJCIrWF91PkVGVyQiKyhRZj5IJkZXJCIrNHYkKmV5RlckIitNRXVINUZKJCIrKHolZmU3RkokIitVIilccTlGSiQiKzc6cmo7RkokIisvcXJPPUZKJCIrZyRIIykpPkZKJCIrWkQ/PEBGSiQiK2dPJEdBI0ZKJCIrSFZjL0JGSiQiKzYpej9PI0ZKJCIrZ2ZJJlIjRkokIiteVVMvQ0ZKJCIrMF13KlEjRkokIitlKikqPk4jRkokIitfcyM+SCNGSiQiK2UqcDBAI0ZKJCIrd0o4NEBGSiQiK2VaKiopKT5GSiQiK0UoKW9ePUZKJCIrUiUpKSkpcCJGSiQiKyllKVFLOkZKJCIrSGgzYThGSiQiKyJIZ2Y7IkZKJCIrNyU9MHEqRlckIitrIT1XbyhGVyQiK1AqKUhLY0ZXJCIrYShHXmMkRlckIit1QGouOkZXJCEreGk8PmBGXF9sJCErUSMqKj1fI0ZXJCErMERgWldGVyQhKzFbMCJIJ0ZXJCErQT8pZS4pRlckISstRHZtJypGVyQhK3dfKXA2IkZKJCErLUNIYDdGSiQhK0ZcYHU4RkokIStpPyMpejlGSiQhK1MoUSVvOkZKJCErTCZbKVI7RkokIStpWnAkcCJGSiQhK3opKnpIPEZKJCErMks7WzxGSiQhKz1yJipbPEZKJCEreDxfSzxGSiQhKyZbZSQqcCJGSiQhK108N107RkokIStSLGgmZSJGSiQhKycqZnYxOkZKJCErPlRoOTlGSiQhK2clXC5KIkZKJCErc1JBJj4iRkokISsjbyNlcTVGSiQhKyVIIlJ5JCpGVyQhKytIZyUpekZXJCErMTleUmxGVyQhK0EtUmVdRlckISsocHpsYiRGVyQhK289TFw/RlckIStMVmU7YkZcX2wkIispeWkiPSMqRlxfbCQiK11jKnBOI0ZXJCIrZUdXU1BGVyQiKyFHLiZmXUZXJCIrJipmVC1qRlckIitpTFRldUZXJCIrJyl5IXleKUZXJCIrNWUxcyUqRlckIis0blFKNUZKJCIrZ15yLjZGSiQiK2tecmo2RkokIislM1Q1QCJGSiQiKyYqKnlhQyJGSiQiK3NmJXBFIkZKJCIraSIpW3Y3RkokIitrd0ZyN0ZKJCIrOyQzWUQiRkokIis4LyplQSJGSiQiK1xXayY9IkZKJCIrXVJcTTZGSiQiK2J1OnQ1RkokIitdKFJDKyJGSiQiK1ZUQUsjKkZXJCIrWXlgayQpRlckIitdVVFKdUZXJCIrKDNpS1cnRlckIispPTY1VCZGVyQiKzgpKXBYVkZXJCIrJlswJmVLRlckIitFOmhnQEZXJCIrIVwmM2o1RlckISpHSkNLI0ZcX2wkIStyU3ooMyJGVyQhK0FNXz9ARlckISs4Ikc9NiRGVyQhK29hc19TRlckISs3PiJcJFxGVyQhKycpPSQzdiZGVyQhK1x3dCRcJ0ZXJCEreiZSeDooRlckIStObiV5dChGVyQhKzx2KipIIylGVyQhK1kvM0onKUZXJCErP14lKlEqKUZXJCErT05TXyIqRlckIStBP0FyIypGVyQhK0NUNSdIKkZXJCErNVttRyMqRlckIStfa1JyISpGVyQhK3JxaUYpKUZXJCEraDhaLCYpRlckISsnNHZ4NClGVyQhK0ZOMEF3RlckISsxV1UhMyhGVyQhKyhRTyZ6a0ZXJCErIlskXEVlRlckISsnUXcoR15GVyQhK2k5OyVSJUZXJCErJD5RMWokRlckIStxZktZR0ZXJCEraDVSXD9GVyQhKyVma3pDIkZXJCErPVNpK1hGXF9sJCIrL1sjXE8kRlxfbCQiKyEqUTkvNkZXJCIrdUNuWD1GVyQiK1UzRWFERlckIitSTGJCS0ZXJCIrUio0eCVRRlckIitcUVhAV0ZXJCIrU0M2U1xGVyQiKyFIXicqUiZGVyQiK142cSd6JkZXJCIrSHhkR2hGVyQiKyNHJUgkUidGVyQiKyUzbiYqZSdGVyQiK2xRIm9yJ0ZXJCIrP2E5dm5GVyQiK0IxTmxuRlckIitMXigpKW8nRlckIitSVHpabEZXJCIrLSZ6Wk0nRlckIitRPDEkMydGVyQiKyFIKFFtZEZXJCIrcDooKilSJkZXJCIrLCRbYSlcRlckIis+ZCIzYCVGVyQiK3MsUVNTRlckIitUeHA+TkZXJCIrc1dedUhGVyQiK05lcTVDRlckIitAZ0BNPUZXJCIrKnAoKjREIkZXJCIrOyNHJnBtRlxfbCQiKztzXSh5KSEjOCQhKzZdMDFbRlxfbCQhK0teMkw1RlckISsjeSZRazpGVyQhK1k0cXA/RlckIStAdWBXREZXJCErMkYhWylIRlckIStSeCNvUSRGVyQhK1daUlpQRlckIStHK3dqU0ZXJCErMTtuTFZGVyQhKyhlIlFiWEZXJCErTExsRlpGVyQhK0ZNd1xbRlckIStcJilcQFxGVyQhKyVSWkolXEZXJCErTnlbOlxGVyQhK0cmcChSW0ZXJCErW0FwPFpGVyQhK0YtUV5YRlckIStQR05WVkZXJCErZz5cJzQlRlckIStJays5UUZXJCErVlJSKlwkRlckIStfM1NjSkZXJCErI0chKSopeSNGVyQhKyUqKls3UyNGVyQhK21NVygqPkZXJCErRGMoPWUiRlckISt4JykpKWU2RlckIStvWzhHdEZcX2wkIStLImYjeklGXF9sJCIrJG9XZjYiRlxfbCQiK0MpUXFAJkZcX2wkIistdUomPSpGXF9sJCIramJUKUgiRlckIitOSiR6bCJGVyQiK0w7JFIqPkZXJCIrIVxgTkkjRlckIispemVVZSNGVyQiKz1UJVEkR0ZXJCIrcXZYXUlGVyQiK0wsaEtLRlckIitOSj16TEZXJCIrRT5WKlskRlckIisleSYpSGMkRlckIitVVCUpKmYkRlckIitJI3AuZyRGVyQiKy5gRmxORlckIisuWWgmXCRGVyQiKyJHZ0ZSJEZXJCIrKVsnUWVLRlckIispKmVXJTQkRlckIitBXDkuSEZXJCIrNXAicG8jRlckIitHT1JbQ0ZXJCIrWGBQIT4jRlckIitAKSp6Oj5GVyQiK1swckY7RlckIitpWEFIOEZXJCIrXSxdQjVGVyQiKyF5WnE4KEZcX2wkIitxYiQpSFNGXF9sJCIrNVomR1cqRmlmbSQhKyE0MyYqMyNGXF9sJCErKFxlRC8mRlxfbCQhK0U/RCgpeUZcX2wkISsvJWYoZjVGVyQhK0RNJVxKIkZXJCErOXAxXzpGVyQhKzk/OXA8RlckISs1Jz1XJz5GVyQhKz9wUk9ARlckIStmdCRRRyNGVyQhKz53dzBDRlckISslKXBbLERGVyQhK1x5Y3FERlckIStPWCZHaCNGVyQhK0MlZiVHRUZXJCErR2x2PEVGVyQhKz5FUCJlI0ZXJCErO2c8P0RGVyQhK0ZLRU5DRlckISsnZlZ6SyNGVyQhK2RAcyo+I0ZXJCErJykzR18/RlckISshZWV1KT1GVyQhK3gmSHNxIkZXJCEraDpvODpGVyQhK3JFKiozOEZXJCErdSIzYTQiRlckIStOI3BAdilGXF9sJCErWyZ5c10nRlxfbCQhKy11WVVVRlxfbCQhK0dWYCEpPkZcX2wkIislUT0/YyNGaWZtJCIraEw2WUNGXF9sJCIreVJZb1hGXF9sJCIreWtmLm1GXF9sJCIrbWU1TCYpRlxfbCQiK28nKipSLiJGVyQiK0RCKTM/IkZXJCIrdDZlXzhGVyQiK0dqKnlbIkZXJCIrXlAiZWciRlckIishZjZicSJGVyQiKyxYTyd5IkZXJCIrR2YleiU9RlckIisjUUcrKj1GVyQiK1c8ZTc+RlckIitbJ3BkIj5GVyQiK0NTJSoqKj1GVyQiKykpeWpsPUZXJCIrKltjTiI9RlckIitBcGNXPEZXJCIrY2pvZjtGVyQiKy0hcCtjIkZXJCIrKSo+KnBXIkZXJCIrUC4lPUsiRlckIitYNTQnPSJGVyQiK1tuSFQ1RlckIitUIylvISopKUZcX2wkIishKUhnNXRGXF9sJCIrWG1bKm8mRlxfbCQiKyQ0KT1XU0ZcX2wkIis7PlkiUiNGXF9sJCIrJEdYIXl1RmlmbSQhKzsmUnFxKUZpZm0kISt3MGJbQ0ZcX2wkIStPPSQ0KFJGXF9sJCErJTQnKVFVJkZcX2wkIStZXVQleidGXF9sJCErYEtmcSEpRlxfbCQhK2dbbVQjKkZcX2wkISslbzcpSDVGVyQhKzopeko3IkZXJCErXHVlLjdGVyQhK0hvXHE3RlckISskUjhOSyJGVyQhKz95UWk4RlckISs6YywoUSJGVyQhK3VfVihSIkZXJCErcl0jUVIiRlckIStEJilcdzhGVyQhKz8nKSplTSJGVyQhK2wzZi04RlckIStzX0RaN0ZXJCErbHRuIT0iRlckIStCJFFQNSJGVyQhK3JWUzw1RlckIStGYTpGIypGXF9sJCErNyRSeD8pRlxfbCQhKyVmOnQ3KEZcX2wkIStyKnp4KmZGXF9sJCErKVx5NyRbRlxfbCQhK3ZDM1NPRlxfbDdhamwiIiEkISsnPm8pMyEqRlckIStvSGYoeSJGSiQhK149eWpFRkokISs6NGRLTkZKJCErZzNNJ1IlRkokIStCJVJsRCZGSiQhK113WjhoRkokIStBNDdtcEZKJCErKFwjKT0ieUZKJCErYVFWWScpRkokIStWTGJqJSpGSiQhK3lUXUQ1RiskIStCUjIsNkYrJCErQHonPTwiRiskIStVN2FPN0YrJCErKCo+bCRIIkYrJCErJFJKPE0iRiskIStPJHAkejhGKyQhKzghNGBTIkYrJCErOjNhPTlGKyQhK0dVUT05RiskIStJbGEvOUYrJCErOyhcclAiRiskISs4O3NPOEYrJCErN0Y6JUciRiskISt6YGk/N0YrJCErW2RfWjZGKyQhKzoqXGoxIkYrJCErKio+Oid5KkZKJCErL2sleSYpKUZKJCErOmg0Iyp5RkokIStjOiM0IXBGSiQhK08ubSUqZUZKJCErX3okPilbRkokIStZUHdwUUZKJCErTGlyakdGSiQhKyJlJj1vPUZKJCErRCdvaycpKUZXJCIqajJjIXlGVyQiK3gyVkI1RkokIisjKTQyWj5GSiQiK2VUVFlHRkokIisyKEcmPVBGSiQiKzN6JCpmWEZKJCIrMCUpW21gRkokIit0PkZMaEZKJCIrNyk0WSZvRkokIitAITNUXyhGSiQiKyFlIXpNIilGSiQiK1c8SnonKUZKJCIrNkpDXSIqRkokIitKQFVTJipGSiQiKyNwUkwlKSpGSiQiKy9PTjA1RiskIittb3A7NUYrJCIrKVtKIj01RiskIislZjwnNDVGKyQiK18pZkoiKipGSiQiK1VPJGVqKkZKJCIrZT1icCMqRkokIismNCM9QCkpRkokIitzI2YoKUgpRkokIitAYjU2eEZKJCIraEBYbnFGSiQiK3ddNHhqRkokIislWzEiXGNGSiQiK3lsNiMqW0ZKJCIrOyVwVDYlRkokIitmKlxFSyRGSiQiKyp6dFVfI0ZKJCIrISl5N0Q8RkokIishKT1cMiQqRlckIitUR0dpOUZXJCErN2l1T2lGVyQhKyQ9N1dQIkZKJCErPjNhLEBGSiQhKyJcRzEhR0ZKJCEreWw9bk1GSiQhK2U8aic0JUZKJCErdUFHJW8lRkokISslPiJSRF9GSiQhKyopcEA6ZEZKJCErOiJIIlxoRkokIStuS3VBbEZKJCErRnoyS29GSiQhK0I1c3RxRkokIStHTipcQyhGSiQhK2hmNFd0RkokIStYV0FxdEZKJCErYl9rQnRGSiQhKyNbP2Q/KEZKJCErIj0jKSk9cUZKJCErdHRibW5GSiQhK0U0MGBrRkokISt1cFIkMydGSiQhK2A7PmpjRkokIStYL1UpPiZGSiQhK09DSCZwJUZKJCErMDs0Z1RGSiQhK0ZCMCpmJEZKJCErJD1qIz1JRkokISspNCtPVSNGSiQhKyx0bz89RkokIStARSlbQCJGSiQhKzVFI0c2J0ZXJCErRW0lUloiRlxfbCQiK114MyxkRlckIisqcDcpUTZGSiQiK1tPMyhvIkZKJCIrLGF4NUFGSiQiKzhtKGVxI0ZKJCIrSXBeb0pGSiQiKyY9KSlcZiRGSiQiK08yeSIpUkZKJCIrKG9MY0slRkokIitnamZCWUZKJCIrJz0oNHRbRkokIis3XyxzXUZKJCIrNSlbKD1fRkokIitZZ0Y3YEZKJCIrVzo/X2BGSiQiK0JjeVFgRkokIitARyZIRiZGSiQiKydmeWk6JkZKJCIrKHphNCpcRkokIitKaHR6WkZKJCIrRzgoZV8lRkokIisjZUJJQiVGSiQiK2lgPTBSRkokIit3RGZZTkZKJCIrIWZVOzskRkokIitnJz5bdiNGSiQiK2wuaklCRkokIitfLmIkKj1GSiQiKy9GKXpXIkZKJCIrTT9GIykqKkZXJCIrYy5mJVsmRlckIis3IVxyLSJGVyQhK2R4Ol5MRlckIStfJ3BIaChGVyQhKyV5XkE8IkZKJCErOypwWGMiRkokIStAJD1dJD5GSiQhK1RqYiFHI0ZKJCErMGROKWYjRkokIStNYCNlKUdGSiQhK0sta1NKRkokIStfSnhnTEZKJCErR3RfV05GSiQhKyM9bzBwJEZKJCErIXpeenokRkokIStdejltUUZKJCErRGYwJipRRkokIStkOCwmKVFGSiQhK0RNeU9RRkokIStmPmNeUEZKJCErKipbJDRqJEZKJCErNXImb1okRkokIStoP2giSCRGSiQhK2Mjb3gyJEZKJCEra0M4UUdGSiQhKzk9cXZERkokIStObWgkSCNGSiQhK0NjNiYqPkZKJCErU1hcJG8iRkokISsmPXA/TyJGSiQhKyUpSDlNNUZKJCErdTR6SHFGVyQhK2NqdzxQRlckIStLQS5sVkZcX2wkIit4QChReSNGVyQiKyZHOVciZkZXJCIrJSkqM3YjKilGVyQiKyJHLSh6NkZKJCIrSWMlKVw5RkokIit2LyI0cSJGSiQiK1dGIzMkPkZKJCIrZUJzUEBGSiQiK3Y1KSo+QkZKJCIrJSkqPmlaI0ZKJCIrYWVLMEVGSiQiKyNbbGtxI0ZKJCIra3I0ekZGSiQiK3UleUgjR0ZKJCIreCdwIlFHRkokIiswQi5ER0ZKJCIrX0VBJXkjRkokIitPLW87RkZKJCIrMzpoQkVGSiQiK0UkcGtdI0ZKJCIrQipHcE8jRkokIislPmhvPyNGSiQiKzdVSUc/RkokIitoUVZMPUZKJCIrSVdgQztGSiQiK24pcFJTIkZKJCIrZ2k6dTZGSiQiKz4kZmBQKkZXJCIrOFFgbHBGVyQiK0U3TE9YRlckIitvTmY2QEZXJCErOixGYUdGXF9sJCErazJDS0VGVyQhKzg4RzJcRlckIStxdjAhNChGVyQhK2UnKUhoIipGVyQhKz8kNC42IkZKJCErWEIhKipHIkZKJCEreCRITVgiRkokISs1X2UqZiJGSiQhKyo+XXNzIkZKJCErVy1dTj1GSiQhKyk+Ok8jPkZKJCErbmAzIio+RkokIStpZGhQP0ZKJCErQXk3aj9GSiQhKypSZXgxI0ZKJCErISplJj0wI0ZKJCErLU4oZiwjRkokIStnKGYzJz5GSiQhK1xrVygpPUZKJCErJUhNb3oiRkokIStQbkYhcCJGSiQhK0VBO3A6RkokIStXYSpcViJGSiQhKyt6UCpHIkZKJCErJUgpKVI4IkZKJCErdylIY3EqRlckISt5J28oMyEpRlckISshW1tzRSdGVyQhK2s4LypcJUZXJCEreFAoPnMjRlckISs8LW5OJipGXF9sJCIrIXlKNip5RlxfbCQiKyRIKWYqWyNGVyQiK3csOktURlckIitTaCk9cSZGVyQiKzFoKls9KEZXJCIrbDxObyYpRlckIit0dGdTKSpGVyQiK1siSCIqNCJGSiQiK24pUTY/IkZKJCIrJyo0TCpHIkZKJCIrbF00ajhGSiQiKzJNKD5VIkZKJCIrKCpSbWw5RkokIitTOy0lXCJGSiQiKzV1MDI6RkokIitBaiRcXSJGSiQiKzBLKHpbIkZKJCIrXW9pYzlGSiQiKzdEXDY5RkokIis5SEhgOEZKJCIrO3onR0ciRkokIitySTssN0ZKJCIrbHNANDZGSiQiKyFwXCIzNUZKJCIrR09ZIioqKUZXJCIrIylSWU15RlckIishUSNHQm1GVyQiKyMpKmY0UCZGVyQiK2AybSE0JUZXJCIraTdhJnojRlckIis7XmkpXCJGVyQiK0xuI283I0ZcX2wkIStTTSkpXDVGVyQhKz5tOXhBRlckIStjIWV4WCRGVyQhK3A7LyJlJUZXJCErQD0zUGNGVyQhK1dNIW9oJ0ZXJCErVFMvN3ZGVyQhK0lxaDokKUZXJCErQWVPQCEqRlckIStrej9DJypGVyQhK3ZvLDc1RkokISt1LmtdNUZKJCEraDQ3eTVGSiQhK0NCUyU0IkZKJCErazNhKjQiRkokISspPi5QNCJGSiQhKyVcaXIyIkZKJCErKHAkSF01RkokISs5eGM4NUZKJCErSWxXdicqRlckISsmUmknRyIqRlckIStFXlotJilGVyQhK2MlelkheUZXJCErdjZpVnFGVyQhK3cxNUdpRlckISsnZnVzTyZGVyQhK0tuYnFXRlckISs4JD52YSRGVyQhKyZmJnoyRUZXJCErSFkpNG0iRlckIStZJiplbHJGXF9sJCIrWjdDaUBGXF9sJCIrSXJTRzZGVyQiK1JCVDY/RlckIiskZjtyJkdGVyQiK2EwI3psJEZXJCIrISk+Im9TJUZXJCIrNUVVKDQmRlckIitqI3lTcyZGVyQiK2M8JT1HJ0ZXJCIrcSNbbHcnRlckIitcRyRbPChGVyQiK1AoXFRdKEZXJCIrX0x5X3hGVyQiK3QyJik+ekZXJCIrS2RIMCEpRlckIisiPXUpNCEpRlckIitzOzhOekZXJCIrc0dQJHkoRlckIitpT2lkdkZXJCIrKlInZWhzRlckIit2KillKipvRlckIitEJUdsWidGVyQiK2wpNHkqZkZXJCIrI2YiR3BhRlckIiskPW5yKltGVyQiK1FhKnpHJUZXJCIrISoqRyZbT0ZXJCIrZUJwJilIRlckIislPitsSSNGVyQiKzdrKXpoIkZXJCIrYSRbOEYqRlxfbCQiKzcnPiIzQ0ZcX2wkISssaihITSVGXF9sJCErRkZ2IjQiRlckIStWTFREPEZXJCErd1hcSEJGVyQhK25naCkqR0ZXJCErSW4leVUkRlckISs8V3U3UkZXJCErTC1SXFZGVyQhK1xzVE1aRlckISsqUkxdMSZGVyQhKyNSUSFSYEZXJCErOVgkW2ImRlckISskSEo5ciZGVyQhKyYqUlczZUZXJCEraGczWWVGVyQhK0BjOkRlRlckISsuaixaZEZXJCErT0VkOGNGVyQhKyNSU3NVJkZXJCErKUg3ND4mRlckIStqJj56IVxGVyQhK2MoKik+ZSVGVyQhK1Y6PzxVRlckIStwbyR6IlFGVyQhK1Y2JCkpUSRGVyQhKyQqPnNNSEZXJCErcXJmZ0NGVyQhK0dAYXI+RlckIStEeW9zOUZXJCErK0JmInAqRlxfbCQhKzkuRGdZRlxfbCQiKzVMNXZKRmlmbSQiKzx6YyQ+JkZcX2wkIitxI2U8IyoqRlxfbCQiK3cpR2VXIkZXJCIrTEBAdz1GVyQiK1BAYHpBRlckIis6cUxfRUZXJCIrO1FhIipIRlckIit3PFklSCRGVyQiK0E4IillTkZXJCIrOChRRnkkRlckIis/aSNbJ1JGVyQiK0B5NC9URlckIis3Ly0rVUZXJCIrcDBdX1VGVyQiKyYqcCg9RSVGVyQiK2MoMypHVUZXJCIrTSZmWjolRlckIisqNHk0LyVGVyQiKytgWiopUUZXJCIrZXlcLVBGVyQiK0AkKmYjWyRGVyQiK0EmM0VCJEZXJCIrUGZmYkhGVyQiK1QkUVtsI0ZXJCIrR0B5TEJGVyQiK1lkKycqPkZXJCIrNTs9WDtGVyQiKyQqeS4mRyJGVyQiKzNlPyQ+KkZcX2wkIisiXGh2XiZGXF9sJCIrdj86Zz1GXF9sJCErV2JBVjxGXF9sJCErcyI9IWVfRlxfbCQhK0U+Ql4nKUZcX2wkIStKRTwqPSJGVyQhK1clW11cIkZXJCErRGQzIXkiRlckISt3eig9LyNGVyQhKzVlSHlBRlckIStpRF0oWyNGVyQhK05uJ3ptI0ZXJCErNTdaPUdGVyQhK3MiPiJRSEZXJCErY25MRUlGVyQhK1VDKEczJEZXJCErZkt6MkpGVyQhKyMpeVosSkZXJCErcm5na0lGVyQhK3clXCIpKkhGVyQhKzYkXEwhSEZXJCErJ2YwPHkjRlckISssUCZcaiNGVyQhK0pKL2xDRlckISs8UjZ1QUZXJCErKWZyVzEjRlckIStrNWNRPUZXJCErYCZSKilmIkZXJCErITRcI1s4RlckIStgJyk9KjMiRlckISsiXGpbQylGXF9sJCErPyg9KG9iRlxfbCQhKygzJlwhKkdGXF9sJCErUFU+bkJGaWZtJCIrZSFcb08jRlxfbCQiK2h4WCYqW0ZcX2wkIitVXWVEdEZcX2wkIis0IVteaipGXF9sJCIrNSNwLj0iRlckIisjNGM3USJGVyQiKzZYXWs6RlckIis/NGxHPEZXJCIrYSZcQyg9RlckIitqLylbKj5GVyQiKyc9Yl40I0ZXJCIrMCtzc0BGVyQiKz5vREZBRlckIitwO29lQUZXJCIrbDQ5bkFGVyQiKy1iK2BBRlckIitnQidvQCNGVyQiKz1bXWZARlckIis4LyM+MyNGVyQiKzJ0RiYpPkZXJCIrSiI0NCg9RlckIisiWylIUzxGVyQiK3YjZl1mIkZXJCIrZCI9cFYiRlckIitmXnBuN0ZXJCIrQFFHKjMiRlckIitOK0pPISpGXF9sJCIrJ2VscjcoRlxfbCQiKzY1TCY9JkZcX2wkIitXMG1JS0ZcX2wkIistPHUjRyJGXF9sJCErInpISFInRmlmbSQhK2tmJXBeI0ZcX2wkISsnKVtfS1ZGXF9sJCErMDlJcGdGXF9sJCErQXVsNnhGXF9sJCErb3VAWCMqRlxfbCQhK0JtcGw1RlckIStSXWAkPiJGVyQhKylST3FJIkZXJCErJ29oYFMiRlckISslZVF5WyJGVyQhK1BbJ1JiIkZXJCErdSk0TWciRlckISshcDlnaiJGVyQhK1Urel47RlckISsiKkgiNGwiRlckIStsPXNMO0ZXJCErQig0MmciRlckISs2a15fOkZXJCErQiQ+KipbIkZXJCErdkgjUVQiRlckIStOdkNEOEZXJCErQmtKRDdGVyQhK1pKQzo2RlckIStzXDxqKipGXF9sJCErVUYiKilwKUZcX2wkIStwbmp0dEZcX2wkISs+X20sZ0ZcX2wkIStDVlwoZiVGXF9sJCErKVteYzwkRlxfbCQhKzk1Yl08RlxfbCQhKyRIUU5PJEZpZm0kIiskXG5KMCJGXF9sJCIrJj5kWlMjRlxfbCQiK1JyejBQRlxfbCQiK0YqPVclXEZcX2wkIisqUjEnNGhGXF9sJCIrJ3AoRyI+KEZcX2wkIisiPTsvPSlGXF9sJCIrTCdSIXAhKkZcX2wkIit6c05dKSpGXF9sJCIrSmsoPTAiRlckIityJCkpcDUiRlckIitLJmUrOiJGVyQiKzojKiozPSJGVyQiK0J1Vyo+IkZX

i:=(33.7-a)/h+1;

NiM+SSJpRzYiJCIlIVEkISIi

op(338,au[1]);

NiMkIis+ZCIzYCUhIzY=

op(338,au[2]);

NiMkIStcc1RNWiEjNg==

plist := []: for i from 1 to n do tmp := [op(i,au[1]), op(i,au[2])]; plist := [op(plist), tmp]; end do:listplot(plist);

LSUlUExPVEc2JC0lJ0NVUlZFU0c2Izdgamw3JCQiIiMiIiEkRixGLDckJCIrKTQjWyYqPiEiKiQhKyc+bykzISohIzY3JCQiK2glSD8pPkYxJCErb0hmKHkiISM1NyQkIistX3dmPkYxJCErXj15akVGOjckJCIrcyR5KEc+RjEkISs6NGRLTkY6NyQkIitxLjgqKT1GMSQhK2czTSdSJUY6NyQkIisqR2ozJT1GMSQhK0IlUmxEJkY6NyQkIiteLSwleSJGMSQhK113WjhoRjo3JCQiK0R5Zz08RjEkIStBNDdtcEY6NyQkIis8MHJXO0YxJCErKFwjKT0ieUY6NyQkIitvc1NpOkYxJCErYVFWWScpRjo3JCQiK2knUj1aIkYxJCErVkxiaiUqRjo3JCQiK2Q4QXQ4RjEkISt5VF1ENUYxNyQkIit3eSZvRSJGMSQhK0JSMiw2RjE3JCQiK2xiO2A2RjEkIStAeic9PCJGMTckJCIrKyMpb0s1RjEkIStVN2FPN0YxNyQkIitcUzRoISpGOiQhKygqPmwkSCJGMTckJCIrRCoqZVV4RjokISskUko8TSJGMTckJCIrM1s2IlEnRjokIStPJHAkejhGMTckJCIrSGB2KClcRjokISs4ITRgUyJGMTckJCIrLClRWmQkRjokISs6M2E9OUYxNyQkIisnZVZeOiNGOiQhK0dVUT05RjE3JCQiKyg0JFFEdUY0JCErSWxhLzlGMTckJCErIm9BVVwnRjQkISs7KFxyUCJGMTckJCErSDBUMj9GOiQhKzg7c084RjE3JCQhK3VHIik9TEY6JCErN0Y6JUciRjE3JCQhKy4qZT9kJUY6JCEremBpPzdGMTckJCErIT5xb3YmRjokIStbZF9aNkYxNyQkIStfclRrb0Y6JCErOipcajEiRjE3JCQhK1RRUSgpeUY6JCErKio+Oid5KkY6NyQkIStoKyYqPikpRjokISsvayV5JikpRjo3JCQhKyR5MXhsKkY6JCErOmg0Iyp5Rjo3JCQhKyxDdlI1RjEkIStjOiM0IXBGOjckJCErYiFSUDUiRjEkIStPLm0lKmVGOjckJCErKUhDdzoiRjEkIStfeiQ+KVtGOjckJCErWil6OD8iRjEkIStZUHdwUUY6NyQkIStxLC9ON0YxJCErTGlyakdGOjckJCErayQqb2U3RjEkISsiZSY9bz1GOjckJCErQDFYczdGMSQhK0Qnb2snKSlGNDckJCErJVt5a0YiRjEkIipqMmMheUY0NyQkISthUiY0RiJGMSQiK3gyVkI1Rjo3JCQhKzVCM2M3RjEkIisjKTQyWj5GOjckJCErKlskNEs3RjEkIitlVFRZR0Y6NyQkIStKWkMqPiJGMSQiKzIoRyY9UEY6NyQkISshPkR5OiJGMSQiKzN6JCpmWEY6NyQkISsmM2kiMzZGMSQiKzAlKVttYEY6NyQkISsqKnlpXTVGMSQiK3Q+RkxoRjo3JCQhK2cxW2MpKkY6JCIrNyk0WSZvRjo3JCQhKz1ZM1AiKkY6JCIrQCEzVF8oRjo3JCQhK24mPk9OKUY6JCIrIWUhek0iKUY6NyQkISsjUUtCXihGOiQiK1c8SnonKUY6NyQkISs+OEA/bUY6JCIrNkpDXSIqRjo3JCQhKyZHeFxvJkY6JCIrSkBVUyYqRjo3JCQhK3YhUV1yJUY6JCIrI3BSTCUpKkY6NyQkIStkSlM+UEY6JCIrL09OMDVGMTckJCErXypmdnEjRjokIittb3A7NUYxNyQkISspRzskKm8iRjokIispW0oiPTVGMTckJCErNUc8WW5GNCQiKyVmPCc0NUYxNyQkIistRl9tS0Y0JCIrXylmSiIqKkY6NyQkIitecydbSSJGOiQiK1VPJGVqKkY6NyQkIis9MCYzRCNGOiQiK2U9YnAjKkY6NyQkIitsJ1FnOiRGOiQiKyY0Iz1AKSlGOjckJCIrI1s7RSwlRjokIitzI2YoKUgpRjo3JCQiKzhbaDhbRjokIitAYjU2eEY6NyQkIitCLShIYiZGOiQiK2hAWG5xRjo3JCQiKyJIW2JBJ0Y6JCIrd100eGpGOjckJCIrR144Rm9GOiQiKyVbMSJcY0Y6NyQkIitxTVNhdEY6JCIreWw2IypbRjo3JCQiK2kxJ1sheUY6JCIrOyVwVDYlRjo3JCQiK01leXciKUY6JCIrZipcRUskRjo3JCQiKytJO3AlKUY6JCIrKnp0VV8jRjo3JCQiKyl5OjtvKUY6JCIrISl5N0Q8Rjo3JCQiK1UjW1YiKSlGOiQiKyEpPVwyJCpGNDckJCIrZFs0bykpRjokIitUR0dpOUY0NyQkIishbyIzVykpRjokISs3aXVPaUY0NyQkIit4IioqUnUpRjokISskPTdXUCJGOjckJCIrOWwpKnAmKUY6JCErPjNhLEBGOjckJCIra2ZpQyQpRjokISsiXEcxIUdGOjckJCIreWAlNCwpRjokISt5bD1uTUY6NyQkIitJZ1VLd0Y6JCErZTxqJzQlRjo3JCQiK21ELCQ+KEY6JCErdUFHJW8lRjo3JCQiKy06NyhwJ0Y6JCErJT4iUkRfRjo3JCQiKzxYa1xoRjokISsqKXBAOmRGOjckJCIrLkglZmImRjokISs6IkgiXGhGOjckJCIrQipIPSNcRjokIStuS3VBbEY6NyQkIismXFFORCVGOiQhK0Z6MktvRjo3JCQiKzVFbmRORjokIStCNXN0cUY6NyQkIitBQDlUR0Y6JCErR04qXEMoRjo3JCQiK0A+MzZARjokIStoZjRXdEY6NyQkIitxeXZ1OEY6JCErWFdBcXRGOjckJCIrYzhqJVInRjQkIStiX2tCdEY6NyQkISpYYCFmKClGNCQhKyNbP2Q/KEY6NyQkIStkXyJRKnpGNCQhKyI9IykpPXFGOjckJCErKik+PSpbIkY6JCErdHRibW5GOjckJCErJnpeMTojRjokIStFNDBga0Y6NyQkISskND56eCNGOiQhK3VwUiQzJ0Y6NyQkISshKVxrbExGOiQhK2A7PmpjRjo3JCQhK2U1MjRSRjokIStYL1UpPiZGOjckJCEreSpcU1MlRjokIStPQ0gmcCVGOjckJCErKT1ncSVbRjokISswOzRnVEY6NyQkISt4cD9OX0Y6JCErRkIwKmYkRjo3JCQhK2dDQG1iRjokISskPWojPUlGOjckJCErOnVSUWVGOiQhKyk0K09VI0Y6NyQkIStYeWxdZ0Y6JCErLHRvPz1GOjckJCErIiopUUM/J0Y6JCErQEUpW0AiRjo3JCQhKyFIM1BIJ0Y6JCErNUUjRzYnRjQ3JCQhK2I1JFxLJ0Y6JCErRW0lUloiISM3NyQkIStyay8oSCdGOiQiK114MyxkRjQ3JCQhK0N6VzZpRjokIisqcDcpUTZGOjckJCErSGInKnBnRjokIitbTzMobyJGOjckJCErTTUmWyhlRjokIissYXg1QUY6NyQkISthU3dHY0Y6JCIrOG0oZXEjRjo3JCQhKy4lZVpMJkY6JCIrSXBeb0pGOjckJCEranFFJypcRjokIismPSkpXGYkRjo3JCQhK2xVMzxZRjokIitPMnkiKVJGOjckJCErU0tNLFVGOiQiKyhvTGNLJUY6NyQkIStYJClbYFBGOiQiK2dqZkJZRjo3JCQhK1QxQ3lLRjokIisnPSg0dFtGOjckJCErM2tiIXkjRjokIis3XyxzXUY6NyQkISsnZnliRSNGOiQiKzUpWyg9X0Y6NyQkISsheiJlUTxGOiQiK1lnRjdgRjo3JCQhK2U8Ilw/IkY6JCIrVzo/X2BGOjckJCErST1AKnAnRjQkIitCY3lRYEY6NyQkIStVWy4qUSJGNCQiK0BHJkhGJkY6NyQkIisqPkooSFFGNCQiKydmeWk6JkY6NyQkIittOUgyKilGNCQiKyh6YTQqXEY6NyQkIitkSGp6OEY6JCIrSmh0elpGOjckJCIrMVBEWD1GOiQiK0c4KGVfJUY6NyQkIitTb10kRyNGOiQiKyNlQklCJUY6NyQkIiswOHAhcCNGOiQiK2lgPTBSRjo3JCQiK3MiPU4xJEY6JCIrd0RmWU5GOjckJCIrSDM4KlIkRjokIishZlU7OyRGOjckJCIrSXQ2JnAkRjokIitnJz5bdiNGOjckJCIrZmReXFJGOiQiK2wuaklCRjo3JCQiKzNQIjM7JUY6JCIrXy5iJCo9Rjo3JCQiKzFJJXpLJUY6JCIrL0YpelciRjo3JCQiK1AzRl1XRjokIitNP0YjKSoqRjQ3JCQiK3ojKWVGWEY6JCIrYy5mJVsmRjQ3JCQiKyY0KDRnWEY6JCIrNyFcci0iRjQ3JCQiKzhhUltYRjokIStkeDpeTEY0NyQkIisib2lNXCVGOiQhK18ncEhoKEY0NyQkIitWU2snUiVGOiQhKyV5XkE8IkY6NyQkIishM00nZlVGOiQhKzsqcFhjIkY6NyQkIitYJ2ZXMyVGOiQhK0AkPV0kPkY6NyQkIitBOVl0UUY6JCErVGpiIUcjRjo3JCQiK0dQRkhPRjokISswZE4pZiNGOjckJCIrdT8hW04kRjokIStNYCNlKUdGOjckJCIrRiEpPmBJRjokIStLLWtTSkY6NyQkIitSNSR5cyNGOiQhK19KeGdMRjo3JCQiK0xxRCNRI0Y6JCErR3RfV05GOjckJCIrYVA9Pz9GOiQhKyM9bzBwJEY6NyQkIisnM0xhayJGOiQhKyF6Xnp6JEY6NyQkIisneSs+RSJGOiQhK116OW1RRjo3JCQiK0lEOU4oKUY0JCErRGYwJipRRjo3JCQiKyVlIio9JVtGNCQhK2Q4LCYpUUY6NyQkIitpbmV5KCpGXV1tJCErRE15T1FGOjckJCErTk5MPkdGNCQhK2Y+Y15QRjo3JCQhK0ZOWDhsRjQkISsqKlskNGokRjo3JCQhKzVXKzI1RjokISs1ciZvWiRGOjckJCEra3FuWDhGOiQhK2g/aCJIJEY6NyQkIStNR1BrO0Y6JCErYyNveDIkRjo3JCQhKyhScS4nPkY6JCEra0M4UUdGOjckJCErRCFSN0IjRjokISs5PXF2REY6NyQkIStZXyZbWiNGOiQhK05taCRII0Y6NyQkISsqKltUKm8jRjokIStDYzYmKj5GOjckJCErdjNXdEdGOiQhK1NYXCRvIkY6NyQkIStJcHlESUY6JCErJj1wP08iRjo3JCQhK0whUWM5JEY6JCErJSlIOU01Rjo3JCQhKydwMkRCJEY6JCErdTR6SHFGNDckJCErSEZCJ0ckRjokIStjanc8UEY0NyQkISsqeWxwSSRGOiQhK0tBLmxWRl1dbTckJCErVWc7JkgkRjokIit4QChReSNGNDckJCErUCIpZV5LRjokIismRzlXImZGNDckJCErXidwczwkRjokIislKSozdiMqKUY0NyQkISs9cV50SUY6JCIrIkctKHo2Rjo3JCQhK3IqKik9JUhGOiQiK0ljJSlcOUY6NyQkISstVz0leSNGOiQiK3YvIjRxIkY6NyQkIStzTlQtRUY6JCIrV0YjMyQ+Rjo3JCQhKypmKHkpUiNGOiQiK2VCc1BARjo3JCQhK242cHZARjokIit2NSkqPkJGOjckJCErYiVmYyQ+RjokIislKSo+aVojRjo3JCQhKyNRXzhvIkY6JCIrYWVLMEVGOjckJCErJUhGYlQiRjokIisjW2xrcSNGOjckJCErby4sVDZGOiQiK2tyNHpGRjo3JCQhKy4ubjEnKUY0JCIrdSV5SCNHRjo3JCQhK3MjPVB4JkY0JCIreCdwIlFHRjo3JCQhKzQ9eVJIRjQkIiswQi5ER0Y6NyQkIStfWiYpRzhGXV1tJCIrX0VBJXkjRjo3JCQiK1hfdT5FRjQkIitPLW87RkY6NyQkIisoUWY+SCZGNCQiKzM6aEJFRjo3JCQiKzR2JCpleUY0JCIrRSRwa10jRjo3JCQiK01FdUg1RjokIitCKkdwTyNGOjckJCIrKHolZmU3RjokIislPmhvPyNGOjckJCIrVSIpXHE5RjokIis3VUlHP0Y6NyQkIis3OnJqO0Y6JCIraFFWTD1GOjckJCIrL3FyTz1GOiQiK0lXYEM7Rjo3JCQiK2ckSCMpKT5GOiQiK24pcFJTIkY6NyQkIitaRD88QEY6JCIrZ2k6dTZGOjckJCIrZ08kR0EjRjokIis+JGZgUCpGNDckJCIrSFZjL0JGOiQiKzhRYGxwRjQ3JCQiKzYpej9PI0Y6JCIrRTdMT1hGNDckJCIrZ2ZJJlIjRjokIitvTmY2QEY0NyQkIiteVVMvQ0Y6JCErOixGYUdGXV1tNyQkIiswXXcqUSNGOiQhK2syQ0tFRjQ3JCQiK2UqKSo+TiNGOiQhKzg4RzJcRjQ3JCQiK19zIz5II0Y6JCErcXYwITQoRjQ3JCQiK2UqcDBAI0Y6JCErZScpSGgiKkY0NyQkIit3Sjg0QEY6JCErPyQ0LjYiRjo3JCQiK2VaKiopKT5GOiQhK1hCISoqRyJGOjckJCIrRSgpb149RjokISt4JEhNWCJGOjckJCIrUiUpKSkpcCJGOiQhKzVfZSpmIkY6NyQkIispZSlRSzpGOiQhKyo+XXNzIkY6NyQkIitIaDNhOEY6JCErVy1dTj1GOjckJCIrIkhnZjsiRjokISspPjpPIz5GOjckJCIrNyU9MHEqRjQkIStuYDMiKj5GOjckJCIrayE9V28oRjQkIStpZGhQP0Y6NyQkIitQKilIS2NGNCQhK0F5N2o/Rjo3JCQiK2EoR15jJEY0JCErKlJleDEjRjo3JCQiK3VAai46RjQkISshKmUmPTAjRjo3JCQhK3hpPD5gRl1dbSQhKy1OKGYsI0Y6NyQkIStRIyoqPV8jRjQkIStnKGYzJz5GOjckJCErMERgWldGNCQhK1xrVygpPUY6NyQkISsxWzAiSCdGNCQhKyVITW96IkY6NyQkIStBPyllLilGNCQhK1BuRiFwIkY6NyQkISstRHZtJypGNCQhK0VBO3A6Rjo3JCQhK3dfKXA2IkY6JCErV2EqXFYiRjo3JCQhKy1DSGA3RjokISsrelAqRyJGOjckJCErRlxgdThGOiQhKyVIKSlSOCJGOjckJCEraT8jKXo5RjokISt3KUhjcSpGNDckJCErUyhRJW86RjokISt5J28oMyEpRjQ3JCQhK0wmWylSO0Y6JCErIVtbc0UnRjQ3JCQhK2lacCRwIkY6JCErazgvKlwlRjQ3JCQhK3opKnpIPEY6JCEreFAoPnMjRjQ3JCQhKzJLO1s8RjokISs8LW5OJipGXV1tNyQkISs9ciYqWzxGOiQiKyF5SjYqeUZdXW03JCQhK3g8X0s8RjokIiskSClmKlsjRjQ3JCQhKyZbZSQqcCJGOiQiK3csOktURjQ3JCQhK108N107RjokIitTaCk9cSZGNDckJCErUixoJmUiRjokIisxaCpbPShGNDckJCErJypmdjE6RjokIitsPE5vJilGNDckJCErPlRoOTlGOiQiK3R0Z1MpKkY0NyQkIStnJVwuSiJGOiQiK1siSCIqNCJGOjckJCErc1JBJj4iRjokIituKVE2PyJGOjckJCErI28jZXE1RjokIisnKjRMKkciRjo3JCQhKyVIIlJ5JCpGNCQiK2xdNGo4Rjo3JCQhKytIZyUpekY0JCIrMk0oPlUiRjo3JCQhKzE5XlJsRjQkIisoKlJtbDlGOjckJCErQS1SZV1GNCQiK1M7LSVcIkY6NyQkISsocHpsYiRGNCQiKzV1MDI6Rjo3JCQhK289TFw/RjQkIitBaiRcXSJGOjckJCErTFZlO2JGXV1tJCIrMEsoelsiRjo3JCQiKyl5aSI9IypGXV1tJCIrXW9pYzlGOjckJCIrXWMqcE4jRjQkIis3RFw2OUY6NyQkIitlR1dTUEY0JCIrOUhIYDhGOjckJCIrIUcuJmZdRjQkIis7eidHRyJGOjckJCIrJipmVC1qRjQkIitySTssN0Y6NyQkIitpTFRldUY0JCIrbHNANDZGOjckJCIrJyl5IXleKUY0JCIrIXBcIjM1Rjo3JCQiKzVlMXMlKkY0JCIrR09ZIioqKUY0NyQkIis0blFKNUY6JCIrIylSWU15RjQ3JCQiK2deci42RjokIishUSNHQm1GNDckJCIra15yajZGOiQiKyMpKmY0UCZGNDckJCIrJTNUNUAiRjokIitgMm0hNCVGNDckJCIrJioqeWFDIkY6JCIraTdhJnojRjQ3JCQiK3NmJXBFIkY6JCIrO15pKVwiRjQ3JCQiK2kiKVt2N0Y6JCIrTG4jbzcjRl1dbTckJCIra3dGcjdGOiQhK1NNKSlcNUY0NyQkIis7JDNZRCJGOiQhKz5tOXhBRjQ3JCQiKzgvKmVBIkY6JCErYyFleFgkRjQ3JCQiK1xXayY9IkY6JCErcDsvImUlRjQ3JCQiK11SXE02RjokIStAPTNQY0Y0NyQkIitidTp0NUY6JCErV00hb2gnRjQ3JCQiK10oUkMrIkY6JCErVFMvN3ZGNDckJCIrVlRBSyMqRjQkIStJcWg6JClGNDckJCIrWXlgayQpRjQkIStBZU9AISpGNDckJCIrXVVRSnVGNCQhK2t6P0MnKkY0NyQkIisoM2lLVydGNCQhK3ZvLDc1Rjo3JCQiKyk9NjVUJkY0JCErdS5rXTVGOjckJCIrOCkpcFhWRjQkIStoNDd5NUY6NyQkIismWzAmZUtGNCQhK0NCUyU0IkY6NyQkIitFOmhnQEY0JCErazNhKjQiRjo3JCQiKyFcJjNqNUY0JCErKT4uUDQiRjo3JCQhKkdKQ0sjRl1dbSQhKyVcaXIyIkY6NyQkIStyU3ooMyJGNCQhKyhwJEhdNUY6NyQkIStBTV8/QEY0JCErOXhjODVGOjckJCErOCJHPTYkRjQkIStJbFd2JypGNDckJCErb2FzX1NGNCQhKyZSaSdHIipGNDckJCErNz4iXCRcRjQkIStFXlotJilGNDckJCErJyk9JDN2JkY0JCErYyV6WSF5RjQ3JCQhK1x3dCRcJ0Y0JCErdjZpVnFGNDckJCEreiZSeDooRjQkISt3MTVHaUY0NyQkIStObiV5dChGNCQhKydmdXNPJkY0NyQkISs8dioqSCMpRjQkIStLbmJxV0Y0NyQkIStZLzNKJylGNCQhKzgkPnZhJEY0NyQkISs/XiUqUSopRjQkISsmZiZ6MkVGNDckJCErT05TXyIqRjQkIStIWSk0bSJGNDckJCErQT9BciMqRjQkIStZJiplbHJGXV1tNyQkIStDVDUnSCpGNCQiK1o3Q2lARl1dbTckJCErNVttRyMqRjQkIitJclNHNkY0NyQkIStfa1JyISpGNCQiK1JCVDY/RjQ3JCQhK3JxaUYpKUY0JCIrJGY7ciZHRjQ3JCQhK2g4WiwmKUY0JCIrYTAjemwkRjQ3JCQhKyc0dng0KUY0JCIrISk+Im9TJUY0NyQkIStGTjBBd0Y0JCIrNUVVKDQmRjQ3JCQhKzFXVSEzKEY0JCIraiN5U3MmRjQ3JCQhKyhRTyZ6a0Y0JCIrYzwlPUcnRjQ3JCQhKyJbJFxFZUY0JCIrcSNbbHcnRjQ3JCQhKydRdyhHXkY0JCIrXEckWzwoRjQ3JCQhK2k5OyVSJUY0JCIrUChcVF0oRjQ3JCQhKyQ+UTFqJEY0JCIrX0x5X3hGNDckJCErcWZLWUdGNCQiK3QyJik+ekY0NyQkIStoNVJcP0Y0JCIrS2RIMCEpRjQ3JCQhKyVma3pDIkY0JCIrIj11KTQhKUY0NyQkISs9U2krWEZdXW0kIitzOzhOekY0NyQkIisvWyNcTyRGXV1tJCIrc0dQJHkoRjQ3JCQiKyEqUTkvNkY0JCIraU9pZHZGNDckJCIrdUNuWD1GNCQiKypSJ2Voc0Y0NyQkIitVM0VhREY0JCIrdiopZSoqb0Y0NyQkIitSTGJCS0Y0JCIrRCVHbFonRjQ3JCQiK1IqNHglUUY0JCIrbCk0eSpmRjQ3JCQiK1xRWEBXRjQkIisjZiJHcGFGNDckJCIrU0M2U1xGNCQiKyQ9bnIqW0Y0NyQkIishSF4nKlImRjQkIitRYSp6RyVGNDckJCIrXjZxJ3omRjQkIishKipHJltPRjQ3JCQiK0h4ZEdoRjQkIitlQnAmKUhGNDckJCIrI0clSCRSJ0Y0JCIrJT4rbEkjRjQ3JCQiKyUzbiYqZSdGNCQiKzdrKXpoIkY0NyQkIitsUSJvcidGNCQiK2EkWzhGKkZdXW03JCQiKz9hOXZuRjQkIis3Jz4iM0NGXV1tNyQkIitCMU5sbkY0JCErLGooSE0lRl1dbTckJCIrTF4oKSlvJ0Y0JCErRkZ2IjQiRjQ3JCQiK1JUelpsRjQkIStWTFREPEY0NyQkIistJnpaTSdGNCQhK3dYXEhCRjQ3JCQiK1E8MSQzJ0Y0JCErbmdoKSpHRjQ3JCQiKyFIKFFtZEY0JCErSW4leVUkRjQ3JCQiK3A6KCopUiZGNCQhKzxXdTdSRjQ3JCQiKywkW2EpXEY0JCErTC1SXFZGNDckJCIrPmQiM2AlRjQkIStcc1RNWkY0NyQkIitzLFFTU0Y0JCErKlJMXTEmRjQ3JCQiK1R4cD5ORjQkISsjUlEhUmBGNDckJCIrc1dedUhGNCQhKzlYJFtiJkY0NyQkIitOZXE1Q0Y0JCErJEhKOXImRjQ3JCQiK0BnQE09RjQkISsmKlJXM2VGNDckJCIrKnAoKjREIkY0JCEraGczWWVGNDckJCIrOyNHJnBtRl1dbSQhK0BjOkRlRjQ3JCQiKztzXSh5KSEjOCQhKy5qLFpkRjQ3JCQhKzZdMDFbRl1dbSQhK09FZDhjRjQ3JCQhK0teMkw1RjQkISsjUlNzVSZGNDckJCErI3kmUWs6RjQkISspSDc0PiZGNDckJCErWTRxcD9GNCQhK2omPnohXEY0NyQkIStAdWBXREY0JCErYygqKT5lJUY0NyQkISsyRiFbKUhGNCQhK1Y6PzxVRjQ3JCQhK1J4I29RJEY0JCErcG8keiJRRjQ3JCQhK1daUlpQRjQkIStWNiQpKVEkRjQ3JCQhK0crd2pTRjQkISskKj5zTUhGNDckJCErMTtuTFZGNCQhK3FyZmdDRjQ3JCQhKyhlIlFiWEY0JCErR0Bhcj5GNDckJCErTExsRlpGNCQhK0R5b3M5RjQ3JCQhK0ZNd1xbRjQkISsrQmYicCpGXV1tNyQkIStcJilcQFxGNCQhKzkuRGdZRl1dbTckJCErJVJaSiVcRjQkIis1TDV2SkZnZ3E3JCQhK055WzpcRjQkIis8emMkPiZGXV1tNyQkIStHJnAoUltGNCQiK3EjZTwjKipGXV1tNyQkIStbQXA8WkY0JCIrdylHZVciRjQ3JCQhK0YtUV5YRjQkIitMQEB3PUY0NyQkIStQR05WVkY0JCIrUEBgekFGNDckJCErZz5cJzQlRjQkIis6cUxfRUY0NyQkIStJays5UUY0JCIrO1FhIipIRjQ3JCQhK1ZSUipcJEY0JCIrdzxZJUgkRjQ3JCQhK18zU2NKRjQkIitBOCIpZU5GNDckJCErI0chKSopeSNGNCQiKzgoUUZ5JEY0NyQkISslKipbN1MjRjQkIis/aSNbJ1JGNDckJCErbU1XKCo+RjQkIitAeTQvVEY0NyQkIStEYyg9ZSJGNCQiKzcvLStVRjQ3JCQhK3gnKSkpZTZGNCQiK3AwXV9VRjQ3JCQhK29bOEd0Rl1dbSQiKyYqcCg9RSVGNDckJCErSyJmI3pJRl1dbSQiK2MoMypHVUY0NyQkIiskb1dmNiJGXV1tJCIrTSZmWjolRjQ3JCQiK0MpUXFAJkZdXW0kIisqNHk0LyVGNDckJCIrLXVKJj0qRl1dbSQiKytgWiopUUY0NyQkIitqYlQpSCJGNCQiK2V5XC1QRjQ3JCQiK05KJHpsIkY0JCIrQCQqZiNbJEY0NyQkIitMOyRSKj5GNCQiK0EmM0VCJEY0NyQkIishXGBOSSNGNCQiK1BmZmJIRjQ3JCQiKyl6ZVVlI0Y0JCIrVCRRW2wjRjQ3JCQiKz1UJVEkR0Y0JCIrR0B5TEJGNDckJCIrcXZYXUlGNCQiK1lkKycqPkY0NyQkIitMLGhLS0Y0JCIrNTs9WDtGNDckJCIrTko9ekxGNCQiKyQqeS4mRyJGNDckJCIrRT5WKlskRjQkIiszZT8kPipGXV1tNyQkIisleSYpSGMkRjQkIisiXGh2XiZGXV1tNyQkIitVVCUpKmYkRjQkIit2PzpnPUZdXW03JCQiK0kjcC5nJEY0JCErV2JBVjxGXV1tNyQkIisuYEZsTkY0JCErcyI9IWVfRl1dbTckJCIrLlloJlwkRjQkIStFPkJeJylGXV1tNyQkIisiR2dGUiRGNCQhK0pFPCo9IkY0NyQkIispWydRZUtGNCQhK1clW11cIkY0NyQkIispKmVXJTQkRjQkIStEZDMheSJGNDckJCIrQVw5LkhGNCQhK3d6KD0vI0Y0NyQkIis1cCJwbyNGNCQhKzVlSHlBRjQ3JCQiK0dPUltDRjQkIStpRF0oWyNGNDckJCIrWGBQIT4jRjQkIStObid6bSNGNDckJCIrQCkqejo+RjQkISs1N1o9R0Y0NyQkIitbMHJGO0Y0JCErcyI+IlFIRjQ3JCQiK2lYQUg4RjQkIStjbkxFSUY0NyQkIitdLF1CNUY0JCErVUMoRzMkRjQ3JCQiKyF5WnE4KEZdXW0kIStmS3oySkY0NyQkIitxYiQpSFNGXV1tJCErIyl5WixKRjQ3JCQiKzVaJkdXKkZnZ3EkIStybmdrSUY0NyQkISshNDMmKjMjRl1dbSQhK3clXCIpKkhGNDckJCErKFxlRC8mRl1dbSQhKzYkXEwhSEY0NyQkIStFP0QoKXlGXV1tJCErJ2YwPHkjRjQ3JCQhKy8lZihmNUY0JCErLFAmXGojRjQ3JCQhK0RNJVxKIkY0JCErSkovbENGNDckJCErOXAxXzpGNCQhKzxSNnVBRjQ3JCQhKzk/OXA8RjQkISspZnJXMSNGNDckJCErNSc9Vyc+RjQkIStrNWNRPUY0NyQkISs/cFJPQEY0JCErYCZSKilmIkY0NyQkIStmdCRRRyNGNCQhKyE0XCNbOEY0NyQkISs+d3cwQ0Y0JCErYCcpPSozIkY0NyQkISslKXBbLERGNCQhKyJcaltDKUZdXW03JCQhK1x5Y3FERjQkISs/KD0ob2JGXV1tNyQkIStPWCZHaCNGNCQhKygzJlwhKkdGXV1tNyQkIStDJWYlR0VGNCQhK1BVPm5CRmdncTckJCErR2x2PEVGNCQiK2UhXG9PI0ZdXW03JCQhKz5FUCJlI0Y0JCIraHhYJipbRl1dbTckJCErO2c8P0RGNCQiK1VdZUR0Rl1dbTckJCErRktFTkNGNCQiKzQhW15qKkZdXW03JCQhKydmVnpLI0Y0JCIrNSNwLj0iRjQ3JCQhK2RAcyo+I0Y0JCIrIzRjN1EiRjQ3JCQhKycpM0dfP0Y0JCIrNlhdazpGNDckJCErIWVldSk9RjQkIis/NGxHPEY0NyQkISt4JkhzcSJGNCQiK2EmXEMoPUY0NyQkIStoOm84OkY0JCIrai8pWyo+RjQ3JCQhK3JFKiozOEY0JCIrJz1iXjQjRjQ3JCQhK3UiM2E0IkY0JCIrMCtzc0BGNDckJCErTiNwQHYpRl1dbSQiKz5vREZBRjQ3JCQhK1smeXNdJ0ZdXW0kIitwO29lQUY0NyQkISstdVlVVUZdXW0kIitsNDluQUY0NyQkIStHVmAhKT5GXV1tJCIrLWIrYEFGNDckJCIrJVE9P2MjRmdncSQiK2dCJ29AI0Y0NyQkIitoTDZZQ0ZdXW0kIis9W11mQEY0NyQkIit5UllvWEZdXW0kIis4LyM+MyNGNDckJCIreWtmLm1GXV1tJCIrMnRGJik+RjQ3JCQiK21lNUwmKUZdXW0kIitKIjQ0KD1GNDckJCIrbycqKlIuIkY0JCIrIlspSFM8RjQ3JCQiK0RCKTM/IkY0JCIrdiNmXWYiRjQ3JCQiK3Q2ZV84RjQkIitkIj1wViJGNDckJCIrR2oqeVsiRjQkIitmXnBuN0Y0NyQkIiteUCJlZyJGNCQiK0BRRyozIkY0NyQkIishZjZicSJGNCQiK04rSk8hKkZdXW03JCQiKyxYTyd5IkY0JCIrJ2VscjcoRl1dbTckJCIrR2YleiU9RjQkIis2NUwmPSZGXV1tNyQkIisjUUcrKj1GNCQiK1cwbUlLRl1dbTckJCIrVzxlNz5GNCQiKy08dSNHIkZdXW03JCQiK1sncGQiPkY0JCErInpISFInRmdncTckJCIrQ1MlKioqPUY0JCEra2YlcF4jRl1dbTckJCIrKSl5amw9RjQkISsnKVtfS1ZGXV1tNyQkIisqW2NOIj1GNCQhKzA5SXBnRl1dbTckJCIrQXBjVzxGNCQhK0F1bDZ4Rl1dbTckJCIrY2pvZjtGNCQhK291QFgjKkZdXW03JCQiKy0hcCtjIkY0JCErQm1wbDVGNDckJCIrKSo+KnBXIkY0JCErUl1gJD4iRjQ3JCQiK1AuJT1LIkY0JCErKVJPcUkiRjQ3JCQiK1g1NCc9IkY0JCErJ29oYFMiRjQ3JCQiK1tuSFQ1RjQkISslZVF5WyJGNDckJCIrVCMpbyEqKSlGXV1tJCErUFsnUmIiRjQ3JCQiKyEpSGc1dEZdXW0kISt1KTRNZyJGNDckJCIrWG1bKm8mRl1dbSQhKyFwOWdqIkY0NyQkIiskNCk9V1NGXV1tJCErVSt6XjtGNDckJCIrOz5ZIlIjRl1dbSQhKyIqSCI0bCJGNDckJCIrJEdYIXl1RmdncSQhK2w9c0w7RjQ3JCQhKzsmUnFxKUZnZ3EkIStCKDQyZyJGNDckJCErdzBiW0NGXV1tJCErNmteXzpGNDckJCErTz0kNChSRl1dbSQhK0IkPioqWyJGNDckJCErJTQnKVFVJkZdXW0kISt2SCNRVCJGNDckJCErWV1UJXonRl1dbSQhK052Q0Q4RjQ3JCQhK2BLZnEhKUZdXW0kIStCa0pEN0Y0NyQkIStnW21UIypGXV1tJCErWkpDOjZGNDckJCErJW83KUg1RjQkIStzXDxqKipGXV1tNyQkISs6KXpKNyJGNCQhK1VGIiopcClGXV1tNyQkIStcdWUuN0Y0JCErcG5qdHRGXV1tNyQkIStIb1xxN0Y0JCErPl9tLGdGXV1tNyQkISskUjhOSyJGNCQhK0NWXChmJUZdXW03JCQhKz95UWk4RjQkISspW15jPCRGXV1tNyQkISs6YywoUSJGNCQhKzk1Yl08Rl1dbTckJCErdV9WKFIiRjQkISskSFFOTyRGZ2dxNyQkIStyXSNRUiJGNCQiKyRcbkowIkZdXW03JCQhK0QmKVx3OEY0JCIrJj5kWlMjRl1dbTckJCErPycpKmVNIkY0JCIrUnJ6MFBGXV1tNyQkIStsM2YtOEY0JCIrRio9VyVcRl1dbTckJCErc19EWjdGNCQiKypSMSc0aEZdXW03JCQhK2x0biE9IkY0JCIrJ3AoRyI+KEZdXW03JCQhK0IkUVA1IkY0JCIrIj07Lz0pRl1dbTckJCErclZTPDVGNCQiK0wnUiFwISpGXV1tNyQkIStGYTpGIypGXV1tJCIrenNOXSkqRl1dbTckJCErNyRSeD8pRl1dbSQiK0prKD0wIkY0NyQkISslZjp0NyhGXV1tJCIrciQpKXA1IkY0NyQkIStyKnp4KmZGXV1tJCIrSyZlKzoiRjQ3JCQhKylceTckW0ZdXW0kIis6IyoqMz0iRjQtJSVGT05URzYkJSpIRUxWRVRJQ0FHIiM1

restart;with(plots,listplot):with(plots,display):runge4 := proc(f1, f2, f3, a, b, h, alpha1, alpha2, alpha3) local i, n, au1, au2, au3, k11, k12, k13, k21, k22, k23, k31, k32, k33, k41, k42, k43, tmp1, tmp2, tmp3; n := (b-a)/h; au1 := [alpha1]; au2 := [alpha2]; au3 := [alpha3]; for i from 0 to n-1 do print(i); k11 := evalf( h*f1(a+i*h, op(-1,au1), op(-1,au2), op(-1,au3))); k12 := evalf( h*f2(a+i*h, op(-1,au1), op(-1,au2), op(-1,au3))); k13 := evalf( h*f3(a+i*h, op(-1,au1), op(-1,au2), op(-1,au3))); k21 := evalf( h*f1( a+(i+1/2)*h,op(-1,au1)+k11/2, op(-1,au2)+k12/2, op(-1,au3)+k13/2)); k22 := evalf( h*f2( a+(i+1/2)*h,op(-1,au1)+k11/2, op(-1,au2)+k12/2, op(-1,au3)+k13/2)); k23 := evalf( h*f3( a+(i+1/2)*h,op(-1,au1)+k11/2, op(-1,au2)+k12/2, op(-1,au3)+k13/2)); k31 := evalf( h*f1( a+(i+1/2)*h,op(-1,au1)+k21/2, op(-1,au2)+k22/2, op(-1,au3)+k23/2)); k32 := evalf( h*f2( a+(i+1/2)*h,op(-1,au1)+k21/2, op(-1,au2)+k22/2, op(-1,au3)+k23/2)); k33 := evalf( h*f3( a+(i+1/2)*h,op(-1,au1)+k21/2, op(-1,au2)+k22/2, op(-1,au3)+k23/2)); k41 := evalf( h*f1( a+(i+1)*h,op(-1,au1)+k31, op(-1,au2)+k32,op(-1,au3)+k33)); k42 := evalf( h*f2( a+(i+1)*h,op(-1,au1)+k31, op(-1,au2)+k32, op(-1,au3)+k33)); k43 := evalf( h*f3( a+(i+1)*h,op(-1,au1)+k31, op(-1,au2)+k32, op(-1,au3)+k33)); tmp1:=op(-1,au1)+(1/6)*(k11+2*k21+2*k31+k41); tmp2:=op(-1,au2)+(1/6)*(k12+2*k22+2*k32+k42); tmp3:=op(-1,au3)+(1/6)*(k13+2*k23+2*k33+k43); au1:=[op(au1),tmp1]; au2:=[op(au2),tmp2]; au3:=[op(au3),tmp3]; end do; return [au1,au2,au3]; end proc:c1 := 0.2: c2 := 0.2: c3 := 5.7: f1 := (t,x,y,z) -> -y-z; f2 := (t,x,y,z) -> x+x1*y; f2 := (t,x,y,z) -> c2+y*(x-c3);

NiM+SSNmMUc2ImYqNiZJInRHRiVJInhHRiVJInlHRiVJInpHRiVGJTYkSSlvcGVyYXRvckdGJUkmYXJyb3dHRiVGJSwmOSYhIiI5J0YxRiVGJUYlNiM+SSNmMkc2ImYqNiZJInRHRiVJInhHRiVJInlHRiVJInpHRiVGJTYkSSlvcGVyYXRvckdGJUkmYXJyb3dHRiVGJSwmOSUiIiIqJkkjeDFHRiVGMTkmRjFGMUYlRiVGJQ==NiM+SSNmMkc2ImYqNiZJInRHRiVJInhHRiVJInlHRiVJInpHRiVGJTYkSSlvcGVyYXRvckdGJUkmYXJyb3dHRiVGJSwmSSNjMkdGJSIiIiomOSZGMSwmOSVGMUkjYzNHRiUhIiJGMUYxRiVGJUYl

a:=0;

NiM+SSJhRzYiIiIh

b:=100;

NiM+SSJiRzYiIiQrIg==

n:=10;

NiM+SSJuRzYiIiM1

h:=(b-a)/n;

NiM+SSJoRzYiIiM1

n:=(b-a)/h;

NiM+SSJuRzYiIiM1

au := runge4(f1, f2, f3 ,a,b,h,0,0,0);

NiMiIiE=NiMiIiI=NiMiIiM=NiMiIiQ=NiMiIiU=

NiMiIiY=

Warning, computation interrupted

au1 := [0]: au2 := [0]: au3 := [0]:for i from 0 to n-1 do print(i); k11 := evalf( h*f1(a+i*h, Last(au1), Last(au2), Last(au3))); k12 := evalf( h*f2(a+i*h, Last(au1), Last(au2), Last(au3))); k13 := evalf( h*f3(a+i*h, Last(au1), Last(au2), Last(au3))); k21 := evalf( h*f1( a+(i+1/2)*h, Last(au1)+k11/2, Last(au2)+k12/2, Last(au3)+k13/2)); k22 := evalf( h*f2( a+(i+1/2)*h, Last(au1)+k11/2, Last(au2)+k12/2, Last(au3)+k13/2)); k23 := evalf( h*f3( a+(i+1/2)*h, Last(au1)+k11/2, Last(au2)+k12/2, Last(au3)+k13/2)); k31 := evalf( h*f1( a+(i+1/2)*h, Last(au1)+k21/2, Last(au2)+k22/2, Last(au3)+k23/2)); k32 := evalf( h*f2( a+(i+1/2)*h, Last(au1)+k21/2, Last(au2)+k22/2, Last(au3)+k23/2)); k33 := evalf( h*f3( a+(i+1/2)*h, Last(au1)+k21/2, Last(au2)+k22/2, Last(au3)+k23/2)); k41 := evalf( h*f1( a+(i+1)*h, Last(au1)+k31, Last(au2)+k32, Last(au3)+k33)); k42 := evalf( h*f2( a+(i+1)*h, Last(au1)+k31, Last(au2)+k32, Last(au3)+k33)); k43 := evalf( h*f3( a+(i+1)*h, Last(au1)+k31, Last(au2)+k32, Last(au3)+k33)); au1 := Append(au1, Last(au1)+(1/6)*(k11+2*k21+2*k31+k41)); au2 := Append(au2, Last(au2)+(1/6)*(k12+2*k22+2*k32+k42)); au3 := Append(au3, Last(au3)+(1/6)*(k13+2*k23+2*k33+k43)); end do: au1; au2; au3;

i:=(88.9-a)/h+1;

NiM+SSJpRzYiJCIlKyopISIi

op(890,au[1]);op(890,au[2]);op(890,au[3]);

plist := []: for i from 1 to n do tmp := [op(i,au[1]), op(i,au[2])]; plist := [op(plist), tmp]; end do:listplot(plist);

LSUlUExPVEc2JC0lJ0NVUlZFU0c2IzdkaW43JCQiJSs1IiIhJCIlKz9GLDckJCIrdDkuTzYhIickIitrX24iPSNGMjckJCIrd2kib0kiRjIkIishZnUhM0JGMjckJCIrbiNvUF4iRjIkIisvPUNeQkYyNyQkIityKT5DdiJGMiQiK00jMzBII0YyNyQkIit1Xzc1P0YyJCIrbDJwQUBGMjckJCIra3dnbUFGMiQiK1goZXknPUYyNyQkIisob3QnKVwjRjIkIisoPmxZYyJGMjckJCIrOzlBKG8jRjIkIitwc1FjN0YyNyQkIit1ZEtBR0YyJCIrWj1ZbSgqISIoNyQkIisjZnpQIUhGMiQiKy5Qekh1RmduNyQkIis4JCo0UUhGMiQiKy1dJnplJkZnbjckJCIreSIpcE1IRjIkIis/LlYhPiVGZ243JCQiKyRRcUkhSEYyJCIrZSM+UjokRmduNyQkIisuSVBeR0YyJCIrTislUlIjRmduNyQkIitFLy0neSNGMiQiK1xxa1E9RmduNyQkIitvMSE9ciNGMiQiKyNbKjRLOUZnbjckJCIrRCR5QGojRjIkIispWyp5SzZGZ243JCQiK1dVaFxERjIkIitldDoyIiohIik3JCQiK09tJGVZI0YyJCIrX0l2V3VGZXE3JCQiK2opXD9RI0YyJCIrdUJ0KT0nRmVxNyQkIitldTMqSCNGMiQiKykpKXk2QiZGZXE3JCQiK3k4XzxBRjIkIitHKG9eXCVGZXE3JCQiKzMhUng4I0YyJCIrKCo+Y0RSRmVxNyQkIitlcioqZj9GMiQiKzxNZCNbJEZlcTckJCIrXypmVyk+RjIkIitARklQSkZlcTckJCIrd2FBNj5GMiQiKzl1aW9HRmVxNyQkIitUWU1TPUYyJCIrWnIuaEVGZXE3JCQiKz9ZJD14IkYyJCIrKWU9Sl0jRmVxNyQkIisqWylvMDxGMiQiKzk3XCdRI0ZlcTckJCIrMEEpPWsiRjIkIit6dzAwQkZlcTckJCIrMih6LmUiRjIkIisib3ZXRCNGZXE3JCQiK2ZtOEA6RjIkIitJeHlKQUZlcTckJCIrX0s1azlGMiQiKzFhO05BRmVxNyQkIisuakE0OUYyJCIrUCZRUEUjRmVxNyQkIitoM1hjOEYyJCIrXCQqWzxCRmVxNyQkIitmOXMwOEYyJCIrZCUpPihSI0ZlcTckJCIrJDQkKXBEIkYyJCIrY1dWL0RGZXE3JCQiK1I/PTU3RjIkIisjZiNlVEVGZXE3JCQiK1JsRWw2RjIkIitzMyI+IkdGZXE3JCQiKzJ1PUE2RjIkIis/V24+SUZlcTckJCIrLScpKjMzIkYyJCIrRDVFcUtGZXE3JCQiKyV6ZDgvIkYyJCIrKipIUXFORmVxNyQkIitdcF8uNUYyJCIrN0RLR1JGZXE3JCQiKzApSFBuKkZnbiQiKz4qUlVOJUZlcTckJCIrLVpvRyQqRmduJCIrT3FiZ1tGZXE3JCQiKyhIRCoqKiopRmduJCIrQitWaWFGZXE3JCQiK0IiPXRvKUZnbiQiK1slSCR5aEZlcTckJCIrXXUiM1IpRmduJCIrPiJbMi4oRmVxNyQkIismUSFbNSIpRmduJCIrVCwxWiEpRmVxNyQkIitdYltZeUZnbiQiK25wYmcjKkZlcTckJCIrO2s6KmYoRmduJCIrPihvNjIiRmduNyQkIitgMSoqb3RGZ24kIis9RCZcQyJGZ243JCQiK2R0cGNyRmduJCIrRyhbTFgiRmduNyQkIis2WUNqcEZnbiQiKzMrYC48RmduNyQkIiskeT8qKnknRmduJCIrTCg+VCsjRmduNyQkIitbVFRRbUZnbiQiKyFRRmFPI0ZnbjckJCIrQ3MiNF4nRmduJCIrKD0lbyp6I0ZnbjckJCIrTllFNWtGZ24kIitfLkRATEZnbjckJCIrO242U2pGZ24kIitFUnVZUkZnbjckJCIrO2Q/MGpGZ24kIitkITRdcCVGZ243JCQiKyZlbjtKJ0ZnbiQiKyZRI3onZSZGZ243JCQiK3kpKltuakZnbiQiKzNhI1JrJ0ZnbjckJCIrIVs5SlsnRmduJCIrW0ByKCl5RmduNyQkIis4JVJBbidGZ24kIisyWjRPJCpGZ243JCQiK1QnZUcmcEZnbiQiKyUqRyEqKjQiRjI3JCQiKyQ0ViZbdEZnbiQiK2IsOShHIkYyNyQkIitCJmUpKil5RmduJCIraicqRyNcIkYyNyQkIityJ2ViaClGZ24kIitqclozPEYyNyQkIitrZG1zJipGZ24kIitAbl1CPkYyNyQkIisvXU8iMyJGMiQiKz1jcT1ARjI3JCQiK1QveVE3RjIkIitcWzRwQUYyNyQkIitFJ2ZCViJGMiQiK3YnempNI0YyNyQkIis6YjxnO0YyJCIrODsoZksjRjI3JCQiK0lweDc+RjIkIisrNDMoPiNGMjckJCIrRSkzQzwjRjIkIisjPiw1KD5GMjckJCIrJGUqPTtDRjIkIitFYnIhbyJGMjckJCIrPyxxQUVGMiQiKy0lZidwOEYyNyQkIit4eEV5RkYyJCIrKmZIajIiRjI3JCQiKyczcyN6R0YyJCIrTTw5VyMpRmduNyQkIispeXcsJEhGMiQiK2g/Pj9pRmduNyQkIisvdHJSSEYyJCIrXFxobFlGZ243JCQiK2pMZjxIRjIkIitFZV8vTkZnbjckJCIrPTZic0dGMiQiK2pZWl1FRmduNyQkIitOMm02R0YyJCIrSm01RT9GZ243JCQiKzJxR1NGRjIkIitybWJwOkZnbjckJCIrOixOaUVGMiQiK2R6Qk03RmduNyQkIiskKTNsIWUjRjIkIis+YzVpKSpGZXE3JCQiK3JNOyhcI0YyJCIrPFpyNiEpRmVxNyQkIisvdEU4Q0YyJCIrKik+Wz1tRmVxNyQkIiswJz4qSEJGMiQiK2EqNClmYkZlcTckJCIrZSZ6eEMjRjIkIishcGAlW1pGZXE3JCQiK1BoSG5ARjIkIitGYCxBVEZlcTckJCIrdyhwKCkzI0YyJCIrbGFlTk9GZXE3JCQiKzJjUjc/RjIkIisjW0FtRCRGZXE3JCQiK0pRSFE+RjIkIisqKVJSaEhGZXE3JCQiK3Y2YG09RjIkIitCJipcS0ZGZXE3JCQiK0BoOCh6IkYyJCIrTypScmIjRmVxNyQkIis7KTQsdCJGMiQiK2ROIWZVI0ZlcTckJCIrTlJWbDtGMiQiKykpPSI+TCNGZXE3JCQiKzZrMi47RjIkIiskXDwtRiNGZXE3JCQiK0tjKkhhIkYyJCIrZzJRUEFGZXE3JCQiK2dOOSZbIkYyJCIrKFt2NkIjRmVxNyQkIitsIW8lSDlGMiQiKy4iKlFdQUZlcTckJCIrKnA5ZlAiRjIkIisqUiVwJUgjRmVxNyQkIitrIUdXSyJGMiQiK150ZGtCRmVxNyQkIisjKUcmXEYiRjIkIitJNkpoQ0ZlcTckJCIrcVtWRjdGMiQiK0YwKHBlI0ZlcTckJCIrajkjPT0iRjIkIiszVFtXRkZlcTckJCIrXUMxUTZGMiQiK2tTc1BIRmVxNyQkIitgMDYnNCJGMiQiKzNmanJKRmVxNyQkIisjKj4jZjAiRjIkIitLQlRfTUZlcTckJCIrZnFYPDVGMiQiKy93cih5JEZlcTckJCIrUnUhbyEpKkZnbiQiK1M2KHA9JUZlcTckJCIrUFVqYiUqRmduJCIrYDdxaFlGZXE3JCQiK25tIjM3KkZnbiQiKzVGKWZBJkZlcTckJCIrST0+LSkpRmduJCIrdVgocCplRmVxNyQkIismb3onKlwpRmduJCIrcnhkJnAnRmVxNyQkIitGcUg4IylGZ24kIitAW0Jad0ZlcTckJCIrP1o8VnpGZ24kIitMbipHeSlGZXE3JCQiK2QieiYqbyhGZ24kIitvJD1TLCJGZ243JCQiK3MiUkhYKEZnbiQiK2NsWnc2RmduNyQkIisoMyEpUUIoRmduJCIrSyUpPnI4RmduNyQkIit0cUVMcUZnbiQiKyUqRyhbZyJGZ243JCQiK20zRV9vRmduJCIrQV9iJik9RmduNyQkIis4KipRI3AnRmduJCIra3kqR0EjRmduNyQkIis6WWtibEZnbiQiKyNbKVFHRUZnbjckJCIrXjlnV2tGZ24kIitEQWQ6SkZnbjckJCIrVXllaWpGZ24kIisjUTwtcSRGZ243JCQiK3hFIVJKJ0ZnbiQiKy9lTitXRmduNyQkIitLSTYvakZnbiQiK3d0M09fRmduNyQkIittWVhTakZnbiQiK0FlKCpHaUZnbjckJCIrJSo+UEtrRmduJCIrUGB4K3VGZ243JCQiKyJRU0FmJ0ZnbiQiK1czNnIoKUZnbjckJCIrVlFKT29GZ24kIisjZmZgLiJGMjckJCIrKj4/Zj0oRmduJCIrb1koXEAiRjI3JCQiKyh6Xilvd0ZnbiQiKyNmLlNUIkYyNyQkIityL3I/JClGZ24kIissZTxGO0YyNyQkIisiM2JiPSpGZ24kIispR1tXJT1GMjckJCIraFZXSjVGMiQiKzNUa1w/RjI3JCQiKyczUGY8IkYyJCIrM1QoKj5BRjI3JCQiK0MkR2ZOIkYyJCIrJj4qeUZCRjI3JCQiK1VFaXI6RjIkIitNSS9ZQkYyNyQkIitNYmk7PUYyJCIrL1llZEFGMjckJCIrKFJYaDIjRjIkIis1PSNcMSNGMjckJCIrODdkR0JGMiQiK1V0dCR6IkYyNyQkIitldzNeREYyJCIrVzFQJlsiRjI3JCQiKyZbLm1zI0YyJCIrRWF4Ij0iRjI3JCQiKzkiKXlaR0YyJCIrZ2M6RiIqRmduNyQkIislPStsIkhGMiQiK09AQDxwRmduNyQkIis8K0xTSEYyJCIrLCkqKVs+JkZnbjckJCIrKyY9KkdIRjIkIitzIjRzKlFGZ243JCQiK0FaXSIqR0YyJCIrPXFUUUhGZ243JCQiK1tNImUkR0YyJCIrYDRbT0FGZ243JCQiK0NQI3l3I0YyJCIrWz5iQjxGZ243JCQiKzhuJj5wI0YyJCIrNndlWjhGZ243JCQiK0VPUjZFRjIkIitiLEhxNUZnbjckJCIrQGtRR0RGMiQiK14zKDRrKUZlcTckJCIrKj44WFcjRjIkIitGZiJRNChGZXE3JCQiK1A9KDNPI0YyJCIrenQ4QWZGZXE3JCQiK0wyQXlBRjIkIitjWiRwLSZGZXE3JCQiK2kpeXE+I0YyJCIrTyVbdkwlRmVxNyQkIisiPSd6PEBGMiQiKyxvQy5RRmVxNyQkIitvSGdTP0YyJCIrLT5PKFEkRmVxNyQkIityV2tsPkYyJCIrOFJAaklGZXE3JCQiKzBgKyQqPUYyJCIrI1x3NyJHRmVxNyQkIitebnNBPUYyJCIrWiQqPjxFRmVxNyQkIistIj5bdiJGMiQiK1hMVXFDRmVxNyQkIitLMEYqbyJGMiQiKy5wP2pCRmVxNyQkIitZTjBFO0YyJCIrQCpSKipHI0ZlcTckJCIrJG9IXmMiRjIkIitGMm1ZQUZlcTckJCIrJipIWDE6RjIkIitHJD0yQiNGZXE3JCQiK15FKCpcOUYyJCIrVTBhU0FGZXE3JCQiKyQqW2omUiJGMiQiKzk3W3ZBRmVxNyQkIitUWFFWOEYyJCIrTCZHZEwjRmVxNyQkIitqaDskSCJGMiQiKydwaUFVI0ZlcTckJCIrPV0jXEMiRjIkIitBZCZvYCNGZXE3JCQiK0V5Zyk+IkYyJCIrJm80QG8jRmVxNyQkIitLTTthNkYyJCIrJFFLOidHRmVxNyQkIispUlY6NiJGMiQiKyR6Xid6SUZlcTckJCIrIWUtMjIiRjIkIitbRzxVTEZlcTckJCIrJm8qZko1RjIkIis/QD1jT0ZlcTckJCIraCh5PiUqKkZnbiQiKy8jMy8uJUZlcTckJCIrJVJhWWUqRmduJCIrISl5YXZXRmVxNyQkIituYHdWIypGZ24kIislb3FZKyZGZXE3JCQiK3Y0Nz4qKUZnbiQiKyJIPFBqJkZlcTckJCIrYCkzMWgpRmduJCIrQys2I1EnRmVxNyQkIitIej89JClGZ24kIitJb1x0c0ZlcTckJCIrXFYrVSEpRmduJCIrYkdsTyQpRmVxNyQkIitLN0AjeShGZ24kIisoeWxsZypGZXE3JCQiKyJSJD5SdkZnbiQiK2dKZDc2RmduNyQkIishXClcOHRGZ24kIit4JW9YSCJGZ243JCQiK1VtKmU1KEZnbiQiKygpeShHXiJGZ243JCQiKzI2VjxwRmduJCIrXzstdjxGZ243JCQiK2NQW1xuRmduJCIrViZIKzQjRmduNyQkIiswLidRZydGZ24kIitNInknb0NGZ243JCQiK0pBIUhbJ0ZnbiQiK3BEc0JIRmduNyQkIitwXGoqUSdGZ24kIitMJiozcU1GZ243JCQiK09jJ3pLJ0ZnbiQiKyJldVw3JUZnbjckJCIrcSZbSEknRmduJCIrJHldeCFcRmduNyQkIitWTTlAakZnbiQiK3A8YVJlRmduNyQkIisoZik0IlInRmduJCIrIXBNQSVwRmduNyQkIitRLCtDbEZnbiQiK3NFZ08jKUZnbjckJCIrVDdgTW5GZ24kIitrUDJSKCpGZ243JCQiKzRiKD4vKEZnbiQiK3dYbFg2RjI3JCQiKypIXDpaKEZnbiQiK3VzJnlMIkYyNyQkIitCOSJlMClGZ24kIitFNWtZOkYyNyQkIispWyZvTikpRmduJCIrN2IiUnciRjI3JCQiK1tZLGcpKkZnbiQiK0xdImYoPkYyNyQkIisieWQiPTZGMiQiK0BHRWlARjI3JCQiK3NPayVHIkYyJCIrIVwybkgjRjI3JCQiK184UChbIkYyJCIrZjMyXkJGMjckJCIrOUJzQTxGMiQiKyV5IUcuQkYyNyQkIiskPXMhej5GMiQiKyFHQHc5I0YyNyQkIisnRygpb0IjRjIkIitLb1MsPkYyNyQkIissJG9IWiNGMiQiK21Ubyw7RjI3JCQiK0spKlFuRUYyJCIrcSQ0P0giRjI3JCQiK0FwLDRHRjIkIitZNW4yNUYyNyQkIitpIylmJypHRjIkIityPFYibyhGZ243JCQiKzRmME9IRjIkIitDN1YjeSZGZ243JCQiK0U8ak9IRjIkIisuRkJPVkZnbjckJCIrR3gheiFIRjIkIitRXFJoS0ZnbjckJCIrb1pCZUdGMiQiKzFDZnNDRmduNyQkIitnUkIleiNGMiQiKyxyPCcqPUZnbjckJCIrIyo+KDNzI0YyJCIrZiZbVloiRmduNyQkIitkTHZURUYyJCIrQzAuazZGZ243JCQiK0RHV2ZERjIkIitXNTpTJCpGZXE3JCQiKzxPdXZDRjIkIisrYzk/d0ZlcTckJCIrVU0iPlIjRjIkIis0RyxBakZlcTckJCIrR04jKTNCRjIkIitTI29MTCZGZXE3JCQiK0hDMkZBRjIkIitwaTl1WEZlcTckJCIrRmoxWkBGMiQiKygqZipwKVJGZXE3JCQiK2hgMnA/RjIkIithJXAwYCRGZXE3JCQiK1BKRiQqPkYyJCIrZlYmWzwkRmVxNyQkIis7WHc+PkYyJCIrXSpIeipHRmVxNyQkIitjX2dbPUYyJCIrZG5xJG8jRmVxNyQkIissZyIpejxGMiQiK0ZnTT9ERmVxNyQkIisxQlI4PEYyJCIrNys6KlIjRmVxNyQkIisoKT1KXDtGMiQiK3MveThCRmVxNyQkIisyKlJ2ZSJGMiQiK2d6cWZBRmVxNyQkIituRy5HOkYyJCIrL2EjUUIjRmVxNyQkIistOnVxOUYyJCIrMDA9TUFGZXE3JCQiK1xIaDo5RjIkIiszZ3pmQUZlcTckJCIrNENmaThGMiQiK0xwYjVCRmVxNyQkIitGV2k2OEYyJCIrPjQ6KFEjRmVxNyQkIis1UmxpN0YyJCIrYEgwIlwjRmVxNyQkIitkcGk6N0YyJCIrYSZcWGkjRmVxNyQkIitdOlxxNkYyJCIrby8hM3ojRmVxNyQkIit5Iik+RjZGMiQiKzMlUVIqSEZlcTckJCIrMS9xJjMiRjIkIitQJjQjUktGZXE3JCQiK1RhJmYvIkYyJCIrYyZmSmAkRmVxNyQkIitBWiN6KyJGMiQiKz1CKFEpUUZlcTckJCIrZGF1OigqRmduJCIrIj5xN0klRmVxNyQkIiteendvJCpGZ24kIitReVsoeiVGZXE3JCQiKzJ4NFEhKkZnbiQiK0lQTChRJkZlcTckJCIrcT1mQigpRmduJCIrQ20mKSkzJ0ZlcTckJCIrZFA+RCUpRmduJCIrTHMuQ3BGZXE3JCQiKyRlWkg5KUZnbiQiKzsjRyc+ekZlcTckJCIrJ287cSh5RmduJCIrUyFvIjMiKkZlcTckJCIrV2lxRndGZ24kIitWJj5IMCJGZ243JCQiK1I9XCZSKEZnbiQiKyIqKm9JQSJGZ243JCQiK11jMCI9KEZnbiQiKz9lMkY5RmduNyQkIitfRUwmKXBGZ24kIis4dCY+biJGZ243JCQiKyQ0cSY0b0ZnbiQiKyhmYWgnPkZnbjckJCIrI28vYWwnRmduJCIrTG53PkJGZ243JCQiK05PJ1xfJ0ZnbiQiK1Z1eVdGRmduNyQkIitdISkqNFUnRmduJCIrRHdKYktGZ243JCQiKyhRZ3FNJ0ZnbiQiKy9vcW5RRmduNyQkIit3QXUyakZnbiQiKzBkYStZRmduNyQkIitWWyoqM2pGZ24kIisoXCZRdWFGZ243JCQiK1M2Y2VqRmduJCIrKT4iKjReJ0ZnbjckJCIrWSVlbFknRmduJCIrbDslPXQoRmduNyQkIitsVkRZbUZnbiQiK3BbWWIiKkZnbjckJCIrKkd4XSJwRmduJCIreGNJejVGMjckJCIrIT5yZUgoRmduJCIrIz1zVEUiRjI3JCQiK24oMyQ9eUZnbiQiK1E8Wm45RjI3JCQiKz9GOT8mKUZnbiQiK3RlJkdvIkYyNyQkIishMyhbWiUqRmduJCIraV0jKSkqPUYyNyQkIisjXFlfMSJGMiQiK19mXCg0I0YyNyQkIitAU2E9N0YyJCIrcWxhYUFGMjckJCIrdlgmeVMiRjIkIitrY3pUQkYyNyQkIis4Vic+aiJGMiQiK3R3dExCRjI3JCQiKyk9L0MpPUYyJCIrXV5XPEFGMjckJCIrRmNKVUBGMiQiK0psVSw/RjI3JCQiKyYqRzcqUSNGMiQiK0okPWxyIkYyNyQkIiszI28zZyNGMiQiKyZmS2RTIkYyNyQkIisjKjR3aUZGMiQiK1F5IikzNkYyNyQkIitkKlIrKEdGMiQiK0w/JVJeKUZnbjckJCIrQXlTRUhGMiQiKyNcKCo+VidGZ243JCQiK3U7SFNIRjIkIit0PyRmI1tGZ243JCQiK0ZGU0BIRjIkIityWkhCT0ZnbjckJCIrOEpseUdGMiQiK050YVBGRmduNyQkIitVKT0kPkdGMiQiKzAobygqMyNGZ243JCQiK3VMJipbRkYyJCIrKSo0QTs7RmduNyQkIitYIUg7biNGMiQiK0x3a283RmduNyQkIitGOkUhZiNGMiQiK2wnKXk2NUZnbjckJCIrNiczcF0jRjIkIitnJSplLiMpRmVxNyQkIitfJTRJVSNGMiQiKy1VeWpuRmVxNyQkIitgTmNSQkYyJCIrJSkpZTNuJkZlcTckJCIrSipmc0QjRjIkIis4JFJTJFtGZXE3JCQiK0cnb2w8I0YyJCIrOWVXKT0lRmVxNyQkIisoKlwheTQjRjIkIitqIz11byRGZXE3JCQiK29aPEA/RjIkIitrRTwoSCRGZXE3JCQiKzRkIW8lPkYyJCIrb3kzJCpIRmVxNyQkIitfKnBaKD1GMiQiKypIQHJ2I0ZlcTckJCIrYSkqNDA9RjIkIitAVyt3REZlcTckJCIrSyQqelA8RjIkIis6dygqUkNGZXE3JCQiK0s/JkduIkYyJCIrZCshPk0jRmVxNyQkIityc0E1O0YyJCIrbUJpd0FGZXE3JCQiKyRSJSkpXDpGMiQiK1lgYVNBRmVxNyQkIiskKmZ4IlwiRjIkIit4Q0pKQUZlcTckJCIrSy4mZVYiRjIkIis9L2daQUZlcTckJCIrb0owI1EiRjIkIislRyQpKilHI0ZlcTckJCIraiJILkwiRjIkIisvWCZlTiNGZXE3JCQiKykpSGkhRyJGMiQiKzJlUlxDRmVxNyQkIis4LSlHQiJGMiQiKz9yZXJERmVxNyQkIitbIVtxPSJGMiQiKypRYV9zI0ZlcTckJCIrJilmMlY2RjIkIisnKVw6OUhGZXE3JCQiKyVROzQ1IkYyJCIrMUU1VkpGZXE3JCQiK01dX2c1RjIkIitJZDg9TUZlcTckJCIrTTwnPS0iRjIkIisjeU5udSRGZXE3JCQiKyMqKjMqWykqRmduJCIrV0w1UVRGZXE3JCQiKypcQWVcKkZnbiQiK2JKXS5ZRmVxNyQkIishKmY2ZiIqRmduJCIrKD4icGNeRmVxNyQkIis/eWhRKSlGZ24kIisyK1c5ZUZlcTckJCIrcnhCTSYpRmduJCIrNzg9KGYnRmVxNyQkIistMilmQylGZ24kIiskMyl5SHZGZXE3JCQiK2tMJ1IoekZnbiQiK0syX1UnKUZlcTckJCIrO2BWPXhGZ24kIit1UjtzKipGZXE3JCQiK1BeISl6dUZnbiQiK2whUWo6IkZnbjckJCIrWk5uZXNGZ24kIis2MC5aOEZnbjckJCIramIoZTAoRmduJCIrJikpUWVkIkZnbjckJCIrMVlgc29GZ24kIitqIVwxJj1GZ243JCQiKydcSyxyJ0ZnbiQiK0ZtIjQ9I0ZnbjckJCIrLDNncWxGZ24kIis9XCF6ZCNGZ243JCQiK0s1V2NrRmduJCIrPjIiXDAkRmduNyQkIitFVCkzUCdGZ24kIisnPVJ1aSRGZ243JCQiK0RPNT1qRmduJCIrc0tFOFZGZ243JCQiK1JCXS5qRmduJCIrbUVFS15GZ243JCQiKyE9KjRNakZnbiQiKzR6KGU1J0ZnbjckJCIrRDQwPmtGZ24kIisneUlmRChGZ243JCQiKyQ0Xi5kJ0ZnbiQiKz4oW0NnKUZnbjckJCIraydcUCFvRmduJCIrciQqKmYsIkYyNyQkIitiKj4qUnJGZ24kIismZSc9JD4iRjI3JCQiK25OJGVnKEZnbiQiKzY7OiFSIkYyNyQkIitJRjhPIylGZ24kIishR3g/ZyJGMjckJCIrIT02UjIqRmduJCIrY3NhPj1GMjckJCIrLi4ncCwiRjIkIitRWDlGP0YyNyQkIit2OGNkNkYyJCIreUoiRz8jRjI3JCQiKzM1TUw4RjIkIitgQEY+QkYyNyQkIismKWUyWDpGMiQiKzUiKSoqW0JGMjckJCIrN2NFKHkiRjIkIis+Uzt0QUYyNyQkIitRJyozWT9GMiQiKyFIXTs0I0YyNyQkIitHbl0rQkYyJCIrbitqRj1GMjckJCIrTHlaRkRGMiQiK2QuTEA6RjI3JCQiK21XJyozRkYyJCIrYmpTOjdGMjckJCIrMyllayRHRjIkIituLDQ5JSpGZ243JCQiKyxlITQiSEYyJCIrT2RrWXJGZ243JCQiK2VcVlJIRjIkIit4cWRxYEZnbjckJCIrKVEsOyRIRjIkIisybEBHU0ZnbjckJCIrO3Z3JypHRjIkIisjUU9aLiRGZ243JCQiK0thJkclR0YyJCIrYFojcEkjRmduNyQkIisoKSpSZ3gjRjIkIitPVzZ2PEZnbjckJCIrJEgxNHEjRjIkIitZU14mUSJGZ243JCQiK0A9dz9FRjIkIitcUlIpNCJGZ243JCQiK3g5JnpgI0YyJCIreHExXikpRmVxNyQkIitOLDdhQ0YyJCIrPUVPX3NGZXE3JCQiK15JVHFCRjIkIitybCFILydGZXE3JCQiK3g/aShHI0YyJCIrcWd0Pl5GZXE3JCQiKywpKkcxQUYyJCIrNmRTNFdGZXE3JCQiK3VHeUVARjIkIitKPEFmUUZlcTckJCIrSU1NXD9GMiQiK0JzN0pNRmVxNyQkIitlXjd1PkYyJCIrOD5YKDQkRmVxNyQkIisyKD03IT5GMiQiK1ptJnokR0ZlcTckJCIreSdwMSQ9RjIkIitfQ3hQRUZlcTckJCIrYDlcaTxGMiQiKy55JmZbI0ZlcTckJCIrPFduJ3AiRjIkIiteWlt1QkZlcTckJCIrSUU+TDtGMiQiK25YXyhII0ZlcTckJCIrOSgzP2QiRjIkIisvWSU0RCNGZXE3JCQiK2J1Mjg6RjIkIitrZyY+QiNGZXE3JCQiKyZbW2pYIkYyJCIrazkoKVFBRmVxNyQkIis5JG88UyJGMiQiKyZ5VjRGI0ZlcTckJCIrcT1HXDhGMiQiK1glbyNHQkZlcTckJCIrNVAkKSlIIkYyJCIrQGR0NkNGZXE3JCQiKyIpKm8uRCJGMiQiK3FmLUJERmVxNyQkIiszVSRRPyJGMiQiK2JOa2tFRmVxNyQkIitwejxmNkYyJCIrLiMpKSpSR0ZlcTckJCIrJ1xeajYiRjIkIissXFlgSUZlcTckJCIrVSM0YDIiRjIkIisuY2o1TEZlcTckJCIrYiU0Zy4iRjIkIis0Ij0lPU9GZXE3JCQiK053OSUpKipGZ24kIitUJ1JgKVJGZXE3JCQiK3ckR1xpKkZnbiQiK1FcJj1VJUZlcTckJCIrSjg8I0cqRmduJCIrKT1UMiVcRmVxNyQkIitcd25iKilGZ24kIitoW2VkYkZlcTckJCIrKlxEYGspRmduJCIrUSV5OEgnRmVxNyQkIitYRjNeJClGZ24kIis9IVtfOyhGZXE3JCQiKyhbQUkyKUZnbiQiK19jTDIjKUZlcTckJCIrUkJNNnlGZ24kIituQic9WCpGZXE3JCQiKzcmKVFtdkZnbiQiK2gpUlM0IkZnbjckJCIrbGhvUXRGZ24kIitrbUxzN0ZnbjckJCIrN3ooKkdyRmduJCIrMyN6aFsiRmduNyQkIisySkZRcEZnbiQiKzsxJEh1IkZnbjckJCIrNDkiencnRmduJCIrIipwVl4/RmduNyQkIitsYWs+bUZnbiQiKzlSRUFDRmduNyQkIit0JFxkXCdGZ24kIisiNEl6J0dGZ243JCQiK3BAOypSJ0ZnbiQiK2RwNS5NRmduNyQkIis5InpPTCdGZ24kIisvMXRXU0ZnbjckJCIrMiUzVUknRmduJCIrKCp5JD4iW0ZnbjckJCIrJkc8ckonRmduJCIrYXFvRGRGZ243JCQiKyU+KnAhUSdGZ24kIisvJVt5IW9GZ243JCQiK2g8IWVdJ0ZnbiQiK1ElWyV6ISlGZ243JCQiKykpKmZtcSdGZ24kIitMLWlkJipGZ243JCQiKykqNCg+KyhGZ24kIis2dTFENkYyNyQkIitQUj87dUZnbiQiK246MTo4RjI3JCQiK3NjLCIpekZnbiQiK0xmREE6RjI3JCQiK3lkTk8oKUZnbiQiK05VNlI8RjI3JCQiK1MmcC10KkZnbiQiKyZ5IWVfPkYyNyQkIitHY2EsNkYyJCIrOyNSSTkjRjI3JCQiK3RaJlJFIkYyJCIrT1l5JUcjRjI3JCQiK2VQZ2k5RjIkIitcWWBcQkYyNyQkIisjUVxZcCJGMiQiKyZHPlNKI0YyNyQkIitSOldcPkYyJCIrJikzP3FARjI3JCQiK2ZPQDNBRjIkIitIM3NLPkYyNyQkIitIdyZ5VyNGMiQiKz9xIXBqIkYyNyQkIis6J1F4ayNGMiQiK1JgUEU4RjI3JCQiK2EsZSZ6I0YyJCIrelshei4iRjI3JCQiK1FcNCopR0YyJCIrV0ZSR3pGZ243JCQiK1dVZUxIRjIkIitWYT11ZkZnbjckJCIrPmozUUhGMiQiK2kzUyFbJUZnbjckJCIrI2VZQSJIRjIkIitEQCN5TyRGZ243JCQiKyV5JWZrR0YyJCIrITM8MGIjRmduNyQkIitPMiY+IUdGMiQiKyZcbUomPkZnbjckJCIrXWRYSEZGMiQiKyZbd2heIkZnbjckJCIrRUEmM2wjRjIkIis4VSRcPiJGZ243JCQiK2R5ISlvREYyJCIrc0xTcSYqRmVxNyQkIitaLj8mWyNGMiQiKyh6J0gkeihGZXE3JCQiK1cvTSxDRjIkIitCJ2VNWCdGZXE3JCQiK3BvOD1CRjIkIitpJGZTViZGZXE3JCQiK0xZQE9BRjIkIismKXAhPmwlRmVxNyQkIis2JykqZjojRjIkIitdNFhaU0ZlcTckJCIrIz5zeDIjRjIkIithJip5eE5GZXE3JCQiKzQlPTwrI0YyJCIrR2UhPUAkRmVxNyQkIitdKVt6Iz5GMiQiK2NveUVIRmVxNyQkIisiNENsJj1GMiQiK0QjcGdxI0ZlcTckJCIrQiNvdXkiRjIkIiswS1JQREZlcTckJCIrayN6MnMiRjIkIitPeHU2Q0ZlcTckJCIrem1WYztGMiQiKyk9ZERLI0ZlcTckJCIrYm9TJWYiRjIkIitZIzNeRSNGZXE3JCQiK0xya006RjIkIiszdTlPQUZlcTckJCIrJHozcloiRjIkIitcO2dMQUZlcTckJCIrLiRSPFUiRjIkIitIKipRY0FGZXE3JCQiK0NTW284RjIkIis4T0kvQkZlcTckJCIrSnZHPDhGMiQiKz04JXpQI0ZlcTckJCIrY1k0bzdGMiQiK0IqKW95Q0ZlcTckJCIrTjgmM0EiRjIkIitcJVIoM0VGZXE3JCQiKzdgXXY2RjIkIitXMDpyRkZlcTckJCIrdG4rSzZGMiQiK1gpUipwSEZlcTckJCIrOiozLjQiRjIkIisjcD8tQCRGZXE3JCQiKzomby4wIkYyJCIrI3ojUSlcJEZlcTckJCIrNWw5NzVGMiQiK0lYS1VRRmVxNyQkIitKaTNjKCpGZ24kIiskKmV1XlVGZXE3JCQiKyF5Z3NTKkZnbiQiK2NMXlFaRmVxNyQkIitbVHd1ISpGZ24kIit2MjY8YEZlcTckJCIrdFlXZSgpRmduJCIrW0o+MGdGZXE3JCQiKy5hQmUlKUZnbiQiK2ZKRUNvRmVxNyQkIisoW29UPClGZ24kIisrTFwreUZlcTckJCIraEhSMXpGZ24kIisrKj1kJyopRmVxNyQkIis2cD5id0ZnbiQiK3RHJ2UuIkZnbjckJCIrJ0hONVUoRmduJCIrKCp5aC03RmduNyQkIisjPmxYPyhGZ24kIitIXl8tOUZnbjckJCIrSjJwMXFGZ24kIitgalhVO0ZnbjckJCIrJSozaUdvRmduJCIrdDxvST5GZ243JCQiK2xXJT5uJ0ZnbiQiKyhwLXJGI0ZnbjckJCIrQiFHKFFsRmduJCIrUSopWyRwI0ZnbjckJCIrLl9rSmtGZ24kIitKZXAkPiRGZ243JCQiK0cmW1ROJ0ZnbiQiK3ptIlF6JEZnbjckJCIrKlsqcDVqRmduJCIrdnA+N1hGZ243JCQiKzckeXFJJ0ZnbiQiKzpAPnBgRmduNyQkIisrQCEzTidGZ24kIitbXFsnUSdGZ243JCQiK0YjKW9ea0ZnbiQiKzF4byZlKEZnbjckJCIrWWtnQW1GZ24kIitpQCRlKSopRmduNyQkIitbOlkhKW9GZ24kIisoPkAqZjVGMjckJCIrKj1BdUMoRmduJCIrTUJcVTdGMjckJCIrMGRLX3hGZ24kIitRbiVSVyJGMjckJCIrblwoPlYpRmduJCIrNiQ9JWU7RjI3JCQiK2UvZkokKkZnbiQiK2tfMHY9RjI3JCQiK3FvR101RjIkIitxXXN3P0YyNyQkIisqcClwKj4iRjIkIitPUndSQUYyNyQkIiteTiNcUSJGMiQiKzxGPU9CRjI3JCQiKzkoKlEwO0YyJCIreUJ0UkJGMjckJCIrbHphYD1GMiQiK055Y05BRjI3JCQiKyY0SE02I0YyJCIrQmhiSD9GMjckJCIraVckR08jRjIkIiskM1wudiJGMjckJCIrKjQiUXpERjIkIis2Wk1TOUYyNyQkIitDJVxzdSNGMiQiKyozVC45IkYyNyQkIis2cGNnR0YyJCIrPyFHeXgpRmduNyQkIitKJWZBI0hGMiQiK0FNR1NtRmduNyQkIiteWVRTSEYyJCIrKG9AVClcRmduNyQkIitOIkdaI0hGMiQiKzYoNDJ1JEZnbjckJCIrNSxFJSlHRjIkIispXChwQkdGZ243JCQiK2BBW0VHRjIkIitfMHdfQEZnbjckJCIrUTg4ZEZGMiQiKy5VUGk7RmduNyQkIitEIkgvbyNGMiQiK0VebC04RmduNyQkIitMWlMqZiNGMiQiKztYL1A1RmduNyQkIis9KSk+O0RGMiQiK292JkdSKUZlcTckJCIrQSk0QlYjRjIkIitPJXBwIXBGZXE3JCQiK1UjeihbQkYyJCIrYXo9IXkmRmVxNyQkIitNZEttQUYyJCIrXztCPVxGZXE3JCQiK20tVyY9I0YyJCIrSVp2YFVGZXE3JCQiKy9JWDFARjIkIitHI2UkUVBGZXE3JCQiKy8tZUg/RjIkIitcXi1QTEZlcTckJCIrTXImXCY+RjIkIitcSkRDSUZlcTckJCIrJUhoRSk9RjIkIitpUE8ieSNGZXE3JCQiKyQpKUdGIj1GMiQiK3NVaSVmI0ZlcTckJCIrbmw7WDxGMiQiK2c6JFJYI0ZlcTckJCIrXSpmKno7RjIkIit4aik9TiNGZXE3JCQiK0goenFoIkYyJCIrPjQ5JEcjRmVxNyQkIisraVtjOkYyJCIrMChSUkMjRmVxNyQkIiskZksiKVwiRjIkIitFIip6SkFGZXE3JCQiK3F2Jz5XIkYyJCIrZSQpR1hBRmVxNyQkIissciR6USJGMiQiKzRqKVFHI0ZlcTckJCIrQGYpZkwiRjIkIitodyp5TSNGZXE3JCQiK2wnZWdHIkYyJCIreHBUUUNGZXE3JCQiK3QyNVE3RjIkIitRVUxkREZlcTckJCIrWiNmPz4iRjIkIitGJnl0cSNGZXE3JCQiKy5MKXk5IkYyJCIrOCMqPiMqR0ZlcTckJCIra1xfMDZGMiQiK0ZqWztKRmVxNyQkIis9J1JcMSJGMiQiK1hSOCdRJEZlcTckJCIrJ2YnM0U1RjIkIistRFgzUEZlcTckJCIrJ3koSCopKSpGZ24kIitDIlFDNCVGZXE3JCQiKyFvI1FNJipGZ24kIitUJTQiXFhGZXE3JCQiK3o7KGU+KkZnbiQiKygpWyM+NCZGZXE3JCQiK1FcZXQpKUZnbiQiKzBqSFBkRmVxNyQkIistRlVuJilGZ24kIisyd0AwbEZlcTckJCIrYSd5dEYpRmduJCIreS0uP3VGZXE3JCQiK1ZqYi4hKUZnbiQiK3IsTjYmKUZlcTckJCIrJmYqPVl4RmduJCIrJCk+PTopKkZlcTckJCIrIzNvY10oRmduJCIrJ3lDdjgiRmduNyQkIit5JnBERyhGZ24kIitSYFhDOEZnbjckJCIrNjVxeHFGZ24kIis1MHNbOkZnbjckJCIrdDU6IypvRmduJCIrbnUvPT1GZ243JCQiK1YjZXRzJ0ZnbiQiKz1xcVRARmduNyQkIitjJSo+JmUnRmduJCIrPDl2SURGZ243JCQiKyVvMSJva0ZnbiQiK0RHQykqSEZnbjckJCIrcDlBempGZ24kIis1alZmTkZnbjckJCIrUThnQWpGZ24kIisvaCY9QiVGZ243JCQiKyg+KFwuakZnbiQiKzRhO05dRmduNyQkIit2dnNHakZnbiQiK2dJbiEqZkZnbjckJCIrZE49MmtGZ24kIit0KVEtNyhGZ243JCQiKyQpZV1dbEZnbiQiKyVHSFVXKUZnbjckJCIrbiR5UngnRmduJCIrOidSIXkqKkZnbjckJCIrSHVtKDQoRmduJCIrO2xsczZGMjckJCIrPyd6eGEoRmduJCIrV11mbjhGMjckJCIrakowZSIpRmduJCIrM21AeTpGMjckJCIrKkddMSgqKUZnbiQiKydIJW8meiJGMjckJCIrXFpgLjVGMiQiKy9rSDA/RjI3JCQiKzEmeS85IkYyJCIrKXo2ZD0jRjI3JCQiKz9bRDc4RjIkIitgODE1QkYyNyQkIitpI1ssXyJGMiQiK1EwW11CRjI3JCQiK2xHW2Y8RjIkIisuZ24nRyNGMjckJCIrX2dQPD9GMiQiKyEpcDM7QEYyNyQkIismPShSdEFGMiQiK1BWVGY9RjI3JCQiKyQpZVEvREYyJCIrYTtqYjpGMjckJCIrLiVmOXAjRjIkIitraCl5QyJGMjckJCIrIypbKVwjR0YyJCIrNHdwJHAqRmduNyQkIitIRCpcIUhGMiQiKz86aHJ0RmduNyQkIitnUzdRSEYyJCIrQXBgVmJGZ243JCQiK18rI1EkSEYyJCIrYF1cZFRGZ243JCQiKzNhYSxIRjIkIisnZiUqKUhKRmduNyQkIitlcFNcR0YyJCIrVTNhd0JGZ243JCQiK10lb1B5I0YyJCIrXHcwRT1GZ243JCQiKzJgUDRGRjIkIit4NydIVSJGZ243JCQiK1kkZidIRUYyJCIrIyoqPWg3IkZnbjckJCIrbnQwWkRGMiQiKy5DOGUhKkZlcTckJCIrKFIhR2pDRjIkIisxeFkzdUZlcTckJCIrI28/JnpCRjIkIitaTnFoaEZlcTckJCIraUtnJ0gjRjIkIitpXSU0QCZGZXE3JCQiK11UNDpBRjIkIit0QCkqeldGZXE3JCQiK2VsUE5ARjIkIisrWDw5UkZlcTckJCIrcVZxZD9GMiQiK2whKTN1TUZlcTckJCIrVCdSQSk+RjIkIishXHI1OCRGZXE3JCQiKyoqKnkhND5GMiQiKzQuPWtHRmVxNyQkIitFQkZRPUYyJCIrKilILmVFRmVxNyQkIitvZyQpcDxGMiQiK0UySixERmVxNyQkIismemlQcSJGMiQiK25OcSZRI0ZlcTckJCIrOyJHK2siRjIkIisnKno7MEJGZXE3JCQiK0tkZnk6RjIkIissUlNiQUZlcTckJCIraDZVPjpGMiQiKzkhKVtMQUZlcTckJCIrMlhYaTlGMiQiK2kzaVBBRmVxNyQkIitcRGsyOUYyJCIrISp6JnBFI0ZlcTckJCIrUS4kXE4iRjIkIitFZF1AQkZlcTckJCIrSkNFLzhGMiQiKzUmbz9TI0ZlcTckJCIrelFlYjdGMiQiKyJvUi1eI0ZlcTckJCIrSTUlKTM3RjIkIiskR0olW0VGZXE3JCQiKzxBKVI7IkYyJCIrRjglKj5HRmVxNyQkIitoJGU0NyJGMiQiK3EkZiFISUZlcTckJCIrVk5zejVGMiQiK2VNQCJHJEZlcTckJCIrd2JCUzVGMiQiK2pOOyRlJEZlcTckJCIrKGVjQysiRjIkIisieltLJVJGZXE3JCQiKylwT05tKkZnbiQiK0l6cHJWRmVxNyQkIisiPScqKj0kKkZnbiQiK3ZgLCIpW0ZlcTckJCIrQnl0ISoqKUZnbiQiK3ZyWCdbJkZlcTckJCIrSDRqeScpRmduJCIrIilbaDFpRmVxNyQkIisuSWojUSlGZ24kIitSeTdrcUZlcTckJCIrWFwhRzUpRmduJCIrJlxabDMpRmVxNyQkIit6KUgkUnlGZ24kIispW3p0SSpGZXE3JCQiK2NtYCNmKEZnbiQiK0dJdHc1RmduNyQkIitSJUhITyhGZ24kIis/JHk6RCJGZ243JCQiK0dfQV5yRmduJCIrWUFEaDlGZ243JCQiK0A9U2VwRmduJCIrZDMoSHIiRmduNyQkIitQJ2ZkeSdGZ24kIiskMzNhLCNGZ243JCQiK1I/K05tRmduJCIrdzckKnlCRmduNyQkIitCMU0zbEZnbiQiK3lpJGUiR0ZnbjckJCIrbFZqM2tGZ24kIis7RmJTTEZnbjckJCIrJSlcZFJqRmduJCIrLylvKHBSRmduNyQkIitOYiRmSSdGZ24kIispKT5SQVpGZ243JCQiK1BhIVJKJ0ZnbiQiK3NBQD5jRmduNyQkIitBRmFyakZnbiQiK1FMMSNvJ0ZnbjckJCIrR0NRKlsnRmduJCIrOXI8S3pGZ243JCQiKyJRWDdvJ0ZnbiQiK1ErSShRKkZnbjckJCIrKG8mR2xwRmduJCIrQSoqcDA2RjI3JCQiKyMpeUZsdEZnbiQiKzwqW05IIkYyNyQkIitlTy43ekZnbiQiK04nUiIqXCJGMjckJCIrRlRkVycpRmduJCIrN2JXOjxGMjckJCIrTGo5NScqRmduJCIrUlAyST5GMjckJCIrJTNFaDMiRjIkIitAazhDQEYyNyQkIishbyJvVzdGMiQiKyI0IVxzQUYyNyQkIitdSFRSOUYyJCIrM2Alb00jRjI3JCQiKyZmIj1vO0YyJCIrKVI+SEsjRjI3JCQiK0MoZjcjPkYyJCIrJ0dnMD4jRjI3JCQiK0J5bCE9I0YyJCIrUjElPSc+RjI3JCQiKzw2V0JDRjIkIisrS0hxO0YyNyQkIisjZXQkR0VGMiQiKyI+KVJmOEYyNyQkIis1MTYjeSNGMiQiKy0nW3MxIkYyNyQkIituaE0iKUdGMiQiK0F1eHAiKUZnbjckJCIrJVxUMiRIRjIkIissN1xpaEZnbjckJCIrP0Q1UkhGMiQiK0NTUUFZRmduNyQkIiskZTloIkhGMiQiKyE9KHpzTUZnbjckJCIrVzZacUdGMiQiK3NJVUZFRmduNyQkIissPz00R0YyJCIrJG8xJTQ/RmduNyQkIitiJGV2dCNGMiQiKyFvTHViIkZnbjckJCIraHpaZkVGMiQiK1FGUkQ3RmduNyQkIis5K3J4REYyJCIreDc4KHoqRmVxNyQkIittZT8lXCNGMiQiKyNbS08nekZlcTckJCIrInlHLlQjRjIkIit3VWsjZSdGZXE3JCQiK0NZLUZCRjIkIitLSiRIYCZGZXE3JCQiKztZJVxDI0YyJCIrJVwzI0daRmVxNyQkIitdQGBrQEYyJCIrPi11MVRGZXE3JCQiK09QMyczI0YyJCIrWHkzQ09GZXE3JCQiK048ejQ/RjIkIitKNy9bS0ZlcTckJCIrJVJ1ZCQ+RjIkIitWVjViSEZlcTckJCIrY3E0az1GMiQiK2FxL0dGRmVxNyQkIitBc3kleiJGMiQiK3BYPmFERmVxNyQkIitdYCV5cyJGMiQiK011QUNDRmVxNyQkIitSRURqO0YyJCIreG9MSkJGZXE3JCQiKyRwdzRnIkYyJCIrI1tIMUYjRmVxNyQkIisnb3Y0YSJGMiQiKy4pMyhRQUZlcTckJCIrXDk/JFsiRjIkIit5V1FMQUZlcTckJCIrJnksd1UiRjIkIispM3ZNRCNGZXE3JCQiK0JBN3U4RjIkIitUPm8pSCNGZXE3JCQiKyFRMkZLIkYyJCIrYSY9JnBCRmVxNyQkIitGP0l0N0YyJCIrN2hHbkNGZXE3JCQiK2A+JmVBIkYyJCIrNCMqMyVmI0ZlcTckJCIrKG8vLj0iRjIkIitFJCopR3YjRmVxNyQkIitILGhPNkYyJCIrYmdmWkhGZXE3JCQiK082cyU0IkYyJCIrOSQpPiQ9JEZlcTckJCIrITMlZmE1RjIkIitKKlJmWSRGZXE3JCQiK0YlKj07NUYyJCIrY2dhLlFGZXE3JCQiK2tPcyV6KkZnbiQiK0AiNGI/JUZlcTckJCIrKSpcOFclKkZnbiQiK0lhVyRvJUZlcTckJCIrYWwqKTQiKkZnbiQiK0RoYF5fRmVxNyQkIitBJVs9eilGZ24kIisiM3BxI2ZGZXE3JCQiK1JTIioqWylGZ24kIisrODVKbkZlcTckJCIrcVU2LyMpRmduJCIrP1ZFKm8oRmVxNyQkIisoeSVlTXpGZ24kIis8eHRLKSlGZXE3JCQiKzt4ZiJvKEZnbiQiK0k8JSo+NUZnbjckJCIrI3opZVh1RmduJCIrMC9gJD0iRmduNyQkIiswPD5Gc0ZnbiQiK0dOaHo4RmduNyQkIisjKT1HRnFGZ24kIitQbyNcaCJGZ243JCQiK01GLlpvRmduJCIrJ1whZSgqPUZnbjckJCIrSCYpKXpvJ0ZnbiQiKyhHIUhQQUZnbjckJCIrLSdmQGInRmduJCIrIlE6Y2sjRmduNyQkIisueTlVa0ZnbiQiK0Qiemg4JEZnbjckJCIrUkpKaGpGZ24kIitENiRbcyRGZ243JCQiK1tuKlJKJ0ZnbiQiK04yb0hXRmduNyQkIitUNCNlSSdGZ24kIit4IyopM0YmRmduNyQkIit0VDRXakZnbiQiKzFiL3FpRmduNyQkIitQcE5Ra0ZnbiQiK0YrJylbdUZnbjckJCIrXzo2LG1GZ24kIitvbXpFKSlGZ243JCQiKyNleShbb0ZnbiQiKycpPnJUNUYyNyQkIitnLSFIPyhGZ24kIitSODJBN0YyNyQkIis/VmAicChGZ24kIisiKT5xQDlGMjckJCIraWdlXSQpRmduJCIrXSd6XmoiRjI3JCQiKy5hU0MjKkZnbiQiK1tMRF89RjI3JCQiKz4wVU81RjIkIiskKTRdYz9GMjckJCIremc8Iz0iRjIkIitqJDRcQSNGMjckJCIrWHJdajhGMiQiKyJlVShIQkYyNyQkIitWIT0vZSJGMiQiKycpKmZUTSNGMjckJCIrIiplQEU9RjIkIitRb2deQUYyNyQkIis+SSFlMyNGMiQiKz1bWWI/RjI3JCQiK1ZNVFBCRjIkIitLZz4jeSJGMjckJCIrKGVVJGVERjIkIitAVVV0OUYyNyQkIitGPyQ9dCNGMiQiK2hQITM8IkYyNyQkIisjek40JkdGMiQiK1NwJ1suKkZnbjckJCIrPV96PEhGMiQiKzlkRVdvRmduNyQkIitoQThTSEYyJCIrV2RkUl5GZ243JCQiK3A5Z0ZIRjIkIiswJj5qJlFGZ243JCQiK1xMUiopR0YyJCIrXHpmM0hGZ243JCQiKylba0okR0YyJCIrdnIlW0AjRmduNyQkIitMKkdbdyNGMiQiKzh5JXlxIkZnbjckJCIrTlh2KW8jRjIkIitfI1JoTCJGZ243JCQiKygpPjMzRUYyJCIrdy4qPTEiRmduNyQkIisyUC5EREYyJCIrSyN5KXkmKUZlcTckJCIrT247VENGMiQiK2BbY1pxRmVxNyQkIis0VGNkQkYyJCIrIUdXdSllRmVxNyQkIiswUShcRiNGMiQiKyIzazIrJkZlcTckJCIrJVwyUj4jRjIkIitIenQ8VkZlcTckJCIrIyoqNFo2I0YyJCIrKXlZIyl5JEZlcTckJCIrZHpnUD9GMiQiKy9bMHdMRmVxNyQkIit2U3VpPkYyJCIrUSQqeWFJRmVxNyQkIitbNj8hKj1GMiQiKypmV14hR0ZlcTckJCIrcyI+KyM9RjIkIitZMCRIaCNGZXE3JCQiK0p2P188RjIkIituaHFuQ0ZlcTckJCIrbFB2J28iRjIkIisuJj09TyNGZXE3JCQiK08qSE9pIkYyJCIrRmFzKkcjRmVxNyQkIispRyh6aTpGMiQiK3cpPXZDI0ZlcTckJCIrIXA0VV0iRjIkIitUO2ZLQUZlcTckJCIrK2kieVciRjIkIitbKTNNQyNGZXE3JCQiKzdJYyRSIkYyJCIrXHhOekFGZXE3JCQiK1lcUlQ4RjIkIit4IWUxTSNGZXE3JCQiKzZtRCJIIkYyJCIrTTlLR0NGZXE3JCQiK0tMNFY3RjIkIitbRTpXREZlcTckJCIrQz4mbz4iRjIkIisoXCp5IXAjRmVxNyQkIitaOFtfNkYyJCIrJCopejwoR0ZlcTckJCIrM0wkKjQ2RjIkIisqKnpwIjQkRmVxNyQkIis/RztwNUYyJCIrWUdJY0xGZXE3JCQiK3koRywuIkYyJCIrcV51c09GZXE3JCQiK3NiJXojKipGZ24kIis/LiQpXFNGZXE3JCQiK15vR3ImKkZnbiQiKzt1TSlcJUZlcTckJCIrJCpmMEojKkZnbiQiK0d2WkpdRmVxNyQkIit2YTEyKilGZ24kIitQakhsY0ZlcTckJCIrO20/KmYpRmduJCIrYigqUT5rRmVxNyQkIitcRVkyJClGZ24kIitgUmc8dEZlcTckJCIrKnlDPi4pRmduJCIrJSpSJiopUSlGZXE3JCQiKyU9N0d4KEZnbiQiKzh0cm8nKkZlcTckJCIrcW9cSXZGZ24kIitnUSgqPjZGZ243JCQiK19fYDB0RmduJCIrJylwUi44RmduNyQkIis5IjMoKTQoRmduJCIrOmhVQjpGZ243JCQiK2s9MjZwRmduJCIrdCpSd3kiRmduNyQkIitKWi1XbkZnbiQiK2FpODBARmduNyQkIit0SFIqZidGZ24kIitxKm9uWyNGZ243JCQiKzpUYXprRmduJCIrayJ6YCVIRmduNyQkIitbcGAoUSdGZ24kIisyXClmXCRGZ243JCQiKyw8S0ZqRmduJCIrdkMoZTolRmduNyQkIitBKjNTSSdGZ24kIit1UFxXXEZnbjckJCIrQEBCQ2pGZ24kIitnLi0kKWVGZ243JCQiK15yaidSJ0ZnbiQiKzNkSyQqcEZnbjckJCIrOyJSRGAnRmduJCIrLy4xJ0gpRmduNyQkIitNNnpZbkZnbiQiKyV6XXQhKSpGZ243JCQiK1V3KillcUZnbiQiKyh5Xkw6IkYyNyQkIittTU4lXChGZ24kIismKVJKWThGMjckJCIrVnovJzMpRmduJCIrbG1mYjpGMjckJCIrOS1EdikpRmduJCIrP3UqR3giRjI3JCQiKzI2MTYqKkZnbiQiKyYzJj0lKT5GMjckJCIrZThpQzZGMiQiK2kkNClvQEYyNyQkIiskSDFFSCJGMiQiK2dITStCRjI3JCQiKyQqW3onXCJGMiQiKyo9JnBdQkYyNyQkIitKeEVMPEYyJCIrTCcpSClII0YyNyQkIis6LS8hKj5GMiQiK1AmUiRRQEYyNyQkIisjZj10QyNGMiQiKyg9NSIqKT1GMjckJCIrWTAiPlsjRjIkIitqTUEpZSJGMjckJCIrdCEqPXVFRjIkIis5bzd6N0YyNyQkIiskcGFNIkdGMiQiK3RNLmwqKkZnbjckJCIrTy4kKSkqR0YyJCIrMF1GImYoRmduNyQkIispM1xrJEhGMiQiKytRNThkRmduNyQkIituNGhOSEYyJCIreEphJUclRmduNyQkIisjNGplIUhGMiQiK1JbX0JLRmduNyQkIitoU1tiR0YyJCIrYFsxWENGZ243JCQiK1oqPTV6I0YyJCIrbz4+dz1GZ243JCQiKzsjcHRyI0YyJCIrQkp6ZjlGZ243JCQiKysnKTNRRUYyJCIraHBPYDZGZ243JCQiK2RWcWJERjIkIis2OFhoIypGZXE3JCQiK3ZGKj5aI0YyJCIrc2FnaHZGZXE3JCQiKy9EPilRI0YyJCIrL3k4eWlGZXE3JCQiK3U4OzBCRjIkIishPm8tSSZGZXE3JCQiK08mKVtCQUYyJCIrLVYwXFhGZXE3JCQiK2hOZFZARjIkIis/RCR6J1JGZXE3JCQiK04+b2w/RjIkIitAIjRoXiRGZXE3JCQiK3FSKSoqKT5GMiQiK3BaJ1I7JEZlcTckJCIrJT4jZTs+RjIkIitQKWYpKilHRmVxNyQkIis3MmBYPUYyJCIrSCMzem4jRmVxNyQkIislKipbb3giRjIkIiskRz9rXiNGZXE3JCQiK3A8YDU8RjIkIitPJDRvUiNGZXE3JCQiK3loYlk7RjIkIisnZUBHSiNGZXE3JCQiKyQqcClbZSJGMiQiKztNKipmQUZlcTckJCIrZS9bRDpGMiQiK0oyRk5BRmVxNyQkIitdckdvOUYyJCIrdWJ1T0FGZXE3JCQiK0RURDg5RjIkIisiZXpNRSNGZXE3JCQiK2tsS2c4RjIkIisncCRSOkJGZXE3JCQiK08hXCU0OEYyJCIrbDBAJFIjRmVxNyQkIis1bGNnN0YyJCIrYzVXKVwjRmVxNyQkIit2Xmk4N0YyJCIrPCxUTEVGZXE3JCQiK05KZG82RjIkIitkLEssR0ZlcTckJCIrQTVPRDZGMiQiK15ETjFJRmVxNyQkIit1RCVSMyJGMiQiKzMsIlFEJEZlcTckJCIrJT52VS8iRjIkIiszYUldTkZlcTckJCIrUzBLMTVGMiQiK2c/Ky9SRmVxNyQkIitnO1grKCpGZ24kIis0Iz1cSyVGZXE3JCQiK2NDQGEkKkZnbiQiK11TSURbRmVxNyQkIitQVEZDISpGZ24kIitBeTU/YUZlcTckJCIrN3hcNSgpRmduJCIrQXdhRmhGZXE3JCQiKyw1JEdUKUZnbiQiKzElNClwcEZlcTckJCIrNU5LSiIpRmduJCIrSlcqUSh6RmVxNyQkIitWWzlteUZnbiQiK1I8a3MiKkZlcTckJCIrandnPHdGZ24kIitTX2ZnNUZnbjckJCIrJGUnPidRKEZnbiQiKyYzQ0FCIkZnbjckJCIrbVpnc3JGZ24kIithSyxROUZnbjckJCIrIm8heXhwRmduJCIreDgvJm8iRmduNyQkIisuJCkpSCFvRmduJCIreDcjPSk+RmduNyQkIishRyYpKVxtRmduJCIrTFdgUUJGZ243JCQiKydHRDFfJ0ZnbiQiKyU0bnN3I0ZnbjckJCIrJT0kKnpUJ0ZnbiQiK1IvQSNHJEZnbjckJCIrT1xlWGpGZ24kIitHbiUpKipRRmduNyQkIitJKFshM2pGZ24kIitHPSQpUVlGZ243JCQiKzcpNDlKJ0ZnbiQiKyUqcHo+YkZnbjckJCIrJSkqM05PJ0ZnbiQiKyNlQlhjJ0ZnbjckJCIrb0dmdWtGZ24kIisoKnpTJXooRmduNyQkIitPKio0ZW1GZ24kIitoJEh4QSpGZ243JCQiK2xAb0pwRmduJCIrNXNeKDMiRjI3JCQiKzcmcCU9dEZnbiQiK0dSSHQ3RjI3JCQiKzZRW1t5RmduJCIrIypcR3g5RjI3JCQiKyFcbClmJilGZ24kIitXPSRIcCJGMjckJCIrKFx0ISpcKkZnbiQiK0wmZSUzPkYyNyQkIispZk09MiJGMiQiK3oxbjBARjI3JCQiK1ojZm5BIkYyJCIraDMqKmZBRjI3JCQiK0dldTw5RjIkIitROD9WQkYyNyQkIis9V0hWO0YyJCIrajw8SUJGMjckJCIrNFRhJSo9RjIkIitjZSQpM0FGMjckJCIrci5HYUBGMiQiK2gkPikpKT5GMjckJCIrZSQ0KSpSI0YyJCIrcW0iPXEiRjI3JCQiK09cUjRFRjIkIis8PysiUiJGMjckJCIrQTRxb0ZGMiQiKyJbNGM0IkYyNyQkIitUIUdNKEdGMiQiKyFHOFpTKUZnbjckJCIrKlImZUZIRjIkIitseGlZakZnbjckJCIrI3pAKFJIRjIkIisrIUc7dyVGZ243JCQiK24tYT5IRjIkIitELyJmZCRGZ243JCQiK0NGKWUoR0YyJCIrKGU4SXEjRmduNyQkIitYJ1JmIkdGMiQiK1tXb2s/RmduNyQkIitYbj1YRkYyJCIrQnInemYiRmduNyQkIis/U2puRUYyJCIrcGVIYjdGZ243JCQiK0IuOidlI0YyJCIrPU4mPisiRmduNyQkIis7JmZGXSNGMiQiK3l1Y0kiKUZlcTckJCIrSG0oKT1DRjIkIis9cTg0bkZlcTckJCIrcVRbTkJGMiQiKzVlbUhjRmVxNyQkIis/I2ZLRCNGMiQiK2tiIUchW0ZlcTckJCIrTk1tc0BGMiQiKzZibmtURmVxNyQkIitpaSslNCNGMiQiKz13SnBPRmVxNyQkIitnIypbPD9GMiQiKydSV01HJEZlcTckJCIrUHNCVj5GMiQiK2AoKXkjKUhGZXE3JCQiK0gtS3I9RjIkIisoPmomXEZGZXE3JCQiKy0icDwhPUYyJCIrY19vcURGZXE3JCQiKz1uZU08RjIkIitySWBPQ0ZlcTckJCIrJylmdnA7RjIkIisvKHArTSNGZXE3JCQiKzdkQzI7RjIkIitGZkF3QUZlcTckJCIrLFwsWjpGMiQiKzBpWVRBRmVxNyQkIitLZiwqWyJGMiQiK1RrW0xBRmVxNyQkIitWcD5MOUYyJCIrKlEzNUQjRmVxNyQkIitXT116OEYyJCIrXXBrJEgjRmVxNyQkIistMil5SyJGMiQiK2FXJD1PI0ZlcTckJCIrU0dGeTdGMiQiK0A8emNDRmVxNyQkIis6ZGlJN0YyJCIrXSlSMGUjRmVxNyQkIitjbSlbPSJGMiQiKyQpKVxmdCNGZXE3JCQiKy5gK1Q2RjIkIis9LCRvI0hGZXE3JCQiKyM+TSopNCJGMiQiKyU0YCFlSkZlcTckJCIrPSRIJ2U1RjIkIis7a3NOTUZlcTckJCIrNDIwPzVGMiQiK19PVG5QRmVxNyQkIisrMGpKKSpGZ24kIispejhDOyVGZXE3JCQiK0tPT3olKkZnbiQiK0ssNktZRmVxNyQkIistIm9NOSpGZ24kIis7PlMhPiZGZXE3JCQiKyhHd1AjKSlGZ24kIitGX0JhZUZlcTckJCIrVEU/PyYpRmduJCIrXUBEV21GZXE3JCQiKyFcZEZCKUZnbiQiK3IwZSZlKEZlcTckJCIrVVNjaHpGZ24kIis7JSp5MygpRmVxNyQkIitdJnpvcShGZ24kIislZS5eKyJGZ243JCQiK1U9N3B1RmduJCIraEl1bDZGZ243JCQiK2BHISpbc0ZnbiQiKyNRayNlOEZnbjckJCIrSDgyWnFGZ24kIiskSHcjKmUiRmduNyQkIitfd3drb0ZnbiQiK3kkUm4nPUZnbjckJCIrLVtcLm5GZ24kIis0ZD4rQUZnbjckJCIrIyo0QGxsRmduJCIrR24rLEVGZ243JCQiK0A7WF9rRmduJCIrKj14RDMkRmduNyQkIitWNlxvakZnbiQiK3FjX2dPRmduNyQkIisqXGZ2SidGZ24kIisjPk5GTiVGZ243JCQiKz1LODBqRmduJCIrM2s8el5GZ243JCQiK2gjRyRRakZnbiQiK0FgTGhoRmduNyQkIitcKEhrVSdGZ24kIitjOio0SyhGZ243JCQiKz81ZyJlJ0ZnbiQiK24keXpuKUZnbjckJCIrKVE6KT5vRmduJCIrInpUWS0iRjI3JCQiK1tcLmlyRmduJCIrIXB6Rz8iRjI3JCQiK1c6Zk53RmduJCIrVEFzKzlGMjckJCIrV2JhdiMpRmduJCIrKT5cSmgiRjI3JCQiKyR5MGE3KkZnbiQiKzBgWUk9RjI3JCQiKypHJ2VCNUYyJCIrJVs8cC4jRjI3JCQiKzA3ImY7IkYyJCIrLSZHLEAjRjI3JCQiK3luYVY4RjIkIit5Om1BQkYyNyQkIitoQixkOkYyJCIrQGZAWkJGMjckJCIrLCEzLyE9RjIkIitaOXNsQUYyNyQkIitUNVtmP0YyJCIrSCFRI3o/RjI3JCQiKyhHUkhKI0YyJCIrLD0xNz1GMjckJCIreCVceWAjRjIkIitwQCJcXSJGMjckJCIrLFNnO0ZGMiQiK0klPSw/IkYyNyQkIiteIkc3JUdGMiQiKz1uOyVHKkZnbjckJCIreCl5SSJIRjIkIitZWDxWcUZnbjckJCIrNTdcUkhGMiQiKydRKG8iSCZGZ243JCQiK0phMElIRjIkIithI3AncFJGZ243JCQiK05bMiUqR0YyJCIreW4iPipIRmduNyQkIitIKHokUkdGMiQiK2MlKnl2QUZnbjckJCIrXVowc0ZGMiQiKzgtWl88RmduNyQkIitEKjRtcCNGMiQiKzRtKCpvOEZnbjckJCIrIj0mSDtFRjIkIisxcUInMyJGZ243JCQiK15LVExERjIkIitPMytoKClGZXE3JCQiK1IpeSZcQ0YyJCIrRzk0Jj0oRmVxNyQkIismKWUiZk8jRjIkIitUN0YjKmZGZXE3JCQiK2M3PyRHI0YyJCIrcHlPIjMmRmVxNyQkIis3cSc+PyNGMiQiK1EtPiFRJUZlcTckJCIrJ1JzRDcjRjIkIitkXSNwJFFGZXE3JCQiKzNWRFg/RjIkIitrbTg5TUZlcTckJCIrSUY7cT5GMiQiK1hrZiUzJEZlcTckJCIrJG8mUSgqPUYyJCIrU0FRR0dGZXE3JCQiKytwJ3AjPUYyJCIrQl8mM2ojRmVxNyQkIitPJT0qZTxGMiQiKyVIVjdbI0ZlcTckJCIrVShIS3AiRjIkIitUYGtyQkZlcTckJCIrSFUoKUg7RjIkIisqZkhqSCNGZXE3JCQiK3dTIilvOkYyJCIrRiZSN0QjRmVxNyQkIiskeS4rXiJGMiQiK3QhW09CI0ZlcTckJCIrNkdSYDlGMiQiK0UqPj5DI0ZlcTckJCIrI2VGKilSIkYyJCIrMVZOdkFGZXE3JCQiKzNJYlk4RjIkIis3KSozTUJGZXE3JCQiKyNvOGlIIkYyJCIrVSNmIT5DRmVxNyQkIitLWyZ5QyJGMiQiK1twKT5gI0ZlcTckJCIrKjRCOT8iRjIkIisieUNhbiNGZXE3JCQiKzNzJ286IkYyJCIrLmIjRyZHRmVxNyQkIit0JlFUNiJGMiQiKyIpW2xvSUZlcTckJCIrZj0+dDVGMiQiK11lYUdMRmVxNyQkIitdYilSLiJGMiQiK0ItXVJPRmVxNyQkIitBYiNbJyoqRmduJCIrJ2ZYLCwlRmVxNyQkIis0JjNsZypGZ24kIislUmQ1WCVGZXE3JCQiK0dFa2sjKkZnbiQiK0ImZV4oXEZlcTckJCIreUouUiopRmduJCIrXj1AKWYmRmVxNyQkIitTSmNIJylGZ24kIitlT1VSakZlcTckJCIrIipmP08kKUZnbiQiK0gjekBBKEZlcTckJCIrWjkvZiEpRmduJCIrN1gkXEYpRmVxNyQkIitzXEYpeihGZ24kIitXKipHSyYqRmVxNyQkIisnPWlVYihGZ24kIisxdGkuNkZnbjckJCIrPCdSdksoRmduJCIrOmF5JEciRmduNyQkIithUCcpPXJGZ24kIitgNigpKlwiRmduNyQkIitPNUVIcEZnbiQiK2ZZS2Y8RmduNyQkIitkQjRnbkZnbiQiK2Q8M3I/RmduNyQkIis9cTg4bUZnbiQiKypvNWVXI0ZnbjckJCIrJypHcSFcJ0ZnbiQiK0NZOScqR0ZnbjckJCIrY0F4JlInRmduJCIrQXgob1YkRmduNyQkIisxaT5LakZnbiQiKzUmcF0zJUZnbjckJCIrTmQmXEknRmduJCIrU3EnKmZbRmduNyQkIislcDgwSydGZ24kIitVOmcjeSZGZ243JCQiK0kpKkcoUSdGZ24kIitELiVbKG9GZ243JCQiKyQqKilIO2xGZ24kIitjL2RkIilGZ243JCQiK3RvKj5zJ0ZnbiQiKyZlZnZrKkZnbjckJCIrMCRvTC0oRmduJCIrKyxDTjZGMjckJCIrXnRDWHVGZ24kIisiKnBHRThGMjckJCIrXihHKD4hKUZnbiQiKzZmQE06RjI3JCQiK20vQih5KUZnbiQiKytrQF48RjI3JCQiK0lEPCd6KkZnbiQiK3lNKVEnPkYyNyQkIis6biMqNDZGMiQiKzglSEE6I0YyNyQkIitQUkx1N0YyJCIrZ29HIUgjRjI3JCQiK1okcFxaIkYyJCIrKDROKVxCRjI3JCQiK0JaZzM8RjIkIit6cj8zQkYyNyQkIitxJDNUJz5GMiQiKytHX2VARjI3JCQiK3JvTEFBRjIkIisxNHY7PkYyNyQkIiszQzlnQ0YyJCIrJz1vIT47RjI3JCQiK00jXHNsI0YyJCIrPiwwNDhGMjckJCIrLVYmPiFHRjIkIisqKjNzQTVGMjckJCIrRygpWyMqR0YyJCIrJylmI1sheUZnbjckJCIrOkFZTUhGMiQiK0QqPSd5ZUZnbjckJCIrRTgsUEhGMiQiKyNHZSkzV0ZnbjckJCIrP2d1NEhGMiQiKy1wRDpMRmduNyQkIitqTTVoR0YyJCIrNypHQV4jRmduNyQkIis1PiF5eiNGMiQiK1BcS0Q+RmduNyQkIisoKlIqW3MjRjIkIitQOChlXCJGZ243JCQiK0xDMFlFRjIkIitHIlErPSJGZ243JCQiKztKKlFjI0YyJCIrUyMqSGclKkZlcTckJCIrQ2tEIVsjRjIkIis7WUM2eEZlcTckJCIrRm1VJ1IjRjIkIisxYyI9UidGZXE3JCQiKzFRSDhCRjIkIitBLVQoUSZGZXE3JCQiK3MpcDlCI0YyJCIrLXpSO1lGZXE3JCQiKyRbcTg6I0YyJCIrIz1CLi0lRmVxNyQkIitHREZ0P0YyJCIrL2IwZE5GZXE3JCQiK1JXTigqPkYyJCIrKmVDZz4kRmVxNyQkIitQWHNCPkYyJCIrOyYzXCJIRmVxNyQkIit3NVdfPUYyJCIrQWxLKHAjRmVxNyQkIisrbF8keSJGMiQiK2RlQUpERmVxLSUlRk9OVEc2JCUqSEVMVkVUSUNBRyIjNQ==

solve( {x*(c1*y-c2) = 0, y*(c3-c4*x) = 0}, {x,y} );

NiQ8JC9JInhHNiIkIiIhRigvSSJ5R0YmRic8JC9GKiQiJCsmRigvRiUkIiUrOkYo

solve( {-y-z=0, x+c1*y=0, c2+z*(x-c3)=0}, {x,y,z});

NiQ8JS9JInlHNiIkISs8QzU4TiEjNi9JInhHRiYkIitNWz9FcSEjNy9JInpHRiYkIis8QzU4TkYpPCUvRiUkISspKm9bWUchIikvRjAkIispKm9bWUdGNy9GKyQiKyZ6dEhwJiEiKg==

au := runge4(f1, f2,a,b,h,1500, 500);

NiM+SSNhdUc2IjckN2VpbiIlKzokRigiIiFGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRik3ZWluIiQrJiRGLEYqRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYtRi1GLUYt

plist := []: for i from 1 to n do tmp := [op(i,au[1]), op(i,au[2])]; plist := [op(plist), tmp]; end do:listplot(plist);

LSUlUExPVEc2JC0lJ0NVUlZFU0c2IzdkaW43JCQiJSs6IiIhJCIkKyZGLEYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKUYpRilGKS0lJUZPTlRHNiQlKkhFTFZFVElDQUciIzU=

restart;with(plots,listplot):with(plots,display):runge4 := proc(f1, f2, a, b, h, alpha1, alpha2) local i, n, au1, au2, k11, k12, k21, k22, k31, k32, k41, k42, tmp1, tmp2; n := (b-a)/h; au1 := [alpha1]; au2 := [alpha2]; for i from 0 to n-1 do k11 := evalf( h*f1(a+i*h, op(-1,au1), op(-1,au2))); k12 := evalf( h*f2(a+i*h, op(-1,au1), op(-1,au2))); k21 := evalf( h*f1( a+(i+1/2)*h,op(-1,au1)+k11/2, op(-1,au2)+k12/2)); k22 := evalf( h*f2( a+(i+1/2)*h,op(-1,au1)+k11/2, op(-1,au2)+k12/2)); k31 := evalf( h*f1( a+(i+1/2)*h,op(-1,au1)+k21/2, op(-1,au2)+k22/2)); k32 := evalf( h*f2( a+(i+1/2)*h,op(-1,au1)+k21/2, op(-1,au2)+k22/2)); k41 := evalf( h*f1( a+(i+1)*h,op(-1,au1)+k31, op(-1,au2)+k32)); k42 := evalf( h*f2( a+(i+1)*h,op(-1,au1)+k31, op(-1,au2)+k32)); tmp1:=op(-1,au1)+(1/6)*(k11+2*k21+2*k31+k41); tmp2:=op(-1,au2)+(1/6)*(k12+2*k22+2*k32+k42); au1:=[op(au1),tmp1]; au2:=[op(au2),tmp2]; end do; return [au1,au2]; end proc:eul_mod := proc(f1, f2, a, b, h, alpha1, alpha2) local i, n, au1, au2, k11, k12, k21, k22, k31, k32, k41, k42, tmp1, tmp2; n := (b-a)/h; au1 := [alpha1]; au2 := [alpha2]; for i from 0 to n-1 do #tmp1 := op(-1,au1)+h/2*(f1(t(i),op(-1,au1),op(-1,au2))+f1(t(i+1),op(-1,au1)+h*f1(t(i),op(-1,au1),op(-1,au2)),op(-1,au2)+h*f1(t(i),op(-1,au1),op(-1,au2)))); #tmp2 := op(-1,au1)+h/2*(f2(t(i),op(-1,au1),op(-1,au2))+f2(t(i+1),op(-1,au1)+h*f2(t(i),op(-1,au1),op(-1,au2)),op(-1,au2)+h*f2(t(i),op(-1,au1),op(-1,au2)))); tmp1 := op(-1,au1)+(h/2)*(f1(t(i),op(-1,au1),op(-1,au2))+f1(t(i+1),op(-1,au1)+h*f1(t(i),op(-1,au1),op(-1,au2)),op(-1,au1)+h*f1(t(i),op(-1,au1),op(-1,au2)))); tmp2 := op(-1,au1)+(h/2)*(f2(t(i),op(-1,au1),op(-1,au2))+f1(t(i+1),op(-1,au1)+h*f2(t(i),op(-1,au1),op(-1,au2)),op(-1,au1)+h*f2(t(i),op(-1,au1),op(-1,au2)))); au1:=[op(au1),evalf(tmp1)]; au2:=[op(au2),evalf(tmp2)]; end do; return [au1,au2]; end proc:################################################################ # y' = u2; # y = u1; # # Gleichung umformen nach u3 und ersetzten ################################################################f1 := (t,u1,u2) -> u2; f2 := (t,u1,u2) -> -sin(u1)-0.2*u2;

NiM+SSNmMUc2ImYqNiVJInRHRiVJI3UxR0YlSSN1MkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlOSZGJUYlRiU=

NiM+SSNmMkc2ImYqNiVJInRHRiVJI3UxR0YlSSN1MkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYtSSRzaW5HNiRJKnByb3RlY3RlZEdGMkkoX3N5c2xpYkdGJTYjOSUhIiIqJiQiIiNGNiIiIjkmRjpGNkYlRiVGJQ==

a:=0;

NiM+SSJhRzYiIiIh

b:=50;

NiM+SSJiRzYiIiNd

n:=50;

NiM+SSJuRzYiIiNd

h:=(b-a)/n;

NiM+SSJoRzYiIiIi

n:=(b-a)/h;

NiM+SSJuRzYiIiNd

au := eul_mod(f1, f2,a,b,h,2,0);

NiM+SSNhdUc2IjckN1UiIiMkIiIkIiIhJCIrdERxIWYnISIqJCIrTSR6RVEiISIpJCIrbi1hYEhGMSQiKyY+LktCJ0YxJCIrSXlaXjghIigkIitTL3AkKUdGOCQiKzMuMCM9J0Y4JCIrVj5NQjghIickIitKIXlDJEdGPyQiK3hbIVwxJ0Y/JCIrNiN5JClIIiEiJiQiK0xaWnpGRkYkIitGck5dZkZGJCIrJW9nUUYiISIlJCIrW00yRkZGTSQiK2N4OFFlRk0kIis5QCQpXDchIiQkIitndGt2RUZUJCIrRydSIUdkRlQkIit5O0VFNyEiIyQiKzB1PURFRmVuJCIrLFksP2NGZW4kIitLYzguNyEiIiQiK1oqeWNkI0ZcbyQiK2ZuLTliRlxvJCIrIXBYLz0iRiskIis5VTVGREYrJCIrQHguNWFGKyQiK2lPPWU2IiIiJCIrXmJXekNGaW8kIisyKTQhM2BGaW8kIit1OU1PNkYoJCIrIm8mb0tDRigkIitSZyF6PyZGKCQiKzM3Ilw2IkYqJCIrZXchb1EjRiokIitJLHA0XkYqJCIrJzQmKVE0IiIiJSQiK1tbelRCRlxxJCIrdmtLOF1GXHEkIis5YkR0NSIiJiQiK000aihII0ZjcSQiK1IseT1cRmNxJCIrJylcLGA1IiInJCIrMioqSGFBRmpxJCIrYm8sRVtGanEkIit1aDpMNSIiKCQiK2dneTZBRmFyJCIrIipIK05aRmFyN1VGKyQiK3REcSE0I0YuJCIreCVSMiVSRi4kIitkRUAmeilGLiQiKyZ6IypHIT1GMSQiKzJOKFw7JUYxJCIrWnB0ayYpRjEkIitbWV5jPUY4JCIrdFJNZ1JGOCQiK2U2bHUlKUY4JCIrSXk9Oz1GPyQiKyJ6Q2spUUY/JCIrdFQyPiQpRj8kIitGXTkieSJGRiQiK2A2MjhRRkYkIitCVSNHOylGRiQiKyVlRnZ1IkZNJCIrK1g2VFBGTSQiKyopPSoqMyEpRk0kIispZW9YciJGVCQiK1Z0YnFPRlQkIit3KVt6JnlGVCQiKyVcTEFvIkZlbiQiK0FXTCxPRmVuJCIrJCpcdjR4RmVuJCIrUiQzMGwiRlxvJCIrZ25UTE5GXG8kIisoeWNWYyhGXG8kIisrOVE+O0YrJCIrUyt5bU1GKyQiK3YwcUB1RiskIisiW1QpKWUiRmlvJCIrSStTLE1GaW8kIiswWnQiRyhGaW8kIis8digpZTpGKCQiKz1JRFBMRigkIitiJTNXOShGKCQiK1YneSVIOkYqJCIrZWRKdUtGKiQiK1Y/bjRxRiokIisuVWorOkZccSQiK0phYzdLRlxxJCIrRm1aeG9GXHEkIitRUExzOUZjcSQiK2Anej46JEZjcSQiK3hVeFpuRmNxJCIrJSpwY1c5RmpxJCIraWtgIzQkRmpxJCIrIip6Xj9tRmpxJCIrLFJLPDlGYXIkIis1VkBNSUZhcg==

i:=(33.7-a)/h+1;

NiM+SSJpRzYiJCIkWiQhIiI=

op(34,au[1]);

NiMkIit1OU1PNiIiIw==

op(338,au[1]);

Error, improper op or subscript selector

op(338,au[2]);

NiMkIStcc1RNWiEjNg==