Моделирование физики

[notHaker]

Познаю азы моделирования физики... Выполнил просчёт (грубый) гравитации, столкновения. Помогите пожалуйста с моделированием трения. Задача проста... Из точки "(x,y)" запускаю объект со скоростью "vi" под заданным углом "ugol". А далше...

procedure TForm1.Button1Click(Sender: TObject);
const
  g = 1; //"грубая" гравитация
  a = 0.3; //трение... сколько скорости потеряет объект при столкновении
var
  vi,ugol : integer;
  vix,viy : real;
  x, y : integer;
begin
  vi := StrToIntDef(Edit2.Text,10);
  ugol := StrToIntDef(Edit1.Text,30);
  vix := vi * cos(ugol/180 * pi); //скорость по координате x
  viy := vi * sin(ugol/180 * pi); //скорость по координате y
  x := 1; //начальные значения
  y := 1;
  while true do
  begin
    if stop then break;
    application.ProcessMessages;
    Canvas.Brush.Color := clWhite;
    Canvas.FillRect(rect(0,0,Self.ClientWidth,Self.ClientHeight));
    Canvas.Draw(x,y,bmp);
    if (x <= 0) or (x >= self.ClientWidth) then //моделирование столкновения
        vix := -vix;
    if (y <= 0) or (y >= self.ClientHeight) then
        viy := -viy;

    x := round(x + vix); //моделирование траектории
    y := round(y + viy);

    if viy < 0 then //моделирование гравитации
      viy := viy + g
    else
      viy := viy - g;

    sleep(25);
  end;
  stop := false;
end;

всё в порядке... но даже без трения, объект при некоторых значениях g и a, просто выпадает за пределы, предписанные ему... бред какой-то... а при значениях в этом листинге, всё в порядке... что за бред?

[notHaker]

откликнитесь... пожалуйста...

[evg_m]

//считаем куда он должен прийти
    x := round(x + vix); //моделирование траектории
    y := round(y + viy);

if x>=self.ClientWidth //выскочили за границу 
  // x+x.with>clientwith если объект имеет размеры
 then begin
  wix:=-wix; //отражаем скорость на будущее(меняем знак)
  x:=self.clientwith-(x-self.clientwith);//отражаем путь пройденный за границей еще здесь надо наверное учесть его размеры
end;

[notHaker]

Сообщение от evg_m

//считаем куда он должен прийти
    x := round(x + vix); //моделирование траектории
    y := round(y + viy);

if x>=self.ClientWidth //выскочили за границу 
  // x+x.with>clientwith если объект имеет размеры
 then begin
  wix:=-wix; //отражаем скорость на будущее(меняем знак)
  x:=self.clientwith-(x-self.clientwith);//отражаем путь пройденный за границей еще здесь надо наверное учесть его размеры
end;

ммм... он выпадает на совсем... это коррекция границ и её я понял, спасибо. НО, скорость абсолютно не зависит от координат, а она не принимает тех значений, которые возвратили бы объект обратно... Когда реализовал трение (просто "отнимал" скорость при столкновении), то вабще сумасброд при гравитации > 1.

[evg_m]

    if viy < 0 then //моделирование гравитации
      viy := viy + g
    else
      viy := viy - g;

считаем что viy>0 //движение вниз viy<0 //движение вверх
viy:=viy+g; всегда
ускорение/замедление зависит от совпадения/несовпадения знаков
viy и g

[Alien41180]

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