说明：  分别给出利用自适应Simpson积分的数值积分方法（QUAD.m）和自适应Lobatto的数值积分算法（QUADL.m）。
(were given the use of adaptive Simpson integral numerical integration methods (QUAD.m) and adaptive Lobatto should the numerical integration algorithm (QUADL.m).)