initial commit / add all books

74 files changed, 904 insertions, 0 deletions

diff --git a/42/CH1/EX1.1/sadiku_1_1.sce b/42/CH1/EX1.1/sadiku_1_1.sce
new file mode 100755
index 000000000..d53165a9a
--- /dev/null
+++ b/42/CH1/EX1.1/sadiku_1_1.sce
@@ -0,0 +1,13 @@
+clear;

+clc;

+format('v',6)

+A=[10,-4,6];

+B=[2,1,0];

+disp(A(1,2),'Component of A along ay : ')

+P=3*A-B;

+disp((P(1,1)^2+P(1,2)^2+P(1,3)^2)^0.5,'magnitude is :')

+C=A+2*B;

+det_C=(C(1,1)^2+C(1,2)^2+C(1,3)^2)^0.5;

+format('v',7)

+ac=C/det_C;

+disp(ac,'Unit Vector along C is :')
\ No newline at end of file
diff --git a/42/CH1/EX1.2/sadiku_1_2.sce b/42/CH1/EX1.2/sadiku_1_2.sce
new file mode 100755
index 000000000..1a76ecaf0
--- /dev/null
+++ b/42/CH1/EX1.2/sadiku_1_2.sce
@@ -0,0 +1,14 @@
+clear;

+clc;

+format('v',6);

+P=[0,2,4];

+Q=[-3,1,5];

+origin=[0,0,0];

+rp=P-origin;

+disp(rp,'Position Vector of P is :')

+rpq=Q-P;

+disp(rpq,'Position Vector from P to Q is :')

+det_rpq=(rpq(1,1)^2+rpq(1,2)^2+rpq(1,3)^2)^0.5;

+disp(det_rpq,'distance between P and Q is :')

+A=10*rpq/det_rpq;

+disp([A;-A],'Vectors parallel to PQ with magnitude of 10 :')
\ No newline at end of file
diff --git a/42/CH1/EX1.3/sadiku_1_3.sce b/42/CH1/EX1.3/sadiku_1_3.sce
new file mode 100755
index 000000000..7c24e6765
--- /dev/null
+++ b/42/CH1/EX1.3/sadiku_1_3.sce
@@ -0,0 +1,9 @@
+format ('v',6);

+vb= [10*cos(%pi/4), -10*sin(%pi/4)]

+vm= [-2*cos(%pi/4), -2*sin(%pi/4)]

+vmg= vb+vm;

+disp (vmg, 'Velocity of man with respect to ground:')

+mod_vmg=(vmg(1,1)^2+vmg(1,2)^2)^.5;

+dir= atand(vmg(1,2)/vmg(1,1))

+disp( mod_vmg,'Absolute velocity of man is:' )

+disp (dir,'Angle with east in radian:')
\ No newline at end of file
diff --git a/42/CH1/EX1.4/sadiku_1_4.sce b/42/CH1/EX1.4/sadiku_1_4.sce
new file mode 100755
index 000000000..9b005ac42
--- /dev/null
+++ b/42/CH1/EX1.4/sadiku_1_4.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+A=[3,4,1];

+B=[0,2,-5];

+det_A=(A(1,1)^2+A(1,2)^2+A(1,3)^2)^0.5;

+det_B=(B(1,1)^2+B(1,2)^2+B(1,3)^2)^0.5;

+theta=acosd((sum(A.*B))/(det_A*det_B));

+disp(theta,'Angle between A and B is :')
\ No newline at end of file
diff --git a/42/CH1/EX1.5/sadiku_1_5.sce b/42/CH1/EX1.5/sadiku_1_5.sce
new file mode 100755
index 000000000..9e4a71d0f
--- /dev/null
+++ b/42/CH1/EX1.5/sadiku_1_5.sce
@@ -0,0 +1,39 @@
+clear;

+clc;

+format('v',7);

+P=[2,0,-1];

+Q=[2,-1,2];

+R=[2,-3,1];

+S=P+Q;

+T=P-Q;

+U1=S(1,2)*T(1,3)-S(1,3)*T(1,2);

+U2=S(1,3)*T(1,1)-S(1,1)*T(1,3);

+U3=S(1,1)*T(1,2)-S(1,2)*T(1,1);

+U=[U1 U2 U3];

+disp(U,'(P+Q)*(P-Q)=')

+V1=R(1,2)*P(1,3)-R(1,3)*P(1,2);

+V2=R(1,3)*P(1,1)-R(1,1)*P(1,3);

+V3=R(1,1)*P(1,2)-R(1,2)*P(1,1);

+V=[V1 V2 V3];

+X=(Q(1,1)*V(1,1)+Q(1,2)*V(1,2)+Q(1,3)*V(1,3));

+disp(X,'Q.R*P')

+W1=Q(1,2)*R(1,3)-Q(1,3)*R(1,2);

+W2=Q(1,3)*R(1,1)-Q(1,1)*R(1,3);

+W3=Q(1,1)*R(1,2)-Q(1,2)*R(1,1);

+W=[W1 W2 W3];

+Y=(W(1,1)*P(1,1)+W(1,2)*P(1,2)+W(1,3)*P(1,3));

+disp(Y,'P.Q*R')

+det_W=(W(1,1)^2+W(1,2)^2+W(1,3)^2)^.5;

+det_Q=(Q(1,1)^2+Q(1,2)^2+Q(1,3)^2)^.5;

+det_R=(R(1,1)^2+R(1,2)^2+R(1,3)^2)^.5

+sineoftheta=(det_W/(det_Q*det_R));

+disp(sineoftheta,'sin of theta=')

+Z1=P(1,2)*W(1,3)-P(1,3)*W(1,2);

+Z2=P(1,3)*W(1,1)-P(1,1)*W(1,3);

+Z3=P(1,1)*W(1,2)-P(1,2)*W(1,1);

+Z=[Z1 Z2 Z3];

+disp(Z,'P* Q*R=')

+disp(W/det_W,'Unit Vector Perpendicular to Q & R')

+q=Q/det_Q;

+C=(P(1,1)*q(1,1)+P(1,2)*q(1,2)+P(1,3)*q(1,3));

+disp(C*q,'Component of P along Q');

diff --git a/42/CH1/EX1.7/sadiku_1_7.sce b/42/CH1/EX1.7/sadiku_1_7.sce
new file mode 100755
index 000000000..2da9cf668
--- /dev/null
+++ b/42/CH1/EX1.7/sadiku_1_7.sce
@@ -0,0 +1,25 @@
+clear;

+clc;

+format('v',6);

+P1=[5 2 -4];

+P2=[1 1 2];

+P3=[-3 0 8];

+P4=[3 -1 0]

+R1=P1-P2;

+R2=P1-P3;

+R3=P2-P3;

+R4=P1-P4;

+U1=R1(1,2)*R2(1,3)-R1(1,3)*R2(1,2);

+U2=R1(1,3)*R2(1,1)-R1(1,1)*R2(1,3);

+U3=R1(1,1)*R2(1,2)-R1(1,2)*R2(1,1);

+U=[U1 U2 U3];

+disp(U)

+disp('Since U is Zero so P1,P2,P3 are in straight line')

+det_R1=(R1(1,1)^2+R1(1,2)^2+R1(1,3)^2)^.5;

+V1=R4(1,2)*R1(1,3)-R4(1,3)*R1(1,2);

+V2=R4(1,3)*R1(1,1)-R4(1,1)*R1(1,3);

+V3=R4(1,1)*R1(1,2)-R4(1,2)*R1(1,1);

+V=[V1 V2 V3];

+det_V=(V(1,1)^2+V(1,2)^2+V(1,3)^2)^.5;

+det_R1=(R1(1,1)^2+R1(1,2)^2+R1(1,3)^2)^.5;

+disp((det_V/det_R1),'Shortest Distance')
\ No newline at end of file
diff --git a/42/CH10/EX10.1/sadiku_10_1.sce b/42/CH10/EX10.1/sadiku_10_1.sce
new file mode 100755
index 000000000..9d22d6b91
--- /dev/null
+++ b/42/CH10/EX10.1/sadiku_10_1.sce
@@ -0,0 +1,22 @@
+clear;

+clc;

+format('v',6);

+disp('Direction of wave propagation is -ax');

+w=10^8,c=3*10^8;

+B=w/c;

+disp(B,'Value of beta=');

+T=2*%pi/w;

+disp(T/2*10^9,'Time taken to travel half of wave length in nS= ');

+t=0

+x=-2*%pi:%pi/16:2*%pi;

+Ey=50*cos(10^8 *t +B*x);

+subplot(2,2,1)

+plot(x,Ey);

+t=T/4;

+Ey=50*cos(10^8 *t +B*x);

+subplot(2,2,2)

+plot(x,Ey);

+t=T/2;

+Ey=50*cos(10^8 *t +B*x);

+subplot(2,2,3)

+plot(x,Ey);
\ No newline at end of file
diff --git a/42/CH10/EX10.10/sadiku_10_10.sce b/42/CH10/EX10.10/sadiku_10_10.sce
new file mode 100755
index 000000000..0951436be
--- /dev/null
+++ b/42/CH10/EX10.10/sadiku_10_10.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+kx=0,ky=.866,kz=.5,Eo=10^-9 /(36*%pi),Uo=4*%pi*10^-7;

+k=sqrt(kx*kx+ky*ky+kz*kz);

+w=k/(sqrt(Uo*Eo));

+disp(w*10^-6,'w im Mrad/sec');

+l=2*%pi/k;

+disp(l,'lamda = ')

diff --git a/42/CH10/EX10.2/sadiku_10_2.sce b/42/CH10/EX10.2/sadiku_10_2.sce
new file mode 100755
index 000000000..fe47ad353
--- /dev/null
+++ b/42/CH10/EX10.2/sadiku_10_2.sce
@@ -0,0 +1,9 @@
+clear;

+clc;

+Ho=10,n=200*%e^(%i*%pi/6),P=atan(3^.5),b=.5,e=10^-9 /(36*%pi);

+Eo=n*Ho;

+disp('a=w*sqrt(u*e/2*(1+(c/(w*e)^2)^.5)-1)');

+disp('b=w*sqrt(u*e/2*(1+(c/(w*e)^2)^.5)+1)');

+a=b*((sqrt(((1+(tan(P))^2)^.5)-1))/(sqrt(((1+(tan(P))^2)^.5)+1)));

+disp(a,'Value of alpha=');

+disp(1/a,'Skin depth = ')
\ No newline at end of file
diff --git a/42/CH10/EX10.3/sadiku_10_3.sce b/42/CH10/EX10.3/sadiku_10_3.sce
new file mode 100755
index 000000000..3e063eb1b
--- /dev/null
+++ b/42/CH10/EX10.3/sadiku_10_3.sce
@@ -0,0 +1,7 @@
+clear;

+clc;

+B=1,n=60*%pi,Ur=1,Eo=10^-9 /(36*%pi),Uo=4*%pi*10^-7;

+Er=Uo*Ur/(n^2 *Eo);

+disp(Er,'Er =');

+w=B/sqrt(Eo*Er*Uo*Ur);

+disp(w*10^-6, 'w in Mrad/sec');
\ No newline at end of file
diff --git a/42/CH10/EX10.4/sadiku_10_4.sce b/42/CH10/EX10.4/sadiku_10_4.sce
new file mode 100755
index 000000000..d1844c6df
--- /dev/null
+++ b/42/CH10/EX10.4/sadiku_10_4.sce
@@ -0,0 +1,5 @@
+clear;

+clc;

+c=3,w=10^8,Ur=20,Eo=10^-9 /(36*%pi),Er=1,Uo=4*%pi*10^-7;

+a=sqrt(Uo*Ur*w*c/2);

+disp(a,'alpha = beta =');//as c/w*E>>1

diff --git a/42/CH10/EX10.6/sadiku_10_6.sce b/42/CH10/EX10.6/sadiku_10_6.sce
new file mode 100755
index 000000000..e7fc6ab92
--- /dev/null
+++ b/42/CH10/EX10.6/sadiku_10_6.sce
@@ -0,0 +1,7 @@
+clear;

+clc;

+a=2*10^-3,b=6*10^-3,t=10^-3,l=2,c=5.8*10^7;

+Ri=l/(c*%pi*a*a);

+Ro=l/(c*%pi*((b+t)^2-b^2));

+Rdc=Ro+Ri;

+disp(Rdc*10^3,'Resistance in mOhm');
\ No newline at end of file
diff --git a/42/CH10/EX10.7/sadiku_10_7.sce b/42/CH10/EX10.7/sadiku_10_7.sce
new file mode 100755
index 000000000..1bcfe4b6f
--- /dev/null
+++ b/42/CH10/EX10.7/sadiku_10_7.sce
@@ -0,0 +1,7 @@
+clear;

+clc;

+a=0,b=.8,Eo=10^-9 /(36*%pi),Uo=4*%pi*10^-7,Ur=1,w=2*%pi*10^7;

+Er=b^2/(Uo*Eo*w*w);

+disp(Er);

+n=sqrt(Uo/(Eo*Er));

+disp(n);
\ No newline at end of file
diff --git a/42/CH11/EX11.1/sadiku_11_1.sce b/42/CH11/EX11.1/sadiku_11_1.sce
new file mode 100755
index 000000000..e33e58c03
--- /dev/null
+++ b/42/CH11/EX11.1/sadiku_11_1.sce
@@ -0,0 +1,9 @@
+clear;

+clc;

+format('v',6);

+R=0,G=0,a=0,Ro=70,B=3,f=100*10^6;

+w=2*%pi*f;

+C=B/(w*Ro);

+disp(C*10^12,'Capacitance per meter of line in pF')

+L=Ro*Ro*C;

+disp(L*10^9,'Inductance per meter in nHz')
\ No newline at end of file
diff --git a/42/CH11/EX11.10/sadiku_11_10.sce b/42/CH11/EX11.10/sadiku_11_10.sce
new file mode 100755
index 000000000..379ba54b4
--- /dev/null
+++ b/42/CH11/EX11.10/sadiku_11_10.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+format('v',6);

+Er=3.8, c=3*10^8;

+r=4.5;//ratio w/h

+Eeff= ((Er+1)/2)+ ((Er-1)/(2*(1+12/r)^.5));

+disp(Eeff, 'The effective relative permittivity = ');

+Zo=(120*%pi)/((r+1.393+ (.667*log(r+1.444)))*((Eeff)^.5));

+disp(Zo, 'Character impedence of line');

+f=10^10;

+l=c/(f*sqrt(Eeff));

+disp(l*1000, 'The wavelength of line at 10 GHz');

diff --git a/42/CH11/EX11.11/sadiku_11_11.sce b/42/CH11/EX11.11/sadiku_11_11.sce
new file mode 100755
index 000000000..325e163c3
--- /dev/null
+++ b/42/CH11/EX11.11/sadiku_11_11.sce
@@ -0,0 +1,13 @@
+clear;

+clc;

+h=1, w=.8, Er=6.6, P= atan(.0001), c= 5.8*10^7,f=10^10,mu=4*%pi*10^(-7),C=3*10^8;

+r=w/h;

+Ee=((Er+1)/2)+ ((Er-1)/(2*(1+12/r)^.5));

+Zo=(120*%pi)/((r+1.393+ (.667*log(r+1.444)))*((Ee)^.5));

+Rs=((%pi*f*mu)/c)^.5;

+ac=8.686*Rs/(w*10^-3*Zo);

+disp(ac,'Conduction Attenuation Constant = ');

+l=C/(f*(Ee)^.5);

+disp(l);

+ad=27.3*(Ee-1)*Er*tan(P)/((Er-1)*Ee*l);

+disp(ad,'Dielectric Attenuation Constant =');
\ No newline at end of file
diff --git a/42/CH11/EX11.2/sadiku_11_2.sce b/42/CH11/EX11.2/sadiku_11_2.sce
new file mode 100755
index 000000000..c4c50b012
--- /dev/null
+++ b/42/CH11/EX11.2/sadiku_11_2.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+Zo=60,a=20*10^-3,u=.6*3*10^8, f=100*10^6;

+R=a*Zo,disp(R,'R=');

+L=Zo/u,disp(L*10^9,'L in nH=');

+G=a*a/R,disp(G*10^6,'G in micro S per meter =');

+C=1/(u*Zo),disp(C*10^12,'C in pF =');

+l=u/f;disp(l,'l=');

diff --git a/42/CH11/EX11.3/sadiku_11_3.sce b/42/CH11/EX11.3/sadiku_11_3.sce
new file mode 100755
index 000000000..9046fa37a
--- /dev/null
+++ b/42/CH11/EX11.3/sadiku_11_3.sce
@@ -0,0 +1,18 @@
+clear;

+clc;

+format('v',6);

+w=10^6,B=1,a=8,Vg=10;

+Zo=60+40*%i,Zg=40,Zl=20+50*%i;

+a=(a/8.686);;//Since 1Np=8.686 dB

+Y=a+B*%i;

+Yl=2*Y;

+h=tanh(Yl);

+Zin=Zo*(Zl+Zo*tanh(Yl))/(Zo+Zl*tanh(Yl));

+disp(Zin,'The input impdence =');

+Io=Vg/(Zin+Zg);//at  z=0

+disp(Io*1000,'Sending end current in mA =');

+Vo=Zin*Io;

+Vop =(Vo+Zo*Io)/2;

+Vom =(Vo-Zo*Io)/2;

+Im= ((Vop * %e^-Y)/Zo)- ((Vom * %e^Y)/Zo);

+disp(Im*1000, 'Current at middle line in mA= ');

diff --git a/42/CH11/EX11.4/sadiku_11_4.sce b/42/CH11/EX11.4/sadiku_11_4.sce
new file mode 100755
index 000000000..bf1c4c41b
--- /dev/null
+++ b/42/CH11/EX11.4/sadiku_11_4.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+format('v',6);

+l=30,Zo=50,f=2*10^6,Zl=60+40*%i,u=.6*3*10^8;

+w=2*%pi*f;

+T=(Zl-Zo)/(Zl+Zo);

+disp(T,'Reflection coefficient  =');

+s=(1+abs(T))/(1-abs(T));

+disp(s,'Standing wave ratio =');

+B=w/u;disp(B*l);

+Zin=Zo*(Zl+Zo*tan(B*l)*%i)/(Zo+Zl*tan(B*l)*%i);

+disp(Zin);
\ No newline at end of file
diff --git a/42/CH11/EX11.5/sadiku_11_5.sce b/42/CH11/EX11.5/sadiku_11_5.sce
new file mode 100755
index 000000000..0d3ddb6a8
--- /dev/null
+++ b/42/CH11/EX11.5/sadiku_11_5.sce
@@ -0,0 +1,23 @@
+clear;

+clc;

+format('v',6);

+Zl=100+150*%i;

+Zo=75;

+zl=Zl/Zo;

+T=(Zl-Zo)/(Zl+Zo);

+disp(T,'T =');

+s=(1+abs(T))/(1-abs(T));

+disp(s,'s =')

+format('v',5);

+Yl=1/Zl;

+disp(Yl*1000,'Load admittance in mS');

+B=2*%pi,l=.4;

+Zin=Zo*(Zl+Zo*tan(B*l)*%i)/(Zo+Zl*tan(B*l)*%i);

+format('v',6);

+disp(Zin,'Zin at .4 l from load')//for .4l

+B=2*%pi,l=.6;

+Zin=Zo*(Zl+Zo*tan(B*l)*%i)/(Zo+Zl*tan(B*l)*%i);

+format('v',6);

+disp(Zin,'Zin at .6 l from load')//for .6l

+

+

diff --git a/42/CH11/EX11.6/sadiku_11_6.sce b/42/CH11/EX11.6/sadiku_11_6.sce
new file mode 100755
index 000000000..e1fb6d8c9
--- /dev/null
+++ b/42/CH11/EX11.6/sadiku_11_6.sce
@@ -0,0 +1,11 @@
+clear;

+clc;

+s=2, l1=11,l2=19,ma=24,mi=16,u=3*10^8,Zo=50;

+l=(l2-l1)*2;

+disp(l,'Lamda =');

+f=u/l;

+disp(f*10^-6,'Frequency im MHz =');

+L=(24-19)/l;//Let us assume load is at 24cm

+zl=1.4+.75*%i; //by smith chart

+Zl=Zo*zl;

+disp(Zl,'Zl =')
\ No newline at end of file
diff --git a/42/CH11/EX11.7/sadiku_11_7.sce b/42/CH11/EX11.7/sadiku_11_7.sce
new file mode 100755
index 000000000..f5be14ec6
--- /dev/null
+++ b/42/CH11/EX11.7/sadiku_11_7.sce
@@ -0,0 +1,14 @@
+clear;

+clc;

+format('v',6);

+Zo=100, Zl=40+30*%i;

+Yo=1/Zo;

+yl=Zo/Zl;

+ys1=1.04*%i, ys2=-1.04*%i; //By smith chart

+Ys1=Yo*ys1, Ys2=Yo*ys2;

+disp([Ys1*1000 Ys2*1000],'Possible values of sub admittance in mS =');

+la=.5 - (62-(-39))/720 ;disp(la, 'distance between load and antenna at A devided by Lamda');

+lb= (62-39)/720;disp(lb, 'distance between load and antenna at B devided by Lamda');//With the help of figure

+da=88/720, db= 272/720;

+format('v',7);

+disp(da,db,'Sub length devided by Lamda');
\ No newline at end of file
diff --git a/42/CH11/EX11.8/sadiku_11_8.sce b/42/CH11/EX11.8/sadiku_11_8.sce
new file mode 100755
index 000000000..9ad5bb7f7
--- /dev/null
+++ b/42/CH11/EX11.8/sadiku_11_8.sce
@@ -0,0 +1,6 @@
+clear;

+clc;

+Zg=100,Zo=50,Zl=200,u=3*10^8,l=100,Vg=12;

+Tg=(Zg-Zo)/(Zg+Zo);

+Tl=(Zl-Zo)/(Zl+Zo);

+t1=l/u;

diff --git a/42/CH12/EX12.1/sadiku_12_1.sce b/42/CH12/EX12.1/sadiku_12_1.sce
new file mode 100755
index 000000000..d6a6943b1
--- /dev/null
+++ b/42/CH12/EX12.1/sadiku_12_1.sce
@@ -0,0 +1,30 @@
+clear;

+clc;

+a=2.5*10^-2, b=1*10^-2,c=0, Ur=1,Er=4,C=3*10^8;

+fc=0,m=0,n=0;

+while(fc*10^-9<15.1)

+fc=(C/(4*a))*sqrt(m^2 + (a*n/b)^2);

+if ((fc*10^-9) < 15.1) then

+n=n+1;

+else disp(n-1,'Max value of n is ='); end

+end 

+fc=0,m=0,n=0;

+while(fc*10^-9<15.1)

+fc=(C/(4*a))*sqrt(m^2 + (a*n/b)^2);

+if ((fc*10^-9) < 15.1) then

+m=m+1;

+else disp(m-1,'Max value of m is ='); end

+end 

+function[p]= modes(m,n);

+p=(C/(4*a))*sqrt(m^2 + (a*n/b)^2);

+if ((p*10^-9) < 15.1) then

+disp([m n],'Transmission mode is possible'); else p=0;end

+endfunction

+for i=1:1:5, for j=1:1:2,modes(i,j);end;

+end

+    

+

+

+

+

+

diff --git a/42/CH12/EX12.3/sadiku_12_3.sce b/42/CH12/EX12.3/sadiku_12_3.sce
new file mode 100755
index 000000000..1f8146b2e
--- /dev/null
+++ b/42/CH12/EX12.3/sadiku_12_3.sce
@@ -0,0 +1,13 @@
+clear;

+clc;

+format('v',7);

+a=1.5*10^-2,b=.8*10^-2,c=0,Uo=4*%pi*10^-7,Ur=1,Eo=10^-9/(36*%pi),Er=4,C=3*10^8,w=%pi*10^11,m=1,n=3;

+u=C/2;

+f=w/(2*%pi);

+fc=u*((m*m)/(a*a) + (n*n)/(b*b) )^.5/2;

+disp(fc*10^-9,'Cutoff frequency = ');

+B=w*sqrt(Uo*Ur*Eo*Er)*sqrt(1-(fc/f)^2);

+disp(B,'Phase constant = ');

+disp(%i*B,'Propagation constant = ');

+n=377/sqrt(Er)*sqrt(1-(fc/f)^2);

+disp(n,'Intrinsic wave impedence = ');
\ No newline at end of file
diff --git a/42/CH12/EX12.4/sadiku_12_4.sce b/42/CH12/EX12.4/sadiku_12_4.sce
new file mode 100755
index 000000000..8640072c1
--- /dev/null
+++ b/42/CH12/EX12.4/sadiku_12_4.sce
@@ -0,0 +1,13 @@
+clear;

+clc;

+a=8.636*10^-2,b=4.318*10^-2,f=4*10^9;

+u=3*10^8;

+fc=u/(2*a);

+disp(fc*10^-9,'Cut off frquency = ');

+if(f>fc) then disp('As f>fc so TE10 mode will propagate')

+else disp('It will not propagate') 

+end

+Up=u/sqrt(1-(fc/f)^2);

+disp(Up*10^-6,'Phase velocity in Mm/sec = ');

+Ug=u*u/Up;

+disp(Ug*10^-6,'Group velocity in Mm/sec = ');
\ No newline at end of file
diff --git a/42/CH12/EX12.5/sadiku_12_5.sce b/42/CH12/EX12.5/sadiku_12_5.sce
new file mode 100755
index 000000000..9d487e2a1
--- /dev/null
+++ b/42/CH12/EX12.5/sadiku_12_5.sce
@@ -0,0 +1,7 @@
+clear;

+clc;

+f=10*10^9,a=4*10^-2,b=2*10^-2,u=3*10^8,Pavg=2*10^-3;

+fc=u/(2*a);

+n=377/sqrt(1-(fc/f)^2);

+E=sqrt(4*n*Pavg/(a*b));

+disp(E,'Peak value of Electric field = ');
\ No newline at end of file
diff --git a/42/CH12/EX12.6/sadiku_12_6.sce b/42/CH12/EX12.6/sadiku_12_6.sce
new file mode 100755
index 000000000..ff01358cf
--- /dev/null
+++ b/42/CH12/EX12.6/sadiku_12_6.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+cc=5.8*10^7, f=4.8*10^9,c=10^-17,Uo=4*%pi*10^-7,Eo=10^-9/(36*%pi),Er=2.55,z=60*10^-2,l=4.2*10^-2,b=2.6*10^-2,P=1.2*10^3;

+n=377/sqrt(Er);

+u=3*10^8 /sqrt(Er);

+fc=u/(2*l);

+ad=c*n/(2*sqrt(1-(fc/f^2)));

+Rs=sqrt(%pi*f*Uo/cc);

+ac=2*Rs*(.5+(b/l)*(fc/f)^2)/(b*n*sqrt(1-(fc/f)^2));

+a=ac;

+Pd=P*(%e^(2*a*z) -1);

+disp(Pd,'Power dissipiated = ');
\ No newline at end of file
diff --git a/42/CH12/EX12.8/sadiku_12_8.sce b/42/CH12/EX12.8/sadiku_12_8.sce
new file mode 100755
index 000000000..462d63fda
--- /dev/null
+++ b/42/CH12/EX12.8/sadiku_12_8.sce
@@ -0,0 +1,7 @@
+clear;

+clc;

+a=5*10^-2,b=4*10^-2,c=10*10^-2,C=5.8*10^7,Uo=4*%pi*10^-7;

+f101=3.335*10^9;

+d=sqrt(1/(%pi*f101*Uo*C));

+Q=(a*a+c*c)*a*b*c/(d*(2*b*(a^3 + c^3)+a*c*(a*a+c*c)));

+disp(Q,'Quality factor of TE101 = ');
\ No newline at end of file
diff --git a/42/CH13/EX13.1/sadiku_13_1.sce b/42/CH13/EX13.1/sadiku_13_1.sce
new file mode 100755
index 000000000..81d97868d
--- /dev/null
+++ b/42/CH13/EX13.1/sadiku_13_1.sce
@@ -0,0 +1,31 @@
+clear;

+clc;

+format('v',5);

+function[P,I]=powerhert(H,P,r,B,dl)

+I=H*4*r*%pi/((B*(dl))*sin(P));    

+P=40*%pi*%pi*I*I*dl*dl;

+disp(P*1000,'Power transmit by Hertizian dipole in mWatt');

+endfunction

+powerhert((5*(10)^-6),%pi/2,2000,(2*%pi),1/25);

+function[P,I]=powerhw(H,P,r)

+I=H*2*r*%pi*sin(P)/(cos((%pi/2)*cos(P)));R=73;    

+P=(I*I*R)/2;

+disp(P*1000,'Power transmit by Half wave dipole in mWatt');

+endfunction

+powerhw((5*(10)^-6),%pi/2,2000);

+function[P,I]=powerqw(H,P,r)

+I=H*2*r*%pi*sin(P)/(cos((%pi/2)*cos(P)));R=36.56;    

+P=(I*I*R)/2;format('v',4);

+disp(P*1000,'Power transmit by Quarterwave monopole in mWatt');

+endfunction

+powerqw((5*(10)^-6),%pi/2,2000);

+function[P,I]=powersingloop(H,r,k);R=192.3;

+I=H*r/(%pi*%pi*10*k*k);

+P=(I*I*R)/2;

+disp(P*1000,'Power transmit by 10 turn loop antena in mWatt');

+endfunction

+powersingloop((5*(10)^-6),2000,1/20);

+

+

+

+

diff --git a/42/CH13/EX13.10/sadiku_13_10.sce b/42/CH13/EX13.10/sadiku_13_10.sce
new file mode 100755
index 000000000..6ab8bf7a9
--- /dev/null
+++ b/42/CH13/EX13.10/sadiku_13_10.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+c=3*(10)^8, f=3*(10)^9,Aet=9,r1=1.852*(10)^5,r2=4*r1,r3=5.556*10^5,Pr=200*(10)^3,a=20;

+l=c/f;

+Gdt=4*%pi*Aet/(l*l);

+P1=Gdt*Pr/(4*%pi*r1*r1);

+P2=Gdt*Pr/(4*%pi*r2*r2);

+disp(P1*1000,'Signal power density at 100nmile in mWatt');

+disp(P2*1000,'Signal power density at 400nmile in mWatt');

+Pr=Aet*a*Gdt*Pr/(4*%pi*r3*r3)^2;

+disp(Pr*10^12,'Power of reflected signal in picoWatt');

+

diff --git a/42/CH13/EX13.2/sadiku_13_2.sce b/42/CH13/EX13.2/sadiku_13_2.sce
new file mode 100755
index 000000000..b4c6cb804
--- /dev/null
+++ b/42/CH13/EX13.2/sadiku_13_2.sce
@@ -0,0 +1,18 @@
+clear;

+clc;

+format('v',6);

+c=3*10^8;

+f=50*10^6;

+disp(c/(2*f),'Length of halfdipole in meter');

+function[P,I]=curpow(E,P,r)

+n=120*%pi; R=73;

+I=E*2*r*%pi*sin(P)/(n*(cos((%pi/2)*cos(P))));   

+P=(I*I*R)/2;

+disp(I*1000,'Current fed to antenna in mA');

+disp(P*1000,'Power radiated by Antenna in mWatt');

+endfunction

+curpow((10*(10)^-6),%pi/2,500*10^3);

+Zl=73+42.5*%i,Zo=75;

+T=(Zl-Zo)/(Zl+Zo);

+s=(1+abs(T))/(1-abs(T));

+disp(s,'Standing wave ratio');
\ No newline at end of file
diff --git a/42/CH13/EX13.3/sadiku_13_3.sce b/42/CH13/EX13.3/sadiku_13_3.sce
new file mode 100755
index 000000000..4a4275e74
--- /dev/null
+++ b/42/CH13/EX13.3/sadiku_13_3.sce
@@ -0,0 +1,5 @@
+clear;

+clc;

+G=(integrate('(sin(P))^3','P',0,%pi))*integrate('1','Q',0,2*%pi)

+Gd=4*%pi/G;

+disp(Gd)

diff --git a/42/CH13/EX13.4/sadiku_13_4.sce b/42/CH13/EX13.4/sadiku_13_4.sce
new file mode 100755
index 000000000..c32f27ba0
--- /dev/null
+++ b/42/CH13/EX13.4/sadiku_13_4.sce
@@ -0,0 +1,10 @@
+clear;

+clc;

+format('v',7);

+G=5;

+r=10*10^3;

+P=20*10^3;

+n=120*%pi;

+Gd=10^(G/10);

+E=sqrt(n*Gd*P/(2*%pi*r*r));

+disp(E,'Electric field intensity at 10 km =');
\ No newline at end of file
diff --git a/42/CH13/EX13.5/sadiku_13_5.sce b/42/CH13/EX13.5/sadiku_13_5.sce
new file mode 100755
index 000000000..239b77680
--- /dev/null
+++ b/42/CH13/EX13.5/sadiku_13_5.sce
@@ -0,0 +1,6 @@
+clear;

+clc;

+Umax=2;

+Uavg=(1/(4*%pi))*2*integrate('(sin (P))^2','P',0,%pi)*integrate('(sin (Q))^3','Q',0,%pi);

+D=Umax/Uavg;

+disp(D,'Directivity of antenna');
\ No newline at end of file
diff --git a/42/CH13/EX13.8/sadiku_13_8.sce b/42/CH13/EX13.8/sadiku_13_8.sce
new file mode 100755
index 000000000..b2b16f08f
--- /dev/null
+++ b/42/CH13/EX13.8/sadiku_13_8.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+c=3*10^8,f=30*10^6,E=2*10^-3;

+l=c/f;

+n=120*%pi,R=73;

+format('v',5);

+Gdmax=n/(%pi*R);

+format('v',6);

+Amax=(l^2 /(4*%pi))*Gdmax;

+disp(Amax,'Maximum effective area');

+Pr=(E*E*Amax)/(2*n);

+disp(Pr*(10^9),'Power rerceived in nWatt')
\ No newline at end of file
diff --git a/42/CH13/EX13.9/sadiku_13_9.sce b/42/CH13/EX13.9/sadiku_13_9.sce
new file mode 100755
index 000000000..19daa2929
--- /dev/null
+++ b/42/CH13/EX13.9/sadiku_13_9.sce
@@ -0,0 +1,6 @@
+clear;

+clc;

+Gt=25,Gr=18,r=200,Pr=5*10^-3 ;

+Gdt=10^(Gt/10),Gdr=10^(Gr/10);

+Pt=Pr*(4*%pi*r)^2 /(Gdr*Gdt);

+disp(Pt,'Minimum power received in Watt =');
\ No newline at end of file
diff --git a/42/CH14/EX14.1/sadiku_14_1.sce b/42/CH14/EX14.1/sadiku_14_1.sce
new file mode 100755
index 000000000..4b5cc8994
--- /dev/null
+++ b/42/CH14/EX14.1/sadiku_14_1.sce
@@ -0,0 +1,10 @@
+clear;

+clc;

+S11=.85*(cosd(-30)+%i*sind(-30));

+S12=.07*(cosd(56)+%i*sind(56));

+S21=1.68*(cosd(120)+%i*sind(120));

+S22=.85*(cosd(-40)+%i*sind(-40));

+Zl=75,Zo=75;

+Tl=(Zl-Zo)/(Zl+Zo);

+Ti=S11+ (S12*S21*Tl)/(1-S22*Tl);

+disp(Ti,'Input reflection coefficient=')
\ No newline at end of file
diff --git a/42/CH14/EX14.2/sadiku_14_2.sce b/42/CH14/EX14.2/sadiku_14_2.sce
new file mode 100755
index 000000000..d80bdc0fa
--- /dev/null
+++ b/42/CH14/EX14.2/sadiku_14_2.sce
@@ -0,0 +1,11 @@
+clear;

+clc;

+format('v',6)

+d=80*(10)^-6;

+n1=1.62,NA=.21,L=8*(10)^-7 ;

+P=asind(NA);

+disp(P,'Acceptance angle');

+n2=sqrt(n1^2 - NA^2);

+disp(n2,'Refractive index');

+V=(%pi*d/L)*sqrt(n1^2 - n2^2);

+disp(V,'No of modes');

diff --git a/42/CH14/EX14.3/sadiku_14_3.sce b/42/CH14/EX14.3/sadiku_14_3.sce
new file mode 100755
index 000000000..6b3a9caab
--- /dev/null
+++ b/42/CH14/EX14.3/sadiku_14_3.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+a=.25;

+P=1-.4;

+l=(10/a)*log10(1/P);

+disp(l,'Distance travelled in Km');

+

+

diff --git a/42/CH2/EX2.1/sadiku_2_1.sce b/42/CH2/EX2.1/sadiku_2_1.sce
new file mode 100755
index 000000000..fa4b595b5
--- /dev/null
+++ b/42/CH2/EX2.1/sadiku_2_1.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+format('v',7);

+x=-2;y=6;z=3;

+r=(x^2+y^2)^.5;

+B=atand(y/x);

+R=sqrt(x^2+y^2+z^2);

+X=atand(r/z);

+disp([r B z ],'Cylindrical acordinate of P:');

+disp([R X B],'Spherical Cordinate of P:');

+A=[cosd(B) sind(B) 0;-sind(B) cosd(B) 0;0 0 1]*[y;x+z;0];

+disp (A,'A in cylindrical cordinates')
\ No newline at end of file
diff --git a/42/CH2/EX2.2/sadiku_2_2.sce b/42/CH2/EX2.2/sadiku_2_2.sce
new file mode 100755
index 000000000..570296bc9
--- /dev/null
+++ b/42/CH2/EX2.2/sadiku_2_2.sce
@@ -0,0 +1,22 @@
+clear;

+clc;

+format('v',6);

+function [X,Y,Z]=sptocart(x,y,z);

+R=sqrt(x^2+y^2+z^2);r=sqrt(x^2+y^2);

+P=asin(r/R);Q=acos(x/r);

+X=(10/R)*sin(P)*cos(Q)+R*(cos(P))^2 *cos(Q)-sin(Q);

+Y=(10/R)*sin(P)*sin(Q)+R*(cos(P))^2 *sin(Q)+cos(Q);

+Z=(10/R)*cos(P)-R*cos(P)*sin(P);

+disp([X Y Z],'B in cartesian cordinate')

+endfunction

+sptocart(-3,4,0);

+function [r,p,z]=sptocylin(r1,p1,z1);

+R=sqrt(r1^2+z1^2);

+P=acos(z1/R);

+r=(10/R)*sin(P)+R*(cos(P))^2 ;

+p=1;

+z=(10/R)*cos(P)-R*cos(P)*sin(P);

+disp([r p z],'B in cylindrical cordinates');

+endfunction

+sptocylin(5,%pi/2,-2);

+

diff --git a/42/CH2/EX2.3/sadiku_2_3.sce b/42/CH2/EX2.3/sadiku_2_3.sce
new file mode 100755
index 000000000..34676fabb
--- /dev/null
+++ b/42/CH2/EX2.3/sadiku_2_3.sce
@@ -0,0 +1,15 @@
+clear;

+clc;

+E=[-5 10 3];ModE=sqrt((-5)^2+10^2+3^2);

+F=[1 2 -6];

+P=[5,%pi/2,3];

+G1=E(1,2)*F(1,3)-E(1,3)*F(1,2);

+G2=E(1,3)*F(1,1)-E(1,1)*F(1,3);

+G3=E(1,1)*F(1,2)-E(1,2)*F(1,1);

+G=[G1 G2 G3];

+disp(sqrt(G1^2+G2^2+G3^2),'Mod of (E*F)');

+ay=[sin(%pi/2) cos(%pi/2) 0];

+Ey=(E(1,1)*ay(1,1)+E(1,2)*ay(1,2)+E(1,3)*ay(1,3));

+disp(Ey,'Component of E parallel to x=2 & z=3');

+P=acosd(3/ModE);

+disp(90-P,'Angle which make E wid Z=3');
\ No newline at end of file
diff --git a/42/CH2/EX2.4/sadiku_2_4.sce b/42/CH2/EX2.4/sadiku_2_4.sce
new file mode 100755
index 000000000..ff0403987
--- /dev/null
+++ b/42/CH2/EX2.4/sadiku_2_4.sce
@@ -0,0 +1,21 @@
+clear;

+clc;

+format('v',6)

+function [R,P,Q]=Posvec(r,p,q);

+R=r*sind(q);P=-sind(p)*cosd(q)/r;Q=r*r;

+D=[R P Q];

+disp(D,'D at P');

+Dn=[r*sind(q) 0 0];

+Dt=D-Dn;

+disp(Dt, 'Tangential component of D at P');

+endfunction

+Posvec(10,150,330);

+D=[-5 .043 100];

+a=[0 1 0];

+U1=D(1,2)*a(1,3)-D(1,3)*a(1,2);

+U2=D(1,3)*a(1,1)-D(1,1)*a(1,3);

+U3=D(1,1)*a(1,2)-D(1,2)*a(1,1);

+U=[U1 U2 U3];

+det_U=sqrt(U1^2+U2^2+U3^2);

+format('v',7);

+disp(U/det_U,'Unit vector P perpendicular to D');

diff --git a/42/CH3/EX3.1/sadiku_3_1.sce b/42/CH3/EX3.1/sadiku_3_1.sce
new file mode 100755
index 000000000..741cb0b8c
--- /dev/null
+++ b/42/CH3/EX3.1/sadiku_3_1.sce
@@ -0,0 +1,7 @@
+clear;

+clc;

+B1=[0,5,0],B2=[5,%pi/2,0],C1=[0 5 10],C2=[5 %pi/2 0],D1=[5 0 10],D2=[5,0,10],p=5;

+BC=integrate('1','Z',0,10);

+disp(BC);   //as dl will be along dz

+CD=integrate('5','Q',0,%pi/2);

+disp(CD);   //dl will be along d(phi)
\ No newline at end of file
diff --git a/42/CH3/EX3.2/sadiku_3_2.sce b/42/CH3/EX3.2/sadiku_3_2.sce
new file mode 100755
index 000000000..018d0aec3
--- /dev/null
+++ b/42/CH3/EX3.2/sadiku_3_2.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+C1=integrate('x^2','x',1,0);//for y=0=z

+C2=0;// as (az.ay)=0

+C3=integrate('x^2 -1','x',0,1);

+C4=integrate('-y-y^2','y',1,0);

+C=C1+C2+C3+C4;

+disp(C);
\ No newline at end of file
diff --git a/42/CH3/EX3.9/sadiku_3_9.sce b/42/CH3/EX3.9/sadiku_3_9.sce
new file mode 100755
index 000000000..371f825f7
--- /dev/null
+++ b/42/CH3/EX3.9/sadiku_3_9.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+ab=integrate('2*sin(P)','P',%pi/3,%pi/6);

+bc=(3^.5 /2)*integrate('p','p',2,5);

+Cd=integrate('5*sin(P)','P',%pi/6,%pi/3);

+da=.5*integrate('p','p',5,2);

+C1=ab+bc+Cd+da;

+disp(C1, 'C1=');

+C2=integrate('sin(Q)','Q',%pi/6,%pi/3)*integrate('(1+p)','p',2,5);

+disp(C2,'C2=');

+disp('Since C1=C2 hence stroke theorem is proved');

+

diff --git a/42/CH4/EX4.1/sadiku_4_1.sce b/42/CH4/EX4.1/sadiku_4_1.sce
new file mode 100755
index 000000000..f93360bff
--- /dev/null
+++ b/42/CH4/EX4.1/sadiku_4_1.sce
@@ -0,0 +1,16 @@
+clear;
+clc;
+format('v',6);
+Q1=1;
+Q2=-2;
+Q=10*10^-9;
+P1=[0 3 1]-[3 2 -1];
+P2=[0 3 1]-[-1 -1 4];
+
+e=10^-9/(36*%pi);
+det1=(P1(1,1)^2+P1(1,2)^2+P1(1,3)^2)^.5;
+det2=(P2(1,1)^2+P2(1,2)^2+P2(1,3)^2)^.5;
+F=[[(Q*Q1)*(P1)]/(4*%pi*e*(det1)^3)]+[[(Q*Q2)*(P2)]/(4*%pi*e*(det2)^3)];
+E=[(10^-6)*(F/Q)];
+disp(F,'F(in mN)=');
+disp(E,'At that point E(in kV)=');
\ No newline at end of file
diff --git a/42/CH4/EX4.10/sadiku_4_10.sce b/42/CH4/EX4.10/sadiku_4_10.sce
new file mode 100755
index 000000000..2c7e2008a
--- /dev/null
+++ b/42/CH4/EX4.10/sadiku_4_10.sce
@@ -0,0 +1,13 @@
+clear;
+clc;
+format('v',6);
+Q1=-4;
+Q2=5;
+R1=[1 0 1]-[2 -1 3];
+R2=[1 0 1]-[0 4 -2];
+e=10^-9/(36*%pi);
+mod_R1=(R1(1,1)^2+R1(1,2)^2+R1(1,3)^2)^.5;
+mod_R2=(R2(1,1)^2+R2(1,2)^2+R2(1,3)^2)^.5;
+C0=0;
+V=10^-6*(([Q1/mod_R1]+[Q2/mod_R2])/(4*%pi*e))+C0;
+disp(V*10^-3,'V(1,0,1)(in kV)=');
\ No newline at end of file
diff --git a/42/CH4/EX4.12/sadiku_4_12.sce b/42/CH4/EX4.12/sadiku_4_12.sce
new file mode 100755
index 000000000..a3ca49a7b
--- /dev/null
+++ b/42/CH4/EX4.12/sadiku_4_12.sce
@@ -0,0 +1,10 @@
+clear;

+clc;

+q=10*10^-6;

+function[V]=pot(r,P,Q);

+V=10*sin(P)*cos(Q)/r^2;

+endfunction

+Va=pot(1,%pi/6,2*%pi/3);

+Vb=pot(4,%pi/2,%pi/3);

+W=q*(Vb-Va);

+disp(W*10^6,'Work done in uJoule');
\ No newline at end of file
diff --git a/42/CH4/EX4.13/sadiku_4_13.sce b/42/CH4/EX4.13/sadiku_4_13.sce
new file mode 100755
index 000000000..221a35756
--- /dev/null
+++ b/42/CH4/EX4.13/sadiku_4_13.sce
@@ -0,0 +1,6 @@
+clear;

+clc;

+p1=-5*10^-9, p2=9*10^-9;

+r1=2,r2=-3,e=10^-9/(36*%pi);

+V=(1/(4*%pi*e))*((p1*abs(r1)/r1^3)+(p2*abs(r2)/r2^3));

+disp(V);
\ No newline at end of file
diff --git a/42/CH4/EX4.14/sadiku_4_14.sce b/42/CH4/EX4.14/sadiku_4_14.sce
new file mode 100755
index 000000000..77309a283
--- /dev/null
+++ b/42/CH4/EX4.14/sadiku_4_14.sce
@@ -0,0 +1,7 @@
+clear;

+clc;

+format('v',6);

+Q1=-1*10^-9 ,Q2=4*10^-9,Q3=3*10^-9,e=10^-9/(36*%pi);

+V1=(1/(4*%pi*e) * (Q2+Q3)),V2=(1/(4*%pi*e)*(Q1+Q3/(2^.5)) ),V3=(1/(4*%pi*e) * (Q1+Q2/(2^.5)));

+W=.5*((V1*Q1)+(V2*Q2)+(V3*Q3));

+disp(W*10^9,'Energy in nJ');
\ No newline at end of file
diff --git a/42/CH4/EX4.6/sadiku_4_6.sce b/42/CH4/EX4.6/sadiku_4_6.sce
new file mode 100755
index 000000000..acb9d7144
--- /dev/null
+++ b/42/CH4/EX4.6/sadiku_4_6.sce
@@ -0,0 +1,15 @@
+clear;

+clc;

+format('v',6);

+p1=10*10^-9;

+p2=15*10^-9;

+pl=10*%pi*10^-9;

+e=(10^-9)/(36*%pi);

+E1=(p1/(2*e))*[-1 0 0];

+E2=(p2/(2*e))*[0 1 0];

+R=[1 0 -3];

+p=(R(1,1)^2+R(1,2)^2+R(1,3)^2);

+a=R/p;

+E3=(pl/(2*%pi*e))*a;

+E=E1+E2+E3;

+disp(E,'E(in V) at (1,1,-1)=');
\ No newline at end of file
diff --git a/42/CH4/EX4.7/sadiku_4_7.sce b/42/CH4/EX4.7/sadiku_4_7.sce
new file mode 100755
index 000000000..2f05d77d3
--- /dev/null
+++ b/42/CH4/EX4.7/sadiku_4_7.sce
@@ -0,0 +1,16 @@
+clear;
+clc;
+format('v',12);
+e=10^-9;
+Q=-5*%pi*10^-3;
+pl=3*%pi*10^-3;
+r=[4 0 3];
+p=(r(1,1)^2+r(1,2)^2+r(1,3)^2)^.5;
+r1=[4,0,0];
+R=r-r1;
+mod_R=(R(1,1)^2+R(1,2)^2+R(1,3)^2)^.5;
+Dq=(Q*R)/(4*%pi*mod_R^3);
+ap=r/p;
+Dl=(pl/(2*%pi*p))*ap;
+D=Dq+Dl;
+disp(D*10^6,'Flux density D(in microC) due to a point charge and a infinite line charge');
\ No newline at end of file
diff --git a/42/CH4/EX4.8/sadiku_4_8.sce b/42/CH4/EX4.8/sadiku_4_8.sce
new file mode 100755
index 000000000..9fc783845
--- /dev/null
+++ b/42/CH4/EX4.8/sadiku_4_8.sce
@@ -0,0 +1,5 @@
+clear;

+clc;

+r1=0,r2=1,z1=-2,z2=2,q1=0,q2=2*%pi;

+Q=integrate('p^2','p',r1,r2)*integrate('(cos(Q)^2)','Q',q1,q2)*integrate('1','z',z1,z2);

+disp(Q,'Total charge is =');
\ No newline at end of file
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)=');
+
diff --git a/42/CH6/EX6.12/sadiku_6_12.sce b/42/CH6/EX6.12/sadiku_6_12.sce
new file mode 100755
index 000000000..f7d4d44b8
--- /dev/null
+++ b/42/CH6/EX6.12/sadiku_6_12.sce
@@ -0,0 +1,11 @@
+clear;

+clc;

+Eo=10^-9 /(36*%pi),Er1=4,Er2=6,d=5*10^-3,S=30*10^-4;

+C1=Eo*Er1*S*2/d;

+C2=Eo*Er2*S*2/d;

+C=C1*C2/(C1+C2);//Since they are in series

+disp(C*10^12,'Capacitance of capacitor in figure a in pF =');

+C1=Eo*Er1*S/(2*d);

+C2=Eo*Er2*S/(2*d);

+C=C1+C2;

+disp(C*10^12,'Capacitance of capacitor in figure b in pF = ')
\ No newline at end of file
diff --git a/42/CH7/EX7.1/sadiku_7_1.sce b/42/CH7/EX7.1/sadiku_7_1.sce
new file mode 100755
index 000000000..4d78bd6ba
--- /dev/null
+++ b/42/CH7/EX7.1/sadiku_7_1.sce
@@ -0,0 +1,5 @@
+clear;

+clc;

+a1=acos(0),a2=acos(2/29^.5),p=5,I=10;

+H=I/(4*%pi*p)*(cos(a1)-cos(a2));

+disp(H*1000,'H at (0,0,5) in mA ');
\ No newline at end of file
diff --git a/42/CH7/EX7.2/sadiku_7_2.sce b/42/CH7/EX7.2/sadiku_7_2.sce
new file mode 100755
index 000000000..e5c48bdf7
--- /dev/null
+++ b/42/CH7/EX7.2/sadiku_7_2.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+a1=acos(0),a2=acos(1),p=5,I=3;

+Hz=I/(4*%pi*p)*(cos(a2)-cos(a1))*[.8 .6 0];

+a2=acos(1),a1=acos(.6),p=4,I=3;

+Hx=I/(4*%pi*p)*(cos(a2)-cos(a1))*[0 0 1];

+H=Hx+Hz;

+disp(H*1000,'H at (0,0,5) in mA ');
\ No newline at end of file
diff --git a/42/CH7/EX7.5/sadiku_7_5.sce b/42/CH7/EX7.5/sadiku_7_5.sce
new file mode 100755
index 000000000..6c15dbcb5
--- /dev/null
+++ b/42/CH7/EX7.5/sadiku_7_5.sce
@@ -0,0 +1,11 @@
+clear;

+clc;

+i0=-10,i4=10;

+H0=.5*i0*-1;// in the positive Y direction

+H4=.5*i4*-1*-1;//in the positive Y direction

+H=H0+H4;

+disp(H, 'H at (1,1,1) =')

+H0=.5*i0*-1;//in the positive Y direction

+H4=.5*i4*-1;//in the negative Y direction

+H=H0+H4;

+disp(H,'H at (0,-3,10 =)');
\ No newline at end of file
diff --git a/42/CH7/EX7.7/sadiku_7_7.sce b/42/CH7/EX7.7/sadiku_7_7.sce
new file mode 100755
index 000000000..b5748756c
--- /dev/null
+++ b/42/CH7/EX7.7/sadiku_7_7.sce
@@ -0,0 +1,6 @@
+clear;

+clc;

+disp('Vector potential A=-p^2/4');

+Q=%pi/2,p1=1,p2=2,z1=0,z2=5

+Y=.5*integrate('p','p',p1,p2)*integrate('1','z',z1,z2);

+disp(Y,'Total magnetic flux= ')
\ No newline at end of file
diff --git a/42/CH8/EX8.1/sadiku_8_1.sce b/42/CH8/EX8.1/sadiku_8_1.sce
new file mode 100755
index 000000000..2d37b2dfe
--- /dev/null
+++ b/42/CH8/EX8.1/sadiku_8_1.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+m=2,q=3,v=[4 0 3],E=[12 10 0],t=1;

+disp(q*E/m,'Acceleration ofthe particle= ');

+u=[22 15 3];

+modofu=sqrt(22*22+15*15+3*3);

+KE=.5*m*(modofu)^2;

+disp(KE, 'Kinetic energy= ')
\ No newline at end of file
diff --git a/42/CH8/EX8.14/sadiku_8_14.sce b/42/CH8/EX8.14/sadiku_8_14.sce
new file mode 100755
index 000000000..29e271009
--- /dev/null
+++ b/42/CH8/EX8.14/sadiku_8_14.sce
@@ -0,0 +1,6 @@
+clear;

+clc;

+p=10*10^-2 ,a=1*10^-2 ,Ur=1000, Uo=4*%pi*10^-7,n=200,phi=.5*10^-3;

+U=Uo*Ur;

+I=phi*2*%pi*p/(Uo*Ur*n*%pi*a*a);

+disp(I);
\ No newline at end of file
diff --git a/42/CH8/EX8.15/sadiku_8_15.sce b/42/CH8/EX8.15/sadiku_8_15.sce
new file mode 100755
index 000000000..82760908c
--- /dev/null
+++ b/42/CH8/EX8.15/sadiku_8_15.sce
@@ -0,0 +1,10 @@
+clear;

+clc;

+Uo=4*%pi*10^-7,Ur=50,l1=30*10^-2,s=10*10^-4,l3=9*10^-2,la=1*10^-2,B=1.5,N=400;

+R1=l1/(Uo*Ur*s);R2=R1;

+R3=l3/(Uo*Ur*s);

+Ra=la/(Uo*s);

+R=R1*R2/(R1+R2);

+Req=R3+Ra+R;

+I=B*s*Req/N;

+disp(I,'Required current= ');

diff --git a/42/CH8/EX8.16/sadiku_8_16.sce b/42/CH8/EX8.16/sadiku_8_16.sce
new file mode 100755
index 000000000..f371b3003
--- /dev/null
+++ b/42/CH8/EX8.16/sadiku_8_16.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+m=400,g=9.8,Ur=3000, Uo=4*%pi*10^-7,S=40*10^-4,la=1*10^-4,li=50*10^-2,I=1;

+B=sqrt(m*g*Uo/S);

+Ra=2*la/(Uo*S);

+Ri=li/(Uo*Ur*S);

+N=(Ra+Ri)/(Ra*Uo)*B*la;

+disp(N,'No of turns= ');
\ No newline at end of file
diff --git a/42/CH8/EX8.8/sadiku_8_8.sce b/42/CH8/EX8.8/sadiku_8_8.sce
new file mode 100755
index 000000000..094dea3fa
--- /dev/null
+++ b/42/CH8/EX8.8/sadiku_8_8.sce
@@ -0,0 +1,9 @@
+clear;

+clc;

+format('v',6);

+H1=[-2 6 4],Uo=4*%pi*10^-7,Ur=5;

+U1=Uo*Ur;

+M1=(Ur-1)*H1;

+disp(M1,'M = ');

+B1=U1*H1;

+disp(B1*10^6,'B in uW/m^2');
\ No newline at end of file
diff --git a/42/CH9/EX9.5/sadiku_9_5.sce b/42/CH9/EX9.5/sadiku_9_5.sce
new file mode 100755
index 000000000..e446f74e4
--- /dev/null
+++ b/42/CH9/EX9.5/sadiku_9_5.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+z3=%i,z4=3+4*%i,z5=-1+6*%i,z6=3+4*%i;

+z1=(z3*z4/(z5*z6));

+disp(z1,'z1 =');

+z7=1+%i, z8=4-8*%i;

+z2=(z7/z8)^.5;

+disp(z2,'z2 = ')
\ No newline at end of file