getf("IDbin.sci");
//
// Representation d'intervalles de dispersion optimaux 
// pour la loi binomiale B(10,p) en fonction de p
//
n=10; p=linspace(0,1,100); alpha=0.1;
inter=[];
for pr=p;
  inter=[inter;IDoptim_bin(n,pr,alpha)];
end;
p=p'*[1,1];
xset("window",0);xbasc();
xset("font",2,3);
plotframe([0,0,1,n],[0,10,0,10],[%f,%f],..
["Intervalles de dispersion loi binomiale","",""]);
plot2d2("gnn",p,inter,[1,1],"000");
//
// Representation d'intervalles de dispersion symetriques 
// pour la loi binomiale B(10,p) en fonction de p
//
n=10; p=linspace(0,1,100); alpha=0.1;
inter=[];
for pr=p;
  inter=[inter;IDsym_bin(n,pr,alpha)];
end;
p=p'*[1,1];
plot2d2("gnn",p,inter,[2,2],"000");

