From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 3369/CH2/EX2.11/Ex2_11.sce | 20 ++++++++++++++++++++ 1 file changed, 20 insertions(+) create mode 100644 3369/CH2/EX2.11/Ex2_11.sce (limited to '3369/CH2/EX2.11') diff --git a/3369/CH2/EX2.11/Ex2_11.sce b/3369/CH2/EX2.11/Ex2_11.sce new file mode 100644 index 000000000..67670ce53 --- /dev/null +++ b/3369/CH2/EX2.11/Ex2_11.sce @@ -0,0 +1,20 @@ +//Chapter 2, Exmaple 11, page 75 +//Calculate the potential within the mesh +clc +clear +//Based on figure 2.38(b) +//equations are obtained using Eq.2.46 +A1 = 1/2*(0.54+0.16) +A2 = 1/2*(0.91+0.14) +S = [0.5571 -0.4571 -0.1;-0.4751 0.828 0.3667;-0.1 0.667 0.4667] +//By obtaining the elements of the global stiffness matrix(Sadiku,1994) +//and by emplying the Eq.2.49(a) +S1 = [1.25 -0.014;-0.014 0.8381] +S2 = [-0.7786 -0.4571;-0.4571 -0.3667] +Phi13 = [0 ;10] +val1 = S2*Phi13 +Phi24 = S1\val1 +disp(-Phi24,"The values of Phi2 and Phi4 are:") + +//Answers may vary due to round of error + -- cgit