/* Copyright Will Wilson 1999 */
/* birth-death simulation (multiple run output) */

#include	<stdlib.h>
#include	<stdio.h>

#define	MAXTIME 	100
#define	PRNTTIME	10
#define	NUMRUNS		10
#define RUNTYPE		1 /* 1-zero stays extinct; 2-reseeding */

int main(void)
{
	int	ttt, event, n, irun;
	int	seed;
	double	alpha, beta;
	double	ave25=0.0, ave50=0.0, ave75=0.0, ave100=0.0;

	alpha = beta = 0.10;

	seed = 1456739853;
	srand48(seed);

	for(irun=0;irun<NUMRUNS;irun++)
	{
		n = 10;
		for(ttt=0;ttt<=MAXTIME;ttt++)
		{
			event = 0;
			if(drand48()<alpha) event = 1;
			if(drand48()<beta) event -= 1;
			n += event;
			if(n<0) n=0;
#if(RUNTYPE==1)
			if(n==0) ttt = MAXTIME;
#endif
			if(ttt==25) ave25 += n;
			if(ttt==50) ave50 += n;
			if(ttt==75) ave75 += n;
			if(ttt==100) ave100 += n;
		}
	}

	printf(" 25 %f\n",ave25/NUMRUNS);
	printf(" 50 %f\n",ave50/NUMRUNS);
	printf(" 75 %f\n",ave75/NUMRUNS);
	printf("100 %f\n",ave100/NUMRUNS);

	return (0);
}