Numerische MathemalJk
|
|
- Millicent McLaughlin
- 6 years ago
- Views:
Transcription
1 Numer. Math. 44, 53-6 (1984) Numerische MathemalJk 9 Springer-Verlag 1984 Discrete Approximations to Spherically Symmetric Distributions* Dedicated to Fritz Bauer on the occasion of his 6th birthday Walter Gautschi Purdue University, Department of Computer Sciences, Mathematical Sciences Building, Room 442, West Lafayette, Indiana 4797, USA Summary. We consider the approximation of spherically symmetric distributions in ~ by linear combinations of Heaviside step functions or Dirac delta functions. The approximations are required to faithfully reproduce as many moments as possible. We discuss stable methods of computing such approximations, taking advantage of the close connection with Gauss- Christoffel quadrature. Numerical results are presented for the distributions of Maxwell, Bose-Einstein, and Fermi-Dirac. Subject Classifications: AMS(MOS): 65D15 CR: Introduction There is some interest among physicists in approximating the distribution functions of statistical mechanics by discrete functions - either linear combinations of Dirac delta functions or linear combinations of Heaviside step functions. For the Maxwell velocity distribution, Laframboise and Stauffer [9] and Calder, Laframboise and Stauffer [1] construct such approximations which are optimal in the sense of matching as many moments as possible. The resulting equations are solved in [9] by what amounts to Prony's method and in [1] by a reduction to an eigenvalue problem involving Hankel matrices. Both methods are classical; for the former, see, e.g. Hildebrand [8, w for the latter, Szeg6 [1, Eq. (2.2.9)]. They are subject to severe ill-conditioning, however, as is well-known. Here we point out that both approximation problems can be formulated in terms of Gauss-Christoffel quadrature, and can therefore be brought into the realm of stable modern methods of constructing orthogonal polynomials; see Gautschi [5]. In particular, algorithmic implementations of Christoffel's theorem (Galant [3], Gautschi [6]) find application here. We use these methods to generate numerical data for the distributions of Maxwell, Bose-Einstein, and Fermi-Dirac. * Work supported in part by the National Science Foundation under Grant MCS A1
2 54 W. Gautschi 2. Approximation by Step Functions We consider a function f which is spherically symmetric in F,. a, hence a function only of the radial distance, f=f(r), <r< ~. We impose the following conditions on f: (i) fe C 11-, (:x3] and f'(r)<o on [, ~]. (ii) The integrals ~f(r)r"dr, ~f'(r)r"dr, m=,1,2... all exist and are finite, o o In particular, f has to be nonnegative on [, ~]. This will not be required in this section. Noting that S f(r)rmdr= 1 f(r)rm+l]~ - 1 T f'(r)w+l dr, m=,1,2,..., o m+l m+l o it follows from (ii) that rm+lf(r) has a limit as r~, which of course must be as follows from (i) and (ii). Thus, equal to zero. Likewise, f(r)~o as r~, lim rmf(r) =, m =, 1, 2... (2.1) We wish to approximate f by a linear combination of Heaviside step functions, tl /(r)= ~ a~h(rv-r), (2.2) ~=1 where H(t)= if t<, H(t)=l if t>. Defining the moments of a function g(r), as in [1], [9], by Sg(r)r~dV, j=o, 1,2... where dv=[2nd/2/f(d/2)]ra-ldr is the volume element of the spherical shell in p,a, d>l, and dv=dr if d=l, we require that f and f have the same moments of orders up to 2n-l, i.e., co ~ avh(rv-r)rj+d-ldr=sf(r)r~+d-ldr, j=,1... 2n-l, v=l or, which is the same, rv av~ rj+a-ldr=~f(r)ri+a-'dr, v=l j=,1... 2n-1. Carrying out the integration on the left, and integrating by parts on the right, we get, upon using (2.1), d j co.. (a~rv)r~= ~ [-rnf'(r)]rjdr, j=, 1,.,2n-1. (2.3) v=l These are precisely the equations for n-point Gauss-Christoffel quadrature relative to the (nonnegative) integration measure d2(r)= -rdf'(r)dr on [, oo]. (2.4)
3 Discrete Approximations to Spherically Symmetric Distributions 55 Hence, r, in (2.2) are the Gaussian abscissas relative to the measure d2(r) in (2.4) and ;t~=avr ~ are the corresponding Christoffel numbers. Once the nodes rv and weights 2~ have been computed, the coefficients a v in (2.2) are simply obtained by a~ = 2,. r;- a, v = 1,2... n. If d= 1, the functions f(r) and f(r) are to be extended to negative values of r symmetrically with respect to the origin. To compute the n-point Gauss-Christoffel formula in question, it suffices to generate the coefficients ak, ilk' k =, 1... n - 1, in the recursion formula 7Ck+l(t):(t--O~k)7~k(t)--flkI~k_l(t), k=o, 1... n-l, (2.5) []i rc l(t)=o, rto(t)=l, for the respective (monic) orthogonal polynomials rck(')= nk(" ; d2). The desired nodes r~ are then the eigenvalues of the Jacobi matrix _ ] Jn = 1 ~1 r.. ' r O~n--1 3 oo and the weights 2v are expressible as 2~-flo -- v~,l, 2 flo = S d2(r), in terms of the o first components v~,l of the corresponding normalized eigenvectors. They all can be readily computed using standard software of linear algebra; see, e.g., Gautschi I-4, w 5.1]. 3. Examples Example3.1. The Maxwell distribution (cf. [9, 1])f(r)=n-a/2e -r2 on [, oc]. In this case, (2.4) yields o 2 dz(r)=~vffra+te-'~dr on [, oo]. (3.1) 7~' The recursion coefficients c~ k, fig in (2.5) that correspond to this measure can be computed in two different ways. We can start with the recursion coefficients co, flo for the measure d2~ on [, oe], which are available in Galant [2] for <k< 19 to 2 significant decimal digits, and then keep multiplying the measure by r (d + 1 times), each time generating the corresponding recursion coefficients by the algorithms described in Galant [3] or Gautschi [6]. The coefficient flo, at the end, is then adjusted to conform to the normalization adopted in (3.1). Alternatively, we may compute the ek, fig directly from (3.1), using the discretized Stieltjes procedure as described in Gautschi [5,
4 56 W. Gautschi Example 4.6]. We have used both these approaches, at the same time extending Galant's table up to k = 49 and recomputing it to 25 decimal places. Using double precision on the CDC 65 (ca. 29 significant decimal digits) we observed agreement to 25 decimal places. The results are listed in Table 1 of the Appendix. They can be used, as described at the end of Sect. 2, to produce Gauss-Christoffel formulae with as many as 5 terms. For reasons of space, we refrain from tabulating any of them. Example 3.2. The Bose-Einstein distribution f(r) = (fle r- 1)- 1 on [, oo], fi > 1, d > 2. The measure of interest now is d2(r)=flr a-2 e-'dr on [, oo]. (3.2) Since for large r the distribution behaves like d2(r)~fl-1 rde-rdr, we generated the recursion coefficients by the discretized Stieltjes procedure, using the Gauss-Laguerre quadrature rule to carry out the discretization (see [7] for a similar application). As a check, we also used a discretization based on the Fej6r quadrature rule applied to each subinterval of the decomposition [, oo] = [, 1] w [1, 1] w [1, 5] w [5, oo] (cf. [5, w In the case fl= 1, d=3 that we computed, the largest discrepancy observed was 1 unit in the 25th decimal place. The results are shown in Table 2 of the Appendix. Example 3.3. The Fermi-Dirac distribution f(r)=(fler2+ 1) 1 on [, oo], fl>. Here, rd+ 1 e--r 2 d2(r)=2fl (fl+e_r2)2 dr on [, ~]. For large r, this measure behaves similarly as the one in Example 3.1. Therefore, we used the same method as in Example 3.1 (the second one described there) to generate the recursion coefficients ~k, ilk' The results for fl= 1, d--3, are shown in Table 3 of the Appendix. 4. Approximation by Dirae Delta Functions We now assume that f~c[, oo], f(r)>o on [, oo], and that the integrals S f(r)r m dr, m =, 1, 2..., all exist and are finite, with f(r) dr >. Approximating f by a linear combination of Dirac delta functions, f(r),~f(r), f(r)= ~ a~(r-rv), (4.1) v=l and using the same moment-matching procedure as in Sect. 2, one is led immediately to the equations oo a-i j_ 1] (a,r,)r,- I [f(r)r"- #dr, j=o, n-1. (4.2) v=l
5 Discrete Approximations to Spherically Symmetric Distributions 57 Thus, r~ are the Gaussian abscissas relative to the (nonnegative) measure dp(r) =rd-lf(r)dr, and a~-14rv - 1-d, v=l,2... n, where Pv are the corresponding Christoffel numbers. For the Maxwell distribution of Example 3.1 this yields dp(r) =z-d/zra-1 e-'2dr, so that, when d=3, we can use the results of Example 3.1 relative to the case d=l, with an obvious modification of rio. For the Bose- Einstein distribution, one has d#(r)= r a- 1 (fie'-1)-1 dr, to which the numerical results of [7] apply, if d = 2 and fl--1. In other cases, one can easily adapt the methods used in Example 3.2. The same holds for the Fermi-Dirac distribution of Example 3.3. Appendix (Tables 1-3) see pages 58-6 References 1. Calder, A.C., Laframboise, J.G., Stauffer, A.D.: Optimum step-function approximation of the Maxwell distribution ~unpublished) oo 2. Galant, D.: Gauss quadrature rules for the evaluation of 2n ~S exp(-x:)f(x)dx. Math. Comput. 23 (1969), Review 42, Loose microfiche suppl. E o 3. Galant, D.: An implementation of Christoffel's theorem in the theory of orthogonal polynomials. Math. Comput (1971) 4. Gautschi, W.: A survey of Gauss-Christoffel quadrature formulae. In: E.B. Christoffel, The Influence of his Work on Mathematics and the Physical Sciences (P.L. Butzer, F. Feh6r, eds.), pp Basel: Birkh~iuser Gautschi, W.: On generating orthogonal polynomials. SIAM J. Sci. Statist. Comput. 3, (1982) 6. Gautschi, W.: An algorithmic implementation of the generalized Christoffel theorem. In: Numerische Integration (G. H~immerlin, ed.), pp Intern. Ser. Numer. Math. 57. Basel: Birkhauser Gautschi, W., Milovanovi6, G.V.: Gaussian quadrature involving Einstein and Fermi functions with an application to convergence acceleration of series, submitted for publication. 8. Hildebrand, F.B.: Introduction to Numerical Analysis, 2nd ed., New York: McGraw-Hill Laframboise, J.G., Stauffer, A.D.: Optimum discrete approximation of the Maxwell distribution. AIAA Journal 7, (1969) 1. Szeg6, G.: Orthogonal Polynomials, AMS Colloq. Publications, Vol. 23, 4th ed. 2nd printing. Providence, R.I.: Amer. Math. Soc Received July 5, 1983/September 7, 1983
6 58 W. Gautschi ~ o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o ~ o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o,,,,, ~++++++g~+ c~ II ~d o o o o o o o o o o o ~ o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o I I II o,...,,.., ~ ,,. O, ~C m ~ o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o VII 2, ~ o o o o o o o o o o ~ o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o oo 2 m,..,,,,,,,,.....,,. ~ ,, ~ 1 7 6
7 Discrete Approximations to Spherically Symmetric Distributions 59 o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o VII o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o ~ ~ ~ o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o ~ I I I ~ o o o o o o o o o o o o o o o o o ~ o o o o o ~ o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o m..,,, ~
8 6 W. Gautschi Table 3. Recursion coefficients ~k, ilk, O< k<49, for the orthogonal polynomials relative to d2(r)=2rn+l(1 +e-'~)-2e-'2dr on [, ~], d=3 k alpha(k) d d+O d d d d d d d+ 9 2.S d+O d+ II d d d d d d d d d d+O d d d S d+O d+OO d+O d d d+O d d d d d S d+O d+O d d d+O d d d Sd d d+O d d+ 48 S.89SO d+O d+ d=3 beta(k) d d d d d d SO d+O d+O 1, d IO id+O d d d d d d d+O d+ 3, d+ 3, d d+ 3, d d d d iI d d d d d d+O d d+O d+ 5, d d Sld+O iOlS d d d d d d+ 7, d+O d d d d d d+O
Numerische MathematJk
Numer. Math. 49, 111-121 (1986) Numerische MathematJk 9 Springer-Verlag 1986 Spline Approximations to Spherically Symmetric Distributions * Walter Gautschi I and Gradimir V. Milovanovi6 1 Department of
More informationSummation of Series and Gaussian Quadratures
Summation of Series Gaussian Quadratures GRADIMIR V. MILOVANOVIĆ Dedicated to Walter Gautschi on the occasion of his 65th birthday Abstract. In 985, Gautschi the author constructed Gaussian quadrature
More informationHOW AND HOW NOT TO CHECK GAUSSIAN QUADRATURE FORMULAE *
BIT 23 (1983), 209--216 HOW AND HOW NOT TO CHECK GAUSSIAN QUADRATURE FORMULAE * WALTER GAUTSCHI Department of Computer Sciences, Purdue University, West Lafayette, Indiana 47907, U.S.A. Abstract. We discuss
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 informationCOMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE
BIT 6-85//41-11 $16., Vol. 4, No. 1, pp. 11 118 c Swets & Zeitlinger COMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE WALTER GAUTSCHI Department of Computer Sciences,
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 informationVariable-precision recurrence coefficients for nonstandard orthogonal polynomials
Numer Algor (29) 52:49 418 DOI 1.17/s1175-9-9283-2 ORIGINAL RESEARCH Variable-precision recurrence coefficients for nonstandard orthogonal polynomials Walter Gautschi Received: 6 January 29 / Accepted:
More informationNumerische Mathematik
Numer. Math. 52, 241-250 (1988) Numerische Mathematik 9 Springer-Verlag 1988 Lower Bounds for the Condition Number of Vandermonde Matrices* Walter Gautschi 1,** and Gabriele Inglese e 1 Department of Computer
More information(1-2) fx\-*f(x)dx = ^^ Z f(x[n)) + R (f), a > - 1,
MATHEMATICS OF COMPUTATION, VOLUME 29, NUMBER 129 JANUARY 1975, PAGES 93-99 Nonexistence of Chebyshev-Type Quadratures on Infinite Intervals* By Walter Gautschi Dedicated to D. H. Lehmer on his 10th birthday
More informationIS THE RECURRENCE RELATION FOR ORTHOGONAL POLYNOMIALS ALWAYS STABLE?* WALTER GAUTSCHI
BIT33 (1993), 277-284. IS THE RECURRENCE RELATION FOR ORTHOGONAL POLYNOMIALS ALWAYS STABLE?* WALTER GAUTSCHI Dept. of Computer Sciences, Purdue University, West Lafayette, IN 47907-1398, USA Abstract.
More informationOn summation/integration methods for slowly convergent series
Stud. Univ. Babeş-Bolyai Math. 6(26), o. 3, 359 375 On summation/integration methods for slowly convergent series Gradimir V. Milovanović Dedicated to Professor Gheorghe Coman on the occasion of his 8th
More informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationGAUSS-LAGUERRE AND GAUSS-HERMITE QUADRATURE ON 64, 96 AND 128 NODES
GAUSS-LAGUERRE AND GAUSS-HERMITE QUADRATURE ON 64, 96 AND 128 NODES RICHARD J. MATHAR Abstract. The manuscript provides tables of abscissae and weights for Gauss- Laguerre integration on 64, 96 and 128
More informationNumerische Mathematik
Numer. Math. (2000) 86: 617 633 Digital Object Identifier (DOI) 10.1007/s002110000186 Numerische Mathematik Quadrature rules for rational functions Walter Gautschi 1, Laura Gori 2, M. Laura Lo Cascio 2
More informationNumerical integration formulas of degree two
Applied Numerical Mathematics 58 (2008) 1515 1520 www.elsevier.com/locate/apnum Numerical integration formulas of degree two ongbin Xiu epartment of Mathematics, Purdue University, West Lafayette, IN 47907,
More informationQuadrature Formulas for Infinite Integrals
Quadrature Formulas for Infinite Integrals By W. M. Harper 1. Introduction. Since the advent of high-speed computers, "mechanical" quadratures of the type (1) f w(x)f(x)dx^ Hif(aj) Ja 3=1 have become increasingly
More informationSzegő-Lobatto quadrature rules
Szegő-Lobatto quadrature rules Carl Jagels a,, Lothar Reichel b,1, a Department of Mathematics and Computer Science, Hanover College, Hanover, IN 47243, USA b Department of Mathematical Sciences, Kent
More information27r---./p. IR.(f)l<=,, l(o) max If(z)l. (1.2) o={z" z= 1/2(u + u-1), u p ei, O<- O <-- 2.a }
SIAM J. NUMER. ANAL. Vol. 27, No. 1, pp. 219-224, February 1990 1990 Society for Industrial and Applied Mathematics 015 A NOTE ON THE CONTOUR INTEGRAL REPRESENTATION OF THE REMAINDER TERM FOR A GAUSS-CHEBYSHEV
More informationA determinant characterization of moment sequences with finitely many mass-points
To appear in Linear and Multilinear Algebra Vol 00, No 00, Month 20XX, 1 9 A determinant characterization of moment sequences with finitely many mass-points Christian Berg a and Ryszard Szwarc b a Department
More informationGENERALIZED GAUSS RADAU AND GAUSS LOBATTO FORMULAE
BIT 0006-3835/98/3804-0101 $12.00 2000, Vol., No., pp. 1 14 c Swets & Zeitlinger GENERALIZED GAUSS RADAU AND GAUSS LOBATTO FORMULAE WALTER GAUTSCHI 1 1 Department of Computer Sciences, Purdue University,
More informationContemporary Mathematicians
Contemporary Mathematicians Joseph P.S. Kung University of North Texas, USA Editor For further volumes: http://www.springer.com/series/4817 Claude Brezinski Ahmed Sameh Editors Walter Gautschi, Volume
More informationarxiv: v1 [physics.comp-ph] 22 Jul 2010
Gaussian integration with rescaling of abscissas and weights arxiv:007.38v [physics.comp-ph] 22 Jul 200 A. Odrzywolek M. Smoluchowski Institute of Physics, Jagiellonian University, Cracov, Poland Abstract
More informationPositivity of Turán determinants for orthogonal polynomials
Positivity of Turán determinants for orthogonal polynomials Ryszard Szwarc Abstract The orthogonal polynomials p n satisfy Turán s inequality if p 2 n (x) p n 1 (x)p n+1 (x) 0 for n 1 and for all x in
More informationCALCULATION OF GAUSS-KRONROD QUADRATURE RULES. 1. Introduction A(2n+ 1)-point Gauss-Kronrod integration rule for the integral.
MATHEMATICS OF COMPUTATION Volume 66, Number 219, July 1997, Pages 1133 1145 S 0025-5718(97)00861-2 CALCULATION OF GAUSS-KRONROD QUADRATURE RULES DIRK P. LAURIE Abstract. The Jacobi matrix of the (2n+1)-point
More informationTHE CIRCLE THEOREM AND RELATED THEOREMS FOR GAUSS-TYPE QUADRATURE RULES. Dedicated to Ed Saff on the occasion of his 60th birthday
THE CIRCLE THEOREM AND RELATED THEOREMS FOR GAUSS-TYPE QUADRATURE RULES WALTER GAUTSCHI Dedicated to Ed Saff on the occasion of his 6th birthday Abstract. In 96, P.J. Davis and P. Rabinowitz established
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationVerified Calculation of the Nodes and Weights for Gaussian Quadrature Formulas
Interval Computations No 4, 1993 Verified Calculation of the Nodes and Weights for Gaussian Quadrature Formulas Ulrike Storck Using Gaussian Quadrature for the verified calculation of an integral assumes
More informationA Note on Extended Gaussian Quadrature
MATHEMATICS OF COMPUTATION, VOLUME 30, NUMBER 136 OCTOBER 1976, PAGES 812-817 A Note on Extended Gaussian Quadrature Rules By Giovanni Monegato* Abstract. Extended Gaussian quadrature rules of the type
More informationSpecial Functions of Mathematical Physics
Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table
More informationBulletin T.CXLV de l Académie serbe des sciences et des arts 2013 Classe des Sciences mathématiques et naturelles Sciences mathématiques, No 38
Bulletin T.CXLV de l Académie serbe des sciences et des arts 213 Classe des Sciences mathématiques et naturelles Sciences mathématiques, No 38 FAMILIES OF EULER-MACLAURIN FORMULAE FOR COMPOSITE GAUSS-LEGENDRE
More informationGaussian interval quadrature rule for exponential weights
Gaussian interval quadrature rule for exponential weights Aleksandar S. Cvetković, a, Gradimir V. Milovanović b a Department of Mathematics, Faculty of Mechanical Engineering, University of Belgrade, Kraljice
More informationMatrix functions that preserve the strong Perron- Frobenius property
Electronic Journal of Linear Algebra Volume 30 Volume 30 (2015) Article 18 2015 Matrix functions that preserve the strong Perron- Frobenius property Pietro Paparella University of Washington, pietrop@uw.edu
More informationGauss-type Quadrature Rules for Rational Functions
Gauss-type Quadrature Rules for Rational Functions arxiv:math/9307223v1 [math.ca] 20 Jul 1993 Walter Gautschi Abstract. When integrating functions that have poles outside the interval of integration, but
More informationƒ f(x)dx ~ X) ^i,nf(%i,n) -1 *=1 are the zeros of P«(#) and where the num
ZEROS OF THE HERMITE POLYNOMIALS AND WEIGHTS FOR GAUSS' MECHANICAL QUADRATURE FORMULA ROBERT E. GREENWOOD AND J. J. MILLER In the numerical integration of a function ƒ (x) it is very desirable to choose
More informationOrthogonal polynomials
Orthogonal polynomials Gérard MEURANT October, 2008 1 Definition 2 Moments 3 Existence 4 Three-term recurrences 5 Jacobi matrices 6 Christoffel-Darboux relation 7 Examples of orthogonal polynomials 8 Variable-signed
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationConstrained Leja points and the numerical solution of the constrained energy problem
Journal of Computational and Applied Mathematics 131 (2001) 427 444 www.elsevier.nl/locate/cam Constrained Leja points and the numerical solution of the constrained energy problem Dan I. Coroian, Peter
More informationOn orthogonal polynomials for certain non-definite linear functionals
On orthogonal polynomials for certain non-definite linear functionals Sven Ehrich a a GSF Research Center, Institute of Biomathematics and Biometry, Ingolstädter Landstr. 1, 85764 Neuherberg, Germany Abstract
More informationOptimally scaled and optimally conditioned Vandermonde and Vandermonde-like matrices
BIT Numer Math (20) 5: 03 25 DOI 0.007/s0543-00-0293- Optimally scaled and optimally conditioned Vandermonde and Vandermonde-like matrices Walter Gautschi Received: July 200 / Accepted: 6 October 200 /
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationNotes on Special Functions
Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.
More informationIntroduction to Numerical Analysis
J. Stoer R. Bulirsch Introduction to Numerical Analysis Translated by R. Bartels, W. Gautschi, and C. Witzgall Springer Science+Business Media, LLC J. Stoer R. Bulirsch Institut fiir Angewandte Mathematik
More informationA lower estimate for the Lebesgue constants of linear means of Laguerre expansions
A lower estimate for the Lebesgue constants of linear means of Laguerre expansions George Gasper 1 and Walter Trebels 2 Dedicated to P. L. Butzer on the occasion of his 70th birthday in gratitude (Oct.
More informationDiscrete Orthogonal Polynomials on Equidistant Nodes
International Mathematical Forum, 2, 2007, no. 21, 1007-1020 Discrete Orthogonal Polynomials on Equidistant Nodes Alfredo Eisinberg and Giuseppe Fedele Dip. di Elettronica, Informatica e Sistemistica Università
More informationNonlinear trigonometric approximation and the Dirac delta function
Journal of Computational and Applied Mathematics 9 (6) 34 45 www.elsevier.com/locate/cam Nonlinear trigonometric approximation and the Dirac delta function Xiubin Xu Department of Information and Computational
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationDisconjugacy and the Secant Conjecture
Disconjugacy and the Secant Conjecture Alexandre Eremenko November 26, 2014 Abstract We discuss the so-called secant conjecture in real algebraic geometry, and show that it follows from another interesting
More informationFAMILIES OF ALGEBRAIC TEST EQUATIONS W. GAUTSCHI (I)
FAMILIES OF ALGEBRAIC TEST EQUATIONS W. GAUTSCHI (I) 1. Introduction. Test matrices have been in use for some time to scrutinize computer algorithms for solving linear algebraic systems and eigenvalue
More informationA Note on the Recursive Calculation of Incomplete Gamma Functions
A Note on the Recursive Calculation of Incomplete Gamma Functions WALTER GAUTSCHI Purdue University It is known that the recurrence relation for incomplete gamma functions a n, x, 0 a 1, n 0,1,2,..., when
More informationWeighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 2, pp. 401 412 (2013) http://campus.mst.edu/adsa Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationCONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION
CONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION Tuncay Aktosun Department of Mathematics and Statistics Mississippi State University Mississippi State, MS 39762 Abstract: For the one-dimensional
More informationFully Symmetric Interpolatory Rules for Multiple Integrals over Infinite Regions with Gaussian Weight
Fully Symmetric Interpolatory ules for Multiple Integrals over Infinite egions with Gaussian Weight Alan Genz Department of Mathematics Washington State University ullman, WA 99164-3113 USA genz@gauss.math.wsu.edu
More informationml then either k + m < n or k + n < m and one of the above scalar
SIAM J. MATH. ANAL. Vol. 23, No. 4, pp. 959-964, July 1992 1992 Society for Industrial and Applied Mathematics 010 ORTHOGONAL POLYNOMIALS AND A DISCRETE BOUNDARY VALUE PROBLEM I* RYSZARD SZWARC" Abstract.
More informationCertain Dual Series Equations Involving Generalized Laguerre Polynomials
International Journal of Computational and Applied Mathematics. ISSN 1819-4966 Volume 11, Number 1 (2016), pp. 55-59 Research India Publications http://www.ripublication.com Certain Dual Series Equations
More informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More informationQuadratures and integral transforms arising from generating functions
Quadratures and integral transforms arising from generating functions Rafael G. Campos Facultad de Ciencias Físico-Matemáticas, Universidad Michoacana, 58060, Morelia, México. rcampos@umich.mx Francisco
More informationASYMPTOTIC DISTRIBUTION OF THE MAXIMUM CUMULATIVE SUM OF INDEPENDENT RANDOM VARIABLES
ASYMPTOTIC DISTRIBUTION OF THE MAXIMUM CUMULATIVE SUM OF INDEPENDENT RANDOM VARIABLES KAI LAI CHUNG The limiting distribution of the maximum cumulative sum 1 of a sequence of independent random variables
More informationExact Solution of the Dirac Equation with a Central Potential
Commun. math. Phys. 27, 155 161 (1972) by Springer-Verlag 1972 Exact Solution of the Dirac Equation with a Central Potential E. J. KANELLOPOULOS, TH. V. KANELLOPOULOS, and K. WILDERMUTH Institut fur Theoretische
More informationA trigonometric orthogonality with respect to a nonnegative Borel measure
Filomat 6:4 01), 689 696 DOI 10.98/FIL104689M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat A trigonometric orthogonality with
More informationGeometry for Physicists
Hung Nguyen-Schafer Jan-Philip Schmidt Tensor Analysis and Elementary Differential Geometry for Physicists and Engineers 4 i Springer Contents 1 General Basis and Bra-Ket Notation 1 1.1 Introduction to
More informationNEEDLET APPROXIMATION FOR ISOTROPIC RANDOM FIELDS ON THE SPHERE
NEEDLET APPROXIMATION FOR ISOTROPIC RANDOM FIELDS ON THE SPHERE Quoc T. Le Gia Ian H. Sloan Yu Guang Wang Robert S. Womersley School of Mathematics and Statistics University of New South Wales, Australia
More informationIntroduction to Numerical Analysis
J. Stoer R. Bulirsch Introduction to Numerical Analysis Second Edition Translated by R. Bartels, W. Gautschi, and C. Witzgall With 35 Illustrations Springer Contents Preface to the Second Edition Preface
More informationESTIMATES OF THE TRACE OF THE INVERSE OF A SYMMETRIC MATRIX USING THE MODIFIED CHEBYSHEV ALGORITHM
ESTIMATES OF THE TRACE OF THE INVERSE OF A SYMMETRIC MATRIX USING THE MODIFIED CHEBYSHEV ALGORITHM GÉRARD MEURANT In memory of Gene H. Golub Abstract. In this paper we study how to compute an estimate
More informationAsymptotics of Integrals of. Hermite Polynomials
Applied Mathematical Sciences, Vol. 4, 010, no. 61, 04-056 Asymptotics of Integrals of Hermite Polynomials R. B. Paris Division of Complex Systems University of Abertay Dundee Dundee DD1 1HG, UK R.Paris@abertay.ac.uk
More informationPOINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES
November 1, 1 POINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES FRITZ KEINERT AND SOON-GEOL KWON,1 Abstract Two-direction multiscaling functions φ and two-direction multiwavelets
More informationComplex symmetric operators
Complex symmetric operators Stephan Ramon Garcia 1 Complex symmetric operators This section is a brief introduction to complex symmetric operators, a certain class of Hilbert space operators which arise
More informationBounds for Eigenvalues of Tridiagonal Symmetric Matrices Computed. by the LR Method. By Gene H. Golub
Bounds for Eigenvalues of Tridiagonal Symmetric Matrices Computed by the LR Method By Gene H. Golub 1. Introduction. In recent years, a number of methods have been proposed for finding the eigenvalues
More informationYALE UNIVERSITY DEPARTMENT OF COMPUTER SCIENCE
fa it) On the Jacobi Polynomials PA ' Gregory Matviyenko Research Report YALETJ/DCS/RR-1076 June 13, 1995 DTI,* QUALITY DEFECTED 3 YALE UNIVERSITY DEPARTMENT OF COMPUTER SCIENCE & k.-. *,«.
More informationOpen Research Online The Open University s repository of research publications and other research outputs
Open Research Online The Open University s repository of research publications and other research outputs A note on the Weiss conjecture Journal Item How to cite: Gill, Nick (2013). A note on the Weiss
More informationSINGULAR VALUE DECOMPOSITION
DEPARMEN OF MAHEMAICS ECHNICAL REPOR SINGULAR VALUE DECOMPOSIION Andrew Lounsbury September 8 No 8- ENNESSEE ECHNOLOGICAL UNIVERSIY Cookeville, N 38 Singular Value Decomposition Andrew Lounsbury Department
More informationarxiv: v1 [nlin.cd] 23 Oct 2009
Lyapunov exponent and natural invariant density determination of chaotic maps: An iterative maximum entropy ansatz arxiv:0910.4561v1 [nlin.cd] 23 Oct 2009 1. Introduction Parthapratim Biswas Department
More informationEnergy spectrum for a short-range 1/r singular potential with a nonorbital barrier using the asymptotic iteration method
Energy spectrum for a short-range 1/r singular potential with a nonorbital barrier using the asymptotic iteration method A. J. Sous 1 and A. D. Alhaidari 1 Al-Quds Open University, Tulkarm, Palestine Saudi
More informationLinear Algebra Practice Problems
Math 7, Professor Ramras Linear Algebra Practice Problems () Consider the following system of linear equations in the variables x, y, and z, in which the constants a and b are real numbers. x y + z = a
More informationGauss Hermite interval quadrature rule
Computers and Mathematics with Applications 54 (2007) 544 555 www.elsevier.com/locate/camwa Gauss Hermite interval quadrature rule Gradimir V. Milovanović, Alesandar S. Cvetović Department of Mathematics,
More informationSpectra of some self-exciting and mutually exciting point processes
Biometrika (1971), 58, 1, p. 83 83 Printed in Great Britain Spectra of some self-exciting and mutually exciting point processes BY ALAN G. HAWKES University of Durham SUMMAKY In recent years methods of
More informationInstructions for Matlab Routines
Instructions for Matlab Routines August, 2011 2 A. Introduction This note is devoted to some instructions to the Matlab routines for the fundamental spectral algorithms presented in the book: Jie Shen,
More informationTheory and Problems of Signals and Systems
SCHAUM'S OUTLINES OF Theory and Problems of Signals and Systems HWEI P. HSU is Professor of Electrical Engineering at Fairleigh Dickinson University. He received his B.S. from National Taiwan University
More informationA NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS. Masaru Nagisa. Received May 19, 2014 ; revised April 10, (Ax, x) 0 for all x C n.
Scientiae Mathematicae Japonicae Online, e-014, 145 15 145 A NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS Masaru Nagisa Received May 19, 014 ; revised April 10, 014 Abstract. Let f be oeprator monotone
More informationarxiv: v1 [math.na] 5 May 2011
ITERATIVE METHODS FOR COMPUTING EIGENVALUES AND EIGENVECTORS MAYSUM PANJU arxiv:1105.1185v1 [math.na] 5 May 2011 Abstract. We examine some numerical iterative methods for computing the eigenvalues and
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationGeneralized Functions Theory and Technique Second Edition
Ram P. Kanwal Generalized Functions Theory and Technique Second Edition Birkhauser Boston Basel Berlin Contents Preface to the Second Edition x Chapter 1. The Dirac Delta Function and Delta Sequences 1
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 informationIndependence of some multiple Poisson stochastic integrals with variable-sign kernels
Independence of some multiple Poisson stochastic integrals with variable-sign kernels Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological
More informationDeviation Measures and Normals of Convex Bodies
Beiträge zur Algebra und Geometrie Contributions to Algebra Geometry Volume 45 (2004), No. 1, 155-167. Deviation Measures Normals of Convex Bodies Dedicated to Professor August Florian on the occasion
More informationAN ITERATIVE METHOD WITH ERROR ESTIMATORS
AN ITERATIVE METHOD WITH ERROR ESTIMATORS D. CALVETTI, S. MORIGI, L. REICHEL, AND F. SGALLARI Abstract. Iterative methods for the solution of linear systems of equations produce a sequence of approximate
More informationSection 6.6 Gaussian Quadrature
Section 6.6 Gaussian Quadrature Key Terms: Method of undetermined coefficients Nonlinear systems Gaussian quadrature Error Legendre polynomials Inner product Adapted from http://pathfinder.scar.utoronto.ca/~dyer/csca57/book_p/node44.html
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 informationA CLASS OF SLOWLY CONVERGENT SEMES AND THEIR SUMMATION BY GAUSSIAN QUADRATURE WALTER GAUTSCHI
mathematics of computation volume 57, number 195 july 1991, pages 309-325 A CLASS OF SLOWLY CONVERGENT SEMES AND THEIR SUMMATION BY GAUSSIAN QUADRATURE WALTER GAUTSCHI Abstract. Series are considered whose
More information1. Introduction. Let µ(t) be a distribution function with infinitely many points of increase in the interval [ π, π] and let
SENSITIVITY ANALYSIS OR SZEGŐ POLYNOMIALS SUN-MI KIM AND LOTHAR REICHEL Abstract Szegő polynomials are orthogonal with respect to an inner product on the unit circle Numerical methods for weighted least-squares
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationA Riemann Lebesgue Lemma for Jacobi expansions
A Riemann Lebesgue Lemma for Jacobi expansions George Gasper 1 and Walter Trebels Dedicated to P. L. Butzer on the occasion of his 65-th birthday (To appear in Contemporary Mathematics) Abstract. A Lemma
More informationGaussian Quadrature Involving Einstein and Fermi Functions With an Application to Summation of Series*
MATHEMATICS OF COMPUTATION VOLUME 44, NUMBER 169 JANUARY 1985, PAGES 177-190 Gaussian Quadrature Involving Einstein and Fermi Functions With an Application to Summation of Series* By Walter Gautschi and
More informationPOLYNOMIAL EQUATIONS OVER MATRICES. Robert Lee Wilson. Here are two well-known facts about polynomial equations over the complex numbers
POLYNOMIAL EQUATIONS OVER MATRICES Robert Lee Wilson Here are two well-known facts about polynomial equations over the complex numbers (C): (I) (Vieta Theorem) For any complex numbers x,..., x n (not necessarily
More informationChapter 10 INTEGRATION METHODS Introduction Unit Coordinate Integration
Chapter 0 INTEGRATION METHODS 0. Introduction The finite element analysis techniques are always based on an integral formulation. At the very minimum it will always be necessary to integrate at least an
More informationSpring 2019 Exam 2 3/27/19 Time Limit: / Problem Points Score. Total: 280
Math 307 Spring 2019 Exam 2 3/27/19 Time Limit: / Name (Print): Problem Points Score 1 15 2 20 3 35 4 30 5 10 6 20 7 20 8 20 9 20 10 20 11 10 12 10 13 10 14 10 15 10 16 10 17 10 Total: 280 Math 307 Exam
More informationOn the existence of higher order polynomial. lattices with large figure of merit ϱ α.
On the existence of higher order polynomial lattices with large figure of merit Josef Dick a, Peter Kritzer b, Friedrich Pillichshammer c, Wolfgang Ch. Schmid b, a School of Mathematics, University of
More informationMATH 690 NOTES. 1. Associated graded rings
MATH 690 NOTES 1. Associated graded rings Question 1.1. Let G = k[g 1 ] = n 0 G n be a one-dimensional standard graded ring over the field k = G 0. What are necessary and sufficient conditions in order
More informationOn Polynomial Cases of the Unichain Classification Problem for Markov Decision Processes
On Polynomial Cases of the Unichain Classification Problem for Markov Decision Processes Eugene A. Feinberg Department of Applied Mathematics and Statistics State University of New York at Stony Brook
More information