A MATHEMATICA PACKAGE FOR COMPUTING ASYMPTOTIC EXPANSIONS OF SOLUTIONS OF P-FINITE RECURRENCE EQUATIONS. 1. The Problem
|
|
- Allan Wilkinson
- 5 years ago
- Views:
Transcription
1 A MATHEMATICA PACKAGE FOR COMPUTING ASYMPTOTIC EXPANSIONS OF SOLUTIONS OF P-FINITE RECURRENCE EQUATIONS MANUEL KAUERS Abstract. We describe a simple package for computig a fudametal system of certai formal series solutios up to a prescribed order of a give P-fiite recurrece equatio. These solutios ca be viewed as describig the asymptotic behavior of sequeces satisfyig the recurrece. 1. The Problem The Mathematica code described below solves the followig problem: Give a liear recurrece equatio of order r with polyomial coefficiets (also kow as a P-fiite recurrece) p 0 ()a + p 1 ()a p r ()a +r = 0 Fid r liearly idepedet solutios of the form γ e u(1/r) ρ α( ( 0 + β 01 + β 0 + β ) + ( 0 + β 11 + β 1 + β ) log() + + ( 1 + β k 11 + β k 1 + β k ) log() k 1) where αγρβ ij are costats rsk are positive itegers ad u is a polyomial. It is well kow ad ot difficult to compute this data [6] ad besides our implemetatio described below there are several others which do the same job [7 1]. The Mathematica code I programmed a couple of years ago was origially iteded for private use oly but over the time I have also give copies to other users ad this ofte raised the questio o how to use the package ad how exactly to iterpret its output. The purpose of this techical report is to aswer these questios. A priori the output is oly correct i some formal algebraic sese. But accordig to Birkhoff ad Trjitzisky [] it is also correct aalytically i the sese that every sequece (a ) satisfyig the iput recurrece has a liear combiatio of the output series as its asymptotic Supported by the Austria FWF grat Y464-N18. 1
2 MANUEL KAUERS expasio. This meas i particular that we ca obtai asymptotic estimates of the form a = c γ e u(1/r) ρ α( ( 0 + β β 0N N + O( 1 N+1 )) + + ( 1 + β k β k 1N N + O ( 1 N+1 )) log() k 1 ) ( ) for ay prescribed N ad some costat c which ca be determied umerically to high accuracy. The also a c γ e u(1/r) ρ α log() k 1 ( ) a where a b meas asymptotic equivalece i the sese that lim b = 1. Ulike the algebraic correctess the aalytic correctess is ot easy to show. It should be remarked that the argumets give by Birkhoff ad Trjitzisky are log ad complicated ad that some people hesitate to believe i them. O cocrete examples however I have ever observed ay mismatch betwee the actual asymptotic behavior of a solutio (a ) ad the formal expasios. A alterative approach whose uderlyig aalytic theory is more widely accepted but which would require more programmig effort is to go via the geeratig fuctio A(z) := =0 a z of the sequece (a ) uder cosideratio. Flajolet ad Sedgewick [3] give a comprehesive accout o the correspodece betwee the asymptotic behavior of A(z) ear its sigularities closest to the origi ad the asymptotic behavior of the sequece (a ). Salvy s Maple package gdev [5] follows this approach.. The Package The package is available for dowload at the URL It was writte ad tested for Mathematica 6 ad 7. Two commads are provided by the package. The first Asymptotics takes a recurrece as iput ad returs the domiat term of all its (formal) asymptotic solutios. I[1]:= << Asymptotics.m Asymptotics Package by Mauel Kauers c RISC Liz V 0.3 ( ) I[]:= Asymptotics[( + 1)f[ + ] (3 + )f[ + 1] + ( + )f[] f[]] { (3 5 ) (1+ ( ) 5)/ (1 5)/ Out[]= For processig this output further it ca at times be more coveiet to get it ot as a expressio. For this case we offer the possibility to receive the output i a iteral format: I[3]:= Asymptotics[(+1)f[+] (3+)f[+1]+(+)f[] f[] Retur iteral ] { { { Out[3]= { 1+ 5 { { 1 5 {1
3 ASYMPTOTIC EXPANSIONS 3 I this format a solutio () γ exp(µ1 1/r + µ /r + + µ r 1 1 1/r )ρ α e ( (0 β β 0N N + O( 1 )) N ( 1 + β k β k 1N N + O ( 1 )) ) log() k 1 N+1 is represeted i the form { γ {µ1... µ r 1 ρα {0β {1β k Here is a example with a logarithmic term: I[4]:= Asymptotics[( + )f[ + ] ( + 3)f[ + 1] + ( + 1)f[] f[]] { 1 Out[4]= 1 + log() I[5]:= Asymptotics[(+)f[+] (+3)f[+1]+(+1)f[] f[] Retur iteral ] Out[5]= {{0 1 { 0 {1 {0 1 { 0 {0 1 {1 0 By default the commad determies oly the domiat terms of the solutio. If more terms of the expasio are desired they ca be requested by explicitly specifyig the order. I[6]:= Asymptotics[( + 1)f[ + ] (3 + )f[ + 1] + ( + )f[] f[]order 3] { (3 5 ) (1+ ( 5)/ ) 3 ( ) (1 ( ) 5)/ I[7]:= Asymptotics[( + 1)f[ + ] (3 + )f[ + 1] + ( + )f[] f[]order 5] { (3 5 ) (1+ ( 5)/ ) 5 ( ) (1 ( ) 5)/ Out[6]= Out[7]= Ofte oly the higher order terms of oe of the solutios are of iterest or oe computes with big effort the first terms of a expasio ad realizes oly afterwards that some additioal terms are eeded. For these cases istead of wastig computatio time ito the computatio of uiterestig data or ito the recomputatio of data which is already kow there is the secod commad of the package FurtherTerms which takes as iput oe trucated series solutio i the iteral format ad completes it to a trucated series solutio of the specified order. I[8]:= sols = Asymptotics[( + 1)f[ + ] (3 + )f[ + 1] + ( + )f[] f[] Order 3 Retur iteral ] { { Out[8]= { 1+ 5 { { { 1 5 { I[9]:= FurtherTerms[( + 1)f[ + ] (3 + )f[ + 1] + ( + )f[] f[]sols[[]]5] { { 1 5 { Out[9]=
4 4 MANUEL KAUERS I[]:= FurtherTerms[( + 1)f[ + ] (3 + )f[ + 1] + ( + )f[] f[]%15] { { 1 5 { Out[]= The Multiplicative Costat The sequece solutios of a P-fiite recurrece equatio p 0 ()a + p 1 ()a p r ()a +r = 0 form a vector space of fiite dimesio. A particular solutio (a ) ca therefore be characterized uiquely by a fiite umber of iitial values a 0 a 1... a N. Give such a particular solutio we may woder about its asymptotic behavior. Clearly if we kow r liearly idepedet asymptotic solutios (b (1) )... (b (r) ) (e.g. from the Asymptotics commad) the we have a c 1 b (1) + + c r b (r) ( ) for certai costats c 1... c r. I geeral these costats caot be computed i closed form but it is possible to obtai very accurate umerical approximatios to them. Here is how A Sigle Domiat Term. First assume for simplicity that the asymptotic solutios (b (1) )... (b (r) ) are such that oe of them asymptotically domiates all the others b say lim (1) = 0 for all k Æ ad i =... r. The the terms (b () )... (b (r) ) are too k b (i) small to cotribute to the asymptotics ad we have i fact a c 1 b (1) ( ) i.e. there is oly a sigle costat c 1 to be determied. Also if fier asymptotic estimates with higher order terms are cosidered oly those comig from (b (1) ) will play a role. I this case the computatio of c 1 is easy. Sice we have lim approximatios for c 1 by computig the quotiet b (1) b (1) a = c 1 we ca obtai decet /a for some large idex. The higher the idex ad the more terms of the asymptotic expasio (b (1) ) are take ito accout the more digits of the quotiet b (1) /a will agree with the digits of the actual costat c 1. I the followig typical example sessio we defie the sequece (a ) by two iitial values ad a recurrece of secod order. This allows to compute a for every specific idex (e.g. for = 50). Next we determie the domiat terms of the two asymptotic solutios of the recurrece. The secod domiates the first because the basis of its expoetial term is greater tha i the first. Next we determie 15 terms of the asymptotic expasio. The
5 ASYMPTOTIC EXPANSIONS 5 approximate value of the costat is fially obtaied by computig the quotiet b[]/a[] for large idices. I[11]:= a[0] = 1; a[1] = ; a[ Iteger] := a[] = ((3( )+)a[ 1] a[ ])/( 1); I[1]:= a[50] Out[1]= I[13]:= Asymptotics[( + 1)f[ + ] (3 + )f[ + 1] + ( + )f[] f[]] { (3 5 ) (1+ ( ) Out[13]= 5)/ (1 5)/ I[14]:= {N[(3 5)/]N[(3 + 5)/] Out[14]= { I[15]:= terms = Last[FurtherTerms[( + 1)f[ + ] (3 + )f[ + 1] + ( + )f[] f[] { { 1 5 {115]]; I[16]:= b[ Iteger] := ( 3+ ) 5 (1 5)/ Sum[terms[[k + 1]] k {k015]; I[17]:= $RecursioLimit = 5 ; I[18]:= N[b[]/a[] / ] Out[18]= I[19]:= N[b[]/a[] / ] Out[19]= Which of these digits ca we trust? It is istructive to compute the quotiet for several idices ad see which digits remai fixed. The calculatio above strogly suggests that at least the first 30 digits are OK. The covergece is very quick because 15 terms of the expasio were take ito accout. This meas the error decays to zero i speed O( 16 ). More terms will lead to eve faster covergece. The figure below illustrates how for the preset example the umber of correct decimal digits (vertical axis) grows with (horizotal axis) whe 15 terms of the expasio are used (left) ad whe 30 terms of the expasio are used (right) As ca be see from this example doublig the umber of terms i the expasio teds to double the the accuracy of the estimate while doublig the evaluatio idex will usually ot double the accuracy. 3.. Several Domiat Terms. Whe o sigle term domiates all others the several costats have to be determied. Two terms are of the same growth for example whe the expoetial parts have the same absolute value (like ad ( ) ) or whe they differ by a polyomial multiple (like ad 3 ; although the latter grows more quickly tha the former both terms have to be take ito accout because the higher order terms of their
6 6 MANUEL KAUERS expasios may iterfere with each other). Suppose that is such that (b (1) )... (b (m) such that lim b (1) k b (i) a c 1 b (1) ) are such that b(i) b (j) + + c rb (r) ( ) = O( k ) for all 1 ij m ad some k ad = 0 for i = m r ad all k. The the first m terms cotribute sigificatly to the asymptotic behavior of (a ) ad the cotributio of the remaiig terms ca be eglected: a c 1 b (1) + + c m b (m) ( ) We have to determie approximatios for the costats c 1... c m. This ca be doe by solvig a suitable system of liear equatios as illustrated i the followig typical example sessio. I this example there are two domiat terms: the secod ad the third. For both of them we compute the first 15 terms of the expasio. The we set up a liear system for the desired costats c c 3 ad solve it. I[0]:= ClearAll[a b]; I[1]:= a[0] = ; a[1] = 1; a[] = 7; I[]:= a[ Iteger] := a[] = (( )a[ 3 + ] + ( )a[ + ] + ( )a[ 1 + ])/( ); I[3]:= a[] Out[3]= I[4]:= rec = 4(1+ )f[]+( 3+0 )f[1+] 3( 3+11 )f[+]+( 1+18 )f[3+]; I[5]:= Asymptotics[rec f[]] Out[5]= { ( 1 ) ( 3 ) (3 71)/6 ( 3 ) (3+ 71)/6 I[6]:= terms = Last[FurtherTerms[rec f[] {0 3 { {115]]; I[7]:= b[ Iteger] := ( 3 ) (3 71)/6 Sum[terms[[k + 1]] k {k015]; I[8]:= terms3 = Last[FurtherTerms[rec f[] {0 3 { {115]]; I[9]:= b3[ Iteger] := ( 3 ) (3+ 71)/6 Sum[terms3[[k + 1]] k {k0 15]; I[30]:= $RecursioLimit = 5 ; I[31]:= A = N[{b[]/a[]b3[]/a[] /. {{ 00 { 00050]; I[3]:= LiearSolve[A {1 1] Out[3]= { I[33]:= A = N[{b[]/a[]b3[]/a[] /. {{ 000 { ]; I[34]:= LiearSolve[A {1 1] Out[34]= { The umber of correct digits is ow less tha before because some accuracy is lost durig the liear system solvig. But still the umber of correct digits ca be icreased by takig ito accout more terms of the expasio. The figure blow shows for this example the umber of correct digits i depedece of the evaluatio idex ( + 00) whe 15 (left) or 30 (right) terms of the expasios are used.
7 ASYMPTOTIC EXPANSIONS Whe you eed more terms tha you ca compute. The best way of gettig a accurate estimate for the multiplicative costat(s) is to use as may terms of the expasio as possible. Sometimes oe would like to use eve more terms tha ca be computed explicitly with reasoable effort. For this situatio there exists a simple way to use higher order terms without eve kowig them explicitly. This is kow as Richardso s covergece acceleratio [4]. The idea is that if ( a c 1 + β k k + β ) k+1 k+1 + ( ) the ad therefore ( a c 1 + β k k k + β ) k+1 k+1 k+1 + k a a k 1 ( c ) k + O( (k+1) ) ( ) ( ). This meas that ( k a a ) coverges to the same limit as (a k 1 ) but oe order of magitude faster. Of course the scheme ca be iterated such as to elimiate several terms at oce. Here is some Mathematica code for doig this. I[35]:= Richardso[expr k Iteger] := Together[ k (expr /. ) expr ]; I[36]:= Richardso[expr {k0 Iteger k1 Iteger] := Fold[Richardso[#1 #]& expr Rage[k0 k1]] k 1 I the followig example we redo the calculatio of Sectio 3.1 by computig oly 1 terms of the expasio explicitly ad elimiatig three more terms with the commad just defied. I[37]:= ClearAll[a b]; I[38]:= a[0] = 1; a[1] = ; a[ Iteger] := a[] = ((3( )+)a[ 1] a[ ])/( 1); I[39]:= terms = Last[FurtherTerms[( + 1)f[ + ] (3 + )f[ + 1] + ( + )f[] f[] { { 1 5 {11]]; I[40]:= b[ Iteger] := ( 3+ ) 5 (1 5)/ Sum[terms[[k + 1]] k {k01]; I[41]:= $RecursioLimit = 5 ; I[4]:= u[ Iteger] = Richardso[b[]/a[] {13 15]; I[43]:= N[u[500] 50] Out[43]= Possible Issues There seems to be a bug related to the costructio of logarithmic terms. Whe a result is retured it seems correct but i some istaces the computatio aborts with a error whe it should ot. This bug will be fixed i a future versio.
8 8 MANUEL KAUERS I some examples Mathematica has trouble hadlig algebraic umbers. Typically these troubles become more likely for higher order terms. I such situatios it ca help to compute oly low order terms first the to rephrase all Root expressios i terms of AlgebraicNumber expressios of a commo umber field ad the to apply FurtherTerms to this. This is ot oly more stable but also more efficiet. Please report other problems to mkauers@risc.jku.at Refereces [1] Cyril Baderier Felix Cher ad Hsie-Kuei Hwag. Asymptotics of D-fiite sequeces. i preparatio. [] G.D. Birkhoff ad W.J. Trjitzisky. Aalytic theory of sigular differece equatios. Acta Mathematica 60: [3] Philippe Flajolet ad Robert Sedgewick. Aalytic Combiatorics. Cambridge Uiversity Press 009. [4] Lewis Fry Richardso. The approximate arithmetic solutio by fiite differece of physical problems ivolvig differetial equatios with a applicatio to the stress i a masory dam. Philosophical Trasactios of the Royal Society Lodo Series A : [5] Bruo Salvy. Examples of automatic asymptotic expasios. SIGSAM Bulleti 5() [6] Jet Wimp ad Doro Zeilberger. Resurrectig the asymptotics of liear recurreces. Joural of Mathematical Aalysis ad Applicatios 111: [7] Doro Zeilberger. Asyrec: A Maple package for computig the asymptotics of solutios of liear recurrece equatios with polyomial coefficiets. The Persoal Joural of Shalsoh B. Ekhad ad Doro Zeilberger 008. Mauel Kauers Research Istitute for Symbolic Computatio J. Kepler Uiversity Liz Austria address: mkauers@risc.ui-liz.ac.at
Math 155 (Lecture 3)
Math 55 (Lecture 3) September 8, I this lecture, we ll cosider the aswer to oe of the most basic coutig problems i combiatorics Questio How may ways are there to choose a -elemet subset of the set {,,,
More informationRecursive Algorithms. Recurrences. Recursive Algorithms Analysis
Recursive Algorithms Recurreces Computer Sciece & Egieerig 35: Discrete Mathematics Christopher M Bourke cbourke@cseuledu A recursive algorithm is oe i which objects are defied i terms of other objects
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
More information1 Generating functions for balls in boxes
Math 566 Fall 05 Some otes o geeratig fuctios Give a sequece a 0, a, a,..., a,..., a geeratig fuctio some way of represetig the sequece as a fuctio. There are may ways to do this, with the most commo ways
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationZeros of Polynomials
Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationarxiv: v1 [cs.sc] 2 Jan 2018
Computig the Iverse Melli Trasform of Holoomic Sequeces usig Kovacic s Algorithm arxiv:8.9v [cs.sc] 2 Ja 28 Research Istitute for Symbolic Computatio RISC) Johaes Kepler Uiversity Liz, Alteberger Straße
More informationNEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE
UPB Sci Bull, Series A, Vol 79, Iss, 207 ISSN 22-7027 NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE Gabriel Bercu We itroduce two ew sequeces of Euler-Mascheroi type which have fast covergece
More informationRecurrence Relations
Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationMath 2784 (or 2794W) University of Connecticut
ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationThe picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled
1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how
More informationCurve Sketching Handout #5 Topic Interpretation Rational Functions
Curve Sketchig Hadout #5 Topic Iterpretatio Ratioal Fuctios A ratioal fuctio is a fuctio f that is a quotiet of two polyomials. I other words, p ( ) ( ) f is a ratioal fuctio if p ( ) ad q ( ) are polyomials
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationTopic 1 2: Sequences and Series. A sequence is an ordered list of numbers, e.g. 1, 2, 4, 8, 16, or
Topic : Sequeces ad Series A sequece is a ordered list of umbers, e.g.,,, 8, 6, or,,,.... A series is a sum of the terms of a sequece, e.g. + + + 8 + 6 + or... Sigma Notatio b The otatio f ( k) is shorthad
More informationSection 11.8: Power Series
Sectio 11.8: Power Series 1. Power Series I this sectio, we cosider geeralizig the cocept of a series. Recall that a series is a ifiite sum of umbers a. We ca talk about whether or ot it coverges ad i
More informationPolynomials with Rational Roots that Differ by a Non-zero Constant. Generalities
Polyomials with Ratioal Roots that Differ by a No-zero Costat Philip Gibbs The problem of fidig two polyomials P(x) ad Q(x) of a give degree i a sigle variable x that have all ratioal roots ad differ by
More informationA Note on the Symmetric Powers of the Standard Representation of S n
A Note o the Symmetric Powers of the Stadard Represetatio of S David Savitt 1 Departmet of Mathematics, Harvard Uiversity Cambridge, MA 0138, USA dsavitt@mathharvardedu Richard P Staley Departmet of Mathematics,
More information3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials
Math 60 www.timetodare.com 3. Properties of Divisio 3.3 Zeros of Polyomials 3.4 Complex ad Ratioal Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered
More informationPhysics 116A Solutions to Homework Set #1 Winter Boas, problem Use equation 1.8 to find a fraction describing
Physics 6A Solutios to Homework Set # Witer 0. Boas, problem. 8 Use equatio.8 to fid a fractio describig 0.694444444... Start with the formula S = a, ad otice that we ca remove ay umber of r fiite decimals
More informationThe Random Walk For Dummies
The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oe-dimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli
More information(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3
MATH 337 Sequeces Dr. Neal, WKU Let X be a metric space with distace fuctio d. We shall defie the geeral cocept of sequece ad limit i a metric space, the apply the results i particular to some special
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More informationSeries: Infinite Sums
Series: Ifiite Sums Series are a way to mae sese of certai types of ifiitely log sums. We will eed to be able to do this if we are to attai our goal of approximatig trascedetal fuctios by usig ifiite degree
More informationWe are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n
Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at
More informationREGRESSION (Physics 1210 Notes, Partial Modified Appendix A)
REGRESSION (Physics 0 Notes, Partial Modified Appedix A) HOW TO PERFORM A LINEAR REGRESSION Cosider the followig data poits ad their graph (Table I ad Figure ): X Y 0 3 5 3 7 4 9 5 Table : Example Data
More informationCALCULATION OF FIBONACCI VECTORS
CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College
More informationx a x a Lecture 2 Series (See Chapter 1 in Boas)
Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio
More informationMa 530 Introduction to Power Series
Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power
More informationComplex Analysis Spring 2001 Homework I Solution
Complex Aalysis Sprig 2001 Homework I Solutio 1. Coway, Chapter 1, sectio 3, problem 3. Describe the set of poits satisfyig the equatio z a z + a = 2c, where c > 0 ad a R. To begi, we see from the triagle
More informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More information1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.
.3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(
More informationDynamic Programming. Sequence Of Decisions
Dyamic Programmig Sequece of decisios. Problem state. Priciple of optimality. Dyamic Programmig Recurrece Equatios. Solutio of recurrece equatios. Sequece Of Decisios As i the greedy method, the solutio
More informationDynamic Programming. Sequence Of Decisions. 0/1 Knapsack Problem. Sequence Of Decisions
Dyamic Programmig Sequece Of Decisios Sequece of decisios. Problem state. Priciple of optimality. Dyamic Programmig Recurrece Equatios. Solutio of recurrece equatios. As i the greedy method, the solutio
More informationMath 312 Lecture Notes One Dimensional Maps
Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,
More informationSECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES
SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES Read Sectio 1.5 (pages 5 9) Overview I Sectio 1.5 we lear to work with summatio otatio ad formulas. We will also itroduce a brief overview of sequeces,
More informationProblem Set 2 Solutions
CS271 Radomess & Computatio, Sprig 2018 Problem Set 2 Solutios Poit totals are i the margi; the maximum total umber of poits was 52. 1. Probabilistic method for domiatig sets 6pts Pick a radom subset S
More informationRandom Matrices with Blocks of Intermediate Scale Strongly Correlated Band Matrices
Radom Matrices with Blocks of Itermediate Scale Strogly Correlated Bad Matrices Jiayi Tog Advisor: Dr. Todd Kemp May 30, 07 Departmet of Mathematics Uiversity of Califoria, Sa Diego Cotets Itroductio Notatio
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationThe Rand and block distances of pairs of set partitions
The Rad ad block distaces of pairs of set partitios Frak Ruskey 1 ad Jeifer Woodcock 1 Dept. of Computer Sciece, Uiversity of Victoria, CANADA Abstract. The Rad distaces of two set partitios is the umber
More informationCSE 1400 Applied Discrete Mathematics Number Theory and Proofs
CSE 1400 Applied Discrete Mathematics Number Theory ad Proofs Departmet of Computer Scieces College of Egieerig Florida Tech Sprig 01 Problems for Number Theory Backgroud Number theory is the brach of
More informationMAT1026 Calculus II Basic Convergence Tests for Series
MAT026 Calculus II Basic Covergece Tests for Series Egi MERMUT 202.03.08 Dokuz Eylül Uiversity Faculty of Sciece Departmet of Mathematics İzmir/TURKEY Cotets Mootoe Covergece Theorem 2 2 Series of Real
More informationLecture 3. Digital Signal Processing. Chapter 3. z-transforms. Mikael Swartling Nedelko Grbic Bengt Mandersson. rev. 2016
Lecture 3 Digital Sigal Processig Chapter 3 z-trasforms Mikael Swartlig Nedelko Grbic Begt Madersso rev. 06 Departmet of Electrical ad Iformatio Techology Lud Uiversity z-trasforms We defie the z-trasform
More informationSequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018
CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical
More information1 Last time: similar and diagonalizable matrices
Last time: similar ad diagoalizable matrices Let be a positive iteger Suppose A is a matrix, v R, ad λ R Recall that v a eigevector for A with eigevalue λ if v ad Av λv, or equivaletly if v is a ozero
More informationMAS111 Convergence and Continuity
MAS Covergece ad Cotiuity Key Objectives At the ed of the course, studets should kow the followig topics ad be able to apply the basic priciples ad theorems therei to solvig various problems cocerig covergece
More information18.440, March 9, Stirling s formula
Stirlig s formula 8.44, March 9, 9 The factorial fuctio! is importat i evaluatig biomial, hypergeometric, ad other probabilities. If is ot too large,! ca be computed directly, by calculators or computers.
More informationPolynomial identity testing and global minimum cut
CHAPTER 6 Polyomial idetity testig ad global miimum cut I this lecture we will cosider two further problems that ca be solved usig probabilistic algorithms. I the first half, we will cosider the problem
More information1 6 = 1 6 = + Factorials and Euler s Gamma function
Royal Holloway Uiversity of Lodo Departmet of Physics Factorials ad Euler s Gamma fuctio Itroductio The is a self-cotaied part of the course dealig, essetially, with the factorial fuctio ad its geeralizatio
More informationAdvanced Course of Algorithm Design and Analysis
Differet complexity measures Advaced Course of Algorithm Desig ad Aalysis Asymptotic complexity Big-Oh otatio Properties of O otatio Aalysis of simple algorithms A algorithm may may have differet executio
More informationSection 7 Fundamentals of Sequences and Series
ectio Fudametals of equeces ad eries. Defiitio ad examples of sequeces A sequece ca be thought of as a ifiite list of umbers. 0, -, -0, -, -0...,,,,,,. (iii),,,,... Defiitio: A sequece is a fuctio which
More informationCLOSED FORM FORMULA FOR THE NUMBER OF RESTRICTED COMPOSITIONS
Submitted to the Bulleti of the Australia Mathematical Society doi:10.1017/s... CLOSED FORM FORMULA FOR THE NUMBER OF RESTRICTED COMPOSITIONS GAŠPER JAKLIČ, VITO VITRIH ad EMIL ŽAGAR Abstract I this paper,
More informationSeptember 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1
September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright
More informationw (1) ˆx w (1) x (1) /ρ and w (2) ˆx w (2) x (2) /ρ.
2 5. Weighted umber of late jobs 5.1. Release dates ad due dates: maximimizig the weight of o-time jobs Oce we add release dates, miimizig the umber of late jobs becomes a sigificatly harder problem. For
More informationAdvanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology
Advaced Aalysis Mi Ya Departmet of Mathematics Hog Kog Uiversity of Sciece ad Techology September 3, 009 Cotets Limit ad Cotiuity 7 Limit of Sequece 8 Defiitio 8 Property 3 3 Ifiity ad Ifiitesimal 8 4
More informationLesson 10: Limits and Continuity
www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals
More informationCALCULUS BASIC SUMMER REVIEW
CALCULUS BASIC SUMMER REVIEW NAME rise y y y Slope of a o vertical lie: m ru Poit Slope Equatio: y y m( ) The slope is m ad a poit o your lie is, ). ( y Slope-Itercept Equatio: y m b slope= m y-itercept=
More informationUnit 6: Sequences and Series
AMHS Hoors Algebra 2 - Uit 6 Uit 6: Sequeces ad Series 26 Sequeces Defiitio: A sequece is a ordered list of umbers ad is formally defied as a fuctio whose domai is the set of positive itegers. It is commo
More informationLecture 3 The Lebesgue Integral
Lecture 3: The Lebesgue Itegral 1 of 14 Course: Theory of Probability I Term: Fall 2013 Istructor: Gorda Zitkovic Lecture 3 The Lebesgue Itegral The costructio of the itegral Uless expressly specified
More information10.2 Infinite Series Contemporary Calculus 1
10. Ifiite Series Cotemporary Calculus 1 10. INFINITE SERIES Our goal i this sectio is to add together the umbers i a sequece. Sice it would take a very log time to add together the ifiite umber of umbers,
More informationDefinition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4.
4. BASES I BAACH SPACES 39 4. BASES I BAACH SPACES Sice a Baach space X is a vector space, it must possess a Hamel, or vector space, basis, i.e., a subset {x γ } γ Γ whose fiite liear spa is all of X ad
More informationis also known as the general term of the sequence
Lesso : Sequeces ad Series Outlie Objectives: I ca determie whether a sequece has a patter. I ca determie whether a sequece ca be geeralized to fid a formula for the geeral term i the sequece. I ca determie
More informationA New Method to Order Functions by Asymptotic Growth Rates Charlie Obimbo Dept. of Computing and Information Science University of Guelph
A New Method to Order Fuctios by Asymptotic Growth Rates Charlie Obimbo Dept. of Computig ad Iformatio Sciece Uiversity of Guelph ABSTRACT A ew method is described to determie the complexity classes of
More informationAPPENDIX F Complex Numbers
APPENDIX F Complex Numbers Operatios with Complex Numbers Complex Solutios of Quadratic Equatios Polar Form of a Complex Number Powers ad Roots of Complex Numbers Operatios with Complex Numbers Some equatios
More informationTHE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.
THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of
More informationAutomated Proofs for Some Stirling Number Identities
Autoated Proofs for Soe Stirlig Nuber Idetities Mauel Kauers ad Carste Scheider Research Istitute for Sybolic Coputatio Johaes Kepler Uiversity Altebergerstraße 69 A4040 Liz, Austria Subitted: Sep 1, 2007;
More informationAnalysis of Algorithms. Introduction. Contents
Itroductio The focus of this module is mathematical aspects of algorithms. Our mai focus is aalysis of algorithms, which meas evaluatig efficiecy of algorithms by aalytical ad mathematical methods. We
More informationNotes on iteration and Newton s method. Iteration
Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f
More information5 Sequences and Series
Bria E. Veitch 5 Sequeces ad Series 5. Sequeces A sequece is a list of umbers i a defiite order. a is the first term a 2 is the secod term a is the -th term The sequece {a, a 2, a 3,..., a,..., } is a
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationSection 6.4: Series. Section 6.4 Series 413
ectio 64 eries 4 ectio 64: eries A couple decides to start a college fud for their daughter They pla to ivest $50 i the fud each moth The fud pays 6% aual iterest, compouded mothly How much moey will they
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationSection A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics
X0/70 NATIONAL QUALIFICATIONS 005 MONDAY, MAY.00 PM 4.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationSummary: Congruences. j=1. 1 Here we use the Mathematica syntax for the function. In Maple worksheets, the function
Summary: Cogrueces j whe divided by, ad determiig the additive order of a iteger mod. As described i the Prelab sectio, cogrueces ca be thought of i terms of coutig with rows, ad for some questios this
More informationChapter 7: Numerical Series
Chapter 7: Numerical Series Chapter 7 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals
More informationChapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics:
Chapter 6 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals (which is what most studets
More information6.003 Homework #3 Solutions
6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the
More informationCHAPTER 1 SEQUENCES AND INFINITE SERIES
CHAPTER SEQUENCES AND INFINITE SERIES SEQUENCES AND INFINITE SERIES (0 meetigs) Sequeces ad limit of a sequece Mootoic ad bouded sequece Ifiite series of costat terms Ifiite series of positive terms Alteratig
More informationRoberto s Notes on Series Chapter 2: Convergence tests Section 7. Alternating series
Roberto s Notes o Series Chapter 2: Covergece tests Sectio 7 Alteratig series What you eed to kow already: All basic covergece tests for evetually positive series. What you ca lear here: A test for series
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS
MASSACHUSTTS INSTITUT OF TCHNOLOGY 6.436J/5.085J Fall 2008 Lecture 9 /7/2008 LAWS OF LARG NUMBRS II Cotets. The strog law of large umbers 2. The Cheroff boud TH STRONG LAW OF LARG NUMBRS While the weak
More informationFall 2013 MTH431/531 Real analysis Section Notes
Fall 013 MTH431/531 Real aalysis Sectio 8.1-8. Notes Yi Su 013.11.1 1. Defiitio of uiform covergece. We look at a sequece of fuctios f (x) ad study the coverget property. Notice we have two parameters
More informationIntroduction to Signals and Systems, Part V: Lecture Summary
EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive
More informationREGULARIZATION OF CERTAIN DIVERGENT SERIES OF POLYNOMIALS
REGULARIZATION OF CERTAIN DIVERGENT SERIES OF POLYNOMIALS LIVIU I. NICOLAESCU ABSTRACT. We ivestigate the geeralized covergece ad sums of series of the form P at P (x, where P R[x], a R,, ad T : R[x] R[x]
More informationOn Generalized Fibonacci Numbers
Applied Mathematical Scieces, Vol. 9, 215, o. 73, 3611-3622 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.5299 O Geeralized Fiboacci Numbers Jerico B. Bacai ad Julius Fergy T. Rabago Departmet
More informationNICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =
AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,
More informationThe Growth of Functions. Theoretical Supplement
The Growth of Fuctios Theoretical Supplemet The Triagle Iequality The triagle iequality is a algebraic tool that is ofte useful i maipulatig absolute values of fuctios. The triagle iequality says that
More informationChapter 3. Strong convergence. 3.1 Definition of almost sure convergence
Chapter 3 Strog covergece As poited out i the Chapter 2, there are multiple ways to defie the otio of covergece of a sequece of radom variables. That chapter defied covergece i probability, covergece i
More informationSection 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations
Differece Equatios to Differetial Equatios Sectio. Calculus: Areas Ad Tagets The study of calculus begis with questios about chage. What happes to the velocity of a swigig pedulum as its positio chages?
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationMATH 21 SECTION NOTES
MATH SECTION NOTES EVAN WARNER. March 9.. Admiistrative miscellay. These weekly sectios will be for some review ad may example problems, i geeral. Attedace will be take as per class policy. We will be
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationInfinite Series and Improper Integrals
8 Special Fuctios Ifiite Series ad Improper Itegrals Ifiite series are importat i almost all areas of mathematics ad egieerig I additio to umerous other uses, they are used to defie certai fuctios ad to
More informationHarmonic Number Identities Via Euler s Transform
1 2 3 47 6 23 11 Joural of Iteger Sequeces, Vol. 12 2009), Article 09.6.1 Harmoic Number Idetities Via Euler s Trasform Khristo N. Boyadzhiev Departmet of Mathematics Ohio Norther Uiversity Ada, Ohio 45810
More informationTEACHER CERTIFICATION STUDY GUIDE
COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra
More information