实验目的:
掌握理查森外推法
实验要求:
1. 给出理查森外推算法
2. 用Matlab实现理查森外推算法
3. 用Matlab实现自适应积分算法
实验内容:
1. 理查森外推算法,数学知识:利用Richardson外推对逐次分半,若记
则有
由Richardson外推方法,可得到左式的误差为
考虑舍入误差,m不能取得太大。
算法描述:
(1)命名函数;(2)如果输入的未知数少于4个,默认精度0.001;(3)描述T表示矩阵坐标;(4)依次赋值计算T表第一列;(5)根据Richardson公式求出T表矩阵的值;(6)若达到精度则运算结束,若未达到则循环计算;(7)输出T表ÿ