next up previous contents index
suivant: Pour s'amuser avec les monter: Factorielles aléatoires précédent: La solution mathématique   Table des matières   Index

La modélisation avec Xcas

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 :

Image factalea

Documentation de giac écrite par Renée De Graeve