xprobks.c

Numerical Recipes Software 提供的算法子程序集

/* Driver for routine probks */

#include <stdio.h>
#define NRANSI
#include "nr.h"

#define NPTS 20
#define EPS 0.1
#define ISCAL 40

int main(void)
{
	int i,j,jmax;
	char txt[ISCAL+1];
	float alam,aval;

	printf("probability function for kolmogorov-smirnov statistic\n\n");
	printf("%7s %10s %13s\n","lambda","value:","graph:");
	for (i=1;i<=NPTS;i++) {
		alam=i*EPS;
		aval=probks(alam);
		jmax=(int) (0.5+(ISCAL-1)*aval);
		for (j=0;j<ISCAL;j++) {
			if (j < jmax)
				txt[j]='*';
			else
				txt[j]=' ';
		}
		txt[ISCAL]='\0';
		printf("%8.6f %10.6f       %s\n",alam,aval,txt);
	}
	return 0;
}
#undef NRANSI
/* (C) Copr. 1986-92 Numerical Recipes Software . */