1 files changed, 31 insertions, 0 deletions
diff --git a/2891/CH2/EX2.22/Ex2_22.sce b/2891/CH2/EX2.22/Ex2_22.sce
new file mode 100755
index 000000000..17fc3b4c3
--- /dev/null
+++ b/2891/CH2/EX2.22/Ex2_22.sce
@@ -0,0 +1,31 @@
+// Exa 2.22
+clc;
+clear;
+close;
+// given :
+f1=60 // frequency in Hz
+omega1=2*%pi*f1 // angular frequency in Hz
+f2=100 // frequency in MHz
+f2=100*10^6 // frequency in Hz
+omega2=2*%pi*f2 // angular frequency in Hz
+sigma=5.8*10^7 // conductivity in mho/m
+epsilon_0=8.854*10^-12 // permittivity in free space in F/m
+mu_0=4*%pi*10^-7 //  permeability in free space in H/m
+epsilon_r=1 // relative permittivity
+mu_r=1 // relative permeability
+epsilon=epsilon_r*epsilon_0 // permittivity
+mu=mu_0*mu_r // permeability
+
+disp("At f=60Hz")
+k1=(sigma)/(omega1*epsilon) // ratio
+disp(k1," ratio k is equal to")
+disp("since k>>1 therefore it is very good conductor at f=60Hz:")
+delta1=(sqrt(2/(omega1*mu*sigma))) // depth of penetration in m
+disp(delta1,"depth of penetration delta1 in m:")
+
+disp("At f=100Hz")
+k2=sigma/(omega2*epsilon) // ratio
+disp(k2,"ratio k is equal to")
+disp("since k2>>1 therefore it is very good conductor at f=100Hz:")
+delta2=(sqrt(2/(omega2*mu*sigma))) // depth of penetration in m
+disp(delta2,"depth of penetration delta2 in m:")