数值分析,选主元的Doolittle矩阵分解,一段C++程序