x=0:0.1:26; %0阶贝塞尔函数的零点
y=besselj(0,x);
LD0=[];
for k=1:260
if y(k)*y(k+1)<0
h=interp1(y(k:k+1),x(k:k+1),0);
LD0=[LD0,h];
end
end
%%%%%%%%一阶贝塞尔函数的零点
x=0:0.1:26;
y=besselj(1,x);
LD1=[];
for k=1:260
if y(k)*y(k+1)<0
h=interp1(y(k:k+1),x(k:k+1),0);
LD1=[LD1,h];
end
end