x1=[0.5404 0.2567 0.1914 0.1988];
y1=.01*[0.0314 0.0431 0.0399 0.0421];
x2=[0.5836 0.2570 0.1947 0.2060];
y2=.01*[0.0370 0.0417 0.0463 0.0578];
x3=[0.1560 0.1134 0.0852 0.0577];
y3=.01*[0.0553 0.1196 0.1725 0.1731];
x4=[0.0758 0.0447 0.0417 0.0479];
y4=.01*[0.0093 0.0153 0.0155 0.0178];
sel=(1:4); %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.2187 0.2349 0.2651 0.3060];
y1=.01*[0.0531 0.0521 0.0535 0.0521];
x2=[0.2329 0.2556 0.2899 0.3377];
y2=.01*[0.0703 0.0675 0.0655 0.0619];
x3=[0.0601 0.0603 0.0811 0.0980];
y3=.01*[0.1784 0.1571 0.1469 0.1310];
x4=[0.0212 0.0246 0.0372 0.0486];
y4=.01*[0.0083 0.0080 0.0077 0.0087];
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:4);
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],'New method','Huber'); xlabel('% bias');
ylabel('standard deviation');
figure(2),q2=plot(x2,y2,'o'); plot(x2,y2); hold off;
legend([p2,q2],'New method','Huber'); xlabel('% bias');
ylabel('standard deviation');
figure(3),q3=plot(abs(x3),y3,'o'); plot(abs(x3),y3); hold off;
legend([p3,q3],'New method','Huber'); xlabel('% bias');
ylabel('standard deviation');
figure(4),q4=plot(abs(x4),y4,'o'); plot(abs(x4),y4); hold off;
legend([p4,q4],'New method','Huber'); xlabel('% bias');
ylabel('standard deviation');

