#include <stdio.h>
//#include <math.h>


double uniform(double a,double b,long int *seed)
{
	double t;
	*seed=16807*(*seed)+1;
	*seed=abs(*seed-(*seed/2147483646)*2147483646);
	t=(*seed)/2147483646.0;
	t=a+(b-a)*t;
	return(t);
}


//main()
//{
//	int i;
//	int c=100;
//	double b;
//	long int s=13579;
//	
//	for(i=0;i<100;i++)
//	{
//		b=uniform(0,1,&s);
//		printf("%13.4f",b);
//	}
//}