一个通用的理查森外推例程，可以执行多个外推步骤，可以容纳任意领先的误差项，可以处理任何统一的细化方案（例如 h/2、h/3、h/5 等），并且可以报告完整的外推表，包括与已知的良好解决方案进行比较。 能够处理矢量输入。 与 MATLAB 和 Octave 兼容。 包括 Octave 风格的单元测试。
一个人为的用法示例：
假设您的近似序列是 [1, 2, 3, 4]，其中每个连续值都将网格细化了 3 倍。假设近似中的三个主要误差项是 O(h^2)、O(h^4 ) 和 O(h^5) 并且您知道确切值是 4.14。 然后你会调用
Octave:78> [result, normtable] = richardson([1,2,3,4],[2,4,5], 3, 4.14) 结果 = 4.14163223140496 常态 =
3.13999999999999968 NaN NaN NaN
- 1
- 2
前往页