
nmin = 0;
bmin = 0;
nmax = 0;
bmax = 0;
np = size(par,1);
for i=1:np
 if par(i,1)<0.001
  bmin = bmin+par(i,2);
  nmin = nmin+1;
 end;
 if par(i,1)>0.999
  bmax = bmax+par(i,2);
  nmax = nmax+1;
 end;
end;
bmin = bmin/nmin;
bmax = bmax/nmax;
plot(0,bmin,'bo');
plot(1,bmax,'bo');