Plot empty lattice modes in lat2d_realfreqsvd.py

[qpms.git] / besseltransforms / 7-4-7

((k^8*(-105 + 16*Sqrt[1 + k^2/(c - I*k0)^2])*(c - I*k0) + 8*k^6*(-105 + 44*Sqrt[1 + k^2/(c - I*k0)^2])*(c - I*k0)^3 + 144*k^4*(-14 + 9*Sqrt[1 + k^2/(c - I*k0)^2])*(c - I*k0)^5 + 320*k^2*(-6 + 5*Sqrt[1 + k^2/(c - I*k0)^2])*(c - I*k0)^7 + 640*(-1 + Sqrt[1 + k^2/(c - I*k0)^2])*(c - I*k0)^9)/(5040*k^7) - (k^8*(-105 + 16*Sqrt[1 + k^2/(2*c - I*k0)^2])*(2*c - I*k0) + 8*k^6*(-105 + 44*Sqrt[1 + k^2/(2*c - I*k0)^2])*(2*c - I*k0)^3 + 144*k^4*(-14 + 9*Sqrt[1 + k^2/(2*c - I*k0)^2])*(2*c - I*k0)^5 + 320*k^2*(-6 + 5*Sqrt[1 + k^2/(2*c - I*k0)^2])*(2*c - I*k0)^7 + 640*(-1 + Sqrt[1 + k^2/(2*c - I*k0)^2])*(2*c - I*k0)^9)/(720*k^7) + (k^8*(-105 + 16*Sqrt[1 + k^2/(3*c - I*k0)^2])*(3*c - I*k0) + 8*k^6*(-105 + 44*Sqrt[1 + k^2/(3*c - I*k0)^2])*(3*c - I*k0)^3 + 144*k^4*(-14 + 9*Sqrt[1 + k^2/(3*c - I*k0)^2])*(3*c - I*k0)^5 + 320*k^2*(-6 + 5*Sqrt[1 + k^2/(3*c - I*k0)^2])*(3*c - I*k0)^7 + 640*(-1 + Sqrt[1 + k^2/(3*c - I*k0)^2])*(3*c - I*k0)^9)/(240*k^7) - (k^8*(-105 + 16*Sqrt[1 + k^2/(4*c - I*k0)^2])*(4*c - I*k0) + 8*k^6*(-105 + 44*Sqrt[1 + k^2/(4*c - I*k0)^2])*(4*c - I*k0)^3 + 144*k^4*(-14 + 9*Sqrt[1 + k^2/(4*c - I*k0)^2])*(4*c - I*k0)^5 + 320*k^2*(-6 + 5*Sqrt[1 + k^2/(4*c - I*k0)^2])*(4*c - I*k0)^7 + 640*(-1 + Sqrt[1 + k^2/(4*c - I*k0)^2])*(4*c - I*k0)^9)/(144*k^7) + (k^8*(-105 + 16*Sqrt[1 + k^2/(5*c - I*k0)^2])*(5*c - I*k0) + 8*k^6*(-105 + 44*Sqrt[1 + k^2/(5*c - I*k0)^2])*(5*c - I*k0)^3 + 144*k^4*(-14 + 9*Sqrt[1 + k^2/(5*c - I*k0)^2])*(5*c - I*k0)^5 + 320*k^2*(-6 + 5*Sqrt[1 + k^2/(5*c - I*k0)^2])*(5*c - I*k0)^7 + 640*(-1 + Sqrt[1 + k^2/(5*c - I*k0)^2])*(5*c - I*k0)^9)/(144*k^7) - (k^8*(-105 + 16*Sqrt[1 + k^2/(6*c - I*k0)^2])*(6*c - I*k0) + 8*k^6*(-105 + 44*Sqrt[1 + k^2/(6*c - I*k0)^2])*(6*c - I*k0)^3 + 144*k^4*(-14 + 9*Sqrt[1 + k^2/(6*c - I*k0)^2])*(6*c - I*k0)^5 + 320*k^2*(-6 + 5*Sqrt[1 + k^2/(6*c - I*k0)^2])*(6*c - I*k0)^7 + 640*(-1 + Sqrt[1 + k^2/(6*c - I*k0)^2])*(6*c - I*k0)^9)/(240*k^7) + (k^8*(-105 + 16*Sqrt[1 + k^2/(7*c - I*k0)^2])*(7*c - I*k0) + 8*k^6*(-105 + 44*Sqrt[1 + k^2/(7*c - I*k0)^2])*(7*c - I*k0)^3 + 144*k^4*(-14 + 9*Sqrt[1 + k^2/(7*c - I*k0)^2])*(7*c - I*k0)^5 + 320*k^2*(-6 + 5*Sqrt[1 + k^2/(7*c - I*k0)^2])*(7*c - I*k0)^7 + 640*(-1 + Sqrt[1 + k^2/(7*c - I*k0)^2])*(7*c - I*k0)^9)/(720*k^7) - (k^8*(-105 + 16*Sqrt[1 + k^2/(8*c - I*k0)^2])*(8*c - I*k0) + 8*k^6*(-105 + 44*Sqrt[1 + k^2/(8*c - I*k0)^2])*(8*c - I*k0)^3 + 144*k^4*(-14 + 9*Sqrt[1 + k^2/(8*c - I*k0)^2])*(8*c - I*k0)^5 + 320*k^2*(-6 + 5*Sqrt[1 + k^2/(8*c - I*k0)^2])*(8*c - I*k0)^7 + 640*(-1 + Sqrt[1 + k^2/(8*c - I*k0)^2])*(8*c - I*k0)^9)/(5040*k^7))/k0^4
SeriesData[k, Infinity, {(1920*c^7)/k0^4, (-93555*c^8)/(2*k0^4) + ((10395*I)*c^7)/k0^3, (480000*c^9)/k0^4 - ((207360*I)*c^8)/k0^3 - (23040*c^7)/k0^2, (-45045*(198*c^10 - (125*I)*c^9*k0 - 27*c^8*k0^2 + (2*I)*c^7*k0^3))/(4*k0^4), 0, (45045*(14391*c^12 - (14457*I)*c^11*k0 - 5940*c^10*k0^2 + (1250*I)*c^9*k0^3 + 135*c^8*k0^4 - (6*I)*c^7*k0^5))/(16*k0^4)}, 5, 11, 1]