With a Pauli matrix, we change the linear respond BdG equations into a matrix.Substitute the numerical result of the ground state wave function;

• 再引入泡利矩阵，把相应的线性响应方程转化为矩阵形式。