#include "uniform.c"

int bernoulli(double p,long int *seed)
{
	int x;

	double u;
	double uniform();
	
	u=uniform(0.0,1.0,seed);
	x=(u<=p)?1:0;
	return(x);
}

//main()
//{int s=13257;
//	//bernoulli(0.5,&s);
//	int i;
//	for(i=0;i<10;i++)
//	{printf("%d\n",bernoulli(0.5,&s));}
//}