время восхода и захода солнца в c

В моем приложении на C я хочу рассчитать время восхода / захода солнца для заданной даты, широты и долготы. Я искал в сети, но не могу найти рабочий образец.

Я попытался реализовать этот образец: http://souptonuts.sourceforge.net/code/sunrise.c.html

Но этот образец работал некорректно.

Есть ли простой исходный код на C или метод, который я могу легко реализовать в своем приложении? ?

Редактировать:
Я использую код по этой ссылке , но это дало мне неправильные значения захода / восхода солнца. Также я попробовал ссылку Саула здесь , но она тоже дала мне неправильный результат.
У меня 41N, 28E локация. Когда я пробую коды, оба образца говорят, что значение восхода солнца составляет приблизительно 10:13, а закат - 23:24. Но правильные значения - 06:06, 20:13.
Не могу понять, в чем проблема.