代码如下：
Hmkj = Table[0, {i, 1, N1}];
HH2[m_, k_, j_] :=
G1[[k, j]]*q[m, k, j]*
BesselI[2,
q[m, k, j]*r[j]]*(Hmkj[[j + 1]]*BesselI[1, q[m, k, j + 1]*r[j]] +
BesselK[1, q[m, k, j + 1]*r[j]]) +
G1[[k, j + 1]]*q[m, k, j + 1]*
BesselI[1,
q[m, k, j]*r[j]]*(BesselK[2, q[m, k, j + 1]*r[j]] -
Hmkj[[j + 1]]*BesselI[2, q[m, k, j + 1]*r[j]])
HH3[m_, k_, j_] :=
G1[[k, j]]*q[m, k, j]*
BesselI[2,
q[m, k, j]*r[j]]*(Hmkj[[j + 1]]*BesselI[1, q[m, k, j + 1]*r[j]] +
BesselK[1, q[m, k, j + 1]*r[j]]) +
G1[[k, j + 1]]*q[m, k, j + 1]*BesselI[1, q[m, k, j]*r[j]]*
BesselK[2, q[m, k, j + 1]*r[j]] -
G1[[k, j + 1]]*q[m, k, j + 1]*BesselI[1, q[m, k, j]*r[j]]*
Hmkj[[j + 1]]*BesselI[2, q[m, k, j + 1]*r[j]];
HH2[24, 1, 1]
求助大佬，帮我看看为啥两种形式计算结果不一致，万分感谢！！！！！！！
希望别被系统删帖了！！！！！！！

不完整的代码等于没有。非要我猜一个的话，你把G1之类的数值表的精度弄高点试试。