由C代码实现的数值计算方法——LU分解