Caída libre con rebotes

function rebote(elasticidad,yo)

v = 0;
dt = 1/10;
tf = 50;
g=-9.81;

times = 0:dt:tf;
x=0*times;
z=0*times;
y = zeros(size(times))

for idx = 1:length(times)
newy = yo + (v+g*dt/2)*dt;
v = v + g*dt;
if newy < 0
newy = 0;
v = -v*elasticidad;
end
y(idx) = newy;
yo = newy;
end

for i=1:length(times)
pause(1/100)  
plot3(x(i),z(i),y(i),'o r','MarkerFaceColor','m')
grid on
axis([-2 2 -2 2 min(y) max(y)])
end

end