Plot empty lattice modes in lat2d_realfreqsvd.py

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

((35*k^8 + 16*k^6*(35 - 12*Sqrt[1 + k^2/(c - I*k0)^2])*(c - I*k0)^2 + 16*k^4*(105 - 64*Sqrt[1 + k^2/(c - I*k0)^2])*(c - I*k0)^4 + 64*k^2*(28 - 23*Sqrt[1 + k^2/(c - I*k0)^2])*(c - I*k0)^6 - 640*(-1 + Sqrt[1 + k^2/(c - I*k0)^2])*(c - I*k0)^8)/(6720*k^6) - (35*k^8 + 16*k^6*(35 - 12*Sqrt[1 + k^2/(2*c - I*k0)^2])*(2*c - I*k0)^2 + 16*k^4*(105 - 64*Sqrt[1 + k^2/(2*c - I*k0)^2])*(2*c - I*k0)^4 + 64*k^2*(28 - 23*Sqrt[1 + k^2/(2*c - I*k0)^2])*(2*c - I*k0)^6 - 640*(-1 + Sqrt[1 + k^2/(2*c - I*k0)^2])*(2*c - I*k0)^8)/(960*k^6) + (35*k^8 + 16*k^6*(35 - 12*Sqrt[1 + k^2/(3*c - I*k0)^2])*(3*c - I*k0)^2 + 16*k^4*(105 - 64*Sqrt[1 + k^2/(3*c - I*k0)^2])*(3*c - I*k0)^4 + 64*k^2*(28 - 23*Sqrt[1 + k^2/(3*c - I*k0)^2])*(3*c - I*k0)^6 - 640*(-1 + Sqrt[1 + k^2/(3*c - I*k0)^2])*(3*c - I*k0)^8)/(320*k^6) - (35*k^8 + 16*k^6*(35 - 12*Sqrt[1 + k^2/(4*c - I*k0)^2])*(4*c - I*k0)^2 + 16*k^4*(105 - 64*Sqrt[1 + k^2/(4*c - I*k0)^2])*(4*c - I*k0)^4 + 64*k^2*(28 - 23*Sqrt[1 + k^2/(4*c - I*k0)^2])*(4*c - I*k0)^6 - 640*(-1 + Sqrt[1 + k^2/(4*c - I*k0)^2])*(4*c - I*k0)^8)/(192*k^6) + (35*k^8 + 16*k^6*(35 - 12*Sqrt[1 + k^2/(5*c - I*k0)^2])*(5*c - I*k0)^2 + 16*k^4*(105 - 64*Sqrt[1 + k^2/(5*c - I*k0)^2])*(5*c - I*k0)^4 + 64*k^2*(28 - 23*Sqrt[1 + k^2/(5*c - I*k0)^2])*(5*c - I*k0)^6 - 640*(-1 + Sqrt[1 + k^2/(5*c - I*k0)^2])*(5*c - I*k0)^8)/(192*k^6) - (35*k^8 + 16*k^6*(35 - 12*Sqrt[1 + k^2/(6*c - I*k0)^2])*(6*c - I*k0)^2 + 16*k^4*(105 - 64*Sqrt[1 + k^2/(6*c - I*k0)^2])*(6*c - I*k0)^4 + 64*k^2*(28 - 23*Sqrt[1 + k^2/(6*c - I*k0)^2])*(6*c - I*k0)^6 - 640*(-1 + Sqrt[1 + k^2/(6*c - I*k0)^2])*(6*c - I*k0)^8)/(320*k^6) + (35*k^8 + 16*k^6*(35 - 12*Sqrt[1 + k^2/(7*c - I*k0)^2])*(7*c - I*k0)^2 + 16*k^4*(105 - 64*Sqrt[1 + k^2/(7*c - I*k0)^2])*(7*c - I*k0)^4 + 64*k^2*(28 - 23*Sqrt[1 + k^2/(7*c - I*k0)^2])*(7*c - I*k0)^6 - 640*(-1 + Sqrt[1 + k^2/(7*c - I*k0)^2])*(7*c - I*k0)^8)/(960*k^6) - (35*k^8 + 16*k^6*(35 - 12*Sqrt[1 + k^2/(8*c - I*k0)^2])*(8*c - I*k0)^2 + 16*k^4*(105 - 64*Sqrt[1 + k^2/(8*c - I*k0)^2])*(8*c - I*k0)^4 + 64*k^2*(28 - 23*Sqrt[1 + k^2/(8*c - I*k0)^2])*(8*c - I*k0)^6 - 640*(-1 + Sqrt[1 + k^2/(8*c - I*k0)^2])*(8*c - I*k0)^8)/(6720*k^6))/k0^4
SeriesData[k, Infinity, {(945*c^7)/k0^4, (-17280*c^8)/k0^4 + ((3840*I)*c^7)/k0^3, (3465*(125*c^9 - (54*I)*c^8*k0 - 6*c^7*k0^2))/(4*k0^4), 0, (-9009*(4819*c^11 - (3960*I)*c^10*k0 - 1250*c^9*k0^2 + (180*I)*c^8*k0^3 + 10*c^7*k0^4))/(16*k0^4)}, 5, 11, 1]