209(9 FD2(X) = ( (-ALOG(l-X) )»»(1/A1) )«B1 



2100 FD3(X)=B2/( (-ALOG(X) )«t(l/A2)) 



2110 PDFl (X)=>EXP(-EXP(-(X-EPSI)/PHI))*EXP(-(X-EPSI)/PHI)/PHI 



2120 PDF2(X)=Al»(X»t(Al-l))»EXP(-(X/Bl)»»Al)/(Bl»»Al) 



2130 PDF3(X)=A2«(B2»*A2)»EXP(-(B2/X)»*A2)/(X»*(A2+1) ) 



2140 CDFl (X)=EXP(-EXP((EPSI-X)/PHI) ) 



2150 CDF2(X)=1.0-EXP(-(X/B1)««A1) 



2160 CDF3(X)=EXP(-(B2/X)*»A2) 



2170 G(X)=W«(1.0-EXP(V2«(X-CUT))) 



2180 70 PRINT, 'SELECT A DISTRIBUTION...' 



2190 PRINT, 'EXTREMAL TYPE I...1' 



2200 PRINT, 'WEIBULL 2' 



2210 PRINT, 'LOG-EXTREMAL 3' 



2220 PRINT, 'SELECT 1, 2, OR 3' 



2230 READ, ID 



2240 IF(ID .LT. 1 .OR. ID .GT. 3) 60 TO 70 



2250 IF(ID .EQ. 1) WRITE(6,104) 



2260 104 F0RMAT(/, IX, "INPUT EXTREMAL TYPE I EPSILON, AND PHI ") 



2270 IF(ID .EQ. I) READ, EPSI , PHI 



2280 IF( ID .EQ. 2) WRITE(6,114) 



2290 114 F0RMAT(/, IX, "INPUT WEIBULL ALPHA AND BETA") 



2300 IF( ID .EQ. 2) READ,A1,B1 



2310 IF( ID .EQ. 3) WRITE(6,124) 



2320 124 FORMAT(/, IX, "INPUT LOG-EXTREMAL ALPHA AND BETA") 



2330 IF( ID .EQ. 3) READ,A2,B2 



2340 WRITE(6,5) 



2350 5 FORMAT(/, IX, "INPUT AVERAGE NUMBER OF EXTREMAL EVENTS PER YEAR,' 



2360 t> /,1X,"THE POISSON 'LAMBDA' PARAMETER") 



2370 READ, LAMBDA 



2380 BL='CUT 



2390 IF(ID .EQ. 1) BU = FD1 (1 -. 0000001 ) 



2400 IFdD .EQ. 2) BU=FD2 (1-. 0000001 ) 



2410 IFdD .EQ. 3) BU = FD3 (1-. 0000001) 



2420 SUM=0 



2430 D=BU-BL 



2440 K=-l 



2450 DO 10 1=1,101 



2460 K=-K 



2470 IF( K .LT. 0) C=4 



2480 IF( K .GT. )C=2 



2490 IF( I .EQ. 1 .OR. I .EQ. 101)C=1 



2500 ADD=FLOAT(I-1)*D/100.0 



2510 XV=BL+ADD „ .^ 



2520 IFdD .EQ. 1 .AND. EXP ( - ( X V-EPS I ) /PH I ) .GT. 82.0) GO TO 10 



0530 IFdD .EQ. 2 .AND. ( ( X V/ B 1 ) »*A1 ) .GT. 82.0) 60 TO 10 



2540 IFdD .EQ. 3 .AND. ( (B2/XV) ♦♦A2) .GT. 82.0) GO TO 10 



2550 IF( ID .EQ. 1) FAC1=PDF 1 ( XV) 



2560 IF( ID .EQ. 2) FAC 1 =PDF2 ( XV) 



2570 IF( ID .EQ. 3) FAC 1 =PDF3 ( XV) 



2580 FAC2=6(XV) 



2590 IF(XV .LT. CUT) FAC2=0 



2600 SUM=SUM+C«FAC1*FAC2 



2610 10 CONTINUE 



A9 



