超越方程根的问题

wuhaimiao

clc;
clear;
syms x;
H=2.6e-2;
a=130;
b=200;
r11=diff(bessely(0,x*b),x);
r22=diff(besselj(0,x*a),x);
r33=diff(besselj(0,x*b),x);
r44=diff(bessely(0,x*a),x);
f=@(x) ((x*r11+H*bessely(0,x*b))*r22-(x*r33+H*besselj(0,x*b))*r44);
x0=[1 2 3];
sol = arrayfun(@(x)fzero(f,x),x0)

Undefined function or method 'isfinite' for input arguments of type 'sym'.
Error in ==> fzero at 323
    elseif ~isfinite(fx) || ~isreal(fx)
Error in ==> @(x)fzero(f,x)

Error in ==> temperature at 13
sol = arrayfun(@(x)fzero(f,x),x0)

wuhaimiao

已经解决了！
clc;
clear;
syms x;
H=2.6e-2;
a=130;
b=200;
f=@(x)((x.*200*bessely(1,200*x)-H*bessely(0,x*b)).*besselj(1,130*x)-(x.*200*besselj(1,200*x)-H*besselj(0,x*b)).*bessely(1,130*x));
x0=[1 2 3];
sol=arrayfun(@(x)fzero(f,x),x0)