小球反弹演示 MATLAB

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