To appear: Math. Comp LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS
|
|
- Sharyl Tate
- 6 years ago
- Views:
Transcription
1 To appear: Math Comp 21 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS GEORGE BOROS AND VICTOR H MOLL Abstract We present a rational version of the classical Landen transformation for elliptic integrals This is employed to obtain explicit closed-form expressions for a large class of integrals of even rational functions and to develop an algorithm for numerical integration of these functions 1 Introduction We consider the space of even rational functions of degree 2p { E 2p := R(z = P (z p 1 P (z := b z 2(p 1 and Q(z := Q(z = } p a z 2(p with positive real coefficients a, b R + normalized by the condition a = a p = 1, the space E := of normalized even rational functions, and the set of 2p 1 parameters p=1 E 2p P 2p := {a 1,, a p 1 ; b,, b p 1 } We describe an algorithm to determine, as a function of the parameter set P 2p, a closed-form expression of the integral (11 I := R(z dz for a large class of functions R E The function R is called symmetric if its denominator Q satisfies Q(1/z = z 2p Q(z This is equivalent to its coefficients being palindromic, ie a j = a p j for 1 j p The class of symmetric functions plays a crucial role in this algorithm Define E s 2p := {R E 2p den(r is symmetric} (where den(r denotes the denominator of R, the class of rational functions with symmetric denominators of degree 2p, and E s := E s 2p Date: April 25, Mathematics Subject Classification Primary 33 Key words and phrases Rational functions, Landen transformation, Integrals The second author was supported in part by NSF Grant DMS p=1 =
2 2 GEORGE BOROS AND VICTOR H MOLL For m N define and E m,s 2p E m 2p := {R E (den(r 1/(m+1 is even, symmetric of degree 2p} := E m 2p E s 2p, so a function R E m,s 2p can be written in the form R(z = P (z Q m+1 (z, where P (z is an even polynomial and Q(z is an even symmetric polynomial of degree 2p The method of partial fractions gives (in principle the value of I in terms of the roots of Q Symbolic computations yield either a closed-form answer, an expression in terms of the roots of an associated polynomial, or the integral returned unevaluated The algorithm described here allows only algebraic operations on elementary functions and changes of variables of the same type In particular, we exclude the solution of algebraic equations of degree higher than 2 We say that a rational function R E is computable if its integral can be evaluated by our algorithm The first step in the algorithm is to consider symmetric rational functions In Section 2 we prove a reduction formula, ie a map F p : E s 2p R p that reduces the computability of the integral of the symmetric function R to that of one of degree 1 2 deg(r Here R p is the space of rational functions with denominator of degree p These new functions are not necessarily symmetric, and this is the main limitation of our algorithm The details of F p require the evaluation of some binomial sums which are presented in Appendix A The classical Wallis formula dz (z m+1 = ( π 2m 2 2m+1 m shows that every R E m 2 is computable The reduction formula now implies that every R E m 4 is computable This is described in Section 3 The computability of R E 4 is also a consequence of the classical theory of hypergeometric functions; the details are given in [2] In Section 4 we compute the integral of every function in E m,s 8, where the reduction method expresses these integrals in terms of functions in E m 4 The algorithm does not, in general, provide a value for the integral of a nonsymmetric function of degree 8 The final piece of the algorithm is described in Section 5: for R E 2p, the symmetrization of its denominator produces a (symmetric rational function in E 4p with the same integral as R The reduction formula now yields a new function in E 2p with the same integral as R We thus obtain a map T 2p : E 2p E 2p such that (12 R(z dz = T 2p (R(z dz In particular, the class of computable rational functions of degree 2p is invariant under forward and bacward iteration of T 2p This map is the rational analog of the original Landen transformation for elliptic integrals described in [4, 13] The map T 2p can also be interpreted as a map on the coefficients Φ 2p : O + 2p O+ 2p where O + 2p = R + p 1 R + p
3 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS 3 The case of Φ 6 is described in detail in Section 6 and is given explicitly by (13 a a 1 + 5a 2 + a 1 a 2 (a 1 + a /3 a 2 a 1 + a (a 1 + a /3 b b + b 1 + b 2 (a 1 + a /3 b 1 b (a b 1 + b 2 (a a 1 + a b + b 2 b 2 (a 1 + a /3 Let x := (a 1, a 2 ; b, b 1, b 2 Then Φ 6 : O + 6 O+ 6 is iterated to produce a sequence x n+1 := Φ 6 (x n of points in O + 6 that yield a sequence of rational functions with constant integral We have proved in [3] the existence of L R +, depending upon the initial point x = (a 1, a 2 ; b, b 1, b 2, such that x n (3, 3; L, 2L, L Thus (14 b z 4 + b 1 z 2 + b 2 z 6 + a 1 z 4 + a 2 z dz = L(x π 2 This establishes a numerical method to compute the integral in (14 Numerical studies on integrals of even degree 2p suggest the existence of a limiting value L = L(x such that the sequence x n := Φ 2p (x n satisfies x n (( ( ( ( ( p p p p 1 p 1,,, ; L, 1 2 p 1 1 L,, ( p 1 p 1 L The integral of the original rational function is thus π 2 L(x The proof of convergence remains open for p 4 Examples are given in Section 7 The most important issues left unresolved in this paper are the convergence of the iteration of the map Φ 2p discussed above and the geometric interpretation of the Landen transformation T 2p Finally, the question of integration of odd functions has not been addressed at all Some history The problem of integration of rational functions R(z = P (z/q(z was considered by J Bernoulli in the 18 th century He completed the original attempt by Leibniz of a general partial fraction decomposition of R(z The main difficulty associated with this procedure is to obtain a complete factorization of Q(z over R Once this is nown the partial fraction decomposition of R(z can be computed The fact is that the primitive of a rational function is always elementary: it consists of a new rational function (its rational part and the logarithm of a second rational function (its transcendental part In his classic monograph [9] G H Hardy states: The solution of the problem (of definite integration in the case of rational functions may therefore be said to be complete; for the difficulty with regard to the explicit solution of algebraical equations is one not of inadequate nowledge but of proved impossibility He goes on to add: It appears from the preceding paragraphs that we can always find the rational part of the integral, and can find the complete integral if we can find the roots of Q(z =
4 4 GEORGE BOROS AND VICTOR H MOLL In the middle of the last century Hermite [1] and Ostrogradsy [15] developed algorithms to compute the rational part of the primitive of R(z without factoring Q(z More recently Horowitz [11] rediscovered this method and discussed its complexity The problem of computing the transcendental part of the primitive was finally solved by Lazard and Rioboo [12], Rothstein [17] and Trager [18] For detailed descriptions and proofs of these algorithms the reader is referred to [5] and [6] 2 The reduction formula In this section we present a map F p : E m,s 2p E m p that is the basis of the integration algorithm described in Section 5 The proof is elementary and the binomial sums discussed in the Appendix are employed Let D p (z be the general symmetric polynomial of degree 4p We express the integral of z 2n /Dp m+1 as a linear combination of integrals of z /Ep m+1 where E p is a polynomial of degree 2p whose coefficients are determined by those of D p Theorem 21 Let m, n, p N Define (21 and (22 D p (d 1, d 2,, d p ; z = E p (d 1, d 2,, d p ; z = + p+1 j=1 p i=1 d j p d p+1 (z 2 + z 4p 2 = z 2p + p i+1 2 2i 1 z 2(p i for d i R + Then for n (m + 1p 1, (m+1p n 1 2 m (23 j= j=1 z 2n dz (D p (d 1,, d p ; z m+1 = ( (m + 1p n 1 + j 4 j j + i 1 i ( j + 2i 2 j 1 z 2((m+1p 1 j ( E p (d 1,, d p ; z d j+i, m+1 dz, and for (m + 1p 1 < n < 2p(m we employ the symmetry rule (24 N n,p = N 2p(m+1 1 n,p Proof First observe that (24 follows from the change of variable z 1/z Now consider z 2n dz (25 N n,p (d 1,, d p ; m := ( p = d p+1 (z 4p 2 + z 2 m+1 for n (m + 1p 1 The substitution z = tan θ yields N n,p = π/2 (1 C 2 n C 4(m+1p 2n 2 dθ ( p = d p+1 {(1 C 2 2p C 2 + (1 C 2 C 4p 2 } m+1,
5 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS 5 where C = cos θ Letting ψ = 2 θ and D = cos ψ = 2C 2 1 then gives π (1 D 2 n (1 + D 2(m+1p 2n 1 dψ N n,p = ( p = d p+1 (1 D 2 {(1 D 2p 2 + (1 + D 2p 2 } m+1 Now observe that the integrals of the odd powers of cosine vanish when we expand (1 + D 2(m+1p 2n 1, producing π/2 (1 D 2 n (m+1p n 1 ( 2(m+1p 2n 1 N n,p = 2 m j= D dθ } m+1 D { p = d p+1 (1 D 2 p j= A second double angle substitution ϕ = 2ψ gives π (1 E n (m+1p n 1 N n,p = 2 m j= ( 2(m+1p 2n 1 { p = d p+1 (1 E p j= ( 2p 2 2 p(m+1 n j 1 (1 + E j dϕ ( 2p 2 } m+1, 2 p j (1 + E j where E = cos ϕ = 2D 2 1 The change of variable z = tan(ϕ/2 then yields z 2n (m+1p n 1 ( 2(m+1p 2n 1 N n,p = 2 m j= (1 + z 2 (m+1p n j 1 dz { p ( = d p ( p+1 z 2 2p 2 } m+1 j= (1 + z2 p j (26 Finally, we modify (26 using Lemma A2 and Lemma A4 with N = (m+1p n 1 to produce (23 Note that the previous theorem associates to each rational function of symmetric denominator R 1 (z = a new rational function (m+1p 1 R 2 (z = 2 m such that n= b s z 2s + b s 1 z 2(s b ( z 4p + d p z 4p d 1 z 2p m+1 (m+1p n 1 b n j= R 1 (z dz = ( (m + 1p n 1 + j 4 j R 2 (z dz z 2((m+1p 1 j (E p (d 1,, d p ; z m+1 3 The quartic case In this section we describe the computability of rational functions R E m 4 These are functions of the form P (z R(z = (z 4 + 2az m+1 where P (z is an even polynomial of degree 4m+2 Observe that the normalization a = a 2 = 1 maes the denominator of R automatically symmetric It suffices to evaluate z 2n dz (31 N n,4 (d 1 ; m := (z 4 + 2d 1 z m+1
6 6 GEORGE BOROS AND VICTOR H MOLL where n 2m+1 is required for convergence From (24 we have N n,4 (d 1 ; m = N 2m 1 n,4 (d 1 ; m, so we may assume n m We now employ Theorem 21 to obtain a closed form expression for N n,4 (d 1 ; m Theorem 31 Let m N and assume n m Then (32 N n,4 (d 1 ; m := z 2n dz (z 4 + 2d 1 z m+1 = m n π 2 3m+3/2 (1 + a ( ( 2m 1 m n + j 2 j (1 + d m+1/2 1 j m j j= For m + 1 n 2m + 1 we have (33 z 2n dz (z 4 + 2d 1 z m+1 = n m 1 π 2 3m+3/2 (1 + d 1 ( ( 2m 1 m n + j 2 j (1 + d m+1/2 1 j m j j= j= ( j ( 1 m j ( j Proof We apply the result of the Theorem 21 with D 1 (d 1 ; z = z 4 + 2d 1 z and E 1 (d 1 ; z = (1 + d 1 z 2 + 2, so that z 2n m n dz ( m n + j (z 4 + 2d 1 z m+1 = 2 m 4 j z 2(m j dz ((1 + d 1 z m+1 The change of variable u = (1 + d 1 z 2 /2 then yields z 2(m j dz ( (1 + d 1 z m+1 = 2 (j+3/2 (1 + d 1 m+j 1/2 where we have used ( 2m = π m j u r 1/2 (1 + u s du = π 2 2(s 1 ( 2r r ( j ( m j ( 2(s r 1 s r 1 ( 1 m j u m j 1/2 du (1 + u m (2m+j+3/2 (1 + d 1 (m j+1/2, ( 1 s 1 r The algorithm also requires a scaled version of N,4 (d 1 ; m Corollary 32 Let b >, c >, a > bc, m N, and n m Define Then for n m, N n,4 (a, b, c; m = π N n,4 (a, b, c; m := z 2n dz (bz 4 + 2az 2 + c m+1 ( c(c/b m n { 8(a + bc} 2m+1 1/2 (34 m n = 2 ( 2m 2 m ( m n + 2 ( 2 ( m 1 ( a + 1, bc
7 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS 7 and for m + 1 n 2m + 1, ( { N n,4 (a, b, c; m = π c(c/b m n 8(a + 2m+1 1/2 bc} (35 n m 1 = 2 ( 2m 2 m ( m n + 2 ( 2 Proof Let n m The substitution u = z(b/c 1/4 yields (36 N n,4 (a, b, c; m = 1 c m n/2+3/4 b n/2+1/4 N n,4 ( m ( a bc ; m 1 ( a + 1 bc so (34 then follows from Theorem 31 From (24 we have N n,4 (a, b, c; m = N 2m+1 n,4 (c, b, a; m for m + 1 n 2m + 1, giving (35 4 The symmetric case of degree 8 In this section we prove the computability of the set E m,s 8 of symmetric rational functions with denominator of degree 8 and establish an explicit formula for the integral N n,8 (a 1, a 2 ; m =, z 2n dz (z 8 + a 2 z 6 + 2a 1 z 4 + a 2 z m+1 where n 4m + 3 is required for convergence Observe that (24 reduces the discussion to the case n 2m + 1 The expression (22, with p = 2, produces E 2 (a 1, a 2 ; z = (1 + a 1 + a 2 z 4 + 2(a 2 + 4z Theorem 41 Every function in E m,s 8 is computable More specifically, define c 1 := a 2 + 4, c 2 := 1 + a 1 + a 2, and t,j (m, n; a 1, a 2 := π2 (3m+2++j/2 c (m j/2 2 (c 1 + 8c 2 j m 1/2 ( ( ( ( ( 1 4m n + 2 2m m + j m n m j j j Then for m + 1 n 2m + 1, 1 + a 1 + a 2 > and a > 8 8(1 + a 1 + a 2, z 2n dz (z 8 + a 2 z 6 + 2a 1 z 4 + a 2 z m+1 = and for n m, m m =n j= =n t,j (m, n; a 1, a 2 + 2m+1 =n m 1 j= z 2n dz (z 8 + a 2 z 6 + 2a 1 z 4 + a 2 z m+1 = 2m+1 =m+1 m 1 j= t,j (m, n; a 1, a 2 t,j (m, n; a 1, a 2, Proof The reduction formula yields z 2n dz (z 8 + a 2 z 6 + 2a 1 z 4 + a 2 z m+1 = 2m+1 ( 4m n m z 2 (41 dz n (c 2 z 4 + 2c 1 z m+1
8 8 GEORGE BOROS AND VICTOR H MOLL We then use Corollary 32 to evaluate (41 5 A sequence of Landen transformations The transformation theory of elliptic integrals was initiated by Landen in 1771 He proved the invariance of the function (51 under the transformation (52 ie that (53 G(a, b := π/2 a 1 = (a + b/2 d θ a2 cos 2 θ + b 2 sin 2 θ G(a 1, b 1 = G(a, b b 1 = ab, Gauss [7] rediscovered this invariance while numerically calculating the length of a lemniscate An elegant proof of (53 is given by Newman in [14] Here, the substitution x = b tan θ converts 2G(a, b into the integral of [ (a 2 + x 2 (b 2 + x 2 ] 1/2 over R; the change of variable t = (x ab/x/2 then completes the proof The Gauss-Landen transformation can be iterated to produce a double sequence (a n, b n such that a n b n < 2 n It follows that a n and b n converge to a common limit, the so-called arithmetic-geometric mean of a and b, denoted by AGM(a, b Passing to the limit in G(a, b = G(a n, b n produces (54 π 2AGM(a, b = π/2 d θ a2 cos 2 θ + b 2 sin 2 θ The reader is referred to [4] and [13] for details The goal of this section is to produce a map T 2p : E 2p E 2p that preserves the integral, ie (55 R(z dz = T 2p (R(z dz This map is the rational analog of the original Landen transformation (52 Theorem 51 Let R(z = P (z/q(z with (56 p 1 P (z = b j z 2(p 1 j and Q(z = j= p a j z 2(p j j= Define a j = for j > p, b j = for j > p 1, j (57 d p+1 j = a p a j for p 1, (58 d 1 = 1 2 = p = a 2 p, (59 c j = 2p 1 = a j b p 1 j+
9 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS 9 for j 2p 1, and also (51 Let (511 α p (i = for 1 i p 1, and (512 { 2 2i 1 p+1 i =1 1 + p +i 1 i =1 d if i = a + i = b + i = Q(1 2i/p+1/p 2 ( +2i 2 1 d+i if 1 i p α p (i 2 2i Q(1 2(1 i/p [ p 1 i ( ] p 1 + i (c + c 2p 1 2i = for i p 1 Finally, define the polynomials (513 p 1 P + (z = b + i z2(p 1 i and Q + (z = = Then T 2p (R(z := P + (z/q + (z satisfies (55, ie (514 P (z Q(z dz = p a + i z2(p i = P + (z Q + (z dz Proof The first step is to convert the polynomial Q(z to its symmetric form: with Then I := P (z Q(z dz = C(z = P (z z 2p Q(1/z := D(z = Q(z z 2p Q(1/z := I = 2p 1 = 2p 1 = C(z D(z dz c z 2 p d p+1 (z 2 + z 2(2p = c z 2 dz Q s (z Now employ the reduction formula in Section 2 to evaluate j= L := z 2 dz Q s (z Observe that one needs to evaluate L only for p 1 Indeed, the usual symmetry rule yields L = L 2p 1 The reduction formula now gives p 1 ( p 1 + j L = 2 z 2(p 1 j dz p i= α p(iz 2(p i = p 1 α p (p j=1 ( p + j 2 2 λ 2(j 1 2p +1 2 with α p (i as in (51 and λ = [α p (p/α p (] 1/2p z 2(p j dz p i= b+ i z2(p i
10 1 GEORGE BOROS AND VICTOR H MOLL Note The extension of this transformation to the case of P (z Q m+1 (z dz requires explicit formulae for the coefficients of P (z ( z 2p Q(1/z m+1 and Q m+1 (z ( z 2p Q(1/z m+1 An algorithm for integration Let x = (a, b with a = (a 1,, a p 1, b = (b,, b p 1, and let O + 2p = R + p 1 R + p We then have a map Φ 2p : O + 2p O + 2p x := (a, b x + := (a +, b + where a + i and b + i are given in (511, 512 Iteration of this map, starting at x, produces a sequence x n+1 := Φ 2p (x n of points in O + 2p The rational functions formed with these parameters have integrals that remain constant along this orbit Numerical studies suggest the existence of a number L = L(x R + such that (( ( ( ( ( ( p p p p 1 p 1 p 1 x n,,, ; L, L,, L 1 2 p 1 1 p 1 Thus the integral of the original rational function is π 2 L 6 The sixth degree case We discuss the map T 2p : E 2p E 2p for the case p = 3 The effect of T 6 on the coefficients P 6 = {b, b 1, b 2, a 1, a 2 } is denoted by Φ 6 : O + 6 O+ 6 and is given explicitly by (61 a a 1 + 5a 2 + a 1 a 2 (a 1 + a /3 a 2 a 1 + a (a 1 + a /3 b b + b 1 + b 2 (a 1 + a /3 b 1 b (a b 1 + b 2 (a a 1 + a b + b 2 b 2 (a 1 + a /3 using Theorem 51 The map Φ 6 preserves the integral (62 U 6 (a 1, a 2, b ; b 1, b 2 := b z 4 + b 1 z 2 + b 2 z 6 + a 1 z 4 + a 2 x dz and the convergence of its iterations has been proved in [3], the main result of which is the following theorem Theorem 61 Let x := (a 1, a 2; b, b 1, b 2 R + 5 Define x n+1 := Φ 6 (x n Then U 6 is invariant under Φ 6 Moreover, the sequence {(a n 1, a n 2 } converges to (3, 3 and
11 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS 11 {(b n, b n 1, b n 2 } converges to (L, 2L, L, where the limit L is a function of the initial data x Therefore b z 4 + b 1 z 2 + b 2 z 6 + a 1 z 4 + a 2 z dz = L(x π 2 This iteration is similar to Landen s transformation for elliptic integrals that has been employed in [4] in the efficient calculation of π Numerical data indicate that the convergence of x n is quadratic The proof of convergence is based on the fact that Φ 6 cuts the distance from (a 1, a 2 to (3, 3 by at least half A sequence of algebraic curves The complete characterization of parameters (a 1, a 2 in the first quadrant that yield computable rational functions R(z := b z 4 + b 1 z 2 + b 2 z 6 + a 1 z 4 + a 2 z of degree 6 remains open The polynomial z 6 + a 1 z 4 + a 2 z factors when a 1 = a 2 so the diagonal := {(a 1, a 2 R + R + : a 1 = a 2 } produces computable functions In view of the invariance of the class of computable functions under iterations by Φ 6, the curves X n := Φ ( n 6 (, with n Z, are also computable The curve X 1 has equation (9 + 5a 1 + 5a 2 + a 1 a 2 3 = (a 1 + a (a 1 + a and consists of two branches meeting at the cusp (3, 3 In terms of the coordinates x = a 1 3 and y = a 2 3 the leading order term is T 1 (x, y = 1728(x y 2 This curve is rational and can be parametrized by (63 a 1 (t = t 2 (t 5 t 4 + 2t 3 t 2 + t + 1 a 2 (t = t 3 (t 5 + t 4 t 3 + 2t 2 t + 1 The rationality of X n for n 1 and its significance for the integration algorithm remains open The complexity of these curves increases with n For example, the curve X 2 := Φ ( 2 3 ( is of total degree 9 in x = a 1 3 and y = a 2 3 with leading term T 2 (x, y := (x y 18 [ 163(x 4 + y xy(x 2 + y 2 174x 2 y 2] The diagonal can be replaced by a 2-parameter family of computable curves X(c, d that are produced from the factorization of the sextic with a 1 = c + d and a 2 = cd + 1/d All the images Φ ( n 6 X(c, d with n Z are computable curves The question of whether these are all the computable parameters remains open 7 Examples In this section we present a variety of closed-form evaluations of integrals of rational functions Example 1 The integral z 2 (z 4 + 4z dz = π is computed by Mathematica 3 using (32 in 1 seconds The direct calculation too 1227 seconds and 64 extra seconds to simplify the answer
12 12 GEORGE BOROS AND VICTOR H MOLL Example 2 The integral of any even rational function with denominator a power of an even quartic polynomial can be computed directly by using Corollary 32 For example: z 6 dz (2z 4 + 2z = 11π( ( /2 Example 3 The case n = in (32 deserves special attention: π (71 N,4 (a; m = 2 m+3/2 (a + 1 P m(a m+1/2 where (72 P m (a = 2 2m m = 2 ( 2m 2 m ( m + The polynomial P m (a has been studied in [1] and [2] Example 4 The case n = m in (32 yields N m,4 (a; m = m (a + 1 z 2m dz (z 4 + 2az m+1 = π 2 3m+3/2 (1 + a m+1/2 The change of variable z z converts this integral to which is [8] N m,4 (a; m = 1 2 z m 1/2 dz (z 2 + 2az + 1 m+1, Example 5 A symmetric function of degree 6 The integral I = x 8 (x 6 + 4x 4 + 4x dx = ( 2m m x 8 [(x 2 + 1(x 4 + 3x 2 + 1] 5 dx can be computed by decomposing the integrand into partial fractions as 1 (x (x (x (x (x (x 4 + 3x x 2 (x 4 + 3x (x 4 + 3x x 2 (x 4 + 3x (x 4 + 3x x 2 (x 4 + 3x (x 4 + 3x x 2 (x 4 + 3x Each of these terms is now computable yielding I = (x 4 + 3x x 2 (x 4 + 3x π Example 6 Non-symmetric functions of degree 6 In this case we can use the scheme (61 to produce numerical approximations to the integral For example,
13 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS 13 the evaluation of 45z z z 6 + z 4 + 3z is shown below: n a n 1 a n 2 b n b n 1 b n Thus L and 45x x x 6 + x 4 + 3x dx π 2 = Example 7 Symmetric functions of degree 8 These integrals can be evaluated using Theorem 41 For example: dz ( π (z 8 + 5z z 4 + 5z = ( /2 Example 8 As in the case of degree 6 we can provide numerical approximations to nonsymmetric integrals of degree 8 The iteration (61 is now replaced by a n+1 1 = an 2 (a n 1 + a n 3 + 4a n 1 a n 3 + 1(a n 1 + a n 3 + 8(a n (a n 1 + an 2 + an /2 a2 n+1 = an 1 a n 3 + 6(a n 1 + a n 3 + 2(a n a n 1 + an 2 + a a n+1 3 = b n+1 = a n 1 + a n (a n 1 + an 2 + an /2 b n + b n 1 + b n 2 + b n 3 (a n 1 + an 2 + an /4 b n+1 1 = bn 3 (3a n 1 + a n b n 2 (a n b n 1 (a n b n (3a n 3 + a n (a n 1 + an 2 + an /4 b n+1 2 = bn 3 (a n b n 2 + b n 1 + b n (a n (a n 1 + an 2 + an /4 b n+1 3 = b n + b n 3 (a n 1 + an 2 + an /4 with initial conditions a 1, a 2, a 3, b, b 1, b 2, b 3 Then (73 U 8 (a 1, a 2, a 3, b, b 1, b 2, b 3 := is invariant under these transformations dz b x 6 + b 1 x 4 + b 2 x 2 + b 3 x 8 + a 1 x 6 + a 2 x 4 + a 3 x dx
14 14 GEORGE BOROS AND VICTOR H MOLL Note Numerical calculations show that (a n 1, a n 2, a n 3 (4, 6, 4 and that (b n, b n 1, b n 2, b n 3 (1, 3, 3, 1 L for some L depending upon the initial conditions Example 9 A symmetric function of degree 12 We use Theorem 21 to evaluate as (74 I := z 18 dz (z z z 8 + 4z z z I = 25π( Here p = 3, n = 9, and m = 2, so n > (m + 1p 1 and we need to apply the transformation z 1/z to reduce the value of n Indeed, we have I = and Theorem 21 now yields z 16 dz (z z z 8 + 4z z z , I = 2 17 z 16 dz 13172(1 + z 2 3 (1 + 4z 2 + z 4 3 The new integrand is expanded in partial fractions in the variable t = z 2 to produce (74 Example 1 We use Theorem 21 to evaluate where I := z 1 dz Q 2 (z Q(z = z 2 + 6z z 16 24z z z z 8 24z z 4 + 6z The factorization with Q(z = (1 + z 2 2 T (zt ( z T (z = z 8 2z 7 + 4z z 5 + 6z 4 14z 3 + 4z 2 + 2z + 1 leads to a partial fraction expansion containing the term 72 51z z z z 4 43z z 6 55z (1 2z + 4z z 3 + 6z 4 14z 5 + 4z 6 + 2z 7 + z 8 2, which we were unable to integrate; furthemore, the roots of T (z = cannot be evaluated by radicals The procedure described in Theorem 21, however, shows that I = z 1 (4 + z 2 (z z z dz (z (z 8 + 3z 6 + 8z 4 + 3z ,
15 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS 15 the integrand of which can be expanded to yield 9 I = dz (z dz z z z z (z 8 + 3z 6 + 8z 4 + 3z dz 91z z z (z 8 + 3z 6 + 8z 4 + 3z dz Every piece is now computable, with the final result I = ( π Example 11 The symmetric functions of degree 16 have denominator D 4 (d 1, d 2, d 3, d 4 ; z = z 16 + d 4 z 14 + d 3 z 12 + d 2 z 1 + 2d 1 z 8 + d 2 z 6 + d 3 z 4 + d 4 z 2 + 1, the integral of which is computed in terms of E 4 (d 1, d 2, d 3, d 4 ; z = (1 + d 1 + d 2 + d 3 + d 4 z 8 + 2(16 + d 2 + 4d 3 + 9d 4 z 6 This new integral is symmetric provided [ ] [ ] d1 15 (75 = + d (2 + d 3 + 6d 4 z (8 + d 4 z [ 3 4] d 3 + [ ] 8 d 7 4 Introduce the new parameters ( ( 8 8 e j = d j for 2 j 4 and e 1 = d j 2 4 Then (75 yields (76 [ e1 e 2 ] = [ ] 3 e Thus, if the original denominator has the form [ ] 8 e 7 4 D 4 (z = (z d 4 (z 14 + z 2 + d 3 (z 12 + z 4 + (112 4d 3 + 7d 4 (z 1 + z 6 the integral +2(15 + 3d 3 8d 4 z 8, P (z ( D 4 (z m+1 dz is reduced to an integral with symmetric denominator of degree 8 and these are computable We can thus evaluate a 2-parameter family of symmetric integrals of degree 16 For example, tae d 3 = d 4 = 1 to obtain (77 R 1 (z = The main theorem yields (78 z 4 (z 16 + z 14 + z z 1 + 2z z 6 + z 4 + z R 2 (z = 124z z z z 1 + 6z 12 + z (16z z z z
16 16 GEORGE BOROS AND VICTOR H MOLL so that R 1 (z dz = Letting f[n] := N n,8 [1, n, 27/32, 9/4] we obtain R 2 (z dz R 1 (z dz = 2 15 (f[] + 6f[1] f[2] + 496f[3] and conclude that z 4 dz (z 16 + z 14 + z z 1 + 2z z 6 + z 4 + z = = ( π Example 12 We classify the symmetric denominators of degree 32 that yield computable integrals These functions depend on 8 parameters (79 D 8 (d 1,, d 8 ; z = 8 d 9 (z 2 + z 2(16 = and the main theorem expresses the integral in terms of E 8 The conditions for E 8 to be symmetric yield and the symmetric E 8 is d 1 = d d 6 312d d 8 d 2 = d 5 11d d d 8 d 3 = d d 6 329d 7 224d 8 d 4 = 496 8d 5 19d 6 + 8d d 8 E 8 (d 5, d 6, d 7, d 8 ; z = 32768(1 + z 16 + ( d 8 (z 2 + z 14 + ( d d 8 (z 4 + z 12 + ( d d d 8 (z 6 + z 1 + d 5 d 6 +( d d d d 8 z 8 The symmetry of E 8 now determines d 5, d 6 in terms of d 7, d 8 and we obtain d d d 3 d 4 = d ( d The function (z is a symmetric polynomial of degree 32 and yields a particular solution to (71 As before let e j = d j ( 16 9 j for 2 j 8 and e 1 = d ( 16 8
17 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS 17 Then (711 e 1 e 2 e 3 e 4 e 5 e 6 = e and as in the case of degree 16 we can compute a 2-parameter family of symmetric integrals of degree 32 e 8 Appendix A Two binomial sums The closed-form evaluation of sums involving binomials coefficients can be obtained by traditional analytical techniques or by using the powerful WZ-method as described in [16] We discuss two sums used to simplify expressions in later sections, presenting one proof in each style Lemma A1 Let, N be positive integers with N Then N ( ( ( 2N + 1 N j 2N (A1 = 4 N j= Proof Multiply the left hand side of (A1 by x and sum over to produce N N ( ( 2N + 1 N j N ( N j 2N + 1 ( N j x = x = j= = j= = N ( 2N + 1 (x + 1 N j j= = (1 + x + 1 2N+1 (1 x + 1 2N+1 2 x + 1 N ( 2N = 4 N x In order to justify the last step we start with the well nown result ( i y = ( 2 + i (A2 y 1 4y 2y (see WILF [19], page 54 Letting x = 4y and i = 2N + 1 in (A2 gives (1 x + 1 2N+1 ( 2N = ( (N+1+ x 2N+1+ ; x + 1 = similarly x = 4y and i = 2N 1 yields (1 + x + 1 2N+1 ( 2N = ( (N+1+ x 2N+1+ x + 1 =2N+1 =
18 18 GEORGE BOROS AND VICTOR H MOLL Thus (1 + x + 1 2N+1 (1 x + 1 2N+1 2 x + 1 = = 2N ( (2N ( 1 4 N x = N ( 2N 4 N x = Lemma A2 Let N N Then N ( 2N + 1 (1 + z 2 N j = j= N ( N + j 4 j z 2(N j j= Proof The coefficient of z 2 on the left hand side is N ( ( 2N + 1 N j, j= and the corresponding coefficient on the right hand side is ( 2N 4 N The result then follows from Lemma A1 Lemma A3 Let, N N with N Then ( ( 2N N j = 22 1 N N j= F (; j = ( + N 1 N Proof This lemma could be proven in the same style as Lemma A1 Instead we use the WZ-method as explained in [16] Indeed, let ( ( 2N N j N N2 2 1( +N 1 N and define, with the pacage EKHAD, the function G(; j = F (; j, j( 1 2(N + ( j + 1 Then F (; j F ( + 1; j = G(; j + 1 G(; j, and summing over j we see that the sum of F (; j over j is independent of The case = N produces 1 as the common value Lemma A4 Let p N, d 1, d 2,, d p be parameters, and define d p+1 := 1 Then p p ( 2p 2 (A3 d p+1 z 2 (1 + z 2 p j = p+1 j=1 d j = z 2p + j= p 2 2i 1 z 2(p i i=1 p+1 i j=1 j + i 1 i ( j + 2i 2 j 1 d j+i
19 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS 19 Proof For fixed i p 1 the coefficient of z 2i on the right hand side of (A3 is [RHS] (2i = 22(p i 1 p+1 ( r + p i 2 (r 1 d r, p i r p + i 1 r=p i+1 and for i = p we have [RHS] (2p = 1 + p j=1 d j Similarly, for the left hand side of (A3, [LHS] (2i = p+1 r=p+1 i d r p i ( 2r 2 j= ( r j 1 r j 1 + i It is easy to chec that the coefficients of z 2p match It suffices to show that for each i such that i p 1 and for each r such that p + 1 i r p + 1 we have p i ( 2r 2 j= ( r 1 j r 1 p + i = 22(p i 1 (r 1 p i This follows from Lemma A1 with = p i and N = r 1 ( r + p i 2 r p + i 1 The suggestions of the referrees and the editor are gratefully acnowledged References [1] BOROS, G - MOLL, V: A criterion for unimodality Elec Jour of Combinatorics, 6 (1999, #R1 [2] BOROS, G - MOLL, V: An integral hidden in Gradshteyn and Ryzhi Jour Comp Appl Math 16, , 1999 [3] BOROS, G - MOLL, V: A rational Landen transformation Contemporary Mathematics 251, 83-91, 2 [4] BORWEIN, J - BORWEIN, P: Pi and the AGM Canadian Mathematical Society Wiley- Interscience Publication [5] BRONSTEIN, M: Symbolic Integration I Transcendental functions Algorithms and Computation in Mathematics, 1 Springer-Verlag, 1997 [6] GEDDES, K - CZAPOR, SR - LABAHN, G: Algorithms for Computer Algebra Kluwer, Dordrecht The Netherlands, 1992 [7] GAUSS, KF: Arithmetische Geometrisches Mittel, 1799 In Were, 3, Konigliche Gesellschaft der Wissenschaft, Gottingen Reprinted by Olms, Hildescheim, 1981 [8] GRADSHTEYN, IS - RYZHIK, IM: Table of Integrals, Series and Products Fifth Edition, ed Alan Jeffrey Academic Press, 1994 [9] HARDY, GH: The Integration of Functions of a Single Variable Cambridge Tracts in Mathematics and Mathematical Physics, 2, Second Edition, Cambridge University Press, 1958 [1] HERMITE, C: Sur l integration des fractions rationelles Nouvelles Annales de Mathematiques ( 2 eme serie 11, , 1872 [11] HOROWITZ, E: Algorithms for partial fraction decomposition and rational function integration Proc of SYMSAM 71, ACM Press, , 1971 [12] LAZARD, D - RIOBOO, R: Integration of Rational Functions: Rational Computation of the Logarithmic Part Journal of Symbolic Computation 9, , 199 [13] MCKEAN, H - MOLL, V: Elliptic Curves: Function Theory, Geometry, Arithmetic Cambridge University Press, 1997 [14] NEWMAN, D: A simplified version of the fast algorithm of Brent and Salamin Math Comp 44, 27-21, 1985 [15] OSTROGRADSKY, MW: De l integration des fractions rationelles Bulletin de la Classe Physico-Mathematiques de l Academie Imperieriale des Sciences de St Petersbourgh, IV, ,
20 2 GEORGE BOROS AND VICTOR H MOLL [16] PETKOVSEK, M - WILF, HS - ZEILBERGER, D: A=B A K Peters, Wellesley, Massachusetts 1996 [17] ROTHSTEIN, M: A new algorithm for the integration of Exponential and Logarithmic Functions, Proc of the 1977 MACSYMA Users Conference, NASA Pub, CP-212, [18] TRAGER, BM: Algebraic factoring and rational function integration Proc SYMSAC 76, [19] WILF, HS: generatingfunctionology Academic Press, 199 Department of Mathematics, Xavier University, New Orleans, Louisiana address: gboros@xulamathedu Department of Mathematics, Tulane University, New Orleans, LA address: vhm@mathtulaneedu
LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS. 1. Introduction We consider the space of even rational functions of degree 2p
MATHEMATICS OF COMPUTATION Volume 71, Number 238, Pages 649 668 S 25-5718(1)1347-3 Article electronically published on November 9, 21 LANDEN TRANSFORMATIONS AND THE INTEGRATION OF RATIONAL FUNCTIONS GEORGE
More informationA GEOMETRIC VIEW OF RATIONAL LANDEN TRANSFORMATIONS
Bull. London Math. Soc. 35 (3 93 3 C 3 London Mathematical Society DOI:./S4693393 A GEOMETRIC VIEW OF RATIONAL LANDEN TRANSFORMATIONS JOHN HUBBARD and VICTOR MOLL Abstract In this paper, a geometric interpretation
More informationA HYPERGEOMETRIC INEQUALITY
A HYPERGEOMETRIC INEQUALITY ATUL DIXIT, VICTOR H. MOLL, AND VERONIKA PILLWEIN Abstract. A sequence of coefficients that appeared in the evaluation of a rational integral has been shown to be unimodal.
More informationA transformation of rational functions
Elem. Math. 58 (23) 73 83 13-618/3/273-11 c Birkhäuser Verlag, Basel, 23 Elemente der Mathematik A transformation of rational functions George Boros, Michael Joyce and Victor H. Moll George Boros was born
More informationSums of arctangents and some formulas of Ramanujan
SCIENTIA Series A: Mathematical Sciences, Vol. 11 (2005), 13 24 Universidad Técnica Federico Santa María Valparaíso, Chile ISSN 0716-8446 c Universidad Técnica Federico Santa María 2005 Sums of arctangents
More informationA PROBABILISTIC PROOF OF WALLIS S FORMULA FOR π. ( 1) n 2n + 1. The proof uses the fact that the derivative of arctan x is 1/(1 + x 2 ), so π/4 =
A PROBABILISTIC PROOF OF WALLIS S FORMULA FOR π STEVEN J. MILLER There are many beautiful formulas for π see for example [4]). The purpose of this note is to introduce an alternate derivation of Wallis
More informationarxiv: v1 [math.ca] 17 Feb 2017
GENERALIZED STIELTJES CONSTANTS AND INTEGRALS INVOLVING THE LOG-LOG FUNCTION: KUMMER S THEOREM IN ACTION OMRAN KOUBA arxiv:7549v [mathca] 7 Feb 7 Abstract In this note, we recall Kummer s Fourier series
More informationTewodros Amdeberhan, Dante Manna and Victor H. Moll Department of Mathematics, Tulane University New Orleans, LA 70118
The -adic valuation of Stirling numbers Tewodros Amdeberhan, Dante Manna and Victor H. Moll Department of Mathematics, Tulane University New Orleans, LA 7011 Abstract We analyze properties of the -adic
More informationRelations and Functions (for Math 026 review)
Section 3.1 Relations and Functions (for Math 026 review) Objective 1: Understanding the s of Relations and Functions Relation A relation is a correspondence between two sets A and B such that each element
More informationA Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)!
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 8 (2005), Article 05.2.3 A Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)! Ira M. Gessel 1 and Guoce Xin Department of Mathematics Brandeis
More informationAn Arithmetic-Geometric Mean Inequality
EXPOSITIONES MATHEMATICAE Mathematical Note Expo. Math. 19 (2001):273-279 Urban & Fischer Verlag www. u rbanfischer.de/jou rnals/expomath An Arithmetic-Geometric Mean Inequality Paul Bracken Centre de
More informationCERTAIN INTEGRALS ARISING FROM RAMANUJAN S NOTEBOOKS
CERTAIN INTEGRALS ARISING FROM RAMANUJAN S NOTEBOOKS BRUCE C. BERNDT AND ARMIN STRAUB Abstract. In his third noteboo, Ramanujan claims that log x dx + π 2 x 2 dx. + 1 In a following cryptic line, which
More informationA COMPUTER PROOF OF A SERIES EVALUATION IN TERMS OF HARMONIC NUMBERS
A COMPUTER PROOF OF A SERIES EVALUATION IN TERMS OF HARMONIC NUMBERS RUSSELL LYONS, PETER PAULE, AND AXEL RIESE Abstract. A fruitful interaction between a new randomized WZ procedure and other computer
More informationThe integrals in Gradshteyn and Ryzhik. Part 10: The digamma function
SCIENTIA Series A: Mathematical Sciences, Vol 7 29, 45 66 Universidad Técnica Federico Santa María Valparaíso, Chile ISSN 76-8446 c Universidad Técnica Federico Santa María 29 The integrals in Gradshteyn
More informationA REMARK ON THE BOROS-MOLL SEQUENCE
#A49 INTEGERS 11 (2011) A REMARK ON THE BOROS-MOLL SEQUENCE J.-P. Allouche CNRS, Institut de Math., Équipe Combinatoire et Optimisation Université Pierre et Marie Curie, Paris, France allouche@math.jussieu.fr
More informationMath Homework 2
Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is
More informationLinear Recurrence Relations for Sums of Products of Two Terms
Linear Recurrence Relations for Sums of Products of Two Terms Yan-Ping Mu College of Science, Tianjin University of Technology Tianjin 300384, P.R. China yanping.mu@gmail.com Submitted: Dec 27, 2010; Accepted:
More informationANALOGUES OF THE TRIPLE PRODUCT IDENTITY, LEBESGUE S IDENTITY AND EULER S PENTAGONAL NUMBER THEOREM
q-hypergeometric PROOFS OF POLYNOMIAL ANALOGUES OF THE TRIPLE PRODUCT IDENTITY, LEBESGUE S IDENTITY AND EULER S PENTAGONAL NUMBER THEOREM S OLE WARNAAR Abstract We present alternative, q-hypergeometric
More information1. Introduction Interest in this project began with curiosity about the Laplace transform of the Digamma function, e as ψ(s + 1)ds,
ON THE LAPLACE TRANSFORM OF THE PSI FUNCTION M. LAWRENCE GLASSER AND DANTE MANNA Abstract. Guided by numerical experimentation, we have been able to prove that Z 8 / x x + ln dx = γ + ln) [cosx)] and to
More informationCombinatorial Analysis of the Geometric Series
Combinatorial Analysis of the Geometric Series David P. Little April 7, 205 www.math.psu.edu/dlittle Analytic Convergence of a Series The series converges analytically if and only if the sequence of partial
More informationA Computer Proof of Moll s Log-Concavity Conjecture
A Computer Proof of Moll s Log-Concavity Conjecture Manuel Kauers and Peter Paule Research Institute for Symbolic Computation (RISC-Linz) Johannes Kepler University Linz, Austria, Europe May, 29, 2006
More informationarxiv:math/ v1 [math.nt] 9 Aug 2004
arxiv:math/0408107v1 [math.nt] 9 Aug 2004 ELEMENTARY RESULTS ON THE BINARY QUADRATIC FORM a 2 + ab + b 2 UMESH P. NAIR Abstract. This paper examines with elementary proofs some interesting properties of
More informationBP -HOMOLOGY AND AN IMPLICATION FOR SYMMETRIC POLYNOMIALS. 1. Introduction and results
BP -HOMOLOGY AND AN IMPLICATION FOR SYMMETRIC POLYNOMIALS DONALD M. DAVIS Abstract. We determine the BP -module structure, mod higher filtration, of the main part of the BP -homology of elementary 2- groups.
More informationAn Approach to Hensel s Lemma
Irish Math. Soc. Bulletin 47 (2001), 15 21 15 An Approach to Hensel s Lemma gary mcguire Abstract. Hensel s Lemma is an important tool in many ways. One application is in factoring polynomials over Z.
More informationExplicit solution of a class of quartic Thue equations
ACTA ARITHMETICA LXIV.3 (1993) Explicit solution of a class of quartic Thue equations by Nikos Tzanakis (Iraklion) 1. Introduction. In this paper we deal with the efficient solution of a certain interesting
More informationLEGENDRE S THEOREM, LEGRANGE S DESCENT
LEGENDRE S THEOREM, LEGRANGE S DESCENT SUPPLEMENT FOR MATH 370: NUMBER THEORY Abstract. Legendre gave simple necessary and sufficient conditions for the solvablility of the diophantine equation ax 2 +
More informationMultiplicative Order of Gauss Periods
Multiplicative Order of Gauss Periods Omran Ahmadi Department of Electrical and Computer Engineering University of Toronto Toronto, Ontario, M5S 3G4, Canada oahmadid@comm.utoronto.ca Igor E. Shparlinski
More informationOn the Length of Lemniscates
On the Length of Lemniscates Alexandre Eremenko & Walter Hayman For a monic polynomial p of degree d, we write E(p) := {z : p(z) =1}. A conjecture of Erdős, Herzog and Piranian [4], repeated by Erdős in
More informationThe Circle Method. Basic ideas
The Circle Method Basic ideas 1 The method Some of the most famous problems in Number Theory are additive problems (Fermat s last theorem, Goldbach conjecture...). It is just asking whether a number can
More informationarxiv: v3 [math.ca] 14 Oct 2015
Symmetry, Integrability and Geometry: Methods and Applications Certain Integrals Arising from Ramanujan s Noteboos SIGMA 11 215), 83, 11 pages Bruce C. BERNDT and Armin STRAUB arxiv:159.886v3 [math.ca]
More informationSome generalizations of a supercongruence of van Hamme
Some generalizations of a supercongruence of van Hamme Victor J. W. Guo School of Mathematical Sciences, Huaiyin Normal University, Huai an, Jiangsu 3300, People s Republic of China jwguo@hytc.edu.cn Abstract.
More informationA COMBINATORIAL PROOF OF A RESULT FROM NUMBER THEORY
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 4 (004), #A09 A COMBINATORIAL PROOF OF A RESULT FROM NUMBER THEORY Shaun Cooper Institute of Information and Mathematical Sciences, Massey University
More informationMath 200 University of Connecticut
IRRATIONALITY OF π AND e KEITH CONRAD Math 2 University of Connecticut Date: Aug. 3, 25. Contents. Introduction 2. Irrationality of π 2 3. Irrationality of e 3 4. General Ideas 4 5. Irrationality of rational
More informationIntroduction to Arithmetic Geometry Fall 2013 Lecture #2 09/10/2013
18.78 Introduction to Arithmetic Geometry Fall 013 Lecture # 09/10/013.1 Plane conics A conic is a plane projective curve of degree. Such a curve has the form C/k : ax + by + cz + dxy + exz + fyz with
More informationCOMBINATORICS OF GENERALIZED q-euler NUMBERS. 1. Introduction The Euler numbers E n are the integers defined by E n x n = sec x + tan x. (1.1) n!
COMBINATORICS OF GENERALIZED q-euler NUMBERS TIM HUBER AND AE JA YEE Abstract New enumerating functions for the Euler numbers are considered Several of the relevant generating functions appear in connection
More informationPART II: Research Proposal Algorithms for the Simplification of Algebraic Formulae
Form 101 Part II 6 Monagan, 195283 PART II: Research Proposal Algorithms for the Simplification of Algebraic Formulae 1 Research Area Computer algebra (CA) or symbolic computation, as my field is known
More information#A31 INTEGERS 18 (2018) A NOTE ON FINITE SUMS OF PRODUCTS OF BERNSTEIN BASIS POLYNOMIALS AND HYPERGEOMETRIC POLYNOMIALS
#A31 INTEGERS 18 (2018) A NOTE ON FINITE SUMS OF PRODUCTS OF BERNSTEIN BASIS POLYNOMIALS AND HYPERGEOMETRIC POLYNOMIALS Steven P. Clar Department of Finance, University of North Carolina at Charlotte,
More informationPOCHHAMMER SYMBOL WITH NEGATIVE INDICES. A NEW RULE FOR THE METHOD OF BRACKETS.
POCHHAMMER SYMBOL WITH NEGATIVE INDICES A NEW RULE FOR THE METHOD OF BRACKETS IVAN GONZALEZ, LIN JIU, AND VICTOR H MOLL Abstract The method of brackets is a method of integration based upon a small number
More informationClassical transcendental curves
Classical transcendental curves Reinhard Schultz May, 2008 In his writings on coordinate geometry, Descartes emphasized that he was only willing to work with curves that could be defined by algebraic equations.
More informationSome theorems on the explicit evaluations of singular moduli 1
General Mathematics Vol 17 No 1 009) 71 87 Some theorems on the explicit evaluations of singular moduli 1 K Sushan Bairy Abstract At scattered places in his notebooks Ramanujan recorded some theorems for
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.nt] 9 Sep 2017
arxiv:179.2954v1 [math.nt] 9 Sep 217 ON THE FACTORIZATION OF x 2 +D GENERALIZED RAMANUJAN-NAGELL EQUATION WITH HUGE SOLUTION) AMIR GHADERMARZI Abstract. Let D be a positive nonsquare integer such that
More informationIntegrals evaluated in terms of Catalan s constant. Graham Jameson and Nick Lord (Math. Gazette, March 2017)
Integrals evaluated in terms of Catalan s constant Graham Jameson and Nick Lord (Math. Gazette, March 7) Catalan s constant, named after E. C. Catalan (84 894) and usually denoted by G, is defined by G
More informationCORRESPONDENCE BETWEEN ELLIPTIC CURVES IN EDWARDS-BERNSTEIN AND WEIERSTRASS FORMS
CORRESPONDENCE BETWEEN ELLIPTIC CURVES IN EDWARDS-BERNSTEIN AND WEIERSTRASS FORMS DEPARTMENT OF MATHEMATICS AND STATISTICS UNIVERSITY OF OTTAWA SUPERVISOR: PROFESSOR MONICA NEVINS STUDENT: DANG NGUYEN
More informationCLOSED FORM CONTINUED FRACTION EXPANSIONS OF SPECIAL QUADRATIC IRRATIONALS
CLOSED FORM CONTINUED FRACTION EXPANSIONS OF SPECIAL QUADRATIC IRRATIONALS DANIEL FISHMAN AND STEVEN J. MILLER ABSTRACT. We derive closed form expressions for the continued fractions of powers of certain
More informationLinear Ordinary Differential Equations Satisfied by Modular Forms
MM Research Preprints, 0 7 KLMM, AMSS, Academia Sinica Vol. 5, December 006 Linear Ordinary Differential Equations Satisfied by Modular Forms Xiaolong Ji Department of Foundation Courses Jiyuan Vocational
More informationAlgebraic. techniques1
techniques Algebraic An electrician, a bank worker, a plumber and so on all have tools of their trade. Without these tools, and a good working knowledge of how to use them, it would be impossible for them
More informationExtend Fermats Small Theorem to r p 1 mod p 3 for divisors r of p ± 1
Extend Fermats Small Theorem to r p 1 mod p 3 for divisors r of p ± 1 Nico F. Benschop AmSpade Research, The Netherlands Abstract By (p ± 1) p p 2 ± 1 mod p 3 and by the lattice structure of Z(.) mod q
More informationON YANG MEANS III. Edward Neuman
BULLETIN OF THE INTERNATIONAL MATHEMATICAL VIRTUAL INSTITUTE ISSN (p) 2303-4874, ISSN (o) 2303-4955 www.imvibl.org /JOURNALS / BULLETIN Vol. 8(2018), 113-122 DOI: 10.7251/BIMVI1801113N Former BULLETIN
More informationDescribing Resolvent Sextics of Quintic Equations. By Frank D. Grosshans
Describing Resolvent Sextics of Quintic Equations By Frank D. Grosshans L. E. Dickson, Bulletin A.M.S. 31, 1925, 515-523 RESOLVENT SEXTICS OF QUINTIC EQUATIONS BY L. E. DICKSON 1. Introduction. The object
More informationCalculus. Integration (III)
Calculus Integration (III) Outline 1 Other Techniques of Integration Partial Fractions Integrals Involving Powers of Trigonometric Functions Trigonometric Substitution 2 Using Tables of Integrals Integration
More informationMathematics 324 Riemann Zeta Function August 5, 2005
Mathematics 324 Riemann Zeta Function August 5, 25 In this note we give an introduction to the Riemann zeta function, which connects the ideas of real analysis with the arithmetic of the integers. Define
More informationA Note on Indefinite Ternary Quadratic Forms Representing All Odd Integers. Key Words: Quadratic Forms, indefinite ternary quadratic.
Bol. Soc. Paran. Mat. (3s.) v. 23 1-2 (2005): 85 92. c SPM ISNN-00378712 A Note on Indefinite Ternary Quadratic Forms Representing All Odd Integers Jean Bureau and Jorge Morales abstract: In this paper
More informationPILLAI S CONJECTURE REVISITED
PILLAI S COJECTURE REVISITED MICHAEL A. BEETT Abstract. We prove a generalization of an old conjecture of Pillai now a theorem of Stroeker and Tijdeman) to the effect that the Diophantine equation 3 x
More informationSome properties and applications of a new arithmetic function in analytic number theory
NNTDM 17 (2011), 3, 38-48 Some properties applications of a new arithmetic function in analytic number theory Ramesh Kumar Muthumalai Department of Mathematics, D.G. Vaishnav College Chennai-600 106, Tamil
More informationPólya s Random Walk Theorem
Pólya s Random Walk Theorem Jonathan ovak Abstract. This note presents a proof of Pólya s random walk theorem using classical methods from special function theory and asymptotic analysis. 1. ITRODUCTIO.
More informationIntegration Using Tables and Summary of Techniques
Integration Using Tables and Summary of Techniques Philippe B. Laval KSU Today Philippe B. Laval (KSU) Summary Today 1 / 13 Introduction We wrap up integration techniques by discussing the following topics:
More informationA Symbolic Operator Approach to Power Series Transformation-Expansion Formulas
A Symbolic Operator Approach to Power Series Transformation-Expansion Formulas Tian- Xiao He Department of Mathematics and Computer Science Illinois Wesleyan University Bloomington, IL 61702-2900, USA
More informationIntroduction to Techniques for Counting
Introduction to Techniques for Counting A generating function is a device somewhat similar to a bag. Instead of carrying many little objects detachedly, which could be embarrassing, we put them all in
More informationMathematics 136 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 19 and 21, 2016
Mathematics 36 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 9 and 2, 206 Every rational function (quotient of polynomials) can be written as a polynomial
More informationQuadratic Transformations of Hypergeometric Function and Series with Harmonic Numbers
Quadratic Transformations of Hypergeometric Function and Series with Harmonic Numbers Martin Nicholson In this brief note, we show how to apply Kummer s and other quadratic transformation formulas for
More informationNew series expansions of the Gauss hypergeometric function
New series expansions of the Gauss hypergeometric function José L. López and Nico. M. Temme 2 Departamento de Matemática e Informática, Universidad Pública de Navarra, 36-Pamplona, Spain. e-mail: jl.lopez@unavarra.es.
More informationPólya s Random Walk Theorem arxiv: v2 [math.pr] 18 Apr 2013
Pólya s Random Wal Theorem arxiv:131.3916v2 [math.pr] 18 Apr 213 Jonathan ova Abstract This note presents a proof of Pólya s random wal theorem using classical methods from special function theory and
More informationAdditive Decompositions in Primitive Extensions
Additive Decompositions in Primitive Extensions KLMM, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China School of Mathematical Sciences, University of Chinese
More informationThe problem of integration in finite terms
The problem of integration in finite terms Symbolic Integration Lecture 14 1. The Freshman Calculus Problem 2. The Artificial Intelligence Problem 3. The Algebraic (Rational Function) Problem 4. The Decision
More informationq-series Michael Gri for the partition function, he developed the basic idea of the q-exponential. From
q-series Michael Gri th History and q-integers The idea of q-series has existed since at least Euler. In constructing the generating function for the partition function, he developed the basic idea of
More informationTHE BAILEY TRANSFORM AND FALSE THETA FUNCTIONS
THE BAILEY TRANSFORM AND FALSE THETA FUNCTIONS GEORGE E ANDREWS 1 AND S OLE WARNAAR 2 Abstract An empirical exploration of five of Ramanujan s intriguing false theta function identities leads to unexpected
More informationYi Wang Department of Applied Mathematics, Dalian University of Technology, Dalian , China (Submitted June 2002)
SELF-INVERSE SEQUENCES RELATED TO A BINOMIAL INVERSE PAIR Yi Wang Department of Applied Mathematics, Dalian University of Technology, Dalian 116024, China (Submitted June 2002) 1 INTRODUCTION Pairs of
More informationMS 3011 Exercises. December 11, 2013
MS 3011 Exercises December 11, 2013 The exercises are divided into (A) easy (B) medium and (C) hard. If you are particularly interested I also have some projects at the end which will deepen your understanding
More informationELLIPSES AND ELLIPTIC CURVES. M. Ram Murty Queen s University
ELLIPSES AND ELLIPTIC CURVES M. Ram Murty Queen s University Planetary orbits are elliptical What is an ellipse? x 2 a 2 + y2 b 2 = 1 An ellipse has two foci From: gomath.com/geometry/ellipse.php Metric
More informationGEOMETRIC CONSTRUCTIONS AND ALGEBRAIC FIELD EXTENSIONS
GEOMETRIC CONSTRUCTIONS AND ALGEBRAIC FIELD EXTENSIONS JENNY WANG Abstract. In this paper, we study field extensions obtained by polynomial rings and maximal ideals in order to determine whether solutions
More informationOn the length of lemniscates
On the length of lemniscates Alexandre Eremenko and Walter Hayman Abstract We show that for a polynomial p(z) =z d +... the length of the level set E(p) :={z : p(z) =1} is at most 9.173 d, which improves
More informationAlg II Syllabus (First Semester)
Alg II Syllabus (First Semester) Unit 1: Solving linear equations and inequalities Lesson 01: Solving linear equations Lesson 02: Solving linear inequalities (See Calculator Appendix A and associated video.)
More informationNewton, Fermat, and Exactly Realizable Sequences
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 8 (2005), Article 05.1.2 Newton, Fermat, and Exactly Realizable Sequences Bau-Sen Du Institute of Mathematics Academia Sinica Taipei 115 TAIWAN mabsdu@sinica.edu.tw
More informationHSC Marking Feedback 2017
HSC Marking Feedback 017 Mathematics Extension 1 Written Examination Question 11 Part (a) The large majority of responses showed correct substitution into the formula x = kx +lx 1 k+l given on the Reference
More informationMCMULLEN S ROOT-FINDING ALGORITHM FOR CUBIC POLYNOMIALS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 130, Number 9, Pages 2583 2592 S 0002-9939(02)06659-5 Article electronically published on April 22, 2002 MCMULLEN S ROOT-FINDING ALGORITHM FOR CUBIC
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 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 information9709 MATHEMATICS. 9709/33 Paper 3, maximum raw mark 75
CAMBRIDGE INTERNATIONAL EXAMINATIONS GCE Advanced Level MARK SCHEME for the May/June 04 series 9709 MATHEMATICS 9709/ Paper, maximum raw mark 75 This mark scheme is published as an aid to teachers and
More informationJUST THE MATHS UNIT NUMBER 1.3. ALGEBRA 3 (Indices and radicals (or surds)) A.J.Hobson
JUST THE MATHS UNIT NUMBER 1 ALGEBRA (Indices and radicals (or surds)) by AJHobson 11 Indices 12 Radicals (or Surds) 1 Exercises 14 Answers to exercises UNIT 1 - ALGEBRA - INDICES AND RADICALS (or Surds)
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationarxiv: v6 [math.nt] 12 Sep 2017
Counterexamples to the conjectured transcendence of /α) k, its closed-form summation and extensions to polygamma functions and zeta series arxiv:09.44v6 [math.nt] Sep 07 F. M. S. Lima Institute of Physics,
More informationCore A-level mathematics reproduced from the QCA s Subject criteria for Mathematics document
Core A-level mathematics reproduced from the QCA s Subject criteria for Mathematics document Background knowledge: (a) The arithmetic of integers (including HCFs and LCMs), of fractions, and of real numbers.
More informationA Note on Extended Binomial Coefficients
A Note on Extended Binomial Coefficients Thorsten Neuschel arxiv:407.749v [math.co] 8 Jul 04 Abstract We study the distribution of the extended binomial coefficients by deriving a complete asymptotic expansion
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 informationApplicable Analysis and Discrete Mathematics available online at ABEL S METHOD ON SUMMATION BY PARTS.
Applicable Analysis and Discrete Mathematics available online at http://pefmathetfrs Appl Anal Discrete Math 4 010), 54 65 doi:1098/aadm1000006c ABEL S METHOD ON SUMMATION BY PARTS AND BALANCED -SERIES
More informationHypergeometric series and the Riemann zeta function
ACTA ARITHMETICA LXXXII.2 (997) Hypergeometric series and the Riemann zeta function by Wenchang Chu (Roma) For infinite series related to the Riemann zeta function, De Doelder [4] established numerous
More informationBinomial coefficients and k-regular sequences
Binomial coefficients and k-regular sequences Eric Rowland Hofstra University New York Combinatorics Seminar CUNY Graduate Center, 2017 12 22 Eric Rowland Binomial coefficients and k-regular sequences
More informationOn Recurrences for Ising Integrals
On Recurrences for Ising Integrals Flavia Stan Research Institute for Symbolic Computation (RISC-Linz) Johannes Kepler University Linz, Austria December 7, 007 Abstract We use WZ-summation methods to compute
More informationA CLOSED-FORM SOLUTION MIGHT BE GIVEN BY A TREE. VALUATIONS OF QUADRATIC POLYNOMIALS
A CLOSED-FORM SOLUTION MIGHT BE GIVEN BY A TREE. VALUATIONS OF QUADRATIC POLYNOMIALS ALYSSA N. BYRNES, JULIE FINK, GARY LAVIGNE, ISABELLE NOGUES, SENTHIL RAJASEKARAN, AMBER YUAN, LEYDA ALMODOVAR, XIAO
More informationAlgebra II. A2.1.1 Recognize and graph various types of functions, including polynomial, rational, and algebraic functions.
Standard 1: Relations and Functions Students graph relations and functions and find zeros. They use function notation and combine functions by composition. They interpret functions in given situations.
More informationPYTHAGOREAN TRIPLES KEITH CONRAD
PYTHAGOREAN TRIPLES KEITH CONRAD 1. Introduction A Pythagorean triple is a triple of positive integers (a, b, c) where a + b = c. Examples include (3, 4, 5), (5, 1, 13), and (8, 15, 17). Below is an ancient
More information2018 Mathematics. Advanced Higher. Finalised Marking Instructions
National Qualifications 08 08 Mathematics Advanced Higher Finalised Marking Instructions Scottish Qualifications Authority 08 The information in this publication may be reproduced to support SQA qualifications
More informationHardy martingales and Jensen s Inequality
Hardy martingales and Jensen s Inequality Nakhlé H. Asmar and Stephen J. Montgomery Smith Department of Mathematics University of Missouri Columbia Columbia, Missouri 65211 U. S. A. Abstract Hardy martingales
More informationarxiv: v1 [math.nt] 4 Jun 2018
On the orbit of a post-critically finite polynomial of the form x d +c arxiv:1806.01208v1 [math.nt] 4 Jun 2018 Vefa Goksel Department of Mathematics University of Wisconsin Madison, WI 53706, USA E-mail:
More informationTHE CONGRUENT NUMBER PROBLEM
THE CONGRUENT NUMBER PROBLEM KEITH CONRAD 1. Introduction A right triangle is called rational when its legs and hypotenuse are all rational numbers. Examples of rational right triangles include Pythagorean
More informationA Diophantine System and a Problem on Cubic Fields
International Mathematical Forum, Vol. 6, 2011, no. 3, 141-146 A Diophantine System and a Problem on Cubic Fields Paul D. Lee Department of Mathematics and Statistics University of British Columbia Okanagan
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 informationHankel Operators plus Orthogonal Polynomials. Yield Combinatorial Identities
Hanel Operators plus Orthogonal Polynomials Yield Combinatorial Identities E. A. Herman, Grinnell College Abstract: A Hanel operator H [h i+j ] can be factored as H MM, where M maps a space of L functions
More information