小编jol*_*ene的帖子

我正在尝试对一个物理问题进行集成，而我编写的代码给了我大约 10 倍的结果。我想知道是否有人可以为我指出正确的方向，无论是我可怕的三重 for 循环还是其他出错的地方。

我正在尝试做这个计算。（如果你有兴趣，它来自这篇关于使用 Hylleraas 坐标计算锂原子基态能量的论文！！）

这是论文的相关部分，我将在下面解释我是如何分解它的。

长话短说，为了获得积分值，无限和 (5) 在 q 的 10 个值之后被截断。公式 T (q) 在 (6) 中给出，并且是三个值的三重嵌套和本身：k_1_2、k_2_3 和 k_3_1。

这是我的 T(q) 代码：

def T_q(j1, j2, j3, j_1_2, j_2_3, j_3_1, alpha, beta, gamma,q):

    '''
    T_q formula for I integral summation 

    '''
    #print("q", q )
    L_1_2 = 1/2 * (j_1_2 +1) #sets adjusted values of j12 etc.  
    #print(L_1_2, "L12")
    L_2_3 = 1/2 * (j_2_3 +1)
    #print(L_2_3, "l23")
    L_3_1 = 1/2 * (j_3_1 +1)
    #print(L_3_1, "L31")

    j_1 …