将上面的数组，代入下面的函数中，其中数组中的数代表lb ，a0是一个常数，求ob

ob (lb^2/Sin[ob]^2) - lb^2 (Cot[ob]) - a0 == 0

{0.0130574, 0.0129021, 0.0127472, 0.0125928, 0.0124388, 0.0122851, \
0.0121317, 0.0119786, 0.0118257, 0.011673, 0.0115205, 0.0113681, \
0.0112159, 0.0110637, 0.0109116, 0.0107594, 0.0106073, 0.0104551, \
0.0103028, 0.0101505, 0.009998, 0.00984536, 0.00969254, 0.00953951, \
0.00938626, 0.00923275, 0.00907896, 0.00892487, 0.00877045, \
0.00861569, 0.00846056, 0.00830504, 0.00814911, 0.00799274, \
0.00783593, 0.00767865, 0.00752089, 0.00736262, 0.00720383, \
0.0070445, 0.00688463, 0.00672418, 0.00656316, 0.00640155, \
0.00623933, 0.0060765, 0.00591304, 0.00574895, 0.00558421, \
0.00541883, 0.00525279, 0.00508609, 0.00491873, 0.00475071, \
0.00458201, 0.00441266, 0.00424263, 0.00407195, 0.00390061, \
0.00372862, 0.00355599, 0.00338273, 0.00320885, 0.00303436, \
0.00285928, 0.00268363, 0.00250742, 0.00233067, 0.00215342, \
0.00197567, 0.00179748}

方法很多，这里介绍个基础的Table用法：
lblst = {0.0130574, 0.0129021, 0.0127472, 0.0125928, 0.0124388, 0.0122851, 0.0121317,
0.0119786, 0.0118257, 0.011673, 0.0115205, 0.0113681, 0.0112159, 0.0110637, 0.0109116,
0.0107594, 0.0106073, 0.0104551, 0.0103028, 0.0101505, 0.009998, 0.00984536,
0.00969254, 0.00953951, 0.00938626, 0.00923275, 0.00907896, 0.00892487, 0.00877045,
0.00861569, 0.00846056, 0.00830504, 0.00814911, 0.00799274, 0.00783593, 0.00767865,
0.00752089, 0.00736262, 0.00720383, 0.0070445, 0.00688463, 0.00672418, 0.00656316,
0.00640155, 0.00623933, 0.0060765, 0.00591304, 0.00574895, 0.00558421, 0.00541883,
0.00525279, 0.00508609, 0.00491873, 0.00475071, 0.00458201, 0.00441266, 0.00424263,
0.00407195, 0.00390061, 0.00372862, 0.00355599, 0.00338273, 0.00320885, 0.00303436,
0.00285928, 0.00268363, 0.00250742, 0.00233067, 0.00215342, 0.00197567, 0.00179748};
eq = (ob lb^2)/Sin[ob]^2 - lb^2 Cot[ob];
eqlst = Table[eq, {lb, lblst}]
这方程显然有无穷多组解，如果是要一组解那就限好范围用FindRoot，需要多解那就用RootSearch（library.wolfram.com/infocenter/Demos/4482/）。