Template Presentasi Skripsi Powerpoint 468x60

Monday 14 May 2012

Matlab Eleminasi Gaus

clc,clear all,close all
disp('==============================================')
disp('  Program Penyelesaian Persamaan Linear Ax=B  ')
disp('             Metode eliminasi Gauss           ')
disp('----------------------------------------------')
disp('==============================================')
disp(' ')
m=input('banyaknya baris adalah ');
n=input('banyaknya kolom adalah ');
for i=1:m
    for j=1:n
        A(i,j)=input('masukan matriks = ');
    end
end
for k=1:n
    B(k)=input('Masukan vektor = ');
end
A
B=B'
if m~=n
    disp (['Baris - ',num2str(m),';Kolom - ',num2str(n)])
    disp ('Bukan Matriks Persegi')
    disp('Matrik [A] Tidak Mempunyai Invers..!')
    return
end
detA=det(A);
if detA==0
    disp('A adalah Matriks singular')
    disp('Matrik (A)tidak mempunyai invers..1')
    return
end
disp('langkah awal : Membuat Matrik M = A[B]')
disp('-----------------------------------------------------')
A=[A B ]
for i=1:m
    disp(['Eliminasi Kolom' num2str(i)])
    disp('-----------------------------------------------------')
    disp(['B' num2str(i) '/a' num2str(i)])
    A(i,:)=A(i,:)./A(i,i)
    io=i-1;
    in=m-i;
    for ii=1:io
        disp(['B' num2str(ii) '-a' num2str(ii) num2str(i) '*B' num2str(i)])
        A(ii,:)=A(ii,:)-A(i,:).*A(ii,i)
    end
     for jj=1:in
        disp(['B' num2str(i+jj) '-a' num2str(i+jj) num2str(i) '*B' num2str(i)])
        A(i+jj,:)=A(i+jj,:)-A(i,:).*A(i+jj,i)
     end
end
disp('SOLUSI PERSAMAAN LINEAR ADALAH ')
disp(' ')
for i=1:m
    disp(['X' num2str(i) '=' num2str(A(i,m+1))])
end

Setelah dimasukkan ke dalam program MATLAB, maka berbentuk seperti gambar dibawah ini :





No comments:

Post a Comment