First of all, thanks for such a robust approach.

In the eigen value there is "_Z1" which you later on replace by "n". What is this?

The eigen value we need to find has to be of the form.

lambda = (a^2+Pi^2*)^(3/2)/a instead of lambda = (a^2+Pi^2*_Z1^2)^(3/2)/a.

Finally, if I understood you correctly, in such situations (in maple for fast execution)

first we need to find the eigenvalue and then use that to solve the system of odes,

instead of going to find both (eigenfunction and eigenvalue) at the same time. 

Thanks

First of all, thanks for such a robust approach.

In the eigen value there is "_Z1" which you later on replace by "n". What is this?

The eigen value we need to find has to be of the form.

lambda = (a^2+Pi^2*)^(3/2)/a instead of lambda = (a^2+Pi^2*_Z1^2)^(3/2)/a.

Finally, if I understood you correctly, in such situations (in maple for fast execution)

first we need to find the eigenvalue and then use that to solve the system of odes,

instead of going to find both (eigenfunction and eigenvalue) at the same time. 

Thanks