In this paper, we visualise and analyse the dynamics of fractals (Julia and Mandelbrot sets) for complex polynomials of the form $$T(z) = z^{n} + mz + r$$
T
(
z
)
=
z
n
+
m
z
+
r
, where $$n \ge 2$$
n
≥
2
and $$m, r \in \mathbb {C}$$
m
,
r
∈
C
, by adopting the viscosity approximation type iteration process which is most widely used iterative method for finding fixed points of non-linear operators. We establish a convergence condition in the form of escape criterion which allows to adapt the escape-time algorithm to the considered iteration scheme. We also present some graphical examples of the Mandelbrot and Julia fractals showing the dependency of Julia and Mandelbrot sets on complex polynomials, contraction mappings, and iteration parameters. Moreover, we propose two numerical measures that allow the study of the dependency of the set shape change on the values of the iteration parameters. Using these two measures, we show that the dependency for the considered iteration method is non-linear.