diff --git a/mp5/Project_5_Maggioni_Claudio/code_template.m b/mp5/Project_5_Maggioni_Claudio/code_template.m index 1ceabd1..0e173ed 100644 --- a/mp5/Project_5_Maggioni_Claudio/code_template.m +++ b/mp5/Project_5_Maggioni_Claudio/code_template.m @@ -7,7 +7,6 @@ load('blur_data/B.mat'); B=double(B); load('blur_data/A.mat'); A=double(A); -ciao = A; % Show Image figure @@ -20,20 +19,18 @@ title('Blured Image') b=B'; b=b(:); -AT = A * A'; +AT = A' * A; bt = A' * b; IL = ichol(AT, struct('type', 'nofill', 'diagcomp', 0.01)); -y = IL \ b; -x0 = IL' \ y; -[x, rvec] = myCG(AT, bt, x0, 200, 1e-6); +[x, rvec] = myCG(AT, bt, zeros(size(b)), 200, 1e-6); [x2, ~, ~, ~, rvec2] = pcg(AT, bt, 1e-6, 200, IL, IL'); figure; -semilogy(rvec); +semilogy(rvec / norm(bt)); hold on; -semilogy(rvec2); +semilogy(rvec2 / norm(bt)); hold off; title('Residual norms over iteration (y is log)') matlab2tikz('showInfo', false, '../res_log.tex'); diff --git a/mp5/img_my-1.png b/mp5/img_my-1.png new file mode 100644 index 0000000..78c7d0e Binary files /dev/null and b/mp5/img_my-1.png differ diff --git a/mp5/img_my.tex b/mp5/img_my.tex index 7d15c02..c32c982 100644 --- a/mp5/img_my.tex +++ b/mp5/img_my.tex @@ -6,15 +6,15 @@ \begin{tikzpicture} \begin{axis}[% -width=4.521in, -height=3.566in, -at={(0.758in,0.481in)}, +width=6.028in, +height=4.754in, +at={(1.011in,0.642in)}, scale only axis, xmin=0, -xmax=250, +xmax=160, ymode=log, -ymin=1e-06, -ymax=118.544185368288, +ymin=1e-09, +ymax=1, yminorticks=true, axis background/.style={fill=white}, title style={font=\bfseries}, @@ -22,236 +22,187 @@ title={Sharp Image (myCG)} ] \addplot [color=mycolor1, forget plot] table[row sep=crcr]{% -1 11.1886780605947\\ -2 37.688016651504\\ -3 20.5612191459901\\ -4 10.1126098486914\\ -5 6.60295496425051\\ -6 4.03660691189902\\ -7 2.83696084221764\\ -8 2.02776259797851\\ -9 1.55623996384361\\ -10 1.26373343652138\\ -11 0.968080560648255\\ -12 0.808014633949812\\ -13 0.666010141343337\\ -14 0.547935882226161\\ -15 0.455778233603678\\ -16 0.379140935249178\\ -17 0.325385471261729\\ -18 0.27468130551472\\ -19 0.241318492631058\\ -20 0.203089896725741\\ -21 0.186402910640978\\ -22 0.164285195168731\\ -23 0.161498606695505\\ -24 0.13554465057248\\ -25 0.108256329622331\\ -26 0.096685060868788\\ -27 0.0901032819163635\\ -28 0.0846560423043152\\ -29 0.0720452421305809\\ -30 0.0638081561682231\\ -31 0.0586938472203618\\ -32 0.0524931123015411\\ -33 0.0500311324114693\\ -34 0.0441377829780386\\ -35 0.0430660141184194\\ -36 0.0399007401033348\\ -37 0.0357338939687774\\ -38 0.0334582949737502\\ -39 0.0301143079798452\\ -40 0.0277450095617363\\ -41 0.0246457162237748\\ -42 0.0226451015644127\\ -43 0.0214585559097522\\ -44 0.0197750620834928\\ -45 0.0184656387170807\\ -46 0.01424126114094\\ -47 0.014925314691824\\ -48 0.0138939777471367\\ -49 0.0124854196790047\\ -50 0.0121443393806818\\ -51 0.0111248201462271\\ -52 0.0117048407914948\\ -53 0.0113323279545552\\ -54 0.0106729739632728\\ -55 0.00944425804326873\\ -56 0.00871066933482761\\ -57 0.00881210600673333\\ -58 0.00795156096327054\\ -59 0.00726613636631538\\ -60 0.00700470814378812\\ -61 0.00621935380743637\\ -62 0.00580029695866593\\ -63 0.00557779006381989\\ -64 0.00510498070037545\\ -65 0.0047009622833176\\ -66 0.00438789538442034\\ -67 0.00435585711109691\\ -68 0.00396640165581629\\ -69 0.00393477324085466\\ -70 0.00371617606606877\\ -71 0.00327886320159879\\ -72 0.00312237367539779\\ -73 0.00299349690875735\\ -74 0.00284074963413765\\ -75 0.00287860679448775\\ -76 0.0027731320103853\\ -77 0.00263952084059243\\ -78 0.00233720240014001\\ -79 0.00212697729708345\\ -80 0.00199788688851914\\ -81 0.00191003161317315\\ -82 0.00189242158319258\\ -83 0.00174836410231261\\ -84 0.0016140985772839\\ -85 0.00153318949743614\\ -86 0.00147245713810495\\ -87 0.00150762346036416\\ -88 0.00141698818640003\\ -89 0.00132612374554337\\ -90 0.0012294024439175\\ -91 0.00124357808130644\\ -92 0.00117351874267278\\ -93 0.00118805545766119\\ -94 0.00114318803197479\\ -95 0.00110915873689051\\ -96 0.00101959224657441\\ -97 0.00102415265668903\\ -98 0.000975302119471475\\ -99 0.000987996366258813\\ -100 0.00095742632740931\\ -101 0.000934265103971271\\ -102 0.000867894211150555\\ -103 0.000857812123034294\\ -104 0.000756099336226914\\ -105 0.000691162784782045\\ -106 0.000695473018163751\\ -107 0.000633925014757231\\ -108 0.000590287819831721\\ -109 0.000535337664369594\\ -110 0.000538252780320825\\ -111 0.000525826256051619\\ -112 0.000497879429315293\\ -113 0.000454616010876659\\ -114 0.000458838481347935\\ -115 0.000424839125494755\\ -116 0.000394335524998984\\ -117 0.000385945051019674\\ -118 0.000369692359961463\\ -119 0.000344267565220491\\ -120 0.000327893992619588\\ -121 0.000325828293258797\\ -122 0.000289237146884203\\ -123 0.000285070703523818\\ -124 0.000254192839215692\\ -125 0.000243444971741076\\ -126 0.000232080954193898\\ -127 0.000232619238643129\\ -128 0.000221784687572175\\ -129 0.000220230466851297\\ -130 0.000211029395043088\\ -131 0.000193736708718616\\ -132 0.000173392267570507\\ -133 0.000177242279228775\\ -134 0.000176692274847942\\ -135 0.000163578052215373\\ -136 0.000154604812916662\\ -137 0.000144545518887022\\ -138 0.00014456301163485\\ -139 0.000142459358743701\\ -140 0.0001308824450235\\ -141 0.000115457857714896\\ -142 0.00011695471803412\\ -143 0.000113103879953932\\ -144 0.000105257482462195\\ -145 0.000106888699118678\\ -146 0.000102374869876372\\ -147 9.76413644339824e-05\\ -148 8.4851357550636e-05\\ -149 7.68259382254786e-05\\ -150 7.45049326964518e-05\\ -151 7.30876640535304e-05\\ -152 7.5125116902824e-05\\ -153 7.84305708028061e-05\\ -154 8.45501238033231e-05\\ -155 8.17331778157062e-05\\ -156 7.55108766404751e-05\\ -157 7.48398655259728e-05\\ -158 7.14089505651308e-05\\ -159 6.99224620714426e-05\\ -160 7.04197347847714e-05\\ -161 6.53525683701177e-05\\ -162 6.08288863097357e-05\\ -163 5.89421929587584e-05\\ -164 5.72484731968775e-05\\ -165 5.18959286578684e-05\\ -166 5.16052819357681e-05\\ -167 5.18373975156787e-05\\ -168 5.1438260693119e-05\\ -169 4.74780120078228e-05\\ -170 4.76692285194498e-05\\ -171 4.63049226435535e-05\\ -172 4.3021574259984e-05\\ -173 4.34691080338747e-05\\ -174 4.07308380659776e-05\\ -175 3.92273577862465e-05\\ -176 3.78686988818524e-05\\ -177 3.75780813333193e-05\\ -178 3.89427911494855e-05\\ -179 3.55021611298483e-05\\ -180 3.21423511450592e-05\\ -181 3.1976175682275e-05\\ -182 3.15208931332488e-05\\ -183 3.02809343645595e-05\\ -184 2.93752192595988e-05\\ -185 2.83798341644248e-05\\ -186 2.84874468500327e-05\\ -187 2.6372748821241e-05\\ -188 2.83572106563901e-05\\ -189 2.73404573666424e-05\\ -190 2.64050982491902e-05\\ -191 2.68346097335579e-05\\ -192 2.66791537496404e-05\\ -193 2.44823861206895e-05\\ -194 2.44167662294569e-05\\ -195 2.42329915293711e-05\\ -196 2.57959950015393e-05\\ -197 2.5330918832603e-05\\ -198 2.43122483048608e-05\\ -199 2.29066675778814e-05\\ -200 2.28422313595056e-05\\ -201 2.2009103500163e-05\\ +1 1\\ +2 0.557596535615506\\ +3 0.0776565625508993\\ +4 0.0198777323861488\\ +5 0.00796998880225869\\ +6 0.00379999756962287\\ +7 0.0020374948726527\\ +8 0.00117378505868276\\ +9 0.000679072640424649\\ +10 0.000460205525227464\\ +11 0.000299176988582771\\ +12 0.000208705837325249\\ +13 0.000154155368421638\\ +14 0.000112123271112138\\ +15 8.51723029158506e-05\\ +16 6.84145894011081e-05\\ +17 5.22773173463893e-05\\ +18 4.19658990284636e-05\\ +19 3.35813556680048e-05\\ +20 2.69909819855066e-05\\ +21 2.16953761012216e-05\\ +22 1.75366202909715e-05\\ +23 1.58613744370523e-05\\ +24 1.07527796769547e-05\\ +25 1.21386374874815e-05\\ +26 1.02818699879239e-05\\ +27 8.94827138078525e-06\\ +28 7.32098660262362e-06\\ +29 6.29855369167726e-06\\ +30 5.75608316894847e-06\\ +31 5.01303204699423e-06\\ +32 4.32711144148808e-06\\ +33 3.75161010629661e-06\\ +34 3.13429056757287e-06\\ +35 2.97535388105766e-06\\ +36 2.59579082920428e-06\\ +37 2.35613630297938e-06\\ +38 2.13570046200869e-06\\ +39 1.91918226507556e-06\\ +40 1.79412386278695e-06\\ +41 1.60497187780807e-06\\ +42 1.50673230165654e-06\\ +43 1.29912483256279e-06\\ +44 1.18346887626955e-06\\ +45 1.20647020748118e-06\\ +46 9.3251240253618e-07\\ +47 1.0443397517104e-06\\ +48 7.81977076368339e-07\\ +49 7.06367946061837e-07\\ +50 6.4448540591806e-07\\ +51 5.8599197172186e-07\\ +52 5.19890971305916e-07\\ +53 4.77565688325081e-07\\ +54 4.45837257417211e-07\\ +55 4.52354751810163e-07\\ +56 4.29038013881601e-07\\ +57 4.19929738903634e-07\\ +58 3.83172910236919e-07\\ +59 3.69121286102551e-07\\ +60 3.7162945951419e-07\\ +61 3.46046785121521e-07\\ +62 3.2757610333777e-07\\ +63 3.15320978192865e-07\\ +64 2.96758711344714e-07\\ +65 2.96817636331721e-07\\ +66 2.71124713000919e-07\\ +67 2.62909013135003e-07\\ +68 2.43462036182122e-07\\ +69 2.28404291071355e-07\\ +70 2.24415864473535e-07\\ +71 2.28468445317059e-07\\ +72 2.16184931431001e-07\\ +73 2.03658004162582e-07\\ +74 1.83416377348938e-07\\ +75 1.77025922497321e-07\\ +76 1.70006600401098e-07\\ +77 1.63491825605367e-07\\ +78 1.57088823572816e-07\\ +79 1.6232157488019e-07\\ +80 1.55517018193347e-07\\ +81 1.44132341985942e-07\\ +82 1.33984748126276e-07\\ +83 1.25196828216566e-07\\ +84 1.21215390780978e-07\\ +85 1.15338855777099e-07\\ +86 1.08276150400162e-07\\ +87 1.04805124156927e-07\\ +88 1.01804326380992e-07\\ +89 1.0010166263465e-07\\ +90 9.48101967192271e-08\\ +91 8.89240656896784e-08\\ +92 8.81989115707906e-08\\ +93 8.26647202805855e-08\\ +94 7.78552030683533e-08\\ +95 8.16807953528358e-08\\ +96 7.66807903209987e-08\\ +97 7.25436978662745e-08\\ +98 6.71702595720625e-08\\ +99 7.01343166575355e-08\\ +100 6.80449977565014e-08\\ +101 6.68192875014292e-08\\ +102 6.22828481503855e-08\\ +103 6.14945471895087e-08\\ +104 6.70743224722449e-08\\ +105 6.88753801346713e-08\\ +106 6.74234915967769e-08\\ +107 6.38752997351369e-08\\ +108 6.10832293660283e-08\\ +109 5.60972305767953e-08\\ +110 5.08646226781961e-08\\ +111 4.95616876939365e-08\\ +112 4.46580704289759e-08\\ +113 4.2687817302539e-08\\ +114 4.10291135987935e-08\\ +115 3.78768817309767e-08\\ +116 3.73532344958446e-08\\ +117 3.51670172127996e-08\\ +118 3.22365221541983e-08\\ +119 3.0914108557954e-08\\ +120 3.02505081670252e-08\\ +121 2.87114558023773e-08\\ +122 2.66853087438055e-08\\ +123 2.60049357892773e-08\\ +124 2.51289922847544e-08\\ +125 2.39005712471205e-08\\ +126 2.24094810183141e-08\\ +127 2.01607019229326e-08\\ +128 1.92164353138885e-08\\ +129 1.82687256142865e-08\\ +130 1.78764123154733e-08\\ +131 1.90420480961642e-08\\ +132 1.79866998160551e-08\\ +133 1.62987103158868e-08\\ +134 1.64054691550642e-08\\ +135 1.53093691487969e-08\\ +136 1.48209887890633e-08\\ +137 1.48119659868498e-08\\ +138 1.39479300275188e-08\\ +139 1.30269784768553e-08\\ +140 1.39128297693094e-08\\ +141 1.35601176412506e-08\\ +142 1.31199716923384e-08\\ +143 1.16605449828788e-08\\ +144 1.1256315915413e-08\\ +145 1.10059046096852e-08\\ +146 1.05336940404569e-08\\ +147 9.87947792918718e-09\\ +148 9.87385221564589e-09\\ +149 9.86634931285779e-09\\ +150 9.87265433744658e-09\\ +151 8.57424740655266e-09\\ +152 7.71781430780305e-09\\ }; \addplot [color=mycolor2, forget plot] table[row sep=crcr]{% -1 118.544185368288\\ -2 0.687337944312083\\ -3 0.0352052378196977\\ -4 0.0177116530504275\\ -5 0.00921115423330148\\ -6 0.00742079353381528\\ -7 0.00435978025623298\\ -8 0.00269077966124769\\ -9 0.00326151218772164\\ -10 0.00207394458931424\\ -11 0.00295893040643868\\ -12 0.00133204038205353\\ -13 0.00183522720801244\\ -14 0.00145996981575274\\ -15 0.000990797743084482\\ -16 0.00126621063572918\\ -17 0.000551496956428083\\ -18 0.000515017395303755\\ -19 0.000232316315215174\\ -20 0.000223232804362511\\ -21 0.000208790455128303\\ -22 0.000151172313680646\\ -23 0.000161588020568097\\ -24 0.000122705865740307\\ -25 0.000118687703674217\\ -26 0.000113287133295243\\ +1 1\\ +2 0.00579815823253321\\ +3 0.000296979878939853\\ +4 0.000149409715840567\\ +5 7.77022863220549e-05\\ +6 6.2599388664747e-05\\ +7 3.67776811885644e-05\\ +8 2.26985377046381e-05\\ +9 2.75130507463434e-05\\ +10 1.74951186586756e-05\\ +11 2.49605697423791e-05\\ +12 1.12366572676273e-05\\ +13 1.54813768580115e-05\\ +14 1.23158281548519e-05\\ +15 8.35804590504644e-06\\ +16 1.06813390449761e-05\\ +17 4.65224806020403e-06\\ +18 4.34451840639611e-06\\ +19 1.95974449943222e-06\\ +20 1.88311897094728e-06\\ +21 1.76128803348428e-06\\ +22 1.275240225499e-06\\ +23 1.36310372428718e-06\\ +24 1.03510657531701e-06\\ +25 1.00121067351798e-06\\ +26 9.55653227048963e-07\\ }; \end{axis} \end{tikzpicture}% \ No newline at end of file diff --git a/mp5/img_rcg-1.png b/mp5/img_rcg-1.png index a5b68fd..35a0f84 100644 Binary files a/mp5/img_rcg-1.png and b/mp5/img_rcg-1.png differ diff --git a/mp5/img_rcg.tex b/mp5/img_rcg.tex index 85b5ede..94c849a 100644 --- a/mp5/img_rcg.tex +++ b/mp5/img_rcg.tex @@ -3,12 +3,12 @@ \begin{tikzpicture} \begin{axis}[% -width=2.604in, -height=2.604in, -at={(0.962in,0.66in)}, +width=1.736in, +height=1.736in, +at={(0.75in,0.653in)}, scale only axis, -point meta min=-0.0235938830592846, -point meta max=0.993939917509565, +point meta min=-0.0128981437433609, +point meta max=0.983438542582913, axis on top, xmin=0.5, xmax=250.5, diff --git a/mp5/res_log.tex b/mp5/res_log.tex index b46fa1b..41f9e2d 100644 --- a/mp5/res_log.tex +++ b/mp5/res_log.tex @@ -6,15 +6,15 @@ \begin{tikzpicture} \begin{axis}[% -width=4.521in, -height=3.566in, -at={(0.758in,0.481in)}, +width=6.028in, +height=4.754in, +at={(1.011in,0.642in)}, scale only axis, xmin=0, -xmax=250, +xmax=160, ymode=log, -ymin=1e-06, -ymax=118.544185368288, +ymin=1e-09, +ymax=1, yminorticks=true, axis background/.style={fill=white}, title style={font=\bfseries}, @@ -22,236 +22,187 @@ title={Residual norms over iteration (y is log)} ] \addplot [color=mycolor1, forget plot] table[row sep=crcr]{% -1 11.1886780605947\\ -2 37.688016651504\\ -3 20.5612191459901\\ -4 10.1126098486914\\ -5 6.60295496425051\\ -6 4.03660691189902\\ -7 2.83696084221764\\ -8 2.02776259797851\\ -9 1.55623996384361\\ -10 1.26373343652138\\ -11 0.968080560648255\\ -12 0.808014633949812\\ -13 0.666010141343337\\ -14 0.547935882226161\\ -15 0.455778233603678\\ -16 0.379140935249178\\ -17 0.325385471261729\\ -18 0.27468130551472\\ -19 0.241318492631058\\ -20 0.203089896725741\\ -21 0.186402910640978\\ -22 0.164285195168731\\ -23 0.161498606695505\\ -24 0.13554465057248\\ -25 0.108256329622331\\ -26 0.096685060868788\\ -27 0.0901032819163635\\ -28 0.0846560423043152\\ -29 0.0720452421305809\\ -30 0.0638081561682231\\ -31 0.0586938472203618\\ -32 0.0524931123015411\\ -33 0.0500311324114693\\ -34 0.0441377829780386\\ -35 0.0430660141184194\\ -36 0.0399007401033348\\ -37 0.0357338939687774\\ -38 0.0334582949737502\\ -39 0.0301143079798452\\ -40 0.0277450095617363\\ -41 0.0246457162237748\\ -42 0.0226451015644127\\ -43 0.0214585559097522\\ -44 0.0197750620834928\\ -45 0.0184656387170807\\ -46 0.01424126114094\\ -47 0.014925314691824\\ -48 0.0138939777471367\\ -49 0.0124854196790047\\ -50 0.0121443393806818\\ -51 0.0111248201462271\\ -52 0.0117048407914948\\ -53 0.0113323279545552\\ -54 0.0106729739632728\\ -55 0.00944425804326873\\ -56 0.00871066933482761\\ -57 0.00881210600673333\\ -58 0.00795156096327054\\ -59 0.00726613636631538\\ -60 0.00700470814378812\\ -61 0.00621935380743637\\ -62 0.00580029695866593\\ -63 0.00557779006381989\\ -64 0.00510498070037545\\ -65 0.0047009622833176\\ -66 0.00438789538442034\\ -67 0.00435585711109691\\ -68 0.00396640165581629\\ -69 0.00393477324085466\\ -70 0.00371617606606877\\ -71 0.00327886320159879\\ -72 0.00312237367539779\\ -73 0.00299349690875735\\ -74 0.00284074963413765\\ -75 0.00287860679448775\\ -76 0.0027731320103853\\ -77 0.00263952084059243\\ -78 0.00233720240014001\\ -79 0.00212697729708345\\ -80 0.00199788688851914\\ -81 0.00191003161317315\\ -82 0.00189242158319258\\ -83 0.00174836410231261\\ -84 0.0016140985772839\\ -85 0.00153318949743614\\ -86 0.00147245713810495\\ -87 0.00150762346036416\\ -88 0.00141698818640003\\ -89 0.00132612374554337\\ -90 0.0012294024439175\\ -91 0.00124357808130644\\ -92 0.00117351874267278\\ -93 0.00118805545766119\\ -94 0.00114318803197479\\ -95 0.00110915873689051\\ -96 0.00101959224657441\\ -97 0.00102415265668903\\ -98 0.000975302119471475\\ -99 0.000987996366258813\\ -100 0.00095742632740931\\ -101 0.000934265103971271\\ -102 0.000867894211150555\\ -103 0.000857812123034294\\ -104 0.000756099336226914\\ -105 0.000691162784782045\\ -106 0.000695473018163751\\ -107 0.000633925014757231\\ -108 0.000590287819831721\\ -109 0.000535337664369594\\ -110 0.000538252780320825\\ -111 0.000525826256051619\\ -112 0.000497879429315293\\ -113 0.000454616010876659\\ -114 0.000458838481347935\\ -115 0.000424839125494755\\ -116 0.000394335524998984\\ -117 0.000385945051019674\\ -118 0.000369692359961463\\ -119 0.000344267565220491\\ -120 0.000327893992619588\\ -121 0.000325828293258797\\ -122 0.000289237146884203\\ -123 0.000285070703523818\\ -124 0.000254192839215692\\ -125 0.000243444971741076\\ -126 0.000232080954193898\\ -127 0.000232619238643129\\ -128 0.000221784687572175\\ -129 0.000220230466851297\\ -130 0.000211029395043088\\ -131 0.000193736708718616\\ -132 0.000173392267570507\\ -133 0.000177242279228775\\ -134 0.000176692274847942\\ -135 0.000163578052215373\\ -136 0.000154604812916662\\ -137 0.000144545518887022\\ -138 0.00014456301163485\\ -139 0.000142459358743701\\ -140 0.0001308824450235\\ -141 0.000115457857714896\\ -142 0.00011695471803412\\ -143 0.000113103879953932\\ -144 0.000105257482462195\\ -145 0.000106888699118678\\ -146 0.000102374869876372\\ -147 9.76413644339824e-05\\ -148 8.4851357550636e-05\\ -149 7.68259382254786e-05\\ -150 7.45049326964518e-05\\ -151 7.30876640535304e-05\\ -152 7.5125116902824e-05\\ -153 7.84305708028061e-05\\ -154 8.45501238033231e-05\\ -155 8.17331778157062e-05\\ -156 7.55108766404751e-05\\ -157 7.48398655259728e-05\\ -158 7.14089505651308e-05\\ -159 6.99224620714426e-05\\ -160 7.04197347847714e-05\\ -161 6.53525683701177e-05\\ -162 6.08288863097357e-05\\ -163 5.89421929587584e-05\\ -164 5.72484731968775e-05\\ -165 5.18959286578684e-05\\ -166 5.16052819357681e-05\\ -167 5.18373975156787e-05\\ -168 5.1438260693119e-05\\ -169 4.74780120078228e-05\\ -170 4.76692285194498e-05\\ -171 4.63049226435535e-05\\ -172 4.3021574259984e-05\\ -173 4.34691080338747e-05\\ -174 4.07308380659776e-05\\ -175 3.92273577862465e-05\\ -176 3.78686988818524e-05\\ -177 3.75780813333193e-05\\ -178 3.89427911494855e-05\\ -179 3.55021611298483e-05\\ -180 3.21423511450592e-05\\ -181 3.1976175682275e-05\\ -182 3.15208931332488e-05\\ -183 3.02809343645595e-05\\ -184 2.93752192595988e-05\\ -185 2.83798341644248e-05\\ -186 2.84874468500327e-05\\ -187 2.6372748821241e-05\\ -188 2.83572106563901e-05\\ -189 2.73404573666424e-05\\ -190 2.64050982491902e-05\\ -191 2.68346097335579e-05\\ -192 2.66791537496404e-05\\ -193 2.44823861206895e-05\\ -194 2.44167662294569e-05\\ -195 2.42329915293711e-05\\ -196 2.57959950015393e-05\\ -197 2.5330918832603e-05\\ -198 2.43122483048608e-05\\ -199 2.29066675778814e-05\\ -200 2.28422313595056e-05\\ -201 2.2009103500163e-05\\ +1 1\\ +2 0.557596535615506\\ +3 0.0776565625508993\\ +4 0.0198777323861488\\ +5 0.00796998880225869\\ +6 0.00379999756962287\\ +7 0.0020374948726527\\ +8 0.00117378505868276\\ +9 0.000679072640424649\\ +10 0.000460205525227464\\ +11 0.000299176988582771\\ +12 0.000208705837325249\\ +13 0.000154155368421638\\ +14 0.000112123271112138\\ +15 8.51723029158506e-05\\ +16 6.84145894011081e-05\\ +17 5.22773173463893e-05\\ +18 4.19658990284636e-05\\ +19 3.35813556680048e-05\\ +20 2.69909819855066e-05\\ +21 2.16953761012216e-05\\ +22 1.75366202909715e-05\\ +23 1.58613744370523e-05\\ +24 1.07527796769547e-05\\ +25 1.21386374874815e-05\\ +26 1.02818699879239e-05\\ +27 8.94827138078525e-06\\ +28 7.32098660262362e-06\\ +29 6.29855369167726e-06\\ +30 5.75608316894847e-06\\ +31 5.01303204699423e-06\\ +32 4.32711144148808e-06\\ +33 3.75161010629661e-06\\ +34 3.13429056757287e-06\\ +35 2.97535388105766e-06\\ +36 2.59579082920428e-06\\ +37 2.35613630297938e-06\\ +38 2.13570046200869e-06\\ +39 1.91918226507556e-06\\ +40 1.79412386278695e-06\\ +41 1.60497187780807e-06\\ +42 1.50673230165654e-06\\ +43 1.29912483256279e-06\\ +44 1.18346887626955e-06\\ +45 1.20647020748118e-06\\ +46 9.3251240253618e-07\\ +47 1.0443397517104e-06\\ +48 7.81977076368339e-07\\ +49 7.06367946061837e-07\\ +50 6.4448540591806e-07\\ +51 5.8599197172186e-07\\ +52 5.19890971305916e-07\\ +53 4.77565688325081e-07\\ +54 4.45837257417211e-07\\ +55 4.52354751810163e-07\\ +56 4.29038013881601e-07\\ +57 4.19929738903634e-07\\ +58 3.83172910236919e-07\\ +59 3.69121286102551e-07\\ +60 3.7162945951419e-07\\ +61 3.46046785121521e-07\\ +62 3.2757610333777e-07\\ +63 3.15320978192865e-07\\ +64 2.96758711344714e-07\\ +65 2.96817636331721e-07\\ +66 2.71124713000919e-07\\ +67 2.62909013135003e-07\\ +68 2.43462036182122e-07\\ +69 2.28404291071355e-07\\ +70 2.24415864473535e-07\\ +71 2.28468445317059e-07\\ +72 2.16184931431001e-07\\ +73 2.03658004162582e-07\\ +74 1.83416377348938e-07\\ +75 1.77025922497321e-07\\ +76 1.70006600401098e-07\\ +77 1.63491825605367e-07\\ +78 1.57088823572816e-07\\ +79 1.6232157488019e-07\\ +80 1.55517018193347e-07\\ +81 1.44132341985942e-07\\ +82 1.33984748126276e-07\\ +83 1.25196828216566e-07\\ +84 1.21215390780978e-07\\ +85 1.15338855777099e-07\\ +86 1.08276150400162e-07\\ +87 1.04805124156927e-07\\ +88 1.01804326380992e-07\\ +89 1.0010166263465e-07\\ +90 9.48101967192271e-08\\ +91 8.89240656896784e-08\\ +92 8.81989115707906e-08\\ +93 8.26647202805855e-08\\ +94 7.78552030683533e-08\\ +95 8.16807953528358e-08\\ +96 7.66807903209987e-08\\ +97 7.25436978662745e-08\\ +98 6.71702595720625e-08\\ +99 7.01343166575355e-08\\ +100 6.80449977565014e-08\\ +101 6.68192875014292e-08\\ +102 6.22828481503855e-08\\ +103 6.14945471895087e-08\\ +104 6.70743224722449e-08\\ +105 6.88753801346713e-08\\ +106 6.74234915967769e-08\\ +107 6.38752997351369e-08\\ +108 6.10832293660283e-08\\ +109 5.60972305767953e-08\\ +110 5.08646226781961e-08\\ +111 4.95616876939365e-08\\ +112 4.46580704289759e-08\\ +113 4.2687817302539e-08\\ +114 4.10291135987935e-08\\ +115 3.78768817309767e-08\\ +116 3.73532344958446e-08\\ +117 3.51670172127996e-08\\ +118 3.22365221541983e-08\\ +119 3.0914108557954e-08\\ +120 3.02505081670252e-08\\ +121 2.87114558023773e-08\\ +122 2.66853087438055e-08\\ +123 2.60049357892773e-08\\ +124 2.51289922847544e-08\\ +125 2.39005712471205e-08\\ +126 2.24094810183141e-08\\ +127 2.01607019229326e-08\\ +128 1.92164353138885e-08\\ +129 1.82687256142865e-08\\ +130 1.78764123154733e-08\\ +131 1.90420480961642e-08\\ +132 1.79866998160551e-08\\ +133 1.62987103158868e-08\\ +134 1.64054691550642e-08\\ +135 1.53093691487969e-08\\ +136 1.48209887890633e-08\\ +137 1.48119659868498e-08\\ +138 1.39479300275188e-08\\ +139 1.30269784768553e-08\\ +140 1.39128297693094e-08\\ +141 1.35601176412506e-08\\ +142 1.31199716923384e-08\\ +143 1.16605449828788e-08\\ +144 1.1256315915413e-08\\ +145 1.10059046096852e-08\\ +146 1.05336940404569e-08\\ +147 9.87947792918718e-09\\ +148 9.87385221564589e-09\\ +149 9.86634931285779e-09\\ +150 9.87265433744658e-09\\ +151 8.57424740655266e-09\\ +152 7.71781430780305e-09\\ }; \addplot [color=mycolor2, forget plot] table[row sep=crcr]{% -1 118.544185368288\\ -2 0.687337944312083\\ -3 0.0352052378196977\\ -4 0.0177116530504275\\ -5 0.00921115423330148\\ -6 0.00742079353381528\\ -7 0.00435978025623298\\ -8 0.00269077966124769\\ -9 0.00326151218772164\\ -10 0.00207394458931424\\ -11 0.00295893040643868\\ -12 0.00133204038205353\\ -13 0.00183522720801244\\ -14 0.00145996981575274\\ -15 0.000990797743084482\\ -16 0.00126621063572918\\ -17 0.000551496956428083\\ -18 0.000515017395303755\\ -19 0.000232316315215174\\ -20 0.000223232804362511\\ -21 0.000208790455128303\\ -22 0.000151172313680646\\ -23 0.000161588020568097\\ -24 0.000122705865740307\\ -25 0.000118687703674217\\ -26 0.000113287133295243\\ +1 1\\ +2 0.00579815823253321\\ +3 0.000296979878939853\\ +4 0.000149409715840567\\ +5 7.77022863220549e-05\\ +6 6.2599388664747e-05\\ +7 3.67776811885644e-05\\ +8 2.26985377046381e-05\\ +9 2.75130507463434e-05\\ +10 1.74951186586756e-05\\ +11 2.49605697423791e-05\\ +12 1.12366572676273e-05\\ +13 1.54813768580115e-05\\ +14 1.23158281548519e-05\\ +15 8.35804590504644e-06\\ +16 1.06813390449761e-05\\ +17 4.65224806020403e-06\\ +18 4.34451840639611e-06\\ +19 1.95974449943222e-06\\ +20 1.88311897094728e-06\\ +21 1.76128803348428e-06\\ +22 1.275240225499e-06\\ +23 1.36310372428718e-06\\ +24 1.03510657531701e-06\\ +25 1.00121067351798e-06\\ +26 9.55653227048963e-07\\ }; \end{axis} \end{tikzpicture}% \ No newline at end of file