---
title: Nonlinear Fractal Histopolation Function
url: https://www.emergentmind.com/topics/nonlinear-fractal-histopolation-function
type: topic
---

# Nonlinear Fractal Histopolation Function

A nonlinear fractal histopolation function is a self-referential bounded function associated with a prescribed histogram such that its integrals over a partition match the histogram areas, while continuity is not required. In the formulation introduced in "Nonlinear Fractal Histopolation Function" [2509.18657], the function is constructed as the fixed point of a nonlinear Read–Bajraktarević operator induced by an iterated function system (IFS) built from affine domain maps and nonlinear vertical maps governed by Rakotch contractions. The resulting object is a bounded, Riemann integrable, generally noncontinuous fractal function whose closed graph is the attractor of the underlying IFS and whose area-matching property is enforced by explicit conditions on the coefficients.

## 1. Histopolation and the underlying data model

The basic distinction is between interpolation and histopolation. Interpolation requires pointwise fitting,
\[
f(t_j)=y_j,\qquad j=0,1,\dots,N,
\]
whereas histopolation requires area matching over prescribed subintervals. For
\[
I=[t_0,t_N],\qquad t_0<t_1<\cdots<t_N,
\]
with
\[
I_j=[t_{j-1},t_j)\quad (j=1,\dots,N-1),\qquad I_N=[t_{N-1},t_N],
\]
and histogram data
\[
F=\{y_1,y_2,\dots,y_N\},
\]
a histopolating function satisfies
\[
\int_{t_{j-1}}^{t_j} f(t)\,dt = y_j\,a_j\,(t_N-t_0),
\]
where the normalized lengths are defined by
\[
|I_j|=a_j(t_N-t_0).
\]
Accordingly, histopolation is an area-preserving analogue of interpolation rather than a point-matching one [2509.18657].

This formulation is technically significant because it removes the requirement that the approximant be continuous or interpolatory at knot points. The 2025 framework treats histogram data directly, rather than recasting the problem as interpolation of cumulative quantities. That distinction is central in the fractal setting, where the relevant fixed points may be noncontinuous and where join-up conditions are not part of the definition of the histopolation problem.

## 2. Self-referential construction

The construction begins with affine domain contractions
\[
l_j:I\to I,\qquad l_j(t)=a_j t+b_j,
\]
chosen so that
\[
l_j(t_0)=t_{j-1},\qquad l_j(t_N)=t_j.
\]
Each \(l_j\) therefore maps the full interval \(I\) onto the subinterval \(I_j\), and the family \(\{I_j\}\) forms a disjoint partition of \(I\) [2509.18657].

The vertical maps are defined on
\[
A=I\times\mathbb{R}
\]
by
\[
F_j:A\to\mathbb{R},\qquad F_j(t,x)=c_j t+\delta_j s_j(x)+d_j,
\]
where \(c_j,d_j,\delta_j\) are constants and \(s_j\) are nonlinear Rakotch contractions. A remark in the same work also allows variable vertical scaling functions \(\delta_j(t)\) under corresponding bounds. The associated IFS maps are
\[
w_j:A\to A,\qquad w_j(t,x)=(l_j(t),F_j(t,x)),
\]
and the IFS is
\[
\mathcal I=\{A;w_j,\ j\in\mathbb N_N\}.
\]

The function itself is obtained through the Read–Bajraktarević operator on the Banach space \(\mathcal B(I)\) of bounded real-valued functions on \(I\), equipped with
\[
\|h\|_\infty=\sup_{t\in I}|h(t)|.
\]
For \(h\in\mathcal B(I)\),
\[
Th(t)=F_j(l_j^{-1}(t),h(l_j^{-1}(t)))
      =c_j l_j^{-1}(t)+\delta_j s_j(h(l_j^{-1}(t)))+d_j,\qquad t\in I_j.
\]
Its unique fixed point \(f\in\mathcal B(I)\) satisfies the self-referential equation
\[
f(t)=c_jl_j^{-1}(t)+\delta_js_j(f(l_j^{-1}(t)))+d_j,\qquad t\in I_j,\ j\in\mathbb N_N.
\]
This fixed point is the nonlinear fractal function [2509.18657].

## 3. Rakotch contractions and fixed-point theory

A central structural feature of the theory is the replacement of Banach contractions by Rakotch contractions. A map \(f\) is a \(\psi\)-contraction if
\[
\rho(f(t),f(s))\le \psi(\rho(t,s)).
\]
It is a Rakotch contraction if, in addition,
\[
\frac{\psi(t)}{t}<1,\qquad \frac{\psi(t)}{t}\le \frac{\psi(s)}{s}\quad \forall\, t>s.
\]
Every Banach contraction is Rakotch, but not conversely. The fixed-point theorem used here states that a Rakotch contraction on a complete metric space has a unique fixed point, and the iterates converge to it [2509.18657].

For the IFS on \(A=I\times\mathbb R\), the metric is
\[
d_\eta((t,x),(t',x'))=|t-t'|+\eta|x-x'|,
\]
with
\[
\eta=\frac{1-\max_{j\in\mathbb N_N}|a_j|}{2(C+1)},\qquad C=\max_{j\in\mathbb N_N}|c_j|.
\]
Under the stated assumptions, each \(w_j\) is a Rakotch contraction, and the IFS therefore has a unique attractor [2509.18657].

The same contractive mechanism governs the Read–Bajraktarević operator. The key estimate is
\[
d_{\mathcal B(I)}(Tg,Tg')\le \delta\,\psi(d_{\mathcal B(I)}(g,g')), \qquad \delta=\max_j|\delta_j|.
\]
Hence \(T\) is a Rakotch contraction on \(\mathcal B(I)\), so it has a unique fixed point. The shift from Banach to Rakotch contractions enlarges the admissible class of nonlinearities and is one of the principal theoretical extensions of the construction.

## 4. Boundedness, graph attractor, and integrability

The fixed point \(f\) is bounded because \(T\) maps bounded functions to bounded functions. Its graph,
\[
G=\{(t,f(t)):t\in I\},
\]
need not be closed, because continuity is not assumed. The closure of the graph is nevertheless characterized exactly by the IFS attractor:
\[
\overline G = B,\qquad B=\bigcup_{j=1}^N w_j(B).
\]
The proof proceeds by establishing both
\[
\overline G\subseteq W(\overline G)
\quad\text{and}\quad
W(\overline G)\subseteq \overline G,
\]
where
\[
W(K)=\bigcup_{j=1}^N w_j(K).
\]
Thus the closed graph attractor is the geometric realization of the nonlinear fractal histopolation function [2509.18657].

Riemann integrability is obtained under the conditions
\[
\alpha=\max_j|a_j|<1,\qquad \delta=\max_j|\delta_j|<1.
\]
The proof uses oscillation estimates and Lebesgue’s criterion. The key estimate is
\[
w_f(l_j(I'))\le C|I'|+\delta w_f(I'),
\]
where \(C=\max_j|c_j|\). Iteration along cylinder sets yields bounds tending to zero as the depth increases, so the oscillation vanishes at almost every point. Equivalently, the discontinuities lie in a countable set \(D\), hence in a measure-zero set. The function is therefore Riemann integrable [2509.18657].

This part of the theory is important because histopolation imposes integral constraints rather than point constraints. The construction does not merely produce a bounded self-referential object; it produces one for which the relevant subinterval integrals are well defined in the Riemann sense.

## 5. Exact histopolation criterion and admissible vertical scaling

The decisive characterization of the area-matching property is the histopolation criterion. The fixed point \(f\) solves the histopolation problem if and only if
\[
d_j=\frac{2y_j(t_N-t_0)-c_j(t_N^2-t_0^2)-2\delta_j\int_I s_j(f(t))\,dt}{2(t_N-t_0)}.
\]
This formula identifies the constants \(d_j\) as the parameters that enforce the histogram-area constraints. The self-referential structure is therefore not independent of the area data: the function is histopolating precisely when the vertical offsets satisfy this integral relation [2509.18657].

A notable novelty is the controlled allowance of vertical scaling factors greater than one. Two cases are distinguished. If
\[
\sup_{t>0}\frac{\psi(t)}{t}<1,
\]
one chooses \(\beta\) such that
\[
\sup_{t>0}\frac{\psi(t)}{t}<\beta<1
\]
and then requires
\[
\delta=\max_j|\delta_j|<\frac{1}{\beta}.
\]
If instead
\[
\sup_{t>0}\frac{\psi(t)}{t}=1,
\]
one requires
\[
\delta=\max_j|\delta_j|<1.
\]
In the first case, \(\delta>1\) is permitted while still ensuring that the composite contraction factor
\[
\zeta(t)=\delta\psi(t)
\]
satisfies
\[
\frac{\zeta(t)}{t}<1.
\]
A corresponding remark allows variable scaling functions \(\delta_j(t)\), provided
\[
\delta=\max_j\sup_{t\in I}|\delta_j(t)|
\]
satisfies the same constraints [2509.18657].

These admissibility conditions distinguish the nonlinear Rakotch-based formulation from many classical fractal interpolation and histopolation schemes, where vertical scales are typically restricted by Banach-type bounds.

## 6. Relation to earlier fractal histopolation and representative examples

Classical fractal interpolation functions are continuous, interpolate data points, typically require join-up conditions, and often use Banach contractions. Earlier fractal histopolation, including the Barnsley–Viswanathan line of work, solves area constraints, is not necessarily continuous, and was developed with Banach contractions. The nonlinear formulation extends that setting by constructing a nonlinear fractal histopolation function, using Rakotch contractions, allowing vertical scaling factors greater than one under precise conditions, and providing a fixed-point characterization together with an explicit integral criterion for histopolation [2509.18657].

The 2015 work "Discontinuous Fractal Functions and Fractal Histopolation" formalized a fractal histopolant as an integrable fractal function satisfying
\[
\int_{x_{i-1}}^{x_i} f(x)\,dx = h_i f_i,\qquad i=1,2,\dots,N,
\]
for a partition \(\Delta=\{x_0<\cdots<x_N\}\), with affine maps
\[
L_i(x)=a_i x+b_i,\qquad F_i(x,y)=\alpha_i y+q_i(x),
\]
and \(|\alpha_i|<1\). It also derived the exact area constraint
\[
\int_I q_i(x)\,dx = \frac{h_i f_i-\alpha_i a_i\sum_{i=1}^N h_i f_i}{a_i},
\]
and established that the discontinuities of the bounded fractal function form a Lebesgue null set, so the function is Riemann integrable [1503.06903]. Relative to that framework, the nonlinear Rakotch-based construction retains the area-preserving objective while enlarging the allowable self-referential dynamics.

The 2022 work "Fractal Interpolation over Nonlinear Partitions" introduced a nonlinear partition setting
\[
X=\coprod_{i=1}^n h_i(X),
\]
with nonlinear diffeomorphisms \(h_i\), and the self-referential equation
\[
\psi(h_i(x))=q_i(x)+s_i(x)\psi(x).
\]
It proved existence and uniqueness in \(B(X,F)\), and also developed \(L^p(X,F)\) and \(C^\alpha(X,F)\) theories under explicit contractivity conditions, but it did not formulate histopolation constraints [2203.07935]. This suggests a broader route toward histopolation over nonlinear partitions, although the 2025 construction itself is formulated on affine interval partitions.

Two examples illustrate the scope of the nonlinear theory. For \(I=[0,1]\) with partition \(\{0,\tfrac12,1\}\), one example takes
\[
l_1(t)=\frac12 t,\qquad l_2(t)=\frac12 t+\frac12,
\]
with
\[
s_1(x)=s_2(x)=\frac12\sin x,
\]
and variable scaling factors
\[
\delta_1(t)=\frac32 t,\qquad \delta_2(t)=\frac74 t.
\]
The resulting fixed point satisfies a piecewise self-referential equation and illustrates that the framework accommodates variable vertical scaling exceeding \(1\) [2509.18657]. A second example, again on the partition \(\{0,\tfrac12,1\}\) but with histogram \(\{5,6\}\), uses
\[
F_1(t,x)=\frac12 t+\frac12\frac{1}{1+x}+d_1,\qquad
F_2(t,x)=\frac14 t+\frac14\frac{1}{1+x}+d_2,
\]
and the fact that \(\frac{1}{1+x}\) is a Rakotch contraction. When \(d_1,d_2\) satisfy the histopolation condition above, the fixed point is a solution to the histopolation problem [2509.18657].

Taken together, these developments place the nonlinear fractal histopolation function at the intersection of fractal interpolation, discontinuous fractal approximation, and area-preserving approximation. Its defining feature is not interpolation of nodal values but exact preservation of prescribed histogram areas through a nonlinear self-referential fixed-point mechanism.

Source: https://www.emergentmind.com/topics/nonlinear-fractal-histopolation-function