PROGRAM MAIN CC THIS PROGRAM LISTS THE GEOMETRIC DISTRIBUTION FOR A MEAN DURATION AND CC CHOSEN MINIMUM PHASE DURATION IMPLICIT DOUBLE PRECISION (A-H,O-Z) OPEN(UNIT=5,FILE='CHIOUT1.TXT') WRITE(*,1200) 1200 FORMAT(1H ,'INPUT AVERAGE DURATION ') READ(*,*)XAVG WRITE(*,1300) 1300 FORMAT(1H ,'INPUT MINIMUM PHASE ') READ(*,*)XMIN WRITE(*,1310) 1310 FORMAT(1H ,'INPUT SAMPLE SIZE ') READ(*,*)XSAMP XAVG = XAVG - XMIN XP = 1.0D0/(XAVG + 1.0D0) WRITE(5,1380) WRITE(5,1390)XAVG,XP,XSAMP SUM = 0.0D0 DO 100 J=1,200 XPROB = XP*(1.0D0-XP)**(J-1) XPRSA = XPROB*XSAMP SUM = XPRSA + SUM JJ = J -1 1380 FORMAT(1H ,'THIS PROGRAM PROVIDES THE EXPECTED NUMBER OF', & ' OBSERVATIONS PER EXIT TIME.'/' WE SUM THE EXPECTED VALUES TO', & ' OBTAIN A MINIMUM OF 6, AND THEN START OVER.'/' USE THIS', & ' TO HELP CONSTRUCT THE CHI-SQUARE TEST, CALCULATED BY CHISQ.FOR' & /) 1390 FORMAT(1H ,'AVERAGE DURATION, HAZARD, N ',3(F10.5,1X)/) WRITE(5,1400)JJ,XPROB,XPRSA,SUM 1400 FORMAT(1H ,'J,PROB,EXPECTED,SUM ',I3,1X,3(F10.5,1X)) IF(SUM.GE.6.0D0)SUM = 0.0D0 100 CONTINUE STOP END