void qr(double *a, double *d, int n) 矩阵的QR分解 假设数组an*n在内存中按行优先次序存放，此函数使用HouseHolder变换将其就地进行QR分解。 d为输出参数，d[0,n)存放QR分解的上三角矩阵对角线元素。 bool householder(double const *qr, double const *d, double *b, int n) 求线性代数方程组的解。 假设矩阵qrn*n为某个矩阵an*n的QR分解，在内