matlab中的命令lu采用的是列主元分解，即PLU；本程序实现的是Doolittle分解。