TUGAS
PRAKTIKUM METODE NUMERIK
Disusun
oleh:
Nama : Indah Novitasari
NIM : 1207045011
Prodi : FISIKA
LABORATORIUM FISIKA
KOMPUTASI DAN PERMODELAN
FAKULTAS
MATEMATIKA ILMU PENGETAHUAN ALAM
UNIVERSITAS MULAWARMAN
SAMARINDA
2013
TUGAS
Carilah
nilai A B C D E dari persamaan berikut ini dengan menggunakan program !!
23A -6B 8C 26D -67E = 28
9A 8B -23.1C
28D -8E = 67.3
3A 11B 26C 32.3D 9E = 89.10
7A -12B 72C -6.7D 9.2E = 70.1
8A 3.4B -81C 8.7D -8.2E =
3
PENYELESAIAN
!PROGRAM AKAR
PERSAMAAN LINIER METODE GAUSS-JORDAN
parameter (n=5)
dimension g(n,n+1),X(n)
open(1,file='a.txt',status='old')
read(1,*)((g(i,j),j=1,n+1),i=1,n)
write (*,*)'*********************************'
write (*,*)'NAMA : INDAH NOVITASARI '
write (*,*)'NIM : 1207045011'
write (*,*)'PRODI : FISIKA'
write (*,*)'*********************************'
write (*,*)''
write (*,*)'MATRIKS TUGAS
PRAKTIKUM METODE NUMERIK'
write (*,*)''
write (*,*)'------------------------------------------------'
print*,'MATRIKS AWAL'
print*,''
write(*,'6f8.1')((g(i,j),j=1,n+1),i=1,n)
a=g(1,1);b=g(2,1);c=g(3,1);d=g(4,1);e=g(5,1)
do j=1,n+1
g(1,j)=g(1,j)/a;g(2,j)=g(2,j)-(g(1,j)*b);g(3,j)=g(3,j)-(g(1,j)*c)
g(4,j)=g(4,j)-(g(1,j)*d);g(5,j)=g(5,j)-(g(1,j)*e)
enddo
write (*,*)'------------------------------------------------'
print*,''
print*,'LANGKAH 1'
print*,''
write(*,'6f8.3')((g(i,j),j=1,n+1),i=1,n)
z=g(2,2);y=g(1,2);q=g(3,2);o=g(4,2);r=g(5,2)
do j=1,n+1
g(2,j)=g(2,j)/z;g(1,j)=g(1,j)-(g(2,j)*y);g(3,j)=g(3,j)-(g(2,j)*q)
g(4,j)=g(4,j)-(g(2,j)*o);g(5,j)=g(5,j)-(g(2,j)*r)
enddo
print*,''
print*,'LANGKAH 2'
print*,''
write(*,'6f8.3')((g(i,j),j=1,n+1),i=1,n)
s=g(4,3);p=g(5,3)
do j=1,n+1
g(3,j)=g(3,j)/54.824;g(1,j)=g(1,j)-(g(3,j)*-0.313)
g(2,j)=g(2,j)-(g(3,j)*-2.535)
g(4,j)=g(4,j)-(g(3,j)*s);g(5,j)=g(5,j)-(g(3,j)*p)
enddo
print*,''
print*,'LANGKAH 3'
print*,''
write(*,'6f8.3')((g(i,j),j=1,n+1),i=1,n)
f=g(4,4);o=g(5,4)
do j=1,n+1
g(4,j)=g(4,j)/f;g(1,j)=g(1,j)-(g(4,j)*1.629)
g(2,j)=g(2,j)-(g(4,j)*2.121);g(3,j)=g(3,j)-(g(4,j)*0.157)
g(5,j)=g(5,j)-(g(4,j)*o)
enddo
print*,''
print*,'LANGKAH 4'
print*,''
write(*,'6f8.3')((g(i,j),j=1,n+1),i=1,n)
t=g(5,5);v=g(1,5);u=g(2,5);w=g(3,5);h=g(4,5)
do j=1,n+1
g(5,j)=g(5,j)/t
g(1,j)=g(1,j)-(g(5,j)*v)
g(2,j)=g(2,j)-(g(5,j)*u)
g(3,j)=g(3,j)-(g(5,j)*w)
g(4,j)=g(4,j)-(g(5,j)*h)
enddo
print*,''
print*,'LANGKAH 5'
print*,''
write(*,'6f8.3')((g(i,j),j=1,n+1),i=1,n)
print*,''
write (*,*)'-------------------'
print*,'MATRIKS HASIL'
print*,''
do i=1,n
write(*,3)i,g(i,n+1)
enddo
write (*,*)'-------------------'
!pembuktian
print*,''
print*,'PEMBUKTIAN'
print*,''
print*,'(23*g(1,j))+(-6*g(2,j))+(8*g(3,j))+(26*g(4,j))+(-67*g(5,j))='
do i=1,n
do j=1,n+1
aa=(23*g(1,j))+(-6*g(2,j))+(8*g(3,j))+(26*g(4,j))+(-67*g(5,j))
enddo
enddo
print*,''
write(*,'f5.1')aa
3 FORMAT (5X,'X(',I2,')=',F8.3)
End