说起编程中的高精度数值,我第一反应就是double类型了.的确,double阶码11位,尾数52位,几乎能应对任何苛刻的要求......然而,当我天真地尝试用double来算泰勒展开式的函数值,离散代数的狰狞性瞬间完全暴露眼前. 我不是数学出身,也并不关心那些奇异的函数形式.我只是想推导一个物理系统的冲击响应函数.那个函数没法求得解析解,我只好用皮卡迭代法求泰勒级数式.经过一番挣扎,级数通项终于确定了,但是函数的绘图却成了隐藏难题. 函数的形式: \[f(x)={\sum_{n=1}^{\inf…