-
Notifications
You must be signed in to change notification settings - Fork 2
/
gauss_seidel.m
51 lines (46 loc) · 839 Bytes
/
gauss_seidel.m
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
39
40
41
42
43
44
45
46
47
48
49
50
function [result] = gauss_seidel(A,b,x)
%GAUSS_SEIDEL Summary of this function goes here
% Detailed explanation goes here
n = length(A);
X = x;
Xold = X;
converged = 0;
iter = 0;
I = eye(3);
L = zeros(3);
U = zeros(3);
for i = 1:3
for j = 1:i
L(i,j) = A(i,j);
end
end
for i = 1:3
for j = i:3
U(i,j) = A(i,j);
end
end
C = -inv(I+L)*U;
disp(L)
disp(U)
while converged == 0
iter = iter+1;
Xold = X;
for i = 1:n
sigma = 0;
for j = 1:n
if j ~= i
sigma = sigma + A(i,j)*X(j);
end
end
X(i) = (1/A(i,i))*(b(i) - sigma);
end
X
for i = 1:n
if abs(X(i) - Xold(i)) < 0.000001
converged = 1;
else
converged = 0;
end
end
end
end