Расчёт текущей фазы и времени ближайшего восхода/заката Луны

Какие есть формулы для сабжа если известны текущие географические координаты, дата и время (UTC, поясное либо местное, не суть важно)? Можно ли, исходя из уже известного времени некоторого события восхода и заката Луны в прошлом, более-менее точно рассчитать время последующих восходов и закатов на довольно продолжительный период, например, пару сотен лет? Гуглом нашёл лишь расчёт текущей фазы в условных целочисленных попугаях от 0 до 7 (первая программа), но хотелось бы побольше градаций, например 29 по числу дней в синодическом месяце, да и как там с точностью непонятно, в конце описывается какой-то глюк с полной Луной.