可以利用solve函数和bessely(0,x)函数、besselj(0,x)函数求解。

解不唯一，我随便用两个数，R和rd你根据自己需要修改。一直到plot这样就是画图看看解的情况。然后下面那个fzero是用数值方式求解，中括号里面的两个数是你想要求解的区间范围，根据你自己的需求修改吧。
R=10;
rd=1;

fun=@(x) bessely(0,x*R)*besselj(1,x*rd)-besselj(0,x*R)*bessely(1,x*rd)

x=linspace(0.1,10,100);
for k=1:length(x)
    y(k)=fun(x(k));
end
plot(x,y)

fzero(fun,[0.1 0.5])，