Polynomial Formula for Sums of Powers of Integers
|
|
- Dwayne O’Neal’
- 5 years ago
- Views:
Transcription
1 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 a Professor in the Electrical Engineering Department at Iowa State University, Ames, Iowa. Keywords Sums of powers, integers, polynomial formula. In this article, it is shown that for any positive integer k, there exist unique real numbers a kr, r=,,...,k +, such that for any integer n k+ n S k,n j k = a kr n r. j= r= The numbers a kr are computed explicitly for r = k +,k,k,...,k 0. This fully determines the polynomials for k =,,...,. The cases k =,, 3 are well known and available in high school algebra books.. Introduction Let k and n be positive integers and S k,n = n j= jk.a well-known story in the history of mathematics is that the great German mathematician Carl Friedrich Gauss, while in elementary school, noted that for k =,S,n n is, when written in reverse order, equal to n+n And when they are added one gets S,n =+n++n+3+n+...+n+ = n +n. This yields S,n = nn+. Shailesh Shirali has told the authors that the great Indian mathematician of ancient times, Aryabhata [], has explicitly mentioned in one of his verses the formulas for the cases k = and k = 3. The three cases, k =,, 3, are now in high school algebra texts. Proofs of these formulas are given using the principle of induction. In this article, we establish that for any integer k there exist unique real numbers a kr,r=,,...,k+, such that for any positive integer n, 7 RESONANCE August 05
2 S k,n GENERAL ARTICLE n k+ j k = a kr n r. j= r= In particular, it is shown that for any k, a kk+ = k +,a kk = ; for k,a kk = k ; for k 3,a kk =0; kk k for k,a kk3 = ; 70 for k 5,a kk =0; kk k k 3k for k,a kk5 = ; 7! for k 7,a kk =0; 3kk...k for k 8,a kk7 = ;!7! for k 9,a kk8 =0; kk...k 8 for k 0,a kk9 =0 ;! and for k,a kk0 =0. It is to be noted that once a formula is proposed, its validity may be verified by the principle of induction. The main problem is to guess what form the formula takes. We guessed that it should be a polynomial based on the known formulas for k =,, 3 and found the explicit polynomial by a recurrence relation. Sury [] has noted that Euler had shown that for any positive integers k and n the sum S k,n can be expressed in terms of Bernoulli polynomials. Sury has also pointed out that Euler s work implies that for each positive integer k the sum S k,n is a polynomial in n. The main contribution of the present article is to determine this polynomial explicitly. We do so from elementary methods, thus making the article accessible to our undergraduate students. RESONANCE August 05 77
3 . An Explicit Formula Fix a positive integer k. For a positive integer n, let n fn S k,n j k j= and set f0 = 0. Then {fn} n 0 satisfies the recurrence relation fn +=fn+n + k,n=0,,, with initial condition f0 = 0. It is clear that {fn} n 0 is uniquely determined by and. Thus, if we find real numbers a kj,j =,...,k + such that g k n k+ r= a kr n r, n 0 3 satisfies and, then g k n hastoequals kn for all n. Since g k 0 = 0 and hence holds, to ensure that g k satisfies it suffices to ensure that the coefficients of powers of n on both sides of the following equation k+ r= are equal. g k n +=g k n+n + k, i.e., a kr n + r = k+ r= This leads to the following conditions: a kr n r +n + k, 5 coeff of n k+ a kk+ = a kk+ k + coeff of n k a kk+ + a kk = a kk RESONANCE August 05
4 k + k coeff of n k a kk+ + a kk k + a kk = a kk + 8 k + k coeff of n k a kk+ + a kk 3 k +a kk + a kk = a kk k + 9 and more generally, k + k coeff of n kr a kk+ + a kk r + r k r a k kr+ k +a kkr = a kkr +.0 r Now provides no information but 7, 8, etc. do. Indeed 7 holds iff k + a kk+ = and 8 holds iff k a kk = k k + a kk+ and 0 holds iff for k r + k k r + k + a kkr+ = a kk+ r r + k a k k... r a kr+ k r +.3 RESONANCE August 05 79
5 for r =0,,,...k. It follows from 7, 8, 9 and 0 that a kr,r=,,, k + are uniquely determined as 0 provides a recursive determination of a kkr+ from the knowledge of a kj for j =k +,...,k r +. Thus, 7 9 yield for any k, a kk+ = k +, a kk = k ; and for k,a kk = k ; and 0 determines a kj for j k. This proves the theorem given in the abstract. Let us call it Theorem A: Theorem A. Let k be a positive integer. Then, there exist unique numbers a kr, r =,,...,k + such that for any integer n>, S k,n n j k = j= k+ r= a kr n r. 3. Another Proof of Theorem A We provide another proof of Theorem A at the end of this equation by first establishing the following: Theorem B. For positive integers k and n let g k n n j= jk.then i g n= nn+ ii For k, n, k +g k n = k k + n r g r n r r= +n k+ +k +n k n. 730 RESONANCE August 05
6 Proof. i This result follows from Gauss s argument: g n= n = n +n implying g n=+n++n n +=n +n. ii For k, n, k + r k k + = r k r= = r= k r= n r g r n k + r n r n j r j= n j r j= where, if n =,weset n j= jr =0forr, n k+ = k + j r k +j k j k+ r j= r=0 n n = j + k+ +n + k + j= n + j= j k+ j= = g k+ n + n + k +g k n n k +g k+ n n k+ = n +k +g k n k +n k n k+ yielding ii. ProofofTheoremA.ByTheoremBi,g n isa polynomial of degree two. Then Theorem B ii implies j k RESONANCE August 05 73
7 that g n is a polynomial of degree three and by induction g k n is a polynomial of degree k +. Further, Theorem B ii also shows that the leading coefficient in g k n isk +.Onecanusethesameformulato show that the coefficient of n k in g k nis. Remark. We now derive explicit expressions for g k n fork =, 3,, 5. From Theorem B ii we deduce that 3 3g n = n g n + n 3 +3n n = 3 n yielding nn + = n 3 + n 3 + n, + n 3 +3n n g n n j= j = n3 3 + n + n. Next, g 3 n = n g n + n g n +n +n 3 n nn + = n n + n n + n +n +n 3 n = n + n 3 + n 3 +n = n +n 3 + n, 73 RESONANCE August 05
8 yielding g 3 n= n + n3 + n. Next, 5 5 5g n = n g n + n g n 5 + n 3 g 3 n + n 5 +5n n 3 nn + = 5 n n +0 n n + n n +0 n 3 + n3 + n + n 5 +5n n = n 5 + n n n n 0 = = n 5 + n 5 + n30 + n 0+n, yielding g n= n5 5 + n + n3 3 n 30. Next, g 5 n = n g n + n g n + n 3 g 3 n + n g n 3 +n +n 5 n n = n + n RESONANCE August
9 n +5 n n + n n +0 n 3 + n3 + n n +5 n n + n3 3 n 30 +n +n 5 n 5 = n + n n 0 0 +n n n = n +3n 5 + n 5 + n3.0 n + n.0, yielding g 5 n= n + n5 + n 5 n. This process can be continued and g k n canbecomputed recursively for all integers k. In Section, we compute g k n fork =,, 3,...,.. A Matrix Method Now that we know that a kr,r=,,...k +are determined uniquely, here is a matrix inversion method to find them. Let Then, since k+ r= ψ k l= l r k,l=,,...,k +. r= a kr l r = S k,l = ψ k l, l=,,...,k +, 5 73 RESONANCE August 05
10 it follows that A k a k a k a kk+ = ψ k ψ k ψ k k +, where A k is the k + k + matrix with entries in the lth row given by l, l,...,l k+, l=,,...,k +. It can be checked that A k x x k+ = 0 0 implies that x =0,x =0,...,x n+ =0. This implies that A k is invertible and hence, for any integer k, a k a k a kk+ 5. An Example = A k ψ k ψ k ψ kk+. 7 We now give an example where the function g k narises in a counting situation. Let S n be the set {0,,,...,n} of integers. How many closed intervals [i,j] are there where i,j S n and i< j? Thereareexactlyn intervals of length one each. These are [i,i +]fori =0,,, n. There are exactly n intervals of length two each. These are [i,i+], i=0,,...,n. And in general there are exactly n k + intervals of length k each. These are RESONANCE August
11 [i,i+ k], i=0,,...,nk. Thus the total number of closed intervals [i,j] withi,j in S n and i<jis n+n +...+nk This is precisely g n. Next, consider the set S n, {i,j i,j S n }. How many squares are there with all vertices in S n, and sides parallel to the axes? There are exactly n such squares whose side length is one. These are with vertices {i,j, i +,j, i,j +, i +,j +} with 0 i n, 0 j n. Next, for any integer l, l n there are exactly n l + squares whose side length is l with vertices all in S n,. These are the ones with vertices {i,j, i+l, j, i,j+l, i+l, j+l} with 0 i n l, 0 j n l. Thus, the total number of squares with all vertices in S n, is n l= n l + = n j= j which is precisely g n. Next, consider the set S n,3 {i,j,k,i,j,k S n }.How many cubes are there with all vertices in S n,3 and faces parallel to the coordinate planes? There are exactly n 3 such cubes whose side length is one. These are the ones with vertices {i,j,k, i +,j,k, i,j +,k, i,j,k +, i+,j+,k,i,j+,k+,i+,j,k+,i+,j+,k+} with 0 i n, 0 j n, 0 k n. Similarly, there are n l + 3 cubes of side length l with vertices in S n,3. Thus, the total number of cubes with vertices in S n,3 is n l= nl+3 = n j= j3 which is precisely g 3 n. Next, for any integer k>3consider the set S n,k {i,i,...,i k i j S n,j =,,...,k}. How many k-dimensional cubes are there with all vertices in S n,k and faces parallel to the coordinate hyperplanes? Arguing as before, the number of such cubes with side length l and vertices in S n,k is n l + k. Thus, the total number of k-dimensional cubes with all vertices in S n,k is n l= n l +k = n j= jk which is precisely g k n. 73 RESONANCE August 05
12 . Determination of a kj for j =k+,k,...,k0 Returning to 7 9 we note that a kk+ = Next, for k 3 k +, a kk = k, a kk = k. a kk = k 3 akk+ k+ = kk 3! k akk 3 k akk k! =0. For k, a kk3 = For k 5, a kk = k ak k+ k+ 5 = kk k = k 5 kk k. 70 akk+ k+ k ark k3 ak k k 3! 5! 3! k k akk 5 akk k ak k3 k = kk k k 3 5!! k3! + 0 =0. RESONANCE August
13 For k, k k + k a kk5 = a kk+ a kk 7 k a kk 5 / k 3 k 5 a k k3 3 = kk k k 3k! 7! 5! ! kk k k 3k =! 7. For k 7, k k + k a kk = a kk+ a kk k k 3 a kk a kk3 / k 5 k a kk5 = kk...k 5 7! 8!! + 70 = kk...k 5 0 = 0.! 7! For k 8, k k + a kk7 = a kk+ 8 9 k k a kk a kk RESONANCE August 05
14 k 3 a kk3 5 / k 5 k 7 a kk5 3 = kk...k 8! 9! 7! ! 3kk...k =. 7!! For k 9, k k + k a kk8 = a kk+ a kk k k 3 a kk a kk3 8 / k 5 k 7 k 5 a kk5 a kk7 = kk...k 7 9! 0! 8! = +!!!! + 3 7!! 5! kk...k ! =0. For k 0, a k k9 = k k + a kk+ 0 RESONANCE August
15 k k a kk a kk 0 9 k 3 k 5 a kk3 a kk 7 5 / k 7 k 9 a kk7 3 = kk...k 8 0!! 70 7!! 5! + 3 7!!3! 9! + kk...k 8 0 =.! For k, k k + k a kk0 = a kk+ a kk k k 3 a kk a kk3 0 8 k 5 k 7 a kk5 a kk7 / k 9 k 0 a kk9 = + kk...k 9! = 0 9 8! ! kk...k 9! 0= RESONANCE August 05
16 Summarizing the above we have: for k a kk+ = a k+ kk = for k a kk = k for k 3 a kk =0 for k a kk3 = kkk 70 for k 5 a kk =0 for k a kk5 = kkkk3k! for k 7 a kk =0 for k 8 a kk7 = 3 kk...k 7!! for k 9 a kk8 =0 for k 0 a kk9 = kk...k80! for k a kk0 =0. 7 This, in turn, yields for integers n, S,n S,n S 3,n S,n S 5,n S,n S 7,n S 8,n S 9,n S 0,n = n + n = n3 3 + n + n = n + n3 + n = n5 5 + n + n3 3 n 30 = n + n5 + n 5 n = n7 7 + n + n5 n3 + n = n8 8 + n7 + n 7 n 7 + n = n9 9 + n8 + n7 3 n n3 9 n 30 = n0 0 + n9 + n83 n n 3 0 n = n + n0 + n95 n7 + n 5 n n RESONANCE August 05 7
17 S,n S,n = n + n + n0 n8 8 + n n 33 + n 5 = n3 3 + n + n n 9 + n7 7 n n n. A natural conjecture from our calculations is that for any k r +,k,rpositive integers, a kkr =0.It may be noted that we have verified it for r =,,...,5 and k r Concluding Remarks 7. Asymptotic Behavior of S k,n For Large n It follows from that lim n S k,n = n k+ n j= n k+ j k a kj n j j= j= n k j k j= = a kk+ + a kjn j j= = j k n k+ = k k + + k +. j= n k+ a kj n j n k+ This also follows from the Riemann sum approximation n j k n = n f k+ k j n n, where f kx=x k, 0 x j= j= which converges to 0 f k xdx = k +. 7 RESONANCE August 05
18 Similarly, the second order behavior is given by n j k = n k+ k + n + k j= j= n n j k n k+ k +. And so on. j= 7. No Induction Involved Suggested Reading a kj n j n k+ For the cases k =,, 3, the proof in elementary texts involves merely the verification of the formula using the principle of induction. Our proofs of the polynomial formula in Theorem A involves no induction but are from first principles. 7.3 Other Treatments There are many papers on this subject. Some of them are listed in the Suggested Reading. [] Walter Eugene Clark, Translation of Aryabhatiya of Aryabhata, The University of Chicago Press, verse, pp.37, 930. [] B Sury, Bernoulli numbers and the Riemann Zeta function, Resonance, No.7, pp.5, 003. [3] A F Beardon, Sums of powers of integers, Am. Math. Monthly, Vol. 03, pp.0 3, 99. [] A W F Edwards, Sums of powers of integers: A Little of the History, The Mathematical Gazette, Vol., No.35, 98. [5] H K Krishnapriyan, Eulerian polynomials and Faulhaber's result on sums of powers of integers, The College Mathematics Journal, Vol., pp.8 3, 995. [] Robert Owens, Sums of powers of integers, Mathematics Magazine, Vol.5, pp.38 0, 99. [7] S Shirali, On sums of powers of integers, Resonance, Indian Academy of Sciences, July 007. Address for Correspondence K B Athreya Department of Mathematics & Statistics Iowa State University Ames, Iowa 500 USA kbathreya@gmail.com RESONANCE August 05 73
Sums of Powers of Integers. Julien Doucet. and A. Saleh-Jahromi
Sums of Powers of Integers by Julien Doucet Louisiana State University in Shreveport Department of Mathematics and Computer Science One University Place Shreveport, LA 75 2399 e-mail: jdoucet@pilot.lsus.edu
More informationMathematical Induction Assignments
1 Mathematical Induction Assignments Prove the Following using Principle of Mathematical induction 1) Prove that for any positive integer number n, n 3 + 2 n is divisible by 3 2) Prove that 1 3 + 2 3 +
More informationGeneralized Akiyama-Tanigawa Algorithm for Hypersums of Powers of Integers
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 16 (2013, Article 1332 Generalized Aiyama-Tanigawa Algorithm for Hypersums of Powers of Integers José Luis Cereceda Distrito Telefónica, Edificio Este
More informationDiscrete Math, Spring Solutions to Problems V
Discrete Math, Spring 202 - Solutions to Problems V Suppose we have statements P, P 2, P 3,, one for each natural number In other words, we have the collection or set of statements {P n n N} a Suppose
More information1 Determinants. 1.1 Determinant
1 Determinants [SB], Chapter 9, p.188-196. [SB], Chapter 26, p.719-739. Bellow w ll study the central question: which additional conditions must satisfy a quadratic matrix A to be invertible, that is to
More informationMATH 1210 Assignment 4 Solutions 16R-T1
MATH 1210 Assignment 4 Solutions 16R-T1 Attempt all questions and show all your work. Due November 13, 2015. 1. Prove using mathematical induction that for any n 2, and collection of n m m matrices A 1,
More informationContinued Fraction Digit Averages and Maclaurin s Inequalities
Continued Fraction Digit Averages and Maclaurin s Inequalities Steven J. Miller, Williams College sjm1@williams.edu, Steven.Miller.MC.96@aya.yale.edu Joint with Francesco Cellarosi, Doug Hensley and Jake
More informationAn Euler-Type Formula for ζ(2k + 1)
An Euler-Type Formula for ζ(k + ) Michael J. Dancs and Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan University Bloomington, IL 670-900, USA Draft, June 30th, 004 Abstract
More informationMath 3000 Section 003 Intro to Abstract Math Homework 6
Math 000 Section 00 Intro to Abstract Math Homework 6 Department of Mathematical and Statistical Sciences University of Colorado Denver, Spring 01 Solutions April, 01 Please note that these solutions are
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 informationMath 324 Summer 2012 Elementary Number Theory Notes on Mathematical Induction
Math 4 Summer 01 Elementary Number Theory Notes on Mathematical Induction Principle of Mathematical Induction Recall the following axiom for the set of integers. Well-Ordering Axiom for the Integers If
More informationHIGHER-ORDER DIFFERENCES AND HIGHER-ORDER PARTIAL SUMS OF EULER S PARTITION FUNCTION
ISSN 2066-6594 Ann Acad Rom Sci Ser Math Appl Vol 10, No 1/2018 HIGHER-ORDER DIFFERENCES AND HIGHER-ORDER PARTIAL SUMS OF EULER S PARTITION FUNCTION Mircea Merca Dedicated to Professor Mihail Megan on
More informationa 2n = . On the other hand, the subsequence a 2n+1 =
Math 316, Intro to Analysis subsequences. This is another note pack which should last us two days. Recall one of our arguments about why a n = ( 1) n diverges. Consider the subsequence a n = It converges
More informationarxiv: v1 [math.nt] 17 Nov 2011
On the representation of k sequences of generalized order-k numbers arxiv:11114057v1 [mathnt] 17 Nov 2011 Kenan Kaygisiz a,, Adem Sahin a a Department of Mathematics, Faculty of Arts Sciences, Gaziosmanpaşa
More informationCombinatorial Proofs and Algebraic Proofs I
Combinatorial Proofs and Algebraic Proofs I Shailesh A Shirali Shailesh A Shirali is Director of Sahyadri School (KFI), Pune, and also Head of the Community Mathematics Centre in Rishi Valley School (AP).
More informationII. Determinant Functions
Supplemental Materials for EE203001 Students II Determinant Functions Chung-Chin Lu Department of Electrical Engineering National Tsing Hua University May 22, 2003 1 Three Axioms for a Determinant Function
More informationSOLVING LINEAR SYSTEMS
SOLVING LINEAR SYSTEMS We want to solve the linear system a, x + + a,n x n = b a n, x + + a n,n x n = b n This will be done by the method used in beginning algebra, by successively eliminating unknowns
More informationLecture Note 2: The Gaussian Elimination and LU Decomposition
MATH 5330: Computational Methods of Linear Algebra Lecture Note 2: The Gaussian Elimination and LU Decomposition The Gaussian elimination Xianyi Zeng Department of Mathematical Sciences, UTEP The method
More informationSome thoughts concerning power sums
Some thoughts concerning power sums Dedicated to the memory of Zoltán Szvetits (929 20 Gábor Nyul Institute of Mathematics, University of Debrecen H 00 Debrecen POBox 2, Hungary e mail: gnyul@scienceunidebhu
More informationCSE 160 Lecture 13. Numerical Linear Algebra
CSE 16 Lecture 13 Numerical Linear Algebra Announcements Section will be held on Friday as announced on Moodle Midterm Return 213 Scott B Baden / CSE 16 / Fall 213 2 Today s lecture Gaussian Elimination
More informationarxiv: v1 [math.co] 18 Oct 2016
ON SOME CLASS OF SUMS ANDREI K. SVININ Abstract. We consider some class of the sums which naturally include the sums of powers of integers. We suggest a number of conjectures concerning a representation
More informationSolving Problems by Inductive Reasoning
Solving Problems by Inductive Reasoning MATH 100 Survey of Mathematical Ideas J. Robert Buchanan Department of Mathematics Fall 2014 Inductive Reasoning (1 of 2) If we observe several similar problems
More information1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )
Direct Methods for Linear Systems Chapter Direct Methods for Solving Linear Systems Per-Olof Persson persson@berkeleyedu Department of Mathematics University of California, Berkeley Math 18A Numerical
More informationMultiple Zeta Values of Even Arguments
Michael E. Hoffman U. S. Naval Academy Seminar arithmetische Geometrie und Zahlentheorie Universität Hamburg 13 June 2012 1 2 3 4 5 6 Values The multiple zeta values (MZVs) are defined by ζ(i 1,..., i
More informationDeterminant formulas for multidimensional hypergeometric period matrices
Advances in Applied Mathematics 29 (2002 137 151 www.academicpress.com Determinant formulas for multidimensional hypergeometric period matrices Donald Richards a,b,,1 andqifuzheng c,2 a School of Mathematics,
More informationCHAPTER 4 SOME METHODS OF PROOF
CHAPTER 4 SOME METHODS OF PROOF In all sciences, general theories usually arise from a number of observations. In the experimental sciences, the validity of the theories can only be tested by carefully
More informationOn generalized Frame-Stewart numbers
On generalized Frame-Stewart numbers Jonathan Chappelon and Akihiro Matsuura August 31, 2010 Abstract For the multi-peg Tower of Hanoi problem with k 4 pegs, so far the best solution is obtained by the
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 informationMathematical Structures Combinations and Permutations
Definitions: Suppose S is a (finite) set and n, k 0 are integers The set C(S, k) of k - combinations consists of all subsets of S that have exactly k elements The set P (S, k) of k - permutations consists
More informationMATRICES. a m,1 a m,n A =
MATRICES Matrices are rectangular arrays of real or complex numbers With them, we define arithmetic operations that are generalizations of those for real and complex numbers The general form a matrix of
More informationProblem Solving in Math (Math 43900) Fall 2013
Problem Solving in Math (Math 43900) Fall 203 Week six (October ) problems recurrences Instructor: David Galvin Definition of a recurrence relation We met recurrences in the induction hand-out. Sometimes
More informationEx. Here's another one. We want to prove that the sum of the cubes of the first n natural numbers is. n = n 2 (n+1) 2 /4.
Lecture One type of mathematical proof that goes everywhere is mathematical induction (tb 147). Induction is essentially used to show something is true for all iterations, i, of a sequence, where i N.
More informationarxiv: v1 [cs.sc] 29 Jul 2009
An Explicit Construction of Gauss-Jordan Elimination Matrix arxiv:09075038v1 [cssc] 29 Jul 2009 Yi Li Laboratory of Computer Reasoning and Trustworthy Computing, University of Electronic Science and Technology
More informationn(n 1) = (n + 1)n(n 1)/ ( 1) n(n 1)(n 2) = (n + 1)n(n 1)(n 2)/4
NOTES FROM THE FIRST TWO CLASSES MIKE ZABROCKI - SEPTEMBER 6 &, 202 main idea of this class + 2 + 3 + + n = n(n + )/2 to r + 2 r + + n r =??? Just to show what we are up against: + 2 + 3 + + n = n(n +
More informationBernoulli Numbers and their Applications
Bernoulli Numbers and their Applications James B Silva Abstract The Bernoulli numbers are a set of numbers that were discovered by Jacob Bernoulli (654-75). This set of numbers holds a deep relationship
More informationMATH 301 INTRO TO ANALYSIS FALL 2016
MATH 301 INTRO TO ANALYSIS FALL 016 Homework 04 Professional Problem Consider the recursive sequence defined by x 1 = 3 and +1 = 1 4 for n 1. (a) Prove that ( ) converges. (Hint: show that ( ) is decreasing
More informationCS1800: Mathematical Induction. Professor Kevin Gold
CS1800: Mathematical Induction Professor Kevin Gold Induction: Used to Prove Patterns Just Keep Going For an algorithm, we may want to prove that it just keeps working, no matter how big the input size
More information. Get closed expressions for the following subsequences and decide if they converge. (1) a n+1 = (2) a 2n = (3) a 2n+1 = (4) a n 2 = (5) b n+1 =
Math 316, Intro to Analysis subsequences. Recall one of our arguments about why a n = ( 1) n diverges. Consider the subsequences a n = ( 1) n = +1. It converges to 1. On the other hand, the subsequences
More informationElementary Row Operations on Matrices
King Saud University September 17, 018 Table of contents 1 Definition A real matrix is a rectangular array whose entries are real numbers. These numbers are organized on rows and columns. An m n matrix
More informationAveraging sums of powers of integers and Faulhaber polynomials
Annales Mathematicae et Informaticae 42 (20. 0 htt://ami.ektf.hu Averaging sums of owers of integers and Faulhaber olynomials José Luis Cereceda a a Distrito Telefónica Madrid Sain jl.cereceda@movistar.es
More informationCOUNTING NUMERICAL SEMIGROUPS BY GENUS AND SOME CASES OF A QUESTION OF WILF
COUNTING NUMERICAL SEMIGROUPS BY GENUS AND SOME CASES OF A QUESTION OF WILF NATHAN KAPLAN Abstract. The genus of a numerical semigroup is the size of its complement. In this paper we will prove some results
More informationA196837: Ordinary Generating Functions for Sums of Powers of the First n Positive Integers
Karlsruhe October 14, 2011 November 1, 2011 A196837: Ordinary Generating Functions for Sums of Powers of the First n Positive Integers Wolfdieter L a n g 1 The sum of the k th power of the first n positive
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 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 informationMATH 341 MIDTERM 2. (a) [5 pts] Demonstrate that A and B are row equivalent by providing a sequence of row operations leading from A to B.
11/01/2011 Bormashenko MATH 341 MIDTERM 2 Show your work for all the problems. Good luck! (1) Let A and B be defined as follows: 1 1 2 A =, B = 1 2 3 0 2 ] 2 1 3 4 Name: (a) 5 pts] Demonstrate that A and
More informationSUMS OF POWERS AND BERNOULLI NUMBERS
SUMS OF POWERS AND BERNOULLI NUMBERS TOM RIKE OAKLAND HIGH SCHOOL Fermat and Pascal On September 22, 636 Fermat claimed in a letter that he could find the area under any higher parabola and Roberval wrote
More informationON A KIND OF GENERALIZED ARITHMETIC-GEOMETRIC PROGRESSION
ON A KIND OF GENERALIZED ARITHMETIC-GEOMETRIC PROGRESSION Leetsch C. Hsu* Dalian University of Technology, Dalian 116024, China (Submitted August 1995) 1. INTRODUCTION Let (t\h) p denote the generalized
More informationDirect Methods for Solving Linear Systems. Simon Fraser University Surrey Campus MACM 316 Spring 2005 Instructor: Ha Le
Direct Methods for Solving Linear Systems Simon Fraser University Surrey Campus MACM 316 Spring 2005 Instructor: Ha Le 1 Overview General Linear Systems Gaussian Elimination Triangular Systems The LU Factorization
More informationChapter 1 Matrices and Systems of Equations
Chapter 1 Matrices and Systems of Equations System of Linear Equations 1. A linear equation in n unknowns is an equation of the form n i=1 a i x i = b where a 1,..., a n, b R and x 1,..., x n are variables.
More informationWeek 15-16: Combinatorial Design
Week 15-16: Combinatorial Design May 8, 2017 A combinatorial design, or simply a design, is an arrangement of the objects of a set into subsets satisfying certain prescribed properties. The area of combinatorial
More information*!5(/i,*)=t(-i)*-, fty.
ON THE DIVISIBILITY BY 2 OF THE STIRLING NUMBERS OF THE SECOND KIND T* Lengyel Occidental College, 1600 Campus Road., Los Angeles, CA 90041 {Submitted August 1992) 1. INTRODUCTION In this paper we characterize
More informationVector analysis and vector identities by means of cartesian tensors
Vector analysis and vector identities by means of cartesian tensors Kenneth H. Carpenter August 29, 2001 1 The cartesian tensor concept 1.1 Introduction The cartesian tensor approach to vector analysis
More informationSolve Linear System with Sylvester s Condensation
International Journal of Algebra, Vol 5, 2011, no 20, 993-1003 Solve Linear System with Sylvester s Condensation Abdelmalek Salem Department of Mathematics University of Tebessa, 12002 Algeria And Faculty
More informationSequences of Real Numbers
Chapter 8 Sequences of Real Numbers In this chapter, we assume the existence of the ordered field of real numbers, though we do not yet discuss or use the completeness of the real numbers. In the next
More informationLinear Systems and Matrices
Department of Mathematics The Chinese University of Hong Kong 1 System of m linear equations in n unknowns (linear system) a 11 x 1 + a 12 x 2 + + a 1n x n = b 1 a 21 x 1 + a 22 x 2 + + a 2n x n = b 2.......
More informationLatin Squares and Orthogonal Arrays
School of Electrical Engineering and Computer Science University of Ottawa lucia@eecs.uottawa.ca Winter 2017 Latin squares Definition A Latin square of order n is an n n array, with symbols in {1,...,
More informationExpected Number of Distinct Subsequences in Randomly Generated Binary Strings
Expected Number of Distinct Subsequences in Randomly Generated Binary Strings arxiv:704.0866v [math.co] 4 Mar 08 Yonah Biers-Ariel, Anant Godbole, Elizabeth Kelley March 6, 08 Abstract When considering
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More information2.1 Convergence of Sequences
Chapter 2 Sequences 2. Convergence of Sequences A sequence is a function f : N R. We write f) = a, f2) = a 2, and in general fn) = a n. We usually identify the sequence with the range of f, which is written
More informationTHREE LECTURES ON QUASIDETERMINANTS
THREE LECTURES ON QUASIDETERMINANTS Robert Lee Wilson Department of Mathematics Rutgers University The determinant of a matrix with entries in a commutative ring is a main organizing tool in commutative
More informationECE608 Midterm 1 Spring 2009 THEN PUT YOUR ADDRESS ON EVERY PAGE OF THE EXAM! NAME: Neatness counts. We will not grade what we cannot read.
ECE608 Midterm Spring 2009 Please fill in the following information. THEN PUT YOUR EMAI ADDRESS ON EVERY PAGE OF THE EXAM! NAME: Email (please put complete address): Neatness counts. We will not grade
More informationEquations in Quadratic Form
Equations in Quadratic Form MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to: make substitutions that allow equations to be written
More informationChapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations
Chapter 1: Systems of linear equations and matrices Section 1.1: Introduction to systems of linear equations Definition: A linear equation in n variables can be expressed in the form a 1 x 1 + a 2 x 2
More informationCS137 Introduction to Scientific Computing Winter Quarter 2004 Solutions to Homework #3
CS137 Introduction to Scientific Computing Winter Quarter 2004 Solutions to Homework #3 Felix Kwok February 27, 2004 Written Problems 1. (Heath E3.10) Let B be an n n matrix, and assume that B is both
More informationUniversity of the Virgin Islands, St. Thomas January 14, 2015 Algorithms and Programming for High Schoolers. Lecture 5
University of the Virgin Islands, St. Thomas January 14, 2015 Algorithms and Programming for High Schoolers Numerical algorithms: Lecture 5 Today we ll cover algorithms for various numerical problems:
More informationThe Truncated Pentagonal Number Theorem
The Truncated Pentagonal Number Theorem George E. Andrews Department of Mathematics The Pennsylvania State University University Park, PA 16802 USA Mircea Merca Doctoral School in Applied Mathematics University
More informationBounds for pairs in partitions of graphs
Bounds for pairs in partitions of graphs Jie Ma Xingxing Yu School of Mathematics Georgia Institute of Technology Atlanta, GA 30332-0160, USA Abstract In this paper we study the following problem of Bollobás
More informationThe cocycle lattice of binary matroids
Published in: Europ. J. Comb. 14 (1993), 241 250. The cocycle lattice of binary matroids László Lovász Eötvös University, Budapest, Hungary, H-1088 Princeton University, Princeton, NJ 08544 Ákos Seress*
More informationLaura Chihara* and Dennis Stanton**
ZEROS OF GENERALIZED KRAWTCHOUK POLYNOMIALS Laura Chihara* and Dennis Stanton** Abstract. The zeros of generalized Krawtchouk polynomials are studied. Some interlacing theorems for the zeros are given.
More information2.1 Gaussian Elimination
2. Gaussian Elimination A common problem encountered in numerical models is the one in which there are n equations and n unknowns. The following is a description of the Gaussian elimination method for
More informationChapter Four Gelfond s Solution of Hilbert s Seventh Problem (Revised January 2, 2011)
Chapter Four Gelfond s Solution of Hilbert s Seventh Problem (Revised January 2, 2011) Before we consider Gelfond s, and then Schneider s, complete solutions to Hilbert s seventh problem let s look back
More informationSeries of Error Terms for Rational Approximations of Irrational Numbers
2 3 47 6 23 Journal of Integer Sequences, Vol. 4 20, Article..4 Series of Error Terms for Rational Approximations of Irrational Numbers Carsten Elsner Fachhochschule für die Wirtschaft Hannover Freundallee
More informationEstimates for probabilities of independent events and infinite series
Estimates for probabilities of independent events and infinite series Jürgen Grahl and Shahar evo September 9, 06 arxiv:609.0894v [math.pr] 8 Sep 06 Abstract This paper deals with finite or infinite sequences
More informationProblem 1A. Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1. (Hint: 1. Solution: The volume is 1. Problem 2A.
Problem 1A Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1 (Hint: 1 1 (something)dz) Solution: The volume is 1 1 4xydz where x = y = 1 z 2 This integral has value 16/3 Problem 2A Let f(x)
More informationACI-matrices all of whose completions have the same rank
ACI-matrices all of whose completions have the same rank Zejun Huang, Xingzhi Zhan Department of Mathematics East China Normal University Shanghai 200241, China Abstract We characterize the ACI-matrices
More informationLargest Values of the Stern Sequence, Alternating Binary Expansions and Continuants
1 3 47 6 3 11 Journal of Integer Sequences, Vol. 0 (017), Article 17..8 Largest Values of the Stern Sequence, Alternating Binary Expansions Continuants Rol Paulin Max-Planck-Institut für Mathematik Vivatsgasse
More informationAssignment 10. Arfken Show that Stirling s formula is an asymptotic expansion. The remainder term is. B 2n 2n(2n 1) x1 2n.
Assignment Arfken 5.. Show that Stirling s formula is an asymptotic expansion. The remainder term is R N (x nn+ for some N. The condition for an asymptotic series, lim x xn R N lim x nn+ B n n(n x n B
More informationCSL361 Problem set 4: Basic linear algebra
CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices
More informationJordan normal form notes (version date: 11/21/07)
Jordan normal form notes (version date: /2/7) If A has an eigenbasis {u,, u n }, ie a basis made up of eigenvectors, so that Au j = λ j u j, then A is diagonal with respect to that basis To see this, let
More informationPrinciple of Mathematical Induction
Advanced Calculus I. Math 451, Fall 2016, Prof. Vershynin Principle of Mathematical Induction 1. Prove that 1 + 2 + + n = 1 n(n + 1) for all n N. 2 2. Prove that 1 2 + 2 2 + + n 2 = 1 n(n + 1)(2n + 1)
More informationFormula for the inverse matrix. Cramer s rule. Review: 3 3 determinants can be computed expanding by any row or column
Math 20F Linear Algebra Lecture 18 1 Determinants, n n Review: The 3 3 case Slide 1 Determinants n n (Expansions by rows and columns Relation with Gauss elimination matrices: Properties) Formula for the
More informationProbabilistic Proofs of Existence of Rare Events. Noga Alon
Probabilistic Proofs of Existence of Rare Events Noga Alon Department of Mathematics Sackler Faculty of Exact Sciences Tel Aviv University Ramat-Aviv, Tel Aviv 69978 ISRAEL 1. The Local Lemma In a typical
More informationProblem Set 5 Solutions
Problem Set 5 Solutions Section 4.. Use mathematical induction to prove each of the following: a) For each natural number n with n, n > + n. Let P n) be the statement n > + n. The base case, P ), is true
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 informationDEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix to upper-triangular
form) Given: matrix C = (c i,j ) n,m i,j=1 ODE and num math: Linear algebra (N) [lectures] c phabala 2016 DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix
More informationMany proofs that the primes are infinite J. Marshall Ash 1 and T. Kyle Petersen
Many proofs that the primes are infinite J. Marshall Ash and T. Kyle Petersen Theorem. There are infinitely many prime numbers. How many proofs do you know that there are infinitely many primes? Nearly
More informationHOMEWORK #2 - MATH 3260
HOMEWORK # - MATH 36 ASSIGNED: JANUARAY 3, 3 DUE: FEBRUARY 1, AT :3PM 1) a) Give by listing the sequence of vertices 4 Hamiltonian cycles in K 9 no two of which have an edge in common. Solution: Here is
More informationALGEBRAIC GEOMETRY I - FINAL PROJECT
ALGEBRAIC GEOMETRY I - FINAL PROJECT ADAM KAYE Abstract This paper begins with a description of the Schubert varieties of a Grassmannian variety Gr(k, n) over C Following the technique of Ryan [3] for
More informationHow to Discover the Rogers Ramanujan Identities
How to Discover the Rogers Ramanujan Identities Gaurav Bhatnagar We examine a method to conjecture two very famous identities that were conjectured by Ramanujan, and later found to be known to Rogers..
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K. R. MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Corrected Version, 7th April 013 Comments to the author at keithmatt@gmail.com Chapter 1 LINEAR EQUATIONS 1.1
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 informationMathematical Induction Again
Mathematical Induction Again James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University January 12, 2017 Outline Mathematical Induction Simple POMI Examples
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 informationThe Heine-Borel and Arzela-Ascoli Theorems
The Heine-Borel and Arzela-Ascoli Theorems David Jekel October 29, 2016 This paper explains two important results about compactness, the Heine- Borel theorem and the Arzela-Ascoli theorem. We prove them
More informationMathematical Induction Again
Mathematical Induction Again James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University January 2, 207 Outline Mathematical Induction 2 Simple POMI Examples
More informationNotes on Determinants and Matrix Inverse
Notes on Determinants and Matrix Inverse University of British Columbia, Vancouver Yue-Xian Li March 17, 2015 1 1 Definition of determinant Determinant is a scalar that measures the magnitude or size of
More informationREVIEW FOR EXAM II. The exam covers sections , the part of 3.7 on Markov chains, and
REVIEW FOR EXAM II The exam covers sections 3.4 3.6, the part of 3.7 on Markov chains, and 4.1 4.3. 1. The LU factorization: An n n matrix A has an LU factorization if A = LU, where L is lower triangular
More informationLinks Between Sums Over Paths in Bernoulli s Triangles and the Fibonacci Numbers
Links Between Sums Over Paths in Bernoulli s Triangles and the Fibonacci Numbers arxiv:1611.09181v1 [math.co] 28 Nov 2016 Denis Neiter and Amsha Proag Ecole Polytechnique Route de Saclay 91128 Palaiseau
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 998 Comments to the author at krm@mathsuqeduau All contents copyright c 99 Keith
More informationP i [B k ] = lim. n=1 p(n) ii <. n=1. V i :=
2.7. Recurrence and transience Consider a Markov chain {X n : n N 0 } on state space E with transition matrix P. Definition 2.7.1. A state i E is called recurrent if P i [X n = i for infinitely many n]
More information