Algorithm 855: Subroutines for the computation of Mathieu characteristic numbers and their general orders

A continued fraction function algorithm is developed to evaluate general-order Mathieu characteristic numbers, and a new technique is presented for evaluating the Mathieu determinant which can be used to compute the order directly. Approximate expressions are developed to estimate the orders and Mathieu characteristic numbers for the root, finding algorithms. The algorithms, with minor modifications, were used for computing Mathieu coefficients of general order. The algorithms can deal with a large range of Mathieu characteristic number c, real and complex order., and parameter h.