Heilbronn矩阵计算