avevar.dem

	PROGRAM D13R5
C	Driver for routine AVEVAR
	PARAMETER(NPTS=1000, EPS=0.1)
	DIMENSION DATA(NPTS)
C	Generate Gaussian distributed data
	IDUM=-5
	WRITE(*,'(1X,T4,A,T14,A,T26,A)') 'Shift','Average','Variance'
	DO 12 I=1,11
		SHIFT=(I-1)*EPS
		DO 11 J=1,NPTS
			DATA(J)=SHIFT+I*GASDEV(IDUM)
11		CONTINUE
		CALL AVEVAR(DATA,NPTS,AVE,VAR)
		WRITE(*,'(1X,F6.2,2F12.2)') SHIFT,AVE,VAR
12	CONTINUE
	END