510 IFi LANBDA .LE. e ) 60 TD 17 



528 28 PRINT, 'INPUT THE NUMBER OF TRANSMITTED WAVE HEIGHT,' 



530 PRINT, 'RETURN PERIOD DATA POINTS YOU HAVE. ' 



548 READ,N2 



SSB IFI N2 .LT. 2) PRINT, 'YOU NEED AT LEAST TWO POINTS' 



S6B IFI N2 .LT. 2) 60 TO 28 



57B IFI N2 .6T. 188) PRINT, 'TOO MANY POINTS-108 IS MAXIMUM' 



SBB IFI N2 .6T. 188) 60 TO 28 



59B PRINT, 'INPUT THE DATA POINTS-ONE AT A TIME.' 



6BB PRINT, 'INPUT TRANSMITTED WAVE HEIGHT, COMMA, THEN RETURN PERIOD IN YEARS ' 



618 DO 38 1=1, N2 



620 READ,1(II),YII) 



630 38 YII)=1.0-1.0/ILA«BDA»YII)) 



640 



650C PUT TRANSMITTED WAVE HEI6HTS IN ORDER 



668 CALL ORDER IN2,X,Y) 



67B CALL WAVTRNSIX,Y, ALPHA, BETA, LAMBDA, N2) 



68B CALL EXPCTIC1,C2,C3, ALPHA, BETA, LAMBDA) 



69fl STOP 



708 END 



71B 



72flC SUBROUTINE TO PUT DATA IN ORDER-LOWEST TO HIGHEST 



738 SUBROUTINE 0RDERIN2,X,Y) 



748 DIMENSION XI100),YI100) 



750 DO 20 K=2,N2 



760 J=N2-K+2 



77B DO 10 1=1, J-1 



7BB IFI XII) .LT. Xd+D) 60 TO 10 



79B T1=XII) 



8BB T2=YII) 



BIB X{I)=XIIH) 



82B YII)=YII+1) 



83B X(I+1)=T1 



840 YII+1)=T2 



850 18 CONTINUE 



868 28 CONTINUE 



878 RETURN 



B6B END 



B9B 



9Bfl SUBROUTINE WAVTRNS IK, Y, ALPHA, BETA, LAMBDA, N) 



91B DIMENSION YACTI200,3) ,YESTI20B,3),DUM1I200),DUM2I200),HSI200),YI200) 



92B DIMENSION YAV6l3),C0RRI3),AI3),BI3),STI3),SBI3),ALPHAi3) ,BETA(3),VARI3) 



93fl DIMENSION RETI5),CHSI5, 3) 



94B REAL MEAN 1 3) 



95fl REAL LAMBDA 



96B F1IX)=EXPI-EXPI-IX-EPSI)/PHI)) 



979 ' F2IX)=1.0-EXPII-IX/SI6MA)*»C)) 



98B F3IX)=EXP(-IISI6MA2/X)HU)) 



998 CHARACTER*20 IFLA6I3) 



1000 CHARACTERtl7 DEF 



1018 CHARACTERI34 F0RMI3) 



1820 CHARACTER»26 TITLEIl) 



1030 CHARACTER*! L06IC 



C7 



