%%
eps =   0.202;
C = eye(2);
[m1,c1,a1,y1,b1,v1] =IptoyRegFunc(eps,C,1);

cols = {'r+','g+','b+','m+','c+','y+','ro','go','bo','mo','co','yo'};

 figure(4); clf;    hold on; 
 loglog(b1,v1,'k');xlabel('bias'); ylabel('variance');
 for j = 1:length(b1)
        loglog(b1(j),v1(j),cols{j});
 end

 %%
%C2 = [1 -1 ; -1 1];
theta = 2*pi/8;
lmd = 0.1;
R = [cos(theta), - sin(theta); sin(theta), cos(theta)];
C2 = R'* diag([1,lmd])*R;
[m2,c2,a2,y2,b2,v2] =IptoyRegFunc(eps,C2,2);


 figure(4);   hold on; 
 loglog(b2,v2,'--k');xlabel('bias'); ylabel('variance');
 for j = 1:length(b2)
        loglog(b2(j),v2(j),cols{j});
 end