cosh 的牛顿迭代
牛顿法的规则是 N(z) = z − f(z) / f′(z)。这里取 f(z) = cosh z − 1、f′(z) = sinh z;每一步先算出修正量 (cosh z − 1) / sinh z,再从 z 里减去它。
方程 cosh z = 1 有无穷多解:z = 2πik,k 为任意整数。它们全是单根:cosh z 在根旁线性增长,牛顿法便在附近二次收敛。远离实轴后,cosh 和 sinh 的指数增长会让修正项忽大忽小;图中的敏感边界便从这里冒出来。
试试 z₀ = 0.2i。由于 i·sinh(0.2) ≈ 0.201i,cosh(0.2) ≈ 1.020,第一次修正很小,轨道很快收敛到 z = 0 处的根。再取 z₀ = 0.5 + 3.0i,它离每个根都远些。cosh 和 sinh 的大指数项会卷进修正,轨道可能跨过几个盆地边界才落定;也可能在根之间混沌地绕许多圈,最后才被某一个捕获。
函数在虚方向以 2πi 为周期,牛顿映射也继承相同的周期性:N(z + 2πi) = N(z) + 2πi。图中的水平带状结构因此反复出现。

