A CLASS OF SLOWLY CONVERGENT SEMES AND THEIR SUMMATION BY GAUSSIAN QUADRATURE WALTER GAUTSCHI
|
|
- Amber Lee
- 5 years ago
- Views:
Transcription
1 mathematics of computation volume 57, number 195 july 1991, pages A CLASS OF SLOWLY CONVERGENT SEMES AND THEIR SUMMATION BY GAUSSIAN QUADRATURE WALTER GAUTSCHI Abstract. Series are considered whose general term is a rational function multiplied by a fractional power. The summation of such series is reduced, via Laplace transformation techniques, to a problem of quadrature, which is then solved by Gaussian quadrature relative to Einstein and Fermi weight functions. A number of examples are worked out in detail. We consider series of the type 1. Introduction 1.0) oo v-l S0 = E^ r(k) k=l orr oo (1.1) Sx=^2(-X)k-lk'/-lr(k), k=l where 0 < v < 1 and r(-) is a rational function.2, ko-eg with p, q real polynomials of degrees degp < deg q. Strict inequality is assumed when necessary for convergence. It is further assumed that the zeros of q all have nonpositive real parts: (1.3) ifq(-a) = 0 then Rea > 0. This condition can always be achieved by a preliminary summation of a few initial terms. The problem can be simplified by first obtaining the partial fraction decomposition of r, (1-4) r(s) = 2 2cpJs+^ym + EEKJs+avym+e?m(s+*yrm], p m=l y m=l Received March 7, 1990; revised September 17, Mathematics Subject Classification (1985 Revision). Primary 40A25; Secondary 44A10, 65D30, 33A65. Key words and phrases. Slowly convergent series, Laplace transformation, summation by quadrature, Gaussian quadrature, orthogonal polynomials. Work supported, in part, by the National Science Foundation under grant CCR American Mathematical Society /91 $1.00+ $.25 per page
2 310 WALTER GAUTSCHI where the first sum is over all real zeros (-a ) of q (with multiplicities m ), and the second sum is over all pairs of conjugate complex zeros (-a, -a) (with multiplicities m ). The coefficients c in the first sum are real, those in the second complex, in general. (We have assumed in (1.4) that degp < degö. If degp = degö, and we are thus dealing with Sx in (1.1), there will be an additional constant term in (1.4). Its contribution to the series in (1.1) is expressible in terms of the Riemann zeta function; cf. (3.9) below for m = 0.) Once the decomposition (1.4) has been obtained (for relevant constructive methods, see, e.g., [8, 7.1]), it clearly suffices to consider (1.5) r(s)= l Rea>0, m>x. (s + a) Without restriction of generality, it may be further assumed that Im a > 0. By interpreting the terms in the series ( 1.0) and ( 1.1 ) as Laplace transforms at integer values, it is possible to express the sum of the series as a weighted integral over R+ of certain special functions related to the incomplete gamma function. The weighting involves the product of a fractional power and either Einstein's function t(e' - 1)_1 (in the case of (1.0)), or Fermi's function (e' + 1)_1 (in the case of ( 1.1 )), both having infinitely many poles on the imaginary axis of the complex plane. Properties of the required special functions are briefly developed in 2. Section 3 discusses the summation of ( 1.0), ( 1.1 ) via Gaussian quadrature. The case v = X of purely rational series is treated in 4 and complements more traditional techniques (e.g., those in [8, 7.211]). Numerical examples for the case v = \ are given in 5, where also comparisons are made with direct summation and accelerated summation using the e-algorithm. Define, for / > 0, 2. Preliminaries r1 (2.1_,) g_ t\v) = g_x(i)= _, 0<v<X, and for t > 0, «= 0,1,2,..., e-at^-l,/ (2.1.) g»it'a>u) = g»{t) = W(r=T)j0eTit-t)x~'/dt' Lemma 2.1. We have Rea>0, lma>0, a^0, 0<v<X. (2.2) g0(t;a,i>) = e~aly*(x -v, -at), where y* is Tricomi's form of the incomplete gamma function (cf. [3, eq ]). Furthermore, (2.3) ^«= ^t{(/ + ^±^)^)-^-1(/)}' n = 0, 1,2,...
3 Proof. By definition, slowly convergent series and their summation 311 which, upon the change of variables u = -ax, gives.. e-at(-atrl rat -u -», g (?) = r(i- ) /o e u du 0-at( nt\-t\-v) e \~at) /i <\ -at */,... = ni-y) y(i -". -«O = «v(i-v,-at) (cf. [3, eqs , 6.5.4]). Next,?n+lW e-atf-x ( +1)!F( t _ e-'f-1 [' at - n+x8"[) (n + xy.t(x-u)j0 e [t,n 1-y j T) x dx. To the last integral we apply integration by parts, letting u = xl "(t - x)n, v = eax, hence u=(x-v)x (t-x) - nx (t-x) = (n+x- i/)(/ - x)"x~u - nt(t - x)"~xx~v, 1 a, v = -e. a This yields (2.3) for n > X. A similar calculation gives 8^={t + ^)go^-añhf)' which, in view of (2.1_,), shows that (2.3) holds also for n = 0. To the author's knowledge, no software seems to be readily available for computing the incomplete gamma function in the domain of interest here (left half plane), and one thus has to rely on standard techniques such as power series and asymptotic expansions. An effort of developing good software for Tricomi's function y*(x - v, z), 0 < v < X, applicable for arbitrary complex z, would certainly be worthwhile. For the case v = 5, however, see Summation of 50 and Sx We employ a technique already used in [7], namely, to interpret the general term of the series as a Laplace transform and thereby converting the series into a suitably weighted integral. We first treat the series S0 in (1.0). Assume r(-) given as in (1.5), and consider first the case a^o. Then 1 f t~v jn-l -at \ (3-1} /_1-(7T^=^{rTr^*(^T)!}' 0<1/<U D
4 312 WALTER GAUTSCHI where i?{-} is the Laplace transform and * means convolution. We thus have = íl ^m_,(i;a, v), rn>x, where gn(t; a,v) Thus, is defined in (2.1n). There follows oo Ku 1 roo Et ^ = EW)w = E/ ^'/(o^ fe=1^ + aj fc=1 fc=1->o = EÍ ^'-V(WV,(í;m)^ OO,V-l»oo L-r 1-1/V^ -(fc-l)i /. s j^ * r E* Sm-l{t\a,v)dt 0 fe=i / ' v-r ;Sm-l(f,a,v)dt. Jo e i (3J) H-nrr^=l ^l/^t)gm_x(t;a,v)dt, m>x,0<u<x, ~ (Ac + a) Jo where «( ) is "Einstein's function" (cf. [7]), (3.4) e(t) = -J, 0<t<oo. e - 1 Since gm_x(-; a, v), by Lemma 2.1, is an entire function, formula (3.3) suggests to apply Gaussian quadrature to the integral on the right, using w(t) = t~ve(t) as a weight function on [0, oo]. The required orthogonal polynomials can be computed by techniques already discussed in [7] (see, in particular, (2.3) and Example 4.4 of that paper). Gauss quadrature will converge quite rapidly, unless Re a and/or Im a is large, in which case "stratified" summation can be employed to regain rapid convergence (see Examples 5.1 and 5.4). For v = \, the first 80 recursion coefficients for the required orthogonal polynomials are given to 25 significant digits in Table 1 of the Appendix. It is easily seen that (3.3) holds also for a = 0 if we define (3.5) 8 {t) = T(2 - u) ' gn+i^= n + 2-vg^t)' n = 0, 1,2,... (fl = 0). The series can then be expressed in terms of Riemann's zeta function, oo (3.6) Efe"(m+1"") = i(m+1-1')' k=\
5 SLOWLY CONVERGENT SERIES AND THEIR SUMMATION 313 and since gm_x is a monomial of degree m - 1, the «-point Gauss formula for the integral in (3.3) gives exact answers (modulo rounding) if «= [(m-\- 1)/2J. For the sum Sx in (1.1), a calculation similar to the one which led to (3.3) now yields, for a ^ 0, oo i'"-1 r (3.7) D-'f'iFnr-JL <">«> * _,«,<.,*>*, where <p(-) is the "Fermi weight function" (cf. [7]), (3.8) (p(t) = -,-, 0<i<oo. e + X m > 0, 0 < v < X, The result (3.7), as noted, holds also for m = 0, if g_x is defined as in (2.1_ ). If a = 0, then (3.7) holds with gm_x, m > 0, defined in (2.1_ ) and (3.5), and represents the series oo (3.9) (-l)*-1fc-(m+1-') = (1-2-m+")Ç(m + 1-1/). k=i Again, Gauss quadrature for the integral in (3.7), in this case, is exact if we take «= [(m + 2)/2J points. The first 80 recursion coefficients for the orthogonal polynomials with respect to the weight function t~"'tp(t), v = \, are listed in Table 2 of the Appendix. Remarks. X. The fractional power kv~ ' in ( 1.0), ( 1.1 ) can easily be generalized to (k + b)v~x, Reb > 0, since this only introduces a factor e~ ' in (3.1) and gives f(t) = tl~ue~ 'gm_x(t; a-b,v) in place of (3.2), hence (3-10) E(^+JS = rt-ve(t)e-btgm_x(t;a-b,v)dt, l~rx (k + a) Jo and the similarly generalized version of (3.7), oo / j u\tj 1 /*oo (3.11) D"1* f^i» =/ t-v<p(t).te-b'gm_x(t;a-b,v)dt. k=l ^ ' ^ 2. As an alternative to Gauss quadrature with weight function Cl'e(t), one could write the integral in (3.3) as / OO /»OO. (3.12) / Jo re(t)gm_x(t;a,u)dt= Jo fv- gm_x(f, a, v)dt l - e and apply Gauss-Laguerre quadrature to the integral on the right. The poles that were previously incorporated in the weight function t~"e(t) then, however, become part of the integrand, which may adversely affect the speed of convergence of the quadrature scheme. This was indeed found to be the case when a is relatively small, but for larger values of a, in particular when used
6 314 WALTER GAUTSCHI in conjunction with stratified summation (cf. Examples 5.1 and 5.4), the Gauss- Laguerre method is competitive and, in Example 5.1, even more efficient. The same alternative approach is possible for the integral in (3.7), if written as / OO /»OO. (3.13)/ fvtp(t)-tgm_x(t;a,v)dt=\ C" e~' --tgm_x(t; a,v)dt, Jo Jo 1 +e although this time the poles affect convergence more severely, since they are half as close to the real axis than before. Still, for sufficiently large values of a, Gauss-Laguerre quadrature here, too, gains the upper hand. Finally, if many different values of v were contemplated, then Gauss- Laguerre would also be preferable, since the recursion coefficients for the respective orthogonal polynomials the generalized Laguerre polynomials are then explicitly known and need not be tabulated. In the case of the weight function w(t) = t~ve~l on (0, oo), the corresponding (monic) orthogonal polynomials 71 k(') = nk('> w} mdeed satisfy f314) nk+l(t) = (t-ak)nk(t)-ßknk_x(t), k = 0, X, 2,..., n_x(t) = 0, n0(t) = X, with the coefficients ak = ak(w), ßk = ßk(w) [ß0(w) = f0 w(t)dt] having the particularly simple form ak(w) = 2k+X-u, k>0; (3.15) k (w = tvet). ßQ(w) = Y(X-u), ßk(w) = k(k-u), k>x 4. Series of purely rational terms So far, we assumed that 0 < v < 1 in (1.0) and (1.1). The same techniques, however, are applicable when v = X, i.e., for series oo (4.1) S; = E^)> S\=^2(-X)k-lr(k). One finds, when k=l oo fc=l (4.2) r(s)= ( * m, Rea>0, for m > 2 that and for m > 1 that with e and cp as given in (3.4) and (3.8), respectively.
7 SLOWLY CONVERGENT SERIES AND THEIR SUMMATION 315 For more general rational functions (4-5) ^) = EEu* + «,) i m=l we get from (4.3) oo ^.OO (4.6) Er(fc)=/ Ws(t)dt, where k=i Jo (4-7) i(o = EEAi', (m-l)! i m=l v ' cim jn-2 -a,t Evidently, g is again an entire function, provided (4.8) Ecn=0, i which is required for the series in (4.6) to converge. Similarly, oo,_, roc (4.9) Y,(-V r(k)= / <p(t)h(t)dt, rr; Jo where k=l (4-10) *<o«ee( rfy i m=l cim jn-\ -a,t an entire function for any choice of the coefficients c/w. The first 40 recursion coefficients for the polynomials orthogonal with respect to e and q> can be found to 25 significant digits in [7, Appendices Al and A2]. The method described in this section provides an alternative to other summation/integration methods, such as those discussed in [8, 7.211]. An advantage of the present method is that it leads to Gaussian quadrature of entire functions, a possible complication, that the interval of integration is infinite. 5. Examples In all of our examples we take v = \. In this case, the function g0 in (2.2) is given by (cf. [3, eq ])... / 1\ -at *(\ \ 2 F(Val) (5.1) 8o{f,a,-2)=e y^,-at)=^-^, where F is Dawson's integral, (5.2) F(z) = e~zl [' / dt. Jo For real z, this can be evaluated with an accuracy of up to about 20 significant digits, using the rational Chebyshev approximations given in [2]. All computations reported below were done in double precision on the Cyber 205 computer (the equivalent of about 29 decimal places). Example 5.1. S0 = ~, k~x,2/(k + a)m.
8 316 WALTER GAUTSCHI This series with a = m = X was communicated to the author by Professor P. J. Davis, who encountered it in his study of spirals [4]. We computed S0 from (3.3) (with v = \), (5.1), and (2.1_1), (2.3), for a =.5, X., 2., 4., 8. and m = 1(1)5. The integral in (3.3) was evaluated by «-point Gaussian quadrature rules, which were generated with the help of the recursion coefficients in Table 1 of the Appendix and well-known eigenvalue techniques (see, e.g., [6, 1.3(iv)]). In Table 5.1 we show only the results for m = X ; those for m > 1 are similar, but exhibit somewhat slower convergence. It is evident that, as a increases, convergence of the Gauss quadrature formula slows down considerably. The reason for this is the behavior of the function g0 on the right of (5.1), which for increasing a approaches a discontinuous function (see Figure 5.1). Table 5.1 n-point quadrature approximations to the integral in (3.3) with _ i2 ' m X, a =.5,X.,2.,4.,% a =.5 a=l. a = l a = 4. a = : To achieve better accuracy, when a is large, we proceed as follows. With ao = LaJ denoting the largest integer < a, and a = a0 + a,, where a0 > X, 0 < a, < 1, the summation over all k > X may be "stratified" by letting k = X + Ka0 and summing over all k > 0 for A= X, 2,..., aq. Thus, *o = E k=l 1/2 (k + a0 + a. EE (5.3) =a (m+l/2) "o oo X=l K- YV tl-t^^ -(m+i/2)^ I V^_(K 0 éííe;<«(k + A/flp) 1/2 'Q(K+X + (X + ax)/a0y + ( + Kapï. -1/2 + KOo + Oo + av (K + À/a0) -1/2 ( K -1/2 ' + (X + ax)/aor (l + (X + ax)/a0yi
9 SLOWLY CONVERGENT SERIES AND THEIR SUMMATION Figure 5.1 The function g0 in (5.1) To the inner sum we now apply (3.10) to obtain (5.4) oo E(C=l (* + Xlap) -1/2 (K+X -Í + (X + ax)/a0) a X lße(t)e-w^gm f=j,i] Jo The "effective" parameter in gm_x is now X + ax/aq, a number close to 1, and the coefficient X/a0 in the exponential is bounded by 1. Gauss quadrature applied to the integral in (5.4) should therefore converge quite rapidly, indeed, more so the larger a0\ This is borne out by the results displayed for m = X and a = a0 = 8., 16., 32. in Table 5.2. The correct number of significant decimal digits produced by direct summation of the series, using 1000 terms, is shown in Table 5.3. The numbers in parentheses are the correct digits obtained by applying the e-algorithm with the same number of terms. For the entries marked by an asterisk, we used (5.3), Table 5.2 Approximations to Sñ in (5.3) for m= X, a, = 0, añ = 8., 16., 32., using n-point quadrature in (5.4) a = 8. a= 16. a =
10 318 WALTER GAUTSCHI Table 5.3 Number of correct significant decimal digits in direct (and accelerated) summation of S0 using 1000 terms m a =.5 a = 1. a = 1. a = 4. 2(2) 4(5) 7(9) 10(13) 14(17) 1(2) 4(5) 7(9) 9(12) 13(15) 2(2) 4(5) 6(8) 9(12) 11(14) 2(2) 4(5) 6(8) 8(11)* 11(13)' 1(1) 3(4) 5(7)* 8(10)* 10(12)* (5.4) to verify the number of correct digits. As is evident from Table 5.3, the e-algorithm is only marginally effective on this particular series. Example 5.2. Sx = k=i\ X)k-xk-l/2/(k + a)m. We now apply Gauss quadrature (obtained from the recursion coefficients in Table 2 of the Appendix) relative to the weight function t~l/2tp(t) to the integral in (3.7) (with v = \), using the same values of a and m as in Example 5.1. The results are similar to those in Table 5.1 of Example 5.1, except that convergence is slightly slower. Stratified summation similar to (5.3), (5.4), on the other hand, gives (5.5) ç _ -(ro+1/2) o, - a0 X=l L (W 1/2 (X + (X + ax)/a0) } where (5-6) sx=< -E -1/2 (K + X/a0) K+ X + (X + a,)/aq) K=l OO / i\k-l yk-ir-l(k+x/a0)-l/2 ^(K+X+(X + ax)/a0) if añ is even, if an is odd, that is, by (3.10) and (3. II), I5.V) s, = l - / f{l2e(t)e-[kla )tgm_x(t; X+ax/a0, X/2)dt, a0 even,./o if -Wa0)t t l/2<p(t)-te gm-i(t;i+a,/a0,x/2)dt, a0 odd. < Jo The use of (5.5), with «-point quadrature applied to (5.7), yields results converging at the same speed as those in Table 5.2. Direct summation using 1000 terms yields 2-3 more correct digits than in Table 5.3, but the e-algorithm is now surprisingly effective, giving full accuracy (20 decimals) with only terms!
11 SLOWLY CONVERGENT SERIES AND THEIR SUMMATION 319 1/2, Example 5.3. S0(b) = E?=i(k + bv /(k + b +!)> 0 < Z>< 1. This is another series of interest in P. J. Davis's work on spirals [4]. It can be readily evaluated with the help of Remark 1 in 3, taking m = X, a = b + X in (3.10), and noting (5.1); one finds (5.8) S0(b) = JO t-l/2e(t)e-btg0(f,l, \dt 2 = / f e -btf(vt) V * -1/2 e(t)dt..,. njo y/t Thus, S0(b) is the Laplace transform of the function (F(y/t)/^/t)t~l/2e(t), which has a square root singularity at the origin and poles at integer multiples of 2ni. The last integral in (5.8) is easily computed by Gaussian quadrature relative to the weight function t~l/2e(t). No more than 24 quadrature points are needed to get 14 correct significant decimal digits in the range 0 < b < X. 1/2, Example 5.4. S0 = ^2k=xk~''"/(k + ia), a>0. Here, a = ia, m = X, hence (5.9) S0 = J ri/2e(t)g0(t;ia,^dt. Use of (5.1) and a simple change of variables gives ( 1\ it'2* c, > r~ gjt; ia,-j = ^ erf(-1z), z = Viat, where erf is the error function. Letting -iz = \y/n( X - i)x, i.e., x = ^/2at/n, one finds (cf. [5, eq ]) (5.10) g0(t; ia, i) = y/2/(al)e "* [C(y/2a~tJ i) + is(s/2a~tfji)]. Here, C(x), S(x) denote the Fresnel integrals [5, eqs , 7.3.2]. They can be computed with an accuracy of up to 18 significant digits from the rational Chebyshev approximations provided (on microfiche cards) in [1]. Gaussian quadrature of (5.9), with g0 given by (5.10), yields the results shown in Table 5.4. For each «, the first entry is the real part, the second the imaginary part of the Gauss approximation to S0. Convergence is seen to deteriorate rapidly with increasing a, which is to be expected in view of the highly oscillatory behavior of g0 in (5.10) when a is large. The device of stratified summation, used successfully in Examples 5.1 and 5.2, however, can also be applied here,
12 320 WALTER GAUTSCHI Table 5.4 n-point quadrature approximations to the integral in (3.3) with " = i m = X, a = ia, a =.5, 1., 2., 4., a = a = 2. a = 4. a = and gives, with aq= [a\, a = aq + ax, aq > 1, 0 < a, < 1, ;5.n) =qo~3/2e{ r ri'e(t)e-^>gat;i\x x=i + ^\,^\dt, KM)3/2 l + ííao + aj/aj Using (5.11) instead of (5.9) yields the results shown in Table 5.5. Direct summation using 1000 terms gives 1-2 correct decimal digits in the real part, and 4-5 in the imaginary part. The epsilon algorithm produces no more than one additional correct digit. Example 5.5. S{=\Z^{(-^k'xl2l(k + ia), a>0. Applying Gauss quadrature to (5.12) -fjo t 1/2<p(t)-tgJt;ia, \)dt, with g0 as in (5.10), or to formulae analogous to those in (5.5) (5.7), yields results similar in quality to those in Table 5.4 (but converging at a slightly slower rate) and to those in Table 5.5. As in Example 5.2, here too, the e-algorithm produces full accuracy (20 decimals) with as few as terms.
13 SLOWLY CONVERGENT SERIES AND THEIR SUMMATION 321 Table 5.5 Approximations to Sc for a, = 0, a0 = 8., 16., 32., using n-point quadrature for the integral in (5.11) a= 16. a = Acknowledgment The stimulus to write this paper came from Professor P. J. Davis who showed the author the two series in Example 5.1 (with a = m = X) and Example 5.3 and their interesting connections with analytic properties of certain spirals. The author thanks Professor Davis for bringing these matters to his attention. Appendix Coefficients ak, ßk in the recurrence relation (3.14) for the (monic) polynomials 7t. ( ; wx) and &.( ; w2) orthogonal on [0, oo] with respect to the - i-1/2 V2, weight functions wx(t) = t e(t) and w2(t) = t ' <p(t), where e and tp are the Einstein and Fermi functions, respectively. Table 1 Recursion coefficients for the polynomials {nk(-;wx)} lb 16! alpha(k) D D D D D D D D D D D D D D D D D D D beta(k) D D D D D D D D D D D D D D D D D D D+03
14 322 WALTER GAUTSCHI Table 1 (continued) k alpha(k) beta(k) 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+03 O D D D D D D D D D D D D D D D D OD 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+04
15 SLOWLY CONVERGENT SERIES AND THEIR SUMMATION 323 Table 2 Recursion coefficients for the polynomials {nk(-; w2)} k alpha(k) beta(k) 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 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+04
16 324 WALTER GAUTSCHI Table 2 (continued) k alpha(k) beta(k) D D D D D D D D D lOD 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+04 Bibliography 1. W. J. Cody, Chebyshev approximations for the Fresnel integrals, Math. Comp. 22 (1968), Loose microfiche suppl. A1-B4. 2. W. J. Cody, K. A. Paciorek, and H. C. Thacher, Jr., Chebyshev approximations for Dawson's integral, Math. Comp. 24 (1970), P. J. Davis, Gamma function and related functions, Handbook of Mathematical Functions, Chapter 6 (M. Abramowitz and I. A. Stegun, eds.), NBS Appl. Math. Series, vol. 55, U.S. Government Printing Office, Washington, D.C., _, Spirals: From Theodorus of Cyrene to Meta-Chaos, The Hedrick Lectures, Math. Assoc. Amer., August, W. Gautschi, Error function and Fresnel integrals, Handbook of Mathematical Functions, Chapter 7 (M. Abramowitz and I. A. Stegun, eds.), NBS Appl. Math. Series, vol. 55, U.S. Government Printing Office, Washington, D.C., _, Computational aspects of orthogonal polynomials, Orthogonal Polynomials Theory and Practice (P. Nevai, ed.), NATO ASI Series, Series C: Mathematical and Physical Sciences, vol. 294, Kluwer, Dordrecht, 1990, pp W. Gautschi and G. V. Milovanovic, Gaussian quadrature involving Einstein and Fermi functions with an application to summation of series, Math. Comp. 44 (1985), P. Henrici, Applied and computational complex analysis, vol. 1, Wiley, New York, Department of Computer Sciences, Purdue University, West Lafayette, Indiana address: wxg@cs.purdue.edu
Summation 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 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 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 informationNumerische MathemalJk
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
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 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 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 informationCalculation of the Ramanujan t-dirichlet. By Robert Spira
MATHEMATICS OF COMPUTATION, VOLUME 27, NUMBER 122, APRIL, 1973 Calculation of the Ramanujan t-dirichlet Series By Robert Spira Abstract. A method is found for calculating the Ramanujan T-Dirichlet series
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 informationw = X ^ = ^ ^, (i*i < ix
A SUMMATION FORMULA FOR POWER SERIES USING EULERIAN FRACTIONS* Xinghua Wang Math. Department, Zhejiang University, Hangzhou 310028, China Leetsch C. Hsu Math. Institute, Dalian University of Technology,
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 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 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 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 informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
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 informationThe Perrin Conjugate and the Laguerre Orthogonal Polynomial
The Perrin Conjugate and the Laguerre Orthogonal Polynomial In a previous chapter I defined the conjugate of a cubic polynomial G(x) = x 3 - Bx Cx - D as G(x)c = x 3 + Bx Cx + D. By multiplying the polynomial
More informationAnalytic Continuation of the Theodorus Spiral
Analytic Continuation of the Theodorus Spiral Jörg Waldvogel Seminar für Angewandte Mathematik, ETH Zürich E-mail: waldvogel@math.ethz.ch October 9, 8 to February 9 Abstract The remarkable classical pattern
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 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 informationSums of Products of Bernoulli Numbers*
journal of number theory 60, 234 (996) article no. 00 Sums of Products of Bernoulli Numbers* Karl Dilcher Department of Mathematics, Statistics, and Computing Science, Dalhousie University, Halifax, Nova
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 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 informationA matrix over a field F is a rectangular array of elements from F. The symbol
Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted
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 informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More informationThe Hardy-Littlewood Function
Hardy-Littlewood p. 1/2 The Hardy-Littlewood Function An Exercise in Slowly Convergent Series Walter Gautschi wxg@cs.purdue.edu Purdue University Hardy-Littlewood p. 2/2 OUTLINE The Hardy-Littlewood (HL)
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 informationWeek 9-10: Recurrence Relations and Generating Functions
Week 9-10: Recurrence Relations and Generating Functions April 3, 2017 1 Some number sequences An infinite sequence (or just a sequence for short is an ordered array a 0, a 1, a 2,..., a n,... of countably
More informationSome Fun with Divergent Series
Some Fun with Divergent Series 1. Preliminary Results We begin by examining the (divergent) infinite series S 1 = 1 + 2 + 3 + 4 + 5 + 6 + = k=1 k S 2 = 1 2 + 2 2 + 3 2 + 4 2 + 5 2 + 6 2 + = k=1 k 2 (i)
More information5.4 Bessel s Equation. Bessel Functions
SEC 54 Bessel s Equation Bessel Functions J (x) 87 # with y dy>dt, etc, constant A, B, C, D, K, and t 5 HYPERGEOMETRIC ODE At B (t t )(t t ), t t, can be reduced to the hypergeometric equation with independent
More informationMATHEMATICAL FORMULAS AND INTEGRALS
MATHEMATICAL FORMULAS AND INTEGRALS ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom Academic Press San Diego New York Boston London
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationA Symbolic Operator Approach to Several Summation Formulas for Power Series
A Symbolic Operator Approach to Several Summation Formulas for Power Series T. X. He, L. C. Hsu 2, P. J.-S. Shiue 3, and D. C. Torney 4 Department of Mathematics and Computer Science Illinois Wesleyan
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 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 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 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 informationApproximate Integration Formulas for Ellipses
322 N. LEE AND A. H. STR0UD 5. Acknowledgements. The author wishes to record his gratitude and sincere thanks to Mr. I. M. Khabaza formerly of the University of London Computer Unit3 for his constant interest
More informationAsymptotic Error Estimates for the Gauss Quadrature Formula
Asymptotic Error Estimates for the Gauss Quadrature Formula By M. M. Chawla and M. K. Jain 1. Introduction. Classical error estimates for the Gauss quadrature formula using derivatives can be used, but
More informationOn the N-tuple Wave Solutions of the Korteweg-de Vnes Equation
Publ. RIMS, Kyoto Univ. 8 (1972/73), 419-427 On the N-tuple Wave Solutions of the Korteweg-de Vnes Equation By Shunichi TANAKA* 1. Introduction In this paper, we discuss properties of the N-tuple wave
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 informationEXISTENCE THEOREMS FOR SOLUTIONS OF DIFFERENTIAL EQUATIONS OF NON-INTEGRAL ORDER*
100 EVERETT PITCHER AND W. E. SEWELL [February EXISTENCE THEOREMS FOR SOLUTIONS OF DIFFERENTIAL EQUATIONS OF NON-INTEGRAL ORDER* EVERETT PITCHER AND W. E. SEWELL 1. Introduction. In this paper we prove
More informationGegenbauer Matrix Polynomials and Second Order Matrix. differential equations.
Gegenbauer Matrix Polynomials and Second Order Matrix Differential Equations Polinomios Matriciales de Gegenbauer y Ecuaciones Diferenciales Matriciales de Segundo Orden K. A. M. Sayyed, M. S. Metwally,
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 informationALTERNATIVE DERIVATION OF SOME REGULAR CONTINUED FRACTIONS
ALTERNATIVE DERIVATION OF SOME REGULAR CONTINUED FRACTIONS R. F. C. WALTERS (Received 21 July 1966, revised 20 February 1967) In this paper we find an expression for g* as the limit of quotients associated
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,
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 informationPolyexponentials. Khristo N. Boyadzhiev Ohio Northern University Departnment of Mathematics Ada, OH
Polyexponentials Khristo N. Boyadzhiev Ohio Northern University Departnment of Mathematics Ada, OH 45810 k-boyadzhiev@onu.edu 1. Introduction. The polylogarithmic function [15] (1.1) and the more general
More informationVALUES OF THE LEGENDRE CHI AND HURWITZ ZETA FUNCTIONS AT RATIONAL ARGUMENTS
MATHEMATICS OF COMPUTATION Volume 68, Number 228, Pages 1623 1630 S 0025-5718(99)01091-1 Article electronically published on May 17, 1999 VALUES OF THE LEGENDRE CHI AND HURWITZ ZETA FUNCTIONS AT RATIONAL
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 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 informationRESEARCH ARTICLE. Spectral Method for Solving The General Form Linear Fredholm Volterra Integro Differential Equations Based on Chebyshev Polynomials
Journal of Modern Methods in Numerical Mathematics ISSN: 29-477 online Vol, No, 2, c 2 Modern Science Publishers wwwm-sciencescom RESEARCH ARTICLE Spectral Method for Solving The General Form Linear Fredholm
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationS \ - Z) *! for j = 1.2,---,k + l. i=i i=i
RESTRICTED COMPOSITIONS S. G. MOHANTY McMaster University, Hamilton, Ontario, Canada As a continuation of [6] and [7], this paper deals with a restricted set of compositions of an integer (to be defined
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 informationa 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2.
Chapter 1 LINEAR EQUATIONS 11 Introduction to linear equations A linear equation in n unknowns x 1, x,, x n is an equation of the form a 1 x 1 + a x + + a n x n = b, where a 1, a,, a n, b are given real
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 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 informationA NOTE ON A BASIS PROBLEM
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 51, Number 2, September 1975 A NOTE ON A BASIS PROBLEM J. M. ANDERSON ABSTRACT. It is shown that the functions {exp xvx\v_. form a basis for the
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 informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More informationA NEW METHOD OF INVERSION OF THE LAPLACE TRANSFORM* ATHANASIOS PAPOULIS. Polytechnic Institute of Brooklyn
405 A NEW METHOD OF INVERSION OF THE LAPLACE TRANSFORM* BY ATHANASIOS PAPOULIS Polytechnic Institute of Brooklyn Introduction. In determining a function r(t) from its Laplace transform R(p) R(p) = [ e"'r(t)
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationswapneel/207
Partial differential equations Swapneel Mahajan www.math.iitb.ac.in/ swapneel/207 1 1 Power series For a real number x 0 and a sequence (a n ) of real numbers, consider the expression a n (x x 0 ) n =
More informationNONNEGATIVE AND SKEW-SYMMETRIC PERTURBATIONS OF A MATRIX WITH POSITIVE INVERSE
mathematics of computation volume 54, number 189 january 1990, pages 189-194 NONNEGATIVE AND SKEW-SYMMETRIC PERTURBATIONS OF A MATRIX WITH POSITIVE INVERSE GIUSEPPE BUFFONI Abstract. Let A be a nonsingular
More informationAdomian Decomposition Method with Laguerre Polynomials for Solving Ordinary Differential Equation
J. Basic. Appl. Sci. Res., 2(12)12236-12241, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal of Basic and Applied Scientific Research www.textroad.com Adomian Decomposition Method with Laguerre
More informationSOME IDENTITIES INVOLVING THE FIBONACCI POLYNOMIALS*
* Yi Yuan and Wenpeeg Zhang Research Center for Basic Science, Xi'an Jiaotong University, Xi'an Shaanxi. P.R. China (Submitted June 2000-Final Revision November 2000) 1. INTROBUCTION AND RESULTS As usual,
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 informationClosed-form Second Solution to the Confluent Hypergeometric Difference Equation in the Degenerate Case
International Journal of Difference Equations ISS 973-669, Volume 11, umber 2, pp. 23 214 (216) http://campus.mst.edu/ijde Closed-form Second Solution to the Confluent Hypergeometric Difference Equation
More informationSome Observations on Interpolation in Higher Dimensions
MATHEMATICS OF COMPUTATION, VOLUME 24, NUMBER 111, JULY, 1970 Some Observations on Interpolation in Higher Dimensions By R. B. Guenther and E. L. Roetman Abstract. This paper presents a method, based on
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 informationRational Chebyshev Approximations for the Bickley Functions Kin(x)
MATHEMATICS OF COMPUTATION, VOLUME, NUMBER JULY, PAGES - Rational Chebyshev Approximations for the Bickley Functions Kin(x) By J. M. Blair, C. A. Edwards and J. H. Johnson Abstract. This report presents
More informationMATHEMATICAL FORMULAS AND INTEGRALS
HANDBOOK OF MATHEMATICAL FORMULAS AND INTEGRALS Second Edition ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom ACADEMIC PRESS A Harcourt
More informationNumerical Quadrature Rules for Some Infinite Range Integrals
MATHEMATICS OF COMPUTATION VOLUME 38, NUMBER 157 JANUARY 1982 Numerical Quadrature Rules for Some Infinite Range Integrals By Avram Sidi Abstract. Recently the present author has given a new approach to
More informationNumerical Quadrature over a Rectangular Domain in Two or More Dimensions
Numerical Quadrature over a Rectangular Domain in Two or More Dimensions Part 2. Quadrature in several dimensions, using special points By J. C. P. Miller 1. Introduction. In Part 1 [1], several approximate
More informationL U C A S SEQUENCES AND FUNCTIONS O F A 3-BY-3 M A T R I X
L U C A S SEQUENCES AND FUNCTIONS O F A 3-BY-3 M A T R I X R. S. Melham School of Mathematical Sciences, University of Technology, Sydney PO Box 123, Broadway, NSW 2007 Australia {SubmittedMay 1997-Final
More informationArithmetic properties of lacunary sums of binomial coefficients
Arithmetic properties of lacunary sums of binomial coefficients Tamás Mathematics Department Occidental College 29th Journées Arithmétiques JA2015, July 6-10, 2015 Arithmetic properties of lacunary sums
More informationMULTISECTION METHOD AND FURTHER FORMULAE FOR π. De-Yin Zheng
Indian J. pure appl. Math., 39(: 37-55, April 008 c Printed in India. MULTISECTION METHOD AND FURTHER FORMULAE FOR π De-Yin Zheng Department of Mathematics, Hangzhou Normal University, Hangzhou 30036,
More informationON THE OSCILLATION OF THE DERIVATIVES OF A PERIODIC FUNCTION
ON THE OSCILLATION OF THE DERIVATIVES OF A PERIODIC FUNCTION BY GEORGE PÖLYA AND NORBERT WIENER 1. Let f(x) be a real valued periodic function of period 2ir defined for all real values of x and possessing
More information-P" -p. and V n = a n + ft\ (1.2)
R. S. Melham School of Mathematical Sciences, University of Technology, Sydney, PO Box 123, Broadway, NSW 27 Australia {Submitted November 1997-Final Revision April 1998) 1. INTRODUCTION Define the sequences
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
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 information/ dr/>(t)=co, I TT^Ti =cf' ; = l,2,...,» = l,...,p.
PROCEEDINGS of the AMERICAN MATHEMATICAL SOCIETY Volume 02, Number, January 988 UNIQUE SOLVABILITY OF AN EXTENDED STIELTJES MOMENT PROBLEM OLAV NJÂSTAD (Communicated by Paul S. Muhly) ABSTRACT. Let ai,...,ap
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 informationPolynomial Formula for Sums of Powers of Integers
Polynomial Formula for Sums of Powers of Integers K B Athreya and S Kothari left K B Athreya is a retired Professor of mathematics and statistics at Iowa State University, Ames, Iowa. right S Kothari is
More informationSeries Solutions. 8.1 Taylor Polynomials
8 Series Solutions 8.1 Taylor Polynomials Polynomial functions, as we have seen, are well behaved. They are continuous everywhere, and have continuous derivatives of all orders everywhere. It also turns
More informationMOMENTS OF THE LIFETIME OF CONDITIONED BROWNIAN MOTION IN CONES BURGESS DAVIS AND BIAO ZHANG. (Communicated by Lawrence F. Gray)
proceedings of the american mathematical society Volume 121, Number 3, July 1994 MOMENTS OF THE LIFETIME OF CONDITIONED BROWNIAN MOTION IN CONES BURGESS DAVIS AND BIAO ZHANG (Communicated by Lawrence F.
More informationEXPANSION OF ANALYTIC FUNCTIONS IN TERMS INVOLVING LUCAS NUMBERS OR SIMILAR NUMBER SEQUENCES
EXPANSION OF ANALYTIC FUNCTIONS IN TERMS INVOLVING LUCAS NUMBERS OR SIMILAR NUMBER SEQUENCES PAUL F. BYRD San Jose State College, San Jose, California 1. INTRODUCTION In a previous article [_ 1J, certain
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
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 informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationPolynomial Solutions of Nth Order Nonhomogeneous Differential Equations
Polynomial Solutions of th Order onhomogeneous Differential Equations Lawrence E. Levine Ray Maleh Department of Mathematical Sciences Stevens Institute of Technology Hoboken,J 07030 llevine@stevens-tech.edu
More informationOn Systems of Diagonal Forms II
On Systems of Diagonal Forms II Michael P Knapp 1 Introduction In a recent paper [8], we considered the system F of homogeneous additive forms F 1 (x) = a 11 x k 1 1 + + a 1s x k 1 s F R (x) = a R1 x k
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 informationExplicit evaluation of the transmission factor T 1. Part I: For small dead-time ratios. by Jorg W. MUller
Rapport BIPM-87/5 Explicit evaluation of the transmission factor T (8,E) Part I: For small dead-time ratios by Jorg W. MUller Bureau International des Poids et Mesures, F-930 Sevres Abstract By a detailed
More informationZeros of the Hankel Function of Real Order and of Its Derivative
MATHEMATICS OF COMPUTATION VOLUME 39, NUMBER 160 OCTOBER 1982, PAGES 639-645 Zeros of the Hankel Function of Real Order and of Its Derivative By Andrés Cruz and Javier Sesma Abstract. The trajectories
More informationNIL DERIVATIONS AND CHAIN CONDITIONS IN PRIME RINGS
proceedings of the american mathematical society Volume 94, Number 2, June 1985 NIL DERIVATIONS AND CHAIN CONDITIONS IN PRIME RINGS L. O. CHUNG AND Y. OBAYASHI Abstract. It is known that in a prime ring,
More information