xprobks.c

来自「< C语言数值算法程序大全>>配套程序」· C语言 代码 · 共 33 行

/* Driver for routine probks */

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

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

main()
{
	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;
}