Формула Винценти для обратной геодезической задачи

Я запрограммировал решение обратной геодезической задачи, следуя статье Винценти. Когда стал тестировать, обнаружил интересную вещь. Формула отлично работает, но лишь до тех пор, пока обе точки не окажутся ровно на экваторе (обе широты равны нулю). В этом случае в ответ мы получаем nan (я проследил по коду, там возникает деление на нуль). Причём если хотя бы одна из широт отлична от нуля (0.00001 достаточно), то всё нормально считается.

Вопрос - я что-то неправильно запрограммировал, или же эти формулы действительно применимы лишь если хотя бы одна из широт ненулевая? Возможно, я что-то пропустил в статье, но найти каких-либо указаний на требования к широтам, помимо того, что нельзя на вход подавать антиподальные точки, не удалось.