Переписал выше написанный алгоритм на Delphi.
Осталось только исправить формулы в цикле и вывести переменную которую искали.
Код:
procedure TForm1.Button1Click(Sender: TObject);
Const
h = 500;
m = 65;
g = 9.8;
dt = 0.1;
Var
ti, xi, ai, ui: Real;
begin
ti:= 0; xi:= 0; ai:= g; ui:= 0;
Repeat
ti:= ti + dt;
xi:= xi + ui*ti + (ai*ti*ti)/2;
ui:= ui + ai*dt;
//ai:=(m*g - k1*u1 - k2*u2)/2; //Непонял тут ничё.
Until xi >= h;
end;