void discsFrame::OnTimer1Trigger(wxTimerEvent& event)
{ wxPaintDC z(Panel1);               // створення об'єкту z для малювання
  z.SetBrush(wxColor(0,0,0));        // призначення чорного кольору для заповнення
  z.SetPen(wxPen(wxColor(0,0,0),1)); // призначення чорного кольору для контуру ширини 1 піксель
  z.DrawRectangle(0,0,600,450);      // заповнення прямокутника з протилежними вершинами (0,0) і (600,450)
if (start)
  { srand(time(NULL));
    for (int j=0; j<n; j++)
    { r[j]  = rand()%10 + 10;
      vx[j] = rand()%11 - 5;
      vy[j] = rand()%11 - 5;
      c[j]  = wxColor(rand()%256, rand()%256, rand()%256);
      if (j==0) { x[j] = rand()%(int) (600-2*r[j]+1) + r[j];
                  y[j] = rand()%(int) (450-2*r[j]+1) + r[j];
                }
      else { bool   find = true;
             while (find)
             { x[j] = rand()%(int) (600-2*r[j]+1) + r[j];
               y[j] = rand()%(int) (450-2*r[j]+1) + r[j];
               int k=0;
               while ((find)&&(k<j-1))
               { find = (x[j]-x[k])*(x[j]-x[k])
                      + (y[j]-y[k])*(y[j]-y[k])
                      > (r[j]-r[k])*(r[j]-r[k]);
                 k++;
               }
               find = !find;
             }
           }
    }
    start = false;
  }
  else
  { for (int j=0; j<n; j++)
    { x[j]+=vx[j]; // рівномірний прямолінійний рух
      y[j]+=vy[j];
                  // пружне відбивання від меж
      if (x[j]<r[j]) {x[j]=2*r[j]-x[j];  vx[j]=-vx[j];}
      if (y[j]<r[j]) {y[j]=2*r[j]-y[j];  vy[j]=-vy[j];}
      if (600-r[j]<x[j]) {x[j]=2*(600-r[j])-x[j];  vx[j]=-vx[j];}
      if (450-r[j]<y[j]) {y[j]=2*(450-r[j])-y[j];  vy[j]=-vy[j];}
    }
    int j=0;
    while (j<n-1)
    { int k=j+1;
      while (k<n)
      { if ((x[j]-x[k])*(x[j]-x[k])
          + (y[j]-y[k])*(y[j]-y[k])
          <=(r[j]+r[k])*(r[j]+r[k]))
        { double rj = r[j]*r[j], // "маса" диску j
                 rk = r[k]*r[k], // "маса" диску k
                 rjk= rj+rk;     // "маса" результату "злипання"

          x[j] = (x[j]*rj+x[k]*rk)/rjk;   // "Злипання дисків"
          y[j] = (y[j]*rj+y[k]*rk)/rjk;   // з наданням властивостям диску j
          vx[j]= (vx[j]*rj+vx[k]*rk)/rjk; // значень властивостей результату злипання
          vy[j]= (vy[j]*rj+vy[k]*rk)/rjk;
          int red   = (int) ((rj*c[j].Red()   + rk*c[k].Red()  )/rjk);
          int green = (int) ((rj*c[j].Green() + rk*c[k].Green())/rjk);
          int blue  = (int) ((rj*c[j].Blue()  + rk*c[k].Blue() )/rjk);
          c[j] = wxColour(red, green, blue);
          r[j] = sqrt(rjk);

          x[k] = x[n-1];   // зменшення кількості дисків
          y[k] = y[n-1];   // з наданням властивостям диску k
          vx[k] = vx[n-1]; // значень властивостей диску (n-1)
          vy[k] = vy[n-1];
          r[k] = r[n-1];
          c[k] = c[n-1];
          k--;
          n--;
        }
        k++;
      }
      j++;
    }
  }
  for (int j=0; j<n; j++)
  { z.SetBrush(c[j]);                      // призначення кольору для заповнення круга j
    z.SetPen(wxPen(c[j],1));               // призначення кольору для контуру круга j ширини 1 піксель
    z.DrawCircle(wxPoint(x[j],y[j]),r[j]); // зображення круга j
  }
}