Relationships between generalized Wiener integrals and conditional analytic Feynman integrals over continuous paths

Let C[0, t] denote a generalized Wiener space, the space of real-valued continuous functions on the interval [0, t], and define a random vector Z _n: C[0, t] → R ⁿ⁺¹ by \({Z_n}\left( x \right) = \left( {x\left( 0 \right) + a\left( 0 \right),\int_o^{{t_1}} {h\left( s \right)dx\left( s \right) + x\left( 0 \right) + a\left( {{t_1}} \right),...,\int_0^{{t_n}} {h\left( s \right)dx\left( s \right) + x\left( 0 \right) + a\left( {{t_n}} \right)} } } \right)\), where a ∈ C[0, t], h ∈ L ₂[0, t], and 0 < t ₁ <... < t _n ≤ t is a partition of [0, t]. Using simple formulas for generalized conditional Wiener integrals, given Z _n we will evaluate the generalized analytic conditional Wiener and Feynman integrals of the functions F in a Banach algebra which corresponds to Cameron-Storvick’s Banach algebra S. Finally, we express the generalized analytic conditional Feynman integral of F as a limit of the non-conditional generalized Wiener integral of a polygonal function using a change of scale transformation for which a normal density is the kernel. This result extends the existing change of scale formulas on the classical Wiener space, abstract Wiener space and the analogue of the Wiener space C[0, t].     

A Neumann problem of Ambrosetti–Prodi type

The aim of this paper is to establish an Ambrosetti–Proditype result for the problem

$$\left\{ \begin{array}{ll}-\Delta{u} = g(x, u,\nabla{u}) + t\varphi \quad {\rm in}\, \Omega,\\ \frac{\partial{u}}{\partial\eta} = 0 \qquad\qquad\qquad\quad {\rm on}\, \partial\Omega ;\end{array} \right.$$

i.e., under appropriate conditions, we will show that there exists a constant t ₀ such that the problem above has no solution if t > t ₀, at least a solution if t =  t ₀ and at least two solutions if t < t ₀. The proof is based on a combination of upper and lower solutions method and the Leray–Schauder degree.

     

Karhunen–Loève expansion for a generalization of Wiener bridge

We derive a Karhunen–Loève expansion of the Gauss process \( {B}_t-g(t){\int}_0^1{g}^{\hbox{'}}(u)\mathrm{d}{B}_u,t\in \left[0,1\right] \), where (B_t)_t?∈?[0,?1] is a standardWiener process, and g?:?[0,?1]?→?? is a twice continuously differentiable function with g(0) = 0 and \( {\int}_0^1{\left(g\hbox{'}(u)\right)}^2\mathrm{d}u=1 \). This process is an important limit process in the theory of goodness-of-fit tests. We formulate two particular cases with the functions \( g(t)=\left(\sqrt{2}/\pi \right)\sin \left(\pi t\right),t\in \left[0,1\right] \), and g(t)?=?t, t?∈?[0,?1]. The latter corresponds to the Wiener bridge over [0, 1] from 0 to 0.     

<Emphasis Type="Italic">L</Emphasis>-Approximation of a linear combination of the poisson kernel and its conjugate kernel by trigonometric polynomials

A linear combination Π _q,α = cos(απ/2)P + sin(απ/2)Q of the Poisson kernel P(t) = 1/2 + q cos t + q ² cos 2t + ... and its conjugate kernel Q(t) = q sin t + q ² sin 2t + ... is considered for α ∈ ? and |q| < 1. A new explicit formula is found for the value E _n?1(Π _q,α ) of the best approximation in the space L = L _2π of the function Π _q,α by the subspace of trigonometric polynomials of order at most n ? 1. More exactly, it is proved that \(E_{n - 1} \left( {\prod _{q,\alpha } } \right) = \left. {\frac{{\left| q \right|^n \left( {1 - q^2 } \right)}}{{1 - q^{4n} }}} \right\|\left. {\frac{{\cos \left( {nt - {{\alpha \pi } \mathord{\left/ {\vphantom {{\alpha \pi } 2}} \right. \kern-\nulldelimiterspace} 2}} \right) - q^{2n} \cos \left( {nt + {{\alpha \pi } \mathord{\left/ {\vphantom {{\alpha \pi } 2}} \right. \kern-\nulldelimiterspace} 2}} \right)}}{{1 - q^2 - 2q \cos t}}} \right\|_L\). In addition, the value E _n?1(Π _q,α ) is represented as a rapidly convergent series.     

Generalized <Emphasis Type="Italic">n</Emphasis>-Laplacian: boundedness of weak solutions to the Dirichlet problem with nonlinearity in the critical growth range

Let n ≥ 2 and let Ω ? ? ⁿ be an open set. We prove the boundedness of weak solutions to the problem

$$u \in W_0^1 L^\Phi \left( \Omega \right) and - div\left( {\Phi '\left( {\left| {\nabla u} \right|} \right)\frac{{\nabla u}}{{\left| {\nabla u} \right|}}} \right) + V\left( x \right)\Phi '\left( {\left| u \right|} \right)\frac{u}{{\left| u \right|}} = f\left( {x,u} \right) + \mu h\left( x \right) in \Omega ,$$

where ? is a Young function such that the space W ₀ ¹ L ^Φ(Ω) is embedded into an exponential or multiple exponential Orlicz space, the nonlinearity f(x, t) has the corresponding critical growth, V(x) is a continuous potential, h ∈ L ^Φ(Ω) is a non-trivial continuous function and µ ≥ 0 is a small parameter. We consider two classical cases: the case of Ω being an open bounded set and the case of Ω = ? ⁿ .

     

Existence of Nonnegative Solutions for a Class of Systems Involving Fractional (p,q)-Laplacian Operators

The authors study the following Dirichlet problem of a system involving fractional (p, q)-Laplacian operators:

$$\left\{ {\begin{array}{*{20}{c}} {\left( { - \Delta } \right)_p^su = \lambda a\left( x \right){{\left| u \right|}^{p - 2}}u + \lambda b\left( x \right){{\left| u \right|}^{\alpha - 2}}{{\left| v \right|}^\beta }u + \frac{{\mu \left( x \right)}}{{\alpha \delta }}{{\left| u \right|}^{\gamma - 2}}{{\left| v \right|}^\delta }uin\Omega ,} \\ {\left( { - \Delta } \right)_q^sv = \lambda c\left( x \right){{\left| v \right|}^{q - 2}}v + \lambda b\left( x \right){{\left| u \right|}^\alpha }{{\left| v \right|}^{\beta - 2}}v + \frac{{\mu \left( x \right)}}{{\beta \gamma }}{{\left| u \right|}^\gamma }{{\left| v \right|}^{\delta - 2}}vin\Omega ,} \\ {u = v = 0on{\mathbb{R}^N}\backslash \Omega ,} \end{array}} \right.$$

where λ > 0 is a real parameter, Ω is a bounded domain in R ^N , with boundary ?Ω Lipschitz continuous, s ∈ (0, 1), 1 < p ≤ q < ∞, sq < N, while (?Δ) _p ^s u is the fractional p-Laplacian operator of u and, similarly, (?Δ) _q ^s v is the fractional q-Laplacian operator of v. Since possibly p ≠ q, the classical definitions of the Nehari manifold for systems and of the Fibering mapping are not suitable. In this paper, the authors modify these definitions to solve the Dirichlet problem above. Then, by virtue of the properties of the first eigenvalue λ₁ for a related system, they prove that there exists a positive solution for the problem when λ < λ₁ by the modified definitions. Moreover, the authors obtain the bifurcation property when λ → λ₁^-. Finally, thanks to the Picone identity, a nonexistence result is also obtained when λ ≥ λ₁.

     

A Spectral Gap Estimate and Applications

We consider the Schrödinger operator

$$ \text{-} \frac{d^{2}}{d x^{2}} + V {\text{on an interval}}~~[a,b]~{\text{with Dirichlet boundary conditions}},$$

where V is bounded from below and prove a lower bound on the first eigenvalue λ ₁ in terms of sublevel estimates: if w _V (y) = |{x ∈ [a, b] : V (x) ≤ y}|, then

$$\lambda_{1} \geq \frac{1}{250} \min\limits_{y > \min V}{\left( \frac{1}{w_{V}(y)^{2}} + y\right)}.$$

The result is sharp up to a universal constant if {x ∈ [a, b] : V(x) ≤ y} is an interval for the value of y solving the minimization problem. An immediate application is as follows: let \({\Omega } \subset \mathbb {R}^{2}\) be a convex domain and let \(u:{\Omega } \rightarrow \mathbb {R}\) be the first eigenfunction of the Laplacian ? Δ on Ω with Dirichlet boundary conditions on ?Ω. We prove

$$\| u \|_{L^{\infty}({\Omega})} \lesssim \frac{1}{\text{inrad}({\Omega})} \left( \frac{\text{inrad}({\Omega})}{\text{diam}({\Omega})} \right)^{1/6} \|u\|_{L^{2}({\Omega})},$$

which answers a question of van den Berg in the special case of two dimensions.

     

Asymptotic behavior of the p-torsion functions as p goes to 1

Let \({\Omega}\) be a Lipschitz bounded domain of \({\mathbb{R}^N}\), \({N\geq2}\), and let \({u_p\in W_0^{1,p}(\Omega)}\) denote the p-torsion function of \({\Omega}\), p > 1. It is observed that the value 1 for the Cheeger constant \({h(\Omega)}\) is threshold with respect to the asymptotic behavior of u_p, as \({p\rightarrow 1^+}\), in the following sense: when \({h(\Omega) > 1}\), one has \({\lim_{p\rightarrow 1^+}\left\|u_{p}\right\| _{L^\infty(\Omega)}=0}\), and when \({h(\Omega) < 1}\), one has \({\lim_{p\rightarrow 1^+}\left\|u_p\right\| _{L^\infty(\Omega)}=\infty}\). In the case \({h(\Omega)=1}\), it is proved that \({\limsup_{p\rightarrow1^+}\left\|u_p\right\|_{L^\infty(\Omega)}<\infty}\). For a radial annulus \({\Omega_{a,b}}\), with inner radius a and outer radius b, it is proved that \({\lim_{p\rightarrow 1^+}\left\|u_p\right\| _{L^\infty(\Omega_{a,b})}=0}\) when \({h(\Omega_{a,b})=1}\).     

Divergent solutions to the <Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript>-supercritical NLS equations

We investigate the nonlinear Schrödinger equation iu _t +Δu+|u| ^p?1 u = 0with 1+ 4/N < p < 1+ 4/N?2 (when N = 1, 2, 1 + 4/N < p < ∞) in energy space H ¹ and study the divergent property of infinite-variance and nonradial solutions. If \(M{\left( u \right)^{\frac{{1 - {s_C}}}{{{s_C}}}}}E\left( u \right) \prec M{\left( Q \right)^{\frac{{1 - {s_C}}}{{{s_C}}}}}E\left( Q \right)\) and \(\left\| {{u_0}} \right\|_2^{\frac{{1 - {s_c}}}{{{s_c}}}}\left\| {\nabla {u_0}} \right\|_2^{\frac{{1 - {s_c}}}{{{s_c}}}}{\left\| {\nabla Q} \right\|_2}\), then either u(t) blows up in finite forward time or u(t) exists globally for positive time and there exists a time sequence t _n → +∞ such that \({\left\| {\nabla u\left( {{t_n}} \right)} \right\|_2} \to + \infty \). Here Q is the ground state solution of ?(1?s _c )Q+ΔQ+Q ^p?1 Q = 0. A similar result holds for negative time. This extend the result of the 3D cubic Schrödinger equation obtained by Holmer to the general mass-supercritical and energy-subcritical case.     

Several extremal approximation problems for the characteristic function of a spherical layer

We discuss three interrelated extremal problems on the set P _n,m of algebraic polynomials of a given degree n on the unit sphere \(\mathbb{S}^{m - 1}\) of the Euclidean space ? ^m of dimension m ≥ 2. (1) Find the norm of the functional \(F\left( \eta \right) = F_h P_n = \int_{\mathbb{G}\left( \eta \right)} {P_n (x)dx}\), which is the integral over the spherical layer \(\mathbb{G}\left( \eta \right) = \left\{ {x = \left( {x_1 , \ldots ,x_m } \right) \in \mathbb{S}^{m - 1} :h' \leqslant x_m \leqslant h''} \right\}\) defined by a pair of real numbers η = (h′, h″), ?1 ≤ h′ < h″ ≤ 1, on the set P _n,m with the norm of the space \(L\left( {\mathbb{S}^{m - 1} } \right)\) of functions summable on the sphere. (2) Find the best approximation in \(L_\infty \left( {\mathbb{S}^{m - 1} } \right)\) of the characteristic function χ _η of the layer \(\mathbb{G}\left( \eta \right)\) by the subspace P _n,m ^⊥ of functions from \(L_\infty \left( {\mathbb{S}^{m - 1} } \right)\) that are orthogonal to the space of polynomials P _n,m . (3) Find the best approximation in the space \(L\left( {\mathbb{S}^{m - 1} } \right)\) of the function χ _η by the space of polynomials P _n,m . We present a solution of all three problems for the values h′ and h″ that are neighboring roots of the polynomial in a single variable of degree n + 1 that deviates least from zero in the space L ₁ ^φ (?1, 1) of functions summable on the interval (?1, 1) with ultraspherical weight φ(t) = (1 ? t ²) ^α , α = (m ? 3)/2. We study the respective one-dimensional problems in the space of functions summable on (?1, 1) with an arbitrary not necessarily ultraspherical weight.     

Estimates for sums of moduli of blocks in trigonometric Fourier series

We consider the following two problems. Problem 1: what conditions on a sequence of finite subsets A _k ? ? and a sequence of functions λ _k : A _k → ? provide the existence of a number C such that any function f ∈ L ₁ satisfies the inequality ‖U _A,Λ(f)‖ _p ≤ C‖f‖₁ and what is the exact constant in this inequality? Here, \(U_{\mathcal{A},\Lambda } \left( f \right)\left( x \right) = \sum\nolimits_{k = 1}^\infty {\left| {\sum\nolimits_{m \in A_k } {\lambda _k \left( m \right)c_m \left( f \right)e^{imx} } } \right|}\) and c _m (f) are Fourier coefficients of the function f ∈ L ₁. Problem 2: what conditions on a sequence of finite subsets A _k ? ? guarantee that the function \(\sum\nolimits_{k = 1}^\infty {\left| {\sum\nolimits_{m \in A_k } {c_m \left( h \right)e^{imx} } } \right|}\) belongs to L _p for every function h of bounded variation?     

Harmonic analysis associated with a discrete Laplacian

It is well known that the fundamental solution of

$${u_t}\left( {n,t} \right) = u\left( {n + 1,t} \right) - 2u\left( {n,t} \right) + u\left( {n - 1,t} \right),n \in \mathbb{Z},$$

with u(n, 0) = δ _nm for every fixed m ∈ Z is given by u(n, t) = e ^?2t I _n?m (2t), where I _k (t) is the Bessel function of imaginary argument. In other words, the heat semigroup of the discrete Laplacian is described by the formal series W _t f(n) = Σ _m∈Z e ^?2t I _n?m (2t)f(m). This formula allows us to analyze some operators associated with the discrete Laplacian using semigroup theory. In particular, we obtain the maximum principle for the discrete fractional Laplacian, weighted ? ^p (Z)-boundedness of conjugate harmonic functions, Riesz transforms and square functions of Littlewood-Paley. We also show that the Riesz transforms essentially coincide with the so-called discrete Hilbert transform defined by D. Hilbert at the beginning of the twentieth century. We also see that these Riesz transforms are limits of the conjugate harmonic functions. The results rely on a careful use of several properties of Bessel functions.

     

On the nonexistence of solutions of some anisotropic elliptic and backward parabolic inequalities

We study the nonexistence of weak solutions of higher-order elliptic and parabolic inequalities of the following types: \(\sum {_{i = 1}^N\sum\nolimits_{{e_i} \leqslant {\alpha _i} \leqslant {m_i}} {D_{{x_i}}^{{\alpha _i}}\left( {{A_{{\alpha _i}}}\left( {x,u} \right)} \right)} \geqslant f\left( {x,u} \right),} x \in {\mathbb{R}^N}\), and \({u_t} + \sum {_{i = 1}^N\sum\nolimits_{{k_i} \leqslant {\beta _i} \leqslant {n_i}} {D_{{x_i}}^{{\beta _i}}\left( {{B_{{\beta _i}}}\left( {x,t,u} \right)} \right)} > g\left( {x,t,u} \right),\left( {x,t} \right)} \in {\mathbb{R}^N} \times {\mathbb{R}_ + }\), where l _i , m _i , k _i , n _i ∈ N satisfy the condition l _i , k _i > 1 for all i = 1,..., N, and A _αi (x, u), B _βi (x, t, u), f(x, u), and g(x, t, u) are some given Carathéodory functions. Under appropriate conditions on the functions A _αi , B _βi , f, and g, we prove theorems on the nonexistence of solutions of these inequalities.     

Some Berezin Number Inequalities for Operator Matrices

The Berezin symbol à of an operator A acting on the reproducing kernel Hilbert space H = H(Ω) over some (nonempty) set is defined by \(\tilde A(\lambda ) = \left\langle {A\hat k_\lambda ,\hat k_\lambda } \right\rangle \), λ ∈ Ω, where \(\hat k_\lambda = k_\lambda /\left\| {k_\lambda } \right\|\) is the normalized reproducing kernel of H. The Berezin number of the operator A is defined by \(ber(A) = \mathop {\sup }\limits_{\lambda \in \Omega } \left| {\tilde A(\lambda )} \right| = \mathop {\sup }\limits_{\lambda \in \Omega } \left| {\left\langle {A\hat k_\lambda ,\hat k_\lambda } \right\rangle } \right|\). Moreover, ber(A) ? w(A) (numerical radius). We present some Berezin number inequalities. Among other inequalities, it is shown that if \(T = \left[ {\begin{array}{*{20}c} A & B \\ C & D \\ \end{array} } \right] \in \mathbb{B}(\mathcal{H}(\Omega _1 ) \oplus \mathcal{H}(\Omega _2 ))\), then

$$ber(T) \leqslant \frac{1}{2}(ber(A) + ber(D)) + \frac{1}{2}\sqrt {(ber(A) - ber(D))^2 + \left( {\left\| B \right\| + \left\| C \right\|} \right)^2 } .$$

     

Estimates for the Torsion Function and Sobolev Constants

Let Ω be an open set in Euclidean space, and let u : Ω → ?^?+? be the expected lifetime of Brownian motion in Ω. It is shown that if u?∈?L ^p (Ω) for some p?∈?[1, ?∞?) then (i) u?∈?L ^q (Ω) for all q?∈?[p,?∞?], and (ii) \({trace}\left(e^{t\Delta_{\Omega}}\right)<\infty\) for all t?>?0, where ??Δ_Ω is the Dirichlet Laplacian acting in L ²(Ω). Pointwise bounds are obtained for u in terms of the first Dirichlet eigenfunction for Ω, assuming that the spectrum of ??Δ_Ω is discrete. It is shown that if Ω is open, bounded and connected in the plane and \(\partial\Omega\) has an interior wedge with opening angle α at vertex v then the first Dirichlet eigenfunction and u are comparable near v if and only if α?≥?π/2. Two sided estimates are obtained for the Sobolev constant

$ C_p(\Omega):= \inf\left\{\Vert \nabla u \Vert_2^2: u \in C_0^{\infty}(\Omega),\ \Vert u\Vert_p = 1\right\}, $

where 0?p?Ω satisfies a strong Hardy inequality, and the distance to the boundary function δ?∈?L ^2p/(2???p)(Ω).

     

Extracting outer function part from Hardy space function

Any analytic signal fa(e~(it)) can be written as a product of its minimum-phase signal part(the outer function part) and its all-phase signal part(the inner function part). Due to the importance of such decomposition, Kumarasan and Rao(1999), implementing the idea of the Szeg?o limit theorem(see below),proposed an algorithm to obtain approximations of the minimum-phase signal of a polynomial analytic signal fa(e~(it)) = e~(iN0t)M∑k=0a_k~(eikt),(0.1)where a_0≠ 0, a_M≠ 0. Their method involves minimizing the energy E(f_a, h_1, h_2,..., h_H) =1/(2π)∫_0~(2π)|1+H∑k=1h_k~(eikt)|~2|fa(e~(it))|~2dt(0.2) with the undetermined complex numbers hk's by the least mean square error method. In the limiting procedure H →∞, one obtains approximate solutions of the minimum-phase signal. What is achieved in the present paper is two-fold. On one hand, we rigorously prove that, if fa(e~(it)) is a polynomial analytic signal as given in(0.1),then for any integer H≥M, and with |fa(e~(it))|~2 in the integrand part of(0.2) being replaced with 1/|fa(e~(it))|~2,the exact solution of the minimum-phase signal of fa(e~(it)) can be extracted out. On the other hand, we show that the Fourier system e~(ikt) used in the above process may be replaced with the Takenaka-Malmquist(TM) system, r_k(e~(it)) :=((1-|α_k|~2e~(it))/(1-α_ke~(it))~(1/2)∏_(j=1)~(k-1)(e~(it)-α_j/(1-α_je~(it))~(1/2), k = 1, 2,..., r_0(e~(it)) = 1, i.e., the least mean square error method based on the TM system can also be used to extract out approximate solutions of minimum-phase signals for any functions f_a in the Hardy space. The advantage of the TM system method is that the parameters α_1,..., α_n,...determining the system can be adaptively selected in order to increase computational efficiency. In particular,adopting the n-best rational(Blaschke form) approximation selection for the n-tuple {α_1,..., α_n}, n≥N, where N is the degree of the given rational analytic signal, the minimum-phase part of a rational analytic signal can be accurately and efficiently extracted out.     

Entropy solution for a nonlinear elliptic problem with lower order term in Musielak–Orlicz spaces

In this paper, we shall be concerned with the existence result of the following problem,

$$\begin{aligned} \left\{ \begin{array}{l} -\text {div}\left( a(x,u,\nabla u)\right) -\text {div}(\Phi (x,u))= f \ \ \mathrm{in}\ \Omega ,\\ u=0 \text { on } \partial \Omega , \end{array} \right. \end{aligned}$$

(0.1)

with the second term f belongs to \(L^1(\Omega )\). The growth and the coercivity conditions on the monotone vector field a are prescribed by a generalized N-function M. We assume any restriction on M, therefore we work with Musielak–Orlicz spaces which are not necessarily reflexive. The lower order term \(\Phi \) is a Carathéodory function satisfying only a growth condition.

     

Compactness and existence results in weighted Sobolev spaces of radial functions. Part II: existence

We apply the compactness results obtained in the first part of this work, to prove existence and multiplicity results for finite energy solutions to the nonlinear elliptic equation

$$-\triangle u + V \left(\left|x\right|\right) u = g \left(\left|x\right|, u\right) \quad {\rm in} \Omega \subseteq \mathbb{R}^{N},\,N \geq 3,$$

where \({\Omega}\) is a radial domain (bounded or unbounded) and u satisfies u =  0 on \({\partial\Omega}\) if \({\Omega \neq\mathbb{R}^{N}}\) and \({u \rightarrow 0}\) as \({\left|x\right| \rightarrow \infty}\) if \({\Omega}\) is unbounded. The potential V may be vanishing or unbounded at zero or at infinity and the nonlinearity g may be superlinear or sublinear. If g is sublinear, the case with a forcing term \({g\left(\left|\cdot\right|, 0\right) \neq 0}\) is also considered. Our results allow to deal with V and g exhibiting behaviours at zero or at infinity which are new in the literature and, when \({\Omega = \mathbb{R}^{N}}\), do not need to be compatible with each other.

     

Divided differences for functions of two variables for irregularly spaced arguments

Divided differences forf (x, y) for completely irregular spacing of points (x _i ,y _i ) are developed here by a natural generalization of Newton's scheme. Existing bivariate schemes either iterate the one-dimensional scheme, thus constraining (x _i ,y _i ) to be at corners of rectangles, or give polynomials Σa _jk x ^j y ^k having more coefficients than interpolation conditions. Here the generalizedn ^th divided difference is defined by (1)\(\left[ {01... n} \right] = \sum\limits_{i = 0}^n {A_i f\left( {x_i , y_i } \right)} \) where (2)\(\sum\limits_{i = 0}^n {A_i x_i^j , y_i^k = 0} \), and 1 for the last or (n+1)^th equation, for every (j, k) wherej+k=0, 1, 2,... in the usual ascending order. The gen. div. diff. [01...n] is symmetric in (x _i ,y _i ), unchanged under translation, 0 forf (x, y) an, ascending binary polynomial as far asn terms, degree-lowering with respect to (X, Y) whenf(x, y) is any polynomialP(X+x, Y+y), and satisfies the 3-term recurrence relation (3) [01...n]=λ{[1...n]?[0...n?1]}, where (4) λ= |1...n|·|01...n?1|/|01...n|·|1...n?1|, the |...i...| denoting determinants inx _i ^j y _i ^k . The generalization of Newton's div. diff. formula is (5)

$$\begin{gathered} f\left( {x, y} \right) = f\left( {x_0 , y_0 } \right) - \frac{{\left| {\alpha 0} \right|}}{{\left| 0 \right|}}\left[ {01} \right] + \frac{{\left| {\alpha 01} \right|}}{{\left| {01} \right|}}\left[ {012} \right] - \frac{{\left| {\alpha 012} \right|}}{{\left| {012} \right|}}\left[ {0123} \right] + \cdots + \hfill \\ + \left( { - 1} \right)^n \frac{{\left| {\alpha 01 \ldots n - 1} \right|}}{{\left| {01 \ldots n - 1} \right|}}\left[ {01 \ldots n} \right] + \left( { - 1} \right)^{n + 1} \frac{{\left| {\alpha 01 \ldots n} \right|}}{{\left| {01 \ldots n} \right|}}\left[ {01 \ldots n} \right], \hfill \\ \end{gathered} $$