Komputasi untuk Sains dan Teknik - Universitas Indonesia

2.6. MENGHITUNG INVERS MATRIK 356 C LANGKAH 1: MEMASUKAN NILAI ELEMEN-ELEMEN MATRIK A7 WRITE (*,’(1X,A)’) ’JUMLAH PERSAMAAN ? ’8 READ (*,*) N9 WRITE (*,*)10 WRITE (*,*) ’MASUKAN ELEMEN-ELEMEN MATRIK A’11 M = N + 112 DO 50 I = 1,N13 DO 60 J = 1,N14 WRITE (*,’(1X,A,I2,A,I2,A)’) ’A(’,I,’,’,J,’) = ’15 READ (*,*) A(I,J)16 60 CONTINUE17 50 CONTINUE18 C LANGKAH 2: MENDEFINISIKAN MATRIK IDENTITAS19 WRITE (*,*) ’MENDEFINISIKAN MATRIK IDENTITAS’20 DO 70 I = 1,N21 DO 80 J = M,N+N22 A(I,J) = 023 IF (I+N .EQ. J) THEN24 A(I,J) = 125 END IF26 80 CONTINUE27 70 CONTINUE28 WRITE (*,*)29 C MENAMPILKAN MATRIK AUGMENT30 WRITE (*,’(1X,A)’) ’MATRIK AUGMENT:’31 DO 110 I = 1,N32 WRITE (*,’(1X,5(F14.8))’) (A(I,J),J=1,N+N)33 110 CONTINUE34 WRITE (*,*)35 C MENGHITUNG JUMLAH TUKAR (TKR) POSISI. MULA2 TKR = 036 TKR = 037 C MENGHITUNG JUMLAH OPERASI BAGI/KALI (BK).38 BK = 039 C MENGHITUNG JUMLAH OPERASI TAMBAH/KURANG (TK).40 TK = 041 C LANGKAH 3: MEMERIKSA ELEMEN2 PIVOT DAN PROSES TUKAR POSISI42 NN = N-143 DO 10 I=1,NN44 C LANGKAH 4: MENDEFINISIKAN P45 P = I46 100 IF (ABS(A(P,I)).GE.1.0E-20 .OR. P.GT.N) GOTO 20047 P = P+148 GOTO 10049 200 IF(P.EQ.N+1)THEN50 C MENAMPILKAN PESAN SINGULAR51 WRITE(*,5)52 GOTO 40053 END IF54 C LANGKAH 5: PROSES TUKAR POSISI55 IF(P.NE.I) THEN56 DO 20 JJ=1,N+N57 C = A(I,JJ)58 A(I,JJ) = A(P,JJ)59 A(P,JJ) = C60 TKR = TKR + 161 20 CONTINUE62 END IF63 C LANGKAH 6: PERSIAPAN PROSES TRIANGULARISASI64 JJ = I+1