自适应数值极限（和残差）估计：仅从函数值对极限（带有误差估计）进行数值外推-matlab开发

LIMEST 将在给定点找到一般函数（仅指定用于评估）的极限。 您可能会认为limest 像quad，但要考虑限制。 虽然对我来说，这似乎比做其他事情更经常作为家庭作业，但对于一般情况，尽可能有力地解决是一个有趣的问题。 作为一个例子，我将使用一个中等难度的，它易于分析，但在数值上更难处理。 乐趣 = @(x) (exp(x)-1-x)./x.^2; 无法在 MATLAB 中在 x = 0 处计算此函数，返回 NaN。 虽然泰勒级数展开式显示极限为 1/2，但对接近零的乐趣的蛮力评估会导致数值垃圾，因为有两个取消顺序。 乐趣(0) 答案 = NaN 乐趣(1e-15) 答案 = 110223024625156 乐趣(1e-10) 答案 = 827.403709626582 乐趣(1e-5) 答案 = 0.500000696482408 乐趣(1e-2) 答案 = 0.50