c语言矩阵分解程序