2 asec - inverse trigonometric secant
8 x real, with absolute value >= 1
9 eps nonzero real, defaults to epsilon()
14 Returns the asec of x to a multiple of eps with error less in
15 absolute value than .75 * eps.
17 v = asec(x) is the number in [0, pi] for which sec(v) = x.
20 > print asec(2, 1e-5), asec(2, 1e-10), asec(2, 1e-15), asec(2, 1e-20)
21 1.0472 1.0471975512 1.047197551196598 1.04719755119659774615
27 NUMBER *qasec(NUMBER *x, NUMBER *eps)
30 asin, acos, atan, acsc, acot, epsilon