Question

Maple求解超越方程，未知数f,n0,n1,n2

eqns := {21.20575042*sqrt(-7.831046235*10^(-16) - 1.302684273*10^(-18)/n0)/sin(0.07539822370*f/sqrt(-7.831046235*10^(-16) - 1.302684273*10^(-18)/n0)) = (21.20575042*sqrt(-7.831046235*10^(-16) - 1.302684273*10^(-18)/n0)*tan(0.03769911185*f/sqrt(-7.831046235*10^(-16) - 1.302684273*10^(-18)/n0)) + 39285.04773*sqrt(1/(1 - 0.34/n1))*tan(0.00004883922020*f/sqrt(1/(1 - 0.34/n1))))*(21.20575042*sqrt(-7.831046235*10^(-16) - 1.302684273*10^(-18)/n0)*tan(0.03769911185*f/sqrt(-7.831046235*10^(-16) - 1.302684273*10^(-18)/n0)) - 3.860999998*10^8/((1 - 0.34/n2)*f) - 39285.04773*sqrt(1/(1 - 0.34/n2))*cot(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))) + 62856.07637*sqrt(1/(1 - 0.34/n2))/sin(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))) - 0.00001590936084*(6.177599996*10^8/((1 - 0.34/n2)*f) - 100569.7222*sqrt(1/(1 - 0.34/n2))*cot(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))))*sin(0.00006104902525*f/sqrt(1/(1 - 0.34/n2)))/((6.177599996*10^8/((1 - 0.34/n2)*f) - 100569.7222*sqrt(1/(1 - 0.34/n2))*cot(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))) + 62856.07637*sqrt(1/(1 - 0.34/n2))/sin(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))))*sqrt(1/(1 - 0.34/n2))))/(42.41150084*sqrt(-7.831046235*10^(-16) - 1.302684273*10^(-18)/n0)*tan(0.03769911185*f/sqrt(-7.831046235*10^(-16) - 1.302684273*10^(-18)/n0)) + 39285.04773*sqrt(1/(1 - 0.34/n1))*tan(0.00004883922020*f/sqrt(1/(1 - 0.34/n1))) - 3.860999998*10^8/((1 - 0.34/n2)*f) - 39285.04773*sqrt(1/(1 - 0.34/n2))*cot(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))) + 62856.07637*sqrt(1/(1 - 0.34/n2))/sin(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))) - 0.00001590936084*(6.177599996*10^8/((1 - 0.34/n2)*f) - 100569.7222*sqrt(1/(1 - 0.34/n2))*cot(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))))*sin(0.00006104902525*f/sqrt(1/(1 - 0.34/n2)))/((6.177599996*10^8/((1 - 0.34/n2)*f) - 100569.7222*sqrt(1/(1 - 0.34/n2))*cot(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))) + 62856.07637*sqrt(1/(1 - 0.34/n2))/sin(0.00006104902525*f/sqrt(1/(1 - 0.34/n2))))*sqrt(1/(1 - 0.34/n2)))), 0.00006600709707*f*BesselJ(0, 0.00006600709707*f/sqrt((-0.4317073171*n0 + 1)/(-0.8634146342*n0 + 0.6707317073)))*(-0.4317073171*n0 + 1)/sqrt((-0.4317073171*n0 + 1)/(-0.8634146342*n0 + 0.6707317073)) - BesselJ(1, 0.00006600709707*f/sqrt((-0.4317073171*n0 + 1)/(-0.8634146342*n0 + 0.6707317073)))*(-0.8634146342*n0 + 0.6707317073) = 0, 0.00004240161354*f*BesselJ(0, 0.00004240161354*f/sqrt((-0.34*n1 + 1)/(-0.68*n1 + 0.66)))*(-0.34*n1 + 1)/sqrt((-0.34*n1 + 1)/(-0.68*n1 + 0.66)) - BesselJ(1, 0.00004240161354*f/sqrt((-0.34*n1 + 1)/(-0.68*n1 + 0.66)))*(-0.68*n1 + 0.66) = 0, 0.00005560919483*f*BesselJ(0, 0.00005560919483*f/sqrt((-0.34*n2 + 1)/(-0.68*n2 + 0.66)))*(-0.34*n2 + 1)/sqrt((-0.34*n2 + 1)/(-0.68*n2 + 0.66)) - BesselJ(1, 0.00005560919483*f/sqrt((-0.34*n2 + 1)/(-0.68*n2 + 0.66)))*(-0.68*n2 + 0.66) = 0, 10000 < f and f < 50000}