:h_az()
:Prgm
:Local le,de,ah,hc,ze
:Local lo,do,tp,az
:Local lns,dns
:©..................objet
:ClrIo
:Output 13,13,"calcule de h
 c & Z avec"
:Output 29,49,"L, d, LHA,"

:Output 45,13,"de la posit
 ion estimée"
:Pause
:©.................saisie
:ClrIO
:Loop
: 0->lns
: Input "latitude Nord est
 imée ?",le
: 1->lns
: If le<>0
: Exit
: Input "latitude Sud esti
 mée ?",le
: 2->lns
: If le<>0
: Exit
:EndLoop
:le->lo
:Loop
: 0->dns
: Input "déclinaison Nord 
 ?",de
: 1->dns
: If de<>0
: Exit
: Input "déclinaison Sud ?
",de
: 2->dns
: If de<>0
: Exit
:EndLoop
:de->do
:Input "angle horaire LHA 
 ?",ah
:©............Dieumegard
:If ah>180 then
: 360-ah->tp
:Else
: ah->tp
:EndIf
:If lns=2 and dns=1 
:360-le->le
:If lns=1 and dns=2
:360-de->de
:
:sin-1(sin(le)*sin(de)+(cos(
 le)*cos(de)*cos(tp))->hc
:
:If hc < 0
:abs(hc)->hc
:©...............Bataille
:
:cos-1((sin(de)-sin(hc)*sin
 (le))/(cos(hc)*cos(le)))->
 ze
:
:If ah<180 then
: 360-ze->az
:Else
: ze->az
:EndIf
:©..............affichage
:ClrIO
:Output -4,1,"latitude est
 imée"
:If lns=1 then
: Output 4,7,"N"
:Else
: Output 4,7,"S"
:EndIf
:Output 4,19,lo>DMS
:Output 12,1,"déclinaison"

:If dns=1 then
: Output 20,7,"N"
:Else
: Output 20,7,"S"
:EndIf
:Output 20,19,do>DMS
:Output 28,1,"angle horair
 e LHA"
:Output 36,19,ah>DMS
:Output 52,1,"hc"
:Output 52,19,hc>DMS
:Output 60,1,"Z"
:Output 60,19,az>DMS
:EndPrgm

// calcul de he (hauteur) et de az (azimut) d'un astre
// comme avec les tables de Dieumegard et Bataille
// les données de départ sont :
// latitude estimée (du point estimé de l'observateur)
// déclinaison de l'astre (des éphémérides)
// angle horaire de l'astre (des éphémérides)