矩阵数值计算中龙贝格积分MATLAB程序分析