% Generacion de valores del tiempo de caida
% valores de H y k
H=10:20:400;
N=length(H);
k=0.01:0.1:0.7;
M=length(k);

% precision en la determinacion de la raiz
epsi=0.001;
g=10;
valores=[];


for i=1:M
	for j=1:N
	   h=k(i)^2*H(j)/g;
	   w=k(i)*sqrt(2*H(j)/g);
	   a=1;
		while a>epsi
			b=(h+1)-exp(-w);
			a=abs(b-w);
			w=b;
		end
		valores(i,j)=w/k(i);	
 	end
end
plot(H,valores,'y')
hold on
plot(H,valores,'xy')
plot(H,sqrt(2*H/g),'o')
xlabel('altura de caida')
ylabel('tiempo de caida')
hold off