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