5 files changed, 101 insertions, 0 deletions
diff --git a/3773/CH25/EX25.1/Ex25_1.sce b/3773/CH25/EX25.1/Ex25_1.sce
new file mode 100644
index 000000000..04cf3d33f
--- /dev/null
+++ b/3773/CH25/EX25.1/Ex25_1.sce
@@ -0,0 +1,19 @@
+//Chapter 25: Sky Wave Propagation
+//Example 25-5.1
+clc;
+
+//Variable Initialization
+muf = 10e6    //Maximum usable frequency (Hz)
+h = 300    //Height of reflection (km)
+n = 0.9      //Maximum value of refractive index (unitless)
+
+//Calculations
+Nmax = (1 - n**2)*(muf**2)/81    //Max. Number of electrons per cubic cm
+fc = 9*sqrt(Nmax)    //Critical frequency (Hz)
+dskip = 2*h*sqrt((muf/fc)**2 - 1)   //Skip distance (km)
+
+
+//Result
+mprintf("The skip distance is %.1f km",dskip)
+
+//An error has been made in the calculation of sqrt((muf/fc)**2 - 1)
diff --git a/3773/CH25/EX25.2/Ex25_2.sce b/3773/CH25/EX25.2/Ex25_2.sce
new file mode 100644
index 000000000..789d4a66d
--- /dev/null
+++ b/3773/CH25/EX25.2/Ex25_2.sce
@@ -0,0 +1,18 @@
+//Chapter 25: Sky Wave Propagation
+//Example 25-5.2
+clc;
+
+//Variable Initialization
+fE = 3e6        //Critical frequency for E layer (Hz)
+fF1 = 5e6       //Critical frequency for F1 layer (Hz)
+fF2 = 9e6       //Critical frequency for F2 layer (Hz)
+
+//Calculations
+N_E = (fE**2)/81    //Concentration of electrons in E layer (per cubic cm)
+N_F1 = (fF1**2)/81 //Concentration of electrons in F1 layer (per cubic cm)
+N_F2 = (fF2**2)/81 //Concentration of electrons in F2 layer (per cubic cm)
+
+//Result
+mprintf( "The concentration of electrons in E layer is %e per cubic cm",N_E)
+mprintf( "\nThe concentration of electrons in F1 layer is %e per cubic cm", N_F1)
+mprintf( "\nThe concentration of electrons in F2 layer is %e per cubic cm", N_F2)
diff --git a/3773/CH25/EX25.3/Ex25_3.sce b/3773/CH25/EX25.3/Ex25_3.sce
new file mode 100644
index 000000000..9113cfc27
--- /dev/null
+++ b/3773/CH25/EX25.3/Ex25_3.sce
@@ -0,0 +1,21 @@
+
+//Chapter 25: Sky Wave Propagation
+//Example 25-5.3
+clc;
+
+//Variable Initialization
+N_E = 0.8*0.111e12  //Concentration of electrons in E layer (per cubic cm)
+N_F1 = 0.8*0.3086e12 //Concentration of electrons in E layer (per cubic cm)
+N_F2 = 0.8*1e12     //Concentration of electrons in E layer (per cubic cm)
+
+//Calculations
+fE = 9*sqrt(N_E)    //Critical frequency in E layer (Hz)
+fF1 = 9*sqrt(N_F1)  //Critical frequency in F1 layer (Hz)
+fF2 = 9*sqrt(N_F2)  //Critical frequency in F2 layer (Hz)
+
+//Result
+disp(fE,"The Critical frequency in E layer in Hz")
+disp(fF1,"The Critical frequency in F1 layer in Hz")
+disp(fF2,"The Critical frequency in F2 layer in Hz")
+
+//The difference appearing for fE,fF1 is a result of approximation
diff --git a/3773/CH25/EX25.4/Ex25_4.sce b/3773/CH25/EX25.4/Ex25_4.sce
new file mode 100644
index 000000000..d7a716549
--- /dev/null
+++ b/3773/CH25/EX25.4/Ex25_4.sce
@@ -0,0 +1,23 @@
+//Chapter 25: Sky Wave Propagation
+//Example 25-6.1
+clc;
+
+//Variable Initialization
+hD = 70     //Height of D layer (km)
+hE = 130     //Height of E layer (km)
+hF1 = 230     //Height of F1 layer (km)
+hF2 = 350     //Height of F2 layer (km)
+theta = 10*%pi/180      //Angle of incidence (radians)
+
+//Calculations
+temp = sqrt((cos(theta))**-2 - 1)
+d1 = 2*hD*temp  //Maximum single hop distance for D layer (km)
+d2 = 2*hE*temp  //Maximum single hop distance for E layer (km)
+d3 = 2*hF1*temp //Maximum single hop distance for F1 layer (km)
+d4 = 2*hF2*temp //Maximum single hop distance for F2 layer (km)
+
+//Result
+mprintf( "The Maximum single hop distance for D layer is %.1f km", d1)
+mprintf( "\nThe Maximum single hop distance for E layer is %.2f km", d2)
+mprintf( "\nThe Maximum single hop distance for F1 layer is %.2f km", d3)
+mprintf( "\nThe Maximum single hop distance for F2 layer is %.1f km", d4)
diff --git a/3773/CH25/EX25.5/Ex25_5.sce b/3773/CH25/EX25.5/Ex25_5.sce
new file mode 100644
index 000000000..295a1bacf
--- /dev/null
+++ b/3773/CH25/EX25.5/Ex25_5.sce
@@ -0,0 +1,20 @@
+//Chapter 25: Sky Wave Propagation
+//Example 25-9.1
+clc;
+clear;
+
+//Variable Initialization
+d = 200     //Height of layer (km)
+bet = 20    //Takeoff angle (degrees)
+R = 6370    //Earth's radius (km)
+
+//Calculations
+phi_0 = 90 - bet   //Take off angle for flat earth (degrees)
+h = (d/2)/(sqrt((cos(phi_0*%pi/180)**-2) - 1))    //Skip distance for case (a) (km)
+
+phi_02 = 90 - bet - 57.2*d/(2*R)                    //Take off angle for spherical earth (degrees)
+h2 = (d/2)/(sqrt((cos(phi_02*%pi/180)**-2) - 1))                   //Skip distance for case (b) (km)
+
+//Result
+mprintf("The skip distance for case (a) is %.3f km", h)
+mprintf("\nThe skip distance for case (b) is %.2f km", h2)