CLS

5 PRINT
PRINT "SADEYA/INFORTECNICA"
PRINT
PRINT "INTENSIDAD SEĽAL RADIO"

INPUT "RECEPTOR 1 LIMITE 2 ORIENTACION 3"; TIP
IF TIP = 3 THEN GOTO 1000
IF TIP = 2 THEN GOTO 100
TIP = 1
10 PRINT
INPUT "DISTANCIA EMISORA EN Km."; DK
DK = DK * 1000
ALTM = (DK * DK) / 12742000




20 INPUT "POTENCIA DE EMISION WAT."; KW
   INPUT "LIMITE POTENCIA DECIBEL (EJEM.50)"; DELIM
   INPUT "EXPONENTE EXTINCION (3 - 3.1)"; EEX

KW = KW * 1000
IT = KW / (DK ^ EEX)

DEC = (LOG(IT) / LOG(10)) * 10 + 125
AA = 10 ^ ((DELIM - 125) / 10)
DK2 = (KW / AA) ^ (1 / EEX)


PRINT "DECIBELIOS RECEP. LIMITE="; DEC

IF DEC > 75 THEN PRINT "SEĽAL EXCELENTE"
IF DEC < 75 AND DEC > 65 THEN PRINT "SEĽAL BUENA"
IF DEC < 65 AND DEC > 55 THEN PRINT "SEĽAL ACEPTABLE"
IF DEC < 55 AND DEC > 45 THEN PRINT "SEĽAL DEBIL"
IF DEC < 45 THEN PRINT "SEĽAL DIFICIL"

PRINT " ALTURA MINIMA METROS="; ALTM
REM POL = 10 ^ (DK ^ EEX) * 10 + 50
REM PRINT "LIMITE POTENCIA"; POL

PRINT "LIMITE DISTANCIA/POT. METROS="; DK2; "PARA DECIBEL="; DELIM


PRINT

GOTO 5
END

100 INPUT " ALTURA EMISOR EN METROS "; ALTM
HORI = SQR(ALTM * 12742000)
LIMIT = 4 / 3 * HORI
PRINT "HORIZONTE METROS ="; HORI
PRINT "MAXIMO INCLUIDA REF. METROS="; LIMIT
DK = LIMIT


GOTO 20

1000  INPUT "LONGITUD Y LATITUD GRADOS DE LA ANTENA"; LON1, lat1
INPUT "LONGITUD Y LATITUD GRADOS DEL RECEPTOR"; LON2, lat2
     rad = 180 / 3.1416
     
     LON1 = LON1 / rad
     LON2 = LON2 / rad
     
     lat1 = lat1 / rad
     lat2 = lat2 / rad
     
     GOTO 500

     ACO = (SIN(LON1) * SIN(LON2)) + (COS(LON1) * COS(LON2) * (COS(lat2 - lat1)))

     ATA = SQR(1 - ACO * ACO) / ACO
     

     D = 60 * ATN(ATA) * (COS((lat2 + lat1) / 2))

     DK = D * 12750 * 3.1416 / 360

     PRINT "DISTANCIA KM.="; DK
     
     ACOO = (SIN(LON2) - SIN(LON1) * COS(D / 60)) / (SIN(D / 60) * COS(LON1))
    
     ATAA = SQR(1 - ACOO * ACOO) / ACOO
     
     H = ATN(ATAA)
    
     PRINT "DISTANCIA GRAD.="; D
     PRINT "AZIMUT="; H * rad
500 PRINT

     INL = (LON2 - LON1) * COS((lat2 + lat1) / 2) * 40000 / 360 * rad
     INB = (lat2 - lat1) * 40075 / 360 * rad
     DDK = SQR(INL * INL + INB * INB)
     hh = ATN(INB / INL)
     snl = ABS(INL) / INL
     snb = ABS(INB) / INB

     PRINT "DIF. LONGITUD KM.="; INL
     PRINT "DIF. LATITUD KM.="; INB
     
    
     PRINT "DISTANCIA Km.="; DDK
     if DDK>100 THEN PRINT "DISTANCIA FUERA DE LIMITE"
     hh = hh * rad

     IF snl > 0 AND snb > 0 THEN hz = 90 - hh
     IF snl > 0 AND snb < 0 THEN hz = 90 + ABS(hh)
     IF snl < 0 AND snb > 0 THEN hz = 270 + ABS(hh)
     IF snl < 0 AND snb < 0 THEN hz = 270 - ABS(hh)
     
     
    REM  PRINT "snl,snb"; snl, snb
    
     PRINT "ANGULO GRAD.="; hh
     hz = 360-hz
     PRINT "AZIMUTH DESDE LA ANTENA AL RECEPTOR.="; hz
     hz2 = 360 - (hz + 180)
     IF hz2 > 360 THEN hz2 = hz2 - 360
     PRINT "AZIMUTH DESDE EL RECEPTOR A LA ANTENA.="; hz2
     
     PRINT
     GOTO 5