matlab program

euclude

ψ are eigenfunctions up to order 2 of the mathieu equation given by {(-d^2ψ/dx^2)-kcos(2x)}=Eψ with k=3.0
and E are the eigen values of the same mathieu equation up to order 2.

please want a matlab program that display eigenvalues and eigenfunctions up to order 2 of this mathieu equation. thank you