1 files changed, 32 insertions, 0 deletions
diff --git a/3588/CH6/EX6.3/EX6_3.sce b/3588/CH6/EX6.3/EX6_3.sce
new file mode 100644
index 000000000..2c5f799b1
--- /dev/null
+++ b/3588/CH6/EX6.3/EX6_3.sce
@@ -0,0 +1,32 @@
+//Clearing console
+clc
+clear
+
+x = poly(0,"x")
+y = poly(0,"y")
+
+//Intializing Variables
+r = 1
+s = 0.5
+
+//Node 2: (3*m + b=1)
+//Node 3: (2.5*m + b=2)
+//constructing matrices for solving b and m
+k = [3 1;2.5 1]
+f = [1;2]
+
+//Solving for b and m
+u(1:2,1)=linsolve(k,-f)
+
+//Calculating x and y for (r, s) = (1, 0.5)
+X =  (1 - r )*(1 - s)/4 + (1 + r )*(1 - s)*(3)/4 + (1 + r )*(1 + s)*(2.5)/4 + (1 - r )*(1 + s)*(1.25)/4
+Y =(1 - r )*(1 - s)/4 + (1 + r )*(1 - s)/4 + (1 + r )*(1 + s)*(2)/4 + (1 - r )*(1 + s)*(1.75)/4
+
+
+
+//Printing Results
+printf('\nResults\n')
+printf('\nelement edge 2-3 is described by y=')
+disp(u(1,1)*x + u(2,1))
+printf('\nFor (r, s) = (1, 0.5) , we obtain')
+printf('\n x=%f y=%f',X,Y)