arxiv:math/ v1 [math.co] 7 Jul 2004

Transcription

1 Newton s Method as a Formal Recurrence arxiv:math/ v1 [math.co] 7 Jul 2004 Hal Canary, Carl Eduist, Samuel Lachterman, Brendan Younger February 1, 2008 Abstract Iterating Newton s method symbolically for the general uadratic ax 2 +bx+c yields a rational function P nx) Q nx), the numerator and denominator of which are polynomials with highly composite coefficients. In particular, the coefficients have no prime factors greater than 2 n after n iterations even though most of the coefficients are much larger than 2 n. 1 Introduction If f : C C is a differentiable function, then Newton s method, applied recursively to an initial value of z 0, yields the seuence of values z 1, z 2,... defined as z n+1 = Nz n ) = z n fz n) f z n ) which, in many cases, converges to a root of f. We restrict our attention to the general uadratic fx) = ax 2 + bx + c for the purposes of this paper. Instead of using euation 1) as a numerical method, we are interested in iterating symbolically by letting the initial guess be z 0 = x where x C. Doing so, we obtain z n = P nx) Q n x) where P n x) and Q n x) are relatively prime polynomials in C[x]. An interesting observation, pointed out to us by Jim Propp, is that the coefficients of these polynomials, though very large, have no prime factors 2 n. On his suggestion that this divisibility property might imply some combinatorial interpretation not immediately apparent in the formulation of euation 1), we have derived an explicit symbolic formula for each iterate. The aim of this paper is to show that the polynomials P n and Q n are given by the following explicit formulas: P n x) = a 2 x 2n 2 n 2 ) 2 2 a + b 2n 2 2 c +1 x 1) 2a) Q n x) = 2 n 1 ) 2 1 a + b 2n 2 1 c x The binomial coefficients immediately explain the maximum size of the prime divisors of the coefficients of these polynomials. Furthermore, euation 2) provides a good starting point for investigating any combinatorial meaning of P n and Q n, which is discussed at the end of this paper. 2b) 1

2 2 Fractional Linear Transformations Let C denote the Riemann sphere C { }. A fractional linear transformation is a map β : C C of the form βτ) = aτ + b cτ + d where a, b, c, d C and ad bc 0. This function is a conformal map which is analytic everywhere except the pole at τ = d c. For our purposes, we define a particular fractional linear transformation ϕτ) as ϕτ) = τ r 1 τ r 2 3) where r 1 and r 2 are the two distinct roots of the uadratic polynomial fx) = ax 2 + bx + c. This particular approach exploits the fact that Newton s method for uadratics with distinct roots) is conugate to z z 2, with respect to ϕ, the fractional linear transformation which sends the roots of the uadratic to 0 and. Ric Kenyon [2] was the first to point this out to us; expositions can be found in both Cayley [1] and McMullen [3]. For distinct roots r 1 and r 2, the fractional linear transformation ϕ 1 τ) exists and the formula from euation 1) can be expressed as Nz n ) = z n az2 n + bz n + c 2az n + b That is, the following diagram of maps is commutative: = ϕ 1 ϕz n ) 2 ) ϕ C N C ϕ C C ) 2 Thus, the formula for the nth iterate of Newton s method is ϕ 1 ϕx) 2n ) = r 1x r 2 ) 2n r 2 x r 1 ) 2n x r 2 ) 2n x r 1 ) 2n = P nx) Q n x) Using this formula, we shall prove the euations in 2). 4) 3 Proof of Explicit Formula Theorem 1. Given a uadratic polynomial fx) = ax 2 + bx + c, with b 2 4ac 0, define polynomials P n and Q n as follows. 1 P n x) = a2n r 1 x r 2 ) 2n r 2 x r 1 ) 2n) 5a) Q n x) = a2 x r 2 ) 2n x r 1 ) 2n) 5b) Then the polynomials given in euation 5) are eual to those in euation 2) which are reproduced below. P n x) = a 2 x 2n Q n x) = 2 n 1 2 n 2 ) 2 2 a + b 2n 2 2 c +1 x ) 2 1 a + b 2n 2 1 c x 2

3 Before we begin the proof of the theorem, we will need the following lemma. Lemma 1. For all x, y R, and all n N, n 1, the following identity holds ) n i 1 x n y n = x y) 1) i x + y) n 2i 1 xy) i 7) i i=0 Proof of Lemma 1. We proceed by induction. For n = 1, 2 the euality is easily verified, so assume the identity is valid for all n. Let Tn) be the right-hand side of euation 7). Since we need only prove that x n+1 y n+1 = x + y)x n y n ) xyx y ) Tn + 1) = x + y)tn) xytn 1) 8) Combining the sums on the right side of euation 8) by shifting indices, we have [ n ) 1 [ ) )] ] n i 1 n i 1 x y) x + y) n + 1) i + x + y) n 2i xy) i 0 i i 1 i=1 Applying Pascal s identity and bringing the leftmost term inside the sum, we obtain ) n i x y) 1) i x + y) n 2i xy) i i i=0 Finally, we note that for n 1 the binomial coefficient is 0 when i = n. Therefore, this last expression is eual to Tn + 1), and so we are done. Proof of Theorem 1. Begin by expanding euation 5) via the binomial theorem to obtain P n x) = a2 = a2 = a2 Q n x) = a2 2 n 2 n [ 2 n)x [ ] r 1 r 2 ) 2n r 2 r 1 ) 2n 2 n)x [ r 1 r 2 r )] 1 ) 2n 1 r 2 ) 2n 1 )x 2n + 2 n)x [ r 1 r 2 r )]] 1 ) 2n 1 r 2 ) 2n 1 2 n)x [ ] r 2 ) 2n r 1 ) 2n We can now apply lemma 1) to the expressions of the form x n y n in both P n and Q n. Doing so and canceling the factor of since b 2 4ac 0) yields ) P n x) = a 2 2 n 2 n 2 x 2n x 2 1) n ) 2 ) 2n 2 2 r 1 r 2 ) +1 ) Q n x) = a 2 2 n 2 n 1 x 2 1) n ) 1 ) 2n 2 1 r 1 r 2 ) 3

4 Replacing r 1 and r 2 with their values in terms of the coefficients a, b, c gives the final form ) 2 P n x) = a 2 n 2 n 2 x 2n x 2 1) n ) 2 a + b 2n 2 2 c +1 ) 2 n 2 n 1 Q n x) = x 2 1) n ) 1 a + b 2n 2 1 c These euations are the same as euations 2a) and 2b), which was to be proved. 4 Further Remars Though this paper does not discuss any combinatorial interpretation of the polynomials P n and Q n we suspect that their may be some fruitful combinatorial euivalence yet to be discovered. To aid further research in this area, we mae a couple of observations about P n and Q n. Simply iterating Newton s method for the general uadratic gives N Pn Q n ) = ap n 2 cq2 n 2aP n Q n + bq 2 = P n+1 9) n Q n+1 where P 0 x) = x and Q 0 x) = 1 so that the initial term is z 0 = x as before. The uestion is whether the numerator and denominator of euation 9) are relatively prime so that we may define P n and Q n recursively in the natural way. This is, in fact, true as the following lemma proves. Lemma 2. The polynomials P n+1 and Q n+1, defined recursively as and P n+1 = ap 2 n cq 2 n Q n+1 = 2aP n Q n + bq 2 n where P 0 x) = x and Q 0 x) = 1, are relatively prime except possibly in the case b 2 4ac = 0. For the duration of the statement and proof of Lemma 2, we are suspending the definition of P n and Q n given earlier, but it will be an immediate conseuence of Lemma 2 that the two definitions agree.) Proof. We proceed inductively by assuming that P i and Q i are relatively prime for all i n. Assume P n+1 and Q n+1 are not relatively prime to derive a contradiction. Then there exists an irreducible polynomial α such that α P n+1 and α Q n+1. If α Q n, then, since α apn 2 cq 2 n, it follows that α P n which contradicts the induction hypothesis that P n and Q n are relatively prime. Hence α Q n. Since α 2aP n Q n + bq 2 n, we now that α 2aP n +bq n. But then α P n + b 2a Q n and also α P n ± c a Q n. Conseuently, α divides their difference, so α b 2a c a) Qn which only occurs when b 2 4ac = 0. In combinatorics it is sometimes useful to consider two formal variables x, y which do not commute with each other but instead obey yx = xy where is another formal variable that commutes with both x and y. This approach, due to Schützenberger [4], is useful in applications such as counting lattice paths. In our case, the polynomials P n and Q n can be easily generalized to the non-commuting case. As a generalization of the usual binomial coefficient, the -binomial coefficient is defined as [ ] n n 1 i+ = 1 i i=1 4

5 and the following -binomial theorem for non-commuting variables x, y due to Schützenberger [4] is n [ ] n x + y) n = x y n Analogously to the above, we can then define non-commuting polynomials P n and Q n recursively as P n+1 = ap 2 n cq 2 n 10a) Q n+1 = ap nq n + aq np n + bq 2 n where P 0 = x, Q 0 = y, and yx = xy. We conecture explicit formulas for both P n and Q n which happen to be the same as euations 2a) and 2b) except for the presence of a -binomial coefficient: P nx, y) = a 2 x 2n 2 Q n x, y) = 5 Conclusion 2 n 2 2 n 1 [ ] 2 2 a + b 2n 2 2 c +1 x y 2n [ ] 2 1 a + b 2n 2 1 c x y 2n We have taen the initial observation that Newton s method, when applied to uadratics, produces polynomials with highly composite coefficients and proved an explicit formula for the nth iterate that explains this compositeness as a conseuence of the inherent compositeness of binomial coefficients. Furthermore, a recursive definition and a conectural non-commutative analogue of the polynomials P n and Q n were noted in hopes of spurring further research into finding a combinatorial interpretation. We believe that a proof of the non-commutative analogue, as well as the larger issue of finding a combinatorial interpretation, are problems which merit further study. It is also worth noting that for higher-degree polynomials, such as cubics, no similar phenomena have been found. In particular, the occurrence of coefficients with large prime factors indicates that no simple product formulas for the coefficients exist, but this does not rule out the existence of more complicated formulas. We heartily than Jim Propp, Ric Kenyon, and the rest of the Spatial Systems Laboratory at UW- Madison for their helpful insight and generous support for this research. In addition, we are indebted to the NSF s Research Experiences for Undergraduates program and the NSA for funding our research, as well as to the computing staff at UW-Madison for providing the computational resources necessary for our investigations. References [1] A. Cayley. Application of the newton-fourier method to an imaginary root of an euation. In Collected Mathematical Papers, volume XI, chapter 736. Cambridge University Press, [2] R. Kenyon. Personal Communication, March [3] C. McMullen. Families of rational maps and iterative root-finding algorithms. Annals of Math., 125: , [4] M. P. Schützenberger. Une interprétation de certaines solutions de l éuation fonctionnelle: Fx + y) = Fx)Fy). C. R. Acad. Sci. Paris, 236: , b) 11a) 11b) 5