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

解不唯一，我隨便用兩個(gè)數(shù)，R和rd你根據(jù)自己需要修改。一直到plot這樣就是畫圖看看解的情況。然后下面那個(gè)fzero是用數(shù)值方式求解，中括號(hào)里面的兩個(gè)數(shù)是你想要求解的區(qū)間范圍，根據(jù)你自己的需求修改吧。
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])，