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) datos complementarios (3) FINALIZAR(1)"; ndp
IF ndp = 1 THEN 1000
IF ndp = 3 THEN 2000
GOTO 11


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

2000 PRINT
PRINT "DATOS COMPLEMENTARIOS"
DIFMA = 15 - MAB
DIFEX = DIFMA / 2.666666667#
EXPL = -4 + DIFEX
LUMI = 10 ^ (EXPL)
PRINT "LUMINOSIDAD SOL=1"; LUMI

MASA = LUMI ^ (1 / 3.3)
PRINT "MASA SOL=1"; MASA
TEMP = 9202.212 * 10 ^ (-.6646191 * ICOL) + 2500
IF ICOL < -.12 THEN TEMP = (ICOL * -10) * 10000



PRINT "TEMPERATURA K."; TEMP
TVID = 1 / (MASA * MASA)
PRINT "VIDA ESTELAR SOL=1"; TVID
PRINT
GOTO 1000