™ = D ™ BESJ680 



TA=0. 



C 



TA = S 



C 



K=K-1 



BESJ6B1 



TB=T0L BESJ682 



KK = 1 B E S J 6 8 3 



605 CONTINUE BESJ684 



BESJ685 



C BACKWARD RECUR UNINDEXED BESJ686 



C BESJ687 



DO 601 1=1, IN BESJ688 



S=TB BESJ689 



TB=TM#TB-TA BESJ690 



TA=S BESJ691 



DTM=DTM-TRX BESJ692 



TM = DTM BESJ69 7 



601 CONTINUE BESJ694 



C NORMALIZATION BESJ695 



IF(KK.NE.l) GO TO 604 BESJ696 



TA=(TA/TB)*TEMP(3) BESJ697 



TB=TEMP(3) BESJ698 



KK=2 BESJ699 



IN=NS BESJ700 



IF(NS.NE.O) 60 TO 605 BESJ701 



604 Y(NN)=TB BESJ702 



615 NZ=N-NN BESJ703 



IF(NN.EQ.l) RETURN BESJ704 



S=TB BESJ705 



TB=TM#TB-TA BESJ706 



BESJ707 



DTM=DTM-TRX BESJ708 



TM=DTM BESJ709 



K=NN_1 BESJ710 



Y(K)=TB BESJ711 



IF(NN.EQ.2) RETURN BESJ712 



KM=K " 1 BESJ713 



BESJ714 



C BACKWARD RECUR INDEXED BESJ715 



BESJ716 



DO 602 1=1, KM BESJ717 



Y(K-1)=TM*Y(K)-Y(K+1> BESJ718 



DTM=DTM-TRX BESJ719 



TM=DTM BESJ720 



BESJ721 



602 CONTINUE BES^'" 



RETURN BESJ72^ 



C BESJ724 



C 



C 



BESJ725 

 BESJ726 



91 CONTINUE BESJ727 

 NZ= " 2 BESJ728 

 RETURN BESJ729 



92 CONTINUE BESJ730 

 NZ = _3 BESJ731 

 RETURN BESJ732 



93 CONTINUE BESJ733 

 NZ=_l BESJ734 



Figure A3. (Sheet 16 of 25) 



A20 



