METODE NUMERIK LANJUTAN
SISTEM PERSAMAAN LINIER DENGAN METODE INVERS
MATRIKS
LA’LA MONICA (18202109005)
UNIVERSITAS SAM RATULANGI
PASCA SARJANA
PROGRAM STUDI TEKNIK SIPIL
2018
PERMASALAHAN
Cari nilai x1, x2, dan x3 pada persamaan (1),(2), dan (3) berikut.
2x1 + 4x2 + x3 = 11 ......(1)
-x1 + 3x2 – 2x3 = -16 ......(2)
2x1 – 3x2 + 5x3 = 21 ......(3)
Dibuat dalam bentuk matriks pada persamaan (4) dan (5).
[𝐴]{𝑥} = {𝐵} ......(4)
2 4 1 𝑥1 −11
[−1 3 −2] {𝑥2} = {−16}.......(5)
2 −3 5 𝑥3 21
DIAGRAM ALIR PEMOGRAMAN
START
INPUT NILAI
MATRIKS A DAN B
MENGHITUNG
KOFAKTOR MATRIKS A
[Ac]
MENGHITUNG
DETERMINAN MATRIKS A
D
MENGHITUNG
ADJOINT MATRIKS A
[Aa]
MENGHITUNG INVERS
MATRIKS A
[Ai]
MENGHITUNG NILAI
X1, X2, DAN X3
END
BAHASA PEMOGRAMAN FORTRAN YANG DIGUNAKAN
PRINT *, "TUGAS 1 METODE NUMERIK LANJUTAN"
PRINT *, "LA'LA MONICA"
PRINT *, "18202109005"
PRINT *, "========================================================================="
PRINT *, "INVERS MATRIKS 3X3"
PRINT *, "Mencari nilai dari X1, X2, dan X3"
PRINT *, "-------------------------------------------------------------------------"
PRINT *, "Input nilai matriks A dan nilai B"
PRINT *, "A11 ="
READ *, A11
PRINT *, "A12 ="
READ *, A12
PRINT *, "A13 ="
READ *, A13
PRINT *, "A21 ="
READ *, A21
PRINT *, "A22 ="
READ *, A22
PRINT *, "A23 ="
READ *, A23
PRINT *, "A31 ="
READ *, A31
PRINT *, "A32 ="
READ *, A32
PRINT *, "A33 ="
READ *, A33
PRINT *, "B1 ="
READ *, B1
PRINT *, "B2 ="
READ *, B2
PRINT *, "B3 ="
READ *, B3
PRINT *, "-------------------------------------------------------------------------"
PRINT *, "Kofaktor Matriks A"
Ac11=(A22*A33)-(A23*A32)
PRINT *, "Ac11 =", Ac11
Ac12=(-1)*((A21*A33)-(A23*A31))
PRINT *, "Ac12 =", Ac12
Ac13=(A21*A32)-(A22*A31)
PRINT *, "Ac13 =", Ac13
Ac21=(-1)*((A12*A33)-(A13*A32))
PRINT *, "Ac21 =", Ac21
Ac22=(A11*A33)-(A13*A31)
PRINT *, "Ac22 =", Ac22
Ac23=(-1)*((A11*A32)-(A12*A31))
PRINT *, "Ac23 =", Ac23
Ac31=(A12*A23)-(A13*A22)
PRINT *, "Ac31 =", Ac31
Ac32=(-1)*((A11*A23)-(A13*A21))
PRINT *, "Ac32 =", Ac32
Ac33=(A11*A22)-(A12*A21)
PRINT *, "Ac33 =", Ac33
PRINT *, "-------------------------------------------------------------------------"
PRINT *, "Adjoint Matriks A"
Aa11=Ac11
Aa12=Ac21
Aa13=Ac31
Aa21=Ac12
Aa22=Ac22
Aa23=Ac32
Aa31=Ac13
Aa32=Ac23
Aa33=Ac33
PRINT *, "Aa11 =", Aa11
PRINT *, "Aa12 =", Aa12
PRINT *, "Aa13 =", Aa13
PRINT *, "Aa21 =", Aa21
PRINT *, "Aa22 =", Aa22
PRINT *, "Aa23 =", Aa23
PRINT *, "Aa31 =", Aa31
PRINT *, "Aa32 =", Aa32
PRINT *, "Aa33 =", Aa33
PRINT *, "-------------------------------------------------------------------------"
PRINT *, "Determinan Matriks A"
D=((A11*A22*A33)+(A12*A23*A31)+(A13*A21*A32))-
((A13*A22*A31)+(A11*A23*A32)+(A12*A21*A33))
PRINT *, D
PRINT *, "-------------------------------------------------------------------------"
PRINT *, "Invers Matriks A"
Ai11=Aa11/D
Ai12=Aa12/D
Ai13=Aa13/D
Ai21=Aa21/D
Ai22=Aa22/D
Ai23=Aa23/D
Ai31=Aa31/D
Ai32=Aa32/D
Ai33=Aa33/D
PRINT *, "Ai11 =", Ai11
PRINT *, "Ai12 =", Ai12
PRINT *, "Ai13 =", Ai13
PRINT *, "Ai21 =", Ai21
PRINT *, "Ai22 =", Ai22
PRINT *, "Ai23 =", Ai23
PRINT *, "Ai31 =", Ai31
PRINT *, "Ai32 =", Ai32
PRINT *, "Ai33 =", Ai33
PRINT *, "-------------------------------------------------------------------------"
PRINT *, "Nilai X1, X2, dan X3"
X1=(Ai11*B1)+(Ai12*B2)+(Ai13*B3)
X2=(Ai21*B1)+(Ai22*B2)+(Ai23*B3)
X3=(Ai31*B1)+(Ai32*B2)+(Ai33*B3)
PRINT *, "X1 =", X1
PRINT *, "X2 =", X2
PRINT *, "X3 =", X3
STOP
END PROGRAM
HASIL DARI PROGRAM YANG DIBUAT
