Численное интегрирование в делфи

[Ира91]

Помогите кто-нибудь....
Искусственный спутник движется по замкнутой орбите вокруг Земли. Траектория спутника описывается уравнениями:
x(t)=a(cosξ-e),
y(t)=a√(1-e^2 ) sin ξ ,
t=1/R0 √(a^3/g)(ξ- esin ξ).
Здесь: ξ — параметр с пределами изменения от 0 до 2π;
a=(rmax + rmin)/2 - большая полуось орбиты;
e=(rmax- rmin)/(rmax+ rmin) - эксцентриситет орбиты;

rmin = 8000 км— перигей орбиты;
rmax = 20000 км — апогей орбиты;
R0 — радиус Земли;
g — ускорение свободного падения.

Построить зависимость длины траектории спутника
S=∫_0^2π Integral(√((dx/dξ)^2+(dy/dξ)^2 ) )dξ
от апогея орбиты rmax в пределах от 8000 до 200000 км.