贡献者: addis
对于可微的二元函数
计算过程中,由于出现了大数加小数,所以实际上 x(i) + h
增加的值可能比 h
有较大不同,可以用 (x(i) + h) - x(i)
来计算实际的增量,这样会使结果更精确。
注意该代码的有效数字估计只包括计算减法的截断误差,不包括用差分代替微分的误差。
另外给出一个变精度计算版本,可以使截断误差精确到最后一位有效数字。显然函数 f
需要支持变精度 sym
类型的变量。
友情链接: 超理论坛 | ©小时科技 保留一切权利