Abstract
In this paper we consider the (inverse) problem of determining the iterated function system (IFS) which produces a shaped fractal interpolant. We develop a new type of rational IFS by using functions of the form , where are cubics and are preassigned quadratics having 3shape parameters. The fixed point of the developed rational cubic IFS is in , but its derivative varies from a piecewise differentiable function to a continuous nowhere differentiable function. An upper bound of the uniform error between the fixed point of a rational IFS and an original function is deduced for the convergence results. The automatic generations of the scaling factors and shape parameters in the rational IFS are formulated so that its fixed point preserves the positive/monotonic features of prescribed data. The presence of scaling factors provides additional freedom to the shape of the fractal interpolant over its classical counterpart in the modeling of discrete data.
1 Introduction
Setting a novel platform for the approximation of natural objects such as trees, clouds, feathers, leaves, flowers, landscapes, glaciers, galaxies, and torrents of water, Mandelbrot [1] introduced the term fractal in the literature. Since fractals capture the nonlinear structures of various objects effectively, the fractal geometry has been successfully used in different problems in applied sciences and engineering [14]. The iterated function system (IFS) was introduced by Hutchinson [5] for the construction of various types of fractal sets, and popularized by Barnsley [6]. An IFS is a dynamical system consisting of a finite collection of continuous maps. Based on the IFS theory, Barnsley [7] constructed a class of functions that are known as FIFs. The graph of a FIF is the fixed point of an IFS. Also a FIF is the fixed point of the ReadBajraktarević operator on a suitable function space. Common features between a FIF and a piecewise polynomial interpolation are that they are geometrical in nature, and they can be computed rapidly, but the main difference is the fractal character, i.e., a FIF satisfies a functional relation related to the selfsimilarity on smaller scales. In the direction of smooth fractal curves, Barnsley and Harington [8] initiated the construction of a restricted class of differentiable FIF or FIF that interpolates the prescribed data if the values of , , at the initial end point of the interval are given, where Φ is the original function. This method is based on the recursive nature of an algorithm, and specifying the boundary conditions similar to the classical splines was found to be quite difficult to handle in this construction. The fractal splines with general boundary conditions have been studied recently [913] by restricting their IFSs parameters suitably.
The motivation of this work is the research on different types of splines by several authors; see, for instance, Schmidt and Heß [14], Fritsch and Carlson [15], Schumaker [16], and Brodlie and Butt [17], and references therein. The uniqueness of spline representation for a given data set turns out to be a disadvantage for shape modification problems. The use of rational functions with the shape parameters was introduced by Späth [18] to preserve different geometric properties attached to a given set of data. Rational interpolants are often used in data visualization problems due to their excellent asymptotic properties, capability to model complicated smooth structures, better interpolation properties, and excellent extrapolating powers. Gregory and Delbourgo [19] introduced the rational cubic spline with one family of shape parameters, and this work inspired a large amount of research in shapepreserving rational spline interpolations, see [20,21] and references therein.
In this paper, we introduce the rational cubic IFS with 3shape parameters in each subinterval of the interpolation domain such that its fixed point generalizes the corresponding classical rational cubic spline functions [20]. The developed rational cubic spline FIF is bounded, and is unique by fixed point theory for a given set of scaling factors and shape parameters. Because of the recursive nature of FIF, the necessary conditions for monotonicity on the derivative values at knots alone may not ensure the monotonicity of a rational cubic fractal interpolant for a given monotonic data. Based on the appropriate condition on the rational IFS parameters: (i) the scaling factors that depend only on given data, and (ii) the shape parameters that depend on both the interpolation data and scaling factors, we construct the shapepreserving rational cubic FIFs for a prescribed positive and/or monotonic data. By varying the scaling factors (within the shapepreserving interval) and shape parameters (according to the conditions derived in our theory), we can make the fixed point of a rational cubic IFS more pleasant and suitable for aesthetic requirements in a modeling problem. The proposed method is suitable for the shapepreserving interpolation problems where a data set originates from an unknown function and its derivative is a continuous nowhere differentiable function, for instance, the motion of single inverted pendulum in nonlinear control theory [22].
Comparison of the proposed rational cubic FIF over some existing schemes:
• When all the scaling factors are zero, the proposed rational cubic FIF reduces to the classical rational cubic interpolant [20], see Remark 1, Section 3.
• To generate shapepreserving interpolants, our construction does not need additional knots in contrast to methods due to Schumaker [16] and Brodlie and Butt [17], which require additional knots for the shapepreserving interpolants.
• The classical interpolants [15,23] are suitable only for monotonicity interpolation whereas the proposed rational cubic FIF is suitable for both monotonicity and positivity interpolation. Moreover, the rational quadratic interpolant [23] is a special case of our rational cubic FIF for the particular choice of the scaling factors and shape parameters, see Remark 2, Section 3.
• For given monotonic data, the monotonic curve generated by the rational quadratic interpolant [23] is unique for fixed shape parameters, whereas for the same monotonic data an infinite number of monotonic curves will be obtained using our rational FIF by suitable modifications in the associated scaling factors. Thus, when the shape parameters are incapable to change the shape of an interpolant in given intervals, then the scaling factors can be used to alter the shape of the interpolant in our method.
• Where monotonicity is concerned, our construction does not need an additional condition on derivatives at knots except for the necessary conditions. But the construction of Fritsch and Carlson [15] needs some restrictions on derivatives at knots apart from the necessary conditions for the same problem.
• The derivatives of the shapepreserving interpolants [1517,23] are piecewise smooth, whereas the derivative of our rational cubic FIF may be piecewise smooth to a nondifferentiable function according to the choice of the scaling factors. Owing to this special feature, the proposed method is preferable over the classical shapepreserving interpolants when the approximation is taken for data originating with an unknown function having a shape with fractality in .
This paper is organized as follows. In Section 2, the general constructions of fractal interpolants and rational cubic FIFs based on IFSs are summarized. Section 3 is devoted to the construction of a suitable rational IFS so that its fixed point is the desired interpolant that can be used for shape preservation. Then we deduce an upper bound of the uniform error bound between the original function and the rational cubic FIF. The fixed point of this rational IFS does not follow any shape constraints. The restrictions on the rational IFS parameters are deduced for a positivity shape in Section 4, and the results are illustrated with suitably chosen examples. In Section 5, the monotonicity problem is considered through the developed rational cubic IFS.
2 IFS for fractal functions
Let be a partition of . Let be the value of original function at , . Denote , , , and let D be a compact subset of ℝ such that , . Let , , be the contractive homeomorphisms such that
It is easy to verify that is a just touching hyperbolic IFS whose unique fixed point is
Let , , , be the continuous realvalued functions on C such that
and , , are the suitable continuous functions. Now define the functions , , as for every . Then is called an IFS related to a given interpolation data . According to [7], the IFS ℐ has a unique fixed point G which is the graph of a continuous function , , . The function ϕ is called a FIF generated by the IFS ℐ, and it takes the form
The existence of a spline FIF based on a polynomial IFS is given in [8]. We have extended this result to the rational IFS with 3shape parameters in the following.
Theorem 1Letbe a given data set, whereare the slope at, and () are thekth derivative values atfor. Consider the rational IFS, wheresatisfies equation (1), is a suitable compact subset of ℝ. , , is a polynomial containingarbitrary constants, andis a nonvanishing quadratic polynomial with 3shape parameters in each subinterval defined onI, and, . Let, whererepresents thederivative ofwith respect tox. With the setting, , , if
then the fixed point of the rational IFSis the graph of therational FIF.
Proof Suppose . Now is a complete metric space, where is the metric on induced by the norm on . Define the ReadBajraktarević operator U on as
Since , the conditions and (4) imply that U is a contractive operator on . The fixed point ψ of U is a fractal function that satisfies the functional equation:
Using equation (4) in equation (7), we get the following system of equations for :
When all arbitrary constants in are determined from equation (8), then exists. By using similar arguments as in [7], it can be shown that IFS has a unique fixed point, and that it is the graph of the rational FIF . □
3 Rational cubic IFS
The construction of the desired rational cubic IFS is given in Section 3.1 such that its fixed point is used for shape preservation in the sequel. The error analysis of the fixed point of rational cubic IFS with an original function is studied in Section 3.2 for convergence results.
3.1 Construction
In the proposed rational cubic IFS, we assume () are the rational functions with 3shape parameters, whose denominators are preassigned quadratics. Based on Theorem 1, with , consider the following fixed point equation:
, , , and are arbitrary constants, and , , and are the shape parameters such that . From this condition, it is easy to see that for all . To make the fixed point ψ a interpolant, the following Hermite interpolatory conditions are imposed:
After evaluation of , , , and using the above Hermite interpolatory conditions, we get the desired rational cubic FIF:
where
Now it is easy to see that rational cubic FIF ψ is the fixed point of the following rational cubic IFS:
where , , are evaluated by using equation (1),
The fixed point ψ of the above rational cubic IFS is unique for every fixed set of scaling factors and shape parameters. Thus by taking different sets of scaling and shape parameters, we can generate an infinite number of fixed points for the above rational cubic IFS. In most applications, the derivatives () are not given, and hence they must be calculated either from the given data or by using numerical approximation methods [24].
Remark 1 If for , then the rational cubic FIF (10) coincides with the corresponding classical rational cubic interpolation function S as
described in the literature [20] with , .
Remark 2 Substituting and for in equation (10), we have
where
After some rigorous calculations, we have found that
where
Now from equation (12), we conclude that for the above choice of , , and , our rational cubic FIF ψ reduces to a monotonicity preserving rational quadratic FIF [25] constructed by our group. Also it is easy to verify that, if , and , , then the rational cubic FIF reduces to the rational quadratic function as in [23].
3.2 Error analysis of fixed point of rational cubic IFS
Theorem 2LetψandS, respectively, be the fixed point of rational cubic IFS (11) and the classical rational cubic function with respect to the dataobtained from the original function. Denote, , . Let the shape parameters satisfy. Then
Proof Since the coefficients of in equation (9) depend on , we can write . From equation (5), the ReadBajraktarević operator (cf. Section 2 with ) is rewritten as for ,
Let ξ and e be the nonzero and zero vectors in respectively. If , is the only function of x for , then the classical rational cubic interpolant is the fixed point of . Let us assume that ψ is a fixed point of a rational cubic IFS (11) associated with a nonzero scale vector ξ. Consequently, ψ is the fixed point of . From equation (14), it is easy to verify that is a contractive operator for a fixed scaling vector ξ:
Now,
Using the mean value theorem for functions of several variables, there exists such that for ,
Using equation (17) in equation (16), we have
Now we wish to calculate the bounds of each term in the righthand side of equation (18). From Remark 1, it is easy to see that
where
Now using equation (20) in equation (19), we get .
Since the above inequality is true for all , we get the following estimation:
Since is independent of , it easy to see that
where
By using similar arguments as used in the estimation of , we have found that
By using equations (21) and (22) in equation (18), we have
Since the above inequality is true for ,
Combining equations (15) and (23) with the inequality
we get
From equation (24), it is evident that for , , the fixed point of rational cubic IFS (11) coincides with the corresponding classical rational cubic interpolant.
Since the original function , it is known that [20]
Therefore, using equations (24)(25) together with the inequality , we get the bound for , and it completes the proof of theorem. □
Corollary 1 (Convergence results)
Assume thatandare bounded as. Then we have the following results:
(i) Since, we conclude from equation (13) that the fixed point of rational cubic IFS equation (11) converges uniformly to the original function Φ as.
(ii) Again from the error estimation (13), () convergence can be obtained if the derivative values are available such that (), and the scaling factors are chosen as () for.
4 Positivity preserving rational cubic FIF
The rational cubic fractal interpolation function developed in Section 3 has deficiencies as far as the positivity preserving issue is concerned. Because of the recursive nature of FIFs, we assume all the scaling factors are nonnegative so that it is easy to derive the sufficient conditions for a positive fixed point of the rational cubic IFS (11). It requires one to assign appropriate restrictions on the scaling factors and shape parameters , and , for , so that the positivity feature of a given set of positive data is preserved in the fixed point of the rational cubic IFS (11). In Section 4.1, the suitable restrictions are developed on the scaling factors and shape parameters for a positivity preserving rational cubic spline FIF. The importance of suitable restrictions on the rational IFS parameters is illustrated in Section 4.2.
4.1 Restrictions on IFS parameters for positivity
Theorem 3Letbe a given positive data. If
(i) the scaling factors, , are selected as
(ii) with, the shape parameters, , , are chosen as
where, , , , then for fixed, , , (), the unique fixed pointψof the rational IFS (11) is positive.
Proof From equation (10), we have
It is easy to verify that using equation (2), if , , the sufficient conditions for for all are for all , . If we assume , , and , then it is easy to see that for any . Thus the initial conditions on the scaling factor and shape parameters are , and , , , respectively, for . Including the initial conditions on the shape parameters, we have , . Thus our problem reduces to finding conditions on the scaling factors and shape parameters for which for all . From equation (10), is rewritten as
where
By substituting in equation (28), for all is equivalent to for all , where , , , .
From [14], we have for all if and only if , where
Let , then we have , , , . Now , if and only if , , respectively. Hence, the restriction on the scaling factor is
If , is true from equation (29), and in this case can be chosen arbitrarily. Otherwise, , we have . Similarly is true when (i) , arbitrary (ii) , . Another set of restrictions on , , , and can be derived if . But we have not considered it here due to the complexity involved in the calculations. The above discussions yield equation (27).
Therefore, whenever equations (29) and (27) are true. Now it is easy to see that the fixed point of the rational cubic IFS (11) is positive if the scaling factors and shape parameters involved in the IFS (11) satisfy equations (26) and (27), respectively. □
4.2 Examples and discussion
In order to demonstrate the positive interpolation using our rational cubic IFS, consider the positive data set , which is taken from a function (see Figure 1(a)) with an irregular derivative function as described in Figure 2(a). Such type of data arises in the motion of a single inverted pendulum in the field of nonlinear control theory [22]. The position of the cart can be taken as smooth and positive in a short time interval with an irregular cart velocity. To approximate such data, we have employed the rational cubic IFS (11). The derivative values at the knots are approximated by the arithmetic mean method [24] as , , , , , , and . The scaling factors are constrained as , , , , , by equation (26) with a choice of .
Figure 1. Illustration of positive rational fractal interpolants with shape parameters.
Figure 2. Derivatives of positive rational fractal interpolants and classical interpolant.
The IFS parameters of the original function are given in Table 1, and aesthetic modifications are illustrated by varying the scaling factors and shape parameters. In order to explain the sensitiveness of a rational cubic FIF with respect to the scaling factors, we have taken a fixed set of shape parameters in the construction of Figures 1(b)(c), see Table 1. By comparing Figure 1(b) with Figure 1(a), we observe that the fractal curve pertaining to the first subinterval converges to a convex shape as , and changes in other subintervals are negligible. By comparing the shapes of Figure 1(a) and Figure 1(c), we notice perceptible variations in the second subinterval, and variations in other subintervals are negligible. By comparing Figures 1(d)(e) with Figure 1(a) and Figure 1(c), respectively, we can observe the sensitivity of the positive FIF with respect to its shape parameters. Finally, we have constructed the classical rational cubic interpolant in Figure 1(f) with the zero scaling vector. From the above discussion, we conclude that the effects due to the scaling factors , , and shape parameters , are very local in nature for the given positive data set.
Table 1. Rational IFS parameters for positive fractal interpolants
From equations (4) and (10), interpolates the data . In this example, the interpolation data for is . The derivative functions of rational cubic FIFs in Figures 1(a)(e) are constructed in Figures 2(a)(e), respectively, and they are typically irregular fractal functions close to a continuous function, but at least they differ from a piecewise differentiable function. We have calculated the uniform errors between this original function Φ in Figure 1(a) and the rational cubic FIFs in Figures 1(b)(f) (see Table 2). Also we have calculated the uniform errors between their derivatives (see Table 2). The effects of the scaling factors are very prominent in the first subinterval of Figure 1(b), but they also render major effects in its derivative (see Figure 2(b) and Table 2). The effects of the scaling factor are prominent in the second and third subintervals of Figure 2(c) in comparison with Figure 2(a). Hence the scaling factor is moderately local in nature in the derivative of rational cubic FIF. The shape parameters and produce similar effects (see the corresponding figures and Table 2). The rational fractal functions in Figures 2(a)(e) are irregular in nature over the interval , whereas the derivative of a classical interpolant is piecewise differentiable in the interval (see Figure 2(f)). Comparing the uniform distances in Table 2, if the original function is smooth and positive but its derivative is very irregular, then our rational cubic IFS is an ideal tool for approximating such a function instead of the classical rational cubic interpolant whose derivative is a piecewise smooth function.
Table 2. Uniform errors betweenΦand rational fractal interpolants, and their derivatives
5 Monotonicity preserving rational cubic FIF
The fixed point ψ of a rational cubic IFS may not preserve the monotonic feature of a given set of monotonic data. For an automatic generation of rational IFS parameters, we restrict them in Section 5.1, and the results are implemented in Section 5.2 through suitable examples.
5.1 Restrictions on IFS parameters for monotonicity
Theorem 4Letbe a given monotonic data. Let the derivative values satisfy the necessary conditions for monotonicity, namely
If (i) the scaling factors, , are chosen as
(ii) the shape parameters, , and, , are selected as
then for a fixedξ, , , (), the unique fixed pointψof the rational cubic IFS (11) is monotonic in nature.
Proof Differentiating equation (10) with respect to x, after some mathematical manipulations, we get
where
Due to the recursive nature of rational fractal function (33), the necessary conditions (30) are not sufficient to ensure the monotonicity of fixed point ψ of the rational cubic IFS (11). We impose additional restrictions on the scaling factors , and shape parameters , , and , , so that these conditions together with the necessary conditions (30) yield the monotonic feature of the fixed point ψ of the IFS (11).
Case I: Monotonically increasing data
Suppose is a given monotonically increasing data set. Due to the recursive nature of IFS and equation (2), it is assumed that all the scaling factors , , are nonnegative for a monotonic fixed point of rational cubic IFS (11). For , , which is monotone on (choose ). Otherwise for , the sufficient conditions for the monotonicity of the fixed point of rational cubic IFS (11) are , . From equation (33),
Similarly,
From equation (33), is rewritten as
Without loss of generality, assume that
i.e., . We search for sufficient conditions that make . For this purpose, we make each term in nonnegative. The selection of with respect to equations (34) and (36) gives and , respectively. Now it remains to make and . In these two inequalities, the product of the shape parameters is involved. Therefore these inequalities are true if we restrict the shape parameters , , and , , respectively, as in equation (32).
Justification for equation (32)
Let be negative, then from equations (32) and (34)(36), we can conclude that is negative. Therefore, and . Similarly, it can be shown that being positive gives similar results.
The above discussion led to the following procedure to make : first choose the scaling factors with respect to equations (34)(36), then select the shape parameters according to equation (32).
Again from equation (33), is rearranged as . Similarly, it is easy to verify that equations (32) and (34)(36) are sufficient for . For simplicity, denote , , . From equation (33), and with the above notations, is rewritten as
Substituting (see equation (32)) in the above expression, we get
From the final expression of , it is easy to verify that equations (32) and (34)(36) are sufficient for . Hence we have proved that the fixed point ψ of the rational cubic IFS (11) is monotonically increasing over , if the scaling factors and shape parameters are chosen according to equation (31) and equation (32), respectively. In the case of , the fixed point of the rational cubic IFS (11) is a constant throughout that subinterval with the value , and .
Case II: Monotonically decreasing data
Suppose is a given monotonically decreasing data set. It is easy to see that the sufficient conditions for monotonicity of equation (10) on are , . As explained in Case I, it is easy to verify that selections of the scaling factors and shape parameters according to equation (31) and (32), respectively, are sufficient for , .
Therefore from the arguments in Case I and Case II, we conclude that if the scaling factors and shape parameters are chosen according to (31) and (32), respectively, then the fixed point ψ of the rational cubic IFS (11) is monotone for given monotonic data. □
Remark 3 Convergence results in Corollary 1 are valid for the shapepreserving rational cubic FIFs.
5.2 Examples and discussion
We construct the rational cubic fractal interpolation functions (RCFIFs) for the standard increasing Akima data [26]. The rational cubic FIFs are generated iteratively (Figures 3(a)(i)) as the fixed points of rational cubic IFS (11). Since for , for , and there is no need to choose and , consequently there is no need to calculate for using equation (32). The derivatives values () are approximated by the arithmetic mean method [24] as , , , , , , , , , , and . Let in equation (31). The scaling factors are restricted as , , , , to capture the monotonicity of the Akima data. A standard rational cubic FIF Φ in Figure 3(a) is generated with a suitable choice of the scaling factors (see Table 3). By comparing Figure 3(b) with Figure 3(a), we observe that the graph of the rational cubic FIF in the subinterval converges to a convex shape as , and the changes in other intervals are negligible. By analyzing the shapes of Figure 3(a) and Figure 3(c), we observe that visually pleasing effects are produced in the ninth subinterval. By analyzing Figure 3(d) with respect to Figure 3(a), we have found excellent variations in . The individual effects of and , respectively, from Figure 3(b) and Figure 3(c) are reflected in Figure 3(d), and thereby we conclude that all the nonzero scaling factors are very much local in nature for this Akima data. Next to visualize the effects of change in the shape parameters, we construct the monotonically increasing rational cubic FIFs in Figures 3(e)(h). By comparing the shapes of Figure 3(e) and Figure 3(a), we notice perceptible variations in the ninth subinterval, and variations in other subintervals are negligible. Again analyzing Figure 3(a), Figure 3(c), and Figure 3(e), we observe that to get appropriate deviations in the rational cubic FIF in the ninth subinterval, one has to vary the scaling parameter , the shape parameters and/or suitably.
Figure 3. Monotonicity preserving rational cubic FIFs and their derivatives.
Table 3. Rational IFS parameters for monotonic fractal interpolants
Next a monotonic rational FIF in Figure 3(f) is constructed as per the data in Table 3 and the variations in of Figure 3(f) with respect to Figure 3(a) are more evident than those of Figure 3(f) with respect to Figure 3(b). So, we can say that the scaling factor is dominant over the shape parameters and in this case at . Hence, for major changes at , one has to modify , and for minor changes (or fine tuning), one has to alter and/or . This observation is useful for aesthetic requirements in various engineering design problems. By analyzing Figures 3(g)(h) with respect to Figure 3(c), we have noticed that the graphs of rational cubic FIFs in the ninth subinterval in Figures 3(g)(h) are concave and convex, respectively. Finally, we construct the classical rational cubic interpolant in Figure 3(i) with respect to the shape parameters of Figure 3(g), and for all i. Since the shape parameters and are the same in Figure 3(g) and Figure 3(i), there is some visual similarity between these two curves in , whereas the same effects are missing in Figure 3(h) due to a variation in . Also, one gets a classical rational cubic interpolant which is similar to Figure 3(h) in our fractal scheme, whenever for all i, and the shape parameters are chosen according to Figure 3(h). The presence of scaling factors in a monotonic rational FIF with shape parameters gives an additional advantage in the choice of interpolant over the classical interpolants with shape parameters. Our construction gives an extra freedom for aesthetic modifications in local shape over the classical rational cubic interpolants to an user. For a qualitative study of the derivatives of monotonic fractal interpolants, the readers are invited to check the effects of rational IFS parameters in Figures 3(a)(g). The uniform errors between monotonic fractal interpolants and their derivatives are given in Table 4 to show the importance of our rational cubic IFS (11).
Table 4. Uniform errors betweenΦand rational fractal interpolants, and their derivatives
From the examples in Sections 45, it is observed that proper interactive adjustments of the scaling factors and shape parameters give us a wide variety of positivity and/or monotonicity preserving fixed points of our rational cubic IFS (11) that can be used in various scientific and engineering problems for aesthetic modifications. In order to get an optimal choice of the fixed point of our rational cubic IFS (11), one can employ a genetic algorithm interactively until the desired accuracy is obtained with the original function.
6 Conclusion
A new type of rational cubic IFS with 3shape parameters is introduced in this work such that its fixed point can be used for shaped data. The developed FIF in this paper includes the corresponding classical rational cubic interpolant as a special case. An upper bound of uniform error between the rational cubic FIF ψ and an original function Φ in is estimated, and consequently we have found that ψ converges uniformly to Φ as . When the accurate derivatives of , are available, and the scaling factors are chosen as , , it is possible to get () convergence for the rational cubic FIF. Automatic data dependent restrictions are derived on the scaling factors and shape parameters of rational cubic IFS so that its fixed point preserves the positivity or monotonicity features of a given set of data. The effects of a change in the scaling factors and shape parameters on the local control of the shape of rational cubic FIF are demonstrated through various examples. Our rational cubic FIFs are more flexible and more suitable for shape related problems in computer graphics, CAD/CAM, CAGD, medical imaging, finance, and engineering applications, and apply equally well to data with or without derivatives. In particular, the proposed method will be an ideal tool in shapepreserving interpolation problems where the data set originates from a positive and/or monotonic function , but its derivative is a continuous and nowhere differentiable function.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed equally to the manuscript and typed, read, and approved the final manuscript.
Acknowledgements
AKBC is thankful to the Department of Science and Technology, Govt. of India for the SERC DST Project No. SR/S4/MS: 694/10. The authors are grateful to the anonymous referees for the valuable comments and suggestions, which improved the presentation of the paper.
References

Mandelbrot, BB: Fractals: Form, Chance and Dimension, Freeman, San Francisco (1977)

West, BJ: Fractals in Physiology and Chaos in Medicine, World Scientific, Singapore (1990)

Yang, YJ, Baleanu, D, Yang, XJ: Analysis of fractal wave equations by local fractional Fourier series method. Adv. Math. Phys.. 2013, (2013) Article ID 632309

Hutchinson, JE: Fractals and self similarity. Indiana Univ. Math. J.. 30, 713–747 (1981). Publisher Full Text

Barnsley, MF: Fractals Everywhere, Academic Press, Orlando, Florida (1988)

Barnsley, MF: Fractal functions and interpolation. Constr. Approx.. 2, 303–329 (1986). Publisher Full Text

Barnsley, MF, Harrington, AN: The calculus of fractal interpolation functions. J. Approx. Theory. 57, 14–34 (1989). Publisher Full Text

Chand, AKB, Kapoor, GP: Generalized cubic spline fractal interpolation functions. SIAM J. Numer. Anal.. 44(2), 655–676 (2006). Publisher Full Text

Chand, AKB, Navascués, MA: Generalized Hermite fractal interpolation. Rev. Acad. Cienc. Exactas, Fís.Quím. Nat. Zaragoza. 64, 107–120 (2009)

Chand, AKB, Navascués, MA: Natural bicubic spline fractal interpolation. Nonlinear Anal.. 69, 3679–3691 (2008). Publisher Full Text

Chand, AKB, Viswanathan, P: Cubic Hermite and cubic spline fractal interpolation functions. AIP Conf. Proc.. 1479, 1467–1470 (2012)

Chand, AKB, Viswanathan, P: A constructive approach to cubic Hermite fractal interpolation function and its constrained aspects. BIT Numer. Math.. 53(4), 841–865 (2013). Publisher Full Text

Schmidt, JW, Heß, W: Positivity of cubic polynomials on intervals and positive spline interpolation. BIT Numer. Math.. 28, 340–352 (1988). Publisher Full Text

Fritsch, FN, Carlson, RE: Monotone piecewise cubic interpolation. SIAM J. Numer. Anal.. 17, 238–246 (1980). Publisher Full Text

Schumaker, LL: On shape preserving quadratic spline interpolation. SIAM J. Numer. Anal.. 20, 854–864 (1983). Publisher Full Text

Brodlie, KW, Butt, S: Preserving convexity using piecewise cubic interpolation. Comput. Graph.. 15, 15–23 (1991). Publisher Full Text

Späth, H: Spline Algorithms for Curves and Surfaces, Utilias Mathematica Pub. Inc., Winnipeg (1974)

Delbourgo, R, Gregory, JA: Shape preserving piecewise rational interpolation. SIAM J. Sci. Stat. Comput.. 6, 967–976 (1985). Publisher Full Text

Sarfraz, M, ALMuhammed, M, Ashraf, F: Preserving monotonic shape of the data using piecewise rational cubic functions. Comput. Graph.. 21, 5–14 (1997). Publisher Full Text

Sarfraz, M, Hussain, MZ, Hussain, M: Shapepreserving curve interpolation. Int. J. Comput. Math.. 89, 35–53 (2012). Publisher Full Text

Kwakernaak, H, Sivan, R: Linear Optimal Control Systems, WileyInterscience, New York (1972)

Gregory, JA, Delbourgo, R: Piecewise rational quadratic interpolation to monotonic data. IMA J. Numer. Anal.. 2, 123–130 (1982). Publisher Full Text

Gregory, JA, Delbourgo, R: Determination of derivative parameters for a monotonic rational quadratic interpolant. IMA J. Numer. Anal.. 5, 397–406 (1985). Publisher Full Text

Chand, AKB, Vijender, N: Monotonicity preserving rational quadratic fractal interpolation functions. Adv. Numer. Anal.. 2014, (2014) Article ID 504825

Akima, H: A new method of interpolation and smooth curve fitting based on local procedures. J. Assoc. Comput. Mach.. 17, 589–602 (1970). Publisher Full Text