При увеличении радиуса, скорость падает, а расстояние все равно увеличивается - Pascal ABC

program ex3;
var T: Text;
x, h, a, Vx, Vy, x1, y1, V, k1_1, k1_2, k1_3, k1_4, k2_1, k2_2, k2_3, k2_4,
k3_1, k3_2, k3_3, k3_4, k4_1, k4_2, k4_3, k4_4: real;
function F(x, Vx, Vy: real; n: byte):real;
begin
case n of
1: f:=1*((0.005033+182.793569*sqrt(Vx*Vx+Vy*Vy))/50)*Vx;
2: f:=-9.81-((0.005033+182.793569*sqrt(Vx*Vx+Vy*Vy))/50)*Vy;
3: f:=Vx;
4: f:=Vy; 
end;
end;
begin
assign(T, 'D:\1.txt');
rewrite(T);
a:=0; h:=0.001; Vx:=6; Vy:=1; 
x1:=0; y1:=10;
x:=a;
while y1>=0 do
begin
writeln(T, x:10:5, V:10:5, x1:10:5, y1:10:5);
k1_1:=h*F(x,Vx, Vy, 1);
k1_2:=h*F(x,Vx, Vy, 2);
k1_3:=h*F(x,Vx, Vy, 3);
k1_4:=h*F(x,Vx, Vy, 4);
k2_1:=h*F(x+h/2, Vx+k1_1/2, Vy+k1_2/2, 1);
k2_2:=h*F(x+h/2, Vx+k1_1/2, Vy+k1_2/2, 2);
k2_3:=h*F(x+h/2, Vx+k1_1/2, Vy+k1_2/2, 3);
k2_4:=h*F(x+h/2, Vx+k1_1/2, Vy+k1_2/2, 4);
k3_1:=h*F(x+h/2, Vx+k2_1/2, Vy+k2_2/2, 1);
k3_2:=h*F(x+h/2, Vx+k2_1/2, Vy+k2_2/2, 2);
k3_3:=h*F(x+h/2, Vx+k2_1/2, Vy+k2_2/2, 3);
k3_4:=h*F(x+h/2, Vx+k2_1/2, Vy+k2_2/2, 4);
k4_1:=h*F(x+h, Vx+k3_1, Vy+k3_2, 1);
k4_2:=h*F(x+h, Vx+k3_1, Vy+k3_2, 2);
k4_3:=h*F(x+h, Vx+k3_1, Vy+k3_2, 3);
k4_4:=h*F(x+h, Vx+k3_1, Vy+k3_2, 4);
x1:=x1+1/6*(k1_1+2*k2_1+2*k3_1+k4_1);
y1:=y1+1/6*(k1_2+2*k2_2+2*k3_2+k4_2);
Vx:=Vx+1/6*(k1_3+2*k2_3+2*k3_3+k4_3);
Vy:=Vy+1/6*(k1_4+2*k2_4+2*k3_4+k4_4);
V:=sqrt(Vx*Vx+Vy*Vy);
x:=x+h;
end;
close(T);
end.