Komputasi untuk Sains dan Teknik - Universitas Indonesia

2.4. ALGORITMA ELIMINASI GAUSS 2719 WRITE (*,’(1X,A)’) ’MATRIK AUGMENT:’20 DO 110 I = 1,N21 WRITE (*,’(1X,5(F14.8))’) (A(I,J),J=1,M)22 110 CONTINUE23 WRITE (*,*)24 C LANGKAH 2: MEMERIKSA ELEMEN-ELEMEN PIVOT DAN PROSES TUKAR POSISI25 NN = N-126 DO 10 I=1,NN27 C LANGKAH 3: MENDEFINISIKAN P28 P = I29 100 IF (ABS(A(P,I)).GE.1.0E-20 .OR. P.GT.N) GOTO 20030 P = P+131 GOTO 10032 200 IF(P.EQ.N+1)THEN33 C MENAMPILKAN PESAN TIDAK UNIK34 WRITE(*,5)35 GOTO 40036 END IF37 C LANGKAH 4: PROSES TUKAR POSISI38 IF(P.NE.I) THEN39 DO 20 JJ=1,M40 C = A(I,JJ)41 A(I,JJ) = A(P,JJ)42 A(P,JJ) = C43 20 CONTINUE44 END IF45 C LANGKAH 5: PERSIAPAN PROSES TRIANGULARISASI46 JJ = I+147 DO 30 J=JJ,N48 C LANGKAH 6: TENTUKAN MJI49 MJI = A(J,I)/A(I,I)50 C LANGKAH 7: MELAKUKAN PROSES TRIANGULARISASI51 DO 40 K=JJ,M52 A(J,K) = A(J,K)-MJI*A(I,K)53 40 CONTINUE54 A(J,I) = 055 30 CONTINUE56 10 CONTINUE57 C MENAMPILKAN HASIL TRIANGULARISASI58 WRITE (*,’(1X,A)’) ’HASIL TRIANGULARISASI:’59 DO 120 I = 1,N60 WRITE (*,’(1X,5(F14.8))’) (A(I,J),J=1,M)61 120 CONTINUE62 C LANGKAH 8: MEMERIKSA ELEMEN A(N,N)63 IF(ABS(A(N,N)).LT.1.0E-20) THEN64 C MENAMPILKAN PESAN TIDAK UNIK65 WRITE(*,5)66 GOTO 40067 END IF68 C LANGKAH 9: MENGHITUNG X(N)69 X(N) = A(N,N+1)/A(N,N)70 C LANGKAH 10: PROSES SUBSTITUSI MUNDUR71 L = N-172 DO 15 K=1,L73 I = L-K+174 JJ = I+175 SUM = 0.076 DO 16 KK=JJ,N77 SUM = SUM+A(I,KK)*X(KK)