Proprietà statistiche di una serie di lanci di dadi!

/*Proprieta' statistiche di una sequenza di lanci di dadi */
#include
#include
#include

int main(void) {
int i, N_lanci=1000,max=3,min=18,a,b,c,sum=0,sum_tot=0;
int hist[19];
double med=0.;


for(i = 0; i < 19; ++i) hist[i]=0;

for(i = 0; i < N_lanci; ++i){
a = 1+((double) rand () /( RAND_MAX + 1.))*6;
b = 1+((double) rand () /( RAND_MAX + 1.))*6;
c = 1+((double) rand () /( RAND_MAX + 1.))*6;

sum = a + b + c ;
sum_tot += sum;
if (sum > max ) max = sum;
if (sum < min ) min = sum;
hist[sum] +=1;
}

med = sum_tot / 1000.;
printf("media di %d lanci di 3 dadi= %f\n",N_lanci,med);
printf("massimo valore su tutti i lanci di 3 dadi: %d\n",max);
printf("minimo valore su tutti i lanci di 3 dadi: %d\n",min);
printf("*** Istogramma\n");
printf("Valore\tNumero di occorrenze\n");
for (i=0; i<19; i++){
if (hist[i] != 0) printf( "%d\t%d\n",i,hist[i] );
}
return 0;
}