[[escafandra]]

Esa función ya te la había puesto aquí.

Empezado por escafandra

Por partes:
1- Según la formula que pones DA = ARCSIN(TAN LT * TAN DECL) pero eso no lo veo en el código que debería ser algo como:

Código Delphi [-]
var   DA, lt, decl : extended; 
begin   
  lt:=    -2.19474036998329;
  decl:= -24.7975765404927;
  DA:= RadToDeg(arcsin(tan(DegToRad(lt))*tan(DegToRad(decl))));
end;

.....

Saludos.