Потом, альфа у тебя а =137 градусов. Приведи к углу 0-90.
a>90, a<180. Угол во второй четверти  синус альфа равен sin(180-а), или нет? Не помню.

А по формулам совпадает, которые ты на рисунке сфоткал?

У меня все совпало. Вот твои цифры:

#include <stdio.h>
#include <math.h>

int main(void)
{
   double result, x, pi;
   pi=3.14;
   x= (3.07327 + (1.33617*sin((180-137.679165)*pi/180))*tan(14.390278*pi/180) )*29.5;
   printf("The x =  %lf\n", x);
   return 0;
}

x=97.4642!

1.у тебя переменные с двойной точностью -мой контроллер их не поддерживает.
2.  ты думаешь , что все углы надо загонять в диапазон от 0 до 90њ ?
ну это можно попробовать

Gute nacht, schlaff gut, bis Morgen

A = (math.cos(Declination)*math.sin(t))/(-math.cos(Latitude)*math.sin(Declination)
     + math.sin(Latitude)*math.cos(Declination)*math.cos(t));//OK

азимут  = math.atan(A);


   a = math.cos(Latitude)*math.cos(t)*math.cos(Declination);
   b = math.sin(Latitude)*math.sin(Declination);
       высота = math.asin(a+b);

если заданные координаты положительные, то расхождение с координатами в проге Старринайт от почти полного совпадения до примерно 20-30 минут.А вот если Declination отрицательная( напр. туманность Орион), то тут разница в несколько градусов!!!!
Как сие можно исправить?

вот и я не пойму, вроде все просто перевел в радианы и гото по формулам.А результатпочти никакой.С такой навигацией не долго и заблудиться в открытом космосе

Вопрос по синтаксису Си++

int& - это какой тип переменной?
D = int(x); - какое действие выполняет это?
D*=-1; - и вот это что такое?
//------------------------------------------------------------------------------
// DMS: Finds degrees, minutes and seconds of arc for a given angle
// Input:
//   Dd        Angle in degrees in decimal representation
// Output:
//   D         Angular degrees
//   M         Minutes of arc
//   S         Seconds of arc
//------------------------------------------------------------------------------
void DMS (double Dd, int& D, int& M, double& S)
{
  //
  // Variables
  //
  double x;


  x = fabs(Dd);   D = int(x);
  x = (x-D)*60.0; M = int(x);  S = (x-M)*60.0;
 
  if (Dd<0.0) { if (D!=0) D*=-1; else if (M!=0) M*=-1; else S*=-1.0; }
}

Народ, а что означает троечка возле второго тангенса?
и как произносится закорючка которая после этого тангенса стоит?

Закорючка похожа на ключ. Я таким стартер снимаю.