diff options
Diffstat (limited to '260/CH3/EX3.2/3_2.sce')
-rw-r--r-- | 260/CH3/EX3.2/3_2.sce | 42 |
1 files changed, 42 insertions, 0 deletions
diff --git a/260/CH3/EX3.2/3_2.sce b/260/CH3/EX3.2/3_2.sce new file mode 100644 index 000000000..f1673aa33 --- /dev/null +++ b/260/CH3/EX3.2/3_2.sce @@ -0,0 +1,42 @@ +//Eg-3.2
+//pg-69
+
+clear
+clc
+
+
+ // Matrices A and B (AX=B)
+ a=[1 1 -1;1 2 -2;-2 1 1];
+ [n,n]=size(a);
+ b=[1;0;1];
+
+ //Augumented matrix of A and B
+ auga=[a b];
+
+ //Algorithm of Naive gauss elimination
+ //Forward elimination
+ for k=1:n-1
+ for i=(k+1):n
+ factr=auga(i,k)/auga(k,k);
+ auga(i,:)=auga(i,:)-factr*auga(k,:);
+ end
+ end
+
+ //Initializing X
+ X=zeros(n,1);
+
+ //Backward substitution
+ X(n)=auga(n,n+1)/auga(n,n);
+ for i=(n-1):-1:1
+ summ=auga(i,n+1);
+ for j=(i+1):n
+ summ=summ-auga(i,j)*X(j);
+ end
+ X(i)=summ/auga(i,i);
+ end
+ //Displaying Solution
+ disp(X(1),"x0=")
+ disp(X(2),"x1=")
+ disp(X(3),"x2=")
+
+
\ No newline at end of file |