return (((1 + sqrt(1 + 8 * h)) / 4. / h) ** 2 - 1) / 16