Root #2: x=O( 1/ )

Iterations:

> H := x -> normal( (-1/x - 3 )/epsilon );

> x.0 := -3/epsilon;

> for k from 0 to 5 do

> x.(k+1) := H(x.k);

> series( x.(k+1), epsilon, k+2 ); # these show how the iterates converge

> od;