diff options
Diffstat (limited to '260/CH3/EX3.17')
-rw-r--r-- | 260/CH3/EX3.17/3_17.sce | 38 |
1 files changed, 38 insertions, 0 deletions
diff --git a/260/CH3/EX3.17/3_17.sce b/260/CH3/EX3.17/3_17.sce new file mode 100644 index 000000000..6f2dd870b --- /dev/null +++ b/260/CH3/EX3.17/3_17.sce @@ -0,0 +1,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)
\ No newline at end of file |