double p = drand48();
double prob[5]={0.08, 0.11, 0.03, 0.16, 0.20};
for(i=1; i<5; i++) prob[i] += prob[i-1];
for(i=0; i<5; i++){
if (p < prob[i]) break;
}
result = i;