1 files changed, 57 insertions, 0 deletions
diff --git a/260/CH3/EX3.12/3_12.sce b/260/CH3/EX3.12/3_12.sce
new file mode 100644
index 000000000..d5cf922be
--- /dev/null
+++ b/260/CH3/EX3.12/3_12.sce
@@ -0,0 +1,57 @@
+//Eg-3.12
+//pg-100
+
+clear
+clc
+
+ A=[6 15 55 ;15 55 225;55 225 979];
+ B=[74.5;262.3;1078.1];
+ 
+
+n=3;
+
+
+matsol=zeros(3,3);//initialising matrix matsol
+Z=A;
+Y=B;
+I=[1 0 0;0 1 0;0 0 1];//creates an identity matrix of size n*n
+X=zeros(3,1);
+inverse=zeros(3,3);
+ABI=zeros(3,7);
+ABI(:,:)=[A(:,:) B(:,:) I(:,:)];//augumented matrix of A,B,I
+  for k=1:n//proceeds from 1st row to the last row
+        u=k;
+        big=abs(ABI(k,k));          
+        for t=k+1:n  //this loop is for selecting the elementhaving max absolute value in a column                   
+        dummy=abs(ABI(t,k));     
+        if dummy>big then
+            big=dummy;
+            u=t;
+        end
+    end
+        if u~=k then
+        for j=1:2*n+1//interchanging rows to make max absolute element as pivot 
+            dummy=ABI(u,j);
+            ABI(u,j)=ABI(k,j);
+            ABI(k,j)=dummy;
+        end
+        end
+        pivot=ABI(k,k);
+       
+        
+        ABI(k,:)=ABI(k,:)/pivot;
+   
+        for i=1:n
+            if i~=k;
+            factor=ABI(i,k);
+            ABI(i,:)=ABI(i,:)-ABI(k,:)*factor;
+        end
+    end
+    X(:,:)=[ABI(:,n+1)];//determining X using augumented matrix
+
+    matsol=inv(Z)*Y;//using matlab inbuilt functions to get X
+
+end
+
+disp("Coefficients of parabola")
+disp(X)
+\ No newline at end of file