'casa9.bas
for n=1 to 21:read hs$,h11$,h12$,hs1$,hs2$,hs3$:next
' Hora Sid Casa X Casa I CASAXI CASA XII CASA II CASA III
data "00:00:00"," 0( 1)","26( 4)10"," 8.6( 2)","22.0( 3)","12.4( 5)"," 2.4( 5)"
data "00:03:40"," 1( 1)","26( 4)51"," 9.6( 2)","22.9( 3)","13.1( 5)"," 3.4( 6)"
data "02:03:08"," 3( 2)","18( 5)12","12.4( 3)","19.2( 4)"," 6.6( 6)"," 0.7( 7)"
data "04:16:27"," 6( 3)","11( 6)39","13.1( 3)","15.2( 5)"," 3.3( 7)"," 1.4( 8)"
data "05:55:38","29( 3)","29( 6)13"," 5.0( 5)"," 5.1( 6)","23.9( 7)","23.2( 8)"
data "06:00:00"," 0( 4)"," 0( 7)00"," 5.9( 5)"," 6.0( 7)","24.0( 7)","24.1( 8)"
data "06:04:22"," 1( 4)"," 0( 7)47"," 6.8( 5)"," 6.9( 6)","24.9( 7)","25.0( 8)"
data "08:12:54"," 1( 5)","23( 7)31"," 5.4( 6)"," 2.6( 7)","20.1( 8)","23.2( 9)"
data "10:23:35"," 4( 6)","16( 8)30"," 5.4( 7)","28.6( 7)","16.3( 9)","24.4(10)"
data "11:56:20","29( 6)"," 3( 9)09","26.8( 7)","16.9( 8)"," 7.1(10)","20.4(11)"
data "12:00:00"," 0( 7)"," 3( 9)50","27.6( 7)","17.6( 8)"," 8.0(10)","21.4(11)"
data "12:03:40"," 1( 7)"," 4( 9)30","28.3( 7)","18.3( 8)"," 8.9(10)","22.6(11)"
data "14:30:21","10( 8)"," 4(10)55"," 0.7( 9)","17.7( 9)","23.6(10)","10.0( 1)"
data "16:46:16","13( 9)","20(11)09"," 0.7(10)","12.9(10)","19.4( 1)","22.0( 2)"
data "17:55:38","29( 9)","27(12)27","17.2(10)","11.2(11)","15.8( 2)","10.5( 3)"
data "18:00:00"," 0(10)"," 0( 1)00","18.3(10)","12.7(11)","17.3( 2)","11.7( 3)"
data "18:04:22"," 1(10)"," 2( 1)33","19.5(10)","14.2(11)","18.8( 2)","12.8( 3)"
data "20:41:41"," 8(11)","12( 3)18"," 4.4(12)","17.5( 1)"," 1.8( 3)","19.0( 3)"
data "22:49:53","11(12)","12( 4)43","16.4( 1)"," 3.1( 3)","28.7( 4)","16.8( 5)"
data "23:56:20","29(12)","25( 4)30"," 7.4( 2)","21.1( 3)","11.7( 5)"," 1.7( 6)"
data "24:00:00"," 0( 1)","26( 4)10"," 8.6( 2)","22.0( 3)","12.4( 5)"," 2.4( 6)"
la$="51N00"
'..................
cls
e$="23:28:00"
pi#=4*atn(1)
ra#=2*pi#/360
'......................
hs=val(left$(hs$,2))+val(mid$(hs$,4,2))/60+val(right$(hs$,2))/3600
if hs=24 then hs=0
la=val(left$(la$,2))+val(right$(la$,2))/60
e#=val(left$(e$,2))+val(mid$(e$,4,2))/60+val(right$(e$,2))/3600
' ................
armc=hs*360/24 'Ascension Recta del Medio Cielo
MC=atn(tan((armc)*ra#)/cos(e#*ra#))/ra# 'Longitud del Medio Cielo
if MC<0 then MC=MC+360
MC=MC+180
if MC>=360 then MC=MC-360
if hs>=0 and hs=<6 then MC=MC-180
if hs>18 and hs=<24 then MC=MC+180
'................
'a=armc+90 'Ascension Oblicua del Ascendente
a=armc+120
if a= 90 then a= 90.0010
if hs=12 and a=270 then a=270.0026
if a=270 then a=270.0001
if a>=360 then a=a-360
AA= (atn((tan(la*ra#))/cos(a*ra#)))/ra#
H1= (atn(tan(a*ra#)*cos(AA*ra#)/cos((AA+e#)*ra#)))/ra# 'Longitud del Ascendente
if H1<0 and hs>12 then H1=H1+180
if H1<0 then H1=H1+360
H1=H1+180
if H1>=360 then H1=H1-360
if hs>=18 then H1=H1-180
senx=tan(e#*ra#)*tan(la*ra#)
cosx=sqr(1-senx^2)
tanx=senx/cosx
x=atn(tanx) ' x en radianes
tanpolo11=sin(2*x/3)*(1/tan(e#*ra#))
AA= (atn(tanpolo11/cos(a*ra#)))/ra#
H11= (atn(tan(a*ra#)*cos(AA*ra#)/cos((AA+e#)*ra#)))/ra# 'Longitud de la casa 11
if H11<0 then H11=H11+360
print "H11";H11
'..................
call trans(H1,H1$)
call trans(MC,MC$)
print "latitud: ";la$;" hora sideral: ";hs$
d$= "\ \ ###,.## \ \ \ \"
print using d$;"Medio Cielo ";MC;MC$;h10$
print using d$;"Ascendente ";H1;H1$;h1$
end
'..................
sub s: while inkey$="":wend: end sub
'..................
sub trans(x,x$)
s=int(x/30)+1
ga=x-(s-1)*30:g=int(ga)
m=cint((ga-g)*60)
if m=60 then m=0:g=g+1
if g=30 then g=0:s=s+1
if s=13 then s=1
x$=str$(g)+"("+str$(s)+")"+str$(m)
end sub
No hay comentarios:
Publicar un comentario