моделирование движения спутника

Ну, моделирование - понятие широкое.
В простейшем случае (круговой орбиты вокруг объекта, находящегося в центре полярной системы координат):
theta := w*t;
r := r0;
где:
r,theta - полярные координаты,
r0 - расстояние от спутника до центра объекта, вокруг которого он вращается,
w - угловая скорость вращения,
t - время.

Значит так. Рисуем circle(320,240), цвет там куда надо и радиус 100 - типа планета:), потом цикл, где счётчиком будет аргумент.Rad:=200. Координата X:=rad*cos(arg), Y:=rad*sin(arg) и, значит, спутник(маааленький такой, типа далеко) вот и процедура Нарисовать_Спутник(X,Y:Word). :)

А в космосе тела не по эллиптическим орбитам вертятся?

Тогда так: X:=A*cos(arg), Y:=B*sin(arg) - A и B - полуоси эллипса.

Просто я помню, что в астрономии были законы Кеплера, из которых, в свою очередь, вытекают законы движения Ньютона... А сам спутник вокруг своей оси (осей) вращаться не может? А масса спутника не важна?
Не исключаю, что могу излишне наворачивать, но глупо предлагать решение задачи небесной механики, не имея никаких исходных данных..

мне кажется Der Meister прав! тут нужны законы астрономии! спутник постепенно будет приближаться к земле, а след-но радиус тоже менятся будет!

Ну ппц. Все физику позабыли? Пишем 2 З.Н. для планеты и для спутника и решаем ОДУ второго порядка любым численным методом. Задача школьная, реализаций полно.

1. Численное решение дифференциальных уравнений находится явно за пределами школьного курса.
2. При решении дифф.ур. численными методаим ВСЕГДА вносятся дополнительные пограшности. Только иногда эти погрешности бывают допустимыми, иногда - нет. Что касается данной задачи, то, во-первых, вероятнее всего эти погрешности недопустимы (т.к. будут приводить к изменению параметров орбиты), а во-вторых, применение численных методов решения ДУ для данной задачи вообще неэффективно.
Опять же, дурной тон решать численными методами то, что имеет аналитическое решение.