2100 



PARA=1.0+1.0/C 



2110 



CALL GAMMA (PARA, HUE) 



2120 



nEAN(2l=SI6MAtMME 



2130 



FACl=SI6MA«»2iHME»*2 



2140 



PARA=1.0+2.0/C 



2150 



CALL 6AMMA(PARA,MV2) 



2160 



FAC2=SI6HAt»2»HV2 



2170 



VAR(2)=FAC2-FAC1 



21B0 



PARA=1.0-1.0/U 



2190 



CALL GAMMA ( PARA, HPC) 



2200 



MEAN(3)=SIGHA2»HPC 



2210 



PARA=1,0-2.B/U 



2220 



CALL 6AMMA(PARA,HPD) 



2230 



VAR(3)=SIGMA2t«2»HPD-MEAN(3)t»2 



2240 





2250C WRITE OUT THE DATA FOR EACH DISTRIBUTION 



2260 



MRITE(6,135) TITLE(l) 



2270 135 



FORMAT(///,20)(,A26,///) 



2280 



DO 150 K=l,3 



229B 



KTEMP=K 



2300 



WRITE(6,160) IFLA6(K),DEF,F0RM(K1 



2310 160 



FORMAT(151(,A30,//,ll(,Ai7,2X,A34) \ 



2320 



IF( K .EQ. 1) «RITE(6,159) EPSI,PHI 



2330 159 



FORMAT(1)1/EPSI="6X,F10.3,/,1X,"PHI=",7K,F10.3) 



2340 



IF( K .GT. 1) WRITE(6,161) ALPHA(K),BETA(K) 



2350 161 



FORMAT(1)(,"ALPHA=",51(,F10.4,/,11(,''BETA=%6X,F10.3) 



2360 



MRITE(6,162) HEAN(K),VAR(K) 



2370 162 



FORMAT(1X,"MEAN=",6X,F10.3,/,1X,»VARIANCE=%21(,F10.3) 



23S0 



IF( LOGIC .EQ. 'N') GO TO 207 



2390 



DO 170 1=1, N 



24B0 



DUM1(I)=YACT(1,K) 



2410 



DUH2(I)=YEST(1,K) 



2420 



H0LD1=C0RR(K) 



2430 



L2=N 



2440 170 



CONTINUE 



2450 



CALL RESIDUAL(HS,DUM1,DUM2,H0LD1,L2) 



2460 207 



«RITE(6,208) 



2470 208 



FORMAT (7X, "RETURN PERIOD TABLE", /,6X, "YEAR", 13X,"H5") 



2480 



DO 211 J=l,5 



2490 



HRITE(6,212) RET(J),CHS(J,K) 



2500 212 



F0RMAT(1X,F9.2,8X,F9.2) 



2510 211 



CONTINUE 



2520 



HRITEi6,165) 



2530 165 



FORMAT(////) 



2540 150 



CONTINUE 



2550 303 



RETURN 



2560 



END 



2570 





2580C SUBROUTINE TO HELP PRINT OUT DATA 



2590 



SUBROUTINE RESIDUAL(X,YACT,YEST,CORR,N) 



2600 



DIMENSION X(N),YACT(N),YEST(N),DIFF(200) 



2610 



SSR=0 



2620 



DO 10 I=1,N 



CIO 



