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 >20 Km/s., media 75Km/s"; kkk
kkk = kkk / 4.7

INPUT "Movimiento propio en A.R. arcosegundos ,,"; xx


INPUT "Movimiento propio Dec. arcosegundos ,,"; yy
INPUT "Paralaje en arcosegundos obervado o experimental"; parae
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"; maga
difmag = -5 + 5 * LOG(dp) / LOG(10)
mab = maga - difmag
PRINT "Difmag.="; difmag

PRINT "Magnitud absoluta transp="; mab
PRINT "Determinar si la posicion es aceptable en el diagrama Her. Rusell"
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