聚类算法中高斯核局部密度函数怎么用python写出来

我看见有人这样写的

    # Gaussian kernel

    for i in range(4):

        for j in range(i + 1, 5):

            rho[i] = rho[i] + np.exp(-(distance_all[i, j] / dc) * (distance_all[i, j] / dc))

            rho[j] = rho[j] + np.exp(-(distance_all[i, j] / dc) * (distance_all[i, j] / dc))

可是自己找了几个数算了下不对。。这个函数还有其他的写法吗？

def gaussian(n):
    ret = 0
    for j in range(1,n+1):
        ret+=np.exp(-(distance_all[i,j]/dc)**2)
    return ret

单看这个公式应该是这么写吧

math库里面好像有

这个参考下算法，算法里面都会提到的