Bounded variation

From Wikipedia, the free encyclopedia

In mathematical analysis, a function of bounded variation refers to a real-valued function whose total variation is bounded (finite): the graph of a function having this property is well behaved in a precise sense. For a continuous function of a single variable, being of bounded variation means that the distance along the direction of y-axis (i.e. the distance calculated neglecting the contribution of motion along x-axis) traveled by an ideal point moving along the graph of the given function (which, under given hypothesis, is also a continuous path) has a finite value. For a continuous function of several variables, the meaning of the definition is the same, except for the fact that the continuous path to be considered cannot be the whole graph of the given function (which is an hypersurface in this case), but can be every intersection of the graph itself with a plane parallel to a fixed x-axis and to the y-axis.

Functions of bounded variation are precisely those with respect to which one may find Riemann-Stieltjes integrals of all continuous functions.

Another characterization states that the functions of bounded variation on a closed interval are exactly those f which can be written as a difference gh, where both g and h are bounded monotone.

In the case of several variables, a function f defined on an open subset Ω of \scriptstyle\mathbb{R}^n is said to have bounded variation if its distributional derivative is a finite vector Radon measure.

One of the most important aspects of functions of bounded variation is that they form a algebra of discontinuous functions whose first derivative exists almost everywhere: due to this fact, they can and frequently are used to define generalized solutions of nonlinear problems involving functionals, ordinary and partial differential equations in mathematics, physics and engineering. Considering the problem of multiplication of distributions or more generally the problem of defining general nonlinear operations on generalized functions, function of bounded variation are the smallest algebra which has to be embedded in every space of generalized functions preserving the result of multiplication.

Contents

[edit] History

According to Golubov, BV functions of a single variable were first introduced by Camille Jordan, in the paper (Jordan 1881) dealing with the convergence of Fourier series. The first step in the generalization of this concept to functions of several variables was due to Leonida Tonelli, who introduced a class of continuous BV functions in 1926 (Cesari 1986, pp. 47-48), to extend his direct method for finding solutions to problems in the calculus of variations in more than one variable. Ten years after, in 1936, Lamberto Cesari changed the continuity requirement in Tonelli's definition to a less restrictive integrability requirement, obtaining for the first time the class of functions of bounded variation of several variables in its full generality: as Jordan did before him, he applied the concept to resolve of a problem concerning the convergence of Fourier series, but for functions of two variables. After him, several authors applied BV functions to study Fourier series in several variables, geometric measure theory, calculus of variations, and mathematical physics: Renato Caccioppoli and Ennio de Giorgi used them to define measure of non smooth boundaries of sets (see voice "Caccioppoli set" for further informations), Edward D. Conway and Joel A. Smoller applied them to the study of a single nonlinear hyperbolic partial differential equation of first order in the paper (Conway & Smoller 1966), proving that the solution of the Cauchy problem for such equations is a function of bounded variation, provided the initial value belongs to the same class. Aizik Isaakovich Vol'pert developed extensively a calculus for BV functions in the paper (Vol'pert 1967) and in the book (Hudjaev & Vol'pert 1986).

[edit] Formal definition

[edit] BV functions of one variable

Definition 1. The total variation of a real-valued function f, defined on a interval  \scriptstyle [a , b] \subset \mathbb{R} is the quantity

 V^a_b(f)=\sup_{P \in \mathcal{P}} \sum_{i=0}^{n_P-1} | f(x_{i+1})-f(x_i) |. \,

where the supremum is taken over the set  \scriptstyle \mathcal{P} =\left\{P=\{ x_0, \dots , x_{n_p}\}|P\text{ is a partition of } [a,b] \right\} of all partitions of the interval considered.

If f is differentiable and its derivative is integrable, its total variation is the vertical component of the arc-length of its graph, that is to say,

 V^a_b(f) = \int _a^b |f'(x)|\, dx.

Definition 2. A real-valued function f on the real line is said to be of bounded variation (BV function) on a chosen interval [a,b] if its total variation is finite, i.e.

 f \in BV([a,b]) \iff V^a_b(f) < +\infty

It can be proved that a real function f is of bounded variation in an interval if and only if it can be written as the difference f = f1f2 of two non decreasing functions (This is known as the Jordan decomposition.)

Through the Stieltjes integral, any function of bounded variation on a closed interval [a,b] defines a bounded linear functional on C([a,b]). In this special case (Kolmogorov & Fomin 1969, pp. 374-376), the Riesz representation theorem states that every bounded linear functional arises uniquely in this way. The normalised positive functionals or probability measures correspond to positive non-decreasing lower semicontinuous functions. This point of view has been important in spectral theory (Riesz & Sz.-Nagy 1990), in particular in its application to ordinary differential equations.

[edit] BV functions of several variables

Functions of bounded variation, BV functions, are functions whose distributional derivative is a finite Radon measure. More precisely:

Definition 1 Let Ω be an open subset of  \scriptstyle\mathbb{R}^n . A locally integrable function u is said of bounded variation (BV function), and write

 u\in BV(\Omega)

if there exists a finite vector Radon measure  \scriptstyle Du\in\mathcal M(\Omega,\mathbb{R}^n) such that the following equality holds

 
\int_\Omega u(x)\,\mathrm{div}\boldsymbol{\phi}(x)\, dx = - \int_\Omega \langle\boldsymbol{\phi}, Du(x)\rangle 
\qquad \forall\boldsymbol{\phi}\in C_c^1(\Omega,\mathbb{R}^n)

that is, u defines a linear functional on the space  \scriptstyle C_c^1(\Omega,\mathbb{R}^n) of continuously differentiable vector functions  \scriptstyle\boldsymbol{\phi} of compact support contained in Ω: the vector measure Du represents therefore the distributional or weak gradient of u.

An equivalent definition is the following.

Definition 2 Given a locally integrable function u, the total variation of u in is defined as

 V(u,\Omega):=\sup\left\{\int_\Omega u\mathrm{div}\boldsymbol{\phi}\colon \phi\in C_c^1(\Omega,\mathbb{R}^n),\ \Vert\boldsymbol{\phi}\Vert_{L^\infty(\Omega)}\le 1\right\}.

where  \scriptstyle \Vert\;\Vert_{L^\infty(\Omega)} is the essential supremum norm.

The space of functions of bounded variation (BV functions) can then be defined as

 BV(\Omega)=\{ u\in L^1(\Omega)\colon V(u,\Omega)<+\infty\}

The two definition are equivalent since if  \scriptstyle V(u,\Omega)<+\infty then

\left|\int_\Omega u(x)\,\mathrm{div}\boldsymbol{\phi}\, dx \right |\leq V(u,\Omega)\Vert\boldsymbol{\phi}\Vert_{L^\infty(\Omega)}
\qquad \forall \boldsymbol{\phi}\in C_c^1(\Omega,\mathbb{R}^n)

therefore \scriptstyle \int_\Omega u(x)\,\mathrm{div}\boldsymbol{\phi}(x) defines a continuous linear functional on the space \scriptstyle C_c^1(\Omega,\mathbb{R}^n). Since \scriptstyle C_c^1(\Omega,\mathbb{R}^n)
\subset C^0(\Omega,\mathbb{R}^n) as a linear subspace, this continuous linear functional can be extended continuously and linearily to the whole \scriptstyle C^0(\Omega,\mathbb{R}^n) by the Hahn–Banach theorem i.e. it defines a Radon measure.

[edit] Basic properties

Only the properties common to functions of one variable and to functions of several variables will be considered in the following, and proofs will be carried on only for functions of several variables since the the proof for the case of one variable is a straightforward adaptation of the several variables case. References (Giusti 1984, pp. 7-9), (Hudjaev & Vol'pert 1986) and (Màlek et al. 1996) are extensively used.

[edit] BV functions have only jump-type singularities

In the case of one variable, the assertion is clear: for each point x0 in the interval  \scriptstyle ]a , b[ \subset \mathbb{R} of definition of the function u, either one of the following two assertions is true

 \lim_{x\rightarrow x_{0^-}}\!\!\!u(x) = \!\!\!\lim_{x\rightarrow x_{0^+}}\!\!\!u(x)
 \lim_{x\rightarrow x_{0^-}}\!\!\!u(x) \neq \!\!\!\lim_{x\rightarrow x_{0^+}}\!\!\!u(x)

while both limits exists and are finite. In the case of functions several variables, there are some premises to understand: first of all, there is a continuum of directions along with is possible to approach a given point x0 belonging to the domain \scriptstyle\Omega\in\mathbb{R}^n. It is necessary to precise a suitable concept of limit: chosing a unit vector \scriptstyle{\boldsymbol\hat a}\in\mathbb{R}^n it is possible to divide Ω in two sets

\Omega_{({\boldsymbol\hat a},\boldsymbol{x}_0)} = \Omega \cap \{\boldsymbol{x}\in\mathbb{R}^n|\langle\boldsymbol{x}-\boldsymbol{x}_0,{\boldsymbol\hat a}\rangle>0\} \qquad 
\Omega_{(-{\boldsymbol\hat a},\boldsymbol{x}_0)} = \Omega \cap \{\boldsymbol{x}\in\mathbb{R}^n|\langle\boldsymbol{x}-\boldsymbol{x}_0,-{\boldsymbol\hat a}\rangle>0\}

Then for each point x0 belonging to the domain \scriptstyle\Omega\in\mathbb{R}^n of the BV function u or one of the following two assertion is true

 \lim_{\overset{\boldsymbol{x}\rightarrow \boldsymbol{x}_0}{\boldsymbol{x}\in\Omega_{({\boldsymbol\hat a},\boldsymbol{x}_0)}}}\!\!\!\!\!\!u(\boldsymbol{x}) = \!\!\!\!\!\!\!\lim_{\overset{\boldsymbol{x}\rightarrow \boldsymbol{x}_0}{\boldsymbol{x}\in\Omega_{(-{\boldsymbol\hat a},\boldsymbol{x}_0)}}}\!\!\!\!\!\!\!u(\boldsymbol{x})
  \lim_{\overset{\boldsymbol{x}\rightarrow \boldsymbol{x}_0}{\boldsymbol{x}\in\Omega_{({\boldsymbol\hat a},\boldsymbol{x}_0)}}}\!\!\!\!\!\!u(\boldsymbol{x}) \neq \!\!\!\!\!\!\!\lim_{\overset{\boldsymbol{x}\rightarrow \boldsymbol{x}_0}{\boldsymbol{x}\in\Omega_{(-{\boldsymbol\hat a},\boldsymbol{x}_0)}}}\!\!\!\!\!\!\!u(\boldsymbol{x})

or x0 belongs to a subset of Ω having zero n − 1-dimensional Hausdorff measure. The quantities

\lim_{\overset{\boldsymbol{x}\rightarrow \boldsymbol{x}_0}{\boldsymbol{x}\in\Omega_{({\boldsymbol\hat a},\boldsymbol{x}_0)}}}\!\!\!\!\!\!u(\boldsymbol{x})=u_{\boldsymbol{\hat a}}(\boldsymbol{x}) \qquad \lim_{\overset{\boldsymbol{x}\rightarrow \boldsymbol{x}_0}{\boldsymbol{x}\in\Omega_{(-{\boldsymbol\hat a},\boldsymbol{x}_0)}}}\!\!\!\!\!\!\!u(\boldsymbol{x})=u_{-\boldsymbol{\hat a}}(\boldsymbol{x})

are called approximate limits of the BV function u at the point x0.

[edit] V(·, Ω) is lower semi-continuous on BV(Ω)

The functional \scriptstyle V(\cdot,\Omega):BV(\Omega)\rightarrow \mathbb{R}^+ is lower semi-continuous: to see this, choose a Cauchy sequence of BV-functions \scriptstyle\{u_n\}_{n\in\mathbb{N}} converging to \scriptstyle u\in L^1_{loc}(\Omega). Then, since all the functions of the sequence and their limit function are integrable and by the definition of lower limit

\liminf_{n\rightarrow\infty}V(u_n,\Omega)\geq\lim_{n\rightarrow\infty}\int_\Omega u_n(x)\,\mathrm{div}\boldsymbol{\phi}\, dx = \int_\Omega u(x)\,\mathrm{div}\boldsymbol{\phi}\, dx \qquad\forall\boldsymbol{\phi}\in C_c^1(\Omega,\mathbb{R}^n),\quad\Vert\boldsymbol{\phi}\Vert_{L^\infty(\Omega)}\leq 1

Now considering the supremum on the set of functions \scriptstyle\boldsymbol{\phi}\in C_c^1(\Omega,\mathbb{R}^n) such that \scriptstyle \Vert\boldsymbol{\phi}\Vert_{L^\infty(\Omega)}\leq 1 then the following inequality holds true

\liminf_{n\rightarrow\infty}V(u_n,\Omega)\geq V(u,\Omega)

which is exactly the definition of lower semicontinuity.

[edit] BV(Ω) is a Banach space

By defintion BV(Ω) is a subset of L1loc(Ω), while linearity follows from the linearity properties of the defining integral i.e.

\begin{align}
\int_\Omega [u(x)+v(x)]\,\mathrm{div}\boldsymbol{\phi}(x) & =\int_\Omega 
                    u(x)\,\mathrm{div}\boldsymbol{\phi}(x)+\int_\Omega v(x)\,\mathrm{div}\boldsymbol{\phi}(x) = \\
                                             & =- \int_\Omega \langle\boldsymbol{\phi}(x), Du(x)\rangle- \int_\Omega \langle \boldsymbol{\phi(x)}, Dv(x)\rangle
                                               =- \int_\Omega \langle \boldsymbol{\phi}(x), [Du(x)+Dv(x)]\rangle 
\end{align}

for all \scriptstyle\phi\in C_c^1(\Omega,\mathbb{R}^n) therefore \scriptstyle u+v\in BV(\Omega)for all \scriptstyle u,v\in BV(\Omega), and

\int_\Omega cu(x)\,\mathrm{div}\boldsymbol{\phi}(x)=c\int_\Omega u(x)\,\mathrm{div}\boldsymbol{\phi}(x)=-c \int_\Omega \langle \boldsymbol{\phi}(x), Du(x)\rangle

for all \scriptstyle c\in\mathbb{R}, therefore \scriptstyle cu\in BV(\Omega) for all \scriptstyle u\in BV(\Omega), and all \scriptstyle c\in\mathbb{R}. The proved vector space properties imply that BV(Ω) is a vector subspace of L1(Ω). Consider now the function \scriptstyle\|\;\|_{BV}:BV(\Omega)\rightarrow\mathbb{R}^+ defined as

\| u \|_{BV} := \| u \|_{L^1} + V(u,\Omega)

where \scriptstyle\| \; \|_{L^1} is the usual L1(Ω) norm: it is easy to prove that this is a norm on BV(Ω). To see that BV(Ω) is complete respect to it, i.e. it is a Banach space, consider a Cauchy sequence \scriptstyle\{u_n\}_{n\in\mathbb{R}} in BV(Ω). By definition it is also a Cauchy sequence in L1(Ω) and therefore has a limit u in L1(Ω): since un is bounded in BV(Ω) for each n, then \scriptstyle \Vert u \Vert_{BV} < +\infty by lower semicontinuity of the variation \scriptstyle V(\cdot,\Omega), therefore u is a BV function. Finally, again by lower semicontinuity, choosing an arbitrary small positive number \scriptstyle\varepsilon

\Vert u_j - u_k \Vert_{BV}<\varepsilon\quad\forall j,k\geq N\in\mathbb{N} \quad\Rightarrow\quad V(u_k-u,\Omega)\leq \liminf_{j\rightarrow +\infty} V(u_k-u_j,\Omega)\leq\varepsilon

[edit] Chain rule for BV functions

Chain rules for nonsmooth functions are very important in mathematics and mathematical physics since there are several important physical models whose behavior is described by functions or functionals with a very limited degree of smoothness.The following version is proved in the paper (Vol'pert 1967, p. 248): all partial derivatives must be intended in a generalized sense. i.e. as generalized derivatives

Theorem. Let \scriptstyle f:\mathbb{R}^p\rightarrow\mathbb{R} be a function of class C1 (i.e. a continuous and differentiable function having continuous derivatives) and let \scriptstyle\boldsymbol{u}(\boldsymbol{x})=(u_1(\boldsymbol{x}),\ldots,u_p(\boldsymbol{x})) be a function in BV(Ω) with Ω being an open subset of  \scriptstyle\mathbb{R}^n . Then \scriptstyle f\circ\boldsymbol{u}(\boldsymbol{x})=f(\boldsymbol{u}(\boldsymbol{x}))\in BV(\Omega) and

\frac{\partial f(\boldsymbol{u}(\boldsymbol{x}))}{\partial x_i}=\sum_{k=1}^p\frac{\partial\bar{f}(\boldsymbol{u}(\boldsymbol{x}))}{\partial u_k}\frac{\partial{u_k(\boldsymbol{x})}}{\partial x_i}
\qquad\forall i=1,\ldots,n

where \scriptstyle\bar f(\boldsymbol{u}(\boldsymbol{x})) is the mean value of the function at the point \scriptstyle x \in\Omega, defined as

\bar f(\boldsymbol{u}(\boldsymbol{x}))=\int_0^1 f\left(\boldsymbol{u}_{\boldsymbol{\hat a}}(\boldsymbol{x})t + \boldsymbol{u}_{-\boldsymbol{\hat a}}(\boldsymbol{x})(1-t)\right)dt

A more general chain rule formula for Lipschitz continuous functions \scriptstyle f:\mathbb{R}^p\rightarrow\mathbb{R}^s has been found by Luigi Ambrosio and Gianni Dal Maso and published in the paper (Ambrosio & Dal Maso 1990). However, even this formula has very important direct consequences: choosing \scriptstyle f(u)=v(\boldsymbol{x})u(\boldsymbol{x}) where \scriptstyle v(\boldsymbol{x}) is a BV function the preceding formula becomes

\frac{\partial v(\boldsymbol{x})u(\boldsymbol{x})}{\partial x_i} = {\bar u(\boldsymbol{x})}\frac{\partial v(\boldsymbol{x})}{\partial x_i} + 
{\bar v(\boldsymbol{x})}\frac{\partial u(\boldsymbol{x})}{\partial x_i}

This implies that the product of two functions of bounded variation is again a function of bounded variation, therefore BV(Ω) is an algebra.

[edit] Generalizations and extensions

[edit] Weighted BV functions

It is possible to generalize the above notion of total variation so that different variations are weighted differently. More precisely, let \scriptstyle \varphi : [0, +\infty)\longrightarrow [0, +\infty) be any increasing function such that \scriptstyle \varphi(0) = \varphi(0+) =\lim_{x\rightarrow 0_+}\varphi(x) = 0 (the weight function) and let \scriptstyle f : [0, T]\longrightarrow X be a function from the interval  \scriptstyle [0 , T] \subset \mathbb{R} taking values in a normed vector space X. Then the \scriptstyle \boldsymbol\varphi-variation of f over [0,T] is defined as

\mathop{\varphi\mbox{-Var}}_{[0, T]} (f) := \sup \sum_{j = 0}^{k} \varphi \left( | f(t_{j + 1}) - f(t_{j}) |_{X} \right),

where, as usual, the supremum is taken over all finite partitions of the interval [0,T], i.e. all the finite sets of real numbers ti such that

0 = t_{0} < t_{1} < \ldots < t_{k} = T.

The original notion of variation considered above is the special case of \scriptstyle \varphi-variation for which the weight function is the identity function: therefore a integrable function f is said to be a weighted BV function (of weight \scriptstyle\varphi) if and only if its \scriptstyle \varphi-variation is finite.

f\in BV_\varphi([0, T];X)\iff \mathop{\varphi\mbox{-Var}}_{[0, T]} (f) <+\infty

The space \scriptstyle BV_\varphi([0, T];X) is a topological vector space with respect to the norm

\| f \|_{BV_\varphi} := \| f \|_{\infty} + \mathop{\varphi \mbox{-Var}}_{[0, T]} (f),

where \scriptstyle\| f \|_{\infty} denotes the usual supremum norm of f. Weighted BV functions were introduced and studied in full generality by Wladislav Orlicz and Julian Musielak in the paper (Musielak & Orlicz 1959): Laurence Chisholm Young studied early the case \scriptstyle\varphi(x)=x^p where p is a positive integer.

[edit] SBV functions

SBV functions i.e. Special functions of Bounded Variation where introduced by Luigi Ambrosio and Ennio de Giorgi in the paper (Ambrosio & De Giorgi 1988), dealing with free discontinuity variational problems: given a open subset Ω of  \scriptstyle\mathbb{R}^n , the space SBV(Ω) is a proper subspace of BV(Ω), since the weak gradient of each function belonging to it const exatcly of the sum of a n-dimensional support and a n − 1-dimensional support measure and no lower-dimensional terms, as seen in the following definition.

Definition. Given a locally integrable function u, then \scriptstyle u\in {S\!BV}(\Omega) if and only if

1. There exist two Borel functions f and g of domain Ω and codomain \scriptstyle \mathbb{R}^n such that

 \int_\Omega\vert f\vert dH^n+ \int_\Omega\vert g\vert dH^{n-1}<+\infty.

2. For all of continuously differentiable vector functions  \scriptstyle\phi of compact support contained in Ω, i.e. for all  \scriptstyle \phi \in 
C_c^1(\Omega,\mathbb{R}^n) the following formula is true:

 \int_\Omega u\mbox{div} \phi dH^n = \int_\Omega \langle \phi, f\rangle dH^n +\int_\Omega \langle \phi, g\rangle dH^{n-1}.

where Hα is the α-dimensional Hausdorff measure.

Details on the properties of SBV functions can be found in works cited in the bibliography section: particularly the paper (De Giorgi 1992) contains a useful bibliography.

[edit] bv sequences

As particular examples of Banach spaces, Dunford & Schwartz (1958, Chapter IV) consider spaces of sequences of bounded variation, in addition to the spaces of functions of bounded variation. The total variation of a sequence x=(xi) of real or complex numbers is defined by

TV(x) = \sum_{i=1}^\infty |x_{i+1}-x_i|.

The space of all sequences of finite total variation is denoted by bv. The norm on bv is given by

\|x\|_{bv} = |x_1| + TV(x) = |x_1| +  \sum_{i=1}^\infty |x_{i+1}-x_i|.

With this norm, the space bv is a Banach space.

The total variation itself defines a norm on a certain subspace of bv, denoted by bv0, consisting of sequences x = (xi) for which

\lim_{n\to\infty} x_n =0.

The norm on bv0 is denoted

\|x\|_{bv_0} = TV(x) = \sum_{i=1}^\infty |x_{i+1}-x_i|.

With respect to this norm bv0 becomes a Banach space as well.

[edit] Examples

The function f(x)=sin(1/x) is not of bounded variation.
The function f(x)=sin(1/x) is not of bounded variation.

The function

f(x) = \begin{cases} 0, & \mbox{if }x =0 \\ \sin(1/x), & \mbox{if } x \neq 0 \end{cases}

is not of bounded variation on the interval [0,2 / π]

The function f(x)=x sin(1/x) is not of bounded variation.
The function f(x)=x sin(1/x) is not of bounded variation.

While it is harder to see, the function

f(x) = \begin{cases} 0, & \mbox{if }x =0 \\ x \sin(1/x), & \mbox{if } x \neq 0 \end{cases}

is not of bounded variation on the interval [0,2 / π] either.

The function f(x)=x2 sin(1/x) is of bounded variation.
The function f(x)=x2 sin(1/x) is of bounded variation.

At the same time, the function

f(x) = \begin{cases} 0, & \mbox{if }x =0 \\ x^2 \sin(1/x), & \mbox{if } x \neq 0 \end{cases}

is of bounded variation on the interval [0,2 / π].

The Sobolev space W1,1(Ω) is a proper subset of BV(Ω). In fact, for each u in W1,1(Ω) it is possible to choose a measure  \scriptstyle \mu:=\nabla u \mathcal L (where  \scriptstyle\mathcal L is the Lebesgue measure on Ω) such that the equality

 \int u\mathrm{div}\phi = -\int \phi\, d\mu = -\int \phi \nabla u \qquad \forall \phi\in C_c^1

holds, since it is nothing more than the definition of weak derivative, and hence holds true. One can easily find an example of a BV function which is not W1,1.

[edit] Applications

[edit] Mathematics

Functions of bounded variation have been studied in connection with the set of discontinuities of functions and differentiability of real functions, and the following results are well-known. If f is a real function of bounded variation on an interval [a, b] then

[edit] Physics and engineering

The ability of BV functions to deal with discontinuities has made their use widespread in the applied sciences: solutions of problems in mechanics, physics, chemical kinetics are very often representable by functions of bounded variation. The book (Hudjaev & Vol'pert 1986) details a very ample set of mathematical physics applications of BV functions. Also there is some modern application which deserves a brief description.

  • The Mumford-Shah functional: the segmentation problem for a two-dimensional image, i.e. the problem of faithful reproduction of contours and grey scales is equivalent to the minimization of such functional.

[edit] See also

[edit] References

[edit] Bibliography

[edit] External links

[edit] Theory

[edit] Other


This article incorporates material from BV function on PlanetMath, which is licensed under the GFDL.