initial commit / add all books

8 files changed, 119 insertions, 0 deletions

diff --git a/42/CH5/EX5.1/sadiku_5_1.sce b/42/CH5/EX5.1/sadiku_5_1.sce
new file mode 100755
index 000000000..e04e8159e
--- /dev/null
+++ b/42/CH5/EX5.1/sadiku_5_1.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+r=.2;

+disp('J=1/r3(2cosP ar + sinP a)')

+I=(2/r)*integrate('sin(P)*cos(P)','P',0,%pi/2)*integrate('1','Q',0,2*%pi);

+disp(I,'Current passing through Hemispherical shell');

+I=(2/r)*integrate('sin(P)*cos(P)','P',0,%pi,10^-10)*integrate('1','Q',0,2*%pi);

+disp(I,' Current through spherical shell= ');
\ No newline at end of file
diff --git a/42/CH5/EX5.10/sadiku_5_10.sce b/42/CH5/EX5.10/sadiku_5_10.sce
new file mode 100755
index 000000000..b99626f49
--- /dev/null
+++ b/42/CH5/EX5.10/sadiku_5_10.sce
@@ -0,0 +1,12 @@
+clear;
+clc;
+format('v',12);
+disp(0,'Point(3,-2,2) is in conductor region hence E=D=');
+ps=2;
+Dn=ps;
+D=[0 Dn 0];
+e=(10^-9)/(36*%pi);
+er=2;
+E=D/(e*er);
+disp(D,'D=');
+disp(E,'E=');
diff --git a/42/CH5/EX5.2/sadiku_5_2.sce b/42/CH5/EX5.2/sadiku_5_2.sce
new file mode 100755
index 000000000..21ac4b7b2
--- /dev/null
+++ b/42/CH5/EX5.2/sadiku_5_2.sce
@@ -0,0 +1,10 @@
+clear;

+clc;

+format('v',12);

+ps=10^-7;

+u=2;

+w=0.1;

+t=5;

+I=ps*u*w;

+Q=I*t*10^9;

+disp(Q,'charge(in nC) collected in 5 sec=');
\ No newline at end of file
diff --git a/42/CH5/EX5.3/sadiku_5_3.sce b/42/CH5/EX5.3/sadiku_5_3.sce
new file mode 100755
index 000000000..c0c3cb9a9
--- /dev/null
+++ b/42/CH5/EX5.3/sadiku_5_3.sce
@@ -0,0 +1,18 @@
+clear;

+clc;

+format('v',12);

+n=10^29;

+e=-1.6*10^-19;

+pv=n*e;

+disp(pv*10^-6,'(a) pv(in MC/m3)=');

+sigma=5*10^7;

+E=10^-2;

+J=sigma*E;

+disp(J*10^-3,'(b) J(in kA/m2)=');

+S=(%pi*10^-6)/4;

+I=J*S;

+format('v',6);

+disp(I,'(c) I(in A)=');

+u=J/pv;

+format('v',12);

+disp(u,'(d) u(in m/s)=');
\ No newline at end of file
diff --git a/42/CH5/EX5.4/sadiku_5_4.sce b/42/CH5/EX5.4/sadiku_5_4.sce
new file mode 100755
index 000000000..47ac6c7c5
--- /dev/null
+++ b/42/CH5/EX5.4/sadiku_5_4.sce
@@ -0,0 +1,10 @@
+clear;

+clc;

+format('v',6);

+l=4;

+d=3;

+r=0.5;

+S=(d^2-(%pi*r^2))*10^-4;

+sigma=5*10^6;

+R=(l*10^6)/(sigma*S);

+disp(R,'R(in microohm)=');        
\ No newline at end of file
diff --git a/42/CH5/EX5.6/sadiku_5_6.sce b/42/CH5/EX5.6/sadiku_5_6.sce
new file mode 100755
index 000000000..9c7a1f946
--- /dev/null
+++ b/42/CH5/EX5.6/sadiku_5_6.sce
@@ -0,0 +1,18 @@
+clear;
+clc;
+format('v',6);
+e0=10^-9/(36*%pi);
+er=2.55;
+E=10^4;
+d=1.5*10^-3;
+D=e0*er*E*10^9;
+disp(D,'D(in nC/m^2)=');
+xe=1.55;
+P=xe*e0*E*10^9;
+disp(P,'P(in nC/m^2)=');
+ps=D;
+disp(ps,'ps(in nC/m^2)=');
+pps=P;
+disp(pps,'pps(in nC/m^2)=');
+V=E*d;
+disp(V,'V(in V)=');
\ No newline at end of file
diff --git a/42/CH5/EX5.7/sadiku_5_7.sce b/42/CH5/EX5.7/sadiku_5_7.sce
new file mode 100755
index 000000000..38d11614b
--- /dev/null
+++ b/42/CH5/EX5.7/sadiku_5_7.sce
@@ -0,0 +1,15 @@
+clear;
+clc;
+format('v',6);
+Q=2*10^-12;
+e0=(10^-9)/(36*%pi);
+er=5.7;
+xr=er-1;
+r=10^-1;
+E=Q*10^12/(4*%pi*e0*er*r^2);
+P=xr*e0*E;
+pps=P*1;
+disp(pps,'(a) pps(in pC/m^2)=');
+Q1=-4*10^-12;
+F=(Q*Q1)*10^12/(4*%pi*e0*er*r^2);
+disp(F,'(b) F(in pN)(in the direction of ar)=');
\ No newline at end of file
diff --git a/42/CH5/EX5.9/sadiku_5_9.sce b/42/CH5/EX5.9/sadiku_5_9.sce
new file mode 100755
index 000000000..5e10dd01e
--- /dev/null
+++ b/42/CH5/EX5.9/sadiku_5_9.sce
@@ -0,0 +1,28 @@
+clear;
+clc;
+format('v',6);
+an=[0 0 1];
+E1=[5 -2 3];
+er1=4;
+er2=3;
+e=(10^-9)/(36*%pi);
+e1n=E1*an';
+E1n=[0 0 e1n];
+E2n=[0 0 E1n*[0;0;1]];
+E1t=E1-E1n;
+E2t=E1t;
+E2n=(er1*E1n)/er2;
+E2=E2t+E2n;
+disp(E2,'E2=');
+theta1=atand(((E1t(1,1)^2+E1t(1,2)^2+E1t(1,3)^2)^0.5)/e1n);
+alpha1=90-theta1;
+disp(alpha1,'Angle of E1 with interface=');
+alpha2=90-atand(((E2t(1,1)^2+E2t(1,2)^2+E2t(1,3)^2)^0.5)/((E2n(1,1)^2+E2n(1,2)^2+E2n(1,3)^2)^0.5));
+disp(alpha2,'Angle of E2 with interface=');
+wE1=0.5*er1*e*10^12*(E1(1,1)^2+E1(1,2)^2+E1(1,3)^2);
+wE2=0.5*er2*e*10^12*(E2(1,1)^2+E2(1,2)^2+E2(1,3)^2);
+disp(wE1,'Energy densities are wE1(in uJ)=');
+disp(wE2,'                     wE2(in uJ)=');
+We=wE2*integrate('1','x',2,4)*integrate('1','y',3,5)*integrate('1','z',-6,-4)*10^-3;
+disp(We,'We(in mJ)=');
+