from __future__ import division
import math
fcr=11*10**6;
D=1000;
h=400;
fmuf=fcr*math.sqrt(1+(D/(2*h))**2);
print("The maximum stable frequency is %g Hz"%fmuf);
from __future__ import division
import math
Nmax=10**11;
phi=(math.pi)/9;
fcr=math.sqrt(81*Nmax);
print("The critical frequency is %g Hz"%fcr);
fmuf=fcr*(1/math.cos(phi));
print("The maximum usable frequency is %g Hz"%fmuf);
from __future__ import division
import math
D=2000;
h=200;
fmuf=30.6*10**6;
fcr=fmuf/math.sqrt(1+(D/(2*h))**2);
print("The critical frequency is %g Hz"%fcr);
from __future__ import division
import math
n=0.9;
fmuf=10*10**6;
f=10*10**6;
h=400*10**3;
Nmax=(1-n**2)*f**2/81;
print("The Nmax value is %g /m^3"%Nmax);
fcr=math.sqrt(81*Nmax);
print("The critical frequency is %g Hz"%fcr);
Dskip=2*h*math.sqrt((fmuf/fcr)**2-1);
print("The skip distance is %g m"%Dskip);
from __future__ import division
import math
ht=150;
hr=2;
Is=9;
d=40*10**3;
f=1.2*10**6;
c=3*10**8;
lamda=c/f;
print("The wavelength is %d m"%lamda);
E=120*(math.pi)*ht*hr*Is/(lamda*d);
print("The electric field is %g V/m"%E);
from __future__ import division
import math
dmax=45*10**3;
ht=(dmax/8.24)**2; #dmax=4.12[sqrt(ht)+sqrt(hr)];ht=hr;
print("The height of transmission is %g m"%ht);
from __future__ import division
import math
fcre=2.5*10**6;
fcrf=8.5*10**6;
Nmaxe=(fcre)**2/81;
Nmaxf=(fcrf)**2/81;
print("The Nmax for e layer is %g /m^3"%Nmaxe);
print("The Nmax for f layer is %g /m^3"%Nmaxf);
from __future__ import division
import math
Nmaxf1=2.5;
Nmaxf2=3.5;
Nmaxf3=1.5;#10^6*10^-6=1;
fcr1=math.sqrt(81*Nmaxf1);
fcr2=math.sqrt(81*Nmaxf2);
fcr3=math.sqrt(81*Nmaxf3);
print("The critical frequencies are %gHz %gHz %gHz"%(fcr1,fcr2,fcr3));
from __future__ import division
import math
fcr1=4.5*10**6;
fcr2=1.5*10**6;
Nmax1=(fcr1/9)**2;
Nmax2=(fcr2/9)**2;
print("The Nmax values are %g m^3 %g m^3"%(Nmax1,Nmax2));
Nmax=Nmax1-Nmax2;
print("The change in electron density is %g m^3"%Nmax);
from __future__ import division
import math
#Note:10^6 is the power and not 10^-6 as mentioned in book
n=0.5;
N=400*10**6;
f=math.sqrt((81*N)/(1-n**2));
print("The frequency is %e Hz"%f);
from __future__ import division
import math
D=1500;
h=250;
fmuf=37.95*10**6;
fcr=fmuf/math.sqrt(1+(D/(2*h))**2);
print("The critical frequency is %e Hz"%fcr);
from __future__ import division
import math
D=2500;
h=200;
fcr=5*10**6;
fmuf=fcr*math.sqrt(1+(D/(2*h))**2);
print("The maximum usable frequency is %g Hz"%fmuf);
from __future__ import division
import math
T=5*10**-3;
c=3*10**8;
h=c*(T/2);
print("The virtual height is given by %g m"%h);
from __future__ import division
import math
ht=40;
hr=25;
f=90*10**6;
p=35;
LOS=4.12*(math.sqrt(ht)+math.sqrt(hr));
print("The line of sight distance is %g m"%LOS);
from __future__ import division
import math
Nmax=1.26*10**12;
fcr=math.sqrt(81*Nmax);
print("The critical frequency is %g Hz"%fcr)
from __future__ import division
import math
Nmax=1.24*10**12;
fcr=math.sqrt(81*Nmax);
print("The critical frequency is %g Hz"%fcr);
from __future__ import division
import math
fcr=6*10**6;
D=200*10**3;
h=200*10**3;
fmuf=fcr*math.sqrt(1+(D/(2*h))**2);
print("The maximum usable frequency is %g Hz"%fmuf);
from __future__ import division
import math
ht=100;
hr=50;
d=1.4142*(math.sqrt(ht)+math.sqrt(hr));
print("The maximum range is %g miles"%d);