blob: 6f2dd870bb737da664c6cb27a5c2c761968e7a8a (
plain)
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
|
//Eg-3.17
//pg-114
clear
clc
A=[-3.5 1 1.5;1 4 -1;-2 -.6 -3.5];
B=[2.5;4;-16];
es=10^-5;
imax=10;
[r,c] = size(A)
n = r;
X=[0;0;0];
iter=1;
lambda=1;
while iter<imax//condition for termination
for i=1:n
summ=B(i);
pivot=A(i,i);
if pivot==0
error('gsie not applicable');//to avoid a/0 forms
end
old=X(i);
for j=1:n
if i~=j
summ=summ-A(i,j)*X(j);
end
end
X(i)=(lambda*summ/pivot)+(1-lambda)*old;//relaxation
end
iter=iter+1;
end
disp("Solution after 10 iterations")
disp(X)
|
