1 CLS
2 PRINT "SADEYA/INFORTECNICA"
PRINT
PRINT "SUPLEMENTO ESTELAR"
PRINT "CALCULO DE LA DISTANCIA, VELOCIDAD TANGENCIAL Y REAL DE UNA ESTRELLA"
PRINT
10 INPUT "Velocidad radial Km/s, si se conoce sino 1"; vr
INPUT "Constante velocidad real estrella. Estandar Vradial*2 (>10 <100)Km/s"; kkk
kkk = kkk / 4.7

INPUT "Movimiento propio en A.R. milesimas de arcosegundos ,,"; xx
xx = xx / 1000
INPUT "Movimiento propio Dec. milesimas de arcosegundos ,,"; yy
yy = yy / 1000



INPUT "Paralaje en milesimas de arcosegundos obervado o experimental"; parae
parae = parae / 1000



dp = 1 / parae


11 PRINT

vts = SQR(xx * xx + yy * yy)

lonal = dp * 3.26156 * 2 * 3.1415926#
lonk = lonal * 365.25 * 24 * 3600 * 300000
segkm = lonk / 1296000
vtana = segkm * vts
vtans = vtana / (365.25 * 24 * 3600)


anpos = ATN(yy / xx) * 180 / 3.1415926#


vtk = vtans

vre = SQR(vtk * vtk + vr * vr)
IF vr < 0 THEN vre = vre * -1

dal = dp * 3.2616
ana = ATN(vtk / vr) * 180 / 3.1416
anb = 90 - ana

cosa = (vr / vre)
rec = cosa * dal

disp = SQR((dal ^ 2) - (rec ^ 2))
anb1 = ATN(rec / disp) * (180 / 3.1416)
anb2 = anb

tperi = rec * (300000 / vre) * -1

vres = (vre / vtk) * vts



parals = vres
dispar = ABS((1 / parals) * kkk)





PRINT "Distancia en parsecs media por mov."; dispar
INPUT "            Si nueva distancia 1 finalizar si paralaje valido 2"; ndis
IF ndis = 1 THEN dp = dispar
IF ndis = 1 THEN 11
IF xx < 0 THEN anpos = anpos + 180
PRINT "Distancia en parsecs, paral. arc.seg"; dispar; 1 / dispar

PRINT "Velocidad tangencial ,,="; vts
PRINT "Velocidad tangencial Km/s.="; vtk
PRINT "Velocidad real Km/s="; vre
PRINT "Velocidad real en ,,"; vres
PRINT "Distancia ans luz"; dal
PRINT "Distancia perigeo"; disp
PRINT "Tiempo perigeo ans"; tperi
PRINT "Angulo del perigeo"; anb1
PRINT "Angulo de posicion del perigeo"; anpos
INPUT "Magnitud aparente visual V "; maga
difmag = -5 + 5 * LOG(dp) / LOG(10)
mab = maga - difmag
PRINT "Difmag.="; difmag

PRINT "Magnitud absoluta transp="; mab
INPUT "Indice mag B-V"; icol
icol = icol * 1.00001
IF icol = 0 THEN icol = .00001


IF icol < (0) THEN PRINT "(Si negativo O), Tipo espectral B"; INT((-.3 - icol) / .3 * 10) * -1
IF icol > (.0001) AND icol < (.3) THEN PRINT "Tipo espectral A"; INT((icol - 0) / .3 * 10)
IF icol > (.300001) AND icol < (.6) THEN PRINT "Tipo espectral F"; INT((icol - .3) / .3 * 10)
IF icol > (.60001) AND icol < (1) THEN PRINT "Tipo espectral G"; INT((icol - .6) / .4 * 10)
IF icol > (1.00001) AND icol < (1.4) THEN PRINT "Tipo espectral K"; INT((icol - 1) / .5 * 10)
IF icol > (1.40001) THEN PRINT "Tipo espectral M"; INT((icol - 1.5) / .5 * 10)



MAPR = icol * 8
MAGG = (icol * 3 - 3) * -1
IF mab < (-3) THEN PRINT "Probable supergigante de magabs"; mab
IF mab > (10) AND icol < (1.5) THEN PRINT "Probable estrella enana blanca de magabs"; mab

PRINT "Si es de la serie principal V magabs med="; MAPR; " magabs cal."; mab
PRINT "Diferencia de magnitudes"; ABS(MAPR - mab)

PRINT "Si es gigante (III) magabs med="; MAGG; " magabs cal."; mab
PRINT "Diferencia de magnitudes"; ABS(MAGG - mab)






PRINT " O, determinar si la posicion es aceptable en el diagrama Her. Rusell"
PRINT
INPUT "Correccion, por nueva magnitud absoluta"; mab
difmag = maga - mab
logadp = (difmag + 5) / 5
dp = 10 ^ logadp
PRINT "Distancia parsecs="; dp
INPUT "Correccion por nueva dist.pars (2) FINALIZAR(1)"; ndp
IF ndp = 1 THEN 1000
GOTO 11


1000 PRINT
INPUT "Si nuevo calculo 1"; nc
IF nc = 1 THEN 10
END