On a generalization of the Gauss formula for numerical integration

Abstract

?

Authors

Original title (in Romanian)

T. Popoviciu, Asupra unei generalizări a formulei de integrare numerică a lui Gauss, Studii şi cerc. şt. Iaşi, 6 (1955) nos. 1-2, pp. 29-57 (in Romanian)

Keywords

Gauss formula, numerical integration

Cite this paper as

T. Popoviciu, Asupra unei generalizări a formulei de integrare numerică a lui Gauss, Studii şi cerc. şt. Iaşi, 6 (1955) nos. 1-2, pp. 29-57 (in Romanian)

PDF

About this paper

Journal

Studii şi cercetari ştiintifice Iaşi

Publisher Name

??

JSTOR

Print ISSN

not available

Online ISSN

not available

??

Paper (preprint) in HTML form

Original text
Rate this translation
Your feedback will be used to help improve Google Translate

ON A GENERALIZATION OF GAUSS'S NUMERICAL INTEGRATION FORMULA

BY
TIBERIU POPOVICIU
Communication presented on May 4, 1955 at the meeting of the 1955 Branch of the Academy of the Romanian Republic

§ 1. Gaussian-type formulas.

  1. 1.

    In practical numerical calculation problems, the value of a linear functional is often required 1 )A[f]\mathrm{A}[f], defined on a vector spaceSSof functionsf=f(x)f=f(x), real, of the real variablexx, defined and continuous on an interval I. In what follows we will assume that the elements ofSSare differentiable a sufficient number of times, at least at the points where these derivatives will occur. We will also assume thatSScontains all polynomials inxxIn the following we will impose the functionalityA[f]\vec{A}[f]and certain restrictive conditions, which we will specify when they occur.

  2. 2.

    Suppose that the values ​​are given

f(j)(xand),j=0,1,,Rand1,and=1,2,,n(f(0)(x)=f(x)),f^{(j)}\left(x_{i}\right),j=0,1,\ldots,r_{i}-1,i=1,2,\ldots,n\quad\left(f^{(0)}(x)=f(x)\right), (1)

of the functionf(x)f(x)and its first derivativesf(x),f"(x),f^{\prime}(x),f^{\prime\prime}(x),\ldotson distinct points
(2)

x1,x2,,xnx_{1},x_{2},\ldots,x_{n}

of interval I. On the pointxandx_{i}the values ​​of the function and its primes are givenRand1r_{i}-1derived, so that the numbersRand,and=1,2,,nr_{i},i=1,2,\ldots,nare entirely positive.

For eachfS,A¯[f]f\in S,\bar{A}[f]can be approximated by a given linear combination of the values ​​(1) of the functionf(x)f(x)and its first derivatives on points (2). We thus obtain the approximation formula

A[f]and=1nj=0Rand1Aand,jf(j)(xand).\mathrm{A}[f]\approx\sum_{i=1}^{n}\sum_{j=0}^{r_{i}-1}a_{i,j}f^{(j)}\left(x_{i}\right). (3)

The points (2) are the nodes of this formula and the numbersRand,and=1,2,,nr_{i},i=1,2,\ldots,nare the multiplicity orders of these nodes. The nodeRandr_{i}has the order to

0 0 footnotetext: 1. By a linear functional we mean an additive and homogenic functional.

multiplicityRand,and=1,2,,nr_{i},i=1,2,\ldots,n. The numbersAand,j,j=0,1,,Rand1,and=1,2,,na_{i,j},j=0,1,\ldots,r_{i}-1,i=1,2,\ldots,ncharacterizes the approximation formula (3) of this type and can be called the coefficients of this formula.

RESTR[f]\mathrm{R}[f]of formula (3) is, by definition, the difference between the first and second members of the formula. If we then add the remainder to the second member of the formula, this approximate equality becomes an ordinary equality. The remainder is, obviously, also a linear functional defined onSS.
3. Another interpretation of the approximation formula (3) can be given. A general method for finding an approximate value forA[f]\mathrm{A}[f]consists of replacing the functionf(x)f(x)through another functionφS\varphi\in Sand take theA[φ]\mathrm{A}[\varphi]as an approximate value forA[f]\mathrm{A}[f], soA[f]A[φ]\mathrm{A}[f]\approx\mathrm{A}[\varphi].

From a theoretical point of view, nothing prevents us from choosing the functionφ(x)\varphi(x)completely arbitrary, with the only restriction that it belongs toSS. But from the point of view of practical applications, the choice of functionsφ(x)\varphi(x)must and can generally be considerably restricted. An important case is whenφ(x)=B[fx]\varphi(x)=\mathrm{B}[f\mid x]is a given linear operator, with values ​​inSS. In this caseA[φ]=A[B[fx]]\mathrm{A}[\varphi]=\mathrm{A}[\mathrm{B}[f\mid x]]is a linear functional offf, well determined and defined onSSThe restR[f]=A[f]A[φ]\mathrm{R}[f]=\mathrm{A}[f]-\mathrm{A}[\varphi]will then also be a well-determined and defined linear functional onSS.

An important class of functionsφ(x)\varphi(x)of the previous form is formed by the generalized linear interpolation functions

φ(x)=B[fx]=and=1nj=0Rand1φand,j(x)f(j)(xand)\varphi(x)=\mathrm{B}[f\mid x]=\sum_{i=1}^{n}\sum_{j=0}^{ri-1}\varphi_{i,j}(x)f^{(j)}\left(x_{i}\right) (4)

corresponding to the nodes (2), with the indicated multiplicity orders and whereφand,j(x)S,j=0,1,,Rand1,and=1,2,,n\varphi_{i,j}(x)\in S,j=0,1,\ldots,r_{i}-1,i=1,2,\ldots,nare some given functions. Then the formulaA[f]A[φ]A[f]\approx A[\varphi]reduces to (3), where the coefficients are given by the formulas

Aand,j=A[φand,j],j=0,1,,Rand1,and=1,2,,n.a_{i,j}=\mathrm{A}\left[\varphi_{i},j\right],j=0,1,\ldots,r_{i}-1,i=1,2,\ldots,n.

In particular, either

φ(x)=IT(x1,x1,,x1R1,x2,x2,,x2R2,,xn,xn,,xnRn;fx)\varphi(x)=\mathrm{L}(\underbrace{x_{1},x_{1},\ldots,x_{1}}_{r_{1}},\underbrace{x_{2},x_{2},\ldots,x_{2}}_{r_{2}},\ldots,\underbrace{x_{n},x_{n},\ldots,x_{n}}_{r_{n}};f\mid x)

the Lagrange-Hermite interpolation polynomial of degree

p=R1+R2++Rn1,p=r_{1}+r_{2}+\cdots+r_{n}-1, (5)

which, together with its firstRand1r_{i}-1derivatives, the respective values ​​(1) on the nodesxand,and=1,2,,nx_{i},i=1,2,\ldots,nThis functionφ(x)\varphi(x)is of the form (4), where

φand,j(x)==itand,j(x),j=0,1,,Rand1,and=1,2,,n\varphi_{i,j}(x)==l_{i,j}(x),j=0,1,\ldots,r_{i}-1,i=1,2,\ldots,n

are the fundamental interpolation polynomials relative to the nodes (2) with their respective multiplicity orders. These polynomials are completely determined. They have some well-known expressions, among which we recall the formulas
(6)itand,Rand1(x)=1(Rand1)!j=1jandn(xandxj)Rjit(x)xxand,and=1,2,,n\quad l_{i},r_{i-1}(x)=\frac{1}{\left(r_{i}-1\right)!\prod_{\begin{subarray}{c}j=1\\ j\neq i\end{subarray}}^{n}\left(x_{i}-x_{j}\right)r_{j}}\cdot\frac{l(x)}{x-x_{i}},i=1,2,\ldots,n
where
(7)

it(x)=(xx1)R1(xx2)R2(xxnit)Rn.l(x)=\left(x-x_{1}\right)^{r_{1}}\left(x-x_{2}\right)^{r_{2}}\ldots\left(x-x_{nl}\right)^{r_{n}}.

In the considered case, formula (3) becomes
(8)

A[f]and=1nj=0RandR1candjf(j)(xand)A[f]\approx\sum_{i=1}^{n}\sum_{j=0}^{r_{i}r^{-1}}c_{ij}f^{(j)}\left(x_{i}\right)

where

cand,j=𝐀𝐧𝐞𝐫[itand,j],j=0,1,,Rand1,and=1,2,,nc_{i,j}=\mathbf{A}\left[l_{i},j\right],j=0,1,\ldots,r_{i}-1,i=1,2,\ldots,n
  1. 4.

    The rest of formula (8) enjoys the important property that it cancels out for any polynomial of degreep,pp,pbeing given by formula (5).

In general if the linear functional𝔸[f]\mathbb{A}[f]is zero for any polynomial of degreem2m^{2}), but is nonzero for at least one polynomial of degreem+1m+1, it is said that it has the degree of accuracymmThis definition naturally extends to the casesm=1m=-1andm=+m=+\inftyMore precisely, the degree of accuracy is an integer1\geq-1or the wrong number++\inftyattached to the functionalityA[f]\mathrm{A}[f]and perfectly characterized by prosperity:
10.m=11^{0}.m=-1, ifA[1]0A[1]\neq 0,
20.A[1]=A[x]==A[xm]=0,A[xn+1]02^{0}.\mathrm{A}[1]=\mathrm{A}[x]=\cdots=\mathrm{A}\left[x^{m}\right]=0,\mathrm{~A}\left[x^{n+1}\right]\neq 0, ifA[1]=0\mathrm{A}[1]=0and at least one of the numbersA[xand],and=1,2,\mathrm{A}\left[x^{i}\right],i=1,2,\ldotsis different from zero,
30.m=+3^{0}.m=+\infty, ifA[xand]=0,and=0,1,\mathrm{A}\left[x^{i}\right]=0,i=0,1,\ldots
Finally, we will say that the functional𝔸[f]\mathbb{A}[f]has the degree of accuracy at leastmm, if it has the degree of accuracym\geq m, so if it cancels out for any polynomial of degreemm.

For simplicity, we will say that an approximation formula (3) has the degree of accuracymmrespectively has the degree of accuracy at leastmm, if the rest of this formula has the degree of accuracymmrespectively has the degree of accuracy at leastmm.

With this convention we can say that formula (8) has the degree of accuracy at leastpp.

We recall the following:
Theorem 1. If formula (3) has the degree of accuracy at least p, this formula necessarily coincides with (8).

I gave the proof of this property in another paper [9]. The proof there referred to a linear functionalA[f]A[f]particular, but it does not depend on the form of this functional.

Theorem 1 tells us that formula (8) plays a special role among the formulas (3) that are obtained by varying the coefficients of this formula. Formula (8) is, among all the formulas (3) corresponding to a functionalA[f]\mathrm{A}[f]and some nodes (2), given together with their respective multiplicity orders, the (unique) one that has the maximum degree of accuracy.
5. In general, the degree of accuracy of formula (8) isppIn particular cases, however, this degree of accuracy may be even greater thanpp.

Definition. We will say that formula (8) is of Gaussian type if its degree of accuracy is at leastp+np+n.

0 0 footnotetext: 2. Through polynomial of degree mm, we understand a polynomial of effective degreem\leqq mA polynomial of degree 0 is a constant, and a polynomial of degree -1, identical to the zero polynomial,

For formula (8) to be Gaussian it is necessary and sufficient that its remainderR[f]\mathrm{R}[f]to cancel out for any polynomial of degreep+np+n.

A polynomialP(x)P(x)of the degreep+np+nis always of the form

P(x)=it(x)Q(x)+Q1(x),P(x)=l(x)Q(x)+Q_{1}(x),

whereit(x)l(x)is the polynomial (7),Q(x)Q(x)a polynomial of degreen1n-1andQ1(x)Q_{1}(x)a polynomial of degreepp. Conversely, any polynomial of this form is of degreep+np+n. ThenR[P]=R[itQ]R[P]=R[lQ]and from (8) it immediately followsR(itQ)=A[itQ]R(lQ)=A[lQ]It follows therefore that the necessary and sufficient condition for formula (4) to be Gaussian is that we have

(10)

A[itQ]=0\mathrm{A}[l\mathrm{Q}]=0

whatever the polynomial isQ(x)\mathrm{Q}(x)of the degreen1n-1.
Two functionsf,gf,gfor which we haveOh[fg]=0\AA [fg]=0can be called orthogonal to the functionalA[f]\mathrm{A}[f]We can therefore state:

Theorem 2. The necessary and sufficient condition for formula (8) to be of Gaussian type is that polynomial (7) is orthogonal to any polynomial of degreen1n-1:

Orthogonality of the polynomialit(x)l(x)with any polynomial of degreen1n-1is equivalent to its orthogonality withnnpolynomials of degreen1n-1linearly independent. Such a system ofnnpolynomials is formed from the firstnnpower1,x,x2,,xn11,x,x^{2},\ldots,x^{n-1}his/hersxxAnother system of this kind is formed by the polynomials

(xx1)(xx2)(xxand1)(xxand+1)(xxn),and=1,2,,n,\left(x-x_{1}\right)\left(x-x_{2}\right)\ldots\left(x-x_{i-1}\right)\left(x-x_{i+1}\right)\ldots\left(x-x_{n}\right),i=1,2,\ldots,n, (11)

because we assume that the nodes are distinct.
This last example shows us that a Gaussian formula (8) is nothing more than a formula of the form (8)

A[f]and=1nj=0Randcand,jf(j)(xand),\mathrm{A}[f]\approx\sum_{i=1}^{n}\sum_{j=0}^{r_{i}}c_{i,j}^{\prime}f^{(j)}\left(x_{i}\right), (12)

corresponding to nodes (2), but of multiplicity orderRand+1r_{i}+1(instead ofRand),and=1,2,,n\left.r_{i}\right),i=1,2,\ldots,nrespectively$1\mathdollar 1in which we havec,andRand=0,and=1,2,.nc^{\prime}{}_{i},r_{i}=0,i=1,2,\ldots.n. This observation is due in principle to AA Markov [4]. The property results from formulas (6); (9) corresponding to formula (12).
6. In the following we will assume that the functionalA[f]\mathrm{A}[f], natural numbernnand the orders of multiplicityR1,R2,,Rnr_{1},r_{2},\ldots,r_{n}of the nodes are given.

The orthogonality condition, based on the observation above, can be interpreted in another way. We can look at

φ=φ(x1,x2,,xn)=A[and=1n(xxand)+and1]\phi=\phi\left(x_{1},x_{2},\ldots,x_{n}\right)=\mathrm{A}\left[\prod_{i=1}^{n}\left(x-x_{i}\right)^{\prime}{}_{i}+1\right] (13)

as a (polynomial) function ofx1,x2,,xnx_{1},x_{2},\ldots,x_{n}. Then

1Rand+1φxand=A[it(x)(xx1)(xx2)(xxand1)(xxand+1)(xxn)]and=1,2,,n\begin{gathered}-\frac{1}{r_{i}+1}\cdot\frac{\partial\phi}{\partial x_{i}}=\mathrm{A}\left[l(x)\left(x-x_{1}\right)\left(x-x_{2}\right)\ldots\left(x-x_{i-1}\right)\left(x-x_{i+1}\right)\ldots\left(x-x_{n}\right)\right]\\ i=1,2,\ldots,n\end{gathered}

It follows that the nodes of a Gaussian formula (8) always form a solution of the algebraic system
(14)

φxand=0,and=1,2,,n.\frac{\partial\phi}{\partial x_{i}}=0,\quad i=1,2,\ldots,n.

Any solution to this system, in which the values ​​of the variablesx1,x2,,xnx_{1},x_{2},\ldots,x_{n}are distinct, it gives us a Gaussian formula. However, system (14) does not always have such a (real) solution.

If multiple numbersRandr_{i}are equal, we do not consider as distinct two Gaussian formulas that differ only by a permutation of the nodes having this multiplicity order. In other words, Gaussian formulas depend only on the distinct values ​​of the multiplicity orders.

Based on this observation, one can easily transform system (14) into the equivalent jump from the point of view of searching for Gaussian-type formulas. Thus ifR1,R2,,Rt(1tn)r_{1}^{\prime},r_{2}^{\prime},\ldots,r_{t}^{\prime}(1\leq t\leq n)are the distinct values ​​of the multiplicity ordersR1,R2,,Rnr_{1},r_{2},\ldots,r_{n}, eithernandn_{i}number of nodes of multiplicity orderRandr_{i}^{\prime}$1σ1(and),σ2(and),,σnand(and)1\sigma_{1}^{(i)},\sigma_{2}^{(i)},\ldots,\sigma_{n_{i}}^{(i)}fundamental symmetric functifiles of these knots,and=1,2,,ti=1,2,\ldots,tThen, from the point of view of searching for Gaussian formulas, system (14) is equivalent to

φσj(and)=0,j=1,2,,nand,and=1,2,,t\frac{\partial\phi}{\partial\sigma_{j}^{(i)}}=0,j=1,2,\ldots,n_{i},i=1,2,\ldots,t (15)

This system is as simple as (14) in the sense that the functionφ\phiis a polynomial with respect toσj(and),j=1,2,,nand,and=1,2,,t\sigma_{j}^{(i)},j=1,2,\ldots,n_{i},i=1,2,\ldots,tThe equivalence of the systems (14', (15), in the specified sense, results from the fact that the functional determinant of the fundamental symmetric functionsΣz1z2..zand,and=1,2,,k\Sigma z_{1}z_{2}\ldots..z_{i},i=1,2,\ldots,kof variablesz1,z2,,zkz_{1},z_{2},\ldots,z_{k}with respect to these variables is different from zero, for any system of different values ​​of these variables (see e.g. [6]).

The same can be said about the system deduced from (14), replacing only in part the nodes corresponding to equal multiplicity orders by their fundamental symmetric functionals. The previous method can also be combined with some linear transformations of some or all of the variables.xandx_{i},
7. We must note that for a given linear functionalOh[f]\AA [f]and for a given system of multiplicity orders, there are not always Gaussian-type formulas.

We will say that the functionalA[f]\mathrm{A}[f]is of the order of positivitykkifA[Q2]>0A\left[Q^{2}\right]>0, for any polynomialQ(x)Q(x)of the degreek1k-1and non-identical null.

We can then observe that for a functional of order positivityk12(p+n+2)k\geq\frac{1}{2}(p+n+2)and if at least one of the multiplicity ordersR1,R2,,Rnr_{1},r_{2},\ldots,r_{n}is even, there is no Gaussian formula. Indeed, let us assume, for the sake of clarity, thatRand,Rand+1,,Rnr_{i},r_{i+1},\ldots,r_{n}are even numbers,0<andn0<i\leqq n, the others (if1<and1<i) being odd. Thenit(x)(xx1)(xx2)(xxand1)l(x)\left(x-x_{1}\right)\left(x-x_{2}\right)\ldots\left(x-x_{i-1}\right)is the square of a polynomial of degree12(p+and)k1\frac{1}{2}(p+i)\leq k-1So we have𝔸[it(x)(xx1)(xx2)(xxand1)]>0,(𝔸[it]>0\mathbb{A}\left[l(x)\left(x-x_{1}\right)\left(x-x_{2}\right)\ldots\left(x-x_{i-1}\right)\right]>0,(\mathbb{A}[l]>0, ifand=1)i=1), which, based on Theorem 2 , proves the property.

On the contrary, we will see that if all the numbersR1,R2,,Rnr_{1},r_{2},\ldots,r_{n}are odd, there is at least one Gaussian formula.

Formulas (14) immediately suggest us to examine the relative extrema of the function (13). A relative extremum reached by a system of different values ​​of the variablesx1,x2,,xnx_{1},x_{2},\ldots,x_{n}, we demonstrate the existence of at least one Gauss-type formula. We will see that based on this observation, this existence property holds in particular if𝔸[f]\mathbb{A}[f]is of the order of positivity12(p+n+1)\geq\frac{1}{2}(p+n+1)and if all multiplicity orders are odd.

In the following we will study an extremum problem which, in particular, will give us the solution to the above problem. The properties obtained in this section should be considered as a generalization of the well-known extremal properties of orthogonal polynomials and their generalizations in the sense of G. Polya [5] and D. Jackson [2, 3].

§ 2. On a minimum problem.

  1. 8.

    Let us consider, in particular, a linear functional of the form

A[f]=and=1kλandf(yand),\mathrm{A}[f]=\sum_{i=1}^{k}\lambda_{i}f\left(y_{i}\right), (16)

wherekkis a natural number,y1,y2,,,,yk,ky_{1},y_{2},,,,y_{k},kdistinct points of the real axis andλ1,λ2,,λk\lambda_{1},\lambda_{2},\ldots,\lambda_{k}SYNTHESISkkgiven positive numbers.

We will note withPnP_{n}the set of (real) polynomials of degreennformalxn+x^{n}+\cdots, so with its coefficientxnx^{n}equal to 1.

They are giventtpositive numbersS1,S2,,Sts_{1},s_{2},\ldots,s_{t}andttnatural numbersn1,n2,,ntn_{1},n_{2},\ldots,n_{t}so that each numberSands_{i}What does a number correspond to?nandn_{i}For the sake of brevity, we will call the numbersSands_{i}powers and numbersnandn_{i}the respective degrees corresponding to these powers. These names are justified by the following.

Whether

μ=μn1,n2,,nt(S1,S2,,St)=childA[and=1t|πand|Sand]\mu=\mu_{n_{1},n_{2},\ldots,n_{t}}^{\left(s_{1},s_{2},\ldots,s_{t}\right)}=\inf\mathrm{A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|s_{i}\right] (17)

where the lower bound is relative to all polynomialsπandPnand\pi_{i}\in P_{n_{i}},and=1,2,,ti=1,2,\ldots,t.

Any particular system of polynomialsπandPnand,and=1,2,,t\pi_{i}\in P_{n_{i}},i=1,2,\ldots,tfor whichA[and=1t|πand|Sand]=μA\left[\prod_{i=1}^{t}\left|\pi_{i}\right|s_{i}\right]=\muit will be called a system of minimizing polynomials or, for short, a minimizing system.

caset=1t=1is well known and has been examined in particular by D. Jackson [3].
9. In the case ofn=n1+n2++ntkn=n_{1}+n_{2}+\cdots+n_{t}\geq kHAVEμ=0\mu=0and because the systemπandPnand,and=1,2,,t\pi_{i}\in P_{n_{i}},i=1,2,\ldots,tto be minimizing it is necessary and sufficient that each pointyj,j=1,2,,ky_{j},j=1,2,\ldots,kto be a root of at least one polynomialπand\pi_{i}In the case ofn<kn<kthe results are less trivial and are given by the following:

Theorem 3.A[f]\mathrm{A}[f]being a linear functional of the form (16) with the pointsyandy_{i}distinct and with numbersλand\lambda_{i}all positive andS1,S2,,Sts_{1},s_{2},\ldots,s_{t}a system of (positive) powers andn1,n2,,ntn_{1},n_{2},\ldots,n_{t}a system of degrees corresponding to given values, with their sumn=n1+n2++nt<kn=n_{1}+n_{2}+\cdots+n_{t}<k:
101^{0}There is at least one minimizing system.
202^{0}IfπandPnand,and=1,2,,t\pi_{i}\in P_{ni},i=1,2,\ldots,tis a minimizing system, the polynomialsπand\pi_{i}they all have real roots.

If the powersS1,S2,,Sts_{1},s_{2},\ldots,s_{t}they are all>1>1, then any minimizing systemπandPnand,and=1,2,,t\pi_{i}\in P_{n_{i}},i=1,2,\ldots,talso check the following properties:
303^{0}. All roots of the polynomial
(18)

π=π(x)=and=1tπand\pi=\pi(x)=\prod_{i=1}^{t}\pi_{i}

are distinct.
404^{0}. We have 3 )
(19)

A|and=1t|πand|Sand1(sgπ)Q]=0,\left.\left.\mathrm{A}\left|\prod_{i=1}^{t}\right|\pi_{i}\right|^{s_{i}-1}(\mathrm{sg}\pi)\mathrm{Q}\right]=0,

whatever the polynomial isQ(x)Q(x)of the degreen1n-1.

55^{\circ}The roots of the polynomial (18) are separated by the pointsy1,y2,,yky_{1},y_{2},\ldots,y_{k}.

We note that in the case of Theorem 3,μ>0\mu>0.

10. For11^{\circ}of Theorem 3 is proved by first showing that there exists a positive numberAaso that if at least one of the coefficients of the polynomialsπand,and=1,2,,t\pi_{i},i=1,2,\ldots,tisA\geq ain absolute value, we have

(20)(20)
A[and=1t|πand|Sand]>μ.\mathrm{A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|^{s_{i}}\right]>\mu.

This property results from the following three lemmas:
Lemma 1. Ifm<km<kandπPm\pi\in P_{m}, there is a positive numberbmb_{m}so that
(21)

|π|bm|\pi|\geq b_{m}

at leastkmk-mfrom the pointsyand,and=1,2,,ky_{i},i=1,2,\ldots,k.
It is for this
(22)

It is(z1,z2,,zm+1)\mathrm{E}\left(z_{1},z_{2},\ldots,z_{m+1}\right)

the best approximation, in Chebisev's sense, ofxnx^{n}in the polis, we are called by the degreem1m-1on distinct pointsz1,z2,,zm+1z_{1},z_{2},\ldots,z_{m+1}. The value of (22) is well known [13] but there is no need to reproduce it here. We only note that this number is positive 4 ).

Let's takebm=minIt is(yand1,yand2,,yandm+1)b_{m}=\min\mathrm{E}\left(y_{i_{1}},y_{i_{2}},\ldots,y_{i_{m+1}}\right), where the minimum refers to all combinationsand1,and2,,andm+1i_{1},i_{2},\ldots,i_{m+1}Citem+1m+1of the indices1,2,,k1,2,\ldots,k. From the properties of best approximation polynomials and the definition of the numberbmb_{m}it turns out that among the firstm+1m+1punctureyandy_{i}there is one, eitheryjy_{j}which|π|bm|\pi|\geq b_{m}Leaving aside the pointyj1y_{j1}, among the firstm+1m+1punctureyandy_{i}left, there is one, eitheryj2y_{j_{2}}, which|π|bm|\pi|\geqq b_{m}. We leave aside the pointyj2y_{j2}and we continue the process. In this way we determine the pointsyj1y_{j1},yj2,,yjm+1y_{j_{2},\ldots,}y_{j_{m+1}}on which inequality (21) is verified.

0 0 footnotetext: 3. We have Sgz=1,0sgz=-1,0resp, 1 asz<z<,==respectively.>0>0IfSsis a positive and even integer, we have|z|S1Sgz=zS1,|z|S1=zS1|z|s-1\mathrm{sg}z=z^{s-1},|z|s-1=z^{s-1}sgzz, whatever it iszz. 4. If ε=min|zandzj|\epsilon=\min\left|z_{i}-z_{j}\right|, we have
It is(z1,z2,,zm+1)>εmm+1\mathrm{E}\left(z_{1},z_{2},\ldots,z_{m+1}\right)>\frac{\epsilon^{m}}{m+1}

We notice that the numberbmb_{m}does not depend on the polynomialπ\pi.
Lemma 2. IfπPm\pi\in P_{m}and if we have|π|M|\pi|\leqq\mathrm{M}, onmmdistinct pointsz1,z2,,zmz_{1},z_{2},\ldots,z_{m}, then there is a positive numberF(z1,z2,,zm)\mathrm{F}\left(z_{1},z_{2},\ldots,z_{m}\right)so that all the coefficients of the polynomialπ\pito be<F(z1,z2,,zm)M<\mathrm{F}\left(z_{1},z_{2},\ldots,z_{m}\right)\mathrm{M}, in absolute valueAˇ\check{a}

The property is well known and we can dispense with giving its proof here 5 ).

Lemma 3. IfπandPnand,and=1,2,,t,n=n1+n2+n<k\pi_{i}\in P_{ni},i=1,2,\ldots,t,n=n_{1}+n_{2}\ldots+n<kand if

and=1t|πand|SandN,\prod_{i=1}^{t}\left|\pi_{i}\right|^{s_{i}}\leqq\mathrm{~N}, (23)

on the pointsy1,y2,,yky_{1},y_{2},\ldots,y_{k}, there is a positive numbercc, independent of the polynomialsπand,and=1,2,,t\pi_{i},i=1,2,\ldots,t, so that all the coefficients of these polynomials are<cN<c\mathrm{~N}in absolute value.

The proof can be done by complete induction onttFort=1t=1the property is true because then (23) becomes|τ1|S1N\left|\tau_{1}\right|^{\mid s_{1}}\leq N
and, based on Lemma 2, the coefficient of the polynomialπ1\pi_{1}ARE<F(y1,y2,,yn1)NS1¯<\mathrm{F}\left(y_{1},y_{2},\ldots,y_{n_{1}}\right)\mathrm{N}^{\overline{s_{1}}}in absolute value.

Suppose the property is true fort1(t>1)t-1(t>1)powers and demonstrate it forttpowers. Let us therefore assume that we have (23). Based on Lemma 1, letyj1,yj2,,yjkn1,kn1y_{j_{1}},y_{j_{2}},\ldots,y_{j_{k-n_{1}}},k-n_{1}between the pointsyandy_{i}which we have|π1|bn1\left|\pi_{1}\right|\geqq b_{n_{1}}. Thenkn1>n2+n3++ntk-n_{1}>n_{2}+n_{3}+\cdots+n_{t}and on these points we havemtand=2|πand|andSNbn1S1\underset{i=2}{\stackrel{{\scriptstyle t}}}\left|\pi_{i}\right|{}^{s}i\leqq Nb_{n_{1}}^{-s_{1}}It follows that there was a numbercc^{\prime}so that all the coefficients of the polynomialsπ2,π3,,πt\pi_{2},\pi_{3},\ldots,\pi_{t}to be<cbn1S1N<c^{\prime}b_{n_{1}}^{-s_{1}}\mathrm{~N}, in absolute value. Similarly, we prove that there is a numberc"c^{\prime\prime}so that all the coefficients of the polynomialsπ1,π3,π4,,πt\pi_{1},\pi_{3},\pi_{4},\ldots,\pi_{t}to be<c"bn2S2N<c^{\prime\prime}b_{n_{2}}^{-s_{2}}\mathrm{~N}, in absolute value. Ifc=MAX(cbn1S1,c"bn2S2)c=\max\left(c^{\prime}b_{n_{1}}^{-s_{1}},c^{\prime\prime}b_{n_{2}}^{-s_{2}}\right), we see that the property is true forttpowers. With this, Lemma 3 is proven.
11. Let us return to pt.101^{0}of theorem 3.

We have obviously

μA[|x|n1S1+n2S2++ntSt]=T.\mu\leqq\mathrm{A}\left[|x|^{\left.\mid n_{1}s_{1}+n_{2}s_{2}+\cdots+n_{t}s_{t}\right]=\mathrm{T}.}\right.

Let us now assume that

and=1t|πand|SandN1=Tmin(λ1,λ2,,λk)\prod_{i=1}^{t}\left|\pi_{i}\right|^{s_{i}}\leqq\mathrm{~N}_{1}=\frac{\mathrm{T}}{\min\left(\lambda_{1},\lambda_{2},\ldots,\lambda_{k}\right)}

and eitherc1c_{1}, the numbercccorresponding toN=N1\mathrm{N}=\mathrm{N}_{1}from Lemma 3.
If thenA=c1N1a=c_{1}N_{1}and if at least one of the coefficients of the polynomialsπ1,π2,,πt\pi_{1},\pi_{2},\ldots,\pi_{t}isA\geqq ain absolute value, we haveand=1|πand|Sand>N1\prod_{i=1}\left|\pi_{i}\right|s_{i}>\mathrm{N}_{1}on
5) Ifδ=MAXand|zand|\delta=\max_{i}\left|z_{i}\right|, we have

F(z1,z2,,zm)<MIt is(z1,z2,,zyouyou)+(δ+1)m.\mathrm{F}\left(z_{1},z_{2},\ldots,z_{m}\right)<\frac{\mathrm{M}}{\mathrm{E}\left(z_{1},z_{2},\ldots,z_{u}u\right)}+(\delta+1)^{m}.

at least one of the pointsyandy_{i}We then haveA|and=1t|πand|Sand]>Tμ\left.\left.A\left|\prod_{i=1}^{t}\right|\pi_{i}\right|^{s_{i}}\right]>T\geq\muand ines
the equality (20) is verified. equality (20) is verified.

By well-known reasoning we can now demonstrate the existence of at least one minimizing system.

From his definitionμ\muand from the previous properties, it follows that we can find an infinite series of polynomial systemsπand(m)Pnand,and=1,2,,t,m=1,2,\pi_{i}{}^{(m)}\in P_{n_{i}},i=1,2,\ldots,t,m=1,2,\ldotsso that, on the one hand, all the coefficients of these polynomials are<A<ain absolute value sí, on the other hand, if for abbreviation we putπ(m)=and=1t|πand(m)|SS\pi^{(m)}=\prod_{i=1}^{t}\left|\pi_{i}^{(m)}\right|s^{s}to have

limm𝔸[Π(m)]=μ.\lim_{m\rightarrow\infty}\mathbb{A}\left[\Pi^{(m)}\right]=\mu.

We can then extract from the string(Π(m))m=1\left(\Pi^{(m)}\right)_{m=1}^{\infty}a partial string(Π(Vm))m=1\left(\Pi^{(vm)}\right)_{m=1}^{\infty}so that

limnπand(γm)=πandPnand,and=1,2,,t\lim_{n\rightarrow\infty}\pi_{i}^{(\gamma m)}=\pi_{i}^{*}\in P_{n_{i}},\quad i=1,2,\ldots,t

uniformly in any finite interval.
It then follows that the polynomialsπand,and=1,2,,t\pi_{i}^{*},i=1,2,\ldots,tform a minimizing system.

With this for101^{0}of Theorem 3 is proved.
12. Letπand,and=1,2,,t\pi_{i},i=1,2,\ldots,ta system of minimizing polynomials. Let us also assume that1you<t1\leq u<t, whereyouuis a natural number and is a polynomialπ1,π2,,πyou\pi_{1},\pi_{2},\ldots,\pi_{u}decomposes into the product of two real factorsπand=πandπand"\pi_{i}=\pi_{i}^{\prime}\pi_{i}^{\prime\prime},and=1,2,,youi=1,2,\ldots,uso thatπandPnandπand"Pn"and\pi_{i}^{\prime}\in P_{n^{\prime}i}\pi_{i}^{\prime\prime}\in P_{n^{\prime\prime}i}so$andnand+nand"=nand\mathdollar\stackreln_{i}^{\prime}+n_{i}^{\prime\prime}=n_{i}We can assumenand>0,and=1,2,,youn_{i}^{\prime}>0,i=1,2,\ldots,uandc\mathrm{c}^{*}, in particular, the factorπand"\pi_{i}^{\prime\prime}can also be reduced to 1 (thennand"=0,nand=nandn_{i}^{\prime\prime}=0,n_{i}^{\prime}=n_{i}).

Let us consider the linear functional

A1[f]=A[(and=1n|πand"|Sand,and=you+1n|πand|Sand)f]\mathrm{A}_{1}[f]=\mathrm{A}\left[\left(\prod_{i=1}^{n}\left|\pi_{i}^{\prime\prime}\right|s_{i},\prod_{i=u+1}^{n}\left|\pi_{i}\right|s_{i}\right)f\right] (24)

This functional is of the formA1[f]=j=1kλjf(yj)A_{1}[f]=\sum_{j=1}^{k}\lambda_{j}^{\prime}f\left(y_{j}\right), where

λj=λjand=1you|πand"(yj)|Sandand=you+1n|πand(yj)|Sand,j=1,2,,k\lambda_{j}^{\prime}=\lambda_{j}\prod_{i=1}^{u}\left|\pi_{i}^{\prime\prime}\left(y_{j}\right)\right|^{s_{i}}\cdot\prod_{i=u+1}^{n}\left|\pi_{i}\left(y_{j}\right)\right|^{s_{i}},j=1,2,\ldots,k

It is seen that at mostn1"+n2"++nyou"+nyou+1nyou+2++ntn_{1}^{\prime\prime}+n_{2}^{\prime\prime}+\cdots+n_{u}^{\prime\prime}+n_{u+1}n_{u+2}+\cdots+n_{t}coefficientsλj\lambda_{j}^{\prime}is canceled and at leastk=k(n+n2++nt+n1+n2+++nyou>n1+n2++nyouk^{\prime}=k-\left(n+n_{2}+\cdots+n_{t}^{\prime}+n_{1}^{\prime}+n_{2}^{\prime}+\cdots+\right.+n_{u}^{\prime}>n_{1}^{\prime}+n_{2}^{\prime}+\cdots+n_{u}^{\prime}are positive. For101^{0}of Theorem 3 applies to the functional (24) and the system of polynomialsπ1,π2,,πyou\pi_{1}^{\prime},\pi_{2}^{\prime},\ldots,\pi_{u}^{\prime}necessarily coincides with a minimizing system for the functional (24), for the powersS1,S2,,Syous_{1},s_{2},\ldots,s_{u}to which the degrees correspond respectivelyn1,n2,,nyoun_{1}^{\prime},n_{2}^{\prime},\ldots,n_{u}^{\prime}To show this, let us assume the opposite and then letπ1,πc,,πyou\pi_{1}^{*},\pi_{c}^{*},\ldots,\pi_{u}^{*}a minimizing system corresponding to the functionalityA1[f]\mathrm{A}_{1}[f]We have
A[and=1you|πandπand"|Sandand=you+1t|πand|Sand]=A1[and=1you|πand|Sand]<A1[and=1you|πand|Sand]=A[and=1t|πand|Sand]\mathrm{A}\left[\prod_{i=1}^{u}\left|\pi_{i}^{*}\pi_{i}^{\prime\prime}\right|^{s_{i}}\prod_{i=u+1}^{t}\left|\pi_{i}\right|^{s_{i}}\right]=\mathrm{A}_{1}\left[\prod_{i=1}^{u}\left|\pi_{i}^{*}\right|^{s_{i}}\right]<\mathrm{A}_{1}\left[\prod_{i=1}^{u}\left|\pi_{i}^{\prime}\right|^{s_{i}}\right]=\mathrm{A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|^{s_{i}}\right]
so

A[and=1you|πandπand"|Sandand=you+1t|πand|Sand]<A[and=1t|πand|Sand]\mathrm{A}\left[\prod_{i=1}^{u}\left|\pi_{i}^{\star}\pi_{i}^{\prime\prime}\right|s_{i}\prod_{i=u+1}^{t}\left|\pi_{i}\right|^{s}i\right]<\mathrm{A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|s_{i}\right]

which contradicts the hypothesis thatπand,and=1,2,,t\pi_{i},i=1,2,\ldots,tit is a minimizing system.

13. Based on the above observations, to demonstrate for22^{\circ}of Theorem 3 it is enough to assumet=1t=1, which simplifies the reasoning.

Be it thenπPn\pi\in P_{n}a minimizing polynomial. Suppose thatπ\piwould not have all real roots. Then this polynomial has a real factor of the form(xA)2+b2(x-a)^{2}+b^{2}, whereb0b\neq 0.

Let's put

π(x)=[(xA)2+b2]Q(x),π1(x)=(xA)2Q(x).\pi(x)=\bigl[(x-a)^{2}+b^{2}\bigr]Q(x),\quad\pi_{1}(x)=(x-a)^{2}Q(x).

Thenπ1Pn\pi_{1}\in P_{n}and|π1||π||\pi_{1}|\leq|\pi|for anythingxxBut the strict inequality|π1|<|π||\pi_{1}|<|\pi|is checked on at least one of the pointsyandy_{i}.

Therefore,

A[|π1|1]<A[|π|1],A\bigl[|\pi_{1}|_{1}\bigr]<A\bigl[|\pi|_{1}\bigr],

which contradicts the hypothesis thatπ\piis a minimizing polynomial.

With this, for202^{0}of Theorem 3 is proven.
14. To prove pt.303^{0}of theorem 3, based on what was established in no. 12, 13, it is enough to consider only the caset=2,n1=n2=1t=2,n_{1}=n_{2}=1. If thenπ1,π2\pi_{1},\pi_{2}is a minimizing system, we must show thatπ1π2\pi_{1}\neq\pi_{2}. Let's assume the opposite, so that we would haveπ1=π2=π\pi_{1}=\pi_{2}=\piand beψ(ε)=A[|π+S2ε||S1πS1ε|S2]\psi(\epsilon)=\mathrm{A}\left[\left|\pi+s_{2}\epsilon\right|{}^{s_{1}}\left|\pi-s_{1}\epsilon\right|s_{2}\right]. Thenψ(ε)\psi(\epsilon)is a continuous function and has a continuous derivative inε(S1,S2>1)\epsilon\left(s_{1},s_{2}>1\right)
We have
(25 ) dψdε=S1S2(S1+S2)εA[|π+S2ε|S11|πS1ε|S21Sg(π+S2ε)(πS1ε)]\frac{d\psi}{d\epsilon}=-s_{1}s_{2}\left(s_{1}+s_{2}\right)\epsilon\mathrm{A}\left[\left|\pi+s_{2}\epsilon\right|^{s_{1}-1}\left|\pi-s_{1}\epsilon\right|^{s_{2}-1}sg\left(\pi+s_{2}\epsilon\right)\left(\pi-s_{1}\epsilon\right)\right]

But

Aπ+S2ε|S11|πS1ε|S21sg(π+S2ε)(πS1ε)>0,\mathrm{A}\|\pi+s_{2}\epsilon\left|s_{1}-1\right|\pi-\left.s_{1}\epsilon\right|^{s_{2}-1}\operatorname{sg}\left(\pi+s_{2}\epsilon\right)\left(\pi-s_{1}\epsilon\right)\mid>0,

for|ε||\epsilon|quite small, since the first term is a continuous function ofε\epsilonwhich forε=0\epsilon=0it reduces toA[|π|S1+S22]>0A\left[|\pi|s_{1}+s_{2}-2\right]>0From (25) it follows that sgdψdε=\frac{d\psi}{d\epsilon}=-sgε\epsilon, for|ε||\epsilon|quite small, which shows us thatΨ(ε)\Psi(\epsilon)has a relatively strict maximum forε=0\epsilon=0For|ε||\epsilon|quite small but0\neq 0we have therefore

𝔸[|π+S2|S1|πS1|S2]<𝔸[|π|+S1S2]\mathbb{A}\left[\left|\pi+s_{2}\in\right|^{s_{1}}\left|\pi-s_{1}\in\right|^{s_{2}}\right]<\mathbb{A}\left[|\pi|{}^{s_{1}}+s_{2}\right]

which contradicts the hypothesis thatπ1=π2=π\pi_{1}=\pi_{2}=\piis a minimizing system. With this we have demonstrated thatπ1π2\pi_{1}\neq\pi_{2}and therefore for303^{0}of Theorem 3.
15. Letx1,x2,,xnx_{1},x_{2},\ldots,x_{n}roots of the polynomialπ1π2πt\pi_{1}\pi_{2}\ldots\pi_{t}, then
(26)

A[and=1t|πand|Sand]=A[and=1n|xxand|Sand]\mathrm{A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|^{s}i\right]=\mathrm{A}\left[\prod_{i=1}^{n}\left|x-x_{i}\right|^{s^{\prime}i}\right]

is a continuous function ofx1,x2,,xnx_{1},x_{2},\ldots,x_{n}. Based on the previous results, at any point where the lower bound (17) is reached, we have a relative minimum of the function. If the powersS1,S2,,Sts_{1},s_{2},\ldots,s_{t}they are all>1>1, these minima are reached only for different values ​​of the variablesxandx_{i}. In this case,
however, the function (26) is differentiable and therefore, at these points, the first-order partial derivatives of the function (26) are zero. We have 6 )

xandA[and=1t|πand|Sand]=SandA[j=1jandn|xxj|Sj|xxand|Sand1sg(xxand)]=\displaystyle\frac{\partial}{\partial x_{i}}\mathrm{~A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|^{si}\right]=s_{i}^{\prime}\mathrm{A}\left[\prod_{\begin{subarray}{c}j=1\\ j\neq i\end{subarray}}^{n}\left|x-x_{j}\right|^{s^{\prime}j}\cdot\left|x-x_{i}\right|^{s^{\prime}i}-1\operatorname{sg}\left(x-x_{i}\right)\right]=
=\displaystyle= SandA[and=1t|πand|Sand1sg(and=1tπand)(xx1)(xxand1)(xxand+1)(xxn)],\displaystyle s_{i}^{\prime}\mathrm{A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|^{s}i^{-1}\cdot\operatorname{sg}\left(\prod_{i=1}^{t}\pi_{i}\right)\cdot\left(x-x_{1}\right)\ldots\left(x-x_{i-1}\right)\left(x-x_{i+1}\right)\ldots\left(x-x_{n}\right)\right],

where does it result from?404^{0}of theorem 3, noting that under the conditions here the polynomials (i1) are linearly independent.
16. The statement for505^{0}of theorem 3 constitutes, in a more complete form, a reciprocal of the previous properties in the sense that for303^{0}results from404^{0}. In fact, the property is somewhat more general and can be stated as follows:

Theorem 4.A[f]\mathrm{A}[f]being a linear functional of the form (16) with the pointsyandy_{i}distinct and with numbersλand\lambda_{i}all positive, andS1,S2,,Sands_{1},s_{2},\ldots,s_{i}a system of all powers>1>1andn1,n2,,ntn_{1},n_{2},\ldots,n_{t}a system of corresponding degrees, with their sumn<kn<k, if the polynomialsπandPnand,and=1,2,,t\pi_{i}\in P_{n_{i}},i=1,2,\ldots,tverify equality (19) for any polynomialQ(x)\mathrm{Q}(x)of the degreen1n-1, then all the roots of the polynomial (18) are real, distinct and separated from the pointsyandy_{i},and=1,2,,ki=1,2,\ldots,k.

The separation property of the statement means that if we assume the pointsyandy_{i}arranged in ascending order, soy1<y2<<yky_{1}<y_{2}<\cdots<y_{k}, we haveπ(y1)0,π(yk)0\pi\left(y_{1}\right)\neq 0,\pi\left(y_{k}\right)\neq 0and the string

π(y1),π(y2).,π(yk)\pi\left(y_{1}\right),\pi\left(y_{2}\right)\ldots.,\pi\left(y_{k}\right) (27)

presents (after removing any zero terms) exactlynnsign variations.

We first observe that the polynomialπ\pibeing of the degreenn, the string (27) represents at mostnnsign variations. In addition, ifπ(y1)=0,π(yk)0\pi\left(y_{1}\right)=0,\pi\left(y_{k}\right)\neq 0or ifπ(y1)0,π(yk)=0\pi\left(y_{1}\right)\neq 0,\pi\left(y_{k}\right)=0at mostn1n-1, and ifπ(y1)=π(yk)=0\pi\left(y_{1}\right)=\pi\left(y_{k}\right)=0at mostn2n-2variations of sign. Finally, we observe that fromn<kn<kit follows that at least one term of the sequence (27) is different from zero.

Let us now suppose that the string would only presentn<nn^{\prime}<nsign variations and either

π(yj1),π(yj2),,π(yjn+1),j1<j2<<jn+1\pi\left(y_{j_{1}}\right),\pi\left(y_{j_{2}}\right),\ldots,\pi\left(y_{j_{n^{\prime}+1}}\right),j_{1}<j_{2}<\cdots<j_{n^{\prime}+1}

a substring of (27) that represents exactlynn^{\prime}sign variations. LetandVi_{v}the smallest natural number such thatπ(yjV)π(yjV+andV)<0\pi\left(y_{jv}\right)\pi\left(y_{jv+iv}\right)<0. Thus1andVjV+1jV1\leq i_{v}\leq j_{v+1}-j_{v},V=1,2,,nv=1,2,\ldots,n^{\prime}Let's take the points.ξ1,ξ2,,ξn\xi_{1},\xi_{2},\ldots,\xi_{n^{\prime}}so thatyjV<ξV<yjV+andVy_{j_{v}}<\xi_{v}<y_{jv+i_{v}},V=1,2,,nv=1,2,\ldots,n^{\prime}and consider the polynomialQ(x)=π(yjn+1)(xξ1)xξ2)(xξn)\left.\mathrm{Q}(x)=\pi\left(y_{jn^{\prime}+1}\right)\left(x-\xi_{1}\right)x-\xi_{2}\right)\ldots\left(x-\xi_{n^{\prime}}\right)which is the degreenn1n^{\prime}\leq n-1and which does not cancel out at any of the pointsyandy_{i}From the way the points were chosenξV\xi_{v}, it follows that we haveSandπ(yj)Q(yt)0,and=1.2,,ksi\pi\left(y_{j}\right)\mathrm{Q}\left(y_{t}\right)\geq 0,i=1.2,\ldots,k, and then we have

  1. 6.

    we
    d|x|Sdx=S|x|S1\frac{d|x|s}{dx}=s|x|s-1sgx forS>1s>1andd2|x|Sdx2=S(S1)|x|S2\frac{d^{2}|x|s}{dx^{2}}=s(s-1)|x|s-2forS2s\geq 2

A[and=1t|πand|Sand1(sgπ)Q]=A[and=1t|πand|Sand1|Q|]>0\mathrm{A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|^{s_{i}-1}(\mathrm{sg}\pi)\mathrm{Q}\right]=\mathrm{A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|s_{i}-1|\mathrm{Q}|\right]>0

in contradiction with equality (19).
With this, Theorem 4 is proven. It generalizes a property established long ago fort=1,S1=2t=1,s_{1}=2[7].
17. Theorem 3 and the previous results show us that the minimum problem treated always returns to the particular case when the degrees corresponding to the powers are all equal to 1.

If for abbreviation we writewithμS1,S2,,St\mathrm{cu}\mu^{s_{1}},s_{2},\ldots,s_{t}the number (17) when all the degreesn1,n2,,ntn_{1},n_{2},\ldots,n_{t}are equal to 1, we have

μS1,S2,,St<μS2+S2,S3,S4,,St(t>3)\mu^{s_{1}},s_{2},\ldots,s_{t}<\mu^{s_{2}}+s_{2},s_{3},s_{4},\ldots,s_{t}\quad(t>3)

and in particularμS1,S2,,St<μS1+S2++St\mu^{s_{1}},s_{2},\ldots,s_{t}<\mu^{s_{1}}+s_{2}+\cdots+s_{t}
In particular cases t=1t=1andS1=S2==Sts_{1}=s_{2}=\cdots=s_{t}are equivalent in the previous sense.

To specify the uniqueness of the minimizing system, we must say that two minimizing polynomial systems in which, for each group of equal powers, the product of the polynomials is not considered distinct are not considered.πand\pi_{i}is the same.

It is known that the minimizing system is unique if the powers are equal and>1>1[3].

If the powersSands_{i}are not all equal, uniqueness no longer occurs in general, as will result from the examples in$4\mathdollar 4These examples also show us that the property expressed by pt.404^{0}of Theorem 3 does not characterize minimizing systems, in other words, there are also polynomials (18) for which the orthogonality property404^{0}of Theorem 3 is verified but which are not formed with a minimizing system.

If the powersS1,S2,,Sts_{1},s_{2},\ldots,s_{t}they are all2\geq 2, so andSand2,and=1,2,,ns_{i}^{\prime}\geq 2,i=1,2,\ldots,nwe can easily prove that any system of polynomialsπand,and=1,2,,t\pi_{i},i=1,2,\ldots,twhich verifies the property44^{\circ}of Theorem 3 , so in particular the minimizing polynomials, correspond to strict relative minima of the function (26). Indeed, in this case, the function (26) also has second-order partial derivatives and we have

2xand2A[and=1t|πand|Sand]=Sand(Sand1)A[(j=1jandn|xxand|Sj)|xxand|Sand2]>0and=1,2.,n\begin{gathered}\frac{\partial^{2}}{\partial x_{i}^{2}}\mathrm{~A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|^{s_{i}}\right]=s_{i}^{\prime}\left(s_{i}^{\prime}-1\right)\mathrm{A}\left[\left(\prod_{\begin{subarray}{c}j=1\\ j\neq i\end{subarray}}^{n}\left|x-x_{i}\right|^{s^{\prime}j}\right)\left|x-x_{i}\right|^{s^{\prime}i-2}\right]>0\\ i=1,2\ldots.,n\end{gathered}

(28)

2xandxjA[and=1t|πand|Sand]=\frac{\partial^{2}}{\partial x_{i}\partial x_{j}}\mathrm{~A}\left[\prod_{i=1}^{t}\left|\pi_{i}\right|^{s_{i}}\right]=

=SandSjA[(V=1Vand,Vjn|xxV|SV)|xxand|Sand1|xxj|Sj1sg(xxand)(xxj)]==s_{i}^{\prime}s_{j}^{\prime}\mathrm{A}\left[\left(\prod_{\begin{subarray}{c}v=1\\ v\neq i,v\neq j\end{subarray}}^{n}\left|x-x_{v}\right|^{s^{\prime}v}\right)\left|x-x_{i}\right|^{s^{\prime}{}_{i}-1}\left|x-x_{j}\right|^{s^{\prime}j^{-1}}\mathrm{sg}\left(x-x_{i}\right)\left(x-x_{j}\right)\right]=
=SandSjA(and=1t|πand|Sand1)(sgand=1tπand)(xx1)(xxand1)(xxand+1)(xxj1)=s_{i}^{\prime}s_{j}^{\prime}\mathrm{A}\mid\left(\prod_{i=1}^{t}\left|\pi_{i}\right|^{s_{i}-1}\right)\left(\mathrm{sg}\prod_{i=1}^{t}\pi_{i}\right)\left(x-x_{1}\right)\ldots\left(x-x_{i-1}\right)\left(x-x_{i+1}\right)\ldots\left(x-xj_{-1}\right).
(xxj+1)(xxn)]and=1,2,,j1,j=2,3,,n\left.\cdot\left(x-x_{j+1}\right)\ldots\left(x-x_{n}\right)\right]\quad i=1,2,\ldots,j-1,j=2,3,\ldots,n.

Based on the orthogonality property (19) at the points considered, all derivatives (28) are zero and the stated property results. In the case whenn>kn>k, it is obvious that, in the above sense, there are an infinity of minimizing systems.

Ifn=kn=k, uniqueness does not occur, in the above sense, unless the powers are equal.

It is noteworthy that for44^{\circ}of Theorem 3 also holds in the casenkn\geq k, even with its reciprocal, in the sense that, if (19) holds for any polynomialQ(x)Q(x)of the degreen1n-1, the systemπand,and=1,2,,t\pi_{i},i=1,2,\ldots,tit is minimizing.

Indeed, let us assume that the system would not be minimizing, so that the polynomial (18) would not vanish at all pointsyandy_{i}Either, for fixing ideas,π(yk)0\pi(y_{k})\neq 0We then have…

A[and=1t|πand|Sand1(sgπ)(xy1)(xy2)(xyk1)]==λkand=1t|πand(yk)|Sand1sg(π(yk))(yky1)(yky2)(ykyk1)0\begin{gathered}A\left[\prod_{i=1}^{t}\left|\pi_{i}\right|^{s_{i}-1}\cdot(\operatorname{sg}\pi)\left(x-y_{1}\right)\left(x-y_{2}\right)\ldots\left(x-y_{k-1}\right)\right]=\\ =\lambda_{k}\prod_{i=1}^{t}\left|\pi_{i}\left(y_{k}\right)\right|^{s_{i}-1}\operatorname{sg}\left(\pi\left(y_{k}\right)\right)\cdot\left(y_{k}-y_{1}\right)\left(y_{k}-y_{2}\right)\ldots\left(y_{k}-y_{k-1}\right)\neq 0\end{gathered}

which contradicts equality (19).

§ 3. The existence of Gauss-type formulas.

  1. 18.

    We first have the following:

Theorem 5. For any linear functionalA[f]\mathrm{A}[f]of the form (16), with the pointsyandy_{i}distinct and with numbersλand\lambda_{i}all positive, relative to any natural numbern<kn<kand to any multiplicity order systemR1,R2,.Rnr_{1},r_{2},\ldots.r_{n}consisting of odd numbers, there is at least one Gaussian formula (8).

The hypothesis that the numbersR1,R2,,Rnr_{1},r_{2},\ldots,r_{n}that they are all odd is essential, as can easily be seen from considerations analogous to those made in No. 7.

Theorem 5 follows from Theorem 3. To see this, it is enough to taket=nt=n, the powersS1,S2,,Sn(2)s_{1},s_{2},\ldots,s_{n}(\geq 2)respectively equal toR1+1,R2+1,,Rn+1r_{1}+1,r_{2}+1,\ldots,r_{n}+1and the corresponding degrees all equal to 1. If (2) are the roots of the polynomial (18) corresponding to a minimizing system, we have

it(x)=and=1n|xxand|Sand1sg(xx1)(xx2)(xxn)l(x)=\prod_{i=1}^{n}\left|x-x_{i}\right|^{s_{i}-1}\operatorname{sg}\left(x-x_{1}\right)\left(x-x_{2}\right)\ldots\left(x-x_{n}\right)

and condition (19) reduces to the orthogonality of the polynomialit(x)l(x)with any polynomial of degreen1n-1.

In this way, each minimizing system corresponds to a Gaussian formula (8).
19. Let us return to a linear functional𝔸[f]\mathbb{A}[f]in order of positivitykk, as it was defined in no. 7. Obviously, if𝔸[f]\mathbb{A}[f]has the order of positivitykk, it also has the order of positivitykk^{\prime}for anythingk<kk^{\prime}<k.

If we introduce the moments

αand=A[xand],and=0,1,\alpha_{i}=\mathrm{A}\left[x^{i}\right],i=0,1,\ldots (29)

and Hankel determinants

Aj=αλ+μλ,μ=0,1,,j,j=0,1,A_{j}=\left\|\alpha_{\lambda+\mu}\right\|_{\lambda,\mu=0,1,\ldots,j},\quad j=0,1,\ldots

appropriate, the necessary and sufficient condition that the functionalA[f]A[f]to have the order of positivitykkit is so that we have…

Aj>0,j=0,1,,k1,A_{j}>0,j=0,1,\ldots,k-1,

or, which is equivalent, as the quadratic formand,j=0k1αand+jξandξj\sum_{i,j=0}^{k-1}\alpha_{i+j}\xi_{i}\xi_{j}to be definite and positive.

It is clear that in general the determination of Gaussian-type formulas depends only on the mutual ratios of the firstp+n+1p+n+1MOMENTSαand\alpha_{i},and=0,1,,p+ni=0,1,\ldots,p+nhis/hersA[f]\mathrm{A}[f]More precisely, in determining Gaussian formulas, one can abstract from a linear transformation of the variablexxand by a non-zero constant factor of the functionalA[f]\mathrm{A}[f]. Otherwise, this observation is generally valid for formulas of the form (3) which, through the indicated transforms s, retain their form and degree of accuracy.

IfA[f]\mathrm{A}[f]is a linear functional of positive orderkk, there is a polynomialρkPk\rho_{k}\in P_{k}and only one that is orthogonal to any polynomial of degreek1k-1This is the orthogonal polynomial of degreekkattached functionallyA[f]7\mathrm{A}[f]^{7}).

The polynomialρk\rho_{k}has all real and distinct roots. Indeed, otherwise, this polynomial should have a divisor of the form(xA)2+b2(x-a)^{2}+b^{2}(A,ba,breal). If thenpk=[(xA)2+b2]Qp_{k}=\left[(x-a)^{2}+b^{2}\right]\mathrm{Q}, Q is a polynomial of degreek2k-2and we haveOh[ρkQ]=A[((xA)Q)2]+b2A[Q2]>0\AA \left[\rho_{k}\mathrm{Q}\right]=\mathrm{A}\left[((x-a)\mathrm{Q})^{2}\right]+b^{2}\mathrm{~A}\left[\mathrm{Q}^{2}\right]>0, which contradicts the orthogonality property.

It is clear that the polynomialρk\rho_{k}depends only on the first2k2kMOMENTSαand,and=0,1,,2k1\alpha_{i},i=0,1,\ldots,2k-1his/hersA[f]\mathrm{A}[f].

In particular, a functional of the form (16), in whichyandy_{i}are distinct andλand\lambda_{i}positive, has the degree of positivitykk. In this caseρk(x)=(xy1)(xy2)(xyk)\rho_{k}(x)=\left(x-y_{1}\right)\left(x-y_{2}\right)\ldots\left(x-y_{k}\right).

But we also have a kind of reciprocal of this property in the following sense. Ifyand,and=1,2,,ky_{i},i=1,2,\ldots,kare the roots of the polynomialρk\rho_{k}and the moments (29) verify the inequalities (30), there is a linear functionalOh(k)[f]\AA ^{(k)}[f]of the form (16), with all coefficientsλand\lambda_{i}positive and so that we have

αand=A(t)[xand],and=0,1,,2k1.\alpha_{i}=\mathrm{A}^{(t)}\left[x^{i}\right],i=0,1,\ldots,2k-1. (31)

Indeed, any polynomialQ(x)\mathrm{Q}(x)of the degreekkis of the form

Q(x)=Aρk(x)+and=1kρk(x)(xyand)ρ(yand)kQyand)\left.\mathrm{Q}(x)=a\rho_{k}(x)+\sum_{i=1}^{k}\frac{\rho_{k}(x)}{\left(x-y_{i}\right)\rho^{\prime}{}_{k}\left(y_{i}\right)}\mathrm{Q}y_{i}\right)

where the constantAais equal to zero if and only ifQ(x)\mathrm{Q}(x)is of the degreek1k-1IfP(x)\mathrm{P}(x)is a polynomial of degree2k12k-1, he is still==one the product of two (real) polynomialsQ(x),Q1(x)\mathrm{Q}(x),\mathrm{Q}_{1}(x), first degreekkand the second in rankk1k-1So we have

P(x)=Q(x)Q1(x)=[Aρk(x)+and=1kρk(x)(xyand)ρk(yand)Q(yand)]\displaystyle\mathrm{P}(x)=\mathrm{Q}(x)\mathrm{Q}_{1}(x)=\left[a\rho_{k}(x)+\sum_{i=1}^{k}\frac{\rho_{k}(x)}{\left(x-y_{i}\right)\rho_{k}^{\prime}\left(y_{i}\right)}\mathrm{Q}\left(y_{i}\right)\right]
[and=1kρk(z)(xyand)ρk(yand)Q1(yand)]=Aand=1kρk(x)ρk(x)(xyand)ρk(yand)Q1(yand)+\displaystyle{\left[\sum_{i=1}^{k}\frac{\rho_{k}(z)}{\left(x-y_{i}\right)\rho_{k}^{\prime}\left(y_{i}\right)}\mathrm{Q}_{1}\left(y_{i}\right)\right]=a\sum_{i=1}^{k}\rho_{k}(x)\frac{\rho_{k}(x)}{\left(x-y_{i}\right)\rho_{k}^{\prime}\left(y_{i}\right)}\mathrm{Q}_{1}\left(y_{i}\right)+}
0 0 footnotetext: 7. Existence and uniqueness of the polynomial ρk\rho_{k}results only fromAk10A_{k-1}\neq 0IfAk1=0A_{k-1}=0such a polynomial either does not exist or it is not uniquely determined.
+and=1k(ρk(x)(xyand)ρkyand)2P(yand)++and,j=1andjkρk(x)ρk(x)(xyand)(xyj)ρk(yand)ρk(yj)Q(yand)Q1(yj)\begin{gathered}+\sum_{i=1}^{k}\left(\frac{\rho_{k}(x)}{\left(x-y_{i}\right)\rho_{k}^{\prime}y_{i}}\right)^{2}\mathrm{P}\left(y_{i}\right)+\\ +\sum_{\begin{subarray}{c}i,j=1\\ i\neq j\end{subarray}}^{k}\rho_{k}(x)\frac{\rho_{k}(x)}{\left(x-y_{i}\right)\left(x-y_{j}\right)\rho_{k}^{\prime}\left(y_{i}\right)\rho_{k}^{\prime}\left(y_{j}\right)}\mathrm{Q}\left(y_{i}\right)\mathrm{Q}_{1}\left(y_{j}\right)\end{gathered}

If we take orthogonality into account, we deduce

A[P]=A(k)[P]=and=1andbλandP(yand)\mathrm{A}[\mathrm{P}]=\mathrm{A}^{(k)}[\mathrm{P}]=\sum_{i=1}^{ib}\lambda_{i}\mathrm{P}\left(y_{i}\right)

where

λand=A[(ρk(x)(xyand)ρk(yand))2]>0,and=1,2,,k\lambda_{i}=\mathrm{A}\left[\left(\frac{\rho_{k}(x)}{\left(x-y_{i}\right)\rho_{k}^{\prime}\left(y_{i}\right)}\right)^{2}\right]>0,i=1,2,\ldots,k

From this, in particular, the formulas (31) result.
From the previous analysis we recall:
Lemma 4. – If the linear functionalA[f]\mathrm{A}[f]has positivity ordinalkk, a linear functional can be foundA(k)[f]\mathrm{A}^{(k)}[f]of the form (16) with all coefficientsλand\lambda_{i}positive and so that we haveA[f]=A(it)[f]\mathrm{A}[f]=\mathrm{A}^{(l)}[f], for any polynomial of degree2k12k-1.

It is easy to see that the linear functionalOh(k)[f]\AA ^{(k)}[f]is unique and is precisely the one determined above.
20. - We assume of course that ifA[f]\mathrm{A}[f]has the order of positís vitakk, intervalu1 I contains the roots of the orthogonal polynomial of degreekkattached to this functionality.

We can observe that ifAais a number such thatA[(xA)Q2]>0\mathrm{A}\left[(x-a)\mathrm{Q}^{2}\right]>0, for any polynomial Q of degreek1k-1, the roots of the polynomialρk\rho_{k}they are all>A>aIndeed, ifρk\rho_{k}would have a rootx0Ax_{0}\leq a, then if we putρk(x)=(xx0)Q(x)\rho_{k}(x)=\left(x-x_{0}\right)\mathrm{Q}(x), we would haveA[ρkQ]=A[(xA)Q2]++(αx0)A[Q2]>0\mathrm{A}\left[\rho_{k}\mathrm{Q}\right]=\mathrm{A}\left[(x-a)\mathrm{Q}^{2}\right]++\left(\alpha-x_{0}\right)\mathrm{A}\left[\mathrm{Q}^{2}\right]>0, which centers the orthogonality. It is also seen that ifA[(xb)Q2]<0\mathrm{A}\left[(x-b)\mathrm{Q}^{2}\right]<0, for any orthogonal Q , its nileρk\rho_{k}they are all<b<b.

Thus, for example, the classical functionals of Jacobí, Laguerre and Hermite

and(α,β)[f]=1+1(1x)α(1+x)βf(x)𝑑x(α,β>1)I^{(\alpha,\beta)}[f]=\int_{-1}^{+1}(1-x)^{\alpha}(1+x)^{\beta}f(x)dx\quad(\alpha,\beta>-1) (32)
IT(α)[f]=0it isxxαf(x)𝑑x(α>1)L^{(\alpha)}[f]=\int_{0}^{\infty}e^{-x}x^{\alpha}f(x)dx\quad(\alpha>-1) (33)
H[f]=+it isx2f(x)𝑑xH[f]=\int_{-\infty}^{+\infty}e^{-x^{2}}f(x)dx (34)

have the order of positivitykk, for anythingkkIn the first case the roots of the orthogonal polynomials are in the interval (1,1-1,1), in the second in the interval(0,)(0,\infty)and in the third case in the interval(,+)(-\infty,+\infty). In these cases it is
therefore sufficient to assume that I coincides with these intervals respectively.
In the case of positive linear functionals it follows in particular that the roots of orthogonal polynomials are inside the interval I .
21. - Returning to our problem, we can now prove

Theorem 6. - For any linear functional A [ff] which has the order of positivitykk, relative to any natural number and any multiplicity order systemR1,R2,,Rnr_{1},r_{2},\ldots,r_{n}consists only of odd numbers so thatk12(p+n+1)k\geq\frac{1}{2}(p+n+1), there is at least one Gaussian formula (8).

To prove this theorem it is enough to consider the functionalOh(k)[f]\AA ^{(k)}[f]determined by Lemma 4. Let (2) then be the nodes of a Gauss formula relative to the functionalOh(t)[f]A\AA ^{(t)}[f]\quad\mathrm{O}such a formula exists becauseknk\geq nThe polynomialit(x)l(x)is orthogonal to any polynomialQ(x)\mathrm{Q}(x)of the degreen1n-1functional sideOh(k)[f]\AA ^{(k)}[f]But the productit(x)Q(x)l(x)\mathrm{Q}(x)is of the degreep+n2k1p+n\leq 2k-1, soA[itQ]=A(itk)[itQ]=0\mathrm{A}[l\mathrm{Q}]=\mathrm{A}^{(lk)}[l\mathrm{Q}]=0The polynomialit(x)l(x)is therefore orthogonal to any polynomial of degreen1n-1compared to the functionalA[f]A[f], which, based on Theorem 2, proves the property.

It is also seen that all Gauss-type formulas related to the functionalA[f]A[f]are obtained in this way.

Equalityk=nk=nit is only possible ifR1=R2==Rn=1r_{1}=r_{2}=\cdots=r_{n}=1Then the Gaussian formula is unique and has as nodes precisely the roots of the polynomialρk\rho_{k}orthogonally attached functionallyA[f]\mathrm{A}[f]Apart from this particular case, in the indicated hypotheses, the nodes of any Gaussian-type formula are separated by the roots of the orthogonal polynomialρk\rho_{k}functional attachmentA[f]\mathrm{A}[f].

We observe that, under the conditions of Theorem 6, we have
(35)

A[and=1n(xxand)Rand+1]==A(k)[and=1n(xxand)Rand+1]+A[xp+n+1]A(k)[xp+n+1].\begin{gathered}\mathrm{A}\left[\prod_{i=1}^{n}\left(x-x_{i}\right)^{r_{i}+1}\right]=\\ =\mathrm{A}^{(k)}\left[\prod_{i=1}^{n}\left(x-x_{i}\right)^{r_{i}+1}\right]+\mathrm{A}\left[x^{p+n+1}\right]-\mathrm{A}^{(k)}\left[x^{p+n+1}\right].\end{gathered}

It is seen that the minimum problem can be posed and solved on expression (35) from$2\mathdollar 2, as in the case of functionals of the form (16). This is of course the minimum problem corresponding to the powersRand+1,and=1,2,,nr_{i}+1,i=1,2,\ldots,nand the respective degrees all equal to1.81.{}^{8}) The problem reduces, on the basis of equality (35), to a corresponding problem for a functional of the form (16). There are therefore, in particular, Gauss-type formulas corresponding to the minimizing system of this problem.

IfA˙[f]\dot{A}[f]has the order of positivityk>12(p+n+1)k>\frac{1}{2}(p+n+1), formula (35) can be replaced by

A[and=1n(xxand)Rand+1]=A(k)[and=1n(xxand)Rand+1]\mathrm{A}\left[\prod_{i=1}^{n}\left(x-x_{i}\right)^{r_{i}+1}\right]=\mathrm{A}^{(k)}\left[\prod_{i=1}^{n}\left(x-x_{i}\right)^{r_{i}+1}\right]

In particular, a positive lymph functional has a positivity orderkkfor anythingkkand we therefore deduce

0 0 footnotetext: 8. The lower bound of this expression is no longer necessarily 0\geq 0.

Corollary 1. - For any positive linear functional, relative to any natural number n and to any multiplicity order systemR1,R2,,Rnr_{1},r_{2},\ldots,r_{n}consisting only of odd numbers, there is at least one Gauss-type formula.

In this case the nodes of a Gaussian formula are inside the interval I and are separated by the roots of any orthogonal polynomial of degreek>12(p+n+1)k>\frac{1}{2}(p+n+1)attached to functionality.

In particular (32), (33), (34) are functionals of this kind.
The existence of Gaussian-type formulas for the caseR1=R2==Rn=1Ms.r_{1}=r_{2}=\ldots=r_{n}=1\mathrm{~ms}par was proved by P. Turán [12].
22.- For the restR[f]\mathrm{R}[f]of the Gauss-type formulas (8) , we have

R[xp+n+1]=R[it(x)and=1n(xxand)]=A[and=1n(xxand)Rand+1]\mathrm{R}\left[x^{p+n+1}\right]=\mathrm{R}\left[l(x)\prod_{i=1}^{n}\left(x-x_{i}\right)\right]=\mathrm{A}\left[\prod_{i=1}^{n}\left(x-x_{i}\right)^{r_{i}+1}\right] (36)

It follows that under the conditions of the theorem6,R[xp+n+1]6,\mathrm{R}\left[x^{p+n+1}\right]is the smallest, for and only for Gaussian formulas that come from minimizing systems.

If the functionalA[f]A[f]has the order of positivityk>12(p+n+1)k>\frac{1}{2}(p+n+1), all Gaussian formulas have the degree of accuracy equal top+np+n, because, in this case from (36) it follows thatR[xp+n+1]>0\mathrm{R}\left[x^{p+n+1}\right]>0.

Based on an observation made in No. 5 and on the well-known expression of the remainder of the Lagrange-Hermite interpolation formula, the remainderR[f]\mathrm{R}[f]of a Gaussian type formula (8) can be written
(37)

R[f]=\mathrm{R}[f]=

=A[and=1n(xxand)Rand+1[x1,x1,,R1+1x1,x2,x2,,x2R2+1,,xn,xn,,xnRn+1,x;f]]=\mathrm{A}[\prod_{i=1}^{n}\left(x-x_{i}\right)^{r_{i}+1}[\underbrace{x_{1},x_{1},\ldots,}_{r_{1}+1}x_{1},\underbrace{x_{2},x_{2},\ldots,x_{2}}_{r_{2}+1},\ldots,\underbrace{x_{n},x_{n},\ldots,x_{n}}_{r_{n}+1},x;f]]
using a convenient notation of divided differences, which can be defined as follows:
Fie
(38)

V(f1,f2,,fm+1z1,z2,,zm+1)=fj(zand)and,j=1,2,,n+1\mathrm{V}\binom{f_{1},f_{2},\ldots,f_{m+1}}{z_{1},z_{2},\ldots,z_{m+1}}=\left\|f_{j}\left(z_{i}\right)\right\|_{i,j=1,2,\ldots,n+1}

determinant of function valuesfand=fand(x),and=1,2,,m+1f_{i}=f_{i}(x),i=1,2,\ldots,m+1on the pointsz1,z2,,zm+1z_{1},z_{2},\ldots,z_{m+1}(andiis the line index andjjof columns), provided that if a group ofVvpuncturezandz_{i}are confused, thoseV(>1)v(>1)The corresponding lines contain the values ​​of the functions and their primes.V1v-1derivatives on this point. In particular

V(z1,z2,,zm+)=V(1,x,,xmz1,z2,,zm+1)\mathrm{V}\left(z_{1},z_{2},\ldots,z_{m+}\right)=\mathrm{V}\binom{1,x,\ldots,x^{m}}{z_{1},z_{2},\ldots,z_{m+1}}

: is the Vandermonde determinant of the numbersz1,z2,,zm+1z_{1},z_{2},\ldots,z_{m+1}and

[z1,z2,zm+1;f]=V(1,x1,,xm1,fz1,z2,,zm+1)V(z1,z2,,zm+1)\left[z_{1},z_{2},\ldots z_{m+1};f\right]=\frac{\mathrm{V}\binom{1,x_{1},\ldots,x^{m-1},f}{z_{1},z_{2},\ldots,z_{m+1}}}{\mathrm{~V}\left(z_{1},z_{2},\ldots,z_{m+1}\right)}

is the divided difference (of the ordermm) of the functionf(x)f(x)on the nodesz1,z2,z_{1},z_{2},\ldots,zm+1z_{m+1}.

In the important case for applications, whenA[f]\mathrm{A}[f]is a positive linear functional, from (37) it follows that we haveR[f]>0\mathrm{R}[f]>0iff(x)f(x)is a convex function of orderp+np+n. It is then known that we have [8],
(39)

R[f]=R[xp+n+1]Dp+n+1[f]\mathrm{R}[f]=\mathrm{R}\left[x^{p+n+1}\right]\mathrm{D}_{p+n+1}[f]

where, for brevity, we denote byDm[f]\mathrm{D}_{m}[f]a difference divided by ordermmof the functionf(x)f(x)onm+1m+1convenient distinct nodes (depending on the functionff) from within the interval I. These nodes can be chosen arbitrarily close to each other[𝟏,𝟏𝟎][\mathbf{1,10}].

If, in particular, the functionf(x)f(x)admits a derivative of the orderp+n+1p+n+1, we have

R[f]=R[xp+n+1](p+n+1)!f(p+n+1)(ξ),\mathrm{R}[f]=\frac{\mathrm{R}\left[x^{p+n+1}\right]}{(p+n+1)!}f^{(p+n+1)}(\xi), (40)

whereξ\xiis a convenient point inside the intervalandIAlso, iff(x)f(x)has a derivative of orderp+np+nwhich verifies an ordinary Lipschitz condition with the constantMM, we have (41)

|R[f]||R[xp+n+1]|(p+n+1)!M|\mathrm{R}[f]|\leqq\frac{\left|\mathrm{R}\left[x^{p+n+1}\right]\right|}{(p+n+1)!}\mathrm{M}

In formulas (39), (40), (41), the coefficientR[xp+n+1]\mathrm{R}\left[x^{p+n+1}\right]can be replaced with its value extracted from (36),

Finally, in this case we observe that in the sense of delimiting the remainder, the most precise Gauss-type formulas are those that correspond to minimizing systems.

In the particular case whenk=nk=nand the functional is of the form (16), the Gaussian (unique) formula reduces to the trivial formula

A[f]and=1nλandf(yand)\mathrm{A}[f]\approx\sum_{i=1}^{n}\lambda_{i}f\left(y_{i}\right)

with the remainder identically zero. The degree of accuracy of this formula is++\infty.

§ 4. - Determination of some Gauss-type formulas

  1. 23.
    • To standardize notations, when it comes to a function: nalà lymphaF[f]\mathrm{F}[f], we denote with the corresponding lowercase Greek letter, affected by the indices, the momentsφand=F[xand],and=0,1,,with\varphi_{i}=\mathrm{F}\left[x^{i}\right],i=0,1,\ldots,\mathrm{cu}corresponding uppercase Greek letter, affected by indices, Hankel's determinantsφand,j=φand+λ+μλ,=10,1,,j,and,j=0,1,\phi_{i,j}=\left\|\varphi_{i+\lambda+\mu}\right\|\lambda,{}^{1}=0,1,\ldots,j,i,j=0,1,\ldotsand, in particular,Φj=Φa,jj=0,1,\Phi_{j}=\Phi_{o,j}j=0,1,\ldots

We also introduce the transformed momentsφand,j(ξ)=F[(ξx)andxj]\varphi_{i,j}(\xi)=\mathrm{F}\left[(\xi-x)^{i}x^{j}\right],and,j=0,1,i,j=0,1,\ldots, whereξ\xiis an independent parameter ofxxWe then have

φand,j(ξ)=V=0and(1)V(andV)ξandVφV+j\varphi_{i,j}(\xi)=\sum_{v=0}^{i}(-1)^{v}\binom{i}{v}\xi^{i-v}\varphi_{v+j}

Hankel's determinantφand,j(ξ)=φand+λ+μμ,0(ξ)λ,μA0,1,j\phi_{i,j}(\xi)=\left\|\varphi_{i+\lambda+\mu^{\mu},0}(\xi)\right\|_{\lambda,\mu_{a-0,1}\ldots,j}which is a polis nom inξ\xi, can be brought, by elementary transformations of lines and columns, to other remarkable forms. Thus, we have (0Rand0\leq r\leq i).
where we put

δandV(ξ)=(andV)ξandV,V,and=0,1,,(δandV(ξ)=0, for and<V)\delta_{i}^{v}(\xi)=\binom{i}{v}\xi^{i-v},v,i=0,1,\ldots,\left(\delta_{i}^{v}(\xi)=0,\text{ pentru }i<v\right)

In particular, forR=andr=iformula (42) becomes

Φand,j(ξ)=|φ0φ1φ1φ2φand+jφand+j+1φjφj+1φand+2jδ00(ξ)δ10(ξ)δand+j0(ξ)δ01(ξ)δ11(ξ)δand+j1(ξ)δ0t1(ξ)δ1and1(ξ1δand+jand1(ξ)|\Phi_{i,j}(\xi)=\left|\begin{array}[]{lllll}\varphi_{0}&\varphi_{1}&\ldots&\ldots&\ldots\\ \varphi_{1}&\varphi_{2}&\ldots&\ldots&\varphi_{i+j}\\ \cdots&\cdots&\cdots&\varphi_{i+j+1}\\ \varphi_{j}&\varphi_{j+1}&&&\varphi_{i+2j}\\ \delta_{0}^{0}(\xi)&\delta_{1}^{0}(\xi)&\ldots&\ldots&\delta_{i+j}^{0}(\xi)\\ \delta_{0}^{1}(\xi)&\delta_{1}^{1}(\xi)&\ldots&\ldots&\delta_{i+j}^{1}(\xi)\\ \cdots&\cdots&\cdots&\cdots&\cdots\\ \delta_{0}^{t-1}(\xi)&\delta_{1}^{i-1}\left(\xi_{1}\right.&\ldots&\delta_{i+j}^{i-1}(\xi)\end{array}\right|

The formula is also valid forR=0r=0in the following form:

Φand,j(ξ)=|φand,j(ξ)φand,1(ξ)φand,j(ξ)φand,1(ξ)φand,2(ξ)φand,j+1(ξ)φand,j(ξ)φand,j+1(ξ)φand,2j(ξ)|\Phi_{i,j}(\xi)=\left|\begin{array}[]{ccccc}\varphi_{i,j}(\xi)&\varphi_{i,1}(\xi)&\ldots&\varphi_{i,j}(\xi)\\ \varphi_{i,1}(\xi)&\varphi_{i,2}(\xi)&\ldots&\varphi_{i,j+1}(\xi)\\ \cdots&\cdots&\cdots&\cdots&\cdots\\ \varphi_{i,j}(\xi)&\varphi_{i,j+1}(\xi)&\ldots&\varphi_{i,2j}(\xi)\end{array}\right|

If we apply a well-known transformation formula to this latter determinant (see e.g., [11]), we deduce:

φand,j(ξ)=1(j+1)!Ft1,t2,,tj+1[V=1j+1(ξtV)and)V2(t1,t2,,tj+1)]\left.\phi_{i,j}(\xi)=\frac{1}{(j+1)!}\mathrm{F}_{t_{1},t_{2},\cdots,t_{j+1}}\left[\prod_{V=1}^{j+1}\left(\xi-t_{V}\right)^{i}\right)V^{2}\left(t_{1},t_{2},\ldots,t_{j+1}\right)\right]

where we notewithFt1,t2,,tj+1\mathrm{cu}\mathrm{F}_{t_{1},t_{2}},\cdots,t_{j+1}successive application of the operator F to the variablest1,t2,,tj+1t_{1},t_{2},\ldots,t_{j+1}.
24. - To effectively determine all Gauss-type formulas, it is sufficient to solve the system obtained from (10) if we replace Q withnnpolynomials of degreen1n-1linearly independent. This comes down to solving system (14) or an equivalent system in the sense of no. 6.

We will deal in particular with the case whenn>1,1you<nn>1,1\leqq u<nandRyou+1=Ryou+2===Rn=1r_{u+1}=r_{u+2}==\ldots=r_{n}=1The other orders of multiplicityR1,R2,.Ryour_{1},r_{2},\ldots.r_{u}are odd numbers (some or all of them can also be equal to 1).

Assuming we have obtained the nodesx1,x2,,xyoux_{1},x_{2},\ldots,x_{u}, the other nodes are uniquely determined as the roots of the orthogonal polynomialρnyou\rho_{n-u}very sadnyoun-uattached to functionality 9 )

C[f]=A[(and=1n(xandxRand+1)f]\mathrm{C}[\mathrm{f}]=\mathrm{A}\left[\left(\prod_{i=1}^{n}\left(x_{i}-xr_{i}+1\right)f\right]\right.

polynomialρnyou\rho_{n-u}can be obtained by solving a linear system. If we put

ρnn=(1)nyouV=0nyoudV(ξ)(ξx)V(dnn(ξ)=1)\rho_{n-n}=(-1)^{n-u}\sum_{v=0}^{n-u}d_{v}(\xi)(\xi-x)^{v}\quad\left(d_{n-n}(\xi)=1\right)

HAVE

V=0nyoudV(ξ)Rand+V,0(ξ)=0,and=0,1,nyou1\sum_{v=0}^{n-u}d_{v}(\xi)r_{i+v,0}(\xi)=0,\quad i=0,1\ldots,n-u-1 (43)

functionA[f]\mathrm{A}[f], having the order of positivityk12(p+n+1)k\geq\frac{1}{2}(p+n+1), functionalC[f]\mathrm{C}[f]will have the order of positivityk12(R1+R2++Ryou+you)nyouk-\frac{1}{2}\left(r_{1}+r_{2}+\ldots+r_{u}+u\right)\geq n-uand therefore the system (43) has a unique well-determined solution. The determinant of this system is independent ofξ\xiand is equal toΓnyou\Gamma_{n-u}.
25.- System (43) is equivalent to the lastnyoun-uequations (14). Taking into account this system, the firstyouuWe will replace equations (14) with others that will contain only the variablesx1,x2,,xyoux_{1},x_{2},\ldots,x_{u}.

For this we consider the functional

C(j)[f]=A[and=1andjyou(xandxandRand+1)f]\mathrm{C}^{(j)}[f]=\mathrm{A}\left[\prod_{\begin{subarray}{c}i=1\\ i\neq j\end{subarray}}^{u}\left(x_{i}-x_{i}r_{i}+1\right)f\right] (44)

ThenjAj^{a}equation (14) can be written
(45)

(j)[(xjx)Rjρyoun2]=0\mathbb{C}^{(j)}\left[\left(x_{j}-x\right)^{r}j\rho_{u-n}^{2}\right]=0

We now observe thatC[f]=C(j)[(xjx)Rj+1f]\mathrm{C}[f]=\mathrm{C}^{(j)}\left[\left(x_{j}-x\right)^{r}j^{+1}f\right]soyand,0(xj)=RR(j)(xj)j+and+and,0y_{i,0}\left(x_{j}\right)=r_{r}^{(j)}{}_{j+i+i,0}\left(x_{j}\right),and=0,1i=0,1\ldots, so forξ=xj\xi=x_{j}, the system (43) becomes

0 0 footnotetext: 9. We can replace (xxand)Rand+1with(xandx)Rand+1\left(x-x_{i}\right)^{r}i^{+1}\mathrm{cu}\left(x_{i}-x\right)^{r}i^{+1}because the numbersRand+1r_{i}+1I am a parent.

(40)

V=0nyoudV(xj)RRj+1+and+V,0(j)(xj)=0,and=1,,nyou1\sum_{v=0}^{n-u}d_{v}\left(x_{j}\right)r_{rj+1+i+v,0}^{(j)}\left(x_{j}\right)=0,\quad i=1,\ldots,n-u-1

If we now take this system into account, equation (45) becomes
(47)

V=0nyoudV(xj)RRj+V,c(j)(xj)=0\sum_{v=0}^{n-u}d_{v}\left(x_{j}\right)r_{rj+v,c}^{(j)}\left(x_{j}\right)=0

Eliminating thedV(xj),V=0,1,,nyou1d_{v}\left(x_{j}\right),v=0,1,\ldots,n-u-1from thosenyou+1n-u+1equations (56 (47), we find
(48)

ΓRj,nyou(j)(xj)=0.\Gamma_{r_{j}^{\prime},n-u}^{(j)}\left(x_{j}\right)=0.

If we do thisj=1,2,,youj=1,2,\ldots,uwe find a system that determines our nodesx1,x2,,xyoux_{1},x_{2},\ldots,x_{u}.

Based on what was said in No. 23, it can be written

ΓRj,nyou(j)(xj)=1(nyou+1)!Ct1,t2,,tnyou+1(j)[n=1nyou+1(xjtn)Rj)\Gamma_{r_{j},n-u}^{(j)}\left(x_{j}\right)=\frac{1}{(n-u+1)!}C_{t_{1},t_{2},\ldots,t_{n-u}+1}^{(j)}\left[\prod_{\nu=1}^{n-u+1}\left(x_{j}-t_{\nu}\right)^{r_{j}}\right)
V2(t1,t2,,tnyou+1)]\left.V^{2}\left(t_{1},t_{2},\ldots,t_{n-u+1}\right)\right]

and if we take into account (44),

ΓRj,nyou(j)(xj)==1(nyou+1)!At1,t2,,tnyou+1[(and=1andjyouyou=1V=1nyou+1(xandtV)Rand+1)nyou+1(xV=1nyou+1tV)Rj).V2(t1,t2,,tnyou+1)]\begin{gathered}\Gamma_{r_{j},n-u}^{(j)}\left(x_{j}\right)=\\ =\frac{1}{(n-u+1)!}\mathrm{A}_{t_{1},t_{2},\ldots,t_{n-u+1}}\left[\left(\prod_{\begin{subarray}{c}i=1\\ i\neq j\end{subarray}}^{u-u=1}\prod_{v=1}^{n-u+1}\left(x_{i}-t_{v}\right)^{r_{i}+1}\right)^{n-u+1}\left(x_{v=1}^{n-u+1}t_{v}\right)^{r}j\right)\\ \left..V^{2}\left(t_{1},t_{2},\ldots,t_{n-u+1}\right)\right]\end{gathered}

Once the nodesx1,x2,,xyoux_{1},x_{2},\ldots,x_{u}determined from the indicated system, the polynomial can be foundρnyou\rho_{n-u}calculating from system (47) the coefficientsdV𝒙j,V=0,1,d_{v}\boldsymbol{x}_{j}{}^{\prime},v=0,1,\ldots…,nyou1n-u-1. We can write this polynomial explicitly using the moments of the functionalC[f]\mathrm{C}[f]We have

we (49)
ρnn(x)=Γ1,nyou1(x)Γnyou1\displaystyle\rho_{n-n}(x)=\frac{\Gamma_{1},n-u-1(x)}{\Gamma_{n-u-1}}

IfA[f]\mathrm{A}[f]has a positive orderR1+R2++Ryou+n\geqq r_{1}+r_{2}+\ldots+r_{u}+n, we can obtain the polynomialρnyou\rho_{n-u}and with the help of a well-known formula of Christoffel (see e.g. [11]). For this let us denote withPm(x)\mathrm{P}_{m}(x)the orthogonal polynomial of degreemmattached to the functionA[f]\mathrm{A}[f]This polynomial is well determined formn+R1+R2++Ryoum\leq n+r_{1}+r_{2}+\ldots+r_{u}.

Then the polynomialρnyou(x)and=1you(xxand)RR+1\rho_{n-u}(x)\prod_{i=1}^{u}\left(x-x_{i}\right)^{r}{}^{r+1}differs only by a constant factor of

StCeroStlasi,VI-1

V(Pnyou,Pnyou+1,.,Pn+R1+R2++Ryoux1,x1,.,x1R1+1,x2,x2,,x2R2+1,,xyou,xyou,,xyou,xRyou+1)\mathrm{V}\binom{\mathrm{P}_{n-u},\mathrm{P}_{n-u+1},\ldots.,\mathrm{P}_{n+r_{1}+r_{2}+\ldots+r_{u}}}{\underbrace{x_{1},x_{1},\ldots.,x_{1}}_{r_{1}+1},\underbrace{x_{2},x_{2},\ldots,x_{2}}_{r_{2}+1},\ldots,\underbrace{x_{u},x_{u},\ldots,x_{u},x}_{r_{u}+1}}
  1. 26.

    In particular, ifyou=1u=1, functionalC(1)[f]C^{(1)}[f]it reduces toA[f]A[f]and it is deduced that the nodex1x_{1}is a root of the polynomial

AR1,n1(x)=\displaystyle A_{r_{1},n-1}(x)= |αR1,0(x)αR1,1(x)αR1,n1(x)αR1,1(x)αR1,2(x)αR1,n(x)αR1,n1(x)αR1,n(x)αR1,2n2(x)|=\displaystyle\left|\begin{array}[]{llll}\alpha_{r_{1},0}(x)&\alpha_{r_{1},1}(x)&\ldots&\alpha_{r_{1},n-1}(x)\\ \alpha_{r_{1},1}&(x)&\alpha_{r_{1},2}(x)&\ldots\\ \cdot&\cdot&\alpha_{r_{1},n}(x)\\ \cdot&\cdot&\cdot&\cdot\\ \alpha_{r_{1},n-1}(x)&\alpha_{r_{1},n}(x)&\ldots&\alpha_{r_{1},2n-2}(x)\end{array}\right|=
1n!At1,t2,,tn[(n=1n(xtn))R11V2(t1,t2,,tn)]\displaystyle\frac{1}{n!}A_{t_{1},t_{2},\ldots,t_{n}}\left[\left(\prod_{\nu=1}^{n}\left(x-t_{\nu}\right)\right)^{\left.r_{1}1V^{2}\left(t_{1},t_{2},\ldots,t_{n}\right)\right]}\right. =\displaystyle=

which, ifR1=1r_{1}=1, differs only by a constant factor from the orthogonal polynomial of degreennattached to the functionalityA[f]A[f].

Regarding the calculation of the polynomialρn1\rho_{n-1}, in this case, we can apply formula (49). IfOh[f]\AA [f]has a positivity ordern+R1\geq n+r_{1}HAVE

V(Pn1,Pn,,Pn+R11x1,x1,,x1R1+1)ρn1(x)(xx1)R1+1=\displaystyle\mathrm{V}(\underbrace{\begin{array}[]{l}\mathrm{P}_{n-1},\mathrm{P}_{n},\ldots,\mathrm{P}_{n+r_{1}-1}\\ x_{1},x_{1},\ldots,x_{1}\end{array}}_{r_{1}+1})\rho_{n-1}(x)\left(x-x_{1}\right)_{r_{1}+1}= (50)
=V(Pn1,Pn,,Pn+R1x1,x1,,x1,xR1+1)\displaystyle=\mathrm{V}\left(\begin{array}[]{l}\mathrm{P}_{n-1},\mathrm{P}_{n},\ldots,\mathrm{P}_{n+r_{1}}\\ x_{1},x_{1},\ldots,x_{1},x\\ r_{1}+1\end{array}\right)\cdot

From the determinant in the second member of formula (50) one can easily extract the factor(xx1)R1+1\left(x-x_{1}\right)r_{1}+1because, first of all, the elementPn+and(x)\mathrm{P}_{n+i}(x)can be replaced, based on the othersR1+1r_{1}+1align with (and=1,0,1,2,,R1i=-1,0,1,2,\ldots,r_{1})
V=R1+1n+and(xx1)VV!Pn+andV(x1)=(xx1)R1+1V=R1+1n+and(wx1)VR11V!Pn+and(V)(x1)\sum_{v=r_{1}+1}^{n+i}\frac{\left(x-x_{1}\right)^{v}}{v!}\mathrm{P}_{n+i}^{v}\left(x_{1}\right)=\left(x-x_{1}\right)^{r_{1}+1}\sum_{v=r_{1}+1}^{n+i}\frac{\left(w-x_{1}\right)^{v}-r_{1}-1}{v!}P_{n+i}^{(v)}\left(x_{1}\right)
27. We can calculate onR[xp+n+1]R\left[x^{p+n+1}\right]relative to the rest of the Gauss-type formulas thus obtained. From formulas (36), (44) we deduce

R[xp+n+1]=C(j)[(xjx)Rj+1ρnyou2]\mathrm{R}\left[x^{p+n+1}\right]=\mathrm{C}^{(j)}\left[\left(x_{j}-x\right)^{r}j+1\rho_{n-u}^{2}\right]

and if we take into account (46) we obtain

Rxp+n+1]=n=0nyoudn(xj)γRj+1+nyou+n,0(and)(xj)\left.\mathrm{R}\mid x^{p+n+1}\right]=\sum_{\nu=0}^{n-u}d_{\nu}\left(x_{j}\right)\gamma_{r^{\prime}j+1+n-u+\nu,0}^{(i)}\left(x_{j}\right) (51)

(51), we obtain
(52)

R[xp+n+1]=ΓRj+1,nyou(j)(xj)ΓRj+1,nyou1(j)(xj)\mathrm{R}\left[x^{p+n+1}\right]=\frac{\Gamma_{r_{j}^{\prime}+1,n-u}^{(j)}\left(x_{j}\right)}{\Gamma_{r_{j}+1,n-u-1}^{(j)}\left(x_{j}\right)}

Of course, in this formula we need to replace the nodesx1,x2,x_{1},x_{2},\ldots,xyoux_{u}with their values ​​calculated from the system that is deduced from (48) if we doj=1,2,,youj=1,2,\ldots,u.

In particular ifyou=1u=1, we have

R[xp+n+11]=AR1+1,n1(x1)AR1+1,n2(x1)\mathrm{R}\left[x^{p+n+1_{1}}\right]=\frac{Ar_{1+1,n-1}\left(x_{1}\right)}{Ar_{1+1,n-2}\left(x_{1}\right)}

In general, the calculation of (51) is complicated. Only ifR1=R2=Ryou=1r_{1}=r_{2}=\ldots r_{u}=1(the Gaussian formula is then unique) we have the well-known value fold (in this casep=n1p=n-1)

R[x2n]=AnAn1\mathrm{R}\left[x^{2n}\right]=\frac{A_{n}}{A_{n-1}}
  1. 28.

    Let's consider the particular casen=2,R1=3,R2=1n=2,r_{1}=3,r_{2}=1and let's assume thatA[f]\mathrm{A}[f]has a positivity order3\geqq 3. The nodex1x_{1}is a root of the equation

(α0α2α12)x6+3(α1α2α0α3)x5+3(α1α3+α0α42α22)x4+(8α2α3\left(\alpha_{0}\alpha_{2}-\alpha_{1}^{2}\right)x^{6}+3\left(\alpha_{1}\alpha_{2}-\alpha_{0}\alpha_{3}\right)x^{5}+3\left(\alpha_{1}\alpha_{3}+\alpha_{0}\alpha_{4}-2\alpha_{2}^{2}\right)x^{4}+\left(8\alpha_{2}\alpha_{3}-\right. (53)

7α1α4α0α5)x3+3(α2α4+α1α52α32)x2+3(α3α4α2α5)x+α3α5α+2=0\left.-7\alpha_{1}\alpha_{4}-\alpha_{0}\alpha_{5}\right)x^{3}+3\left(\alpha_{2}\alpha_{4}+\alpha_{1}\alpha_{5}-2\alpha_{3}^{2}\right)x^{2}+3\left(\alpha_{3}\alpha_{4}-\alpha_{2}\alpha_{5}\right)x+\alpha_{3}\alpha_{5}-\alpha_{+}^{2}=0
Nodex2x_{2}is given by the formula

x2=α1x133α2x12+3α3x1α4α0x133α1x12+3α2x1α3x_{2}=\frac{\alpha_{1}x_{1}^{3}-3\alpha_{2}x_{1}^{2}+3\alpha_{3}x_{1}-\alpha_{4}}{\alpha_{0}x_{1}^{3}-3\alpha_{1}x_{1}^{2}+3\alpha_{2}x_{1}-\alpha_{3}}

and the numberR[x6]\mathrm{R}\left[x^{6}\right]of

R[x6]=α2x144α3x13+6α4x124α5x1+α6\displaystyle R\left[x^{6}\right]=\alpha_{2}x_{1}^{4}-4\alpha_{3}x_{1}^{3}+6\alpha_{4}x_{1}^{2}-4\alpha_{5}x_{1}+\alpha_{6}-
(α1x144α2x13+6α3x124α4x1+α5)2α0x144α1x13+6α2x124α3x1+α4\displaystyle\quad-\frac{\left(\alpha_{1}x_{1}^{4}-4\alpha_{2}x_{1}^{3}+6\alpha_{3}x_{1}^{2}-4\alpha_{4}x_{1}+\alpha_{5}\right)^{2}}{\alpha_{0}x_{1}^{4}-4\alpha_{1}x_{1}^{3}+6\alpha_{2}x_{1}^{2}-4\alpha_{3}x_{1}+\alpha_{4}}

Equation (53) has at least two real roots, because, based on the assumptions made, it is of even degree and has at least one real root 10. We therefore have at least two Gauss-type formulas, of which, however, in general, only one corresponds to the minimizing system.

To show this, it is enough to conveniently customize the moments α0=1,α1=0,α2=4,α3=24,α4=200,α5=408\alpha_{0}=1,\ \alpha_{1}=0,\ \alpha_{2}=4,\ \alpha_{3}=24,\ \alpha_{4}=200,\ \alpha_{5}=-408.

Then equation (53) reduces to

(x2+2x8)(x420x3+174x2214x+1556)=0,(x^{2}+2x-8)(x^{4}-20x^{3}+174x^{2}-214x+1556)=0,

which has only two real roots,22and4-4These are the node valuesx1x_{1}in the two corresponding Gauss-type formulas.

Node valuesx2x_{2}are13-13,55, and its valuesR[x6]R[x^{6}]are12920+α6-12920+\alpha_{6}, respectively10760+α6-10760+\alpha_{6}.

The two Gaussian formulas can be written

A[f]13375[(3367f(2)6630f(2)+12600f"(2)+8f(13)]\displaystyle A[f]\approx\frac{1}{3375}\left[\left(3367f(2)-6630f^{\prime}(2)+12600f^{\prime\prime}(2)+8f(-13)\right]\right.
A[f]1729[593f(4)+1696f(4)+1782f"(4)+136f(5)]\displaystyle A[f]\approx\frac{1}{729}\left[593f(-4)+1696f^{\prime}(-4)+1782f^{\prime\prime}(-4)+136f(5)\right]

The first formula alone corresponds to a minimizing system.
29. We will say that the functionalA[f]\mathrm{A}[f]is symmetric of the orderkkif by a linear transformation of the variablexxwe can make moments with odd indicesα2and1,and=12,,k\alpha_{2i-1},i=12,\ldots,kto become null. Such functionals are, for example, those of the formjρ(x)f(x)𝑑x\int^{j}\rho(x)f(x)dx, whereA,ba,bare finite and the functionρ(x)\rho(x)check the propertyρ(x)=ρ(b+Ax)\rho(x)=\rho(b+a-x)Also functions of the formρ(x)f(x)𝑑x\int_{-\infty}^{\infty}\rho(x)f(x)dx, whereρ(x)\rho(x)is a para function.

Returning to the casen=2n=2.R1=3,R2=1r_{1}=3,r_{2}=1studied above, let's assume that the functionalA[f]\mathrm{A}[f]has a positivity order3\geq 3andα5=0\alpha_{5}=0yes an order of symmetry3\geqq 3We can then assumeα1=α3=α5=0\alpha_{1}=\alpha_{3}=\alpha_{5}=0and equation (53) becomes

α0α2x6+3(α0α42α22)x4+3α2α4x2α42=0\alpha_{0}\alpha_{2}x^{6}+3\left(\alpha_{0}\alpha_{4}-2\alpha_{2}^{2}\right)x^{4}+3\alpha_{2}\alpha_{4}x^{2}-\alpha_{4}^{2}=0 (54)

In the conditions we are in (α0>0,α2>0,α0α4α22>0\alpha_{0}>0,\alpha_{2}>0,\alpha_{0}\alpha_{4}-\alpha_{2}^{2}>0) it is immediately seen that this equation has only two real roots unequal and equal in absolute value 11 ). We therefore have two Gauss-type formulas with the sameR[x6]\mathrm{R}\left[x^{6}\right].

In particular, for functionalityJ(0,0)[f]J^{(0,0)}[f], we haveα0=2,α2=23\alpha_{0}=2,\alpha_{2}=\frac{2}{3},α4=25\alpha_{4}=\frac{2}{5}sj the roots of equation (54) are15,15\frac{1}{\sqrt{5}},-\frac{1}{\sqrt{5}}. Doing the calculations, we obtain the Gauss-type formula

f(0,0)[f]=1128[175f(15)405f(15)+323f"(15)+81f(53)]++1281575D6[f]\begin{gathered}f^{(0,0)}[f]=\frac{1}{128}\left[175f\left(\frac{1}{\sqrt{5}}\right)-\frac{40}{\sqrt{5}}f^{\prime}\left(\frac{1}{\sqrt{5}}\right)+\frac{32}{3}f^{\prime\prime}\left(\frac{1}{\sqrt{5}}\right)+81f\left(-\frac{\sqrt{5}}{3}\right)\right]+\\ +\frac{128}{1575}\mathrm{D}_{6}[f]\end{gathered}

and a second Gauss-type formula, with the same remainder, which is deduced from this by replacing5with5\sqrt{5}\mathrm{cu}-\sqrt{5}.
30. We will also consider the case whenn=3,R1=3,R2=R3=1n=3,r_{1}=3,r_{2}=r_{3}=1, assuming thatA[f]A[f]has a positivity order4\geq 4and an order of symmetry4\geq 4Then we can assumeα1=α3=α5=α7=0\alpha_{1}=\alpha_{3}=\alpha_{5}=\alpha_{7}=0and the knotx1x_{1}is 0 root of the equation

0 0 footnotetext: 11. For the discussion it is enough to assume α0=α2=1\alpha_{0}=\alpha_{2}=1. Then ifα42\alpha_{4}\geq 2, the property results from Descartes' rule of signs and if1<α4<21<\alpha_{4}<2from the fact that the derivative of the equation inx2x^{2}it has no real roots.

(55)

α2(α0α4α22)x9+3(α22α42α0α42+α0α2α6)x7+3(3α2α424α22α6+\displaystyle\alpha_{2}\left(\alpha_{0}\alpha_{4}-\alpha_{2}^{2}\right)x^{9}+3\left(\alpha_{2}^{2}\alpha_{4}-2\alpha_{0}\alpha_{4}^{2}+\alpha_{0}\alpha_{2}\alpha_{6}\right)x^{7}+3\left(3\alpha_{2}\alpha_{4}^{2}-4\alpha_{2}^{2}\alpha_{6}+\right.
+α0α4α6)x5+(11α2α4α610α43α0α62)x3+3α6(α42α2α6)x=0.\displaystyle\left.+\alpha_{0}\alpha_{4}\alpha_{6}\right)x^{5}+\left(11\alpha_{2}\alpha_{4}\alpha_{6}-10\alpha_{4}^{3}-\alpha_{0}\alpha_{6}^{2}\right)x^{3}+3\alpha_{6}\left(\alpha_{4}^{2}-\alpha_{2}\alpha_{6}\right)x=0.

Becauseα2(α0α4α22)>0\alpha_{2}(\alpha_{0}\alpha_{4}-\alpha_{2}^{2})>0,α6(α42α2α6)<0\alpha_{6}(\alpha_{4}^{2}-\alpha_{2}\alpha_{6})<0, this equation has at least 3 real roots, namely the root0and two other different and equal in absolute value 12. So we have at least three Gauss-type formulas.

The other two nodes are the roots of the orthogonal polynomial of degree22attached to the functionalC[f]C[f]appropriate.

Apart from a constant factor different from zero, this polynomial, based on formula (49), can be written in the form

[(α0x14+6α2x12+α4)x+4(α2x13+α4x1)][(α2x14+6α4xand3+α6)x+4(α4x19++α6x1)][4(α2x13+α4x1)x+α2x14+6α4x12+α6]2\begin{gathered}{\left[\left(\alpha_{0}x_{1}^{4}+6\alpha_{2}x_{1}^{2}+\alpha_{4}\right)x+4\left(\alpha_{2}x_{1}^{3}+\alpha_{4}x_{1}\right)\right]\left[\left(\alpha_{2}x_{1}^{4}+6\alpha_{4}x_{i}^{3}+\alpha_{6}\right)x+4\left(\alpha_{4}x_{1}^{9}+\right.\right.}\\ \left.\left.+\alpha_{6}x_{1}\right)\right]-\left[4\left(\alpha_{2}x_{1}^{3}+\alpha_{4}x_{1}\right)x+\alpha_{2}x_{1}^{4}+6\alpha_{4}x_{1}^{2}+\alpha_{6}\right]^{2}\end{gathered}

numberR[x8]R\left[x^{8}\right]relative to the rest of the formula is equal to

|α0x14+6α2x12+α44α2x13+4α4x1α2x14+6α4x12+α64α2x13+4α4x1α2x14+6α4x12+α64α4x13+4α6x1α2x14+6α4x12+α64α4x13+4α6x1α4x14+6α6x12+α8||α0x14+6α2x12+α44α2x13+4α4x14α2x13+4α4x1α2x14+6α4x12+α6|\frac{\left|\begin{array}[]{ccr}\alpha_{0}x_{1}^{4}+6\alpha_{2}x_{1}^{2}+\alpha_{4}&4\alpha_{2}x_{1}^{3}+4\alpha_{4}x_{1}&\alpha_{2}x_{1}^{4}+6\alpha_{4}x_{1}^{2}+\alpha_{6}\\ 4\alpha_{2}x_{1}^{3}+4\alpha_{4}x_{1}&\alpha_{2}x_{1}^{4}+6\alpha_{4}x_{1}^{2}+\alpha_{6}&4\alpha_{4}x_{1}^{3}+4\alpha_{6}x_{1}\\ \alpha_{2}x_{1}^{4}+6\alpha_{4}x_{1}^{2}+\alpha_{6}&4\alpha_{4}x_{1}^{3}+4\alpha_{6}x_{1}&\alpha_{4}x_{1}^{4}+6\alpha_{6}x_{1}^{2}+\alpha_{8}\end{array}\right|}{\left|\begin{array}[]{cc}\alpha_{0}x_{1}^{4}+6\alpha_{2}x_{1}^{2}+\alpha_{4}&4\alpha_{2}x_{1}^{3}+4\alpha_{4}x_{1}\\ 4\alpha_{2}x_{1}^{3}+4\alpha_{4}x_{1}&\alpha_{2}x_{1}^{4}+6\alpha_{4}x_{1}^{2}+\alpha_{6}\end{array}\right|}

If the nodex1x_{1}is equal to 0, the other two nodes are
α6α4,α6α4\sqrt{\frac{\alpha_{6}}{\alpha_{4}}},-\sqrt{\frac{\alpha_{6}}{\alpha_{4}}}, andR[x8]=α4α8α62α4\mathrm{R}\left[x^{8}\right]=\frac{\alpha_{4}\alpha_{8}-\alpha_{6}^{2}}{\alpha_{4}}so we find the Gaussian formulaOh[f]12α62[2(α0α62α43)f(0)+(α2α6α42)α0f"(0)+α43[f(α6α4)+f(α6α4)]]\AA [f]\approx\frac{1}{2\alpha_{6}^{2}}\left[2\left(\alpha_{0}\alpha_{6}^{2}-\alpha_{4}^{3}\right)f(0)+\left(\alpha_{2}\alpha_{6}-\alpha_{4}^{2}\right)\alpha_{0}f^{\prime\prime}(0)+\alpha_{4}^{3}\left[f\left(\sqrt{\frac{\alpha_{6}}{\alpha_{4}}}\right)+f\left(-\sqrt{\frac{\alpha_{6}}{\alpha_{4}}}\right)\right]\right]

In particular we have the formulas

J(α,α)[f]=22α+2Γ(α+1)Γ(α+3)25Γ(2α+0)[8(α+1,(2α+57,f(0)+20(α+1)f"(0)++3(2α+7)2(f(52α+7)+f(52α+7))]++1522α+7Γ(α+1)Γ(α+5)(2α+7)Γ(2α+10)D8[f](α>1)13)and(0,0)[f]=1375[456f(0)+20f"(0)+147(f(57)+f(57))]++8441D8[f]\begin{gathered}\begin{aligned} J^{(\alpha,\alpha)}[f]=&\frac{2^{2\alpha+2}\Gamma(\alpha+1)\Gamma(\alpha+3)}{25\Gamma(2\alpha+0)}\left[8\left(\alpha+1,\left(2\alpha+57,f(0)+20(\alpha+1)f^{\prime\prime}(0)+\right.\right.\right.\\ &\left.+3(2\alpha+7)^{2}\left(f\left(\sqrt{\frac{5}{2\alpha+7}}\right)+f\left(-\sqrt{\frac{5}{2\alpha+7}}\right)\right)\right]+\\ &\left.+15\frac{2^{2\alpha+7}\Gamma(\alpha+1)\Gamma(\alpha+5)}{(2\alpha+7)\Gamma(2\alpha+10)}D_{8}[f]\quad(\alpha>-1)^{13}\right)\end{aligned}\\ \begin{aligned} I^{(0,0)}[f]=&\frac{1}{375}\left[456f(0)+20f^{\prime\prime}(0)+147\left(f\left(\sqrt{\frac{5}{7}}\right)+f\left(-\sqrt{\frac{5}{7}}\right)\right)\right]+\\ &+\frac{8}{441}D_{8}[f]\end{aligned}\end{gathered}
  1. 12.

    It would also be interesting here to demonstrate that, under the conditions of the problem, equation (55) has only three real roots.

  2. 13.

    In this formulaΓ(x)\Gamma(x)is the second-order Eulerian function.
    Hf]=π50[44f(0)+5f"(0)+3(f(52)+f(52))]+15π8D8[f]\mathrm{H}\mid f]=\frac{\sqrt{\pi}}{50}\left[44f(0)+5f^{\prime\prime}(0)+3\left(f\left(\sqrt{\frac{5}{2}}\right)+f\left(-\sqrt{\frac{5}{2}}\right)\right)\right]+\frac{15\sqrt{\pi}}{8}\mathrm{D}_{8}[f].

  1. 31.

    To show that not all three formulas correspond to a minimizing system in the sense of$2\mathdollar 2, let us consider the particular case whenα0=1,α2=2,α4=8,α6=2885\alpha_{0}=1,\alpha_{2}=2,\alpha_{4}=8,\alpha_{6}=\frac{288}{5}Equation (55) then becomes

x(x24)(25x6+100x4+400x2+6912)=0x\left(x^{2}-4\right)\left(25x^{6}+100x^{4}+400x^{2}+6912\right)=0

which has real roots 0,2 and - 2. These being the possible values ​​of the nodex1x_{1}, the other two nodesx2,x3x_{2},x_{3}are respectively the roots of the equations5x236=0,5x2+20x+16=05x^{2}-36=0,5x^{2}+20x+16=0, and5x220x+16=05x^{2}-20x+16=0The corresponding values ​​ofR[x8]\mathrm{R}\left[x^{8}\right]SYNTHESISα827.8152\alpha_{8}-\frac{2^{7}.81}{5^{2}}forx=0x=0andα827.8952\alpha_{8}-\frac{2^{7}.89}{5^{2}}for

The Gauss-type formulas thus obtained are
A[f]1324¯[274f(0)+144f"(0)+25(f(65)+f(65))]\displaystyle\qquad\mathrm{A}[f]\approx\frac{1}{3\overline{24}}\left[274f(0)+144f^{\prime\prime}(0)+25\left(f\left(\frac{6}{\sqrt{5}}\right)+f\left(-\frac{6}{\sqrt{5}}\right)\right)\right]
+252(291+1075)f(12.198[6443f(2)3686f(2)+1444f"(2)+\displaystyle+\frac{25}{2}(91+07\sqrt{5})f\left(\frac{1}{2.19^{8}}\left[6443f(2)-3686f^{\prime}(2)+1444f^{\prime\prime}(2)+\right.\right.
A[f]12.193[6443f(2)+3686f(2)+1444f"(2)+\displaystyle\mathrm{A}[f]\approx\frac{1}{2.19^{3}}\left[6443f(-2)+3686f^{\prime}(-2)+1444f^{\prime\prime}(-2)+\right.
+252(2911075f(10+255)+252(291+1075)f(10255)]\displaystyle+\frac{25}{2}\left(291-107\sqrt{5}f\left(\frac{10+2\sqrt{5}}{5}\right)+\frac{25}{2}(291+107\sqrt{5})f\left(\frac{10-2\sqrt{5}}{5}\right)\right]
  1. 32.

    For the case when all multiplicity orders are equal to each other, we have the following property due to P. Turán [12].

Theorem 7. For any linear functionalA[f]\mathrm{A}[f]which has the order of positivitykk, relative to any natural number n and to any system ofnnmultiplicity orders, all equal to the same odd numberRrso thatk12n(R+1)k\geqq\frac{1}{2}n(r+1), there is one and only one Gaussian formula (8).

P. Turá's proof consists in observing that if, in this case, the orthogonality condition (10) is verified, the polynomialπ=(xx1)(xx2)(xxn)\pi=\left(x-x_{1}\right)\left(x-x_{2}\right)\ldots\left(x-x_{n}\right)is minimizing. Indeed, eitherQεPn\mathrm{Q}\epsilon P_{n}a polynomial different fromπ\piWe then have

QR+1πR+1=(R+1)(Qπ)πR+R+12[(Qπ)πR12]2+\mathrm{Q}^{r+1}-\pi^{r+1}=(r+1)(\mathrm{Q}-\pi)\pi^{r}+\frac{r+1}{2}\left[(\mathrm{Q}-\pi)\pi^{\frac{r-1}{2}}\right]^{2}+
+n=1R1n[(Q2π2)QR32n+1πn1]2+\sum_{\nu=1}^{r-1}\nu\left[\left(\mathrm{Q}^{2}-\pi^{2}\right)\mathrm{Q}^{\frac{r-3}{2}-\nu+1}\pi^{\nu-1}\right]^{2}

and finally taking into account orthogonality, the fact thatQπQ-\piis a polynomial of degreen1n-1, and

(Qπ)R12,(Q2π2)QR32n+1πn1,n=1.2,,R12(\mathrm{Q}-\pi)^{\frac{r-1}{2}},\left(\mathrm{Q}^{2}-\pi^{2}\right)\mathrm{Q}^{\frac{r-3}{2}-\nu+1}\pi^{\nu-1},\nu=1.2,\ldots,\frac{r-1}{2}

are polynomials of degree12n(R+1)1k1\frac{1}{2}n(r+1)-1\leqq k-1, we deduce

A[QR+1]A[πR+1]=R+12A[((Qπ)πR12)2]++n=1R12nA[((Q2π2)QR32n+1π1)2]R+12A[((Qπ)πR12)2]>0\begin{gathered}\mathrm{A}\left[\mathrm{Q}^{r+1}\right]-\mathrm{A}\left[\pi^{r+1}\right]=\frac{r+1}{2}\mathrm{~A}\left[\left((\mathrm{Q}-\pi)\pi\frac{r-1}{2}\right)^{2}\right]+\\ +\sum_{\nu=1}^{\frac{r-1}{2}}\nu\mathrm{~A}\left[\left(\left(\mathrm{Q}^{2}-\pi^{2}\right)\mathrm{Q}^{\frac{r-3}{2}-\nu+1}\pi^{-1}\right)^{2}\right]\\ \geq\frac{r+1}{2}\mathrm{~A}\left[\left((\mathrm{Q}-\pi)\pi^{\frac{r-1}{2}}\right)^{2}\right]>0\end{gathered}

soA[QR+1]>A[πR+1]\mathrm{A}\left[Q^{r+1}\right]>\mathrm{A}\left[\pi^{r+1}\right], which proves the theorem. It is seen that the proof remains valid forR=1r=1.

BIBLIOGRAPHY

  1. 1.

    Cauchy A., Sur les fonctions Interpolaires. Comptes rendus Ac. Sci. Parls, 11, pp. 775,789, 1840.

  2. 2.

    Jackson D., The theory of approximation, 1930.

  3. 3.

    Ditto thons. Year hons. A. Diff Ma1cm., (2), 25, pp. 184,152, 1924.

  4. 4.

    Mar a G. Sur, Dalgorithme toulour,

  5. 5.

    Po1 ya G., Jur un algorithme toujours convergent pourobtenir les polynomes de mention continuous, quel coque. Fully rendered Ac. Sci. Parls, 157, pp, 840-843, 1913.

  6. 6.

    Popové timatic, Burresti, 1926, pr 35.38, algebraic, Bul. Soc. Studentilor in ma*

  7. 7.

    Ditto Sur a is bem Gis.

  8. 8.

    Idem Acad. Roumaine, XVI, pp. 214=217, 1934.

  9. 9.

    Id Work on the formal remainder in some approximation formulas of analysis, sess. $t., Acad RPR, 1950, pp. 183:186.

  10. 10.

    Ditto Circle. Math., III, pp. 53:122, 1952.

1955

Related Posts