diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /260/CH3/EX3.4/3_4.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '260/CH3/EX3.4/3_4.sce')
| -rw-r--r-- | 260/CH3/EX3.4/3_4.sce | 45 |
1 files changed, 45 insertions, 0 deletions
diff --git a/260/CH3/EX3.4/3_4.sce b/260/CH3/EX3.4/3_4.sce new file mode 100644 index 000000000..dc05f6003 --- /dev/null +++ b/260/CH3/EX3.4/3_4.sce @@ -0,0 +1,45 @@ +//Eg-3.4
+//pg-74
+
+clear
+clc
+
+ // Matrices A and B (AX=B)...given by 3 sets of material balance equations
+ a=[.8 .02 .06;.1 .83 .12;.1 .15 .82];
+ [n,n]=size(a);
+ b=[50;30;20];
+
+ //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
+
+//Resuts
+P=X(1);
+Q=X(2);
+R=X(3);
+
+//Displaying results
+disp(P,"P=")
+disp(Q,"Q=")
+disp(R,"R=")
\ No newline at end of file |