FractalPark
NewtonHard

Newton Cosh

Newton's method applied to cosh z = 1: infinitely many evenly spaced double roots turn the plane into repeating convergence bands.

Open in Explorer
Newton fractal for cosh z = 1, with repeating horizontal bands and warped boundaries stacked along the imaginary axis
FractalPark’s Newton view of cosh z = 1: repeating convergence bands along the imaginary axis, separated by folded, threadlike borders.

Overview

Newton's usual root-finding step changes character when f(z) = cosh z − 1. Its roots are 2πik for every integer k, and each is a double root. A polynomial gives an orbit only finitely many possible destinations; these roots keep marching along the imaginary axis at equal intervals. The fractal therefore arranges itself into bands running upward and downward, rather than settling into a few basins.

The correction (cosh z − 1) / sinh z simplifies, away from zeros, to tanh(z/2). The denominator sinh z vanishes at every πik: the even multiples are the roots and are removable in the simplified quotient, while the genuine poles and sensitive singularities are the odd multiples (2k+1)πi. Near those odd multiples, a tiny change in the starting point can greatly change the step.

You can see the periodicity along the imaginary axis at a glance — no math background needed. The horizontal bands repeat above and below one another, one for each period of the roots. At their joins, fine vertical threads record where one basin hands an orbit to the next.

The Mathematics

Newton’s cosh iteration

z(n+1) = z(n) - (cosh(z(n)) - 1) / sinh(z(n))

Newton's rule is N(z) = z − f(z) / f′(z). Here f(z) = cosh z − 1 and f′(z) = sinh z. Each step subtracts the correction (cosh z − 1) / sinh z, which is tanh(z/2) away from the zeros of sinh z.

The equation cosh z = 1 has infinitely many solutions: z = 2πik for every integer k. Each is a double root, so nearby Newton errors are approximately halved on each step: convergence is linear, not quadratic. The denominator vanishes at every πik. At even multiples the quotient has a removable value because those points are roots; the genuine poles are the odd multiples (2k+1)πi.

Try z₀ = 0.2i. The correction is tanh(0.1i) = i tan(0.1) ≈ 0.1003i, so z₁ ≈ 0.0997i: the error is visibly about halved. Near an odd multiple such as πi, the correction can become very large. In the implementation, the real component is clamped for numerical stability, and the step is left unchanged when |sinh(z)|² < 1e-10; those are safeguards, not claims about the unsafeguarded map.

The function repeats with period 2πi in the imaginary direction, and the Newton map does too: N(z + 2πi) = N(z) + 2πi. That repeated rule is the source of the horizontal banding.

Visual Characteristics

Start with the horizontal bands along the imaginary axis. A period of roots makes one band, and the convergence pattern returns in the next. The borders do not run straight: they bend and fold, then pull out fine vertical filaments like taut threads.

Within each band, the basin geometry resembles its neighbours, with small distortions. Think of a ribbon folded repeatedly. The warm and cool zones come from the colouring; the sharp breaks gather near the odd multiples of πi, where the Newton correction has genuine poles.

Remix and Examples

Start from the canonical Document

Open the same approved formula state used by this guide, then change the view, coloring, transforms, or animation in the interactive Explorer.

Frequently Asked Questions

What equation is being solved?

It seeks the complex solutions of cosh(z) = 1: z = 2πik for every integer k. These equally spaced roots line the imaginary axis, and each is a double root. Their repeating arrangement creates the horizontal bands.

Why are some regions extremely sensitive?

The correction divides by sinh(z), which vanishes at z = πik for integer k. The even multiples are roots and are removable in the simplified quotient; the genuine poles are the odd multiples (2k+1)πi. Close to an odd multiple, the correction can become huge, sending an orbit across several basin boundaries and pulling fine vertical filaments and unsettled border zones into view.

Formula Record

newtonCosh

This canonical FractalPark implementation is published and runnable.

Deterministic preview of newtonCosh

Canonical source

Language
frm-like/1
Standard library
1
Profile mode
parameter-plane
Profile center
0, 0
Profile zoom
0.25
Profile iterations
96
Profile evidence
mechanical

Parameters

No declared parameters

Canonical source

This verified canonical revision is read-only. Remix creates a separate editable fork.

Loading canonical source…

Source and implementation

Formula ID
a89891b1-8ccb-5d58-9fbb-05944b85ce3c
Canonical name
newtonCosh
Original name
newtonCosh
Historical source
FractalPark
Current implementation
Project-owned

This is a FractalPark project-owned implementation. The linked file is the pinned canonical Definition.

Rights or attribution concern?

Send the Formula ID and supporting evidence. The maintainer can hold or withdraw the implementation while the claim is reviewed. contact@fractalpark.com