x1=[0.2038 0.1292 0.0546 0.0320 0.3674 0.2548]*100;
y1=.01*[0.0413 0.0418 0.0415 0.0431 0.0306 0.0331];
x2=[0.2180 0.1348 0.0531 0.0296 0.4231 0.2917]*100;
y2=.01*[0.0432 0.0431 0.0478 0.0481 0.0318 0.0345];
x3=[-0.0980 -0.0811 -0.0593 -0.0573 -0.10078 -0.0851]*100;
y3=.01*[0.1310 0.1475 0.1597 0.1842 0.1138 0.1319];
x4=[-0.0337 -0.0195 -0.0046 0.0004 -0.08206 -0.0550]*100;
y4=.01*[0.0129 0.0130 0.0132 0.0136 0.0109 0.0115];
sel=[4 2 5];
x1=x1(sel); y1=y1(sel);
x2=x2(sel); y2=y2(sel);
x3=x3(sel); y3=y3(sel);
x4=x4(sel); y4=y4(sel);
x5=[0.0150 0.0144 0.0145 0.0152];
y5=[0.0477 0.0473 0.0470 0.0471];

figure(1),p1=plot(x1,y1,'*'); hold on; plot(x1,y1);
figure(2),p2=plot(x2,y2,'*'); hold on; plot(x2,y2);
figure(3),p3=plot(abs(x3),y3,'*'); hold on; plot(abs(x3),y3);
figure(4),p4=plot(abs(x4),y4,'*'); hold on; plot(abs(x4),y4);

x1=[0.1025 0.0269 0.0523 0.1106 0.2088 0.5289 1.0783]*100;
y1=.01*[0.0826 0.0600 0.0567 0.0541 0.0498 0.0347 0.0147];
x2=[0.0976 0.0202 0.0481 0.1119 0.2190 0.5836 1.1933]*100;
y2=.01*[0.0673 0.0663 0.0612 0.0560 0.0513 0.0371 0.0156];
x3=[-0.0098 -0.0286 -0.0535 -0.0851 -0.1134  -0.1560 -0.2049]*100;
y3=.01*[0.2450 0.2405 0.2224  0.1739 0.1196 0.0553 0.0126];
x4=[0.0130 0.0075 0.0042 -0.0020 -0.0142 -0.0620 -0.1665]*100;
y4=.01*[0.0200 0.0196 0.0197 0.0195 0.0184 0.0134 0.0063];
x5=[0.0229 0.0085 0.0042 0.0070 0.0132 0.0124];
y5=[1.8208 0.5879 0.3173 0.2200 0.1525 0.0132];
sel=(1:6);
x1=x1(sel); y1=y1(sel); x2=x2(sel); y2=y2(sel);
x3=x3(sel); y3=y3(sel); x4=x4(sel); y4=y4(sel);

figure(1),q1=plot(x1,y1,'o'); plot(x1,y1); hold off;
legend([p1,q1],'Proposed penalty','Huber penalty'); xlabel('% bias');
ylabel('standard deviation');
figure(2),q2=plot(x2,y2,'o'); plot(x2,y2); hold off;
legend([p2,q2],'Proposed penalty','Huber penalty'); xlabel('% bias');
ylabel('standard deviation');
figure(3),q3=plot(abs(x3),y3,'o'); plot(abs(x3),y3); hold off;
legend([p3,q3],'Proposed penalty','Huber penalty'); xlabel('% bias');
ylabel('standard deviation');
figure(4),q4=plot(abs(x4),y4,'o'); plot(abs(x4),y4); hold off;
legend([p4,q4],'Proposed penalty','Huber penalty'); xlabel('% bias');
ylabel('standard deviation');

