Игра "Солнечная система" Опции

[Account]

Это не совсем игра, но так написано в моем задании.
Вот что мне надо постараться создать:

Моделируется движение планет в солнечной системе и прохождение комет через нее. В упрощенной модели не следует соблюдать реальные пропорции размеров планет, Солнца и расстояний между ними. В данной модели все планеты движутся вокруг Солнца в одной плоскости но с разной скоростью. У планет, Солнца и комет есть определенные размеры и масса, пропорциональная кубу линейного размера. В одном из углов "Солнечной системы" в случайные моменты времени генерируются кометы, летящие в произвольном направлении. Комета и планета престают существовать при столкновении. При столкновении кометы о Солнцем исчезает только комета. На кометы действуют силы притяжения планет и Солнца в соответствии с законом всемирного тяготения, что изменяет их траекторию. Комета вышедшая за пределы "Солнечной системы* исчезает на всегда. У комет есть хвост, направленный в сторону от Солнца, Управление моделью : при помощи клавиш: начать, закончить.

Итак, делаться будет на с++, с помощью псевдографики как в этой теме.
Как изначально представляю себе это. Будет базовый класс фигура, который будет иметь такие показатели. как радиус планеты, масса, сила притяжения, радиус орбиты. Каждая планета наследник от этого класса.
Не понял вот этого в задании

масса, пропорциональная кубу линейного размера

Так же пока не представляю как в случайном какого-либо конца экрана запустить комету в каком то направлении, рисовать ей хвост(какми способом) от солнца, ну и физику притяжения с изменением траектории.

Как наверстаю начальный код по планетам выложу, далее нужна будет конечно помощь. Заранее благодарю всех кто откликнется, за советы и помощь.

Сообщение отредактировано: Account - 17.07.2011 21:53

[IUnknown]

Только не пойму вот что, в чем вот разница объявлений(присваивание 0),

virtual void Show(int) = 0;
virtual void Recalc()

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

Забыл еще спросить по частицам, массив размерностью определенный переменный const int parrSize = 1500;
а используется только 100 ?

maxParticles(100)

Ага... Тут дело вот в чем: если взять просто массив из 100 частичек, инициализировать их, и потом, когда частички удалятся от кометы, переинициализировать опять рядом с кометой - то будет очень своеобразное поведение? то есть, хвост будет "фонтанировать": он будет создаваться, потом частички будут удаляться все дальше и дальше (при этом ни одна из них еще не уйдет слишком далеко, следовательно, новых-то частиц рядом с кометой не будет), пока не начнут уходить из области влияния комету, тогда она начнет испускать новую порцию частичек. Я попробовал - выглядит ужасно. Поэтому я решил пойти другим путем: я на каждой итерации добавляю по StepParticles - Count новых частичек рядом с кометой (где Count - счетчик удалившихся старых частиц), чтобы выброс шел постоянно. Естественно, что общее число частиц начинает возрастать. И возрастает до тех пор, пока длина хвоста не установится, потом возрастание прекращается. Так вот, для этого и нужно место под 1500 частичек, чтобы массив не переполнялся. Хотя, возможно, это я чересчур загнул, надо гонять программу и проверять, сколько максимально частичек присутствует в хвосте. Хочешь - сделай такую проверку. Мне лень было, я выделил место под 1500 элементов, так сказать, с запасом. Все равно, если б я делал для себя, я б не стал использовать массив, написал бы класс связного списка, и все... Сколько засунешь в него - столько он и будет хранить. Новые будут прибывать, старые - удаляться... Но тебе, чтоб не перегружать программу, я этого делать не стал.

Ну и нужен совет как дальше делать, столкновение кометы с планетами и солнцем.

Я тебе уже говорил выше: проверяй комету с каждым телом, если Расстояние(Комета, НебесноеТело) меньше чем R_Комета + R_{НебесноеТело}, значит, произошло столкновение... Но тебе до этого надо еще принимать меры, когда комета слишком близко подходит к планете - это будет менять ее траекторию. Принцип - то же самый: проверять на каждом шаге расстояние, и если оно меньше какого-то расстояния, вычисленного по массе планеты (там Шека тебе чего-то писал) - то менять Vx и Vy кометы...

Попробовал на нетбуке, тоже вылет, пишется Abnormal program termination.

Значит, запускай трассировку, и смотри, на какой именно строке программа вылетает. Возможно для этого придется перед каждой строкой в Comet::Recalc, Comet::Show и всех вызываемых оттуда функциях выводить тестовое сообщение, и потом смотреть, какое именно из них вывелось последним, значит, на последующей строке произошел вылет. Возможно - будет достаточно просто запустить отладку программы (F7 в IDE), поставить курсор на строку, скажем, с closegraph(), и нажать F4. При вылете программы курсор может быть установлен на строку с ошибкой. Так что ищи, где ошибка, коли она только у тебя воспроизводится...

[Account]

IUnknown, а можешь скинуть свой компилятор, на какой-нито файлообменник?
В нем проверю.
По F4 ошибок не выдает, по F7 что то при создании частиц слишком долго) но постараюсь сделать

Сообщение отредактировано: Account - 23.07.2011 19:44

[IUnknown]

Вот тут мой компилятор валяется: http://vlady.uzelok.net/progs/Tc30.exe
Просто запускай его, он сам распакуется куда скажешь...

[Account]

Вот тут мой компилятор валяется: http://vlady.uzelok.net/progs/Tc30.exe
Просто запускай его, он сам распакуется куда скажешь...

Спасибо. Блин вылет имеется. Буду пробовать по F7.

[IUnknown]

   for(i = 0; i <= stepParticles - count; i++)
   {
      // Вот отсюда ...
      if(maxParticles + i >= parrSize)
      {
         cout << "Oops..." << endl;
      }
      // ... до сюда
      trace[maxParticles + i] = new Particle((int)(getX()) - 1 + random(2), (int)(getY()) - 1 + random(2),
				    Vx + Cvt((int)((getX() - sun->getX()) / Ex)) - 1 + random(2),
				    Vy + Cvt((int)((getY() - sun->getY()) / Ey)) - 1 + random(2));
   // ...

Попробуй добавить вот этот кусок, поставить на вывод Oops курсор и нажать на F4... Если оно туда придет, значит, как я и предполагаю, у тебя по какой-то причине генерируется слишком много частичек и происходит вылет за границы массива. Если нет - то я вообще не понимаю, в чем дело, больше просто негде ошибаться... Вообще-то обычно есть сообщение, по какой причине Abnormal programm termination. Либо Floating point error: Domain, либо еще что-то. Почему у тебя этого не показывает - непонятно...

Тебе сейчас главное - локализовать место, где происходит вылет. Через F7 это на самом деле будет долго, для первого раза выполняй методы не заходя в них, то есть, по F8 на имени метода... Потом, когда будешь знать, в каком именно методе происходит вылет - зайдешь в него и там уже жми на F8... Так каждый раз будешь подбираться все ближе и ближе, пока не найдешь, что именно глючит... Найдешь - говори...

[Account]

код (Показать/Скрыть)

   
for(i = 0; i <= stepParticles - count; i++)
   {
      // Вот отсюда ...
      if(maxParticles + i >= parrSize)
      {
         cout << "Oops..." << endl;
      }
      // ... до сюда
      trace[maxParticles + i] = new Particle((int)(getX()) - 1 + random(2), (int)(getY()) - 1 + random(2),
				    Vx + Cvt((int)((getX() - sun->getX()) / Ex)) - 1 + random(2),
				    Vy + Cvt((int)((getY() - sun->getY()) / Ey)) - 1 + random(2));
   // ...

Попробуй добавить вот этот кусок, поставить на вывод Oops курсор и нажать на F4... Если оно туда придет, значит, как я и предполагаю, у тебя по какой-то причине генерируется слишком много частичек и происходит вылет за границы массива. Если нет - то я вообще не понимаю, в чем дело, больше просто негде ошибаться... Вообще-то обычно есть сообщение, по какой причине Abnormal programm termination. Либо Floating point error: Domain, либо еще что-то. Почему у тебя этого не показывает - непонятно...

Тебе сейчас главное - локализовать место, где происходит вылет. Через F7 это на самом деле будет долго, для первого раза выполняй методы не заходя в них, то есть, по F8 на имени метода... Потом, когда будешь знать, в каком именно методе происходит вылет - зайдешь в него и там уже жми на F8... Так каждый раз будешь подбираться все ближе и ближе, пока не найдешь, что именно глючит... Найдешь - говори...

Итак сидел и держал F8, вроде на апдейте кометы, подозреваю на ее recalc(), НО хочу дать потестить кому-ни то из знакомыx, так же было бы не плохо если бы кто-ни то из здешних у себя запустил и протестил. Вдруг может что из-за видео дров, у меня тут периодически вылет их был, когда возвращался из IDE (дос режима)

[Account]

Итак, дал погонять еще двоим, у одного вылет через 3 секунды, у другого не так быстро но также вылет, с тем же сообщением что я выше писал(((

Что то в комете. Буду искать по возможности.

[IUnknown]

Ау!!! Почему у меня (ноут с WinXP SP3, запускается через DosBox 0.74, мне на фиг не нужно смотреть на это все в полноэкранном режиме) нет никаких вылетов при том же коде, работающем с кометой ( видео прилагается: http://www.youtube.com/watch?v=PzSfcH3jBEE )? Это что, карма?

[Account]

Да уж, возможно и карма, а скинь досбокс, я у себя через него запущу. Так же если что установлю uvScreenCamera_4.6.0.100, если будет ошибка попробую записать.

[IUnknown]

Где-то вот тут: http://www.dosbox.com/download.php?main=1

[Account]

Где-то вот тут: http://www.dosbox.com/download.php?main=1

Все запустил)) Вот скрин и ВИДЕО

Сообщение отредактировано: Account - 24.07.2011 0:36

[IUnknown]

Не понимаю... Причины такого поведения программы - не понимаю. Знаешь, почему? Потому что я, чтоб было проще отлаживать, безо всяких DosBox-ов, даже без Windows, переписал программу с использованием GLUT (там переделывать-то ничего не пришлось, только чуть-чуть перегруппировать методы и добавить условную компиляцию для двух графических примитивов - для putpixel и для circle - типа, это для ДОС, а это для glut) и запустил ее под Code::Blocks у себя на компьютере (не на ноуте, а именно на десктопе). Ты будешь смеяться, но оно не вылетает. Если уж и GCC ошибки не ловит - я не знаю, что еще делать...

[Account]

Не понимаю... Причины такого поведения программы - не понимаю. Знаешь, почему? Потому что я, чтоб было проще отлаживать, безо всяких DosBox-ов, даже без Windows, переписал программу с использованием GLUT (там переделывать-то ничего не пришлось, только чуть-чуть перегруппировать методы и добавить условную компиляцию для двух графических примитивов - для putpixel и для circle - типа, это для ДОС, а это для glut) и запустил ее под Code::Blocks у себя на компьютере (не на ноуте, а именно на десктопе). Ты будешь смеяться, но оно не вылетает. Если уж и GCC ошибки не ловит - я не знаю, что еще делать...

У меня есть стационарный комп, но на нем стоит фряха, проверить пока больше нигде не могу. Как представиться возможность где-нито еще погонять. отпишусь по результатам. А пока займусь столкновением.
Я так понимаю, как ты предлагал что планеты буду как массив типа класса планета, чтобы можно было в цикле потом сравнить растояние каждой с кометой, только вот вычисления эти производить в Run() или как лучьше?

[Lapp]

[offtop]
Посмотрел видео, а после него, как обычно, выскакивает еще несколько "похожих" видео. Среди них песня некоей CINDY SANCHEZ, COMETI UN ERROR.avi. Меня привлекло слово cometi, которое проассоциировалось у меня с кометами, и я послушал.. Оказывается, название песни переводится (с испанского) как "я делаю ошибку" . Я только потом понял, что слово "комета" действительно похоже на англ. "to commit" - совершать. Интересно также, что в испанской его версии только одно m, и за ним следует e, что делает его еще больше похожим на русское "комета". Я не знаком с этимологией слова "комета" (может, кто-нить что-нить нароет?) - возможно, тут есть смысловая связь.. Но! Но я не нашел никаких упоминаний о кометах в видео Account'а.. Что-это - совпадение? Если да, то сразу по нескольким пунктам (ибо есть возможность, что ошибка все же делается..)) Не знаю уж, есть ли тут какой-то скрытый смысл
[/offtop]

[Lapp]

Кометы выглядят замечательно.. Но я не понимаю, почему они не притягиваются ни к Солнцу, ни к планетам? Некоторые проходят прямо через Солнце - и ничего! )
Притяжение еще не реализовано?

[IUnknown]

Притяжение еще не реализовано?

Нет еще, ищем пока ошибку, которая приводит к вылету программы у Account-а.

Присоединяйся, если есть желание У тебя ж есть компилятор который умеет работать с GLUT-ом? Присоединяю получившийся у меня "комбайн", который компилируется и в GCC, и в Турбо-Сях, и ведет себя одинаково, ибо отличается только сама отрисовка. Где может сбоить - ума не приложу.

 main.cpp ( 8.89 килобайт ) Кол-во скачиваний: 589

(сохранено в UTF8)

[Lapp]

Присоединяйся, если есть желание У тебя ж есть компилятор который умеет работать с GLUT-ом?

Спасибо - боюсь, поздновато, но я попробую (кстати, с той эмуляцией bgi-графики я не особо преуспел, хотя и пытался - в лом стало тратить на это время). GLUT не пользовал, сейчас пробую SDL - его должно с лихвой хватать для 2D.

Вечером посижу, попробую вникнуть..

Добавлено через 12 мин.

Нет еще, ...

Мне это немного странно.. Я бы с этого начал, потому что если его сделать - ВСЕ объекты будут двигаться единым способом. Ведь, если посмотреть с общей точки зрения - это просто система нескольких гравитирующих тел. В общем случае их движение представляет полный хаос, но в рассматриваемом случае (супермассивная звезда, массивные планеты, легкие кометы) мы просто знаем (глядя на небо)), что порядок возможен . То есть, не нужно отдельно делать движение планет по кругу - просто нужно задать подходящую массу и вектор скорости. Проблема будет только в том, как "отладить" всю систему, чтоб она не разваливалась из-за межпланетного взаимодействия (впрочем, вот его-то можно и "выключить" для простоты)).

[Account]

IUnknown, сделал вот что. Взяв за пример функцию dist() вот вставил вот такую boom() для проверки расположения кометы и каждой планеты

int boom(Celestial &a, Celestial &b)
{
   double dx =abs (a.getX() - b.getX());
   double dy =abs( a.getY() - b.getY());
   if ((dx<(a.getSize()+3)) || (dy<(a.getSize()+3)))
   return 1;
   return 0;
}

Далее переправил такие места

//.....
Planet *planets[2];

//....
planets[0]->Update();
planets[1]->Update();
//....
planets[0] = new Planet(maxx / 2.+ 3 * 20, maxy / 2., 3,60, 360./365., "earth");
planets[1] = new Planet(maxx / 2.+ 4 * 20, maxy / 2., 3,80, 360./687., "mars");
//.....

Пока думаю как лучше организовать проверку на столкновение и выход кометы за систему, ведь в одном случе уничтожаются два объекта в другом только комета и при обоих комета заново создается.

Сообщение отредактировано: Account - 24.07.2011 17:49

[IUnknown]

но в рассматриваемом случае (супермассивная звезда, массивные планеты, легкие кометы) мы просто знаем (глядя на небо)), что порядок возможен . То есть, не нужно отдельно делать движение планет по кругу - просто нужно задать подходящую массу и вектор скорости.

Да, да... Теоретически - так оно и есть. Как только дело доходит до практики... Попробовал я реализовать то же самое с использованием сил гравитации. Вот так и есть, реальная орбитальная скорость, реальная масса Солнца и планет, реальное расстояние между планетами и Солнцем + масштабирование, чтоб не было необходимости в мониторе с диагональю около 30 а.е.

Так вот, взаимодействие Солнца и планет отрабатывают на ура, все движется, да не по окружностям, а по эллиптическим орбитам. Красота... Вот, думаю, сейчас добавлю сюда комету, и посмотрю, как она будет двигаться. Так вот комете (при реальной же массе, которая на 16 порядков меньше массы Солнца - скажем, масса кометы Галлея = 2.2×10¹⁴ кг. против 2.0×10³⁰ кг. солнечной массы) крышу сносит на первой же секунде, её солнце притягивает и выстреливает ей куда-то в неизвестном направлении как из пращи, со скоростью порядка 1000 км/сек (при том, что реальная орбитальная скорость Меркурия, как самой быстрой планеты, не превышает 50 км/сек). Начинаю увеличивать массу кометы - все постепенно начинает приходить в норму, скорость выброса из системы уменьшается, а когда масса кометы становится соизмерима с массами планет - то комета переходит на гелиоцентрическую орбиту... Я явно что-то упускаю, а вот что - пока сообразить не могу...

[Account]

Вот что пытался сделать для столкновения, но что то не работает так как надо

код (Показать/Скрыть)

#define MSGS_OFF
#include <iostream.h>
#include <math.h>
#include <graphics.h>
#include <stdlib.h>
#include <conio.h>
#include <dos.h>

int maxx;
int maxy;
class Celestial
{
protected:
   double X, Y;
   int Size;
   virtual void Show(int) = 0;
   virtual void Recalc()
   {
   }
public:
   Celestial(double px, double py, int psize);
   double getX() const { return X; }
   double getY() const { return Y; }
   int getSize() const { return Size; }
   void Update();
};
Celestial :: Celestial(double px, double py, int psize) :
	     X(px), Y(py), Size(psize)
{
#ifdef MSGS
   cout << "Celestial created (" << X << " " << Y
	<< " " << Size << " )" << endl;
#endif
}
void Celestial :: Update()
{
   Show(0);
   Recalc();
   Show(1);
}
double dist(Celestial &a, Celestial &b)
{
   double dx = a.getX() - b.getX();
   double dy = a.getY() - b.getY();
   return sqrt(dx * dx + dy * dy);
}

int boom(Celestial &a, Celestial &b)
{
   double dx =abs (b.getX()-a.getX());
   double dy =abs(b.getY()-a.getY());
   if ( (dx<(a.getSize()/2+b.getSize()/2)) || (dy<(a.getSize()/2+b.getSize()/2)) )
   return 1; else return 0;
}
class Star : public Celestial
{
   int Shown;
public:
   Star(int psize) : Celestial(maxx / 2, maxy / 2, psize),
                     Shown(0)
   {
   #ifdef MSGS
      cout << "star created" << endl;
   #endif
   }
   virtual void Show(int);
};
void Star :: Show(int bShow)
{
   if(!Shown)
   {
      setcolor( bShow ? (Shown = 1, YELLOW) : getbkcolor() );
      circle((int)X, (int)Y, Size);
   }
}
Star *sun;
class Planet : public Celestial
{
private:
   int R;
   const char *name;
   double Angle, DeltaAngle; //<----Изменения, естественно в конструкторе тип поменял

protected:
   virtual void Recalc()
   {
      Angle -= DeltaAngle;
      if(Angle < 0) Angle += 360;

      X = maxx/2 + R*cos(Angle*M_PI/180); //<----Изменения
      Y = maxy/2 + R*sin(Angle*M_PI/180); //
   }

public:
   Planet(double px, double py, int psize,
	  int pr, double pdeltaangle, const char *pname);
   virtual void Show(int);
};
Planet :: Planet(double px,double py, int psize,
		 int pr, double pdeltaangle, const char *pname) :
		 Celestial(px, py, psize),
		 R(pr), Angle(360), name(pname),
                 DeltaAngle(pdeltaangle)
{
   #ifdef MSGS
      cout << "planet " << name << " created ("
	   << R << " " << Angle << " )" << endl;
   #endif
}


void Planet :: Show(int bShow)
{
   setcolor( bShow ? LIGHTGREEN : getbkcolor() );
   circle((int)X, (int)Y, Size);
}
class Particle : public Celestial
{
private:
   double Vx, Vy;
public:
   Particle(double px, double py, double pvx, double pvy);
   ~Particle();
   virtual void Show(int);
   virtual void Recalc();
};
Particle :: Particle(double px, double py, double pvx, double pvy) :
	    Celestial(px, py, 30),
	    Vx(pvx), Vy(pvy)
{
   #ifdef MSGS
      cout << "particle created : Vx = " << Vx << " Vy = " << Vy << endl;
   #endif
}
Particle :: ~Particle()
{
   Show(0);
}
void Particle :: Show(int bShow)
{
   putpixel( (int)X, (int)Y, bShow ? LIGHTGRAY : getbkcolor() );
}
void Particle :: Recalc()
{
   Vx += ((getX() - sun->getX()) / dist(*sun, *this));
   Vy += ((getY() - sun->getY()) / dist(*sun, *this));
   X += Vx;
   Y += Vy;
}
const int parrSize = 1500;
const int stepParticles = 15;//1;
const int maxPartSpeed = 3;
class Comet : public Celestial
{
private:
   double Ex, Ey;
   int Vx, Vy;
   Particle *trace[parrSize];
   int maxParticles;
private:
   int Cvt(int curr);
protected:
   virtual void Recalc();
public:
   Comet(double px, double py, int psize, int pvx, int pvy);
   virtual void Show(int);
 int OutOfSystem()
   {
      return (X < -5 || X > maxx + 5 || Y < -5 || Y > maxy + 5) ? 1 : 0;
   }
};
Planet *planets[2];
Comet *c_01;
int Comet :: Cvt(int curr)
{
   int sign;
   if(!curr) sign = 0;
   else sign = curr / abs(curr);
   return sign * random(maxPartSpeed - abs(curr));
}
Comet :: Comet(double px, double py, int psize, int pvx, int pvy) :
	 Celestial(px, py, psize),
	 Vx(pvx), Vy(pvy), maxParticles(1)
{
   Ex = (sun->getX() - sun->getSize()) / maxPartSpeed;
   Ey = (sun->getY() - sun->getSize()) / maxPartSpeed;
   for(int i = 0; i < maxParticles; i++)
   {
      double dx = (px - sun->getX());
      double dy = (py - sun->getY());
      #ifdef MSGS
	 cout << "i = " << i << " Dx = " << dx << " Dy = " << dy << endl;
      #endif
      trace[i] = new Particle(px - 1 + random(2), py - 1 + random(2),
			      Vx + Cvt((int)(dx / Ex)) - 1 + random(2),
			      Vy + Cvt((int)(dy / Ey)) - 1 + random(2));
   }
   #ifdef MSGS
      cout << "comet created (" << Vx << " " << Vy << " )" << endl;
      cout << "dist to sun = " << dist(*sun, *this) << endl;
   #endif
}
void Comet :: Recalc()
{
   X += Vx;
   Y += Vy;
   #ifdef MSGS
      cout << "comet : dist to sun = " << dist(*sun, *this) << endl;
   #endif
   int count = 0;
   for(int i = 0; i < maxParticles; i++)
   {
      trace[i]->Recalc();
      double d = dist(*trace[i], *this);
      #ifdef MSGS
	 cout << "i = " << i << " Dist = " << d << endl;
      #endif
      if(d > 25 * Size)
      {
	 delete trace[i];
	 count += 1;
	 trace[i] = new Particle((int)(getX()) - 1 + random(2), (int)(getY()) - 1 + random(2),
				 Vx + Cvt((int)((getX() - sun->getX()) / Ex)) - 1 + random(2),
				 Vy + Cvt((int)((getY() - sun->getY()) / Ey)) - 1 + random(2));
      }
   }
   #ifdef MSGS
      cout << count << " new particles created" << endl;
   #endif
   for(i = 0; i <= stepParticles - count; i++)
   {
//--------------------------------------
 if(maxParticles + i >= parrSize)
      {
         cout << "Oops..." << endl;
      }
//---------------------------
      trace[maxParticles + i] = new Particle((int)(getX()) - 1 + random(2), (int)(getY()) - 1 + random(2),
				    Vx + Cvt((int)((getX() - sun->getX()) / Ex)) - 1 + random(2),
				    Vy + Cvt((int)((getY() - sun->getY()) / Ey)) - 1 + random(2));
   }
   if(stepParticles - count > 0)
      maxParticles += (stepParticles - count);
   #ifdef MSGS
      cout << maxParticles << " currently. Max = " << parrSize << endl;
   #endif
}
void Comet :: Show(int bShow)
{  setcolor( bShow ? WHITE : getbkcolor() );
   circle((int)X, (int)Y, Size);
   for(int i = 0; i < maxParticles; i++)
      trace[i]->Show(bShow);
}
class Observer
{
public:
   Observer();
   ~Observer();
   void Run();
};
Observer :: Observer()
{
   int gdriver = DETECT, gmode, errorcode;
   initgraph(&gdriver, &gmode, "../BGI");
   errorcode = graphresult();
   if(errorcode != grOk)
   {
      cerr << "Graphics error: " << grapherrormsg(errorcode) << endl;
      exit(-1);
   }
   maxx = getmaxx(); maxy = getmaxy();
   randomize();
}
Observer :: ~Observer()
{
   closegraph();
}
void Observer :: Run()
{
   while(1)
   {
      sun->Update();
      c_01->Update();
      if(kbhit())
      {
	 if(getch() == 27) break;
      }

      if (!(c_01->OutOfSystem())) // Проверяем, если вышли за пределы - то пересоздаем комету...
         { for(int i=0;i<2;i++)
         {if (boom(*planets[i],*c_01))
           {delete planets[i];delete c_01;c_01 = new Comet(rand()%maxx, maxy, 2, 5, -5);}
          else planets[i]->Update();}
         }
      else
      {
         delete c_01;
	 c_01 = new Comet(rand()%maxx, maxy, 2, 5, -5);
      }

      delay(50);
   }
}
int main()
{
   Observer obj;
   sun = new Star(10);
   planets[0] = new Planet(maxx / 2.+ 3 * 20, maxy / 2., 3,60, 360./365., "earth");
   planets[1] = new Planet(maxx / 2.+ 4 * 20, maxy / 2., 3,80, 360./687., "mars");
   c_01 = new Comet(rand()%maxx, maxy, 2, 5, -5);
   obj.Run();
   return 0;
}