PROGRAM D13R10
C	Driver for routine KSONE
	PARAMETER(NPTS=1000,EPS=0.1)
	DIMENSION DATA(NPTS)
	EXTERNAL FUNC
	IDUM=-5
	WRITE(*,'(/1X,T5,A,T24,A,T44,A/)')
     *			'Variance Ratio','K-S Statistic','Probability'
	DO 12 I=1,11
		VAR=1.0+(I-1)*EPS
		FACTR=SQRT(VAR)
		DO 11 J=1,NPTS
			DATA(J)=FACTR*ABS(GASDEV(IDUM))
11		CONTINUE
		CALL KSONE(DATA,NPTS,FUNC,D,PROB)
		WRITE(*,'(1X,F14.6,F18.6,E20.4)') VAR,D,PROB
12	CONTINUE
	END

	FUNCTION FUNC(X)
	Y=X/SQRT(2.0)	
	FUNC=ERF(Y)
	END