mp5: corrected ex44 matlab code

This commit is contained in:
Claudio Maggioni 2020-11-30 20:26:06 +01:00
parent 5bbd92b438
commit 0b1939a9d2
6 changed files with 377 additions and 478 deletions

View file

@ -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');

BIN
mp5/img_my-1.png Normal file

Binary file not shown.

After

Width:  |  Height:  |  Size: 56 KiB

View file

@ -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}%

Binary file not shown.

Before

Width:  |  Height:  |  Size: 58 KiB

After

Width:  |  Height:  |  Size: 58 KiB

View file

@ -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,

View file

@ -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}%