diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /25/CH2/EX2.3 | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '25/CH2/EX2.3')
-rwxr-xr-x | 25/CH2/EX2.3/2_3.sce | 23 |
1 files changed, 23 insertions, 0 deletions
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 |