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.1:0.5;		%valores de gama

for j=1:length(ga)

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

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

end;