用于解决广义特征值问题的包

Question

H_ij 是矩阵元素。运算符 det() 代表矩阵的行列式。在这里，n=15, 16, ..., 1000。我想设计代码来解决关于变量 E 的这个方程。我想知道可以使用哪些python包或Mathematica包。