A Note on the Recursive Calculation of Incomplete Gamma Functions
|
|
- Eileen Walsh
- 5 years ago
- Views:
Transcription
1 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 x is positive, is unstable more so the larger x. Nevertheless, the recursion can be used in the range 0 n x practically without error growth, and in larger ranges 0 n N with a loss of accuracy that can be controlled by suitably limiting N. Categories and Subject Descriptors: G.1.0 [Numerical Analysis]: General stability (and instability); G.1.2 [Numerical Analysis]: Approximation General Terms: Algorithms, Reliability Additional Key Words and Phrases: Incomplete gamma functions, recursive calculation 1. INTRODUCTION Our concern is with the incomplete gamma function Pa, x a, x a, a, x 0 x e t t a1 dt a a, x (1.1) for nonnegative values of the parameter a, and positive x. More precisely, we are interested in the generation of P n Pa n, x, 0 a 1, for n 0,1,2,... by means of the well-known recurrence relation P n P n1 xan1 e x For a 0 we define P 0 lim a30 Pa, x 1. a n, n 1, 2,... ; P 0 Pa, x. (1.2) Author s address: Department of Computer Sciences, Purdue University, West Lafayette, IN 47907; wxg@cs.purdue.edu. Permission to make digital / hard copy of part or all of this work for personal or classroom use is granted without fee provided that the copies are not made or distributed for profit or commercial advantage, the copyright notice, the title of the publication, and its date appear, and notice is given that copying is by permission of the ACM, Inc. To copy otherwise, to republish, to post on servers, or to redistribute to lists, requires prior specific permission and / or a fee ACM /99/ $5.00 ACM Transactions on Mathematical Software, Vol. 25, No. 1, March 1999, Pages
2 102 W. Gautschi It is known (cf., e.g., Gautschi [1997, Sect. 2.4] or Gautschi [1998, Sect. 6.2]) that, for x 0, this recursion is technically unstable in the sense of rounding errors being amplified by factors that grow monotonically to infinity as n 3. Since, however, serious error growth is delayed when x is large, the recursion can still be used without significant loss of accuracy, or with a controlled loss of accuracy, if the range of n is suitably limited. This is studied in Section 2 by means of uniform asymptotics and is illustrated in Section 3 by numerical examples. 2. ANALYSIS OF ERROR PROPAGATION The error propagation in the recursion (1.2) can be described by amplification factors n that determine the amplification s3t at n t of a small error introduced at n s, assuming no additional errors. Indeed (cf. Gautschi [1997, Eq. (2.11)] or Gautschi [1998, Eq. (6.2)]), s3t t s, (2.1) where, if t s, we are dealing with forward recursion and, if t s, with backward recursion. In the case of (1.2) we have 1 Pa, x n n a, x, n 0,1,2,..., (2.2) Pa n, x which, for 0 a 1, x 0, increases monotonically from 0 1 to. This means that forward recurrence is unstable, and backward recurrence very stable, so much so that Pa n, x for any fixed n can be obtained to any prescribed accuracy by recurring backward from some sufficiently large down to n, using an arbitrary starting value, say zero. On the other hand, when x is large, the growth of n is initially very slow, so that we may get away with recurring forward up to some n 0 n 0 a, x without suffering any significant loss of accuracy. We will show that n 0 x for any a with 0 a The Case a 0 If a 0, then n 0, x 1 Pn, x n, x. (2.3) 1 In Gautschi [1997; 1998] different normalizations are used for the incomplete gamma function which, however, do not affect the amplification factors, since they are invariant with respect to scaling.
3 A Note on the Recursive Calculation of Incomplete Gamma Functions 103 We are interested in estimating n such that n c 10 d, where 1 c 10. This means that, recurring from 0 to n, we must expect a loss of about d or d 1 decimal digits of accuracy. To obtain a realistic estimate of n, we make use of Temme s uniform asymptotic approximation [Temme 1975] in terms of the complementary error function erfc [Abramowitz et al. 1965, Chap. 7], a, x 1 erfc a/2, (2.4) a 2 where 21ln, x/a, and the plus or minus sign holds depending on whether 1 or 1. (The validity of (2.4) is uniform in x/a 0 as a 3 and/or x 3.) The equation n c 10 d then becomes Fn 0 (2.5) with Fn erfc x n nlnx/n 2 c 10d if n x, erfc x n nlnx/n 2 c 10d if n x. (2.6) We solve this equation numerically by the method of bisections, for x and d 0, c 2; d 0, c 5; and d 1110, c 1, using rational approximations of Cody [1969] to evaluate the complementary error function. The results are shown graphically in the upper left frame of Figure 1, where the higher lines correspond to larger values of d respectively c. The plot suggests that n 1 is relatively small for n x. We show indeed that n 0, x 2 whenever n x. (2.7) Since n, x is decreasing with increasing x, it suffices to prove (2.7) for x n. In this case, by (2.3), n 0, n n, n n, n 1, n, n 1 and the known inequalities n, n 1, n 1, 2, 3,..., (2.8) 2 attributed by Vietoris [1983] to G. Lochs, establish (2.7).
4 104 W. Gautschi Fig. 1. Lines of constant amplification factors for a 0,0.25,0.5, The Case 0 a 1 Rewriting (2.2) as a, x a n n n a, x a a n, x, (2.9) we can use (2.4) again to estimate the value of n for which n c 10 d. The relevant equation now is (2.5) with erfc x a n a nlnx/a n erfc x a alnx/a 1 c 10d if a n x, Fn (2.10) erfc x a n a nlnx/a n erfc x a alnx/a 1 c 10d if a n x. The solution, for a , is shown in the last three frames of Figure 1 for the same values of x, c, and d as before in Section 2.1. The
5 A Note on the Recursive Calculation of Incomplete Gamma Functions 105 curves are all quite similar and in particular suggest that forward recursion is safe if n does not exceed x. We can show this analytically if we use the fact that a n/a n, x in (2.9) is increasing as a function of a [Tricomi 1950] and a, x as a function of x. Then, for a 1, whenever n 1 x, we have n a, x n 1 n 1, x n 1 n 1, n 1 n 1 n 1 1, n 1 1 2, n 1, n 1 1 n 1 where (2.8) has been used again. Thus, similarly as in (2.7), n a, x 2 if 0 a 1, n x 1. (2.11) We can say, therefore, that for large x forward recursion in (1.2) can safely be used when n x and backward recursion, as indicated earlier, if n x. The graphs in Figure 1 moreover show the extent accuracy deteriorates when forward recursion is used up to an n that is significantly larger than x. 3. NUMERICAL EXAMPLES We conclude with some numerical examples that illustrate the results obtained in Section 2. For a and x , weletn 0 x and first generate P b n for n 0, 1,..., n 0, and for n n k 1 k/10n 0, k 1, 2,..., 5, by backward recursion, P b b 0, P n1 P b n xan1 e x, n, 1,..., 1, (3.1) a n b choosing large enough so that P nk (and hence also all P b n, n n k ) approximate the true values accurately to 12 significant decimal digits. We then generate P f n for the same values of n by forward recursion (1.2) and compute the relative errors e n P n f P n b /P n b (3.2) for n 0,1,..., n 0 and for n n k, k 1,2,..., 5. The results are similar for all values of a and are shown in Table I for a 0. The starting value Pa, x required in (1.2) is computed by a double-precision version of the routine gam in Gautschi [1979]. The fractional term on the right of (1.2)
6 106 W. Gautschi Table I. Errors in Forward Recursion Beyond n x, x 20(20)100 x n k P nk e nk max 0nn0 e n D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D D Table II. Errors in Forward Recursion Beyond n x, x 500, 1000 x n k P nk e nk max 0nn0 e n D D D D D D D D D D D D D D D D and (3.1), when x and n are large, is prone to overflow if computed as written. It is better to first take logarithms and then exponentiate, i.e., to compute this term as x an1 e x expa n 1lnx x lna n. (3.3) a n
7 A Note on the Recursive Calculation of Incomplete Gamma Functions 107 One can see from Table I that the maximum relative error in the range 0 n n 0 is consistent with the one requested, whereas the errors for n n k, k 1,2,....5, eventually become significantly larger, as predicted by the graphs in Figure 1. For still larger values of x, the effect is more pronounced, as is shown in Table II for x 500 and x It is interesting to observe that the paradigm of forward recursion for n x and backward recursion for n x is valid, though for different reasons, also when x is negative and Tricomi s definition of the incomplete gamma function is used (cf. Gautschi [1998, Sect. 6.2]). REFERENCES ABRAMOWITZ, M. AND STEGUN, I. A., Eds Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables. Dover Publications, Inc., Mineola, NY. CODY, W. J Rational Chebyshev approximations for the error function. Math. Comput. 23, GAUTSCHI, W Algorithm 542: Incomplete gamma functions. ACM Trans. Math. Softw. 5, 4 (Dec.), GAUTSCHI, W The computation of special functions by linear difference equations. In Advances in Difference Equations, S. Elaydi, I. Győri, and G. Ladas, Eds. Gordon and Breach Science Publishers, Inc., Langhorn, PA, GAUTSCHI, W The incomplete gamma functions since Tricomi. In Tricomi s Ideas and Contemporary Applied Mathematics. Atti dei Convegni Lincei, vol Accademia Nazionale dei Lincei, Rome, Italy, TEMME, N. M Uniform asymptotic expansions of the incomplete gamma functions and the incomplete beta function. Math. Comput. 29, TRICOMI, F. G Sulla funzione gamma incompleta. Ann. Mat. Pura Appl. 4, 31, VIETORIS, L Dritter Beweis der die unvollständige Gammafunktion betreffenden Lochsschen Ungleichungen. Österreich. Akad. Wiss. Math. Natur. Kl. Sitzungsber. 2, 192, Received: September 1998; revised: October 1998; accepted: November 1998
Defining Incomplete Gamma Type Function with Negative Arguments and Polygamma functions ψ (n) ( m)
1 Defining Incomplete Gamma Type Function with Negative Arguments and Polygamma functions ψ (n) (m) arxiv:147.349v1 [math.ca] 27 Jun 214 Emin Özc. ağ and İnci Ege Abstract. In this paper the incomplete
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 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 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 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 information1. Introduction The incomplete beta function I(a, b, x) is defined by [1, p.269, Eq and p.944, Eq ]
MATHEMATICS OF COMPUTATION Volume 65, Number 215 July 1996, Pages 1283 1288 AN ASYMPTOTIC EXPANSION FOR THE INCOMPLETE BETA FUNCTION B.G.S. DOMAN Abstract. A new asymptotic expansion is derived for the
More informationy(a, x) = f f "V' dt (a > 0) (1) ^ (-1)"
mathematics of computation, volume 26, number 119, july, 1972 On the Zeros of the Incomplete Gamma Function By K. S. Kölbig Abstract. Some asymptotic formulae given elsewhere for the zeros of the incomplete
More informationSUB- AND SUPERADDITIVE PROPERTIES OF EULER S GAMMA FUNCTION
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 135, Number 11, November 27, Pages 3641 3648 S 2-9939(7)957- Article electronically published on August 6, 27 SUB- AND SUPERADDITIVE PROPERTIES OF
More informationStability Analysis of Predator- Prey Models via the Liapunov Method
Stability Analysis of Predator- Prey Models via the Liapunov Method Gatto, M. and Rinaldi, S. IIASA Research Memorandum October 1975 Gatto, M. and Rinaldi, S. (1975) Stability Analysis of Predator-Prey
More informationMaximum Achievable Gain of a Two Port Network
Maximum Achievable Gain of a Two Port Network Ali Niknejad Siva Thyagarajan Electrical Engineering and Computer Sciences University of California at Berkeley Technical Report No. UCB/EECS-2016-15 http://www.eecs.berkeley.edu/pubs/techrpts/2016/eecs-2016-15.html
More informationOn the stirling expansion into negative powers of a triangular number
MATHEMATICAL COMMUNICATIONS 359 Math. Commun., Vol. 5, No. 2, pp. 359-364 200) On the stirling expansion into negative powers of a triangular number Cristinel Mortici, Department of Mathematics, Valahia
More informationCalculation of Cylindrical Functions using Correction of Asymptotic Expansions
Universal Journal of Applied Mathematics & Computation (4), 5-46 www.papersciences.com Calculation of Cylindrical Functions using Correction of Asymptotic Epansions G.B. Muravskii Faculty of Civil and
More informationPUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes
PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES Notes. x n+ = ax n has the general solution x n = x a n. 2. x n+ = x n + b has the general solution x n = x + (n )b. 3. x n+ = ax n + b (with a ) can be
More informationUniversity of Warwick institutional repository:
University of Warwick institutional repository: http://go.warwick.ac.uk/wrap This paper is made available online in accordance with publisher policies. Please scroll down to view the document itself. Please
More informationTitle: Closed-form approximations to the Error and Complementary Error Functions and their applications in atmospheric science
Title: Closed-form approximations to the Error and Complementary Error Functions and their applications in atmospheric science Short title: Closed-form of Error and Complementary Error Functions Authors:
More informationPeriodicity and Solution of Rational Recurrence Relation of Order Six
Applied Mathematics 3 79-733 http://d.doi.org/.436/am..377 Published Online July (http://www.scirp.org/journal/am) Periodicity and Solution of Rational Recurrence Relation of Order Si Tarek F. Ibrahim
More information7.4* General logarithmic and exponential functions
7.4* General logarithmic and exponential functions Mark Woodard Furman U Fall 2010 Mark Woodard (Furman U) 7.4* General logarithmic and exponential functions Fall 2010 1 / 9 Outline 1 General exponential
More informationINEQUALITIES FOR THE GAMMA FUNCTION
INEQUALITIES FOR THE GAMMA FUNCTION Received: 16 October, 26 Accepted: 9 February, 27 Communicated by: XIN LI AND CHAO-PING CHEN College of Mathematics and Informatics, Henan Polytechnic University, Jiaozuo
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 informationComputation of the error functions erf and erfc in arbitrary precision with correct rounding
Computation of the error functions erf and erfc in arbitrary precision with correct rounding Sylvain Chevillard Arenaire, LIP, ENS-Lyon, France Sylvain.Chevillard@ens-lyon.fr Nathalie Revol INRIA, Arenaire,
More informationBlame and Coercion: Together Again for the First Time
Blame and Coercion: Together Again for the First Time Supplementary Material Jeremy Siek Indiana University, USA jsiek@indiana.edu Peter Thiemann Universität Freiburg, Germany thiemann@informatik.uni-freiburg.de
More informationA Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions
mathematics of computation volume 39, number 159 july 1982, pages 201-206 A Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions By John P. Boyd Abstract. The theorem proved here extends
More informationInteger Least Squares: Sphere Decoding and the LLL Algorithm
Integer Least Squares: Sphere Decoding and the LLL Algorithm Sanzheng Qiao Department of Computing and Software McMaster University 28 Main St. West Hamilton Ontario L8S 4L7 Canada. ABSTRACT This paper
More information6.3 Exponential Integrals. Chi-Square Probability Function
6.3 Exponential Integrals 5 Chi-Square Probability Function P (χ ν is defined as the probability that the observed chi-square for a correct model should be less than a value χ. (We will discuss the use
More informationClosed form expressions for two harmonic continued fractions
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part B Faculty of Engineering and Information Sciences 07 Closed form expressions for two harmonic continued
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 informationSharp inequalities and complete monotonicity for the Wallis ratio
Sharp inequalities and complete monotonicity for the Wallis ratio Cristinel Mortici Abstract The aim of this paper is to prove the complete monotonicity of a class of functions arising from Kazarinoff
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 informationFractional part integral representation for derivatives of a function related to lnγ(x+1)
arxiv:.4257v2 [math-ph] 23 Aug 2 Fractional part integral representation for derivatives of a function related to lnγ(x+) For x > let Mark W. Coffey Department of Physics Colorado School of Mines Golden,
More informationOptimization of Rational Approximations by Continued Fractions
Bergische Universität Wuppertal Optimization of Rational Approximations by Continued Fractions Frithjof Blomquist Preprint BUW-WRSWT 2006/3 Wissenschaftliches Rechnen/ Softwaretechnologie Impressum Herausgeber:
More informationCopyright 1972, by the author(s). All rights reserved.
Copyright 1972, by the author(s). All rights reserved. Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are
More informationOscillation Criteria for Delay and Advanced Difference Equations with General Arguments
Advances in Dynamical Systems Applications ISSN 0973-531, Volume 8, Number, pp. 349 364 (013) http://campus.mst.edu/adsa Oscillation Criteria for Delay Advanced Difference Equations with General Arguments
More informationIntroductions to ExpIntegralEi
Introductions to ExpIntegralEi Introduction to the exponential integrals General The exponential-type integrals have a long history. After the early developments of differential calculus, mathematicians
More informationConditions for convexity of a derivative and some applications to the Gamma function
Aequationes Math. 55 (1998) 73 80 0001-9054/98/03073-8 $ 1.50+0.0/0 c Birkhäuser Verlag, Basel, 1998 Aequationes Mathematicae Conditions for conveity of a derivative and some applications to the Gamma
More informationSkill 6 Exponential and Logarithmic Functions
Skill 6 Exponential and Logarithmic Functions Skill 6a: Graphs of Exponential Functions Skill 6b: Solving Exponential Equations (not requiring logarithms) Skill 6c: Definition of Logarithms Skill 6d: Graphs
More informationThe Generating Functions for Pochhammer
The Generating Functions for Pochhammer Symbol { }, n N Aleksandar Petoević University of Novi Sad Teacher Training Faculty, Department of Mathematics Podgorička 4, 25000 Sombor SERBIA and MONTENEGRO Email
More informationRecursive Computation of the Repeated Integrals of the Error Function
Recursive Computation of the Repeated Integrals of the Error Function By Walter Gautschi 1. This paper is concerned with a special technique, originated by J. C. P. Miller [1, p. xvii], of computing a
More informationarxiv:math-ph/ v1 10 Jan 2005
Asymptotic and eact series representations for the incomplete Gamma function arxiv:math-ph/5119v1 1 Jan 5 Paolo Amore Facultad de Ciencias, Universidad de Colima, Bernal Díaz del Castillo 34, Colima, Colima,
More informationHorst Alzer A MEAN VALUE INEQUALITY FOR THE DIGAMMA FUNCTION
Rendiconti Sem. Mat. Univ. Pol. Torino Vol. 75, 2 (207), 9 25 Horst Alzer A MEAN VALUE INEQUALITY FOR THE DIGAMMA FUNCTION Abstract. A recently published result states that for all ψ is greater than or
More informationCOMPLETE MONOTONICITIES OF FUNCTIONS INVOLVING THE GAMMA AND DIGAMMA FUNCTIONS. 1. Introduction
COMPLETE MONOTONICITIES OF FUNCTIONS INVOLVING THE GAMMA AND DIGAMMA FUNCTIONS FENG QI AND BAI-NI GUO Abstract. In the article, the completely monotonic results of the functions [Γ( + 1)] 1/, [Γ(+α+1)]1/(+α),
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES Infinite sequences and series were introduced briefly in A Preview of Calculus in connection with Zeno s paradoxes and the decimal representation
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES We have: Seen how to interpret derivatives as slopes and rates of change Seen how to estimate derivatives of functions given by tables of values Learned how
More informationPreliminary Results on Social Learning with Partial Observations
Preliminary Results on Social Learning with Partial Observations Ilan Lobel, Daron Acemoglu, Munther Dahleh and Asuman Ozdaglar ABSTRACT We study a model of social learning with partial observations from
More informationInequalities and Monotonicity For The Ratio of Γ p Functions
Int. J. Open Problems Compt. Math., Vol. 3, No., March 200 ISSN 998-6262; Copyright c ICSRS Publication, 200 www.i-csrs.org Inequalities and Monotonicity For The Ratio of Γ p Functions Valmir Krasniqi
More information6.4 Incomplete Beta Function, Student s Distribution, F-Distribution, Cumulative Binomial Distribution
6.4 Incomplete Beta Function, Student s Distribution, F-Distribution, Cumulative Binomial Distribution 29 CITED REFERENCES AND FURTHER READING: Stegun, I.A., and Zucker, R. 974, Journal of Research of
More informationDetermining Appropriate Precisions for Signals in Fixed-Point IIR Filters
38.3 Determining Appropriate Precisions for Signals in Fixed-Point IIR Filters Joan Carletta Akron, OH 4435-3904 + 330 97-5993 Robert Veillette Akron, OH 4435-3904 + 330 97-5403 Frederick Krach Akron,
More information2.2. Limits Involving Infinity. Copyright 2007 Pearson Education, Inc. Publishing as Pearson Prentice Hall
2.2 Limits Involving Infinity Copyright 2007 Pearson Education, Inc. Publishing as Pearson Prentice Hall Finite Limits as x ± What you ll learn about Sandwich Theorem Revisited Infinite Limits as x a End
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 informationThe Complexity of Forecast Testing
The Complexity of Forecast Testing LANCE FORTNOW and RAKESH V. VOHRA Northwestern University Consider a weather forecaster predicting the probability of rain for the next day. We consider tests that given
More informationGlobal Asymptotic Stability of a Nonlinear Recursive Sequence
International Mathematical Forum, 5, 200, no. 22, 083-089 Global Asymptotic Stability of a Nonlinear Recursive Sequence Mustafa Bayram Department of Mathematics, Faculty of Arts and Sciences Fatih University,
More informationThe numerical evaluation of a challenging integral
The numerical evaluation of a challenging integral Walter Gautschi Abstract Standard numerical analysis tools, combined with elementary calculus, are deployed to evaluate a densely and wildly oscillatory
More informationAlgorithm 808: ARfit A Matlab package for the estimation of parameters and eigenmodes of multivariate autoregressive models
Algorithm 808: ARfit A Matlab package for the estimation of parameters and eigenmodes of multivariate autoregressive models TAPIO SCHNEIDER New York University and ARNOLD NEUMAIER Universität Wien ARfit
More informationObjectives. Use the number e to write and graph exponential functions representing realworld
Objectives Use the number e to write and graph exponential functions representing realworld situations. Solve equations and problems involving e or natural logarithms. natural logarithm Vocabulary natural
More informationConvergence, Steady State, and Global Gains of Agent-Based Coalition Formation in E-markets
Convergence, Steady State, and Global Gains of Agent-Based Coalition Formation in E-markets KWANG MONG SIM The Chinese University of Hong Kong, Hong Kong YUANSHI WANG and HONG WU Zhongshan University,
More informationSEQUENTIAL ESTIMATION OF DYNAMIC DISCRETE GAMES: A COMMENT. MARTIN PESENDORFER London School of Economics, WC2A 2AE London, U.K.
http://www.econometricsociety.org/ Econometrica, Vol. 78, No. 2 (arch, 2010), 833 842 SEQUENTIAL ESTIATION OF DYNAIC DISCRETE GAES: A COENT ARTIN PESENDORFER London School of Economics, WC2A 2AE London,
More informationVerbatim copying and redistribution of this document. are permitted in any medium provided this notice and. the copyright notice are preserved.
New Results on Primes from an Old Proof of Euler s by Charles W. Neville CWN Research 55 Maplewood Ave. West Hartford, CT 06119, U.S.A. cwneville@cwnresearch.com September 25, 2002 Revision 1, April, 2003
More informationERROR ANALYSIS FOR THE COMPUTATION OF ZEROS OF REGULAR COULOMB WAVE FUNCTION AND ITS FIRST DERIVATIVE
MATHEMATICS OF COMPUTATION Volume 70, Number 35, Pages 1195 104 S 005-5718(00)0141- Article electronically published on March 4, 000 ERROR ANALYSIS FOR THE COMPUTATION OF ZEROS OF REGULAR COULOMB WAVE
More informationJUST THE MATHS UNIT NUMBER 1.4. ALGEBRA 4 (Logarithms) A.J.Hobson
JUST THE MATHS UNIT NUMBER 1.4 ALGEBRA 4 (Logarithms) by A.J.Hobson 1.4.1 Common logarithms 1.4.2 Logarithms in general 1.4.3 Useful Results 1.4.4 Properties of logarithms 1.4.5 Natural logarithms 1.4.6
More informationSimulating from a gamma distribution with small shape parameter
Simulating from a gamma distribution with small shape parameter arxiv:1302.1884v2 [stat.co] 16 Jun 2013 Ryan Martin Department of Mathematics, Statistics, and Computer Science University of Illinois at
More informationRADICAL AND RATIONAL FUNCTIONS REVIEW
RADICAL AND RATIONAL FUNCTIONS REVIEW Name: Block: Date: Total = % 2 202 Page of 4 Unit 2 . Sketch the graph of the following functions. State the domain and range. y = 2 x + 3 Domain: Range: 2. Identify
More informationSOME EQUATIONS WITH FEATURES OF DIGIT REVERSAL AND POWERS
SOME EQUATIONS WITH FEATURES OF DIGIT REVERSAL AND POWERS GEOFFREY B CAMPBELL AND ALEKSANDER ZUJEV Abstract. In this paper we consider integers in base 10 like abc, where a, b, c are digits of the integer,
More informationInvestigating Geometric and Exponential Polynomials with Euler-Seidel Matrices
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 14 (2011), Article 11.4.6 Investigating Geometric and Exponential Polynomials with Euler-Seidel Matrices Ayhan Dil and Veli Kurt Department of Mathematics
More informationInformation Storage Capacity of Crossbar Switching Networks
Information Storage Capacity of Crossbar Switching etworks ABSTRACT In this work we ask the fundamental uestion: How many bits of information can be stored in a crossbar switching network? The answer is
More informationON THE OSCILLATION OF THE SOLUTIONS TO LINEAR DIFFERENCE EQUATIONS WITH VARIABLE DELAY
Electronic Journal of Differential Equations, Vol. 008(008, No. 50, pp. 1 15. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp ON THE
More informationA proof of a partition conjecture of Bateman and Erdős
proof of a partition conjecture of Bateman and Erdős Jason P. Bell Department of Mathematics University of California, San Diego La Jolla C, 92093-0112. US jbell@math.ucsd.edu 1 Proposed Running Head:
More informationApproximations for the Bessel and Struve Functions
MATHEMATICS OF COMPUTATION VOLUME 43, NUMBER 168 OCTOBER 1984. PAGES 551-556 Approximations for the Bessel and Struve Functions By J. N. Newman Abstract. Polynomials and rational-fraction approximations
More informationMath 118B Solutions. Charles Martin. March 6, d i (x i, y i ) + d i (y i, z i ) = d(x, y) + d(y, z). i=1
Math 8B Solutions Charles Martin March 6, Homework Problems. Let (X i, d i ), i n, be finitely many metric spaces. Construct a metric on the product space X = X X n. Proof. Denote points in X as x = (x,
More informationGeometric interpretation of signals: background
Geometric interpretation of signals: background David G. Messerschmitt Electrical Engineering and Computer Sciences University of California at Berkeley Technical Report No. UCB/EECS-006-9 http://www.eecs.berkeley.edu/pubs/techrpts/006/eecs-006-9.html
More informationMathematical preliminaries and error analysis
Mathematical preliminaries and error analysis Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan September 12, 2015 Outline 1 Round-off errors and computer arithmetic
More informationExponentials of Symmetric Matrices through Tridiagonal Reductions
Exponentials of Symmetric Matrices through Tridiagonal Reductions Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A simple and efficient numerical algorithm
More informationExit Criteria for Simpson's Compound Rule*
MATHEMATICS OF COMPUTATION, VOLUME 26, NUMBER 119, JULY, 1972 Exit Criteria for Simpson's Compound Rule* By J. H. Rowland and Y. L. Varol Abstract. In many automated numerical algorithms, the calculations
More informationSkill 6 Exponential and Logarithmic Functions
Skill 6 Exponential and Logarithmic Functions Skill 6a: Graphs of Exponential Functions Skill 6b: Solving Exponential Equations (not requiring logarithms) Skill 6c: Definition of Logarithms Skill 6d: Graphs
More informationS = {x 1 2 _< x n, x prime} for any given n 2. (2.1)
An Exercise in Program Explanation JAYADEV MISRA The University of Texas at Austin A combination of program-proving ideas and stepwise refinement is used to develop and explain an algorithm which uses
More informationWKB solution of the wave equation for a plane angular sector
PHYSICAL REVIEW E VOLUME 58, NUMBER JULY 998 WKB solution of the wave equation for a plane angular sector Ahmad T. Abawi* Roger F. Dashen Physics Department, University of California, San Diego, La Jolla,
More informationLatter research on Euler-Mascheroni constant. 313, Bucharest, Romania, Târgovişte, Romania,
Latter research on Euler-Mascheroni constant Valentin Gabriel Cristea and Cristinel Mortici arxiv:3.4397v [math.ca] 6 Dec 03 Ph. D. Student, University Politehnica of Bucharest, Splaiul Independenţei 33,
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 informationSOME INEQUALITIES FOR THE q-digamma FUNCTION
Volume 10 (009), Issue 1, Article 1, 8 pp SOME INEQUALITIES FOR THE -DIGAMMA FUNCTION TOUFIK MANSOUR AND ARMEND SH SHABANI DEPARTMENT OF MATHEMATICS UNIVERSITY OF HAIFA 31905 HAIFA, ISRAEL toufik@mathhaifaacil
More informationRevision notes for Pure 1(9709/12)
Revision notes for Pure 1(9709/12) By WaqasSuleman A-Level Teacher Beaconhouse School System Contents 1. Sequence and Series 2. Functions & Quadratics 3. Binomial theorem 4. Coordinate Geometry 5. Trigonometry
More information2.1 Limits, Rates of Change and Slopes of Tangent Lines
2.1 Limits, Rates of Change and Slopes of Tangent Lines (1) Average rate of change of y f x over an interval x 0,x 1 : f x 1 f x 0 x 1 x 0 Instantaneous rate of change of f x at x x 0 : f x lim 1 f x 0
More informationPartition of Integers into Distinct Summands with Upper Bounds. Partition of Integers into Even Summands. An Example
Partition of Integers into Even Summands We ask for the number of partitions of m Z + into positive even integers The desired number is the coefficient of x m in + x + x 4 + ) + x 4 + x 8 + ) + x 6 + x
More informationA VERTICAL LOOK AT KEY CONCEPTS AND PROCEDURES ALGEBRA I. Texas Education Agency
A VERTICAL LOOK AT KEY CONCEPTS AND PROCEDURES ALGEBRA I Texas Education Agency The materials are copyrighted (c) and trademarked (tm) as the property of the Texas Education Agency (TEA) and may not be
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 informationAsymptotic Notation. such that t(n) cf(n) for all n n 0. for some positive real constant c and integer threshold n 0
Asymptotic Notation Asymptotic notation deals with the behaviour of a function in the limit, that is, for sufficiently large values of its parameter. Often, when analysing the run time of an algorithm,
More informationS. R. Tate. Stable Computation of the Complex Roots of Unity, IEEE Transactions on Signal Processing, Vol. 43, No. 7, 1995, pp
Stable Computation of the Complex Roots of Unity By: Stephen R. Tate S. R. Tate. Stable Computation of the Complex Roots of Unity, IEEE Transactions on Signal Processing, Vol. 43, No. 7, 1995, pp. 1709
More informationHighly accurate evaluation of the few-body auxiliary functions and four-body integrals
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSB: ATOMIC, MOLECULAR AND OPTICAL PHYSICS J. Phys. B: At. Mol. Opt. Phys. 36 (2003) 1857 1867 PII: S0953-4075(03)58918-0 Highly accurate evaluation of the
More informationTotal Ordering on Subgroups and Cosets
Total Ordering on Subgroups and Cosets Alexander Hulpke Department of Mathematics Colorado State University 1874 Campus Delivery Fort Collins, CO 80523-1874 hulpke@math.colostate.edu Steve Linton Centre
More information5.6 Logarithmic and Exponential Equations
SECTION 5.6 Logarithmic and Exponential Equations 305 5.6 Logarithmic and Exponential Equations PREPARING FOR THIS SECTION Before getting started, review the following: Solving Equations Using a Graphing
More informationThe School District of Palm Beach County Algebra 1 Honors Unit A: Data Analysis
Unit A: Data Analysis MAFS.912.S ID.1.1 MAFS.912.S ID.1.2 MAFS.912.S ID.1.3 MAFS.912.S ID.2.5 Calculator: Yes Mathematics Florida Represent data with plots on the real number line (dot plots, histograms,
More informationLIMITS AND DERIVATIVES
2 LIMITS AND DERIVATIVES LIMITS AND DERIVATIVES 2.2 The Limit of a Function In this section, we will learn: About limits in general and about numerical and graphical methods for computing them. THE LIMIT
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
On Runs of Residues Author(s): D. H. Lehmer and Emma Lehmer Source: Proceedings of the American Mathematical Society, Vol. 13, No. 1 (Feb., 1962), pp. 102-106 Published by: American Mathematical Society
More informationCURRICULUM MAP. Course/Subject: Honors Math I Grade: 10 Teacher: Davis. Month: September (19 instructional days)
Month: September (19 instructional days) Numbers, Number Systems and Number Relationships Standard 2.1.11.A: Use operations (e.g., opposite, reciprocal, absolute value, raising to a power, finding roots,
More informationLogarithmic Functions
Logarithmic Functions Definition 1. For x > 0, a > 0, and a 1, y = log a x if and only if x = a y. The function f(x) = log a x is called the logarithmic function with base a. Example 1. Evaluate the following
More information8.5 Taylor Polynomials and Taylor Series
8.5. TAYLOR POLYNOMIALS AND TAYLOR SERIES 50 8.5 Taylor Polynomials and Taylor Series Motivating Questions In this section, we strive to understand the ideas generated by the following important questions:
More informationAn Algorithm for Numerical Reference Generation in Symbolic Analysis of Large Analog Circuits
An Algorithm for Numerical Reference Generation in Symbolic Analysis of Large Analog Circuits Ignacio García-Vargas, Mariano Galán, Francisco V. Fernández and Angel Rodríguez-Vázquez IMSE, Centro Nacional
More informationarxiv: v2 [math.gm] 24 May 2016
arxiv:508.00533v [math.gm] 4 May 06 A PROOF OF THE RIEMANN HYPOTHESIS USING THE REMAINDER TERM OF THE DIRICHLET ETA FUNCTION. JEONWON KIM Abstract. The Dirichlet eta function can be divided into n-th partial
More informationMATHEMATICS OF COMPUTATION, VOLUME 26, NUMBER 118, APRIL An Algorithm for Computing Logarithms. and Arctangents. By B. C.
MATHEMATICS OF COMPUTATION, VOLUME 26, NUMBER 118, APRIL 1972 An Algorithm for Computing Logarithms and Arctangents By B. C. Carlson Abstract. An iterative algorithm with fast convergence can be used to
More informationOptimality of an Algorithm Solving the Bottleneck Tower of Hanoi Problem
25 Optimality of an Algorithm Solving the Bottleneck Tower of Hanoi Problem YEFIM DINITZ AND SHAY SOLOMON Ben-Gurion University of the Negev Abstract. We study the Bottleneck Tower of Hanoi puzzle posed
More information6.3 Exponential Integrals. Chi-Square Probability Function
6.3 Exponential Integrals 5 Chi-Square Probability Function P (χ ν is defined as the probability that the observed chi-square for a correct model should be less than a value χ. (We will discuss the use
More informationDifference Equations for Multiple Charlier and Meixner Polynomials 1
Difference Equations for Multiple Charlier and Meixner Polynomials 1 WALTER VAN ASSCHE Department of Mathematics Katholieke Universiteit Leuven B-3001 Leuven, Belgium E-mail: walter@wis.kuleuven.ac.be
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