.... кроме того, в самой этой таблице даныые довольно странные- между полнолуниями время разное...

Ничего странного - лунная орбита не круговая, поэтому орбитальная скорость Луны (и видимая скорость перемещения Луны по земному небу) слегка меняется (видимая - от 11^o до 15^o в сутки), и разные полнолуния приходятся на разные точки лунной орбиты. К тому же и сами точки перигея и апогея смещаются, и в результате интервал между полнолуниями в разные месяцы получается неодинаковым - от 29д 6ч 15м до 29д 19ч 12м. А чтобы учесть все эти (и не только эти) эффекты, нужны довольно сложные вычисления. Гораздо проще воспользоваться готовой программой, благо их сейчас много, в том числе и бесплатных - можно за считанные минуты скачать из сети, установить и пользоваться. Например, есть замечательная программа StarCalc.

Уважаемый ilal, на странице http://www.moonphasecalendar.com/questions.htm приведена формула (дает моменты лунных фаз с точностью до 12 часов):

Ф = 0.20439731 + t x 0.03386319269 + t^2 x 0.0301 x 10^{-14}
t - время в сутках с момента 2001 января 1, 0 ч всемирного времени UTC
После расчета нужно отбросить от Ф целую часть, оставить только дробную:
Ф = 0 - новолуние
Ф = 0.25 - первая четверть
Ф = 0.5 - полнолуние
Ф = 0.75 - последняя четверть.

Число суток с 01.01.2001 проще всего подсчитать при помощи юлианских дней (JD, http://www.astronet.ru/db/msg/1173749).
01.01.2001 - JD = 2451911
Удобный способ расчета юлианских дней на любую дату - на странице http://m31.spb.ru/astronomy/time/julianday.shtm

При такой формуле фаза Луны +_12 часов, если нужно точнее, все сильно усложняется. Нужно брать громоздкие формулы для эклиптических долгот Солнца и Луны, например, из Астрономического Ежегодника. Там уже учитываются все эффекты, в том числе эллиптичность орбит Земли и Луны, а также движение линии узлов лунной орбиты.

^ это знак возведения в степень: t^2 значит t в квадрате. 10^{-14} десять в минус 14-й степени.