#include "k.h"
/* 
*	Generate outliers using a relative risk (gamma) which perturbs
*   the 1/2,1/2 distribution somewhat.                     
*
*   Input m - number of alleles, gamma
* 	 *idum - pointer to random number work field.
*
*	Output rslt - vector of length 3: counts of each genotype 
*/
void datdis(int m, double gamma, long *idum, int *rslt) {

double p1;

p1=gamma/(1.+gamma);						/* calc p1				*/
datmult(m,p1,idum,rslt); 
return;
}