ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE OF ANALYTIC FUNCTIONS

MATHEMATICS OF COMPUTATION Volume, Number, Pages S 5-57(XX)- ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE OF ANALYTIC FUNCTIONS GRADIMIR V. MILOVANOVIĆ AND MIODRAG M. SPALEVIĆ This paper is dedicated to Professor Walter Gautschi on the occasion of his 75th birthday Abstract. We study the kernels of the remainder term R n,s(f) of Gauss- Turán quadrature formulas n s f(t)w(t) dt = A i,ν f (i) (τ ν) + R n,s(f) (n N; s N ) ν= i= for classes of analytic functions on elliptical contours with foci at ±, when the weight w is one of the special Jacobi weights w (α,β) (t) = ( t) α ( + t) β (α = β = /; α = β = / + s; α = /, β = / + s; α = / + s, β = /). We investigate the location on the contour where the modulus of the kernel attains its maximum value. Some numerical examples are included.. Introduction Quadrature formulae with multiple nodes appeared more than years after the famous Gaussian quadratures. Starting from the Hermite interpolation formula and taking any system of n distinct nodes {τ,...,τ n } with arbitrary multiplicities m ν (ν =,...,n), Chakalov [] in 9 obtained such a general quadrature, which is exact for all algebraic polynomials of degree at most m + +m n. Taking all the multiplicities to be equal, Turán [] was the first who introduced the corresponding quadrature formula of Gaussian type. Let w be an integrable weight function on the interval (,). In this paper we consider the Gauss-Turán quadrature formula with multiple nodes, (.) f(t)w(t)dt = n s A i,ν f (i) (τ ν ) + R n,s (f), ν= i= where A i,ν = A (n,s) i,ν, τ ν = τ ν (n,s) (i =,,...,s; ν =,...,n), which is exact for all algebraic polynomials of degree at most (s + )n. The nodes τ ν in (.) Received by the editor February 7,. 99 Mathematics Subject Classification. Primary A55; Secondary 65D3, 65D3. Key words and phrases. Gauss-Turán quadrature, s-orthogonality, zeros, multiple nodes, weight, measure, degree of exactness, remainder term for analytic functions, error estimate, contour integral representation, kernel function. The authors were supported in part by the Serbian Ministry of Science, Technology and Development (Project: Applied Orthogonal Systems, Constructive Approximation and Numerical Methods). c 997 American Mathematical Society

G. V. MILOVANOVIĆ AND M. M. SPALEVIĆ must be the zeros of the (monic) polynomial π n,s (t) which minimizes the integral Φ(a,a,...,a n ) = π n (t) s+ w(t)dt, where π n (t) = t n + a n t n + + a t + a. In order to minimize Φ we must have the following orthogonality conditions (.) π n (t) s+ t k w(t)dt =, k =,,...,n. R Polynomials π n = π n,s which satisfy this type of orthogonality (.) (so called power orthogonality ) are known as s-orthogonal (or s-self associated) polynomials with respect to the measure dλ(t) = w(t)dt. For s = this reduces to the standard case of orthogonal polynomials. Several classes of s-orthogonal polynomials, as well as their generalizations known as σ-orthogonal polynomials, were investigated mainly by Italian mathematicians, e.g., Ossicini [6, 7], Ghizzetti and Ossicini [, 3], Ossicini and Rosati [9, 3, 3], Gori Nicolò-Amati [5] (see the survey paper [] for details and references). A generalization of the Gauss-Turán quadrature formula (.) to rules having nodes with arbitrary multiplicities was derived independently by Chakalov [3, ] and Popoviciu [36]. Important theoretical progress on this subject was made by Stancu [, ] (see also [3]). Methods for constructing the nodes τ ν and/or coefficients A i,ν in the Gauss- Turán quadratures, as well as in the generalized Chakalov-Popoviciu-Stancu formulas, can be found in [,,, 3,, 5, 39,, 3, 5]. The remainder term in formulas with multiple nodes was studied by Chakalov [3], Ionescu [9], Ossicini [7], Pavel [3, 33, 3], and Milovanović and Spalević [5]. The case of holomorphic functions f in the Gauss-Turán quadrature (.) was considered by Ossicini and Rosati [9]. In this paper we consider the remainder term R n,s (f) of Gauss-Turán quadrature formulas for classes of analytic functions on elliptical contours, when the weight function w in (.) is one of the special Jacobi weights w (α,β) (t) = ( t) α ( + t) β, with parameters α = β = ; α = β = + s; α =, β = + s; α = + s, β =, where s N. The reason for these choices is explained near the end of Section. The paper is organized as follows. The remainder term of Gauss-Turán formulas for analytic functions and some properties of the kernels in the contour representations of the remainder terms are given in Section. The cases of elliptic contours with foci at the points ±, when w is any one of the four Jacobi weight functions, are studied in Section 3. More precisely, the location on the contour where the modulus of the kernel attains its maximum value is investigated. Some numerical examples are included.. The remainder term in Gauss-Turán quadrature formulae Let Γ be a simple closed curve in the complex plane surrounding the interval [,] and D be its interior. If the integrand f is analytic in D and continuous on D, R

ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE 3 then the remainder term R n,s (f) in (.) admits the contour integral representation (cf. [9]) (.) R n,s (f) = K n,s (z)f(z)dz. πi The kernel is given by (.) K n,s (z) = where (.3) n,s (z) = Γ n,s (z), z / [,], [π n,s (z)] s+ w(t) [π n,s(t)] s+ z t dt, n N, and π n,s (t) is the monic s-orthogonal polynomial with respect to the measure dλ(t) = w(t)dt on (,). For s =, the formulas (.) and (.) reduce to the corresponding formulas for Gaussian quadratures. An alternative representation for K n,s (z) is ( ) K n,s (z) = R n,s = z w(t) n z t dt s ν= i= i!a i,ν (z τ ν ) i+. Put N = (s + )n and n t N = (t τ ν ) s+ + q(t) = π n,s (t) s+ + q(t) (q P N ), ν= where P m denotes the set of all algebraic polynomials of degree at most m. Expanding the integrand of (.3) in descending powers of z, and using (.) and R n,s (t N ) = R n,s (π s+ n,s ) = π n,s (t) dµ(t) = π n,s dµ, where dµ(t) = π n,s (t) s w(t)dt, we conclude that K n,s (z) = R n,s(t N ) z N+ + = π n,s ( dµ z N+ + C z + C ) z +, where C,C,... are constants. The integral representation (.) leads to the error estimate (.) R n,s (f) l(γ) ( ) ( ) max π K n,s(z) max f(z), z Γ z Γ where l(γ) is the length of the contour Γ. We thus have to study the magnitude of K n,s (z) on Γ. It seems that the first unified approach described above was taken by Donaldson and Elliott [5]. They applied it to several kinds of interpolatory and noninterpolatory quadrature rules. Error bounds for Gaussian quadratures of analytic functions were studied by Gautschi and Varga [] (see also []). In particular, they investigated some cases with special Jacobi weights with parameters ±/ (Chebyshev weights). The cases of Gaussian rules with Bernstein-Szegő weight functions and with some symmetric weights including especially the Gegenbauer weight were studied by Peherstorfer [35] and Schira [3], respectively. Some of the results have been extended to Gauss-Radau and Gauss-Lobatto formulas (cf. Gautschi [6], Gautschi and Li [7], Schira [37], Hunter and Nikolov []).

G. V. MILOVANOVIĆ AND M. M. SPALEVIĆ In the sequel we give some properties of the kernel (.). Lemma.. Let the kernel K n,s (z) be given by (.) and (.3). Then, for each z C \ [,], (.5) K n,s (z) = K n,s (z) Moreover, if the weight function in (.) is even, i.e., w( t) = w(t), then (.6) K n,s ( z) = K n,s (z). Proof. According to (.) it is clear that K n,s (z) = n,s (z) [π n,s (z)] s+ = K n,s(z), implying (.5). If w is an even function, i.e., w( t) = w(t), we have π n,s ( z) = ( ) n π n,s (z) and i.e., so that n,s (z) = n,s (z) = ( ) n(s+) K n,s ( z) = w(t) [π n,s(t)] s+ z t dt = w(t) [π n,s(t)] s+ z + t w( t) [π n,s( t)] s+ z + t dt, dt = ( ) n(s+) n,s ( z), n,s ( z) [π n,s ( z)] s+ = ( )n(s+) n,s (z) ( ) n(s+) [π n,s (z)] = K n,s(z) = K s+ n,s (z). Thus, in this case we get (.6). A particularly interesting case is the Chebyshev measure dλ (t) = ( t ) / dt. In 93, S. Bernstein [] showed that the monic Chebyshev polynomial ˆTn (t) = T n (t)/ n minimizes all integrals of the form π n (t) k+ t dt (k ). This means that the Chebyshev polynomials T n are s-orthogonal on (,) for each s. Ossicini and Rosati [9] found three other measures dλ k (t) (k =,3,) for which the s-orthogonal polynomials can be identified as Chebyshev polynomials of the second, third, and fourth kind: U n, V n, and W n, which are defined by U n (cos θ) = sin(n + )θ, V n (cos θ) = cos(n + )θ sin θ cos θ, W n (cos θ) = sin(n + )θ sin θ, respectively (cf. Gautschi and Notaris [9]). However, these measures depend on s, dλ (t) = ( t ) /+s dt, dλ 3 (t) = ( + t)/+s ( t) / dt, dλ (t) = ( t)/+s dt. ( + t) / It is easy to see that W n ( t) = ( ) n V n (t), so that in the investigation it is sufficient to study only first three Jacobi measures dλ k (t), k =,,3. Recently, Ossicini, Martinelli, and Rosati [] have proved the convergence as n + (alternatively, as s + ), of the Gauss-Turán quadrature formula

ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE 5 (.) for the cases dλ (t) and dλ (t), on the basis of results from [9], by taking f to be a holomorphic function on int Γ, where the contour Γ is an ellipse with foci at ± and sum of semiaxes >. Using estimates obtained for R n,s (f), they proved the convergence and rate of convergence of the quadrature formulae, R n,s (f) = O ( ρ n(s+)), n +. 3. The maximum modulus of the kernel on confocal ellipses In this section we take as the contour Γ an ellipse with foci at the points ± and sum of semiaxes >, { (3.) E = z C: z = ( e iθ + e iθ) }, θ < π. When, then the ellipse shrinks to the interval [,], while with increasing it becomes more and more circle-like. Since the ellipse E has length l(e ) = ε E(ε), where ε is the eccentricity of E, i.e., ε = /( + ), and E(ε) = π/ ε sin θ dθ is the complete elliptic integral of the second kind, the estimate (.) reduces to (3.) R n,s (f) E(ε) ( ) max K n,s (z) f, ε = πε z E +, where f = max f(z). As we can see, the bound on the right in (3.) is a z E function of, so that it can be optimized with respect to >. In this section we study the magnitude of K n,s (z) on the contour E. More precisely, for the measures dλ k (k =,,3) defined at the end of the previous section, we investigate the locations on the confocal ellipses where the modulus of the corresponding kernels attain their maximum values. Because of (.5), i.e., symmetry with respect to the real axis, the consideration of K n,s (z), when z = ( e iθ + e iθ) E, may be restricted to the interval θ π. Moreover, if the weight function is even, as in the cases of dλ (t) and dλ (t) (symmetry with respect to both coordinate axes), the consideration may be restricted to the first quarter of E, i.e., to the interval θ π/ (see (.6)). In the sequel we give explicit representations of the kernels K n,s (ν) on the ellipse E for the measures dλ ν (t), ν =,,3, and discuss the maximum points on this ellipse in order to get the exact value of max K n,s(z) (ν) or some estimate. z E 3.. The measure dλ (t) = ( t ) / dt. According to (.3), in this case we have n,s (z) = ( t ) / [ n T n (t)] s+ z t dt, n N, z / [,],

6 G. V. MILOVANOVIĆ AND M. M. SPALEVIĆ where T n (t) is the Chebyshev polynomial of the first kind of degree n. By substituting t = cos θ we obtain n,s (z) = ( n)(s+) π = ( n)(s+) s π [cos nθ] s+ z cos θ z cos θ dθ ( ) s + cos(s + k)nθ dθ, k where for the transformation of [cosnθ] s+ we used a formula from [9, p. 3]. Now, the kernel (.) has the form ( ) s + π s cos(s + k)nθ dθ k K n,s(z) () z cos θ = [T n (z)] s+. Furthermore, using [7, Eq. 3.63.], one finds (see also [, p. 76]) π ( (3.3) z m z ), m N, and we obtain K () n,s(z) = cos mθ z cos θ dθ = s π z π z ( ) s + (z z ) (s+ k)n k [T n (z)] s+. It is well-known that (3.) T n (z) = [( z + ) n ( z + z ) n ] z, z C. i.e., Putting z = (u + u ), we get where K n,s(z) () π = (u u )u n [u n + u n ] s+ (3.5) Z () n,s(u) = Introducing ( ) s + k u, (s k)n K n,s(z) () πz n,s(u) () = (u u )u n [u n + u n ] s+, ( ) s + u nk. s + k + (3.6) a j = a j ( ) = ( j + j ), j N, >, we have u u = (a cos θ) and u n + u n = (a n + cos nθ), when u = e iθ, so that (3.7) K n,s(z) () = s π Z n,s( e () iθ ), z E, n (a cos θ) / (a n + cos nθ) s+/

ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE 7 6 5 5 5 6 3 3 5. 3.5 3.5 3 Figure. The functions θ K (), (z) and θ K() 5, (z) (z = (u + u ), u = e iθ ) for Chebyshev weight of the first kind and =. (top) and =.5 (bottom) where Z n,s(u) () is given by (3.5) and the ellipse E by (3.). Note that the case s =, for which Z () n, (u) =, was analyzed in [, Eq. (5.)]. 6 5 3 6 3.5 3.5.5.5 3....6.. 3 3 3 Figure. The function θ K (),s (z), z = (u+u ), u = e iθ, for s = (dashed line), s = (dot-dashed line), s = 3 (solid line), when =.5, =. (top), and =., =. (bottom)

G. V. MILOVANOVIĆ AND M. M. SPALEVIĆ Table. Maximum value of K () n,s(z), z E, and the bound (3.) for n =,5,, s =,,3, and =.,.5, and. =. =.5 =. s n max bound max bound max bound 9.996(3).73() 5.996().55() 3.69() 5.5() 5.779() 5.35().9( ).56( ).( 6).( 6).97().36().9( 6).9( 6) 5.76( 5) 5.76( 5) 7.53(5).66(6).() 3.993().6() 3.57() 5.6().769(3) 5.7( ) 5.( ).5( ).56( ).73().733().9( ).9( ) 9.6( 3) 9.6( 3) 3 6.369(7).65() 5.9().7(3).6().() 5.55(3) 6.3(3).3( 5).573( 5) 6.35( ) 6.39( ).7().7() 5.37( ) 5.( ).77( 3).77( 3) An analysis of (3.7) shows that the point of the maximum of K n,s(z) () for a given depends on n. The graphics θ K () n, (z) (z = (u+u )/, u = e iθ ) for n = and n = 5 are displayed in Fig., when =. and =.5. The cases for s =,,3, when n = and =.5,.,., and., are presented in Fig.. Using the inequality (see [, Proof of Thm. 5.]) (a cos θ)(a n + cos nθ) (a )(a n + ), θ π/, a simple estimate of (3.7) can be given in the form (3.) K () n,s(z) πz () n,s( ) n( n n ) s ( )( n + n ) ) = K n,s( ( n () ( + + n ) s ), n n where Z n,s(u) () is defined by (3.5). By the crude inequality Z n,s( ) () < Z n,s() () = s ( > ), the inequality (3.) can be simplified to ( ) s (3.9) K n,s(z) () π n( n + n )( ) n n. Numerical values of the actual maximum of K () n,s(z), when z E, and the corresponding bounds (3.) for some selected values of n, s, and are presented in Table. (Numbers in parenthesis indicate decimal exponents.) Based on the previous calculation we can state the following conjecture: Conjecture 3.. For each fixed > and s N there exists n = n (,s) N such that ) max K n,s(z) () = K n,s( () z E ( + ) for each n n. Using (3.9), the estimate (3.) becomes ( M (3.) R n,s (f) n( n + n ) n n ) s,

ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE 9 where M = E(ε) f ( + )/( ). On the basis of (3.) we conclude that the corresponding Gauss-Turán quadrature formulae converge if s is a fixed integer and n +, since lim R n,s(f) =. n + Moreover, we conclude that R n,s (f) = O( n(s+) ), when n +. Assuming that /( n n ) <, we see that also R n,s (f), when s + and n is fixed. This condition is satisfied if n n >, i.e., if n > +. The same conclusion was obtained in [, Eq. (5.5)]. Remark 3.. Recently, Gori and Micchelli [6] have introduced for each n a class of weight functions defined on [, ] for which explicit Gauss-Turán quadrature formulas can be found for all s. Indeed, these classes of weight functions have the peculiarity that the corresponding s-orthogonal polynomials, of the same degree, are independent of s. This class includes certain generalized Jacobi weight functions w n,µ (t) = U n (t)/n µ+ ( t ) µ, where U n (cos θ) = sinnθ/sin θ (Chebyshev polynomial of the second kind) and µ >. In this case, the Chebyshev polynomials T n appear as s-orthogonal polynomials. Since U n (cos θ) = sin nθ sin θ n, i.e., U n (t)/n µ+ ( t ) µ+/, by arguing, for example, in an analogous way as in [] we can obtain in this case that lim R n,s(f) =, under the previous s + condition n > +, where n is a fixed positive integer. 3.. The measure dλ (t) = ( t ) s+/ dt, s N. In this case we have n,s (z) = ( t ) /+s [ n U n (t)] s+ dt, z t n N, z / [,], where U n (t) is the Chebyshev polynomial of the second kind for the weight function w(x) = t, for which U n (cos θ) = By substituting t = cos θ we obtain n,s (z) = n(s+) π = n(s+) s π sin(n + )θ. sin θ sin θ z cos θ [sin(n + )θ]s+ dθ sinθ z cos θ ( ) s + ( ) s+k sin(s + k)(n + )θ dθ, k where for the transformation of [sin(n+)θ] s+ we used a formula from [9, p. 3]. By using the well-known representation [ U n (z) = (z + z ) n+ (z z ) n+], z the substitution z = (u + u ) and the formula (see [7, Eq. 9.63.3]) π sin(m + )θ sin θ z cos θ dθ = π u m+, m N,

G. V. MILOVANOVIĆ AND M. M. SPALEVIĆ yield K () n,s(z) = [ π s u n+ u u u n+ u (n+) ] s+ ( ) s + ( ) s+k k u (s k)(n+). The sum on the right-hand side we denote by Z n,s(u) () and rewrite it in the following form ( ) s + (3.) Z n,s(u) () = ( ) k u (n+)k, s + k + so that i.e., K () n,s(z) = (3.) K () n,s(z) = where ( π s n+ π s n+ a cos θ a n+ cos(n + )θ (a cos θ) s+/ (a n+ cos(n + )θ) s z = ( eiθ + e iθ ) E and a j is defined by (3.6). Now, we consider the last factor in (3.) when n is odd. ) s+/ Z () n,s( e iθ ), Z () n,s( e iθ ) (a n+ cos(n + )θ) /, Lemma 3.3. Let a j and Z () n,s(u) be defined by (3.6) and (3.), respectively. If n is odd, then Z n,s( e () iθ ) (a n+ cos(n + )θ) Z n,s(i ) (), θ π/, / (a n+ ) / with equality for θ = π/. Proof. First we note that (3.3) Z () n,s(u) = where ζ n,s (u): = = [(s )/] ν= ( ν+ k=ν ) { if s is odd, u (n+)s as well as ζ n,s ( e iθ ) = ζ n,s (i ). Putting ν+ ( ) s + S ν (u): = = u ν(n+) s + ν + where k=ν (3.) α = s ν s + ν + = if s is even, + ζ n,s (u), ( ) s + u (ν+)(n+) s + ν + ( ) s + u ν(n+)( αu (n+)), s + ν + and α <,

ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE we see that S ν ( e iθ ) = ( s + ) ν(n+) s + ν + α (n+) cos(n + )θ + α (n+). Now, we consider the quotient F ν (,θ,n): = S ν ( e iθ ) (a n+ cos(n + )θ) /, when n is odd, and for θ π/ we wish to prove the inequality F ν (,θ,n) F ν (,π/,n), i.e., (3.5) i.e., S ν ( e iθ ) (a n+ cos(n + )θ) S ν (i ), θ π/. / (a n+ ) / Using the previous facts, inequality (3.5) reduces to (a n+ )( q cos(n + )θ + q ) (a n+ cos(n + )θ)( q), ( cos(n + )θ)( qa n+ + q ), where q = α (n+). Using (3.6), we find qa n+ + q = ( α)( α (n+) ), so that the previous inequality becomes (3.6) ( cos(n + )θ)( α)( α (n+) ). Since cos(n + )θ, >, and α < (see (3.)) we conclude that inequality (3.6) is true. This also proves inequality (3.5). According to (3.3) and (3.5) we have Z () n,s( e iθ ) (a n+ cos(n + )θ) / [(s )/] ν= F ν (,θ,n) + ζ n,s ( e iθ ) (a n+ cos(n + )θ) / [(s )/] ν= S ν (i ) (a n+ ) + ζ n,s (i ) / (a n+ ) / = with equality holding for θ = π/. Z () n,s(i ) (a n+ ) /, Theorem 3.. If dλ(t) = ( t ) s+/ dt on (,), s N, and n is odd, then ( (3.7) max K n,s(z) () i = K n,s () ( )), z E i.e., the maximum of K () n,s(z) (n odd) on E is attained on the imaginary axis. Proof. For the second factor in (3.), it is obvious that (a cos θ) s+/ (a n+ cos(n + )θ) s (a + ) s+/, for all θ, all n, (a n+ ) s with equality holding when θ = π/ and n is odd. Now, this inequality and Lemma 3.3 give the desired result.

G. V. MILOVANOVIĆ AND M. M. SPALEVIĆ When n is even in Theorem 3., computation shows that the maximum of K n,s(z) () on the ellipse E is attained slightly off the imaginary axis. The graphics θ K n,s(z) () (z = (u + u )/, u = e iθ ) for n = and s =,,3 are displayed in Fig. 3, when =.5,.,., and.. 3..6... 3.5 3.3.5...75..5.5 3 3 Figure 3. The function θ K (),s (z), z = (u+u ), u = e iθ, for s = (dashed line), s = (dot-dashed line), s = 3 (solid line), when =.5, =. (top), and =., =. (bottom) As in the case of the measure dλ (t) we can get here a simple crude bound of the remainder, (3.) R n,s (f) M ( + ) s+, M = E(ε) f ( + ), n+ n+ (n+) which holds for each n N. According to (3.) we may conclude that the corresponding Gauss-Turán quadrature formulae converge if s is a fixed integer and n +. Moreover, R n,s (f) = O( n(s+) ), when n +. Assuming that ( + )/( n+ n ) <, we see that also R n,s (f), when s + and n is fixed. This condition is satisfied if n+ ( + ) n >, i.e., n > + + + 6 +. 3.3. The measure dλ 3 (t) = ( t) / ( + t) /+s, s N. In this case it was shown (see [9, ]) that the monic s-orthogonal polynomials are the monic Jacobi orthogonal polynomials with parameters α = /, β = /, i.e., (n + )θ cos π n,s (t) = n n V n (t) = cos θ, t = cos θ.

ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE 3 Therefore, by (.3), where w(t) = ( t) / (+t) /+s, substituting t = cos θ, we have π cos θ [ ] s+ (n + )θ cos n,s (z) = s+ n(s+) dθ, n N, z / [,]. z cos θ Using the representation (see [9]) we obtain V n (z) = π cos θ s+ K n,s(z) (3) z cos θ = ) ( T n+ ( + z), ( + z) [ ] s+ (n + )θ cos dθ ) ( T n+ ( + z) ( + z) The numerator of the last fraction has the form s+ π cos θ [ cos z cos θ = s+ π ] s+ (n + )θ dθ cos θ z cos θ s ( s + = ( ) s + π cos θ s k z cos θ = ( ) s + π s k = = π s z π s (u u ) k [ cos(m + )θ z cos θ ) cos cos (m + )θ + s+. (s + k)(n + )θ dθ ] cos mθ dθ z cos θ ( ) s + [(z z ) m+ ( + z z ) ] m k ( ) s + u + k u m+, where (3.3) has been used and we have put m = (n + )(s k) + n and z = (u + u ). On the other hand, according to (3.) and the fact that u + (z + ) = u, z = (u + u ), dθ

G. V. MILOVANOVIĆ AND M. M. SPALEVIĆ we get Therefore, K (3) n,s(z) = s π(u + ) u n+ (u u ) Using (3.6) and putting ) ( T n+ ( + z) u n+ + u n =. (z + ) u + ( (3.9) Z (3) n,s(u) = we obtain u + u n+ + u n ) s+ ( ) s + u (n+)k, s + k + ( ) s + u (n+)k. s + k + (3.) K n,s(z) (3) = s π (a + cos θ) s+ Z n,s( e (3) iθ ) n+/ (a cos θ) / (a n+ + cos(n + )θ), s+/ when z = ( eiθ + e iθ ) E. An analysis of (3.) shows that the point of the maximum of K (3) n,s(z) for a given depends on n as in the case of the measure dλ (t). If there exists a sequence of the local maxima, numerical experiments show that it decreases when θ runs over [,π]. Because of that and better clarity in the following figures, the graphics of the function θ K (3) n,s(z) (z = (u + u )/, u = e iθ ) for some selected n, s, and, are presented only for θ [,π/]. The case n =, s =, is given in Fig. for =.,.5,., and.5. The graphics for s =,,3, when n = and =. and.5, are presented in Fig. 5. On the basis of numerical experiments a similar conjecture for K n,s(z) (3) on the ellipse E as Conjecture 3. can be stated. A useful estimate of (3.) can be given by using the fact that Z (3) Z (3) n,s( ) < Z (3) n,s() = s and the inequality (see [, p. 79]) n,s( e iθ ) a + cos θ (a cos θ)(a n+ + cos(n + )θ) a + (a )(a n+ + ), θ π, which is equivalent to a + cos θ (a cos θ)(a n+ + cos(n + )θ) a + (a )(a n+ + ), θ π, because of a = a and a cos θ = (a cos θ)/(a + cos θ). In this way, we get K n,s(z) (3) π ( ) s a + (a + cos θ), n+/ (a )(a n+ + ) a n+ + cos(n + )θ i.e., ( ) s (3.) K n,s(z) (3) π( + ) ( / + / ) n( )( n+ + n, ) n+/ n /

ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE 5 35 3 5 5 5 5 5 3 3 6...6. 3 3. Figure. The function θ K (3), (z) (z = (u+u ), u = e iθ ) for =., =.5 (top) and =., =.5 (bottom) 7.5 5.5 7.5 5.5...6.. 3 3 Figure 5. The function θ K (3),s (z), z = (u+u ), u = e iθ, for s = (dashed line), s = (dot-dashed line), s = 3 (solid line), when =. (left) and =.5 (right) and then, the estimate of the corresponding remainder (3.) becomes ( ) s M ( + ) R n,s (f) n( n+ + n ) n+ n, M = E(ε) f ( + ) +. It is clear that R n,s (f) = O( n(s+) ) when n +. Assuming that + n+ <, n which is equivalent to n > + + + +, we conclude that R n,s (f) when s + and n is fixed.

6 G. V. MILOVANOVIĆ AND M. M. SPALEVIĆ Acknowledgments The authors would like to thank the referee for a careful reading of the manuscript and for his valuable comments. References. S. Bernstein, Sur les polynomes orthogonaux relatifs à un segment fini, J. Math. Pures Appl. 9 (93), 7 77.. L. Chakalov, Über eine allgemeine Quadraturformel, C.R. Acad. Bulgar. Sci. (9), 9. 3. L. Chakalov, General quadrature formulae of Gaussian type, Bulgar. Akad. Nauk Izv. Mat. Inst. (95), 67 (Bulgarian) [English transl. East J. Approx. (995), 6 76].. L. Chakalov, Formules générales de quadrature mécanique du type de Gauss, Colloq. Math. 5 (957), 69 73. 5. J. D. Donaldson and D. Elliott, A unified approach to quadrature rules with asymptotic estimates of their remainders, SIAM J. Numer. Anal. 9 (97), 573 6. 6. W. Gautschi, On the remainder term for analytic functions of Gauss-Lobatto and Gauss- Radau quadratures, Rocky Mountain J. Math. (99) 9 6. 7. W. Gautschi and S. Li, The remainder term for analytic functions of Gauss-Radau and Gauss- Lobatto quadrature rules with multiple points, J. Comput. Appl. Math. 33 (99), 35 39.. W. Gautschi and G. V. Milovanović, S orthogonality and construction of Gauss-Turán type quadrature formulae, J. Comput. Appl. Math. 6 (997), 5. 9. W. Gautschi and S. E. Notaris, Gauss-Kronrod quadrature formulae for weight function of Bernstein-Szegő type, J. Comput. Appl. Math. 5 (99), 99.. W. Gautschi and R. S. Varga, Error bounds for Gaussian quadrature of analytic functions, SIAM J. Numer. Anal. (93), 7 6.. W. Gautschi, E. Tychopoulos and R. S. Varga, A note on the contour integral representation of the remainder term for a Gauss-Chebyshev quadrature rule, SIAM J. Numer. Anal. 7 (99), 9.. A. Ghizzetti and A. Ossicini, Quadrature formulae, Akademie Verlag, Berlin, 97. 3. A. Ghizzetti and A. Ossicini, Sull esistenza e unicità delle formule di quadratura gaussiane, Rend. Mat. (6) (975), 5.. G. H. Golub and J. Kautsky, Calculation of Gauss quadratures with multiple free and fixed knots, Numer. Math. (93), 7 63. 5. L. Gori Nicolò-Amati, On the behaviour of the zeros of some s-orthogonal polynomials, in Orthogonal Polynomials and Their Applications, nd Int. Symp. (Segovia, 96), Monogr. Acad. Cienc. Exactas, Fis., Quim., Nat., Zaragoza, 9, pp. 7 5. 6. L. Gori and C. A. Micchelli, On weight functions which admit explicit Gauss-Turán quadrature formulas, Math. Comp. 65 (996), 567 5. 7. I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products, Sixth Edition (A. Jeffrey and D. Zwillinger, eds.), Academic Press, San Diego,.. D. Hunter and G. Nikolov, On the error term of symmetric Gauss-Lobatto quadrature formulae for analytic functions, Math. Comp. 69 (999), 69. 9. D. V. Ionescu, Restes des formules de quadrature de Gauss et de Turán, Acta Math. Acad. Sci. Hungar. (967), 3 95.. S. Karlin and A. Pinkus, Gaussian quadrature with multiple nodes, in Studies in Spline Functions and Approximation Theory (S. Karlin, C. A. Micchelli, A. Pinkus, and I. J. Schoenberg, eds.), Academic Press, New York, 976, pp. 3.. G. V. Milovanović, Construction of s-orthogonal polynomials and Turán quadrature formulae, in Numerical Methods and Approximation Theory. III (Niš, 97), (G. V. Milovanović, ed.), Univ. Niš, Niš, 9, pp. 3 3.. G. V. Milovanović, Quadratures with multiple nodes, power orthogonality, and momentpreserving spline approximation, in Numerical Analysis in th Century, Vol. 5 (W. Gautschi, F. Marcellán, and L. Reichel, eds.), J. Comput. Appl. Math. 7 (), 67 6. 3. G. V. Milovanović and M. M. Spalević, A numerical procedure for coefficients in generalized Gauss Turán quadratures, Filomat (Niš) 9() (995),.. G. V. Milovanović and M. M. Spalević, Construction of Chakalov Popoviciu s type quadrature formulae, Rend. Circ. Mat. Palermo 5 (99), 65 636.

ERROR BOUNDS FOR GAUSS-TURÁN QUADRATURE FORMULAE 7 5. G. V. Milovanović and M. M. Spalević, Quadrature formulae connected to σ orthogonal polynomials, J. Comput. Appl. Math. (), 69 637. 6. A. Ossicini, Costruzione di formule di quadratura di tipo Gaussiano, Ann. Mat. Pura Appl. () 7 (966), 3 37. 7. A. Ossicini, Le funzioni di influenza nel problema di Gauss sulle formule di quadratura, Matematiche (Catania) 3 (96), 7 3.. A. Ossicini, M. R. Martinelli, and F. Rosati, Funzioni caratteristiche e polinomi s-ortogonali, Rend. Mat. (99), 355 366. 9. A. Ossicini and F. Rosati, Funzioni caratteristiche nelle formule di quadratura gaussiane con nodi multipli, Boll. Un. Mat. Ital. () (975), 37. 3. A. Ossicini and F. Rosati, Comparison theorems for the zeros of s-orthogonal polynomials, Calcolo 6 (979), 37 3 (Italian). 3. A. Ossicini and F. Rosati, s-orthogonal Jacobi polynomials, Rend. Mat. Appl. (7) (99), 399 3 (Italian). 3. P. Pavel, On the remainder of some Gaussian formulae, Studia Univ. Babeş-Bolyai Ser. Math.-Phys. (967), 65 7. 33. P. Pavel, On some quadrature formulae of Gaussian type, Studia Univ. Babeş-Bolyai Ser. Math.-Phys. 3 (96), 5 5 (Romanian). 3. P. Pavel, On the remainder of certain quadrature formulae of Gauss-Christoffel type, Studia Univ. Babeş-Bolyai Ser. Math.-Phys. 3 (96), 67 7 (Romanian). 35. F. Peherstorfer, On the remainder of Gaussian quadrature formulas for Bernstein-Szegő weight functions, Math. Comp. 6 (993), 37 35. 36. T. Popoviciu, Sur une généralisation de la formule d integration numérique de Gauss, Acad. R. P. Romîne Fil. Iaşi Stud. Cerc. Şti. 6 (955), 9 57. 37. T. Schira, The remainder term for analytic functions of Gauss-Lobatto quadratures, J. Comput. Appl. Math. 76 (996), 7 93. 3. T. Schira, The remainder term for analytic functions of symmetric Gaussian quadratures, Math. Comp. 66 (997), 97 3. 39. M. M. Spalević, Product of Turán quadratures for cube, simplex, surface of the sphere, E r n, En r, J. Comput. Appl. Math. 6 (999), 99 5.. M. M. Spalević, Calculation of Chakalov Popoviciu s quadratures of Radau and Lobatto type, ANZIAM J. (formerly as J. Aust. Math. Soc. B) 3 (), 9 7.. D. D. Stancu, On a class of orthogonal polynomials and on some general quadrature formulas with minimum number of terms, Bull. Math. Soc. Sci. Math. Phys. R.P. Romîne (N.S) (9) (957), 79 9.. D. D. Stancu, On certain general numerical integration formulas, Acad. R.P. Romîne. Stud. Cerc. Mat. 9 (95), 9 6 (Romanian). 3. A. H. Stroud and D. D. Stancu, Quadrature formulas with multiple Gaussian nodes, J. SIAM Numer. Anal. Ser. B (965), 9 3.. P. Turán, On the theory of the mechanical quadrature, Acta Sci. Math. Szeged (95), 3 37. 5. G. Vincenti, On the computation of the coefficients of s-orthogonal polynomials, SIAM J. Numer. Anal 3 (96), 9 9. Department of Mathematics, University of Niš, Faculty of Electronic Engineering, P.O. Box 73, Niš, Serbia Yugoslavia E-mail address: grade@gauss.elfak.ni.ac.yu Faculty of Science, Department of Mathematics and Informatics, P.O. Box 6, 3 Kragujevac, Serbia - Yugoslavia E-mail address: spale@knez.uis.kg.ac.yu