Download NASA Technical Memorandum X-2417
Transcript
C - F P H A S ALL ONES MOVED C.:'fi\lOW S E T U P TM OA K E P llPPER LEFT nF DIAGnl\lAL. T OE X T R E M E P = D ( K l 3 +1) PERMIITATION MATRIX C.- Z E R O OIJT P , WILL R E ZEROS AND ONES D O 3 9 I = 1 , ONT C - Z E RH OO I J S E K E E P I NAGR R AOYN LNYE EFDI R SCTO L I I M N . JL(I1 = IZ K = K D + I 39 D ( K ) = OFZO L =1 M =1 n L =1 7803 J = CHECK FOR IF( I = JL(0L) TEEDS ROW OF D I A G . T H AN J .EO. I Z ) G OT O7 R 6 (K-l)*OS + J + YD = 1.D0 nL = O L + 1 C / F / A R KT H I S 1 AS ZRO T O B EF I L L E D - 7901 JL(M) = K M = Y +1 A 1 MOVED I N T OI T D(I) I TI S MOVEn IJP AND OIJTQF HERE Grl TO 780 + L P A R T OF I D E N T I T YA T 7 8 6 D C N -N T EED TO MOVE 1 T n A HnLE. O ( J ) = 1.00 780 L = OR1 + L C. R E T U R NM. A T R I C E S COMPLETED E W I T HI R li3S D E L I R E S A T F L YL E F T OUT. RETIJRN 6 9 1 CALL L N C N T ( 1) IdRITE (6,1691 1 OSIONT 7H6 T,.MAKE c- n io J = KD