Adrian Lawrenz F-Matlab
Adrian Lawrenz F-Matlab
Pesebre
Meeting No. 1
MATLAB CODE:
a= [ 2 4 6 28
3 5 7 34
6 1 5 23 ];
%gauss elimination
[m,n]=size(a);
for j=1:m-1
for z=2:m
if a(j,j)==0
t=a(j,:);a(j,:)=a(z,:);
a(z,:)=t;
end
end
for i=j+1:m
a(i,:)=a(i,:)-a(j,:)*(a(i,j)/a(j,j));
end
end
x=zeros(1,m);
for s=m:-1:1
c=0;
for k=2:m
c=c+a(s,k)*x(k);
end
x(s)= (a(s,n)-c)/a(s,s);
end
disp('Gauss elimination method:');
a
x'
display:
a=
2 4 6 28
0 -1 -2 -8
0 0 9 27
ans =
1
2
3