316B C(2>=.5772156649B15329 



3170 C (3) =-.6558780715282538 



3180 C (4)=-. 0420026350340952 



3190 C(5)=. 1665386113822915 



3200 C(6)=-.8421977345555443 



3210 C(7)=-. 009621971527887 



3220 C (8) =.087218943246663 



3230 C(9)=-. 0011651675918591 



3240 C(10)=-.0002152416741149 



3250 C(ll)=. 0001280582823882 



3268 C(12)=-. 0000201348547807 



3270 C(13)=-. 0080012504934821 



3280 C(14)=. 0000011330272320 



3290 C(15)=-. 0000002056338417 



3300 C(16)=6.116095E-e9 



3310 C(17)=5.0020075E-09 



3328 C(18)=-1.1812746E-89 



3338 C(19)=1.843427E-10 



3340 C(20)=7.7823E-12 



3350 C(21)=-3.69680E-12 



3360 C(22)=5.1E-13 



3378 C(23)=-2.06E-14 



3380 C(24)=-5.4E-15 



3390 C(25)=1.4E-15 



3488 



34i0C sun SERIES 



3420 SUH=0.0 



3430 DO 50 K=l,25 



3440 SUI1=SUM+C(K)»IALPHAHK) 



3458 58 AREA=6F/SUn 



3468 200 RETURN 



3470 END 



3480 



3498 SUBROUTINE EXPCT(M,V2, CUT, ALPHA, BETA, LAMBDA) 



3508 DIMENSION ALPHA (3), BETA (3) 



351B DOUBLE PRECISION BU 



3528 REAL LAHBDA 



3538 FD1(1()=-(AL06(-AL06(X))»PHI)+EPSI 



3548 FD2(X) = ((-AL0B(1-I())*t(l/A1))»B1 



355B FD3(l()=B2/((-AL06(X))»*(l/A2)) 



3S6B PDF1(1()=E1(P(-EXP(-(X-EPSI)/PHI))»EXP(-(X-EPSI)/PHI)/PHI 



357B PDF2(X)=Al*(X»*(Al-l))tEXP(-(X/Bl)*»Al)/(BltfAl) 



358B PDF3(X)=A2»(B2»»A2)»EXP(-(B2/X)«A2)/(Xt*(A2+l)) 



359B CDF1(X)=EXP(-EXP((EPSI-X)/PHI)) 



36BB CDF2(X)=1.0-EXP(-(X/B1)"A1) 



3618 CDF3(X)=EXP(-(B2/X)t*A2) 



3620 G(X)=««(l.0-EXP(V2*(X-CUT))) 



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



3648 PRINT, 'EXTREMAL TYPE I...r 



3658 PRINT, 'WEIBULL 2' 



3668 PRINT, 'L06-EXTREMAL 3' 



367B PRINT, 'SELECT 1, 2, OR 3' 



36BB READ, ID 



C12 



