1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
//Example 3.6
//Gauss-Jordan Method
//Page no. 57
clc;clear;close;
A=[10,1,1,12;1,10,1,12;1,1,10,12]; //augmented matrix
for i=1:3
j=i
while (A(i,i)==0 & j<=3)
for k=1:4
B(1,k)=A(j+1,k)
A(j+1,k)=A(i,k)
A(i,k)=B(1,k)
end
disp(A)
j=j+1
end
for k=4:-1:i
A(i,k)=A(i,k)/A(i,i)
end
disp(A)
for k=1:3
if(k~=i) then
l=A(k,i)/A(i,i)
for m=i:4
A(k,m)=A(k,m)-l*A(i,m)
end
end
end
disp(A)
end
for i=1:3
printf('\nx(%i) = %g\n',i,A(i,4))
end
|
