  We consider two non-linear generalizations of fractal interpolating functions
generated from iterated function systems. The first corresponds to fitting data
using a Kth-order polynomial, while the second relates to the freedom of adding
certain arbitrary functions. An escape-time algorithm that can be used for such
systems to generate fractal images like those associated with Julia or
Mandelbrot sets is also described.
