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 --- 25/CH2/EX2.3/2_3.sce | 23 +++++++++++++++++++++++ 1 file changed, 23 insertions(+) create mode 100755 25/CH2/EX2.3/2_3.sce (limited to '25/CH2/EX2.3') diff --git a/25/CH2/EX2.3/2_3.sce b/25/CH2/EX2.3/2_3.sce new file mode 100755 index 000000000..a949b68fb --- /dev/null +++ b/25/CH2/EX2.3/2_3.sce @@ -0,0 +1,23 @@ +// example:-2.3,page no.-31. +//program to find the resulting fields by assumibg plane waves on either side of the current sheet and enforcing the boundary conditions. +syms E x E1 E2 H1 H2 z Jo A B c N n d ko y; +sym('n*(E2-E1)=0'); //boundary condition to be satisfied at z=0 +sym('z*(E2-E1)=0'); // " " " " +sym('n*(H2-H1)=Jo'); // " " " " +sym('z*(H2-H1)=Jo'); // " " " " +E1=A*N*exp(%i*ko*z)*x; // x component of electric field (region z<0). +H1=A*N*exp(%i*ko*z)*(-y); // -y component of magnetic field (region z<0). +E2=B*N*exp(-%i*ko*z)*x; // x component of electric field (region z>0). +H2=B*N*exp(-%i*ko*z)*y; // y component of electric field (region z>0). +disp(E1,'for z<0, E1=') +disp(H1,'for z<0, H1=') +disp(E2,'for z>0, E2=') +disp(H2,'for z>0, H2=') +//from boundary conditions imposed.we get:- +c=[-1 -1;1 -1]; +d=[A;B]; +c*d==[Jo;0]; +d=inv(c)*[Jo;0]; +//result +// A=-Jo/2; B=-Jo/2. +disp(d) \ No newline at end of file -- cgit