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.