这是上一个问题的后续解答：Square Root of Positive-Semidefinite Matrix, worried about bad eigenvalue
当我第一次运行矩阵平方根命令时，它很快就完成了。但是随后的每一次，它的速度都异常缓慢。很奇怪。我只能以为第一次出了点问题，但结果却很快但不正确。
无论如何，我认为我应该利用问题的结构。它不仅仅是任何对称的正半定矩阵；这是特定的：所有对角线条目都是相同的值（称为a），而所有非对角线条目都是相同的值（称为b）。此外，a>0，b<0和a=(n-1)|b|。
我假设平方根矩阵还将具有所有具有相同值的对角线条目（称为x）和所有具有相同值的非对角线条目（称为y）。那么我要求解的方程组是：
x^2+(n-1)y^2=a
2xy+(n-2)y^2=b
已知n,a,b。
numpy可以用数字方式求解这样的系统吗？
谢谢