/* Driver for routine bnldev */#include <stdio.h>#define NRANSI#include "nr.h"#define N 20#define NPTS 1000#define ISCAL 200#define NN 100#define LLEN 50int main(void){	char txt[LLEN+1];	int i,j,k,klim,dist[N+1];	long idum=(-133);	float pp,xm,dd;	for (;;) {		for (j=0;j<=N;j++) dist[j]=0;		do {			printf("Mean of binomial distribution (0.0 to %d.0)",N);			printf(" - Negative to end:\n");			scanf("%f",&xm);		} while (xm > N);		if (xm < 0.0) break;		pp=xm/NN;		for (i=1;i<=NPTS;i++) {			j=bnldev(pp,NN,&idum);			if (j >= 0 && j <= N) ++dist[j];		}		printf("Binomial-distributed deviate, mean %5.2f of %6d points\n",			xm,NPTS);		printf("%4s %8s %10s\n","x","p(x)","graph:");		for (j=0;j<N;j++) {			for (k=0;k<=LLEN;k++) txt[k]=' ';			dd=(float) dist[j]/NPTS;			klim=(int) (ISCAL*dd+1);			if (klim > LLEN) klim=LLEN;			for (k=1;k<=klim;k++) txt[k]='*';			txt[LLEN]='\0';			printf("%4d %9.4f   %s\n",j,dd,txt);		}	}	return 0;}#undef NRANSI