Автоматическая коррекция времени в таймере ут1-pic

Подскажите, пожалуйста, по какой формуле вычисляется астрономическое время восхода(захода) солнца в опорную дату в местности с заданной географической широтой? Или, что еще лучше, приведите пример конкретного алгоритма для pic16. В руководстве по эксплуатации таймера Овен УТ1-pic (см. раздел Программирование прибора п.5.3.9.) сказано, что для местности с широтой 56 градусов будет вычислено, что 20 февраля восход произойдет в 7ч. 14мин. Напишите, пожалуйста, подробнее, как это вычисляется. В интернете я нашел только формулу для вычисления часового угла восхода и захода солнца: cos t = — tg ф tg d, где t собственно часовой угол, ф широта местности, d склонение светила. Но архитектура микроконтроллера pic16 не поддерживает подобные вычисления и из исходных данных известна только широта.