initial commit / add all books

1 files changed, 27 insertions, 0 deletions

diff --git a/3760/CH2/EX2.19/Ex2_19.sce b/3760/CH2/EX2.19/Ex2_19.sce
new file mode 100644
index 000000000..baf6ef5e9
--- /dev/null
+++ b/3760/CH2/EX2.19/Ex2_19.sce
@@ -0,0 +1,27 @@
+clc;

+// L1=3+(1/(2*X)) -self inductance of coil 1

+// L1=2+(1/(2*X)) -self inductance of coil 2

+// M=1/(2*X) -mutual inductance 

+// X is displacement 

+i1=10; // current of coil 1

+i2=-5; // current of coil 2

+// from expression W=(L1*i1^2)/2 + (L2*i2^2)/2 + (i1*i2*M) where W is work done and is equal to 175+25/(4*x);

+// fe=-25/(4*x^2) where fe is electromagnetic force and is equal to rate of change of work done with respect to x

+disp('case a');

+x1=0.5;

+x2=1; 

+Wm=integrate('-25/(4*x^2)','x',x1,x2);

+printf('Mechanical work done for given increment in displacement is %f watt-sec\n',Wm);

+disp('case b')

+// psi1=(3+(1/(2*x)))*i1+(1/(2*x))*i2 psi1 is flux linkage with coil 1  

+// psi2=(2+(1/(2*x)))*i2+(1/(2*x))*i1 psi2 is flux linkage with coil 2 

+We1=(((3+(1/(2*x2)))*i1+(1/(2*x2))*i2)-((3+(1/(2*x1)))*i1+(1/(2*x1))*i2))*i1;

+We2=(((2+(1/(2*x2)))*i2+(1/(2*x2))*i1)-((2+(1/(2*x1)))*i2+(1/(2*x1))*i1))*i2;

+printf('Electrical energy supplied by source 1 is %f watt-sec\n',We1);

+printf('Electrical energy supplied by source 2 is %f watt-sec\n',We2);

+disp('case c');

+// W=175+(25/(4*x));

+dw=175+(25/(4*x2))-(175+(25/(4*x1)));

+printf('Change in field energy is %f Watt-sec\n',dw);

+

+