x1=[0.1888 0.1835 0.1839 0.1293 0.2312 0.3317]*100;
y1=.01*[0.0640  0.0650 0.07 0.0778 0.0618 0.0585];
x2=[0.2149 0.2114 0.2161 0.1554 0.2657 0.3857]*100;
y2=.01*[0.0767 0.0813 0.0886 0.0991 0.0728 0.0658];
x3=[-0.0811 -0.0826 -0.0843 -0.0755 -0.0876 -0.1028]*100;
y3=.01*[0.2384 0.2774 0.3017 0.3292 0.2227 0.1917];
x4=[-0.0277 -0.0272 -0.0274 -0.0154 -0.0401 -0.0670]*100;
y4=.01*[0.0223 .0215 0.0215 0.0232 0.0217 0.0206];
sel=[1 5 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),p1=plot(x1,y1,'*'); hold on; plot(x1,y1);
%for i=1:length(x1),
%  plot([x1(i)/100*t1-y1(i)/sqrt(50),x1(i)/100*t1+y1(i)/sqrt(50)]/t1*100,...
%		[y1(i),y1(i)],'r'); end
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.9172 0.4916 0.3064 0.1822 0.0999 0.1038]*100;
y1=.01*[0.0369 0.0627 0.0707 0.0754 0.0789 0.0924];
x2=[1.0311 0.5449 0.3319 0.1935 0.1062 0.1034]*100;
y2=.01*[0.0428 0.0709 0.0804 0.0855 0.0944 0.1094];
x3=[-0.1743 -0.1537 -.1303 -0.1014 -0.0698 -0.0349]*100;
y3=.01*[0.0443 0.1134 0.1491 0.1932 0.2794 0.3360];
x4=[-0.1292 -0.0574 -.0301 -0.0134 -0.0041 0.0025]*100;
y4=.01*[0.0144 0.0238 0.0265 0.0278 0.0299 0.0329];
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'); axis([0 110 4e-4 12e-4])
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'); axis([0 18 0 3.5e-3]);
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');

