两个问题 一个是你的文件名字mulStablePoint function中间不要有空格 我下载后中间有个加号 你最好连着写 或者写成mulStablePoint_function
另一个你的程序中
f = [p^2 + beta^2 - k^2*n1^2; p^2 + a0^2 - k^2*(n1^2 - n3^2); p^2 + a^2 - k^2*(n1^2 - n2^2);(a*Dbsli*besselj(m,p*b)-p*Dbsli*Dbslj) * (a0*Dbslk*bessely(m,pc)-p*besselk(m,a0*c)*Dbsly) + (a*Dbsli*bessely(m,p*b)-p*besseli(m,a*b)*Dbslym) * (p*Dbsljm*besselk(m,a0*c)-a0*besselj(m,p*c)*Dbslk)]
这个有个pc 我觉得应该p*c改过来后程序调试成功了
结果为
f =
beta^2 + p^2 - 5384511084535211/128
a0^2 + p^2
a^2 + p^2 - 5719360367604753/1024
(a*bessely(0, p/50000)*(besseli(m - 1, a*b)/2 + besseli(m + 1, a*b)/2) - p*besseli(0, a/50000)*(bessely(m - 1, b*p)/2 - bessely(m + 1, b*p)/2))*(a0*besselj(0, (3*p)/100000)*(besselk(m - 1, a0*c)/2 - besselk(m + 1, a0*c)/2) + p*besselk(0, (3*a0)/100000)*(besselj(m - 1, c*p)/2 - besselj(m + 1, c*p)/2)) + (p*(besseli(m - 1, a*b)/2 + besseli(m + 1, a*b)/2)*(besselj(m - 1, b*p)/2 - besselj(m + 1, b*p)/2) - a*besselj(0, p/50000)*(besseli(m - 1, a*b)/2 + besseli(m + 1, a*b)/2))*(a0*bessely(0, (3*p)/100000)*(besselk(m - 1, a0*c)/2 - besselk(m + 1, a0*c)/2) + p*besselk(0, (3*a0)/100000)*(bessely(m - 1, c*p)/2 - bessely(m + 1, c*p)/2))