On the Chebyshev quadrature rule of finite part integrals
|
|
- Lillian Merritt
- 5 years ago
- Views:
Transcription
1 Journal of Computational and Applied Mathematics 18 25) On the Chebyshev quadrature rule of finite part integrals Philsu Kim a,,1, Soyoung Ahn b, U. Jin Choi b a Major in Mathematics, Dong-A University, 84, Hadan-2 Dong, Saha-Ku, Pusan , South Korea b Department of Mathematics, KAIST, 373-1, Guseong-dong, Yuseong-gu, Daejon 35-71, South Korea Received 5 October 24; received in revised form 16 October 24 Abstract The stability and the convergence of the Chebyshev quadrature rule of one-sided finite part integrals J. Approx. Theory ) ) is considered with some numerical examples. The three-term recurrence relation is suggested for computing the quadrature weights. 24 Elsevier B.V. All rights reserved. Keywords: Finite part integrals; Chebyshev interpolation 1. Introduction The following finite part integral defined by Qf ) is considered. Qf ) := fτ) wτ) dτ, 1 τ 1.1) wτ) = 1 τ) α 1 + τ) β, β > 1, < α <, α + β = ) and fτ) is assumed to be smooth on [, 1]. The integral is to be interpreted as the Hadamard finite part. For the basic properties of it, we refer to [4] and references therein. Corresponding author. Tel.: ; fax: addresses: kimps@donga.ac.kr P. Kim), ckskan@kaist.ac.kr S. Ahn), ujinchoi@numer.kaist.ac.kr U.J. Choi). 1 This work was supported by Korea Research Foundation Grant KRF-22-3-C /$ - see front matter 24 Elsevier B.V. All rights reserved. doi:1.116/j.cam
2 148 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) In [2], Criscuolo has developed a Gaussian type quadrature rule for 1.1) based on the Lagrange interpolation polynomial at the zeros of the th orthogonal polynomial with respect to the weight w, obtaining the optimal stability result O 2α ) and rectifying an analogous result proved by Monegato in [9]. Furthermore, in [2] the more general case α < is analyzed. In [4], Elliott has analyzed two algorithms for the approximate evaluation of Qf ) α is negative, not an integer. He showed that the quadrature sums are closely related to finite-part integrals of the Bernstein polynomials. Then, Delbourgo and Elliott have extended the result of [4] to the case of integer exponents see [3]). In [6], Elliott has developed a numerical evaluation scheme for approximating 1.1) with β = using the Chebyshev interpolation polynomial at the classical node points and obtained the stability result O 2α log ), is the degree of the Chebyshev interpolation polynomial. By subtracting out the singularity, 1.1) is equivalently defined by Qf ) = If)+ f1)q, 1.2) If) = fτ) f1) wτ) dτ, q = 1 τ wτ) 1 τ dτ. Using the ewton s binomial theorem, we write 1 + τ) β = 2 β 1 1 τ ) β = 2 β 2 j= a) j denotes the Pochhammer symbol β) j j! ) 1 τ j, 1.3) 2 a) j := aa + 1) a + j 1) = Γa + j), a) Γa) = ) From the definition of Hadamard finite-part integral defined in [4, 1.3)] and 1.4), we have 1 τ) j+α dτ = 2j+α j + α = 2j+α α) j α1 + α) j, j. 1.5) Combining 1.3) and 1.5), q becomes q = 2 α+β 1 α j= β) j α) j j!1 + α) j. Hence, from the definition of hypergeometric function given in [1, ) and )], we see that q = 2 α+β 1 α Fα, β; 1 + α; 1) = 2α+β Γα)Γβ + 1) Γα + β + 1). 1.6) In [7] the author developed the Chebyshev quadrature rule based on the approximation to f by the Chebyshev interpolation polynomial p f at τ j = cos πj/, which are the zeroes of the polynomial
3 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) T τ) T +1 τ), commonly named the practical abscissae, T k t) is the Chebyshev polynomial of the first kind of degree k. It is defined with the remainder term R f ) as follows: Qf ) = I f ) + f1)q + R f ), 1.7) I f ) = Ip f ) = j= fτ j )w j, w j = 2 k= T k τ j )J k, J k = IT k ). 1.8) ote that the summation with primes denotes a sum with the first and last terms halved. For very restricted cases, β = and < α < 2 1, [7] shows that the above rule has the same stability result as the Gaussian quadrature rule. In this paper, the three-term recurrence relation is suggested as the method for computing the quadrature weights in Section 2. Most of all, we continue the analysis of the stability factor Λ defined by Λ = wj j= and obtain the same stability as that of the Gaussian quadrature rule of [9] without any restriction for β and α through a much simpler idea than that of [7] see Theorem 4). In Section 3, the convergence rate of the remainder R f ) is improved see Theorem 5) and some numerical examples are given in the final section, Section ) 2. The stability analysis 2.1. Computation of the quadrature weights To begin with the analysis of the stability, we first concern a practical method for computing the quadrature weights wj in 1.8). First, we note that J = IT ) =. Using the following expression for T k τ): T k τ) = k j= k) j k) j 1/2) j j!2 j 1 τ)j. We obtain J k in terms of hypergeometric function as follows. α+β Γα)Γβ + 1) J k = IT k ) = 2 Γα + β + 1) g k 1), k =, 1,..., g = 1 and g k k 1) is defined by [ ] k k) j k) j α) j k, k, α; 1 g k = = 1/2) j= j j!α + β + 1) 3 F 2 1 j 2, α + β + 1.
4 15 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) Then, using the results of [1], we see that g k satisfies the following three-term recurrence relation: β + α + k + 1)g k+1 + 2α β 1)g k + β + α k + 1)g k = 2.1) with the initial values g = 1, g 1 = β α + 1 α + β + 1. Therefore, using the three-term recurrence relation 2.1), we can easily calculate the quadrature weights wj in 1.8). For the numerical stability of the recurrence relation of 2.1), we refer to see [1] Stability analysis Throughout this paper, we assume that is an even positive integer. For estimating Λ, we need the following identity for wj of 1.8) ) wj = )j 1 1 T τ) T +1 τ) 2 τ j 1 T τ) T +1 τ) 1 τ τ j τ wτ) dτ, which is proved in [7, 2.9)]. Hence we can split w j of 1.8) as wj = w j1 + w j2 + w j3, j 1, 2.2) wj1 = )j 1 wτ j )1 + τ j ) wj2 = )j 1 τ j ) U τ) τ dτ, τ j wτ)1 τ 2 ) wτ j )1 τ j )2 ) τ j τ U τ) dτ, w j3 = )j 2 J J +1 1 τ. j Here, U k τ) denotes the second kind Chebyshev polynomial of order k. Before the estimation of Λ,we consider the following lemma, which is a basis technique of the estimation for w j2 and w j3. Lemma 1. Assume that g L 1 [, π]) is increasing or decreasing) on [, θ] and decreasing or increasing) on [ θ, π], θ [i/)π, [i + 1)/]π] for a fixed i =,..., 1. Then we have π ) gθ) sin θ dθ π max gθ) sin θ + max gπ θ) sin θ + mi g, θ [,π/] θ [,π/] {, i =, 1, m i g = 3 max gθ), i = 1,..., 2. θ [i/,i+1)/]
5 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) Proof. We first denote by A the integral in the left-hand side of the above formula and note the integrand function highly oscillates. For θ [i/)π, [i + 1)/]π], i=,..., 1, let A be split as follows: [i+1)/]π π ) [i+1)/]π A = + gθ) sin θ dθ = A i 1 + Ai 2 Ai 3. i/)π i/)π ote that { A A = 2, i =, A 1, i = 1. First, by replacing the variable θ by θ + π/ in the integral A i 1 and taking the mean value of the newand old integrals, we can find that π/) ) [i+1)/]π 2A i 1 = i + gθ) sin θ dθ + C k, i/)π k= k+1)π)/ C k = gθ) g θ + π )) sin θ dθ. kπ/ Replacing the variable θ by θ + kπ/ in each integral C k and from monotonicity conditions of g,weget i π/ i C k sin θ g θ + kπ ) g θ + k + 1 ) k= π dθ k= π/ [i+1)/]π gθ) sin θ dθ + gθ) sin θ dθ. i/)π Therefore, π/ [i+1)/]π A i 1 gθ) sin θ dθ + gθ) sin θ dθ. i/)π Similarly, by replacing the variable θ by θ π/ in A i 2,wehave π/ [i+1)/]π A i 2 gπ θ) sin θ dθ + gθ) sin θ dθ. i/)π Finally, we complete the proof using the fact sin θ sin θ. For the estimation of wj2,1 j 1, we let g α,β θ,j)= γ α,βθ) γ α,β πj/) cos θ cos πj/, γ α,βθ) = sin 2α+2 θ 2 cos2β+2 θ )
6 152 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) Then for < α, β, g α,β π θ, j)= g β,α θ, j), g γ α,β θ) + sin θg α,β θ,j) α,β θ,j)= cos θ cos πj/, 2.4) the prime in g α,β θ,j) and γ α,β θ) denotes the derivative of the function about θ. We note that the equation γ α,β θ) + sin θg α,β θ,j)= has always πj/ as a solution in, π). Lemma 2. The equation γ α,β θ)+sin θg α,β θ,j)=has the other solution except πj/, say θ, only in the case β > in, π), and then there exists i=1,..., 2dependent of α, β such that θ [i/,i+1)/] and for enough large we have ) ) i 2α max θ [i/,i+1/)] g α,β θ,j) = g α,β θ,j) =O. 2.5) Proof. Putting fx)= 1 x) α x α β + β + 1)x 2 β + α + β + 2) cos 2 πj 2 + cos 2 πj ) 2α+2 πj πj β + 1) sin cos2β , we obtain the following formula for x = cos 2 θ/2: fx)= f cos 2 θ ) 2 ) x = cos θ cos πj/)γ α,β θ) + sin θg α,β θ,j)). sin θ Then we know that fx)= has a double root at x = cos 2 πj/2. From f x) = 1 x) α x β gx) cos 2 πj ) 2 x, gx) = 1 + α + β)2 + α + β)x β)1 + α + β)x + β1 + β)). It is easy to check that the equation gx) = has only one solution in, 1) only in the case β > and otherwise no solution in, 1). Hence, only in the case β >, fx) has two extreme points in, 1) and the other simple root x = cos 2 θ/2, since f)< and limx 1 fx) =+ for β >. We know γ α,β θ) + sin θg α,β θ,j) =. Therefore, there exists i = 1,..., 2 dependent of α, β such that θ [i/,[i + 1)/]]. In particular, fx) is negative for cos 2 i + 1)/2 <x<cos 2 θ/2 and positive for cos 2 θ/2 <x<cos 2 i/2. Therefore, we know that g α,β θ,j)is a local maximum,
7 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) since fcos 2 θ/2)g α,β θ,j)>. From g α,β θ,j)= 1 2α θ θ 2 sin cos2β α + β + 2) sin 2 θ ) α + 1), for enough large, we get the desired behavior 2.5). Lemma 3. Assume that < α, β and 1 j 1. Then we have max θ [,π/] sin θ g α,β θ,j) =OM α,β), 2.6) ) j 2α+1 O 1 j ) ) 2β+1, < α 1 2, < β 1 2, ) ) j 2α+1 O, < α 1 M α,β = 2, 2 < β, O 1 j ) ) 2β+1, 1 2 < α, < β 1 2, O 1), 1 2 < α, 2 < β, and max θ [,π/] sin θ g α,β π θ,j) has the same behavior by 2.4). Proof. We first note that for any j = 1, 2,..., 1 and θ [, π/], sin θ θ πj/ 1, j = 1, 2,..., 1. cos θ cos πj/ Define ζ α θ,j)= sin 2α+2 θ/2 sin 2α+2 πj/2 θ πj/, ξ β θ,j)= cos 2β+2 θ/2 cos 2β+2 πj/2 θ πj/. Then the mean value theorem yields ξ β θ,j)= β + 1) cos 2β+1 χ 2 sin χ 1 j/) 2β+1, < β 1 c 1 2, 2 1, 1 2 < β, c 1 = β + 1), χ [θ, πj/] and θ [, π/]. On the other hand, if < α 1 2, sin2α+2 θ/2is increasing and convex. So, /πj sin 2α+2 πj/2 1 2 πj/2)2α+1. Hence, again the mean value theorem gives ζ α θ,j) { 12 πj/2) 2α+1, < α 1 2, α + 1, 1 2 < α.
8 154 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) Since sin θ g α,β θ,j) sin 2α+2 θ/2 cos 2β+2 θ/2 sin 2α+2 πj/2 cos 2β+2 πj/2 θ πj/ ζ πj α θ,j)cos2β πj ξ β θ,j)sin2α ζ α θ,j)ξ β θ,j) ζ π α θ,j) 2 πj ) 2β+2 πj + ξ β 2 θ,j) 2 ) πj θ ) 2α+2 + ζ α θ,j)ξ β θ,j)πj, summarizing the above inequalities, we obtain 2.6). Based on the above Lemmas 2 and 3, finally we prove the following asymptotic behavior of the stability Λ. Theorem 4. Let Λ be the stability factor defined in 1.9) of the quadrature rule 1.7). Then we have that Λ = O 2α ), < α <. 2.7) Proof. By using the facts [8, 2.2),2.8),2.1)] and the change of variable τ = cos θ, we first rewrite w j1 as follows: for 1 j 1, Since w j1 = )j 2 α+β+1 = 2α+β+2 2α πj πj sin cos2β πj πj sin2α cos2β π /2 k= sin θ cos πj/ cos θ dθ sin2k + 1)πj/). 2k + 1 /2 k= sin2k + 1)πj/) 2k + 1 πj/ = sin θ 2 sin θ dθ has maximum value at j = 1, we obtain the following estimation: j=1 wj1 2α+β+1 π/ = O 1 /2 j=1 sin θ sin θ dθ 2α πj πj sin cos2β j=1 ) 2β+1 ) j 2α + j=/2+1 1 j = O 2α ), < α <, 2.8) we used the Riemann lower sum.
9 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) For the estimation of wj2, by the change of variable τ = cos θ in w j2, we note that wj2 = 1 2 α+β+2 π 1 cos πj/) g α,β θ,j)sin θ dθ = O 1 j ) 2 π ) g α,β θ,j)sin θ dθ, g α,β θ,j)is defined in 2.3). Then using Lemma 1 and the asymptotic behaviors 2.5) and 2.6), we obtain j=1 w j2 =O 2α ). 2.9) For estimating w j3, we use the identity T τ) T +1 τ) = 2U τ)1 τ 2 ). Then J J +1 = 2 U τ)1 + τ)wτ) dτ = 2 π gθ) sin θ dθ, gθ) = 1 + cos θ) β+1 1 cos θ) α. ote that gθ) is a decreasing function. Thus, the result of Lemma 1 yields J J +1 =O 2α ). From this asymptotic behavior, we obtain wj3 = 1 J J τ = O 2α ). 2.1) j j=1 j=1 For estimating w, we need the following asymptotic behavior for J k: J k =Ok 2α ), k 1, 2.11) which is derived from the asymptotic analysis of [1] and the three-term recurrence relation 2.1). From 1.8) and 2.11), we get ) w 2 J k + J = O 1 k 2α = O 2α ) 2.12) and also, w = 1 k=1 k=1 U τ)wτ) dτ wτ) dτ = O1). 2.13) Hence, by summarizing the asymptotic behaviors 2.8) 2.1), 2.12) and 2.13), we complete the asymptotic behavior 2.7).
10 156 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) The convergence rate Using this Theorem 4 and recalling [7, Lemma 4.1], we improve the convergence rate of the rule as follows: Theorem 5. Let us consider the quadrature rule 1.7). Suppose the function fτ) possesses continuous derivatives up to order p 1 and the derivative f p) τ) satisfies Hölder continuity of order ρ. Then the remainder term R f ) satisfies R f ) =O p ρ 2α ). Proof. Let p τ) be a polynomial of degree satisfying max τ [,1] fτ) p τ) M p ρ with a constant M see [7, Lemma 4.1]). Putting r τ) = fτ) p τ), [7, Theorem 4 and Lemma 4.1] showthat R f ) Ir ) +C 1 2α p ρ, Ir ) = wτ) r τ) r 1) dτ 3.1) 1 τ for some constant C 1. For the estimation of Ir ), we use the following result of [5]: inf p sup x<y 1 r x) r y) x y λ = O p ρ+λ ) for < λ 1 and λ p + ρ. Thus, if we take λ > α, then Ir ) can be estimated as follows: since inf p Ir ) inf 1 τ) α+λ 1 + τ) β r τ) r 1) p 1 τ λ dτ = O p ρ+λ ), 1 τ) α+λ 1 + τ) β dτ <. Then we choose λ as small as we can, viz. λ = α + ε for any sufficiently small ε >. Therefore, we obtain R f ) =O p ρ 2α + α+ε )) = O p ρ 2α ). Remark6. We remark that in [2] Criscuolo proved that the Gaussian type quadrature converges with order O p ρ ) under the same assumptions of Theorem 5, with ρ = 1 when < α 2 1. So that, we deduce that the method in [2] performs better than the present method. 4. Some numerical examples In order to test estimation 2.7), according to formulas 1.8) and 1.9), we calculate the values of Λ / 2α with α =.3 varying the node number from 8 to 2 by step size 2 and increasing the index β from.6 to.6 by step size.2. In Table 1, we list these values. For each fixed index β, Table 1 shows
11 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) Table 1 Evaluation of Λ / 2α with α =.3 β = 8 = 1 = 12 = 14 = 16 = 18 = Table R f ) with α =.9 β Present method Elliott [6] the 2nd the 3rd = = = = = = = that the quantities Λ / 2α are almost constant independent of = 2k, k = 4, 5,...,1. On the other hand, it also shows that the parameter index β hardly influences the behavior of Λ / 2α. We nowtest the rule Q f ) = I f ) + f1)q with the improper integral Qf ) = 1 τ) τ) 1.1 dτ = 2Γ2.1)Γ.9). To do this, we set wτ)=1 τ).9 1+τ) β,fτ)=1 τ) τ) 1.1 β with β=.5+.2k, k=, 1, 3, 4. In Table 2, we display 1 2 R f ) :=1 2 Qf ) Q f ), varying both and β. We also list the corresponding errors to the second and third quadrature rules in [6]. For these rules in [6], we take fτ) = 1 τ) τ) 1.1, wτ) = 1 τ).9. From Table 2, it appears that the present method is converging better than the second rule of [6], which is similar to our method and based on the Chebyshev interpolation polynomial with the classical node points cos2j )π/2), j =1,...,. The last column of Table 2 is the result obtained by the third rule of [6], which is based on the Muneer s approximation and the Chebyshev interpolation. The results show
12 158 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) Table R f ) with α =.9 Present method The 2nd rule of [6] 1 5 R f ) 1 5 R f ) log 1 5 R f ) = = = = = that this method is converging a little faster than the present method. We also see that the convergence of the present method does not depend on the parameter β. Finally, we conclude this section by the test of the rule Q f ) = I f ) + f1)q with the finite part integral Qf ) = 1 + τ) 2.1 Γ.9)Γ3.1) dτ = τ) Γ2.2) We let wτ) = 1 τ).9,fτ) = 1 + τ) 2.1.InTable 3, we display the error 1 5 R f ) :=1 5 Qf ) Q f ) and also list the corresponding errors to the second rule in [6]. Table 3 shows that our method s results are better than that of the method in [6] by log. Aswe remarked, it is caused from the fact that our method makes the growth of the stability factor in [6] to slow by log, and it means our method is more stable. Acknowledgements The authors would like to express their sincere gratitude to the reviewers for many helpful comments and valuable suggestions on the first draft of this paper. References [1] M. Abramowitz, I.A. Stegun, Handbook of Mathematical Functions, Applied Mathematics Series, vol. 55, ational Bureau of Standard, Washington, DC, [2] G. Criscuolo, On the numerical evaluation of certain Hadamard finite-part integrals, Suppl. Rend. Circ. Matem. Palermo Ser. II ) [3] D. Delbourgo, D. Elliott, On the approximate evaluation of Hadamard finite-part integrals, IMA J. umer. Anal ) [4] D. Elliott, An asymptotic analysis of two algorithms for certain Hadamard finite-part integrals, IMA J. umer. Anal ) [5] D. Elliott, On the Hölder semi-norm of the remainder in polynomial approximation, Bull. Austral. Math. Soc ) [6] D. Elliott, Three algorithms for Hadamard finite-part integrals and fractional derivatives, J. Comp. Appl. Math )
13 P. Kim et al. / Journal of Computational and Applied Mathematics 18 25) [7] P. Kim, A Chebyshev quadrature rule for one sided finite part integrals, J. Approx. Theory ) [8] P. Kim, U.J. Choi, A quadrature rule of interpolatory type for Cauchy integrals, J. Comp. Appl. Math ) [9] G. Monegato, The numerical evaluation of a 2-D Cauchy principal value integral arising in boundary integral equation methods, Math. Comp ) [1] R. Piessens, M. Branders, The evaluation and application of some modified moments, BIT )
On the convergence of interpolatory-type quadrature rules for evaluating Cauchy integrals
Journal of Computational and Applied Mathematics 49 () 38 395 www.elsevier.com/locate/cam On the convergence of interpolatory-type quadrature rules for evaluating Cauchy integrals Philsu Kim a;, Beong
More informationA quadrature rule of interpolatory type for Cauchy integrals
Journal of Computational and Applied Mathematics 126 (2000) 207 220 www.elsevier.nl/locate/cam A quadrature rule of interpolatory type for Cauchy integrals Philsu Kim, U. Jin Choi 1 Department of Mathematics,
More informationOn the remainder term of Gauss Radau quadratures for analytic functions
Journal of Computational and Applied Mathematics 218 2008) 281 289 www.elsevier.com/locate/cam On the remainder term of Gauss Radau quadratures for analytic functions Gradimir V. Milovanović a,1, Miodrag
More informationON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT
ON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT Received: 31 July, 2008 Accepted: 06 February, 2009 Communicated by: SIMON J SMITH Department of Mathematics and Statistics La Trobe University,
More informationBivariate Lagrange interpolation at the Padua points: The generating curve approach
Journal of Approximation Theory 43 (6) 5 5 www.elsevier.com/locate/jat Bivariate Lagrange interpolation at the Padua points: The generating curve approach Len Bos a, Marco Caliari b, Stefano De Marchi
More informationMean Convergence of Interpolation at Zeros of Airy Functions
Mean Convergence of Interpolation at Zeros of Airy Functions D. S. Lubinsky Abstract The classical Erdős-Turán theorem established mean convergence of Lagrange interpolants at zeros of orthogonal polynomials.
More informationNumerische Mathematik
Numer. Math. (28) 9:43 65 DOI.7/s2-7-25-7 Numerische Mathematik The superconvergence of Newton Cotes rules for the Hadamard finite-part integral on an interval Jiming Wu Weiwei Sun Received: 2 May 27 /
More informationThe Hardy-Littlewood Function: An Exercise in Slowly Convergent Series
The Hardy-Littlewood Function: An Exercise in Slowly Convergent Series Walter Gautschi Department of Computer Sciences Purdue University West Lafayette, IN 4797-66 U.S.A. Dedicated to Olav Njåstad on the
More informationA derivative-free nonmonotone line search and its application to the spectral residual method
IMA Journal of Numerical Analysis (2009) 29, 814 825 doi:10.1093/imanum/drn019 Advance Access publication on November 14, 2008 A derivative-free nonmonotone line search and its application to the spectral
More informationRiemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,
Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson
Advanced Computational Fluid Dynamics AA5A Lecture Approximation Theory Antony Jameson Winter Quarter, 6, Stanford, CA Last revised on January 7, 6 Contents Approximation Theory. Least Squares Approximation
More informationComparison of Clenshaw-Curtis and Gauss Quadrature
WDS'11 Proceedings of Contributed Papers, Part I, 67 71, 211. ISB 978-8-7378-184-2 MATFYZPRESS Comparison of Clenshaw-Curtis and Gauss Quadrature M. ovelinková Charles University, Faculty of Mathematics
More informationξ,i = x nx i x 3 + δ ni + x n x = 0. x Dξ = x i ξ,i = x nx i x i x 3 Du = λ x λ 2 xh + x λ h Dξ,
1 PDE, HW 3 solutions Problem 1. No. If a sequence of harmonic polynomials on [ 1,1] n converges uniformly to a limit f then f is harmonic. Problem 2. By definition U r U for every r >. Suppose w is a
More informationConvergence and Stability of a New Quadrature Rule for Evaluating Hilbert Transform
Convergence and Stability of a New Quadrature Rule for Evaluating Hilbert Transform Maria Rosaria Capobianco CNR - National Research Council of Italy Institute for Computational Applications Mauro Picone,
More informationOn the solution of integral equations of the first kind with singular kernels of Cauchy-type
International Journal of Mathematics and Computer Science, 7(202), no. 2, 29 40 M CS On the solution of integral equations of the first kind with singular kernels of Cauchy-type G. E. Okecha, C. E. Onwukwe
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationOn Some Estimates of the Remainder in Taylor s Formula
Journal of Mathematical Analysis and Applications 263, 246 263 (2) doi:.6/jmaa.2.7622, available online at http://www.idealibrary.com on On Some Estimates of the Remainder in Taylor s Formula G. A. Anastassiou
More informationSpectral gradient projection method for solving nonlinear monotone equations
Journal of Computational and Applied Mathematics 196 (2006) 478 484 www.elsevier.com/locate/cam Spectral gradient projection method for solving nonlinear monotone equations Li Zhang, Weijun Zhou Department
More informationBURGESS INEQUALITY IN F p 2. Mei-Chu Chang
BURGESS INEQUALITY IN F p 2 Mei-Chu Chang Abstract. Let be a nontrivial multiplicative character of F p 2. We obtain the following results.. Given ε > 0, there is δ > 0 such that if ω F p 2\F p and I,
More informationA note on the bounds of the error of Gauss Turán-type quadratures
Journal of Computational and Applied Mathematic 2 27 276 282 www.elevier.com/locate/cam A note on the bound of the error of Gau Turán-type quadrature Gradimir V. Milovanović a, Miodrag M. Spalević b, a
More informationThe ABC of hyper recursions
Journal of Computational and Applied Mathematics 90 2006 270 286 www.elsevier.com/locate/cam The ABC of hyper recursions Amparo Gil a, Javier Segura a,, Nico M. Temme b a Departamento de Matemáticas, Estadística
More informationSome asymptotic properties of solutions for Burgers equation in L p (R)
ARMA manuscript No. (will be inserted by the editor) Some asymptotic properties of solutions for Burgers equation in L p (R) PAULO R. ZINGANO Abstract We discuss time asymptotic properties of solutions
More informationA Nodal Spline Collocation Method for the Solution of Cauchy Singular Integral Equations 1
European Society of Computational Methods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics (JNAIAM) vol. 3, no. 3-4, 2008, pp. 211-220 ISSN 1790 8140
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 9, 1999, pp. 39-52. Copyright 1999,. ISSN 1068-9613. ETNA QUADRATURE FORMULAS FOR RATIONAL FUNCTIONS F. CALA RODRIGUEZ, P. GONZALEZ VERA, AND M. JIMENEZ
More informationTHE ZERO-DISTRIBUTION AND THE ASYMPTOTIC BEHAVIOR OF A FOURIER INTEGRAL. Haseo Ki and Young One Kim
THE ZERO-DISTRIBUTION AND THE ASYMPTOTIC BEHAVIOR OF A FOURIER INTEGRAL Haseo Ki Young One Kim Abstract. The zero-distribution of the Fourier integral Q(u)eP (u)+izu du, where P is a polynomial with leading
More informationarxiv: v1 [math.cv] 18 Aug 2015
arxiv:508.04376v [math.cv] 8 Aug 205 Saddle-point integration of C bump functions Steven G. Johnson, MIT Applied Mathematics Created November 23, 2006; updated August 9, 205 Abstract This technical note
More informationON CERTAIN CLASSES OF UNIVALENT MEROMORPHIC FUNCTIONS ASSOCIATED WITH INTEGRAL OPERATORS
TWMS J. App. Eng. Math. V.4, No.1, 214, pp. 45-49. ON CERTAIN CLASSES OF UNIVALENT MEROMORPHIC FUNCTIONS ASSOCIATED WITH INTEGRAL OPERATORS F. GHANIM 1 Abstract. This paper illustrates how some inclusion
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 35, pp. 148-163, 2009. Copyright 2009,. ISSN 1068-9613. ETNA SPHERICAL QUADRATURE FORMULAS WITH EQUALLY SPACED NODES ON LATITUDINAL CIRCLES DANIELA
More informationSOME RESULTS FOR MAX-PRODUCT OPERATORS VIA POWER SERIES METHOD. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXXXVII, 2 (208), pp. 9 98 9 SOME RESULTS FOR MAX-PRODUCT OPERATORS VIA POWER SERIES METHOD T. YURDAKADIM and E. TAŞ Abstract. In this paper, we obtain an approximation
More informationMath Homework 2
Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is
More informationSimultaneous Gaussian quadrature for Angelesco systems
for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced
More informationInequalities for the numerical radius, the norm and the maximum of the real part of bounded linear operators in Hilbert spaces
Available online at www.sciencedirect.com Linear Algebra its Applications 48 008 980 994 www.elsevier.com/locate/laa Inequalities for the numerical radius, the norm the maximum of the real part of bounded
More informationON A SURVEY OF UNIFORM INTEGRABILITY OF SEQUENCES OF RANDOM VARIABLES
ON A SURVEY OF UNIFORM INTEGRABILITY OF SEQUENCES OF RANDOM VARIABLES S. A. Adeosun and S. O. Edeki 1 Department of Mathematical Sciences, Crescent University, Abeokuta, Nigeria. 2 Department of Industrial
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Mathematics 233 2010) 1554 1561 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: wwwelseviercom/locate/cam
More informationBackground and Definitions...2. Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and Functions...
Legendre Polynomials and Functions Reading Problems Outline Background and Definitions...2 Definitions...3 Theory...4 Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and
More informationequations in the semiclassical regime
Adaptive time splitting for nonlinear Schrödinger equations in the semiclassical regime Local error representation and a posteriori error estimator Workshop, Numerical Analysis of Evolution Equations Vill,
More informationDirect method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
Journal of Computational and Applied Mathematics 22 (28) 51 57 wwwelseviercom/locate/cam Direct method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
More informationDiscrete uniform limit law for additive functions on shifted primes
Nonlinear Analysis: Modelling and Control, Vol. 2, No. 4, 437 447 ISSN 392-53 http://dx.doi.org/0.5388/na.206.4. Discrete uniform it law for additive functions on shifted primes Gediminas Stepanauskas,
More informationJournal of Mathematical Analysis and Applications
J. Math. Anal. Appl. 35 29) 276 282 Contents lists available at ScienceDirect Journal of Mathematical Analysis and Applications www.elsevier.com/locate/jmaa A Turán-type inequality for the gamma function
More informationIterative algorithms based on the hybrid steepest descent method for the split feasibility problem
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (206), 424 4225 Research Article Iterative algorithms based on the hybrid steepest descent method for the split feasibility problem Jong Soo
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationMarkov Sonin Gaussian rule for singular functions
Journal of Computational and Applied Mathematics 169 (2004) 197 212 www.elsevier.com/locate/cam Markov Sonin Gaussian rule for singular functions G.Mastroianni, D.Occorsio Dipartimento di Matematica, Universita
More information= 1 2 x (x 1) + 1 {x} (1 {x}). [t] dt = 1 x (x 1) + O (1), [t] dt = 1 2 x2 + O (x), (where the error is not now zero when x is an integer.
Problem Sheet,. i) Draw the graphs for [] and {}. ii) Show that for α R, α+ α [t] dt = α and α+ α {t} dt =. Hint Split these integrals at the integer which must lie in any interval of length, such as [α,
More informationExplicit polynomial expansions of regular real functions by means of even order Bernoulli polynomials and boundary values
Journal of Computational and Applied Mathematics 176 (5 77 9 www.elsevier.com/locate/cam Explicit polynomial expansions of regular real functions by means of even order Bernoulli polynomials and boundary
More informationA convolution test equation for double delay integral equations
Journal of Computational and Applied Mathematics 228 (2009) 589 599 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationThe generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations.
The generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations. Johannes Tausch Abstract An extension of the Euler-Maclaurin formula to singular integrals was introduced
More informationRemarks on Bronštein s root theorem
Remarks on Bronštein s root theorem Guy Métivier January 23, 2017 1 Introduction In [Br1], M.D.Bronštein proved that the roots of hyperbolic polynomials (1.1) p(t, τ) = τ m + m p k (t)τ m k. which depend
More informationThe Expansion of the Confluent Hypergeometric Function on the Positive Real Axis
Applied Mathematical Sciences, Vol. 12, 2018, no. 1, 19-26 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.712351 The Expansion of the Confluent Hypergeometric Function on the Positive Real
More informationNON-EXTINCTION OF SOLUTIONS TO A FAST DIFFUSION SYSTEM WITH NONLOCAL SOURCES
Electronic Journal of Differential Equations, Vol. 2016 (2016, No. 45, pp. 1 5. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu NON-EXTINCTION OF
More informationIMPROVEMENTS OF COMPOSITION RULE FOR THE CANAVATI FRACTIONAL DERIVATIVES AND APPLICATIONS TO OPIAL-TYPE INEQUALITIES
Dynamic Systems and Applications ( 383-394 IMPROVEMENTS OF COMPOSITION RULE FOR THE CANAVATI FRACTIONAL DERIVATIVES AND APPLICATIONS TO OPIAL-TYPE INEQUALITIES M ANDRIĆ, J PEČARIĆ, AND I PERIĆ Faculty
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationA COMPLETELY MONOTONIC FUNCTION INVOLVING THE TRI- AND TETRA-GAMMA FUNCTIONS
ao DOI:.2478/s275-3-9-2 Math. Slovaca 63 (23), No. 3, 469 478 A COMPLETELY MONOTONIC FUNCTION INVOLVING THE TRI- AND TETRA-GAMMA FUNCTIONS Bai-Ni Guo* Jiao-Lian Zhao** Feng Qi* (Communicated by Ján Borsík
More informationSome Approximation Properties of Szasz-Mirakyan-Bernstein Operators
EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol 7, No 4, 04, 49-48 ISSN 307-5543 wwwejpamcom Some Approximation Properties of Szasz-Mirakyan-Bernstein Operators Tuncay Tunç, Ersin Şimşek, Department
More informationLARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS. S. G. Bobkov and F. L. Nazarov. September 25, 2011
LARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS S. G. Bobkov and F. L. Nazarov September 25, 20 Abstract We study large deviations of linear functionals on an isotropic
More informationAN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION
J. KSIAM Vol.19, No.4, 409 416, 2015 http://dx.doi.org/10.12941/jksiam.2015.19.409 AN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION MORAN KIM 1 AND CHOHONG MIN
More informationApplied Mathematics Letters
Applied Mathematics Letters 25 (2012) 974 979 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml On dual vector equilibrium problems
More informationINTEGRAL FORMULAS FOR CHEBYSHEV POLYNOMIALS AND THE ERROR TERM OF INTERPOLATORY QUADRATURE FORMULAE FOR ANALYTIC FUNCTIONS
MATHEMATICS OF COMPUTATION Volume 75 Number 255 July 2006 Pages 27 23 S 0025-578(06)0859-X Article electronically published on May 2006 INTEGRAL FORMULAS FOR CHEBYSHEV POLYNOMIALS AND THE ERROR TERM OF
More informationSZEGÖ ASYMPTOTICS OF EXTREMAL POLYNOMIALS ON THE SEGMENT [ 1, +1]: THE CASE OF A MEASURE WITH FINITE DISCRETE PART
Georgian Mathematical Journal Volume 4 (27), Number 4, 673 68 SZEGÖ ASYMPOICS OF EXREMAL POLYNOMIALS ON HE SEGMEN [, +]: HE CASE OF A MEASURE WIH FINIE DISCREE PAR RABAH KHALDI Abstract. he strong asymptotics
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationCh. 03 Numerical Quadrature. Andrea Mignone Physics Department, University of Torino AA
Ch. 03 Numerical Quadrature Andrea Mignone Physics Department, University of Torino AA 2017-2018 Numerical Quadrature In numerical analysis quadrature refers to the computation of definite integrals. y
More informationQuadrature Formula for Computed Tomography
Quadrature Formula for Computed Tomography orislav ojanov, Guergana Petrova August 13, 009 Abstract We give a bivariate analog of the Micchelli-Rivlin quadrature for computing the integral of a function
More informationApplied Mathematics Letters. A reproducing kernel method for solving nonlocal fractional boundary value problems
Applied Mathematics Letters 25 (2012) 818 823 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml A reproducing kernel method for
More informationA NEW CLASS OF MEROMORPHIC FUNCTIONS RELATED TO CHO-KWON-SRIVASTAVA OPERATOR. F. Ghanim and M. Darus. 1. Introduction
MATEMATIQKI VESNIK 66, 1 14, 9 18 March 14 originalni nauqni rad research paper A NEW CLASS OF MEROMORPHIC FUNCTIONS RELATED TO CHO-KWON-SRIVASTAVA OPERATOR F. Ghanim and M. Darus Abstract. In the present
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematics ON A HYBRID FAMILY OF SUMMATION INTEGRAL TYPE OPERATORS VIJAY GUPTA AND ESRA ERKUŞ School of Applied Sciences Netaji Subhas Institute of Technology
More informationTAYLOR SERIES FOR THE ASKEY-WILSON OPERATOR AND CLASSICAL SUMMATION FORMULAS arxiv:math/ v2 [math.ca] 26 Aug 2005
TAYLOR SERIES FOR THE ASKEY-WILSON OPERATOR AND CLASSICAL SUMMATION FORMULAS arxiv:math/0348v [math.ca] 6 Aug 005 BERNARDO LÓPEZ, JOSÉ MANUEL MARCO AND JAVIER PARCET Abstract. An analogue of Taylor s formula,
More informationInternational Journal of Pure and Applied Sciences and Technology
Int. J. Pure Appl. Sci. Technol., 15(2) (2013), pp. 68-72 International Journal of Pure and Applied Sciences and Technology ISSN 2229-6107 Available online at www.ijopaasat.in Research Paper Generalization
More informationA gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
More informationA NEW COLLOCATION METHOD FOR SOLVING CERTAIN HADAMARD FINITE-PART INTEGRAL EQUATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 6, Number, Pages 4 54 c 9 Institute for Scientific Computing and Information A NEW COLLOCATION METHOD FOR SOLVING CERTAIN HADAMARD FINITE-PART
More informationPart II NUMERICAL MATHEMATICS
Part II NUMERICAL MATHEMATICS BIT 31 (1991). 438-446. QUADRATURE FORMULAE ON HALF-INFINITE INTERVALS* WALTER GAUTSCHI Department of Computer Sciences, Purdue University, West Lafayette, IN 47907, USA Abstract.
More informationRiemann integral and Jordan measure are generalized to unbounded functions. is a Jordan measurable set, and its volume is a Riemann integral, R n
Tel Aviv University, 214/15 Analysis-III,IV 161 1 Improper integral 1a What is the problem................ 161 1b Positive integrands................. 162 1c Newton potential.................. 166 1d Special
More informationA Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations
Acta Appl Math (21) 19: 861 873 DOI 1.17/s144-8-9351-y A Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations M. Rasty M. Hadizadeh Received:
More informationSharp estimates for a class of hyperbolic pseudo-differential equations
Results in Math., 41 (2002), 361-368. Sharp estimates for a class of hyperbolic pseudo-differential equations Michael Ruzhansky Abstract In this paper we consider the Cauchy problem for a class of hyperbolic
More informationOn Newton-type methods with cubic convergence
Journal of Computational and Applied Mathematics 176 (2005) 425 432 www.elsevier.com/locate/cam On Newton-type methods with cubic convergence H.H.H. Homeier a,b, a Science + Computing Ag, IT Services Muenchen,
More informationJacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation
Wei et al. SpringerPlus (06) 5:70 DOI 0.86/s40064-06-3358-z RESEARCH Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation Open Access
More informationInterpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials
Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials Lawrence A. Harris Abstract. We extend the definition of Geronimus nodes to include pairs of real numbers where
More informationMarkov s Inequality for Polynomials on Normed Linear Spaces Lawrence A. Harris
New Series Vol. 16, 2002, Fasc. 1-4 Markov s Inequality for Polynomials on Normed Linear Spaces Lawrence A. Harris This article is dedicated to the 70th anniversary of Acad. Bl. Sendov It is a longstanding
More informationNONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian UNDER NONHOMOGENEOUS NEUMANN BOUNDARY CONDITION
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 210, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian
More informationOn the Lebesgue constant of subperiodic trigonometric interpolation
On the Lebesgue constant of subperiodic trigonometric interpolation Gaspare Da Fies and Marco Vianello November 4, 202 Abstract We solve a recent conjecture, proving that the Lebesgue constant of Chebyshev-like
More informationNonnegative linearization for little q-laguerre polynomials and Faber basis
Journal of Computational and Applied Mathematics 199 (2007) 89 94 www.elsevier.com/locate/cam Nonnegative linearization for little q-laguerre polynomials and Faber basis Josef Obermaier a,, Ryszard Szwarc
More informationON BLEIMANN-BUTZER-HAHN OPERATORS FOR EXPONENTIAL FUNCTIONS ULRICH ABEL AND MIRCEA IVAN
BULL. AUSTRAL. MATH. SOC. VOL. 75 (2007) [409-415] 41A36, 05A19, 05A20 ON BLEIMANN-BUTZER-HAHN OPERATORS FOR EXPONENTIAL FUNCTIONS ULRICH ABEL AND MIRCEA IVAN Some inequalities involving the binomial coefficients
More informationSTRONG CONVERGENCE OF AN IMPLICIT ITERATION PROCESS FOR ASYMPTOTICALLY NONEXPANSIVE IN THE INTERMEDIATE SENSE MAPPINGS IN BANACH SPACES
Kragujevac Journal of Mathematics Volume 36 Number 2 (2012), Pages 237 249. STRONG CONVERGENCE OF AN IMPLICIT ITERATION PROCESS FOR ASYMPTOTICALLY NONEXPANSIVE IN THE INTERMEDIATE SENSE MAPPINGS IN BANACH
More informationZERO DISTRIBUTION OF POLYNOMIALS ORTHOGONAL ON THE RADIAL RAYS IN THE COMPLEX PLANE* G. V. Milovanović, P. M. Rajković and Z. M.
FACTA UNIVERSITATIS (NIŠ) Ser. Math. Inform. 12 (1997), 127 142 ZERO DISTRIBUTION OF POLYNOMIALS ORTHOGONAL ON THE RADIAL RAYS IN THE COMPLEX PLANE* G. V. Milovanović, P. M. Rajković and Z. M. Marjanović
More informationExistence and multiple solutions for a second-order difference boundary value problem via critical point theory
J. Math. Anal. Appl. 36 (7) 511 5 www.elsevier.com/locate/jmaa Existence and multiple solutions for a second-order difference boundary value problem via critical point theory Haihua Liang a,b,, Peixuan
More informationComplete monotonicity of a function involving the p-psi function and alternative proofs
Global Journal of Mathematical Analysis, 2 (3) (24) 24-28 c Science Publishing Corporation www.sciencepubco.com/index.php/gjma doi:.449/gjma.v2i3.396 Research Paper Complete monotonicity of a function
More informationRate of convergence for certain families of summation integral type operators
J Math Anal Appl 296 24 68 618 wwwelseviercom/locate/jmaa Rate of convergence for certain families of summation integral type operators Vijay Gupta a,,mkgupta b a School of Applied Sciences, Netaji Subhas
More informationINTEGRABILITY OF SUPERHARMONIC FUNCTIONS IN A JOHN DOMAIN. Hiroaki Aikawa
INTEGRABILITY OF SUPERHARMONIC FUNCTIONS IN A OHN OMAIN Hiroaki Aikawa Abstract. The integrability of positive erharmonic functions on a bounded fat ohn domain is established. No exterior conditions are
More informationOn the least values of L p -norms for the Kontorovich Lebedev transform and its convolution
Journal of Approimation Theory 131 4 31 4 www.elsevier.com/locate/jat On the least values of L p -norms for the Kontorovich Lebedev transform and its convolution Semyon B. Yakubovich Department of Pure
More informationMathematical Tripos Part III Michaelmas 2017 Distribution Theory & Applications, Example sheet 1 (answers) Dr A.C.L. Ashton
Mathematical Tripos Part III Michaelmas 7 Distribution Theory & Applications, Eample sheet (answers) Dr A.C.L. Ashton Comments and corrections to acla@damtp.cam.ac.uk.. Construct a non-zero element of
More informationEXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69 4 (217 271 28 December 217 research paper originalni nauqni rad EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM Saeid Shokooh and Ghasem A.
More informationNumerical Integration for Multivariable. October Abstract. We consider the numerical integration of functions with point singularities over
Numerical Integration for Multivariable Functions with Point Singularities Yaun Yang and Kendall E. Atkinson y October 199 Abstract We consider the numerical integration of functions with point singularities
More informationMeasuring Ellipsoids 1
Measuring Ellipsoids 1 Igor Rivin Temple University 2 What is an ellipsoid? E = {x E n Ax = 1}, where A is a non-singular linear transformation of E n. Remark that Ax = Ax, Ax = x, A t Ax. The matrix Q
More informationDETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL
More informationPseudo-Differential Operators associated with Bessel Type Operators - I
Thai Journal of Mathematics Volume 8 (21 Number 1 : 51 62 www.math.science.cmu.ac.th/thaijournal Online ISSN 1686-29 Pseudo-Differential Operators associated with Bessel Type Operators - I B. B. Waphare
More informationThe Gauss hypergeometric function F (a, b; c; z) for large c
The Gauss hypergeometric function F a, b; c; z) for large c Chelo Ferreira, José L. López 2 and Ester Pérez Sinusía 2 Departamento de Matemática Aplicada Universidad de Zaragoza, 5003-Zaragoza, Spain.
More information(x x 0 )(x x 1 )... (x x n ) (x x 0 ) + y 0.
> 5. Numerical Integration Review of Interpolation Find p n (x) with p n (x j ) = y j, j = 0, 1,,..., n. Solution: p n (x) = y 0 l 0 (x) + y 1 l 1 (x) +... + y n l n (x), l k (x) = n j=1,j k Theorem Let
More informationConstructing optimal polynomial meshes on planar starlike domains
Constructing optimal polynomial meshes on planar starlike domains F. Piazzon and M. Vianello 1 Dept. of Mathematics, University of Padova (Italy) October 13, 2014 Abstract We construct polynomial norming
More informationSEMI-INNER PRODUCTS AND THE NUMERICAL RADIUS OF BOUNDED LINEAR OPERATORS IN HILBERT SPACES
SEMI-INNER PRODUCTS AND THE NUMERICAL RADIUS OF BOUNDED LINEAR OPERATORS IN HILBERT SPACES S.S. DRAGOMIR Abstract. The main aim of this paper is to establish some connections that exist between the numerical
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 39, pp. 1-11, 1. Copyright 1,. ISSN 168-9613. ETNA TRIGONOMETRIC GAUSSIAN QUADRATURE ON SUBINTERVALS OF THE PERIOD GASPARE DA FIES AND MARCO VIANELLO
More information