Open access peer-reviewed chapter

# How Are Fractal Interpolation Functions Related to Several Contractions?

Written By

SongIl Ri and Vasileios Drakopoulos

Submitted: December 23rd, 2019 Reviewed: April 27th, 2020 Published: August 21st, 2020

DOI: 10.5772/intechopen.92662

From the Edited Volume

## Mathematical Theorems

Edited by Lyudmila Alexeyeva

Chapter metrics overview

View Full Metrics

## Abstract

This chapter provides an overview of several types of fractal interpolation functions that are often studied by many researchers and includes some of the latest research made by the authors. Furthermore, it focuses on the connections between fractal interpolation functions resulting from Banach contractions as well as those resulting from Rakotch contractions. Our aim is to give theoretical and practical significance for the generation of fractal (graph of) functions in two and three dimensions for interpolation purposes that are not necessarily associated with Banach contractions.

### Keywords

• attractor
• contraction
• fixed point
• iterated function system
• fractal interpolation

## 1. Introduction

Interpolation is a method of constructing new data points within the range of a discrete set of known data points or the process of estimating the value of a function at a point from its values at nearby points. Although a large number of interpolation schemes are available in the mathematical field of numerical analysis, the majority of these conventional interpolation methods produce interpolants, i.e., functions used to generate interpolation, that are differentiable a number of times except possibly at a finite set of points. Taking into account that the smoothness of a function is a property measured by the number of continuous derivatives it has over some domain, the aforementioned interpolants are considered smooth.

On the other hand, many real-world and experimental signals are intricate and rarely show a sensation of smoothness in their traces. Consequently, to model these signals, we require interpolants that are nondifferentiable in dense sets of points in the domain. To address this issue, interpolation by fractal (graph of) functions is introduced in [1, 2], which is based on the theory of iterated function system. A fractal interpolation function can be considered as a continuous function whose graph is the attractor, a fractal set, of an appropriately chosen iterated function system. If this graph has a Hausdorff-Besicovitch dimension between 1 and 2, the resulting attractor is called fractal interpolation curved line or fractal interpolation curve. If this graph has a Hausdorff-Besicovitch dimension between 2 and 3, the resulting attractor is called fractal interpolation surface. Various types of fractal interpolation functions have been constructed, and some significant properties of them, including calculus, dimension, smoothness, stability, perturbation error, etc., have been widely studied [3, 4, 5].

Fractal interpolation is an advanced technique for analysis and synthesis of scientific and engineering data, whereas the approximation of natural curves and surfaces in these areas has emerged as an important research field. Fractal functions are currently being given considerable attention due to their applications in areas such as Metallurgy, Earth Sciences, Surface Physics, Chemistry and Medical Sciences. In the development of fractal interpolation theory, many researchers have generalised the notion in different ways [6, 7, 8, 9]. Two key issues should be addressed in constructing fractal interpolation functions. They regard to ensuring continuity and the existence of the contractivity, or vertical scaling, factors; see [10, 11]. In [12], nonlinear fractal interpolation surfaces resulting from Rakotch or Geraghty contractions together with some continuity conditions were introduced as well as explicit illustrative examples were given.

The concept of iterated function system was originally introduced as a generalisation of the well-known Banach contraction principle. Since it has become a powerful tool for constructing and analysing fractal interpolation functions, one can use the well-known fixed point results obtained in the fixed point theory in order to construct them in a more general sense. A comparison of various definitions of contractive mappings as well as fixed point theorems that can be used to construct iterated function systems can be found in [13, 14, 15]. In [14], the authors proposed some iterated function systems by using various fixed point theorems, but unfortunately, one does not know whether fractal interpolation functions correspond to those may exist or not. As far as we know, the first significant generalisation of Banach’s principle was obtained by Rakotch [16] in 1962. Recently, a method to generate nonlinear fractal interpolation functions by using the Rakotch or Geraghty fixed point theorem instead of Banach fixed point theorem was presented in [12, 17, 18].

The aim of our article is to provide the connections between several fractal interpolation functions and the contractions used to generate them; it is organised as follows. In Section 2, we recall the results obtained in construction of fractal interpolation curved lines and fractal interpolation surfaces by using Rakotch contractions (or Geraghty contractions) instead of Banach contractions. In Section 3, we only present the connection between fractal interpolation functions by using the Banach contractions and fractal interpolation functions by using the Rakotch contractions because in the case of Geraghty contractions, the existence of fractal interpolation curved lines and fractal interpolation surfaces is similar to the case of Rakotch contractions.

## 2. Preliminaries

Let X ρ and Y σ be metric spaces. A mapping T : X Y is called a Hölder mapping of exponent or order a, if

σ T x T y c ρ x y a

for x , y X , a 0 and for some constant c . Note that, if a > 1 , the functions are constants. Obviously, c 0 . The mapping T is called a Lipschitz mapping, if a may be taken to be equal to 1 . If c = 1 , T is said to be nonexpansive. A Lipschitz function is a contraction with contractivity factor c , if c < 1 . We call T contractive, if for all x , y X and x y , we have σ T x T y < ρ x y . Note that ‘contraction contractive nonexpansive Lipschitz’.

An iterated function system, or IFS for short, is a collection of a complete metric space X ρ together with a finite set of continuous mappings, f n : X X , n = 1 , 2 , , N . It is often convenient to write an IFS formally as X f 1 f 2 f N or, somewhat more briefly, as X f 1 N . The associated map of subsets W : H X H X is given by:

W E = n = 1 N f n E for all E H X ,

where H X is the metric space of all nonempty, compact subsets of X with respect to some metric, e.g., the Hausdorff metric. The map W is called the Hutchinson operator or the collage map to alert us to the fact that W E is formed as a union or ‘collage’ of sets.

If w n are contractions with corresponding contractivity factors s n for n = 1 , 2 , , N , the IFS is termed hyperbolic and the map W itself is then a contraction with contractivity factor s = max s 1 s 2 s N ([2], Theorem 7.1, p. 81). In what follows, we abbreviate by f k the k -fold composition f f f .

Definition 2.1.Let X be a set. A self-map on X or a transformation is a mapping from X to itself.

1. A self-map f on a metric space X ρ is called a φ -contraction, if there exists a function φ : 0 + 0 + with ϕ 0 = 0 and ϕ t < t for all t > 0 such that for all x , y X , ρ f x f y φ ρ x y .

2. We say that f is a Rakotch contraction, if f is a φ -contraction such that for any t > 0 , α t φ t t < 1 and the function 0 + t φ t t is nonincreasing.

3. If f is a φ -contraction for some function φ : 0 + 0 + such that for any t > 0 , α t φ t t < 1 and the function 0 + t φ t t is nonincreasing (or nondecreasing, or continuous), then we call such a function a Geraghty contraction.

From [14], we have the following.

Theorem 2.1.Let X be a complete metric space and X f 1 N be an IFS consisting of Rakotch or Geraghty contractions. Then there is a unique nonempty compact set K H X such that

K = n = 1 N f n K .

### 2.1 Fractal interpolation in R

Let N be a positive integer greater than 1 and I = x 0 x N R . Let a set of interpolation points x i y i I × R : i = 0 1 N be given, where x 0 < x 1 < < x N and y 0 , y 1 , , y N R . Set I n = x n 1 x n I and define, for all n = 1 , 2 , , N , contractive homeomorphisms L n : I I n by

L n x a n x + b n ,

where the real numbers a n , b n are chosen to ensure that L n I = I n .

Let φ : 0 + 0 + be a nondecreasing continuous function such that for any t > 0 , α t φ t t < 1 and the function 0 + t φ t t is nonincreasing. Let d n : I R be a continuously differentiable function such that

max x I d n x 1 .

Now, consider an IFS of the form I × R w n n = 1 2 N in which the maps are nonlinear transformations of the special structure

w n x y = L n x F n x y = a n x + b n c n x + d n x s n y + e n ,

where the transformations are constrained by the data according to

w n x 0 y 0 = x n 1 y n 1 , w n x N y N = x n y n

for n = 1 , 2 , , N , and s n are some Rakotch or Geraghty contractions.

Let us denote by C D the linear space of all real-valued continuous functions defined on D , i.e., C D = f : D R   |   f continuous . Let C I C I denote the set of continuous functions f : I R such that f x 0 = y 0 and f x N = y N , that is,

C I f C I : f x 0 = y 0 f x N = y N .

Let C I C I C I be the set of continuous functions that pass through the given data points x i y i I × R : i = 0 1 N , that is,

C I f C I : f x i = y i i = 0 1 N .

Define a metric d C I on the space C I by

d C I g h max x x 0 x N g x h x

for all g , h C I . Define a mapping T : C I C I for all f C I by

Tf x F n L n 1 x f L n 1 x = c n L n 1 x + d n L n 1 x s n f L n 1 x + e n

for x x n 1 x n and n = 1 , 2 , , N . From [17], we have the following.

Theorem 2.2.Let I × R w n n = 1 2 N denote the IFS defined above. Let each s n be a bounded Rakotch or Geraghty contraction. Then,

1. there is a unique continuous function f : I R which is a fixed point of T ;

2. f x i = y i for all i = 0 , 1 , , N ;

3. if G I × R is the graph of f , then

G = n = 1 N w n G .

An extremely explicit simple example is the following; cf. [12].

Example 1.Let φ t t 1 + t for t 0 + . Let a set of data x i y i : i = 0 1 N be given, where 0 = x 0 < x 1 < < x N = 1 and y i 0 1 for all i = 0 , 1 , , N . Let for all n = 1 , 2 , , N , d n x x n . Let for y 0 + and n = 1 , 2 , , N , s n y y 1 + ny . That is, each s n is a Rakotch contraction (with the same function φ ) that is not a Banach contraction on 0 + . Let for all n = 1 , 2 , , N ,

w n x y a n x + b n c n x + d n x s n y + e n ,

where

a n = x n x n 1 , b n = x n 1 , c n = y n y n 1 , e n = y n 1 .

Then, there exists a continuous function f : 0 1 R that interpolates the given points x i y i : i = 0 1 N . Moreover, the graph G of f is invariant with respect to 0 1 × R w 1 w 2 w N , i.e.,

G = n = 1 N w n G .

### 2.2 Fractal interpolation in R 2

Let M , N be two positive integers greater than 1 . Let us represent the given set of interpolation points as x i y j z i , j K : i = 0 1 M j = 0 1 N , where x 0 < x 1 < < x M , y 0 < y 1 < < y N and z i , j a b for all i = 0 , 1 , , M and j = 0 , 1 , , N . Set I = x 0 x M R and J = y 0 y N R . Throughout this section, we will work in the complete metric space K = D × R , where D = I × J , with respect to the Euclidean, or to some other equivalent, metric.

Set I m = x m 1 x m , J n = y n 1 y n , D m , n = I m × J n and let u m : I I m , v n : J J n , L m , n : D D m , n be defined for m = 1 , 2 , , M and n = 1 , 2 , , N , by

L m , n x y = u m x v n y = a m x + b m c n y + d n .

Thus, for m = 1 , 2 , , M and n = 1 , 2 , , N ,

a m = x m x m 1 x M x 0 , b m = x m 1 x m x m 1 x M x 0 x 0 , c n = y n y n 1 y N y 0 , d n = y n 1 y n y n 1 y N y 0 y 0 .

Furthermore, for m = 1 , 2 , , M and n = 1 , 2 , , N , let mappings F m , n : K R be continuous with respect to each variable. We consider an IFS of the form K w m , n m = 1 2 M n = 1 2 N in which maps w m , n : D × R D m , n × R are transformations of the special structure

w m , n x y z L m , n x y F m , n x y z ,

where the transformations are constrained by the data according to

w m , n x 0 y 0 z 0 , 0 = x m 1 y n 1 z m 1 , n 1 , w m , n x 0 y N z 0 , N = x m 1 y n z m 1 , n ,
w m , n x M y 0 z M , 0 = x m y n 1 z m , n 1 , w m , n x M y N z M , N = x m y n z m , n

for m = 1 , 2 , , M and n = 1 , 2 , , N .

Let B D denote the set of bounded functions f : D R and

B D = { f B D : f x 0 y 0 = z 0 , 0 , f x 0 y N = z 0 , N , f x M y 0 = z M , 0 , f x M y N = z M , N } .

Let B D B D be the set of bounded functions that pass through the given interpolation points x i y j z i , j K = D × a b : i = 0 1 M j = 0 1 N , that is,

B D = f B D : f x i y j = z i , j i = 0 1 M j = 0 1 N .

Define an operator T : B D B D for all f B D by

Tf x y = F m , n u m 1 x v n 1 y f u m 1 x v n 1 y

for x y D m , n , m = 1 , 2 , , M and n = 1 , 2 , , N . In [18], we see the following.

Theorem 2.3.Let D × R w m , n m = 1 2 M n = 1 2 N denote the IFS defined above. Assume that the maps F m , n are Rakotch or Geraghty contractions with respect to the third variable, and uniformly Lipschitz with respect to the first and second variable. Then,

1. there is a unique bounded function f : D R which is a fixed point of T ;

2. f x i y j = z i , j for i = 0 , 1 , , M and j = 0 , 1 , , N ;

3. if G D × R is the graph of f , then

G = m = 1 M n = 1 N w m , n G .

Let for all i = 0 , 1 , , M and j = 0 , 1 , , N , z 0 , j = z i , 0 = z M , j = z i , N and define

F m , n x y z = e m , n x + f m , n y + g m , n xy + s m , n z + h m , n ,

where s m , n are Rakotch or Geraghty contractions. Let

C D = { f C D : f x 0 y 0 = z 0 , 0 , f x 0 y N = z 0 , N , f x M y 0 = z M , 0 , f x M y N = z M , N }

and

C D = f C D : f x i y j = z i , j i = 0 1 M j = 0 1 N .

Let C 0 D C D be the set of continuous functions f : D R such that

f x 0 1 λ y 0 + λ y N = z , , f x M 1 λ y 0 + λ y N = z , , f 1 λ x 0 + λ x M y 0 = z , , f 1 λ x 0 + λ x M y N = z ,

for all λ 0 1 , where for all i = 0 , 1 , , M and j = 0 , 1 , , N ,

z z 0 , j = z i , 0 = z M , j = z i , N .

Let C 0 D f C 0 D : f x i y j = z i , j i = 0 1 M j = 0 1 N C D . For f C 0 D , we define T : C 0 D B D by

Tf x y = F m , n u m 1 x v n 1 y f u m 1 x v n 1 y = e m , n u m 1 x + f m , n v n 1 y + g m , n u m 1 x v n 1 y + s m , n f u m 1 x v n 1 y + h m , n

for x y D m , n , m = 1 , 2 , , M and n = 1 , 2 , , N .

Corollary 2.1(see [18]) Let D × R w m , n m = 1 2 M n = 1 2 N denote the IFS defined above. Then,

1. there is a unique continuous function f : D R which is a fixed point of T ;

2. f x i y j = z i , j for all i = 0 , 1 , , M and j = 0 , 1 , , N ;

3. if G D × R is the graph of f , then

G = m = 1 M n = 1 N w m , n G .

The most simple example is the following; cf. [12].

Example 2.Let φ t t 1 + t for t 0 + . Let a set of data x i y j z i , j : i = 0,1,2 j = 0,1,2 be given, where 0 = x 0 < x 1 < x 2 = 1 , 0 = y 0 < y 1 < y 2 = 1 and z i , j 0 1 for all i = 0,1,2 ; j = 0,1,2 . Let for all i = 0,1,2 and j = 0,1,2 ,

z 0 , j = z i , 0 = z 2 , j = z i , 2 = 0 .

Let for z 0 + ,

s 1 , 1 z z 1 + z , s 1 , 2 z z 1 + 2 z ,
s 2 , 1 z z 1 + 3 z , s 2 , 2 z z 1 + 4 z .

Then, s 1 , 1 , s 1 , 2 , s 2 , 1 , s 2 , 2 are Rakotch contractions (with the same function φ ) that are not Banach contractions on 0 + . So, there exists a continuous function f : 0 1 × 0 1 R that interpolates the given data x i y j z i , j : i = 0,1,2 j = 0,1,2 .

Let d m , n : D R be a function such that max x y D d m , n x y 1 ,

d m , n x 0 y = d m , n x M y = d m , n x y 0 = d m , n x y N = 0

and for some L 1 , L 2 > 0 ,

d m , n x y d m , n x y L 1 x x + L 2 y y .

Let

F m , n x y z = e m , n x + f m , n y + g m , n xy + d m , n x y s m , n z + h m , n ,

where s m , n is a Rakotch or Geraghty contraction. For f C D , we define T : C D B D by

Tf x y = F m , n u m 1 x v n 1 y f u m 1 x v n 1 y = e m , n u m 1 x + f m , n v n 1 y + g m , n u m 1 x v n 1 y + d m , n u m 1 x v n 1 y s m , n f u m 1 x v n 1 y + h m , n

for x y D m , n , m = 1 , 2 , , M and n = 1 , 2 , , N .

For the next, see [12] for details.

Corollary 2.2.Let D × R w m , n m = 1 2 M n = 1 2 N denote the IFS defined above. If each s m , n be a bounded function, then.

1. there is a unique continuous function f : D R which is a fixed point of T ;

2. f x i y j = z i , j for all i = 0 , 1 , , M and j = 0 , 1 , , N ;

3. if G D × R is a graph of f , then

G = m = 1 M n = 1 N w m , n G .

An especially simple example is the following; see [12].

Example 3.Let φ t t 1 + t for t 0 + . Let a set of data x i y j z i , j : i = 0,1,2 j = 0,1,2 be given, where 0 = x 0 < x 1 < x 2 = 1 , 0 = y 0 < y 1 < y 2 = 1 and z i , j 0 1 for all i = 0,1,2 ; j = 0,1,2 . Here, a set of data points is not necessarily the case that z 0 , j = z i , 0 = z 2 , j = z i , 2 for all i = 0,1,2 ; j = 0,1,2 . Let for all i = 1 , 2 ; j = 1 , 2 and x y 0 1 × 0 1 ,

d m , n x y 2 2 m + n x m 1 x m y n 1 y n .

Let for z 0 + ,

s 1 , 1 z 1 1 + z , s 1 , 2 z z 1 + z ,
s 2 , 1 z z 1 + 2 z , s 2 , 2 z z 1 + 3 z .

Then, s 1 , 1 , s 1 , 2 , s 2 , 1 , s 2 , 2 are Rakotch contractions (with the same function φ ) that are not Banach contractions on 0 + . So, there exists a continuous function f : 0 1 × 0 1 R that interpolates the given data x i y j z i , j : i = 0,1,2 j = 0,1,2 .

## 3. Interconnections between FIFs and contractions

In this section, we only present the interconnections between FIFs resulting from Banach contractions and FIFs resulting from Rakotch contractions because in the case of Geraghty contractions, the existence of FICs and FISs is derived similarly to the case of Rakotch contractions.

### Connection 1

1. Each Banach contraction is a Rakotch contraction, since a self-map is a Banach contraction if and only if it is a φ -contraction for a function φ t = αt , for some 0 α < 1 . There exist examples of Rakotch contraction maps that are not Banach contraction maps on X R with respect to the Euclidean metric (see [13]).

2. The Rakotch’s functional condition for convergence of a contractive iteration in a complete metric space can be replaced by an equivalent (or another) functional condition; for instance, a map is a Rakotch contraction if and only if it is a φ -contraction for some nondecreasing function φ : 0 + 0 + such that additionally φ t < t for t > 0 and the map t φ t t is nonincreasing (see [19]).

### Connection 2

1. C I d C I , C I d C I and C I d C I are complete metric spaces, where

d C I f g max x I f x g x

for all f , g C I (see [2]).

1. B D d B D , B D d B D and B D d B D are complete metric spaces, where

d B D f g sup x y D f x y g x y

for all f , g B D [10].

1. C 0 D , C 0 D , C D , C D and C D are closed subspaces of B D with C 0 D C 0 D C D C D B D and C 0 D C D C D C D B D , and so they are complete metric spaces.

### Connection 3

Let d n : I R be a continuously differentiable function such that

max x I d n x 1 .

Then, by the Differential Mean Value Theorem and the extreme value theorem, we can see that for some L d n > 0 ,

d n x d n x L d n x x ,

where x , x I . Hence, d n is Lipschitz continuous function defined on I satisfying max x I d n x 1 , but the converse is not true in general.

### Connection 4

1. The function d n x s n y is a generalisation of the bivariable function d n x y with vertical scaling factors as (continuous) ‘contraction functions’. In fact, in the case when 0 < max x I d n x < 1 (see [20], p. 3), obviously,

d n x y = d n x max x I d n x max x I d n x y .

Let s n y = max x I d n x y and d n x = d n x max x I d n x . Then d n x y = d n x s n y , max x I d n x = 1 and s n is a Banach (or Rakotch) contraction.

1. The functional condition max x I d n x 1 is essential in order to show the difference between Banach contractibility of F n y and Rakotch contractibility of F n y ; compare with [20]. In fact, since φ t < t for any t > 0 ,

F n x y F n x y = d n x s n y s n y max x I d n x s n y s n y max x I d n x φ y y max x I d n x y y ,

where x y , x y R 2 . Hence, if max x I d n x < 1 , as can be seen, notwithstanding each s n is a Rakotch contraction that is not a Banach contraction, each F n is Banach contraction with respect to the second variable because

F n x y F n x y max x I d n x y y .

On the other hand, if max x I d n x = 1 , then we conclude that each F n is Rakotch contraction with respect to the second variable whenever each s n is a Rakotch contraction because

F n x y F n x y max x I d n x φ y y .

1. In Theorem 2.2, for all x y , x y I × R ,

F n x y F n x y = d n x s n y s n y s n y s n y φ y y .

That is, each w n x y is chosen so that function F n x y is Rakotch contraction with respect to the second variable.

1. Even though s n : R R are Rakotch contractions, w n : I × R are not in general Rakotch contractions on the metric space I × R d 0 , and thus, the IFSs defined above are not IFSs of [14] (cf. second and third line in p. 215 of [2]).

### Connection 5

In the case where the vertical scaling factors are constants, in [1], the existence of affine FIFs by using the Banach fixed point theorem was investigated, whereas in [20], a generalisation of affine FIFs by using vertical scaling factors as (continuous) ‘contraction functions’ and Banach’s fixed point theorem was introduced. Theorem 2.2 gives the existence of fractal interpolation curves by using the Rakotch fixed point theorem and vertical scaling factors as (continuous) ‘contraction functions’.

### Connection 6

The boundedness of s n is the essential condition to establish a unique invariant set of an iterated function system. In the fractal interpolation curve with vertical scaling factors as ‘contraction function’, 0 < max x I d n x < 1 (see [20]). Let M max x I c n x + f n and h M 1 max x I d n x . Then for all y h h ,

F n x y = c n x + d n x y + f n M + max x I d n x y M + max x I d n x h h .

So, for all x y I × h h , we can see that F n x y h h . That is, an IFS of the form I × h h w 1 N has been constructed (cf. [21], p. 1897). Thus D s n = h h and s n y max x I d n x y is bounded in D s n . Hence the boundedness of s n in D s n is the essential condition to establish a unique invariant set of an IFS (cf. [21], p. 1897).

### Connection 7

In view of a φ -contraction, the connections between the coefficients of y variable are obtained as follows:

1. In the affine FIF (cf. [1], p. 308, Example 1), for all t > 0 ,

φ t max n = 1 , 2 , , N d n t ,

where d n < 1 for all i = 1 , 2 , , N .

1. In the FIF with vertical scaling factors as (continuous) ‘contraction functions’ (cf. [20], p. 3), for all t > 0 ,

φ t max i = 1 , 2 , , N max x I d n x t ,

where d n x is Lipschitz function defined on I satisfying sup x I d n x < 1 for all n = 1 , 2 , , N .

### Connection 8

We refer to f of Theorem 2.2 as a nonlinear FIF. The reason is that the functions F n take the form

F n x y = c n x + d n x s n y + e n ,

where max x I d n x 1 and each s n is Rakotch contraction. That is, each F n , in general, is nonlinear with respect to the second variable (cf. [17]). In fact, in [2] or [20], since 0 < d n x d n < 1 or 0 < max x I d n x < 1 and

d n x y = d n x max x I d n x max x I d n x y ,

we can see that

F n x y = c n x + d n x y + e n = c n x + d n x s n y + e n ,

where d n x d n x max x I d n x and s n y max x I d n x y , and thus, each s n is a special Banach contraction and linear with respect to the second variable. Obviously, we can say that nonlinear FIFs may have more flexibility and applicability.

### Connection 9

1. The well-known FIS in theory and applications is generated by an IFS of the form K w m , n : m = 1 2 M n = 1 2 N under some conditions, where the maps are transformations of the special structure

w m , n x y z = u m x v n y F m , n x y z = a m x + b m c n y + d n e m , n x + f m , n y + g m , n xy + d m , n x y z + h m , n ,

where d m , n x y < 1 for all x y D R 2 . Then for all x y z , x y z K ,

F m , n x y z F m , n x y z = d m , n x y z d m , n x y z max x y D d m , n x y z z .

That is, each w m , n x y z is chosen so that function F m , n x y z is a Banach contraction with respect to the third variable. So, the existence of bivariable FIFs follows from Banach’s fixed point theorem. In fact, in [22], since for all x y D R 2 , d m , n x y s m , n and 0 s m , n < 1 , we can see that each w m , n x y z is chosen so that function F m , n x y z is Banach contraction with respect to the third variable. Also in [21], since

d m , n x y = λ m , n x x 0 x M x y y 0 y N y

and

λ m , n < 16 x M x 0 2 y N y 0 2 ,

we can see that max x y D d m , n x y < 1 , and so each w m , n x y z is chosen so that function F m , n x y z is a Banach contraction with respect to the third variable.

1. In Theorem 2.3, for all x y z , x y z K R 3 ,

F m , n x y z F m , n x y z = d m , n x y s m , n z s m , n z s m , n z s m , n z φ z z .

That is, each F m , n x y z is Rakotch contraction with respect to the third variable. So, each w m , n x y z is chosen so that the function F m , n x y z is a Rakotch contraction with respect to the third variable.

### Connection 10

In view of a φ -contraction, the connections between the coefficients of variable z are obtained as follows:

1. In the affine FIS (cf. [22]), for all t > 0 ,

φ t max m = 1 , 2 , , M max n = 1 , 2 , , N d m , n t ,

where d m , n < 1 for all m = 1 , 2 , , M and n = 1 , 2 , , N .

1. In the FIS with vertical scaling factors as function (cf. [21]), for all t > 0 ,

φ t max m = 1 , 2 , , M max n = 1 , 2 , , N max x I d m , n x t ,

where d m , n x is Lipschitz function defined on I satisfying sup x I d m , n x < 1 for all m = 1 , 2 , , M and n = 1 , 2 , , N .

### Connection 11

The continuity of bivariable FIFs differ from the continuity of univariable FIFs.

1. The graphs of linear univariable FIFs are always continuous curves.

2. There are bivariable discontinuous functions that interpolate the given data; (see for instance [23], p. 630, 631).

3. Theorem 2.3 ensures that attractors of constructed IFSs are graphs of some bounded functions which interpolate the given data, but these graphs (i.e., the graphs of bivariable FIFs) are not always continuous surfaces. Some continuity conditions of bivariable FIFs are given explicitly by Corollary 2.1 and Corollary 2.2.

### Connection 12

The key difficulty in constructing fractal interpolation surfaces (or volumes) involves ensuring continuity. Another important element necessary in modelling complicated surfaces of this type is the existence of the contractivity, or vertical scaling, factors.

1. In order to ensure continuity of a fractal interpolation surface, in [22], the interpolation points on the boundary was assumed collinear, whereas in [21], vertical scaling factors as (continuous) ‘contraction functions’ were used.

2. A new bivariable fractal interpolation function by using the Matkowski fixed point theorem and the Rakotch contraction is presented in [18]. In order to ensure the continuity of nonlinear FIS, the coplanarity of all the interpolation points on the boundaries instead of collinearity of interpolation points on the boundary was assumed in [18], whereas in [12], vertical scaling factors as (continuous) ‘contraction functions’ were used.

### Connection 13

1. In Theorem 2.2, we can see that

a n = x n x n 1 x N x 0 , b n = x N x n 1 x 0 x n x N x 0 c n = y n y n 1 x N x 0 d n x N s n y N d n x 0 s n y 0 x N x 0 , f n = x N y n 1 x 0 y n x N x 0 x N d n x 0 s n y 0 x 0 d n x N s n y N x N x 0 .

1. In Corollary 2.1, we can see that

a m = x m x m 1 x M x 0 , b m = x M x m 1 x 0 x m x M x 0 , c n = y n y n 1 y N y 0 , d n = y N y n 1 y 0 y n y N y 0 , g m , n = z m , n z m 1 , n z m , n 1 z m 1 , n 1 x M x 0 y N y 0 , e m , n = y N z m , n 1 z m 1 , n 1 y 0 z m , n z m 1 , n x M x 0 y N y 0 , f m , n = x M z m 1 , n z m 1 , n 1 x 0 z m , n z m , n 1 x M x 0 y N y 0 , h m , n = x 0 y 0 z m , n x 0 y N z m , n 1 x M y 0 z m 1 , n + x M y N z m 1 , n 1 x M x 0 y N y 0 s m , n z M , N .

1. In Corollary 2.2, we can see that (compare with above coefficients)

g m , n = z m , n z m 1 , n z m , n 1 z m 1 , n 1 x M x 0 y N y 0 , e m , n = y N z m , n 1 z m 1 , n 1 y 0 z m , n z m 1 , n x M x 0 y N y 0 , f m , n = x M z m 1 , n z m 1 , n 1 x 0 z m , n z m , n 1 x M x 0 y N y 0 , h m , n = x 0 y 0 z m , n x 0 y N z m , n 1 x M y 0 z m 1 , n + x M y N z m 1 , n 1 x M x 0 y N y 0 .

Figures 1(a) and 2(a) are associated with Banach contractions, whereas Figures 1(b) and 2(b) are not necessarily associated with Banach contractions.

## 4. Conclusions and further work

We reviewed nonlinear fractal interpolation functions by using the Geraghty fixed point theorem instead of the Banach fixed point theorem (or the Rakotch fixed point theorem) since Banach contraction (or Rakotch contraction) is a special case of Geraghty contraction. Theorems 2.1, 2.2 and 2.3 ensure that attractors of constructed nonlinear iterated function systems are graphs of some continuous functions which interpolate the given data. In particular, Examples 1, 2 and 3 show that our results remain still true under essentially weaker conditions on the maps of iterated function systems. The methods presented here can be directly extended to piecewise fractal interpolation functions that are based on recurrent IFS. A premise for future work is to extend these methods to hidden-variable fractal interpolation surfaces as well as to identify the parameters of such surfaces.

## References

1. 1. Barnsley MF. Fractal functions and interpolation. Constructive Approximation. 1986;2:303-329
2. 2. Barnsley MF. Fractals Everywhere. 3rd ed. New York: Dover Publications, Inc.; 2012
3. 3. Chand AKB, Kapoor GP. Generalized cubic spline fractal interpolation functions. SIAM Journal on Numerical Analysis. 2006;44(2):655-676
4. 4. Viswanathan P, Chand AKB, Navascués MA. Fractal perturbation preserving fundamental shapes: Bounds on the scale factors. Journal of Mathematical Analysis and Applications. 2014;419:804-817
5. 5. Vijender N. Bernstein fractal trigonometric approximation. Acta Applicandae Mathematicae. 2018;159:11-27
6. 6. Barnsley MF, Elton JH, Hardin DP, Massopust P. Hidden variable fractal interpolation functions. SIAM Journal on Mathematical Analysis. 1989;20:1218-1242
7. 7. Chand AKB, Kapoor GP. Hidden variable bivariate fractal interpolation surfaces. Fractals. 2003;11(3):277-288
8. 8. Chand AKB, Kapoor GP. Stability of affine coalescence hidden variable fractal interpolation functions. Nonlinear Analysis. 2008;68:3757-3770
9. 9. Barnsley MF, Massopust PR. Bilinear fractal interpolation and box dimension. Journal of Approximation Theory. 2015;192:362-378
10. 10. Massopust PR. Fractal Functions, Fractal Surfaces and Wavelets. San Diego: Academic Press; 1994
11. 11. Drakopoulos V, Manousopoulos P. On non-tensor product bivariate fractal interpolation surfaces on rectangular grids. Mathematics. 2020;8(4):525. DOI: 10.3390/math8040525
12. 12. Ri S. New types of fractal interpolation surfaces. Chaos, Solitons and Fractals. 2019;119:291-297
13. 13. Rhoades BE. A comparison of various definitions of contractive mappings. Transactions of the American Mathematical Society. 1977;226:257-290
14. 14. Strobin F, Swaczyna J. On a certain generalisation of the iterated function system. Bulletin of the Australian Mathematical Society. 2013;87(1):37-54
15. 15. Xu S, Cheng S, Zhou Z. Reich’s iterated function systems and well-posedness via fixed point theory. Fixed Point Theory and Applications. 2015;2015(1):71. DOI: 10.1186/s13663-015-0320-7
16. 16. Rakotch E. A note on contractive mappings. Proceedings of the American Mathematical Society. 1962;13:459-465
17. 17. Ri S. A new nonlinear fractal interpolation function. Fractals. 2017;25(6):1750063 (12 pages)
18. 18. Ri S. A new nonlinear bivariate fractal interpolation function. Fractals. 2018;26(4):1850054 (14 pages)
19. 19. Jachymski J, Jóźwik I. Nonlinear contractive conditions: A comparison and related problems. Banach Center Publications. 2007;77:123-146
20. 20. Wang H-Y, Jia-Shan Y. Fractal interpolation functions with variable parameters and their analytical properties. Journal of Approximation Theory. 2013;175:1-18
21. 21. Feng Z, Feng Y, Yuan Z. Fractal interpolation surfaces with function vertical scaling factors. Applied Mathematics Letters. 2012;25(11):1896-1900
22. 22. Dalla L. Bivariate fractal interpolation functions on grids. Fractals. 2002;10(1):53-58
23. 23. Xie H, Sun H. The study on bivariate fractal interpolation functions and creation of fractal interpolated surfaces. Fractals. 1997;5(4):625-634

Written By

SongIl Ri and Vasileios Drakopoulos

Submitted: December 23rd, 2019 Reviewed: April 27th, 2020 Published: August 21st, 2020