1 CLS
2 PRINT "SADEYA/INFORTECNICA"
PRINT
PRINT "SUPLEMENTO ESTELAR"
PRINT "CALCULO DE LA VELOCIDAD TANGENCIAL Y REAL DE UNA ESTRELLA"
PRINT
10 INPUT "Velocidad radial Km/s"; vr
INPUT "Constante  8-14"; kkk


INPUT "Movimiento propio en A.R. ,,"; xx
INPUT "Movimiento propio Dec.,,"; yy
INPUT "Distancia en parsecs"; dp
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)


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 2"; ndis
IF ndis = 1 THEN dp = dispar
IF ndis = 1 THEN 11

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; anb2


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