suivant: Pour s'amuser avec les
monter: Factorielles aléatoires
précédent: La solution mathématique
Table des matières
Index
On tape le programme factalea(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 un produit) 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
factalea(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);
f:=1
while (alea!=1){
f:=f*alea;
alea:=(rand(r)+1);
}
l:=l+f;
n:=n+1;
}
LdivN:=LdivN,evalf(l/n);
}
return [LdivN];
}:;
On tape :
F3:=factalea(3,10,100);
plotlist(F3)
On obtient :
Documentation de giac écrite par Renée De Graeve