initial commit / add all books

1 files changed, 32 insertions, 0 deletions

diff --git a/803/CH6/EX6.1/ex6_1.sce b/803/CH6/EX6.1/ex6_1.sce
new file mode 100644
index 000000000..3e250b662
--- /dev/null
+++ b/803/CH6/EX6.1/ex6_1.sce
@@ -0,0 +1,32 @@
+clc

+f=1.2*(10^6);...............................//Assigning values to parameters

+Yn=56;

+deltalat=37.45;

+x=23.44*sind(0.9856*(Yn-80.7));

+deltah=acosd(tand(x)/tand(deltalat));

+y=acosd((sind(deltalat)*sind(x))+(cosd(deltalat)*cosd(x)*cosd(deltah)));

+disp("degree",x,"Solar declination");

+disp("degree",deltah,"hour angle of the sun");

+disp("degree",y,"Solar zenith angle");

+m=4.8*(10^(-26));................//mean molecular mass of air

+k=1.38*(10^(-23));...............//Boltzmann constant

+g=9.8;...........................//gravitational constant

+e=1.6*(10^(-19));....//electron charge

+me=9.11*(10^(-31));

+epsilon=8.85*(10^(-12));

+hmax=[100 200 300 120 250];

+T=[341 1360 1710 341 1540];

+f=1.2*(10^6);

+Nm=[1.5*10^11 3*10^11 12.5*10^11 0.8*1010 4*1011];

+h=[122 256 335 132 276];

+for i=1:5

+    H(i)=(g * m * [h(i)-hmax(i)])/ (k * T(i));

+    No(i)=(Nm(i)*(secd(x)^0.5));

+    Ne(i)=((No(i))*exp(0.5*(1-(H(i)-secd(x)*exp(-H(i))))));

+    fc(i)=(1/(2*3.14))*(sqrt(((e^2)*(Ne(i)))/(epsilon*me)));

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

+    disp(Ne(i),"Electron density of each layer:");

+    disp(fc(i),"Critical frequency of each layer:");

+    disp(n(i),"Refractive index of each layer:");       

+end

+