代码如下：
L = 5;
c = 0.2;
FF = Tan[L*b] - (c*b)/(b^2 - 0.01);
Plot[FF, {b, 200, 210}]
q1 = b /. FindRoot[FF, {b, 0.3}]
q2 = b /. FindRoot[FF, {b, 0.5}]
q3 = b /. FindRoot[FF, {b, 7}];
q4 = b /. FindRoot[FF, {b, 7.5}];
q4 - q3
m = 50;
B = Table[0, {j, 1, m}];
B[[1]] = B[[1]] + b /. FindRoot[FF, {b, 0.3}];
B[[2]] = B[[2]] + b /. FindRoot[FF, {b, 0.5}];
NE = 0.6;
Table[
x1 = B[[2]] + (i - 2)*NE1;
B[[i]] =
B[[i]] + b /.
FindRoot[FF, {b, B[[i - 1]] + NE}, MaxIterations -> Infinity]
, {i, 3, m}];
B[[30]]

在已知L,c时，可以通过调试得到一系类的解，但是如果c是变化的数，如何得到任意c下的B[[x]]呢？

求助求助！！！！！

Tan以Pi为周期且每个周期内取值从负无穷到正无穷，所以每个周期内至少有一个根，指定初值时指定在不同周期里就可以得到不同的解了
fun[L_, c_] := Module[{init, sol},
(*初始值，代表在0.2到40范围内设置初始值，可以换成你需要的范围*)
init = Range[0.2, 40, \[Pi]/L];
sol = FindRoot[Tan[L*b] - (c*b)/(b^2 - 0.01), {b, init}];
DeleteDuplicates[b /. sol]
]
fun[5, 0.2]

RootSearch试过了没：https://tieba.baidu.com/p/5891416543?pid=122314916403&cid=0#122314916403