Determinacion de cuspides. Sistema de Morinus. Programa en Basic

'morinus
cls

pi=3.14159296535
cdr= pi/180

obl=23.45*cdr

dim z(12)
dim p(12)
dim v(12)
dim x(12)
dim y(12)
dim l(12)
dim signo(12)
dim grado1(12)
dim minuto(12)

lat=43*cdr
st=226.02
for n=1 to 12
 z(n)=st+60 + n*30
 if z(n)>360 then z(n)=z(n)-360
 l(n)=(atn(sin(z(n)*cdr) * cos(obl) / cos (z(n)*cdr)))/cdr
 if n<4 then l(n)=l(n)+180
 if n>3 and n<10 then l(n)=l(n)+270
 if l(n)>360 then l(n)=l(n)-360
 if l(n)<0 then l(n)=l(n)+180
 signo(n)=int(l(n)/30)+1
 grado=l(n)-30*(signo(n)-1)
 grado1(n)=int(grado)
 minuto(n)=int((grado - grado1(n))*60)
next


for n=1 to 12
  print using "## \\ ###,.####  \\ ###,.####  ##(##)##"_
  ; n;"  p";p(n);"  l";_
        l(n);grado1(n);signo(n);minuto(n)
next