This lecture presents a modern approach for the computation of Mathieu functions. The authors review the traditional approach used for these functions, show its limitations, and provide an alternative "tuned" approach enabling improved accuracy and convergence.