clear; clf
g=9.8;		% g en m/s2	
m=1; 		% masa en Kg
H=10;		% altura en metros

dt=0.001;	% dt en segundos

x0=0; v0=0;	% CI

x(1)=x0;
x(2)=x0+v0*dt;

N=ceil(sqrt(2*H/g)/dt)*10;  %estimacion de tiempo de caida y num de puntos

t=0:dt:(N-1)*dt;


ga=0.0;		%valor de gama

for i=2:N-1
  x(i+1)=g*dt^2+(2-ga/m*dt)*x(i)+(ga/m*dt-1)*x(i-1);
end;

ix=max(find((10-x)>0));		%encuentra momento de paso por h=10-x=0
tc=t(ix)			% tiempo anterior
vc=(x(ix+1)-x(ix-1))/(2*dt)	% vel de encuentro
