  In a previous paper, dealing with "Applications in $\mathbb{R}^1$," the
authors developed a new approach to the computation of the Hausdorff dimension
of the invariant set of an iterated function system or IFS and studied some
applications in one dimension. The key idea, which has been known in varying
degrees of generality for many years, is to associate to the IFS a parametrized
family of positive, linear, Perron-Frobenius operators $L_s$. In our context,
$L_s$ is studied in a space of $C^m$ functions and is not compact.
Nevertheless, it is has a strictly positive $C^m$ eigenfunction $v_s$ with
positive eigenvalue $\lambda_s$ equal to the spectral radius of $L_s$. Under
appropriate assumptions on the IFS, the Hausdorff dimension of the invariant
set of the IFS is the value $s=s_*$ for which $\lambda_s =1$. To compute the
Hausdorff dimension of an IFS associated to complex continued fractions, (which
may arise from an infinite iterated function system), we again approximate the
eigenvalue problem by a collocation method, but now using continuous piecewise
bilinear functions. Using the theory of positive linear operators and explicit
a priori bounds on the partial derivatives of the strictly positive
eigenfunction $v_s$, we are able to give rigorous upper and lower bounds for
the Hausdorff dimension $s_*$, and these bounds converge to $s_*$ as the mesh
size approaches zero. We also demonstrate by numerical computations that
improved estimates can be obtained by the use of higher order piecewise tensor
product polynomial approximations, although the present theory does not
guarantee that these are strict upper and lower bounds. An important feature of
our approach is that it also applies to the much more general problem of
computing approximations to the spectral radius of positive transfer operators,
which arise in many other applications.
