How do I generate Mathieu functions in Matlab ?

The besselj function in Matlab has been a great tool for using bessel functions. It appears there are no equivalent tools for Mathieu functions.

I am trying to reproduce the figure below (Source: http://optica.mty.itesm.mx/pmog/Papers/P001.pdf), which is the even modified Mathieu function of 0-order as a function of q (for fixed epsilon_o = 0.693).

How can I go about achieving this on Matlab ? Ultimately, I would like to adapt the script for odd functions and higher orders.

(for context, s=0:0.1:15;plot(s,besselj(s,0.693)) would be all I'd need if it was a first order bessel function)