suivant: Factorielles aléatoires
monter: Nombres triangulaires aléatoires
précédent: La solution mathématique
Table des matières
Index
On tape le programme trialea(r,q,p) qui tire au hasard des nombres entre
1 et r. On fait p fois des échantillons de taille q, et on dessine
les résultats intermédiaires obtenus :
l contient les sommes cumulées des résultats (ici une somme) c'est
à dire la
somme d'un échantillon de taille n=k+1+j*q avec k=0..q-1 et
j=0..p-1. Dans Ldiv on met evalf(l/n) lorsque
n=q,2*q...p*q
trialea(r,q,p):={
local j,k,l,n,LdivN,alea;
LdivN:=NULL;
l:=0;
n:=0;
for (j:=0;j<p;j++){
for (k:=0;k<q;k++){
alea:=(rand(r)+1);
while (alea!=1){
l:=l+alea;
alea:=(rand(r)+1);
}
l:=l+1;
n:=n+1;
}
LdivN:=LdivN,evalf(l/n);
}
return LdivN;
}:;
On tape :
L10:=trialea(10,100,10000);
plotlist(L10)
On obtient :
Documentation de giac écrite par Renée De Graeve