function time = falling_time( y, yd, xs, xe )
% �~‰ºŽžŠÔ
  f = @(x) sqrt((1+yd(x).^2)./y(x));
  time = integral(f, xs, xe, 'RelTol', 0, 'AbsTol', 1e-12);
end
