Une des difficulté, lorsqu'on effectue des traitements statistiques avec Scilab est :

• de définir une densité de probabilité,
• de transformer un histogramme en densité de probabilité.

Pour cela, on est souvent confronté au problème de la normalisation d'une tableau de valeurs de telle sorte que ce tableau puisse être interprété comme l'échantillonné d'une densité de probabilité.

On se souvient alors qu'une condition nécessaire pour qu'une fonction puisse être interprétée comme une densité de probabilité est que son intégrale calculée sur son support vaille 1 :

La version discrète de cette intégrale est alors :

On considère alors une portion de script Scilab dans lequel on a construit un tableau devant être normalisé pour pouvoir être considéré comme une densité de probabilité. On a, par exemple, construit un histogramme qu'il faut normaliser :

  ...
  du = ...;
  ab = a:du:b;                      // la variable u est définie dans l'intervalle [a,b]
  [ind,h,info] = dsearch(U,ab,...); // le tableau h doit être normalisé
  H = sum(h)*du;                    // facteur de normalisation
  f = h/H;                          // normalisation de h
  ...

La fontion Scilab sum(x) est équivalente à l'expression :

  x*ones(x)'