next up previous contents index
suivant: Factorielles aléatoires monter: Nombres triangulaires aléatoires précédent: La solution mathématique   Table des matières   Index

La modélisation avec Xcas

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 :

Image trialea

Documentation de giac écrite par Renée De Graeve