double dbinom(long,long);
double dlgama(double);
void ierm1(STRING,long,long,STRING,STRING,long,byte);
void ierv1(STRING,long,byte);