Question: Plot Mathieu Solutions

Hi,

 

I am new to Maple and asked a question a few days a go and got a great answer (thanks for that!). I am getting more familiar with the software.

 

However, I am struggling to plot a 'stability diagram' with the Mathieu solutions. 

 

In an attempt to solve the equations from (link):

http://www.qmg700.com/quadinfo/Literature/Quadrupole%20basics%20mass%20filter%20neutral%20v4.

In Maple:

de1:=diff(x(t),t$2)+(a+2*q*cos(omega*t))*x(t)=0;
de2:=diff(y(t),t$2)-(a+2*q*cos(omega*t))*y(t)=0;
de3:=diff(z(t),t$2)=0;
sys:=[de1,de2,de3]:
sol:=dsolve(sys);

Now my issue is plotting the stability boundary curves in a,q space. Any ideas?

 

Can anyone please help?

 

Thanks

Please Wait...