PROGRAM GAUSS JORDAN pada fortran f90

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