On the computation of Mathieu functions