Decomposition of multivariate function using the Heaviside step function

Whereas the Dirac delta function introduced by P. A. M. Dirac in 1930 to develop his theory of quantum mechanics has been well studied, a not famous formula related to the delta function using the Heaviside step function in a single-variable form, also given by Dirac, has been poorly studied. Following Dirac’s method, we demonstrate the decomposition of a multivariate function into a sum of integrals in which each integrand is composed of a derivative of the function and a direct product of Heaviside step functions. It is an extension of Dirac’s single-variable form to that for multiple variables.

1. Introduction P. A. M. Dirac introduced in 1930 a function, now called the Dirac delta function, to develop his theory of quantum mechanics (Dirac 1958). It takes value zero at x ≠ 0 and its integral is unity. A fundamental property derivable from the definition of the Dirac delta function is that any multivariate real function can be expressed with delta functions δ and integrals as follows, The importance of this property is analogous to the Fourier transform (Bracewell 1965) for its ability to yield an alternative representation of any multivariate function in which the variables of the function are changed. The delta function can be seen in applications from physics to engineering: such as quantum mechanical states (Lee 1992); quantum similarity integrals (Safouhi and Berlu 2006); pseudopotential (Derevianko 2003); a spin system with a classical environment (Calvani et al. 2013); and generally, numbers of formulae in the Fourier and Laplace transforms, and differential equations (Schwartz 1966;Kreyszig 2011). A more rigorous mathematical theory for the delta function has also been developed and expanded under the branch in pure mathematics called the theory of distributions by L. Schwartz (Schwartz 1945). Further development is the generalized delta impulse (Corinthios 2003), which is an extension of the Dirac delta function to that on the complex plane and is applied to theories of generalized Laplace, z, Hilbert, and Fourier-related transforms (Corinthios 2005(Corinthios , 2007.
Transforming the integral expression in one variable using the Dirac delta function δ into one using the Heaviside step function σ, is essentially described in Dirac's quantum mechanics text (Dirac 1958). It is derived from the relation between the Dirac delta function and the derivative of the Heaviside step function. As these two expressions, the left-and right-hand sides of (1.2), are mathematically equivalent, the step-function expression would be expected to find potent applications in physics and engineering. One general example must be approximation theory (Milovanovic and Rassias 2014). The other example is replacement of the Dirac delta functions with the Heaviside step functions by which (divergent) delta functions can be hidden from integrand. It helps a rigorous formalism by using only bounded functions without advanced Schwartz distributions. Nevertheless, the step-function expression has not been extended so far compared with the application of the delta-function expression. Here we demonstrate a unified formula that extends this step-function expression for singlevariable functions to multiple-variable functions. It can be interpreted as the decomposition of any multivariate function with respect to the Heaviside step function.
2. Decomposition of multivariate functions using the Heaviside step function

Definition
Let R(X 1 , X 2 , ⋅ ⋅⋅, X N ) be a continuous real function defined for 0 ≤ X i < ∞ and satisfies: Whose derivatives exist and continuous where α = α 1 + ⋯ + α N and α i ≥ 0 is a natural number.
can be integrated with respect to a given X i while the other variables are held fixed. Sequences of functions 1 for X j and j ≠ i.

Premise
We follow Dirac's method: The Dirac delta function is regarded as a function (not a Schwartz distribution) The derivative of the Heaviside step function is regarded as the Dirac delta function

Theory
We demonstrate that the defined function R(X 1 , X 2 , ⋯, X N ) can be decomposed into: defines a set of Heaviside step functions for each and α i = 0 or 1.

Proof
The formula (2.1) may be expressed as: We use mathematical induction.
Using the definition of the Dirac delta-function δ, then for any real function R(X 1 , The expression (2.2) for single-variable functions (N =1), ð2:6Þ holds for some N. Using (2.6), the following, holds because R(X 1 , ⋯, X N , μ N + 1 ) can be regarded as one of the R(X 1 , ⋯, X N ) appending a parameter μ N + 1 . Multiplying both sides of (2.7) by δ(μ N + 1 − X N + 1 ) and then integrating each term with respect to μ N + 1 , one obtains on the left-hand side, ð2:8Þ and on the right-hand side,

ð2:9Þ
The order of the integrations can be changed because each integrand can be integrated with respect to its corresponding μ i while holding other variables fixed. Therefore, (2.9) can be transformed into The terms enclosed in braces in (2.10) can be transformed using Lemma 2.1 as follows, ð2:11Þ where Lemma 2.2 was also used.

Finally (2.11) becomes
ð2:12Þ Thus, assuming expression (2.6) for N leads to the same expression for N +1. The formula for N =1 also holds as described above. Therefore, (2.1) holds for any natural number N.
The key is transforming integrands with Dirac delta functions to ones with Heaviside step functions. It is represented by Lemma 2.1, which is intuitively understandable with equality between sums of vertical stripes and sums of horizontal stripes under the target integrand F (Figure 1).

Lemmas
where F is R or its derivatives.