Question: How to solve this equation?

Hi all

How I use "solve" or "fsolve" for this equation ?

M2 := evalf[4](Matrix(4, 4, {(1, 1) = BesselJ(0, 0.5e-1*sqrt(0.1111111111e-16*omega^2-25.00027778)), (1, 2) = -BesselJ(0, 0.5e-1*sqrt(0.4444444445e-16*omega^2-25)), (1, 3) = -BesselY(0, 0.5e-1*sqrt(0.4444444445e-16*omega^2-25)), (1, 4) = 0, (2, 1) = (0.1111111111e-16*I)*omega*(1-25000000000000/omega^2)*BesselJ(1, 0.5e-1*sqrt(0.1111111111e-16*omega^2-25.00027778))/sqrt(0.1111111111e-16*omega^2-25.00027778), (2, 2) = -(0.4444444444e-16*I)*BesselJ(1, 0.5e-1*sqrt(0.4444444445e-16*omega^2-25))/sqrt(0.4444444445e-16*omega^2-25), (2, 3) = (0.4444444444e-16*I)*BesselY(1, 0.5e-1*sqrt(0.4444444445e-16*omega^2-25))/sqrt(0.4444444445e-16*omega^2-25), (2, 4) = 0, (3, 1) = 0, (3, 2) = BesselJ(0, 0.60e-1*sqrt(0.4444444445e-16*omega^2-25)), (3, 3) = BesselY(0, 0.60e-1*sqrt(0.4444444445e-16*omega^2-25)), (3, 4) = -BesselY(0, 0.60e-1*sqrt(0.1111111111e-16*omega^2-25)), (4, 1) = 0, (4, 2) = (0.4444444444e-16*I)*BesselJ(1, 0.60e-1*sqrt(0.4444444445e-16*omega^2-25))/sqrt(0.4444444445e-16*omega^2-25), (4, 3) = (0.4444444444e-16*I)*BesselY(1, 0.60e-1*sqrt(0.4444444445e-16*omega^2-25))/sqrt(0.4444444445e-16*omega^2-25), (4, 4) = -(0.1111111111e-16*I)*omega*BesselY(1, 0.60e-1*sqrt(0.1111111111e-16*omega^2-25))/sqrt(0.1111111111e-16*omega^2-25)})):

with(LinearAlgebra):
DETM2 := Determinant(M2):
solve(DETM2 = 0, omega);

Error, (in solve) cannot solve for an unknown function with other operations in its arguments

Is this Error because of combination of bessel functions? if I use asymptatic forms, does it work?

Thanks