//Conditions de bonne utilisation : commencer par réduire la console scilab de façon à ce que elle ne chevauche pas la fenetre du dessous (que l'on ne met amais au premier plan, pour des raisons de manque de surface sur l'écran).

clear;
stacksize(150000000);
driver("X11");
theta0=%pi*0.1;
theta1=%pi*0.4;
theta0u=%pi*0.23;
theta1u=%pi*0.28;	
Rcondu=0.95;
ncolhisto=20;
n_sample_histo=10000;


%symbol=1; 
%times=2; 
%times_italic=3; 
%times_bold=4; 
%helvetica= 6; 
%helvetica_italic=7; 
%helvetica_bold=8; 
%helvetica_bold_italique=9; 
%8pts=0;%10pts=1;%12pts=2; 
%14pts=3;%18pts=4;%24pts=5; 

deff('y=changecouleur(t)','y=modulo(t+1,35)+1');
deff('y=densite_r_gauss(x)','y=x.*exp(-x.*x./2)');

function y=densite_r_unif(r)
	un=ones(r);
	g=min(1,un./r);
	y=r.*max(0,asin(g)-acos(g));
endfunction

function y=densite_t_unif(theta)
	c=cos(theta);
	s=sin(theta);
	y=ones(theta)./(8*max(c.*c,s.*s));
endfunction
	

securite=0.95;
couleur=1;
npoints_cercles=1000;
npoints_densites=1000;
coeff_duree=1; //pour une valeur 1, les parametres sont choisis de façon à ce que l'animation se déroule en un temps adapté à l'amphi
nbre_histo_fluct=floor(10*coeff_duree);

Rmax=1.5*sqrt(-2*log(1-securite));// la proba qu'un couple de gaussiennes standard definisse un pt à distance <Rmax de l'origine est securite
Rmaxu=sqrt(2);
Rmaxcarre=Rmax^2;
vecteur_nombres=[1,20,50,100,1000,10,10,10,10,10,10,10,10,10];
paquet_monochrome=25;
abs_densite_r_gauss=linspace(0,1.5*Rmax,npoints_densites);
ord_densite_r_gauss=densite_r_gauss(abs_densite_r_gauss);
deuxpi=2*%pi;
abs_densite_t_gauss=linspace(-%pi,%pi,npoints_densites);
ord_densite_t_gauss=ones(1,npoints_densites)./(deuxpi);

abs_densite_r_unif=linspace(0.01,Rmaxu,npoints_densites);
ord_densite_r_unif=densite_r_unif(abs_densite_r_unif);
abs_densite_t_unif=linspace(-%pi,%pi,npoints_densites);
ord_densite_t_unif=densite_t_unif(abs_densite_t_unif);

theta_cercles=linspace(0,deuxpi,1000);

function affiche_cercles(r1,r2)
	abs_c=cos(theta_cercles);
	ord_c=sin(theta_cercles);
	abs_pol=[r1*abs_c,r2*abs_c];
	ord_pol=[r1*ord_c,r2*ord_c];	
	xpoly(abs_pol,ord_pol);
endfunction

function affiche_cercle(r1)
	abs_pol=r1*cos(theta_cercles);
	ord_pol=r1*sin(theta_cercles);
	xpoly(abs_pol,ord_pol);
endfunction


function cadre_haut_gauss() 
	scf(1); 
	clf(1);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%24pts;)
	a.title.font_style = 6; 
	a.title.text="Points aleatoires independants de loi gaussienne";
	a.thickness=2; 
	xselect();
	show_pixmap();
endfunction 

function cbgrg() //cadre bas gauche rayon gaussien 
	scf(2); 
	clf(2);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Histogramme des rayons";
	a.thickness=2; 
	plot2d(abs_densite_r_gauss,ord_densite_r_gauss,style=5, rect=[0,0,Rmax,1], frameflag=1);
endfunction 

function cbdtg() //cadre bas droite theta gaussien 
	scf(3); 
	clf(3);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Histogramme des angles";
	a.thickness=2; 
	plot2d(abs_densite_t_gauss,ord_densite_t_gauss,style=5, rect=[-%pi,0,%pi,0.4], frameflag=1);
endfunction 

function cbgtg() //cadre bas gauche theta gaussien 
	scf(2); 
	clf(2);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Histogramme des angles";
	a.thickness=2; 
	plot2d(abs_densite_t_gauss,ord_densite_t_gauss,style=5, rect=[-%pi,0,%pi,0.4], frameflag=1);
endfunction 

function cbdrg() //cadre bas droite rayon gaussien 
	scf(3); 
	clf(3);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Rayons conditionnes (a<theta<b)";
	a.thickness=2; 
	plot2d(abs_densite_r_gauss,ord_densite_r_gauss,style=5, rect=[0,0,Rmax,1], frameflag=1);
endfunction 

function cbdtcg() //cadre bas droite theta conditionne gaussien 
	scf(3); 
	clf(3);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Angles conditionnes (A<R<B)";
	a.thickness=2; 
	plot2d(abs_densite_t_gauss,ord_densite_t_gauss,style=5, rect=[-%pi,0,%pi,0.4], frameflag=1);
endfunction 

function cadre_haut_unif() 
	scf(1); 
	clf(1);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Points aleatoires independants uniformes sur le carre";
	a.thickness=2; 
	xselect();
	show_pixmap();
endfunction 


function cbgru() //cadre bas gauche rayon unif
	scf(2); 
	clf(2);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Histogramme des rayons";
	a.thickness=2; 
	plot2d(abs_densite_r_unif,ord_densite_r_unif,style=5, rect=[0,0,Rmaxu,2], frameflag=1);
endfunction 

function cbdtu() //cadre bas droite theta unif
	scf(3); 
	clf(3);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Histogramme des angles";
	a.thickness=2; 
	plot2d(abs_densite_t_unif,ord_densite_t_unif,style=5, rect=[-%pi,0,%pi,0.35], frameflag=1);
endfunction 

function cbgtu() //cadre bas gauche theta unif
	scf(2); 
	clf(2);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Histogramme des angles";
	a.thickness=2; 
	plot2d(abs_densite_t_unif,ord_densite_t_unif,style=5, rect=[-%pi,0,%pi,0.35], frameflag=1);
endfunction 

function cbdru() //cadre bas droite rayon unif 
	scf(3); 
	clf(3);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Rayons conditionnes (a<theta<b)";
	a.thickness=2; 
	plot2d(abs_densite_r_unif,ord_densite_r_unif,style=5, rect=[0,0,Rmaxu,2], frameflag=1);
endfunction 

function cbdtcu() //cadre bas droite theta conditionne unif
	scf(3); 
	clf(3);
	a = gca(); // a est le handle sur ces axes par defaut 
	a.title.font_size = 5; //(=%18pts;)
	a.title.font_style = 6; //(= %times_bold;)
	a.title.text="Angles conditionnes (R>A)";
	a.thickness=2; 
	plot2d(abs_densite_t_unif,ord_densite_t_unif,style=5, rect=[-%pi,0,%pi,0.35], frameflag=1);
endfunction 


function trois_figures()
	//scf(1);
	//xdel(1);
	scf(1);
	clf(1);
	xset("wpos",205,0)
	xselect();
	f=gcf();
	f.pixmap="on";
	xset("wdim",600,300);
	
	//scf(2); 
	xdel(2);
	scf(2);
	clf(2);
	xset("wpos",10,310)
	f=gcf();
	f.pixmap="on";
	xset("wdim",490,340);

	
	//scf(3); 
	xdel(3);
	scf(3);
	clf(3);
	xset("wpos",514,310)
	f=gcf();
	f.pixmap="on";
	xset("wdim",490,340);

	
	scf(1);
endfunction





while %T
	 select x_choose(["gaussienne";"uniforme"],"type de loi","Quitter"),
	
	 	case 0 then	
	 		break;
	
		case 1 then
			
			nbre_pts_gauss=0;
			echantrgauss=[];
			echantthetagauss=[];


			trois_figures();
			
			rand("normal");
			cadre_haut_gauss();
			
			while %T
			
				xselect();
				nbre_col_histo=min(nbre_pts_gauss,ncolhisto);

				Choix=x_choose(["Un point";...
					     "+20";"+50";"+100";"+1000";... 
					    "Histos";"Cond. par l''angle";"Cond. par le rayon"],...
					    string(nbre_pts_gauss)+" points gaussiens","Retour" );
				
				select Choix,
				
				case 0 then
					clear nbre_pts_gauss;
					clear echantrgauss;
					clear echantthetagauss;
					break;
		
				case 1 then
					x=rand(1,1);
					y=rand(1,1);
					plot2d(0,0,style=1,rect=[-Rmax,-Rmax,Rmax,Rmax], axesflag=5,frameflag=1);
					a=gca();
					a.isoview="on";
				    xpoly(x,y);
				    p=get("hdl");
				    p.mark_foreground=couleur;
				    p.mark_style=1;
				    p.thickness=1;
					echantrgauss=[echantrgauss,sqrt(x^2+y^2)];
					echantthetagauss=[echantthetagauss,atan(y,x)];
					couleur=changecouleur(couleur);	
					nbre_pts_gauss=nbre_pts_gauss+1;
					xselect();
					show_pixmap();
					
				case 5 then
					plot2d(0,0,style=1,rect=[-Rmax,-Rmax,Rmax,Rmax], axesflag=5,frameflag=1);
					a=gca();
					a.isoview="on";		
					laux=floor(vecteur_nombres(Choix)/paquet_monochrome);
					newr=[];
					newtheta=[];
					for i=1:laux
						x=rand(1,paquet_monochrome);
						y=rand(1,paquet_monochrome);
						xpoly(x,y);
					    p=get("hdl");
					    p.mark_foreground=couleur;
					    p.mark_style=1;
					    p.thickness=1;
						p.line_mode="off";
						newr=[newr,sqrt(x^2+y^2)];
						newtheta=[newtheta,atan(y,x)];
						couleur=changecouleur(couleur);
						xselect();
						show_pixmap();
					end	
					echantrgauss=[echantrgauss,newr];
					echantthetagauss=[echantthetagauss,newtheta];
					nbre_pts_gauss=nbre_pts_gauss+laux*paquet_monochrome;

				case 6 then
					if nbre_pts_gauss<n_sample_histo
						n_missing=n_sample_histo-nbre_pts_gauss;
						x=rand(1,n_missing);
						y=rand(1,n_missing);
						echantrgauss=[echantrgauss,sqrt(x^2+y^2)];
						echantthetagauss=[echantthetagauss,atan(y,x)];
						nbre_pts_gauss=n_sample_histo;
					end
					cbgrg();
					show_pixmap();	
					for j=1:nbre_histo_fluct
						cbgrg();
						n=max(floor(nbre_pts_gauss*j/nbre_histo_fluct),2);
						histplot(linspace(0,1.5*Rmax,nbre_col_histo),echantrgauss(1:n),rect=[0,0,1.5*Rmax,1], frameflag=1);
						show_pixmap();
					end
					cbdtg();
					show_pixmap();	
					for j=1:nbre_histo_fluct
						cbdtg();
						n=max(floor(nbre_pts_gauss*j/nbre_histo_fluct),2);
						histplot(nbre_col_histo,echantthetagauss(1:n),rect=[-%pi,0,%pi,0.4], frameflag=1);
						show_pixmap();
					end
					scf(1);
					
				case 7 then
					if nbre_pts_gauss<n_sample_histo
						n_missing=n_sample_histo-nbre_pts_gauss;
						x=rand(1,n_missing);
						y=rand(1,n_missing);
						echantrgauss=[echantrgauss,sqrt(x^2+y^2)];
						echantthetagauss=[echantthetagauss,atan(y,x)];
						nbre_pts_gauss=n_sample_histo;
					end
					plot2d(0,0,style=1,rect=[-Rmax,-Rmax,Rmax,Rmax], axesflag=5,frameflag=3);
					a=gca();
					a.isoview="on";		
					uaux=[5*Rmax,0,5*Rmax];
					vaux=[5*Rmax*tan(theta0),0,5*Rmax*tan(theta1)];
					xpoly(uaux,vaux);
					p=get("hdl");
					p.thickness=2;
					xarcs([-2 -2.5; 2 2.5; 4 5; 4 5; 0 0; theta0*180*64/%pi theta1*180*64/%pi ]);
					theta1mid=(theta1+theta0)/2;
					theta0mid=theta0/3;
					xstring( 2.6*cos(theta1mid), 2.6*sin(theta1mid),"b");
					t=get("hdl");  //get the handle of the newly created object
  					//t.font_foreground=6;
					t.font_size=5;
					t.font_style=5;
					xstring( 2.1*cos(theta0mid),2.1*sin(theta0mid),"a");
					t=get("hdl");  //get the handle of the newly created object
  					//t.font_foreground=6;
					t.font_size=5;
					t.font_style=5;
					xselect();
					show_pixmap();
					xclick();
					echantrgausscond=echantrgauss((echantthetagauss>theta0)&(echantthetagauss<theta1));
					ncondgr=length(echantrgausscond);
					cbdrg();
					show_pixmap();	
					cbgrg();
					histplot(linspace(0,1.5*Rmax,nbre_col_histo),echantrgauss,rect=[0,0,1.5*Rmax,1], frameflag=1);
					show_pixmap();
					for j=1:nbre_histo_fluct
						cbdrg();
						n=max(floor(ncondgr*j/nbre_histo_fluct),2);
						histplot(linspace(0,1.5*Rmax,nbre_col_histo),echantrgausscond(1:n),rect=[0,0,1.5*Rmax,1], frameflag=1);
						show_pixmap();
					end				
					scf(1);
					a=gca();
					xclick();
					delete(a.children([1,2,3,4]));
					clear echantrgausscond;
					show_pixmap();
					
				case 8 then
					if nbre_pts_gauss<n_sample_histo
						n_missing=n_sample_histo-nbre_pts_gauss;
						x=rand(1,n_missing);
						y=rand(1,n_missing);
						echantrgauss=[echantrgauss,sqrt(x^2+y^2)];
						echantthetagauss=[echantthetagauss,atan(y,x)];
						nbre_pts_gauss=n_sample_histo;
					end
					plot2d(0,0,style=1,rect=[-Rmax,-Rmax,Rmax,Rmax], axesflag=5,frameflag=1);
					a=gca();
					a.isoview="on";		
					affiche_cercles(0.5,1.5);
					xstring( 0.6*cos(%pi/6), 0.6*sin(%pi/6),"A");
					t=get("hdl");  //get the handle of the newly created object
  					//t.font_foreground=6;
					t.font_size=4;
					t.font_style=5;
					xstring( 1.6*cos(%pi/6), 1.6*sin(%pi/6),"B");
					t=get("hdl");  //get the handle of the newly created object
  					//t.font_foreground=6;
					t.font_size=5;
					t.font_style=5;
					xselect();
					show_pixmap();
					xclick();
					echantthetagausscond=echantthetagauss((echantrgauss>0.5)&(echantrgauss<1.5));
					ncondgt=length(echantthetagausscond);
					cbdtcg();
					show_pixmap();	
					cbgtg();
					histplot(nbre_col_histo,echantthetagauss,rect=[-%pi,0,%pi,0.4], frameflag=1);
					show_pixmap();
					for j=1:nbre_histo_fluct
						cbdtcg();
						n=max(floor(ncondgt*j/nbre_histo_fluct),2);
						histplot(nbre_col_histo,echantthetagausscond(1:n), rect=[-%pi,0,%pi,0.4], frameflag=1);
						show_pixmap();
					end
					scf(1);
					a=gca();
					xclick();
					delete(a.children([1,2,3]));
					clear echantthetagausscond;
					show_pixmap();
					
				else 
				   
						plot2d(0,0,style=1,rect=[-Rmax,-Rmax,Rmax,Rmax], axesflag=5,frameflag=1);
						a=gca();
						a.isoview="on";
						newr=[];
						newtheta=[];
						laux=vecteur_nombres(Choix);
						bigx=rand(1,laux);
						bigy=rand(1,laux);
						for i=1:laux
							x=bigx(i);
							y=bigy(i);
						    xpoly(x,y);
						    p=get("hdl");
						    p.mark_foreground=couleur;
						    p.mark_style=1;
						    p.thickness=1;
							newr=[newr,sqrt(x^2+y^2)];
							newtheta=[newtheta,atan(y,x)];
							couleur=changecouleur(couleur);	
							xselect();
							show_pixmap();
						end
						nbre_pts_gauss=nbre_pts_gauss+laux;	
						echantrgauss=[echantrgauss,newr];
						echantthetagauss=[echantthetagauss,newtheta];
				
				end
			end //(pour le while)
		
		case 2 then
			
			trois_figures();
			rand("uniform");
			cadre_haut_unif();

			nbre_pts_unif=0;
			echantrunif=[];
			echantthetaunif=[];
			
			
			while %T
				xselect();
				nbre_col_histo=min(nbre_pts_unif,ncolhisto);

				Choix=x_choose(["Un point";...
					      "+20";"+50";"+100";"+1000";...  
					    "Histos";"Cond. par l''angle";"Cond. par le rayon"],...
					    string(nbre_pts_unif)+" points uniformes","Retour" );
				
				select Choix,
				
				case 0 then
					clear nbre_pts_unif;
					clear echantrunif;
					clear echantthetaunif;
					break;
		
				case 1 then
					x=2*(rand(1,1)-0.5);
					y=2*(rand(1,1)-0.5);
					plot2d(0,0,style=1,rect=[-Rmaxu,-Rmaxu,Rmaxu,Rmaxu], axesflag=5,frameflag=1);
					a=gca();
					a.isoview="on";
				    xpoly(x,y);
				    p=get("hdl");
				    p.mark_foreground=couleur;
				    p.mark_style=1;
				    p.thickness=1;
					echantrunif=[echantrunif,sqrt(x^2+y^2)];
					echantthetaunif=[echantthetaunif,atan(y,x)];
					couleur=changecouleur(couleur);	
					nbre_pts_unif=nbre_pts_unif+1;
					xselect();
					show_pixmap();
					
				case 5 then
					if nbre_pts_unif<n_sample_histo
						n_missing=n_sample_histo-nbre_pts_unif;
						x=2*(rand(1,n_missing)-0.5);
						y=2*(rand(1,n_missing)-0.5);
						echantrunif=[echantrunif,sqrt(x^2+y^2)];
						echantthetaunif=[echantthetaunif,atan(y,x)];
						nbre_pts_unif=n_sample_histo;
					end
					plot2d(0,0,style=1,rect=[-Rmaxu,-Rmaxu,Rmaxu,Rmaxu], axesflag=5,frameflag=1);
					a=gca();
					a.isoview="on";		
					laux=floor(vecteur_nombres(Choix)/paquet_monochrome);
					newr=[];
					newtheta=[];
					for i=1:laux
						x=2*(rand(1,paquet_monochrome)-0.5);
						y=2*(rand(1,paquet_monochrome)-0.5);
						xpoly(x,y);
					    p=get("hdl");
					    p.mark_foreground=couleur;
					    p.mark_style=1;
					    p.thickness=1;
						p.line_mode="off";
						newr=[newr,sqrt(x^2+y^2)];
						newtheta=[newtheta,atan(y,x)];
						couleur=changecouleur(couleur);
						xselect();
						show_pixmap();
					end	
					echantrunif=[echantrunif,newr];
					echantthetaunif=[echantthetaunif,newtheta];
					nbre_pts_unif=nbre_pts_unif+laux*paquet_monochrome;

				case 6 then
					if nbre_pts_unif<n_sample_histo
						n_missing=n_sample_histo-nbre_pts_unif;
						x=2*(rand(1,n_missing)-0.5);
						y=2*(rand(1,n_missing)-0.5);
						echantrunif=[echantrunif,sqrt(x^2+y^2)];
						echantthetaunif=[echantthetaunif,atan(y,x)];
						nbre_pts_unif=n_sample_histo;
					end
					cbgru();
					show_pixmap();	
					for j=1:nbre_histo_fluct
						cbgru();
						n=max(floor(nbre_pts_unif*j/nbre_histo_fluct),2);
						histplot(nbre_col_histo,echantrunif(1:n),rect=[0,0,Rmaxu,2], frameflag=1);
						show_pixmap();
					end
					cbdtu();
					show_pixmap();	
					for j=1:nbre_histo_fluct
						cbdtu();
						n=max(floor(nbre_pts_unif*j/nbre_histo_fluct),2);
						histplot(nbre_col_histo,echantthetaunif(1:n),rect=[-%pi,0,%pi,0.35], frameflag=1);
						show_pixmap();
					end
					scf(1);
					
				case 7 then
					if nbre_pts_unif<n_sample_histo
						n_missing=n_sample_histo-nbre_pts_unif;
						x=2*(rand(1,n_missing)-0.5);
						y=2*(rand(1,n_missing)-0.5);
						echantrunif=[echantrunif,sqrt(x^2+y^2)];
						echantthetaunif=[echantthetaunif,atan(y,x)];
						nbre_pts_unif=n_sample_histo;
					end
					plot2d(0,0,style=1,rect=[-Rmaxu,-Rmaxu,Rmaxu,Rmaxu], axesflag=5,frameflag=1);
					a=gca();
					a.isoview="on";	
					uaux=[5*Rmaxu,0,5*Rmaxu];
					vaux=[5*Rmaxu*tan(theta0u),0,5*Rmaxu*tan(theta1u)];
					xpoly(uaux,vaux);
					p=get("hdl");
					p.thickness=2;
					xarcs([-0.4 -0.8; 0.4 .8; .8 1.6; .8 1.6; 0 0; theta0u*180*64/%pi theta1u*180*64/%pi ]);
					theta1mid=theta1u/3;
					theta0mid=theta0u/3;
					xstring( 0.45*cos(theta0mid), 0.45*sin(theta0mid),"a");
					t=get("hdl");  //get the handle of the newly created object
  					//t.font_foreground=6;
					t.font_size=5;
					t.font_style=5;
					xstring( .85*cos(theta1mid),.85*sin(theta1mid),"b");
					t=get("hdl");  //get the handle of the newly created object
  					//t.font_foreground=6;
					t.font_size=5;
					t.font_style=5;
					xselect();
					show_pixmap();
					xclick();
					echantrunifcond=echantrunif((echantthetaunif>theta0u)&(echantthetaunif<theta1u));
					ncondur=length(echantrunifcond);
					cbdru();
					show_pixmap();	
					cbgru();
					histplot(nbre_col_histo,echantrunif,rect=[0,0,Rmaxu,2], frameflag=1);
					show_pixmap();
					for j=1:nbre_histo_fluct
						cbdru();
						n=max(floor(ncondur*j/nbre_histo_fluct),2);
						histplot(nbre_col_histo,echantrunifcond(1:n),rect=[0,0,Rmaxu,2], frameflag=1);
						show_pixmap();
					end
					scf(1);
					a=gca();
					xclick();
					delete(a.children([1,2,3,4]));
					clear echantrunifcond;
					show_pixmap();
					
				case 8 then
					if nbre_pts_unif<n_sample_histo
						n_missing=n_sample_histo-nbre_pts_unif;
						x=2*(rand(1,n_missing)-0.5);
						y=2*(rand(1,n_missing)-0.5);
						echantrunif=[echantrunif,sqrt(x^2+y^2)];
						echantthetaunif=[echantthetaunif,atan(y,x)];
						nbre_pts_unif=n_sample_histo;
					end				
					plot2d(0,0,style=1,rect=[-Rmaxu,-Rmaxu,Rmaxu,Rmaxu], axesflag=5,frameflag=1);
					a=gca();
					a.isoview="on";		
					affiche_cercle(Rcondu);
					xstring( 1.05*Rcondu*cos(%pi/15), 1.05*Rcondu*sin(%pi/15),"A");
					t=get("hdl");  //get the handle of the newly created object
  					//t.font_foreground=6;
					t.font_size=4;
					t.font_style=5;
					p=get("hdl");
					xselect();
					show_pixmap();
					xclick();
					echantthetaunifcond=echantthetaunif((echantrunif>Rcondu));
					ncondut=length(echantthetaunifcond);
					cbdtcu();
					show_pixmap();	
					cbgtu();
					histplot(nbre_col_histo,echantthetaunif,rect=[-%pi,0,%pi,0.35], frameflag=1);
					show_pixmap();
					for j=1:nbre_histo_fluct
						cbdtcu();
						n=max(floor(ncondut*j/nbre_histo_fluct),2);
						histplot(nbre_col_histo,echantthetaunifcond(1:n), rect=[-%pi,0,%pi,0.35], frameflag=1);
						show_pixmap();
					end
					scf(1);
					a=gca();
					xclick();
					delete(a.children([1,2]));
					clear echantthetaunifcond;
					show_pixmap();

				else 
				   
						plot2d(0,0,style=1,rect=[-Rmaxu,-Rmaxu,Rmaxu,Rmaxu], axesflag=5,frameflag=1);
						a=gca();
						a.isoview="on";
						newr=[];
						newtheta=[];
						laux=vecteur_nombres(Choix);
						bigx=2*(rand(1,laux)-0.5);
						bigy=2*(rand(1,laux)-0.5);
						for i=1:laux
							x=bigx(i);
							y=bigy(i);
						    xpoly(x,y);
						    p=get("hdl");
						    p.mark_foreground=couleur;
						    p.mark_style=1;
						    p.thickness=1;
							newr=[newr,sqrt(x^2+y^2)];
							newtheta=[newtheta,atan(y,x)];
							couleur=changecouleur(couleur);	
							xselect();
							show_pixmap();
						end
						nbre_pts_unif=nbre_pts_unif+laux;	
						echantrunif=[echantrunif,newr];
						echantthetaunif=[echantthetaunif,newtheta];
				
				end
			end //(pour le while)		
	end	//pour le select
end //(pour le while)

scf(1);	
f=gcf();
clf();
f.pixmap="off";
scf(2);
f=gcf();
clf();
f.pixmap="off";
scf(3);
f=gcf();
clf();
f.pixmap="off";
driver("Rec");