Площадь правофокусного сектора мне очень помогла найти книга Д.Е.Охоцимский, Ю.Г.Сихарулидзе ''Основы механики космического полёта''. Во 2-й главе ''Задача двух тел'' в §2.5 ''Связь времени с положением на орбите'' подробно выводится решение указанного ниже интеграла.
Для времени выводится решение такое:
t=(E-e*sinE)*a*sqrt[a/(GM)].
Е-эксцентрическая аномалия.
Для площади сектора отсюда следует:
S=(E-e*sinE)*a*b/2.
Но удобнее выражать площадь через угол v при фокусе(истинную аномалию).
v=arccos[(cosE-e)/(1-ecosE)].
Тогда площадь сектора равна
S=0.5b2[(a/b)arccos[(e+cosv)/(1+ecosv)]-(esinv)/(1+ecosv)].
Ещё легко (через табличный интеграл) можно вывести формулу площади центрального сектора эллипса
S=0.5a*b*arctg[(a/b)*tg(beta)],
где beta-угол между большой полуосью и радиус-вектором,проведённым из центра эллипса .