clear?all;
g=9.8;
v0=0;
y0=1.2;
m=0.4;
t0=0;
K=0.85;
n=700;
dt=0.005;
v=v0;
y=y0;
vx=0;
vz=0;
sx=0;
sz=0;

tic;
for?k=1:n
????if?y>0
????????v=v-g*dt;
????????y=y+v*dt;
????else
????????y=y-K.*v*dt;
????????v=-K.*v-g*dt;
????end
????Fx=r