用1stOpt试了一下,不知对否?
Constant a=0.04,b=0.1841;
Function ((k^3+(k-1)/a^2)*(k^3-(k-1)/a^2)*BessI(1,k*a)*bessj(1,k*a))/(k^2*(k^3-(k-1)/a^2)*bessi(1,k*a)*bessj(0,k*a)+k^2*(k^3+(k-1)/a^2)*bessi(0,k*a)*bessj(1,k*a)+2*k^4*(b-1)/a*bessi(1,k*a)*bessj(1,k*a));
(k^3+625*k-625)*(k^3-625*k+625)*bessi(1,1/25*k)*bessj(1,1/25*k)/(k^2*(k^3-625*k+625)*bessi(1,1/25*k)*bessj(0,1/25*k)+k^2*(k^3+625*k-625)*bessi(0,1/25*k)*bessj(1,1/25*k)-8159/200*k^4*bessi(1,1/25*k)*bessj(1,1/25*k));
结果:k =24.4841433328478 |