Pienimm"an neli"osumman sovitus V2 19.3. 2002, V3 lokakuu -02 2013kevat/maple/

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkjbWlHRiQ2JVEocmVzdGFydEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJw== with(LinearAlgebra):alias(Tr=Transpose): LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEld2l0aEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JKG1mZW5jZWRHRiQ2JC1GIzYjLUYsNiVRJnBsb3RzRidGL0YyL0YzUSdub3JtYWxGJy1JI21vR0YkNi1RIjpGJ0Y9LyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0ZFLyUpc3RyZXRjaHlHRkUvJSpzeW1tZXRyaWNHRkUvJShsYXJnZW9wR0ZFLyUubW92YWJsZWxpbWl0c0dGRS8lJ2FjY2VudEdGRS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRlQ= xd:=[-1.3,-0.1,0.2,1.3]; yd:=[0.103,1.099,0.808,1.897]; n:=nops(xd); plot([seq([xd[i],yd[i]],i=1..n)],style=point,symbol=circle,symbolsize=15); datakuva:=%: C:=<<1,1,1,1>|<op(xd)>>; M:=Tr(C).C; B:=Tr(C).Vector(yd); ab:=LinearSolve(M,B); suora:=ab[1]+ab[2]*x; display(plot(suora,x=-2..2,color=blue),datakuva);

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYrLUkjbWlHRiQ2JVEpJXJlc3RhcnRGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSI7RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0YxLyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLUY2Ni1RIn5GJ0Y5RjsvRj9GPUZARkJGREZGRkhGSi9GTkZMLUYsNiVRJXdpdGhGJ0YvRjItSShtZmVuY2VkR0YkNiQtRiM2Ji1GLDYlUSZwbG90c0YnRi9GMi8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJStleGVjdXRhYmxlR0Y9RjlGOS1GNjYtUSI6RidGOUY7RlNGQEZCRkRGRkZIL0ZLRk9GTUZqbkZdb0Y5

Muistellaan aluksi interpolaatioteht. Esim:

xd:=[$1..10];yd:=[1.3,3.5,4.2,5.0,7.0,8.8,10.1,12.5,13.0,15.6]; LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVElcGxvdEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JKG1mZW5jZWRHRiQ2JC1GIzYlLUYsNiVRI3hkRidGL0YyLUkjbW9HRiQ2LVEiLEYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGMS8lKXN0cmV0Y2h5R0ZFLyUqc3ltbWV0cmljR0ZFLyUobGFyZ2VvcEdGRS8lLm1vdmFibGVsaW1pdHNHRkUvJSdhY2NlbnRHRkUvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR1EsMC4zMzMzMzMzZW1GJy1GLDYlUSN5ZEYnRi9GMkZBLUY+Ni1RIjtGJ0ZBRkNGRkZIRkpGTEZORlBGUi9GVlEsMC4yNzc3Nzc4ZW1GJw== LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2JVEpZGF0YWt1dmFGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSomY29sb25lcTtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRj0vJSlzdHJldGNoeUdGPS8lKnN5bW1ldHJpY0dGPS8lKGxhcmdlb3BHRj0vJS5tb3ZhYmxlbGltaXRzR0Y9LyUnYWNjZW50R0Y9LyUnbHNwYWNlR1EsMC4yNzc3Nzc4ZW1GJy8lJ3JzcGFjZUdGTC1GLDYlUSVwbG90RidGL0YyLUkobWZlbmNlZEdGJDYkLUYjNjEtRiw2JVEjeGRGJ0YvRjItRjY2LVEiLEYnRjlGOy9GP0YxRkBGQkZERkZGSC9GS1EmMC4wZW1GJy9GTlEsMC4zMzMzMzMzZW1GJy1GLDYlUSN5ZEYnRi9GMkZaLUYsNiVRJnN0eWxlRidGL0YyLUY2Ni1RIj1GJ0Y5RjtGPkZARkJGREZGRkhGSkZNLUYsNiVRJnBvaW50RidGL0YyRlotRiw2JVEnc3ltYm9sRidGL0YyRmJvLUYsNiVRJmNyb3NzRidGL0YyRlotRiw2JVErc3ltYm9sc2l6ZUYnRi9GMkZiby1JI21uR0YkNiRRIzE4RidGOUY5LUY2Ni1RIjpGJ0Y5RjtGPkZARkJGREZGRkhGSkZN plot([seq([xd[i],yd[i]],i=1..10)],style=point,symbol=cross,symbolsize=18); # Vaihtoehtoinen tapa: pisteiden lista. p:=interp(xd,yd,x); # Ei tarvitse ladata. display(datakuva,plot(p,x=0..10,color=blue)); Interpolaatio pakottaa polynomin usein voimakkaaseen heilahteluun, jolloin "mets"a katoaa n"akyvist"a puilta". Luovutaan vaatimuksesta, ett"a polynomin tulisi kulkea datapisteiden kautta. Valitaan alempiasteinen polynomi, joka sopivassa mieless"a approksimoi parhaiten annettua dataa. Kriteeriksi valitaan usein funktiomallin (t"ass polynomi) ja datapisteiss"a laskettujen erotuksten neli"osumman minimointi. T"at"a kautta menetelm"a liittyy funktioiden "a"ariarvoteemaan. Polynomi-PNS-teht"av"ass"a on annettu datapisteit"a ( LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUklbXN1YkdGJDYlLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYkLUYvNiVRImlGJ0YyRjUvRjZRJ25vcm1hbEYnLyUvc3Vic2NyaXB0c2hpZnRHUSIwRictSSNtb0dGJDYtUSIsRidGPS8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGNC8lKXN0cmV0Y2h5R0ZILyUqc3ltbWV0cmljR0ZILyUobGFyZ2VvcEdGSC8lLm1vdmFibGVsaW1pdHNHRkgvJSdhY2NlbnRHRkgvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR1EsMC4zMzMzMzMzZW1GJy1GLDYlLUYvNiVRInlGJ0YyRjVGOEY/Rj0=) (paljon) ja teht"av"an"a on sovittaa esim. polynomi LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbW9HRiQ2LVEhRicvJSxtYXRodmFyaWFudEdRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y0LyUpc3RyZXRjaHlHRjQvJSpzeW1tZXRyaWNHRjQvJShsYXJnZW9wR0Y0LyUubW92YWJsZWxpbWl0c0dGNC8lJ2FjY2VudEdGNC8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHRkMtRiM2KEYrLUYjNiYtSSNtaUdGJDYlUSJwRicvJSdpdGFsaWNHUSV0cnVlRicvRjBRJ2l0YWxpY0YnLUYsNi1RMCZBcHBseUZ1bmN0aW9uO0YnRi9GMkY1RjdGOUY7Rj1GP0ZBRkQtSShtZmVuY2VkR0YkNiQtRiM2JC1GSzYlUSJ4RidGTkZRRi9GL0YvLUYsNi1RIj1GJ0YvRjJGNUY3RjlGO0Y9Rj8vRkJRLDAuMjc3Nzc3OGVtRicvRkVGXG8tRiM2Ky1JJW1zdWJHRiQ2JS1GSzYlUSJiRidGTkZRLUYjNiQtSSNtbkdGJDYkUSIwRidGL0YvLyUvc3Vic2NyaXB0c2hpZnRHRltwLUYsNi1RIitGJ0YvRjJGNUY3RjlGO0Y9Rj8vRkJRLDAuMjIyMjIyMmVtRicvRkVGYnAtRiM2Ji1GYW82JUZjby1GIzYkLUZpbzYkUSIxRidGL0YvRlxwLUYsNi1RMSZJbnZpc2libGVUaW1lcztGJ0YvRjJGNUY3RjlGO0Y9Rj9GQUZERmVuRi9GXnAtRiM2Jy1GYW82JUZjby1GIzYkLUZpbzYkUSIyRidGL0YvRlxwRl1xLUYjNiQtSSVtc3VwR0YkNiVGZW5GZnEvJTFzdXBlcnNjcmlwdHNoaWZ0R0ZbcEYvRitGL0ZecC1GIzYnLUZhbzYlRmNvLUYjNiQtRmlvNiRRIjNGJ0YvRi9GXHBGXXEtRiM2JC1GXHI2JUZlbkZmckZeckYvRitGL0YrRi9GK0YvRitGLw== annettuun dataan niin, ett"a LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbW9HRiQ2LVEhRicvJSxtYXRodmFyaWFudEdRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y0LyUpc3RyZXRjaHlHRjQvJSpzeW1tZXRyaWNHRjQvJShsYXJnZW9wR0Y0LyUubW92YWJsZWxpbWl0c0dGNC8lJ2FjY2VudEdGNC8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHRkMtRiM2KC1JK211bmRlcm92ZXJHRiQ2Jy1GLDYtUSYmU3VtO0YnRi9GMkY1L0Y4USV0cnVlRidGOS9GPEZPL0Y+Rk9GP0ZBL0ZFUSwwLjE2NjY2NjdlbUYnLUYjNiYtSSNtaUdGJDYlUSJpRicvJSdpdGFsaWNHRk8vRjBRJ2l0YWxpY0YnLUYsNi1RIj1GJ0YvRjJGNUY3RjlGO0Y9Rj8vRkJRLDAuMjc3Nzc3OGVtRicvRkVGXG8tSSNtbkdGJDYkUSIxRidGL0YvLUZXNiVRIm5GJ0ZaRmZuRj8vJSxhY2NlbnR1bmRlckdGNEYrLUknbXNwYWNlR0YkNiYvJSdoZWlnaHRHUSYwLjBleEYnLyUmd2lkdGhHUSQ1LjBGJy8lJmRlcHRoR0ZccC8lKmxpbmVicmVha0dRJWF1dG9GJy1GIzYkLUklbXN1cEdGJDYlLUkobWZlbmNlZEdGJDYkLUYjNictSSVtc3ViR0YkNiUtRlc2JVEieUYnRlpGZm4tRiM2JEZWRi8vJS9zdWJzY3JpcHRzaGlmdEdRIjBGJy1GLDYtUSgmbWludXM7RidGL0YyRjVGN0Y5RjtGPUY/L0ZCUSwwLjIyMjIyMjJlbUYnL0ZFRl5yLUYjNiYtRlc2JVEicEYnRlpGZm4tRiw2LVEwJkFwcGx5RnVuY3Rpb247RidGL0YyRjVGN0Y5RjtGPUY/RkFGRC1GW3E2JC1GIzYkLUZgcTYlLUZXNiVRInhGJ0ZaRmZuRmVxRmdxRi9GL0YvRitGL0YvLUZfbzYkUSIyRidGLy8lMXN1cGVyc2NyaXB0c2hpZnRHRmlxRi9GK0YvRitGLw== minimoituu. Yleens"a datapisteiden lkm =n on paljon suurempi kuin polynomin asteluku m (t"ass"a 3). Jos A on matriisi , jonka alkiot ovat LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbW9HRiQ2LVEhRicvJSxtYXRodmFyaWFudEdRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y0LyUpc3RyZXRjaHlHRjQvJSpzeW1tZXRyaWNHRjQvJShsYXJnZW9wR0Y0LyUubW92YWJsZWxpbWl0c0dGNC8lJ2FjY2VudEdGNC8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHRkMtRiM2Ji1JJW1zdWJHRiQ2JS1JI21pR0YkNiVRImFGJy8lJ2l0YWxpY0dRJXRydWVGJy9GMFEnaXRhbGljRictRiM2Ji1GTDYlUSJpRidGT0ZSLUYsNi1RIixGJ0YvRjIvRjZGUUY3RjlGO0Y9Rj9GQS9GRVEsMC4zMzMzMzMzZW1GJy1GTDYlUSJqRidGT0ZSRi8vJS9zdWJzY3JpcHRzaGlmdEdRIjBGJy1GLDYtUSI9RidGL0YyRjVGN0Y5RjtGPUY/L0ZCUSwwLjI3Nzc3NzhlbUYnL0ZFRmNvLUkobXN1YnN1cEdGJDYnLUZMNiVRInhGJ0ZPRlItRiM2JEZWRi9GaW4vJTFzdXBlcnNjcmlwdHNoaWZ0R0Zeb0Zcb0YvRitGLw== , i=1..n, j=0..3, saadaan b-kertoimien m"a"ar"a"amiseksi ylim"a"ar"aytyv"a yht"al"osysteemi A b = y, miss"a siis A on n x 3 - matriisi. PNS-ratkaisuun p"a"ast"a"an kertomalla t"am"a systeemi A:n transpoosilla, jolloin ratkaistavaksi tulee systeemi (1) (A^T)A b = A^T y (miss"a A^T tarkoittaa A:n transpoosia) (T"am"a johdetaan joko osittaisderivoimalla tai ortogonaaliprojektioon perustuvalla geometrisluonteisella argumentilla.) Jos tuntemattomia kertoimia on yht"a monta kuin datapisteit"a (polynomin asteluku = n-1), niin systeemill"a A b = y on yksik"asitteinen ratkaisu, mik"ali x-pisteet ovat erillisi"a. T"all"oin on kyse interpolaatiopolynomista. Yll"a oleva ratkaisu voidaan Maplella tehd"a rakentamalla matriisi A ja ratkaisemalla systeemi (1) . Maplen LinearAlgebra-kirjastossa on my"os LeastSquares- funktio, jolle annetaan pelkk"a A - matriisi ja datapisteet, joka siis ratkaisee ylim"a"ar"aytyv"an systeemin A b = y PNS-mieless"a. Jotta valinta ei olisi liian helppoa, Maplessa on my"os kirjasto CurveFitting ja siell"a funktio LeastSquares , joka suorittaa sovituksen suoraan datan perusteella, ja k"aytt"aj"all"a on valittavanaan eri menetelmi"a. restart:with(LinearAlgebra): #with(linalg): with(plots):#with(stats): alias(Tr=Transpose,Van=VandermondeMatrix):

Sovita 2. asteen polynomi annettuun dataan xd:=[8,10,12,16,20,30,40,60,100]; yd:=[0.88,1.22,1.64,2.72,3.96,7.66,11.96,21.56,43.16]; n:=nops(xd); C:=Van(xd,n,3); # tai # Matrix((i,j)->xd[i]^(j-1),n,3); M:=Tr(C).C; b:=Tr(C).Vector(yd); a:=LinearSolve(M,b); p:=add(a[i]*x^(i-1),i=1..3); datakuva:=plot([seq([xd[i],yd[i]],i=1..n)],style=point,symbol=cross,symbolsize=18): display(plot(p,x=7..100,color=blue),datakuva); Tarkistetaan LinearAlgebra- ja CurveFitting-pakkauksten LeastSquares-funktioilla. Huomaa, ett"a samannimisten funktioiden tapauksessa on syyt"a k"aytt"a"a pitk"a"a muotoa, muuten toiminta riipuu kirjastojen latausj"arjestyksest"a (onneksi menee virheeseen, eik"a laske jotain p"a"at"ont"a). QyQtJkkuTGluZWFyQWxnZWJyYUc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSS1MZWFzdFNxdWFyZXNHRik2JEkiQ0dGKS1JJ1ZlY3RvckdGJjYjSSN5ZEdGKSIiIg== LSZJLUN1cnZlRml0dGluZ0c2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSS1MZWFzdFNxdWFyZXNHRig2JkkjeGRHRihJI3lkR0YoSSJ4R0YoL0kmY3VydmVHRigsKComSSNjMEdGKCIiIkYuIiIjRjQqJkkjYzFHRihGNEYuRjRGNEkjYzJHRihGNA== JSFH

LUkmYWxpYXNHJSpwcm90ZWN0ZWRHNiMvSSVjb25kRzYiSTBDb25kaXRpb25OdW1iZXJHRig= n:=5: m:=n-2:x:=[$1..n]; C:=Van(x,n,m);cond(Tr(C).C): evalf(%); n:=6: m:=n-2:x:=[$1..n]:C:=Van(x,n,m):cond(Tr(C).C): evalf(%); n:=10: m:=n-2:x:=[$1..n]:C:=Van(x,n,m):cond(Tr(C).C): evalf(%); H"airi"oluku cond mittaa sensitiivisyytt"a virheille (datan virheille ja py"oristysvirheille). Normaaliyht"al"oratkaisun varjopuoli on suuri h"ari"oalttius, jos teht"av"an koko on v"ah"ank"a"an suuri. n:=6: m:=2:x:=[$1..n]:C:=Van(x,n,m):cond(Tr(C).C): evalf(%); n:=20: m:=2:x:=[$1..n]:C:=Van(x,n,m):cond(Tr(C).C): evalf(%); n:=20: m:=3:x:=[$1..n]:C:=Van(x,n,m):cond(Tr(C).C): evalf(%); Paremmat menetelm"at perustuvat ortogonalisointiin tai ns. singulaariarvohajoitelmaan.

Sovita a) 2. asteen PNS-polynomi ja b) interpolaatiopolynomi P:=[[1.4,7400],[1.8,7500],[2.4,7600],[3,7500],[4,7200]];

Nyc3JCQiIzkhIiIiJSt1NyQkIiM9RiYiJSt2NyQkIiNDRiYiJSt3NyQiIiRGKzckIiIlIiUrcw==

xd:=map(lis->lis[1],P); yd:=map(lis->lis[2],P);

NyckIiM5ISIiJCIjPUYlJCIjQ0YlIiIkIiIl

NyciJSt1IiUrdiIlK3dGJCIlK3M=

Jos tuo n"aytt"a"a oudolta, voitaisiin tehd"a funktiot poimieka:=x->x[1]; poimitoka:=x->x[2]; ja suorittaa xd:=map(poimieka,P);yd:=map(poimitoka,P); T"ass"a on k"atev"a ja yleisp"atev"a tapa muodostaa matriisi. C:=Matrix(5,3,(i,j)->xd[i]^(j-1));

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciNXknKW8sInlTdVklPQ==

Tai, kuten edell"a, suoraan LinearAlgebra-kirjaston funktiolla (joka yll"a "aliasoitiinkin" nimelle Van). VandermondeMatrix(xd):%[1..5,1..3]; #

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciNT0qKW8sInlTdVklPQ==

K"aytet"a"an nyt sit"a keskimm"aist"a "harmaata laatikkoa" kertoimet:=LinearAlgebra[LeastSquares](C,Vector(yd));

LSZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSSdjb2x1bW5HRig2Iy9JJCVpZEdGKCI1US1wLCJ5U3VZJT0=

PkkieEc2Ii5GIw==

SSJ4RzYi

poly2:=Tr(kertoimet).<1,x,x^2>;

LChdNDBCOUVGOTczRjE3MThDNyIiIkkieEc2Il00MDg3RDczQTc2Rjg1MzY4KiRGJSIiI11DMDYzODBCMTY2OUM1REQ3

display(plot(P,style=point,symbol=circle,color=black),plot(poly2,x=1.3..4.2));

LSUlUExPVEc2KC0lJ0NVUlZFU0c2JjcnNyQkIiM5ISIiJCImK1MoISIiNyQkIiM9ISIiJCImK10oISIiNyQkIiNDISIiJCImK2coISIiNyQkIiNJISIiJCImK10oISIiNyQkIiNTISIiJCImKz8oISIiLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSdTWU1CT0xHNiMlJ0NJUkNMRUctJSZTVFlMRUc2IyUmUE9JTlRHLSUnQ1VSVkVTRzYkN2R3NyQkIiM4ISIiJCIxXnJDPykpb250ISM3NyQkIjJWJ0dkYXFDOjghIzskIjImNCkqbzwmKjR0dCEjODckJCIyZTpKX1Y4Jkc4ISM7JCIxJCkpeSM9JFt4UCghIzc3JCQiMlYlKW94J0dWVjghIzskIjFHKikzcjIiSFEoISM3NyQkIjAmKiopZj5eJWU4ISM5JCIwUiNIJUhQISlRKCEjNjckJCIxJGY9WjkpUnQ4ISM6JCIxWjdkciZwSVIoISM3NyQkIjJMbUsmemVEKFEiISM7JCIxb10peiZHbihSKCEjNzckJCIybEloN3cvO1MiISM7JCIxSiFlczl3QlMoISM3NyQkIjJvTnJfUldrVCIhIzskIjFCdFg6RjwydSEjNzckJCIyWiRwUVFrQko5ISM7JCIxPFluOWIpPVQoISM3NyQkIjJFXy9IKD1YWTkhIzskIjBUdz8oPm07dSEjNjckJCIyeWE0ZmpgKWY5ISM7JCIxWi16UyM0M1UoISM3NyQkIjI6RmFvJTMlXFoiISM7JCIxOnQpPi02YVUoISM3NyQkIjJNb090KyE0IVwiISM7JCIxMy8lPkJnKkh1ISM3NyQkIjJObUswLipvLzohIzskIjFQTU16a0ZNdSEjNzckJCIyLS4xQU1Zel4iISM7JCIxXFVuRSVRIlF1ISM3NyQkIjI7SWchbzByTDohIzskIjEkW2JgRmZFVyghIzc3JCQiMXpkOk5cMVo6ISM6JCIxZiFcbE5Ha1coISM3NyQkIjJgLzQpUSMqZmk6ISM7JCIxRCR6b25VMlgoISM3NyQkIjJMa0coSChcamQiISM7JCIxd3hiYSgpXGF1ISM3NyQkIjJeLC53Pk85ZiIhIzskIjE+Nz4hKj1iZXUhIzc3JCQiMjpIZUUzLWVnIiEjOyQiMU5EaD5hTWl1ISM3NyQkIjJmO0xtUCJ6PzshIzskIjF5Zih6IlxCbXUhIzc3JCQiMmtFYHdBY1hqIiEjOyQiMVFdIilRXHVwdSEjNzckJCIxRV8vek1TXDshIzokIjEoelwycGtNWighIzc3JCQiMnZbKFwrYyNbbSIhIzskIjFeVilRamJzWighIzc3JCQiMSp5ZERnXSN5OyEjOiQiMWgwNVFfXCFbKCEjNzckJCIyWCdHZHMqXEZwIiEjOyQiMSJvXyVHNCRSWyghIzc3JCQiMlUiR2NTI0h4cSIhIzskIjFCTGBROVQoWyghIzc3JCQiMiRbJ0g+YCRRQTwhIzskIjFcRSVbb1syXCghIzc3JCQiMk9vT0JKaWx0IiEjOyQiMSg0O3kmUSJSXCghIzc3JCQiMTlHY0lhSV88ISM6JCIwdXNsdWF0XCghIzY3JCQiMkJZI1xBOFhtPCEjOyQiMXl5OnIwUSt2ISM3NyQkIjJ5YDI6M2I6eSIhIzskIjImKSlcYmRCYS52ISM4NyQkIjJmPE4/TFRfeiIhIzskIjFHYzxIZk0xdiEjNzckJCIxc1YoKT1TPzU9ISM6JCIxMEFXIzRXJDR2ISM3NyQkIjEuMTdURkdDPSEjOiQiMW5UW1E4NTd2ISM3NyQkIjFQdFlpISoqKlE9ISM6JCIxUmV4UnQiXF4oISM3NyQkIjJAVCNbY3BRYD0hIzskIjBQVj47MHdeKCEjNjckJCIxeWI2WyhcJW89ISM6JCIxKVJqKlwpXC5fKCEjNzckJCIyIm9OcjFxJkgpPSEjOyQiMV5RXiFRRUhfKCEjNzckJCIxKnpmUngjeigqPSEjOiQiMC0zW0wkXER2ISM2NyQkIjJVIkdjJHAwRCI+ISM7JCIxbHJTLjcoel8oISM3NyQkIjJrRWAnUWAtRT4hIzskIjFYKVFaYik9SXYhIzc3JCQiMjI0PXdQPzolPiEjOyQiMWVnIVt1ZkVgKCEjNzckJCIyKEhmPUQkemAmPiEjOyQiMTFibE1sIVtgKCEjNzckJCIyQ1UlKVEmZTpxPiEjOyQiMVpXXFUlSHFgKCEjNzckJCIydVgiSHgpKUglKT4hIzskIjFONGBNSzRSdiEjNzckJCIwKWY+KClHKys/ISM5JCIyWGElPSFvNjhhKCEjODckJCIxc1YoZVMsTywjISM6JCIxJVI2NV9xSmEoISM3NyQkIjAuMWklZS5IPyEjOSQiMU5rXTIvQFh2ISM3NyQkIjJgLTBJbS1KLyMhIzskIjEtLU8mei9xYSghIzc3JCQiMkNaJSopXCYqW2U/ISM7JCIyJlxGYGtvKilbdiEjODckJCIxRmAxSmV3cj8hIzokIjFpJz54LXEvYighIzc3JCQiMi4vM2MnZilvMyMhIzskIjFFLCh5byU+X3YhIzc3JCQiMiNHYzdwN1wsQCEjOyQiMXcoengpR3pgdiEjNzckJCIxZjxOSnEzO0AhIzokIjEoZV5gYUJgYighIzc3JCQiMCwtJSk0SDE4IyEjOSQiMSZ5WzxZI3ljdiEjNzckJCIxI1F3QVoqZldAISM6JCIxeFVsdj03ZXYhIzc3JCQiMjIxN2siPXFmQCEjOyQiMXJPczg4XWZ2ISM3NyQkIjFoQFZYYTR1QCEjOiQiMTxJZnkoXDJjKCEjNzckJCIxMDU/aVNDKj0jISM6JCIxeCQ+YyFSKj5jKCEjNzckJCIybEloIzNzJkg/IyEjOyQiMVtCMCpRZUljKCEjNzckJCIxWCopeSxHNj1BISM6JCIxb1xaemw7a3YhIzc3JCQiMWhAVj4jSEVCIyEjOiQiMjpLeiRIMztsdiEjODckJCIxViYzUEk1ckMjISM6JCIxdmhQVnIzbXYhIzc3JCQiMTQ9T19lQmlBISM6JCIxKTNhRyJbKXBjKCEjNzckJCIxWSNcWy5uaEYjISM6JCIyOl4nW1UlW3hjKCEjODckJCIxJGU7OFpNL0gjISM6JCIxXyRSRnRuJW92ISM3NyQkIjI8SmlrVCg9MUIhIzskIjE+WmBJIik9cHYhIzc3JCQiMVwpcHpUXS9LIyEjOiQiMXIzOCdmdChwdiEjNzckJCIxIjQ9JzNxLk5CISM6JCIxKCk9eChvMS5kKCEjNzckJCIyWSJIZSY+cilcQiEjOyQiMUc0Pk4yeXF2ISM3NyQkIjJqQVgrcTFOTyMhIzskIjFjb2k0ZjpydiEjNzckJCIxYzdEMTwveUIhIzokIjFBaC1cPlxydiEjNzckJCIya0FYNTJtQ1IjISM7JCIxSEctZC13cnYhIzc3JCQiMk9wUXhsJip5UyMhIzskIjF5InpVUnY+ZCghIzc3JCQiMCV6ZSpmNjpVIyEjOSQiMWlqelJONXN2ISM3NyQkIjInKXBSKnA6QVBDISM7JCIxLHddNyZ6QGQoISM3NyQkIjF2XCpwZylSXkMhIzokIjF3JzRVJz49c3YhIzc3JCQiMiVHYzdXZ1VsQyEjOyQiMXYzPWlFN3N2ISM3NyQkIjEmKXBSWCteIVsjISM6JCIxJXoicCJRISo+ZCghIzc3JCQiMlciR2NrXWwmXCMhIzskIjF5elNPaHlydiEjNzckJCIyS2hBJip5TSU0RCEjOyQiMVx2NTwiUTpkKCEjNzckJCIyT29PYENvU18jISM7JCIxdTpaZyk0N2QoISM3NyQkIjIwKyshbyFcJFFEISM7JCIxQFptKjNEM2QoISM3NyQkIjJkOEYvakhRYiMhIzskIjFgKSlHNmhMcXYhIzc3JCQiMlsjXCk0aUNzYyMhIzskIjJsaHc8bV8pcHYhIzg3JCQiMjM0PVlENEVlIyEjOyQiMTB4bUAkRyNwdiEjNzckJCIxXiwuPjsyKGYjISM6JCIyOmtANjF1Jm92ISM4NyQkIjAoUXgjSCVSNkUhIzkkIjE0QGgmemh5YyghIzc3JCQiMWU6SmckR2RpIyEjOiQiMScpKlFbKFszbnYhIzc3JCQiMUFWJz0/Wy5rIyEjOiQiMEc9MVVFaWMoISM2NyQkIjInXCoqKSopKVwlZWwjISM7JCIxS2FdY09DbHYhIzc3JCQiMkNWJ0dEUzBxRSEjOyQiMW5QcC1uRmt2ISM3NyQkIjEiPlEnbywzJW8jISM6JCIxPnF2MC9FanYhIzc3JCQiMjk/U10pUTIqcCMhIzskIjE3KTRPNjFAYyghIzc3JCQiMjg7SzknXC45RiEjOyQiMXhOQTdXKTNjKCEjNzckJCIxXChcZlRqdXMjISM6JCIxaSlISFFHKGZ2ISM3NyQkIjJ4WSRwKVFxTHUjISM7JCIwL2VzPCdHZXYhIzY3JCQiMSc+UmVnIW9jRiEjOiQiMUtcaEooPXFiKCEjNzckJCIxWSJIKSpSW0F4IyEjOiQiMTFpUnBoWWJ2ISM3NyQkIjF4XzBeajUoeSMhIzokIjEsR05uUSJSYighIzc3JCQiMiRvTnJKRlArRyEjOyQiMWxqZUwncENiKCEjNzckJCIybEReVTsjSDpHISM7JCIwMjg/JCl6MmIoISM2NyQkIjJCT3NDXDUuJEchIzskIjF0Smt6JzMhXHYhIzc3JCQiMTA1P1R1RFhHISM6JCIxSm16a2c8WnYhIzc3JCQiMmIyOmc8OiJmRyEjOyQiMVhHVzxaVFh2ISM3NyQkIjE+UHVkU1l0RyEjOiQiMUVNNSh6Rk5hKCEjNzckJCIxcFB2InAuJCkpRyEjOiQiMSxVNGgoMzphKCEjNzckJCIxWiRwW3QmNC5IISM6JCIxOUAkUSJ5VVJ2ISM3NyQkIjFOcFFwNko9SCEjOiQiMXJBI3oxO3NgKCEjNzckJCIwJz5SS0hySkghIzkkIjFpYU5CIjNfYCghIzc3JCQiMk5vT0w5K28lSCEjOyQiMUBvRCNmISlHYCghIzc3JCQiMSc0PlFJXD4nSCEjOiQiMW5AXFA/Wkl2ISM3NyQkIjJkMjpxS1tsKEghIzskIjEmPlokKj0kM0d2ISM3NyQkIjFVJSlvUWMhKSopSCEjOiQiMSl5Z3ZFY2VfKCEjNzckJCIxOUZhaylwYiskISM6JCIxQE8kPidvOEJ2ISM3NyQkIjEiPlE7QkMqPUkhIzokIjFkQS0pXHMyXyghIzc3JCQiMnhYIkhOJmVXLiQhIzskIi90ZmpAJnpeKCEjNTckJCIyYjA2aS00I1tJISM7JCIxX20qRydHUjp2ISM3NyQkIjJ0VSYzJVwmSGpJISM7JCIvWV5GcV43diEjNTckJCIyUHFTIno4bXhJISM7JCIxXi0oemA3KDR2ISM3NyQkIjIkeWI2dDFsIzQkISM7JCIxcDFzJnA8bl0oISM3NyQkIjF6YzhDYlQxSiEjOiQiMTwlPUd1MFJdKCEjNzckJCIyJFF3X3ZGRUBKISM7JCIxJz0vP1IxM10oISM3NyQkIi8qenAqW29PSiEjOCQiMFhjUT85dlwoISM2NyQkIjI5P1MhKiopNCw6JCEjOyQiMSpHcDkjemUlXCghIzc3JCQiMmpGYiFwIzRZOyQhIzskIjF6KHBgR2s4XCghIzc3JCQiMUVfL1AmKWV6SiEjOiQiMEw1J2ZdJ3pbKCEjNjckJCIyMzE3JUdHQyU+JCEjOyQiMSJbJGZFPWQlWyghIzc3JCQiMSc0PikzO1UzSyEjOiQiMWVTdTlcQSJbKCEjNzckJCIybkFYcXNrVEEkISM7JCIxcFxlVl9WeHUhIzc3JCQiMlYoWygqPTFKUUshIzskIjF6T1wpM2tSWighIzc3JCQiMCYqKil6UDlNRCQhIzkkIjFkT0xpKik9cXUhIzc3JCQiMiQpZTwmRzE1bkshIzskIjEuJ1wsbTFuWSghIzc3JCQiMSV5Y2BKaj9HJCEjOiQiMmIhMy4uRiRHWSghIzg3JCQiMmIsLncuVWhIJCEjOyQiMSkqcE0+UTdmdSEjNzckJCIyJ1woXCplJGUzSiQhIzskIjFZWSU+JjM9YnUhIzc3JCQiMUNbJ0hEWV9LJCEjOiQiMXNueVkwRV51ISM3NyQkIjIxKnpmVyE0Lk0kISM7JCIxbWshejgoM1p1ISM3NyQkIjApZj4uaiJbTiQhIzkkIjF4bCEqPTIrVnUhIzc3JCQiMjpAVS8yXydwTCEjOyQiMSVSTVwlUXZRdSEjNzckJCIxRl8vISpcTyVRJCEjOiQiMWlhPzJWWk11ISM3NyQkIjIneWM4TlkpeVIkISM7JCIxMyg9UkkjW0l1ISM3NyQkIjJEXStUbnpMVCQhIzskIjI6KnknRyhvJGVVKCEjODckJCIyRU1vO2lRc1UkISM7JCIxOTomPlo9O1UoISM3NyQkIjFNb09dXixVTSEjOiQiMVtBJCllWjA8dSEjNzckJCIyJXBReHQiZWhYJCEjOyQiMiY0JWUnKSlHaTd1ISM4NyQkIjJCUnlPPWk9WiQhIzskIjE+bkBPKEd3UyghIzc3JCQiMSV5Y0JxZ2FbJCEjOiQiMXgsMSgqPkMudSEjNzckJCIyQlcpb1VeKjNdJCEjOyQiMSpmNVc0JD4pUighIzc3JCQiMnZWKFtmPidcXiQhIzskIjE2Iz0wI29fJFIoISM3NyQkIjElKW9QWSlbLmAkISM6JCIxZFMyVD5OKVEoISM3NyQkIjFSeGFGXmlWTiEjOiQiMWAqZkZfRlFRKCEjNzckJCIyRFghNGlfdWVOISM7JCIwJXpMXHhneXQhIzY3JCQiMVQhM2NjXUxkJCEjOiQiMTkkXCNwelx0dCEjNzckJCIxclQkeUtZemUkISM6JCIxLlIpKkddS290ISM3NyQkIjFBVylbUilbLU8hIzokIjFwMygpM101anQhIzc3JCQiMSV6ZShvKGVraCQhIzokIjFCeHFiIUchZXQhIzc3JCQiMXRXKkc2aDpqJCEjOiQiMVI3MEo3Wl90ISM3NyQkIjF0WCI+dWFmayQhIzokIjJOKytZJCozck0oISM4NyQkIjE9TW9lTDVoTyEjOiQiMXV1aVliUlR0ISM3NyQkIjIoPVB1L2wiW24kISM7JCIxQT9kLj07T3QhIzc3JCQiMnROciMpNHMqKm8kISM7JCIwOCRvdSM0LkwoISM2NyQkIjJQZDlmXilbL1AhIzskIjExY1ZvampDdCEjNzckJCIxYjQ+KydwKj1QISM6JCIxaF9LYTwiKj10ISM3NyQkIjI4QVcpW140TVAhIzskIjA8PE9daUdKKCEjNjckJCIyJWU7TEpqLVtQISM7JCIybCU+Om14QTJ0ISM4NyQkIjJjKiopem5QSGlQISM7JCIxM2FHYFVSLHQhIzc3JCQiMk5zV0hyWSF5UCEjOyQiMSpmPnlkeltIKCEjNzckJCIxaUFYOSg0QnokISM6JCIxcDw/eFUiKilHKCEjNzckJCIxLjA1MGoqbyFRISM6JCIxXEpPJCl6dSNHKCEjNzckJCIxRmAxI1xJPCNRISM6JCIxa05KRyoza0YoISM3NyQkIjImUXdfJypmT05RISM7JCIxWEh6TjpfcXMhIzc3JCQiMiVvT3QtNSEqXFEhIzskIjEibzg7KT09a3MhIzc3JCQiMidRd19uYEtrUSEjOyQiMT1gJT1KRHlEKCEjNzckJCIxMTZBYVx2elEhIzokIjFyWToqKVEmNEQoISM3NyQkIjFfLjInKjNQJCpRISM6JCIxN3ZdXyVHW0MoISM3NyQkIjI4NkFrJzMzNFIhIzskIjE0amUqMypvUHMhIzc3JCQiMndReE0hekRCUiEjOyQiMTkpZXJEIT1KcyEjNzckJCIyMS8zMU0mR1BSISM7JCIwWUgpeiR5WUEoISM2NyQkIjJ4Uno9TXBCJlIhIzskIjEyZT5oIz13QCghIzc3JCQiMUVfL2hWXm5SISM6JCIxY1NCYCJlL0AoISM3NyQkIjJhLTBnMyVIIilSISM7JCIwMUspRzgpUT8oISM2NyQkIjJqND49YUZmKlIhIzskIjFkRklpQCRvPighIzc3JCQiMThDW2skMy0sJSEjOiQiMUZISlMlKSkpKj0oISM3NyQkIjFbJjRwIyopb0RTISM6JCIxZiplJypbKkcjPSghIzc3JCQiMmxMbnUiUjNSUyEjOyQiMS88TVBSbHZyISM3NyQkIjEuMDVeJm9XMCUhIzokIjEuLl0kb2p6OyghIzc3JCQiMk5jN2IiNCQqb1MhIzskIjJsTnlENG4xOyghIzg3JCQiMSNHYyMqZWBLMyUhIzokIjJEZ2dWbHdMOighIzg3JCQiMChSemN3ZSg0JSEjOSQiMUFQKnlLO2c5KCEjNzckJCIxTW5NKVwyQTYlISM6JCIxMCtmO0tXUXIhIzc3JCQiMWlCWiZHL3g3JSEjOiQiMTFyWEdJTUlyISM3NyQkIjFXKW88Szg+OSUhIzokIjFVN3kkNF1HNyghIzc3JCQiMS4xN2wlUmY6JSEjOiQiMSJwYmQiPVI6ciEjNzckJCIxOUVfIj1MNDwlISM6JCIxMGglPTpedDUoISM3NyQkIjJOZDl6RCUqZT0lISM7JCIvcHNwIWUjKjQoISM1NyQkIiNVISIiJCIxcCopKlJxajo0KCEjNy0lJkNPTE9SRzYnJSRSR0JHJCIpQyllcSUhIikkIiIhISIiJCIoJz4hXCYhIiktJStfQVRUUklCVVRFRzYjLyUnc291cmNlRyUsbWF0aGRlZmF1bHRHLSYlJl9BWElTRzYjIiIiNictJStfR1JJRExJTkVTRzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNictJStfR1JJRExJTkVTRzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUrQVhFU0xBQkVMU0c2JS1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRIng2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRITYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1LSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIkSSYhIiItJSlCT1VORFNfWUc2IyQiJD8iISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlK00hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlcU0hIiItJSlDSElMRFJFTkc2Ig==

b) Nyt haetaan interpolaatioplynomia. Sen asteluku = datapisteiden lkm -1 = 4. No siihen l"oytyy taas ihan valmiskin funktio: (Jos ratkaistaisiin LSQ-teht"av"an"a, matriisina olisi koko Vandermonde (transposilla kertominen voitaisiin "supistaa pois").) intpol:=interp(xd,yd,x);

LCwqJEkieEc2IiIiJSQiKy1QX3QhKSEiKSokRiQiIiQkIStoIjMmZSIpISIoKiRGJCIiIyQiK0FPMnhGISInRiQkIStNO3EibyRGMyQiK1QsJylSISpGMyIiIg==

display(plot(intpol,x=1.3..4.2,color=blue),plot(P,style=point,symbol=circle,color=black),plot(poly2,x=1.3..4.2,color=gray));

LSUlUExPVEc2KS0lJ0NVUlZFU0c2JDdkdzckJCIjOCEiIiQiMUA2WTBoMSZRKCEjNzckJCIyVidHZGFxQzo4ISM7JCIxUkFqbmYnb1EoISM3NyQkIjJlOkpfVjgmRzghIzskIjF1Vnl4P2UpUSghIzc3JCQiMlYlKW94J0dWVjghIzskIjFIM05QPG4hUighIzc3JCQiMCYqKilmPl4lZTghIzkkIjEwTDVdJlFIUighIzc3JCQiMSRmPVo5KVJ0OCEjOiQiMXlrO0wnXGBSKCEjNzckJCIyTG1LJnplRChRIiEjOyQiMWU0QlBtcihSKCEjNzckJCIybEloN3cvO1MiISM7JCIxJ2ZXW18lSCt1ISM3NyQkIjJvTnJfUldrVCIhIzskIjFgI0dKbiozLnUhIzc3JCQiMlokcFFRa0JKOSEjOyQiMU8kXGBgKipmUyghIzc3JCQiMkVfL0goPVhZOSEjOyQiMXhdPUdVNjR1ISM3NyQkIjJ5YTRmamApZjkhIzskIjE+OiZwKlEmPlQoISM3NyQkIjI6RmFvJTMlXFoiISM7JCIxImZHXy9fX1QoISM3NyQkIjJNb090KyE0IVwiISM7JCIxTW9bOlxtPXUhIzc3JCQiMk5tSzAuKm8vOiEjOyQiMSRSQFtAVj9VKCEjNzckJCIyLS4xQU1Zel4iISM7JCIwNHFvdSI9RHUhIzY3JCQiMjtJZyFvMHJMOiEjOyQiMWVDeEVbKipHdSEjNzckJCIxemQ6TlwxWjohIzokIjImZWQ+ImUpR0t1ISM4NyQkIjJgLzQpUSMqZmk6ISM7JCIxQzFkKilwPU91ISM3NyQkIjJMa0coSChcamQiISM7JCIxPC0qXFkjcFJ1ISM3NyQkIjJeLC53Pk85ZiIhIzskIjAxbk9vImZWdSEjNjckJCIyOkhlRTMtZWciISM7JCIxeCZRdHpddFcoISM3NyQkIjJmO0xtUCJ6PzshIzskIjE8YTRlXkpedSEjNzckJCIya0Vgd0FjWGoiISM7JCIxSjs+TCYpKVxYKCEjNzckJCIxRV8vek1TXDshIzokIjEuWXdeMSkqZXUhIzc3JCQiMnZbKFwrYyNbbSIhIzskIjFzTSVlSmFKWSghIzc3JCQiMSp5ZERnXSN5OyEjOiQiMT14PkheIW9ZKCEjNzckJCIyWCdHZHMqXEZwIiEjOyQiMTw1JCkqb2gyWighIzc3JCQiMlUiR2NTI0h4cSIhIzskIjF5SzB2KWVbWighIzc3JCQiMiRbJ0g+YCRRQTwhIzskIjFhYHpwOigpeXUhIzc3JCQiMk9vT0JKaWx0IiEjOyQiMWEnUV9dYEZbKCEjNzckJCIxOUdjSWFJXzwhIzokIjEnb1ZQcGRxWyghIzc3JCQiMkJZI1xBOFhtPCEjOyQiMTB5YXRdIjRcKCEjNzckJCIyeWAyOjNiOnkiISM7JCIxM0dWWSI9XVwoISM3NyQkIjJmPE4/TFRfeiIhIzskIjEvWiE9Pz0oKVwoISM3NyQkIjFzVigpPVM/NT0hIzokIjFZdGJhJFJGXSghIzc3JCQiMS4xN1RGR0M9ISM6JCIxJXkqKW8pZlwxdiEjNzckJCIxUHRZaSEqKipRPSEjOiQiMWEseWkwUjV2ISM3NyQkIjJAVCNbY3BRYD0hIzskIjJEJEhbZkI7OXYhIzg3JCQiMXliNlsoXCVvPSEjOiQiMSI0ZWZ5byE9diEjNzckJCIyIm9OcjFxJkgpPSEjOyQiMV9qcSc9J3lAdiEjNzckJCIxKnpmUngjeigqPSEjOiQiMiZ5Y28nM1FiXyghIzg3JCQiMlUiR2MkcDBEIj4hIzskIjImKTMjcHNgP0h2ISM4NyQkIjJrRWAnUWAtRT4hIzskIjFma0g7XF9LdiEjNzckJCIyMjQ9d1A/OiU+ISM7JCIxUWs9dmpFT3YhIzc3JCQiMihIZj1EJHpgJj4hIzskIjF4Qm1NQmJSdiEjNzckJCIyQ1UlKVEmZTpxPiEjOyQiMSJvZWxcKilIYSghIzc3JCQiMnVYIkh4KSlIJSk+ISM7JCIxYHlwdD5AWXYhIzc3JCQiMClmPigpRysrPyEjOSQiMURPV2YjNChcdiEjNzckJCIxc1YoZVMsTywjISM6JCIxPidIRCNmbV92ISM3NyQkIjAuMWklZS5IPyEjOSQiMW9tRkx4JGZiKCEjNzckJCIyYC0wSW0tSi8jISM7JCIxcldTKiopUSlldiEjNzckJCIyQ1olKilcJipbZT8hIzskIjFeXz9JNCM+YyghIzc3JCQiMUZgMUpld3I/ISM6JCIxcDZrbTRda3YhIzc3JCQiMi4vM2MnZilvMyMhIzskIjFrM0RGck1udiEjNzckJCIyI0djN3A3XCxAISM7JCIxOFJNYy8rcXYhIzc3JCQiMWY8TkpxMztAISM6JCIxTyZSc0NiRGQoISM3NyQkIjAsLSUpNEgxOCMhIzkkIjImPiN6LiE+K3Z2ISM4NyQkIjEjUXdBWipmV0AhIzokIjEmUU03SmRzZCghIzc3JCQiMjIxN2siPXFmQCEjOyQiMXJeIiopRyllenYhIzc3JCQiMWhAVlhhNHVAISM6JCIxW3ElPSwvPGUoISM3NyQkIjEwNT9pU0MqPSMhIzokIjBeaGZ6O1FlKCEjNjckJCIybEloIzNzJkg/IyEjOyQiMTMqeklzRWNlKCEjNzckJCIxWCopeSxHNj1BISM6JCIxXipSOF02dmUoISM3NyQkIjFoQFY+I0hFQiMhIzokIjImZUZJXDY/KmUoISM4NyQkIjFWJjNQSTVyQyMhIzokIjFBOFBXQHghZighIzc3JCQiMTQ9T19lQmlBISM6JCIxIikzNHghKkcjZighIzc3JCQiMVkjXFsubmhGIyEjOiQiMUg/YjhFZCRmKCEjNzckJCIxJGU7OFpNL0gjISM6JCIxPCozP3dzWmYoISM3NyQkIjI8SmlrVCg9MUIhIzskIjI6ISlmZiE9J2ZmKCEjODckJCIxXClwelRdL0sjISM6JCIxWCQzbSJRInBmKCEjNzckJCIxIjQ9JzNxLk5CISM6JCIxIkdHZm9qeGYoISM3NyQkIjJZIkhlJj5yKVxCISM7JCIxUipwd1UpXClmKCEjNzckJCIyakFYK3ExTk8jISM7JCIxPWFScHYwKmYoISM3NyQkIjFjN0QxPC95QiEjOiQiMWAqKjNQK2AqZighIzc3JCQiMmtBWDUybUNSIyEjOyQiMU9LR2o+KCkqZighIzc3JCQiMk9wUXhsJip5UyMhIzskIjE0cWQsbzQrdyEjNzckJCIwJXplKmY2OlUjISM5JCIxTzolXF90LGcoISM3NyQkIjInKXBSKnA6QVBDISM7JCIxKClRYFsmPixnKCEjNzckJCIxdlwqcGcpUl5DISM6JCIwX0V4P1IqKmYoISM2NyQkIjIlR2M3V2dVbEMhIzskIjEiKj1QKFJQJypmKCEjNzckJCIxJilwUlgrXiFbIyEjOiQiMV8sajJiPCpmKCEjNzckJCIyVyJHY2tdbCZcIyEjOyQiMWJEI2Y/byYpZighIzc3JCQiMktoQSYqeU0lNEQhIzskIjE3d2VrLyp5ZighIzc3JCQiMk9vT2BDb1NfIyEjOyQiMlhSRkArU3FmKCEjODckJCIyMCsrIW8hXCRRRCEjOyQiMS0leSJcLDMnZighIzc3JCQiMmQ4Ri9qSFFiIyEjOyQiMCRIJTN4JSpbZighIzY3JCQiMlsjXCk0aUNzYyMhIzskIjFKeiZ5YlpQZighIzc3JCQiMjM0PVlENEVlIyEjOyQiMS4/O1I1SCNmKCEjNzckJCIxXiwuPjsyKGYjISM6JCIxKDQuM3IneSFmKCEjNzckJCIwKFF4I0glUjZFISM5JCIxWylwb0FvIiplKCEjNzckJCIxZTpKZyRHZGkjISM6JCIxZHQib28/dWUoISM3NyQkIjFBVic9P1suayMhIzokIjIwZSJcJj0yYmUoISM4NyQkIjInXCoqKSopKVwlZWwjISM7JCIxIz16YjBOTGUoISM3NyQkIjJDVidHRFMwcUUhIzskIjFmU2ZpU0AiZSghIzc3JCQiMSI+USdvLDMlbyMhIzokIjFiIm9udioqKnl2ISM3NyQkIjI5P1NdKVEyKnAjISM7JCIwOHUtNSxsZCghIzY3JCQiMjg7SzknXC45RiEjOyQiMSlwK0RSdFFkKCEjNzckJCIxXChcZlRqdXMjISM6JCIxdnJ1QjlTcnYhIzc3JCQiMnhZJHApUXFMdSMhIzskIjEiPWNuak4kb3YhIzc3JCQiMSc+UmVnIW9jRiEjOiQiMXN4djFvbGx2ISM3NyQkIjFZIkgpKlJbQXgjISM6JCIyYilwSUxRUml2ISM4NyQkIjF4XzBeajUoeSMhIzokIjE3bV1pMTpmdiEjNzckJCIyJG9OckpGUCtHISM7JCIxXyF6anhcaGIoISM3NyQkIjJsRF5VOyNIOkchIzskIjAmM1teIWVFYighIzY3JCQiMkJPc0NcNS4kRyEjOyQiMUpBNVgqPiFcdiEjNzckJCIxMDU/VHVEWEchIzokIjEtciM9SHlfYSghIzc3JCQiMmIyOmc8OiJmRyEjOyQiMXNjJltPLjxhKCEjNzckJCIxPlB1ZFNZdEchIzokIjF5MXNmaCp5YCghIzc3JCQiMXBQdiJwLiQpKUchIzokIjFVc0YveSVRYCghIzc3JCQiMVokcFt0JjQuSCEjOiQiMWg6UHlCcUh2ISM3NyQkIjFOcFFwNko9SCEjOiQiMVV5Oi1oS0R2ISM3NyQkIjAnPlJLSHJKSCEjOSQiMW1AZidIejhfKCEjNzckJCIyTm9PTDkrbyVIISM7JCIxW1YkKlJfJG9eKCEjNzckJCIxJzQ+UUlcPidIISM6JCIxKj1Gc25uQF4oISM3NyQkIjJkMjpxS1tsKEghIzskIjF5JTR5KkdkMnYhIzc3JCQiMVUlKW9RYyEpKilIISM6JCIxJypHcEIvSy52ISM3NyQkIjE5RmFrKXBiKyQhIzokIjFrcSsjcG4iKVwoISM3NyQkIjEiPlE7QkMqPUkhIzokIjEkPTtdWEJQXCghIzc3JCQiMnhYIkhOJmVXLiQhIzskIjFvOjwoXGwlKVsoISM3NyQkIjJiMDZpLTQjW0khIzskIjFqXSd6Klt0JFsoISM3NyQkIjJ0VSYzJVwmSGpJISM7JCIxKFFyR0xsJXl1ISM3NyQkIjJQcVMiejhteEkhIzskIjFmKUhsVHRMWighIzc3JCQiMiR5YjZ0MWwjNCQhIzskIjFGIyp6cHApelkoISM3NyQkIjF6YzhDYlQxSiEjOiQiMHVPZlZ4SFkoISM2NyQkIjIkUXdfdkZFQEohIzskIjEob0VmbzR2WCghIzc3JCQiLyp6cCpbb09KISM4JCIxQUQqXFtqPFgoISM3NyQkIjI5P1MhKiopNCw6JCEjOyQiMXUuYz4kNG5XKCEjNzckJCIyakZiIXAjNFk7JCEjOyQiMSVvSiRIKCo+VHUhIzc3JCQiMUVfL1AmKWV6SiEjOiQiMTs9QXhnWE51ISM3NyQkIjIzMTclR0dDJT4kISM7JCIxWmMlND0hekh1ISM3NyQkIjEnND4pMztVM0shIzokIjFXc2JZeUVDdSEjNzckJCIybkFYcXNrVEEkISM7JCIxakhMI1IlND11ISM3NyQkIjJWKFsoKj0xSlFLISM7JCIxR1kqUSdRXjd1ISM3NyQkIjAmKiopelA5TUQkISM5JCIwTSRIeV1fMXUhIzY3JCQiMiQpZTwmRzE1bkshIzskIjE3Y18jPXY1UyghIzc3JCQiMSV5Y2BKaj9HJCEjOiQiMCpvdVxpNCZSKCEjNjckJCIyYiwudy5VaEgkISM7JCIxPD9jalxYKlEoISM3NyQkIjInXChcKmUkZTNKJCEjOyQiMUc0bzxtYSRRKCEjNzckJCIxQ1snSERZX0skISM6JCIxS08xSEt3eHQhIzc3JCQiMjEqemZXITQuTSQhIzskIjEjb1RxRDE8UCghIzc3JCQiMClmPi5qIltOJCEjOSQiMVV3OWdgKGVPKCEjNzckJCIyOkBVLzJfJ3BMISM7JCIxN3I/ei0jKmZ0ISM3NyQkIjFGXy8hKlxPJVEkISM6JCIxZi8pSF9GU04oISM3NyQkIjIneWM4TlkpeVIkISM7JCIyWFJ2PmpHJ1t0ISM4NyQkIjJEXStUbnpMVCQhIzskIjFvYVZnXllVdCEjNzckJCIyRU1vO2lRc1UkISM7JCIxKSlwcCpRenBMKCEjNzckJCIxTW9PXV4sVU0hIzokIjFcJioqXHpqNkwoISM3NyQkIjIlcFF4dCJlaFgkISM7JCIxNDFdR1xqRHQhIzc3JCQiMkJSeU89aT1aJCEjOyQiMTxCbkVcYT50ISM3NyQkIjEleWNCcWdhWyQhIzokIjEpSDEuVD1WSighIzc3JCQiMkJXKW9VXiozXSQhIzskIjFgLjppX1czdCEjNzckJCIydlYoW2Y+J1xeJCEjOyQiMFtiXXVfSkkoISM2NyQkIjElKW9QWSlbLmAkISM6JCIxekMib3hOdUgoISM3NyQkIjFSeGFGXmlWTiEjOiQiMWknXG8nKW9ESCghIzc3JCQiMkRYITRpX3VlTiEjOyQiMTAwQHdrNShHKCEjNzckJCIxVCEzY2NdTGQkISM6JCIyRClHKzJ3Ij5HKCEjODckJCIxclQkeUtZemUkISM6JCIxJlwjeipSQ29GKCEjNzckJCIxQVcpW1IpWy1PISM6JCIxSkw2dXolPUYoISM3NyQkIjElemUobyhla2gkISM6JCIxQih5PkdsckUoISM3NyQkIjF0VypHNmg6aiQhIzokIjFFM0MlKnpAaXMhIzc3JCQiMXRYIj51YWZrJCEjOiQiMTIjUlFXP3dEKCEjNzckJCIxPU1vZUw1aE8hIzokIjB2XDUmRyJIRCghIzY3JCQiMig9UHUvbCJbbiQhIzskIi8iWzpgdChbcyEjNTckJCIydE5yIyk0cyoqbyQhIzskIjJidFgmejRNV3MhIzg3JCQiMlBkOWZeKVsvUCEjOyQiMUNZQDRCQ1NzISM3NyQkIjFiND4rJ3AqPVAhIzokIjJiWWBzKFFJT3MhIzg3JCQiMjhBVylbXjRNUCEjOyQiMWtpOyIpek5LcyEjNzckJCIyJWU7TEpqLVtQISM7JCIxaDhvJz0jKSlHcyEjNzckJCIyYyoqKXpuUEhpUCEjOyQiMSxJeSYpb1tEcyEjNzckJCIyTnNXSHJZIXlQISM7JCIxOSVlO2JSPkEoISM3NyQkIjFpQVg5KDRCeiQhIzokIjF4c25XeCIqPXMhIzc3JCQiMS4wNTBqKm8hUSEjOiQiMSZHZkp5QGdAKCEjNzckJCIxRmAxI1xJPCNRISM6JCIydjlwOTcnRzhzISM4NyQkIjImUXdfJypmT05RISM7JCIxVDhXdVonNEAoISM3NyQkIjIlb090LTUhKlxRISM7JCIxXlh6UzpxM3MhIzc3JCQiMidRd19uYEtrUSEjOyQiMUJ4TjQneW0/KCEjNzckJCIxMTZBYVx2elEhIzokIjFOYyo9W3BaPyghIzc3JCQiMV8uMicqM1AkKlEhIzokIjEqKW81dTRKLnMhIzc3JCQiMjg2QWsnMzM0UiEjOyQiMWV4M1QyIT4/KCEjNzckJCIyd1F4TSF6REJSISM7JCIxKEhfMCpwKTM/KCEjNzckJCIyMS8zMU0mR1BSISM7JCIwaVR2WEwsPyghIzY3JCQiMnhSej1NcEImUiEjOyQiMXlKd0goMycqPighIzc3JCQiMUVfL2hWXm5SISM6JCIxRE5LKCkpKVEqPighIzc3JCQiMmEtMGczJUgiKVIhIzskIjFIIykpPkVrJSo+KCEjNzckJCIyajQ+PWFGZipSISM7JCIyOnd5UjFTKSo+KCEjODckJCIxOENbayQzLSwlISM6JCIxLjVdYykzMD8oISM3NyQkIjFbJjRwIyopb0RTISM6JCIxXGhqLzFlLHMhIzc3JCQiMmxMbnUiUjNSUyEjOyQiMTdvJ2Y3Mkc/KCEjNzckJCIxLjA1XiZvVzAlISM6JCIxS2BhPHljL3MhIzc3JCQiMk5jN2IiNCQqb1MhIzskIjJYeThhR3ZsPyghIzg3JCQiMSNHYyMqZWBLMyUhIzokIjFsXWEpKikzKjNzISM3NyQkIjAoUnpjd2UoNCUhIzkkIjJYbnA3XShmNnMhIzg3JCQiMU1uTSlcMkE2JSEjOiQiMVVEJ3onR3I5cyEjNzckJCIxaUJaJkcveDclISM6JCIwI3kuN3JWPXMhIzY3JCQiMVcpbzxLOD45JSEjOiQiMTQoeipvSkNBcyEjNzckJCIxLjE3bCVSZjolISM6JCIxM2hOaGpQRXMhIzc3JCQiMTlFXyI9TDQ8JSEjOiQiMSpccyR6eEBKcyEjNzckJCIyTmQ5ekQlKmU9JSEjOyQiMXBwLWtbXE9zISM3NyQkIiNVISIiJCIyJil5Ils7IikpPUMoISM4LSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIzUhIiItJSdDVVJWRVNHNiY3JzckJCIjOSEiIiQiJitTKCEiIjckJCIjPSEiIiQiJitdKCEiIjckJCIjQyEiIiQiJitnKCEiIjckJCIjSSEiIiQiJitdKCEiIjckJCIjUyEiIiQiJis/KCEiIi0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUnU1lNQk9MRzYjJSdDSVJDTEVHLSUmU1RZTEVHNiMlJlBPSU5URy0lJ0NVUlZFU0c2JDdkdzckJCIjOCEiIiQiMV5yQz8pKW9udCEjNzckJCIyVidHZGFxQzo4ISM7JCIyJjQpKm88Jio0dHQhIzg3JCQiMmU6Sl9WOCZHOCEjOyQiMSQpKXkjPSRbeFAoISM3NyQkIjJWJSlveCdHVlY4ISM7JCIxRyopM3IyIkhRKCEjNzckJCIwJioqKWY+XiVlOCEjOSQiMFIjSCVIUCEpUSghIzY3JCQiMSRmPVo5KVJ0OCEjOiQiMVo3ZHImcElSKCEjNzckJCIyTG1LJnplRChRIiEjOyQiMW9dKXomR24oUighIzc3JCQiMmxJaDd3LztTIiEjOyQiMUohZXM5d0JTKCEjNzckJCIyb05yX1JXa1QiISM7JCIxQnRYOkY8MnUhIzc3JCQiMlokcFFRa0JKOSEjOyQiMTxZbjliKT1UKCEjNzckJCIyRV8vSCg9WFk5ISM7JCIwVHc/KD5tO3UhIzY3JCQiMnlhNGZqYClmOSEjOyQiMVotelMjNDNVKCEjNzckJCIyOkZhbyUzJVxaIiEjOyQiMTp0KT4tNmFVKCEjNzckJCIyTW9PdCshNCFcIiEjOyQiMTMvJT5CZypIdSEjNzckJCIyTm1LMC4qby86ISM7JCIxUE1NemtGTXUhIzc3JCQiMi0uMUFNWXpeIiEjOyQiMVxVbkUlUSJRdSEjNzckJCIyO0lnIW8wckw6ISM7JCIxJFtiYEZmRVcoISM3NyQkIjF6ZDpOXDFaOiEjOiQiMWYhXGxOR2tXKCEjNzckJCIyYC80KVEjKmZpOiEjOyQiMUQkem9uVTJYKCEjNzckJCIyTGtHKEgoXGpkIiEjOyQiMXd4YmEoKVxhdSEjNzckJCIyXiwudz5POWYiISM7JCIxPjc+ISo9YmV1ISM3NyQkIjI6SGVFMy1lZyIhIzskIjFORGg+YU1pdSEjNzckJCIyZjtMbVAiej87ISM7JCIxeWYoeiJcQm11ISM3NyQkIjJrRWB3QWNYaiIhIzskIjFRXSIpUVx1cHUhIzc3JCQiMUVfL3pNU1w7ISM6JCIxKHpcMnBrTVooISM3NyQkIjJ2WyhcK2MjW20iISM7JCIxXlYpUWpic1ooISM3NyQkIjEqeWREZ10jeTshIzokIjFoMDVRX1whWyghIzc3JCQiMlgnR2RzKlxGcCIhIzskIjEib18lRzQkUlsoISM3NyQkIjJVIkdjUyNIeHEiISM7JCIxQkxgUTlUKFsoISM3NyQkIjIkWydIPmAkUUE8ISM7JCIxXEUlW29bMlwoISM3NyQkIjJPb09CSmlsdCIhIzskIjEoNDt5JlEiUlwoISM3NyQkIjE5R2NJYUlfPCEjOiQiMHVzbHVhdFwoISM2NyQkIjJCWSNcQThYbTwhIzskIjF5eTpyMFErdiEjNzckJCIyeWAyOjNiOnkiISM7JCIyJikpXGJkQmEudiEjODckJCIyZjxOP0xUX3oiISM7JCIxR2M8SGZNMXYhIzc3JCQiMXNWKCk9Uz81PSEjOiQiMTBBVyM0VyQ0diEjNzckJCIxLjE3VEZHQz0hIzokIjFuVFtRODU3diEjNzckJCIxUHRZaSEqKipRPSEjOiQiMVJleFJ0IlxeKCEjNzckJCIyQFQjW2NwUWA9ISM7JCIwUFY+OzB3XighIzY3JCQiMXliNlsoXCVvPSEjOiQiMSlSaipcKVwuXyghIzc3JCQiMiJvTnIxcSZIKT0hIzskIjFeUV4hUUVIXyghIzc3JCQiMSp6ZlJ4I3ooKj0hIzokIjAtM1tMJFxEdiEjNjckJCIyVSJHYyRwMEQiPiEjOyQiMWxyUy43KHpfKCEjNzckJCIya0VgJ1FgLUU+ISM7JCIxWClRWmIpPUl2ISM3NyQkIjIyND13UD86JT4hIzskIjFlZyFbdWZFYCghIzc3JCQiMihIZj1EJHpgJj4hIzskIjExYmxNbCFbYCghIzc3JCQiMkNVJSlRJmU6cT4hIzskIjFaV1xVJUhxYCghIzc3JCQiMnVYIkh4KSlIJSk+ISM7JCIxTjRgTUs0UnYhIzc3JCQiMClmPigpRysrPyEjOSQiMlhhJT0hbzY4YSghIzg3JCQiMXNWKGVTLE8sIyEjOiQiMSVSNjVfcUphKCEjNzckJCIwLjFpJWUuSD8hIzkkIjFOa10yL0BYdiEjNzckJCIyYC0wSW0tSi8jISM7JCIxLS1PJnovcWEoISM3NyQkIjJDWiUqKVwmKltlPyEjOyQiMiZcRmBrbyopW3YhIzg3JCQiMUZgMUpld3I/ISM6JCIxaSc+eC1xL2IoISM3NyQkIjIuLzNjJ2YpbzMjISM7JCIxRSwoeW8lPl92ISM3NyQkIjIjR2M3cDdcLEAhIzskIjF3KHp4KUd6YHYhIzc3JCQiMWY8TkpxMztAISM6JCIxKGVeYGFCYGIoISM3NyQkIjAsLSUpNEgxOCMhIzkkIjEmeVs8WSN5Y3YhIzc3JCQiMSNRd0FaKmZXQCEjOiQiMXhVbHY9N2V2ISM3NyQkIjIyMTdrIj1xZkAhIzskIjFyT3M4OF1mdiEjNzckJCIxaEBWWGE0dUAhIzokIjE8SWZ5KFwyYyghIzc3JCQiMTA1P2lTQyo9IyEjOiQiMXgkPmMhUio+YyghIzc3JCQiMmxJaCMzcyZIPyMhIzskIjFbQjAqUWVJYyghIzc3JCQiMVgqKXksRzY9QSEjOiQiMW9cWnpsO2t2ISM3NyQkIjFoQFY+I0hFQiMhIzokIjI6S3okSDM7bHYhIzg3JCQiMVYmM1BJNXJDIyEjOiQiMXZoUFZyM212ISM3NyQkIjE0PU9fZUJpQSEjOiQiMSkzYUciWylwYyghIzc3JCQiMVkjXFsubmhGIyEjOiQiMjpeJ1tVJVt4YyghIzg3JCQiMSRlOzhaTS9IIyEjOiQiMV8kUkZ0biVvdiEjNzckJCIyPEppa1QoPTFCISM7JCIxPlpgSSIpPXB2ISM3NyQkIjFcKXB6VF0vSyMhIzokIjFyMzgnZnQocHYhIzc3JCQiMSI0PSczcS5OQiEjOiQiMSgpPXgobzEuZCghIzc3JCQiMlkiSGUmPnIpXEIhIzskIjFHND5OMnlxdiEjNzckJCIyakFYK3ExTk8jISM7JCIxY29pNGY6cnYhIzc3JCQiMWM3RDE8L3lCISM6JCIxQWgtXD5ccnYhIzc3JCQiMmtBWDUybUNSIyEjOyQiMUhHLWQtd3J2ISM3NyQkIjJPcFF4bCYqeVMjISM7JCIxeSJ6VVJ2PmQoISM3NyQkIjAlemUqZjY6VSMhIzkkIjFpanpSTjVzdiEjNzckJCIyJylwUipwOkFQQyEjOyQiMSx3XTcmekBkKCEjNzckJCIxdlwqcGcpUl5DISM6JCIxdyc0VSc+PXN2ISM3NyQkIjIlR2M3V2dVbEMhIzskIjF2Mz1pRTdzdiEjNzckJCIxJilwUlgrXiFbIyEjOiQiMSV6InAiUSEqPmQoISM3NyQkIjJXIkdja11sJlwjISM7JCIxeXpTT2h5cnYhIzc3JCQiMktoQSYqeU0lNEQhIzskIjFcdjU8IlE6ZCghIzc3JCQiMk9vT2BDb1NfIyEjOyQiMXU6WmcpNDdkKCEjNzckJCIyMCsrIW8hXCRRRCEjOyQiMUBabSozRDNkKCEjNzckJCIyZDhGL2pIUWIjISM7JCIxYCkpRzZoTHF2ISM3NyQkIjJbI1wpNGlDc2MjISM7JCIybGh3PG1fKXB2ISM4NyQkIjIzND1ZRDRFZSMhIzskIjEweG1AJEcjcHYhIzc3JCQiMV4sLj47MihmIyEjOiQiMjprQDYxdSZvdiEjODckJCIwKFF4I0glUjZFISM5JCIxNEBoJnpoeWMoISM3NyQkIjFlOkpnJEdkaSMhIzokIjEnKSpRWyhbM252ISM3NyQkIjFBVic9P1suayMhIzokIjBHPTFVRWljKCEjNjckJCIyJ1wqKikqKSlcJWVsIyEjOyQiMUthXWNPQ2x2ISM3NyQkIjJDVidHRFMwcUUhIzskIjFuUHAtbkZrdiEjNzckJCIxIj5RJ28sMyVvIyEjOiQiMT5xdjAvRWp2ISM3NyQkIjI5P1NdKVEyKnAjISM7JCIxNyk0TzYxQGMoISM3NyQkIjI4O0s5J1wuOUYhIzskIjF4TkE3VykzYyghIzc3JCQiMVwoXGZUanVzIyEjOiQiMWkpSEhRRyhmdiEjNzckJCIyeFkkcClRcUx1IyEjOyQiMC9lczwnR2V2ISM2NyQkIjEnPlJlZyFvY0YhIzokIjFLXGhKKD1xYighIzc3JCQiMVkiSCkqUltBeCMhIzokIjExaVJwaFlidiEjNzckJCIxeF8wXmo1KHkjISM6JCIxLEdOblEiUmIoISM3NyQkIjIkb05ySkZQK0chIzskIjFsamVMJ3BDYighIzc3JCQiMmxEXlU7I0g6RyEjOyQiMDI4PyQpejJiKCEjNjckJCIyQk9zQ1w1LiRHISM7JCIxdEpreiczIVx2ISM3NyQkIjEwNT9UdURYRyEjOiQiMUptemtnPFp2ISM3NyQkIjJiMjpnPDoiZkchIzskIjFYR1c8WlRYdiEjNzckJCIxPlB1ZFNZdEchIzokIjFFTTUoekZOYSghIzc3JCQiMXBQdiJwLiQpKUchIzokIjEsVTRoKDM6YSghIzc3JCQiMVokcFt0JjQuSCEjOiQiMTlAJFEieVVSdiEjNzckJCIxTnBRcDZKPUghIzokIjFyQSN6MTtzYCghIzc3JCQiMCc+UktIckpIISM5JCIxaWFOQiIzX2AoISM3NyQkIjJOb09MOStvJUghIzskIjFAb0QjZiEpR2AoISM3NyQkIjEnND5RSVw+J0ghIzokIjFuQFxQP1pJdiEjNzckJCIyZDI6cUtbbChIISM7JCIxJj5aJCo9JDNHdiEjNzckJCIxVSUpb1FjISkqKUghIzokIjEpeWd2RWNlXyghIzc3JCQiMTlGYWspcGIrJCEjOiQiMUBPJD4nbzhCdiEjNzckJCIxIj5RO0JDKj1JISM6JCIxZEEtKVxzMl8oISM3NyQkIjJ4WCJITiZlVy4kISM7JCIvdGZqQCZ6XighIzU3JCQiMmIwNmktNCNbSSEjOyQiMV9tKkcnR1I6diEjNzckJCIydFUmMyVcJkhqSSEjOyQiL1leRnFeN3YhIzU3JCQiMlBxUyJ6OG14SSEjOyQiMV4tKHpgNyg0diEjNzckJCIyJHliNnQxbCM0JCEjOyQiMXAxcyZwPG5dKCEjNzckJCIxemM4Q2JUMUohIzokIjE8JT1HdTBSXSghIzc3JCQiMiRRd192RkVASiEjOyQiMSc9Lz9SMTNdKCEjNzckJCIvKnpwKltvT0ohIzgkIjBYY1E/OXZcKCEjNjckJCIyOT9TISoqKTQsOiQhIzskIjEqR3A5I3plJVwoISM3NyQkIjJqRmIhcCM0WTskISM7JCIxeihwYEdrOFwoISM3NyQkIjFFXy9QJillekohIzokIjBMNSdmXSd6WyghIzY3JCQiMjMxNyVHR0MlPiQhIzskIjEiWyRmRT1kJVsoISM3NyQkIjEnND4pMztVM0shIzokIjFlU3U5XEEiWyghIzc3JCQiMm5BWHFza1RBJCEjOyQiMXBcZVZfVnh1ISM3NyQkIjJWKFsoKj0xSlFLISM7JCIxek9cKTNrUlooISM3NyQkIjAmKiopelA5TUQkISM5JCIxZE9MaSopPXF1ISM3NyQkIjIkKWU8JkcxNW5LISM7JCIxLidcLG0xblkoISM3NyQkIjEleWNgSmo/RyQhIzokIjJiITMuLkYkR1koISM4NyQkIjJiLC53LlVoSCQhIzskIjEpKnBNPlE3ZnUhIzc3JCQiMidcKFwqZSRlM0okISM7JCIxWVklPiYzPWJ1ISM3NyQkIjFDWydIRFlfSyQhIzokIjFzbnlZMEVedSEjNzckJCIyMSp6ZlchNC5NJCEjOyQiMW1rIXo4KDNadSEjNzckJCIwKWY+LmoiW04kISM5JCIxeGwhKj0yK1Z1ISM3NyQkIjI6QFUvMl8ncEwhIzskIjElUk1cJVF2UXUhIzc3JCQiMUZfLyEqXE8lUSQhIzokIjFpYT8yVlpNdSEjNzckJCIyJ3ljOE5ZKXlSJCEjOyQiMTMoPVJJI1tJdSEjNzckJCIyRF0rVG56TFQkISM7JCIyOip5J0cobyRlVSghIzg3JCQiMkVNbztpUXNVJCEjOyQiMTk6Jj5aPTtVKCEjNzckJCIxTW9PXV4sVU0hIzokIjFbQSQpZVowPHUhIzc3JCQiMiVwUXh0ImVoWCQhIzskIjImNCVlJykpR2k3dSEjODckJCIyQlJ5Tz1pPVokISM7JCIxPm5ATyhHd1MoISM3NyQkIjEleWNCcWdhWyQhIzokIjF4LDEoKj5DLnUhIzc3JCQiMkJXKW9VXiozXSQhIzskIjEqZjVXNCQ+KVIoISM3NyQkIjJ2VihbZj4nXF4kISM7JCIxNiM9MCNvXyRSKCEjNzckJCIxJSlvUFkpWy5gJCEjOiQiMWRTMlQ+TilRKCEjNzckJCIxUnhhRl5pVk4hIzokIjFgKmZGX0ZRUSghIzc3JCQiMkRYITRpX3VlTiEjOyQiMCV6TFx4Z3l0ISM2NyQkIjFUITNjY11MZCQhIzokIjE5JFwjcHpcdHQhIzc3JCQiMXJUJHlLWXplJCEjOiQiMS5SKSpHXUtvdCEjNzckJCIxQVcpW1IpWy1PISM6JCIxcDMoKTNdNWp0ISM3NyQkIjElemUobyhla2gkISM6JCIxQnhxYiFHIWV0ISM3NyQkIjF0VypHNmg6aiQhIzokIjFSNzBKN1pfdCEjNzckJCIxdFgiPnVhZmskISM6JCIyTisrWSQqM3JNKCEjODckJCIxPU1vZUw1aE8hIzokIjF1dWlZYlJUdCEjNzckJCIyKD1QdS9sIltuJCEjOyQiMUE/ZC49O090ISM3NyQkIjJ0TnIjKTRzKipvJCEjOyQiMDgkb3UjNC5MKCEjNjckJCIyUGQ5Zl4pWy9QISM7JCIxMWNWb2pqQ3QhIzc3JCQiMWI0PisncCo9UCEjOiQiMWhfS2E8Iio9dCEjNzckJCIyOEFXKVteNE1QISM7JCIwPDxPXWlHSighIzY3JCQiMiVlO0xKai1bUCEjOyQiMmwlPjpteEEydCEjODckJCIyYyoqKXpuUEhpUCEjOyQiMTNhR2BVUix0ISM3NyQkIjJOc1dIclkheVAhIzskIjEqZj55ZHpbSCghIzc3JCQiMWlBWDkoNEJ6JCEjOiQiMXA8P3hVIiopRyghIzc3JCQiMS4wNTBqKm8hUSEjOiQiMVxKTyQpenUjRyghIzc3JCQiMUZgMSNcSTwjUSEjOiQiMWtOSkcqM2tGKCEjNzckJCIyJlF3XycqZk9OUSEjOyQiMVhIek46X3FzISM3NyQkIjIlb090LTUhKlxRISM7JCIxIm84Oyk9PWtzISM3NyQkIjInUXdfbmBLa1EhIzskIjE9YCU9SkR5RCghIzc3JCQiMTE2QWFcdnpRISM6JCIxclk6KilRJjREKCEjNzckJCIxXy4yJyozUCQqUSEjOiQiMTd2XV8lR1tDKCEjNzckJCIyODZBayczMzRSISM7JCIxNGplKjMqb1BzISM3NyQkIjJ3UXhNIXpEQlIhIzskIjE5KWVyRCE9SnMhIzc3JCQiMjEvMzFNJkdQUiEjOyQiMFlIKXokeVlBKCEjNjckJCIyeFJ6PU1wQiZSISM7JCIxMmU+aCM9d0AoISM3NyQkIjFFXy9oVl5uUiEjOiQiMWNTQmAiZS9AKCEjNzckJCIyYS0wZzMlSCIpUiEjOyQiMDFLKUc4KVE/KCEjNjckJCIyajQ+PWFGZipSISM7JCIxZEZJaUAkbz4oISM3NyQkIjE4Q1trJDMtLCUhIzokIjFGSEpTJSkpKSo9KCEjNzckJCIxWyY0cCMqKW9EUyEjOiQiMWYqZScqWypHIz0oISM3NyQkIjJsTG51IlIzUlMhIzskIjEvPE1QUmx2ciEjNzckJCIxLjA1XiZvVzAlISM6JCIxLi5dJG9qejsoISM3NyQkIjJOYzdiIjQkKm9TISM7JCIybE55RDRuMTsoISM4NyQkIjEjR2MjKmVgSzMlISM6JCIyRGdnVmx3TDooISM4NyQkIjAoUnpjd2UoNCUhIzkkIjFBUCp5SztnOSghIzc3JCQiMU1uTSlcMkE2JSEjOiQiMTArZjtLV1FyISM3NyQkIjFpQlomRy94NyUhIzokIjExclhHSU1JciEjNzckJCIxVylvPEs4PjklISM6JCIxVTd5JDRdRzcoISM3NyQkIjEuMTdsJVJmOiUhIzokIjEicGJkIj1SOnIhIzc3JCQiMTlFXyI9TDQ8JSEjOiQiMTBoJT06XnQ1KCEjNzckJCIyTmQ5ekQlKmU9JSEjOyQiL3BzcCFlIyo0KCEjNTckJCIjVSEiIiQiMXAqKSpScWo6NCghIzctJSZDT0xPUkc2JiUkUkdCRyQiKDclSHYhIigkIig3JUh2ISIoJCIoNyVIdiEiKC0mJSZfQVhJU0c2IyIiIjYnLSUrX0dSSURMSU5FU0c2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiIzYnLSUrX0dSSURMSU5FU0c2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lK0FYRVNMQUJFTFNHNiUtSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJ4NiIvJSdmYW1pbHlHUShERUZBVUxUNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiUSE2Ii0lJUZPTlRHNiUlKERFRkFVTFRHJShERUZBVUxURyIjNS0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJEkmISIiLSUpQk9VTkRTX1lHNiMkIiQ1IiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStNISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJXFNISIiLSUpQ0hJTERSRU5HNiI=

Hyv"a tiet"a"a: K"aytt"aj"a voi valita menetelm"an CurveFitting-pakkauksen LeastSquares-funktiossa.

Polynomin asteluvun nostaminen lis"a"a dramaattisesti h"airi"oalttiutta.

TTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2NzQ0MDc4MTAxNjg4Njc4WCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjMCImIiQkIiM1ISIiRiZGJiIiIkYpJCIjOUYoJCIjPUYoJCIjQ0YoIiIkIiIlJCIkJz4hIiMkIiRDJEY0JCIkdyZGNCIiKiIjO0YlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2NzQ0MDc4MTAxNjg4OTE4WCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjMCImIiQkIiM1ISIiRiZGJiIiIkYpJCIjOUYoJCIjPUYoJCIjQ0YoIiIkIiIlJCIkJz4hIiMkIiRDJEY0JCIkdyZGNCIiKiIjO0YlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2NzQ0MDc4MTAxNjkwMjM4WColKWFueXRoaW5nRzYiRiVbZ2wnIyUhISEiJCIkNDBCOUVGOTczRjE3MThDNzQwODdENzNBNzZGODUzNjhDMDYzODBCMTY2OUM1REQ3RiU=