//
//  programme 7
//
write(%io(2),'programme 7: trace de la frequence de pannes entre 0 et T d''une machine')
write(%io(2),'et de la repartition autour de la valeur moyenne ce nombre') 
write(%io(2),'//')
//
Ness=evstr(x_dialog('nombre d''appareils testes?','100'));
//
T=evstr(x_dialog('duree du test (s)?','1000'));
//
xbasc() 
//
// parametres du calcul 
//
lambda=0.5;   // parametre de la loi de Poisson
//
R=[];Z=[]; 
//
for i=1:Ness 
  t=0; 
  Ncont=-1; 
  while(t<T)  
    Ncont=Ncont+1; 
    y=(-1/lambda)* log(rand()); 
    t=t+y; 
  end; 
 u=sqrt(T)*(Ncont/T-lambda); 
 R=[R,Ncont/T]; 			// nombre moyen de panne observe
 Z=[Z,u]; 				// ecart a la moyenne theorique 
end; 
//
// trace des resultats
//
xsetech([0,0,1,0.5]) 
plot2d(1:Ness,R,-1) 
xtitle('frequence de pannes','appareil no','freq')
xsetech([0,0.5,1,0.5]) 
histplot(-2:0.1:2,Z) 
xtitle('fonction de repartition de l''ecart','ecart',' ')
//
 
 
 
 
 
