Bernoulli, Ramanujan, Toeplitz e le matrici triangolari
|
|
- Derrick Bailey
- 6 years ago
- Views:
Transcription
1 Due Giori di Algebra Lieare Numerica dibeede/gg/home.html Geova, 6 7 Febbraio Beroulli, Ramauja, Toeplitz e le matrici triagolari Carmie Di Fiore, Fracesco Tudisco, Paolo Zellii Speaker: Carmie Di Fiore By usig oe of the defiitios of the Beroulli umbers, we observe that they solve particular odd ad eve lower triagular Toeplitz (l.t.t.) systems. I a paper Ramauja writes dow a sparse lower triagular system solved by Beroulli umbers; we observe that such system is equivalet to a sparse l.t.t. system. The attempt to obtai the sparse l.t.t. Ramauja system from the l.t.t. odd ad eve systems, leads us to study efficiet methods for solvig geeric l.t.t. systems.
2 Beroulli umbers are the ratioal umbers satisfyig the followig idetity t + e t So, they satisfy the followig liear equatios j ] [ j + B () t +! t + B k () (k)! tk. k ( ) j B k k (), j, 3, 4,..., k ( ) j eve : ( ) ( ) 4 4 ( ) ( ) ( ) ( ) ( ) ( 4 ) ( ) ( ) j odd : ( ) ( ) 3 3 ( ) ( ) ( ) ( ) ( ) ( 4 ) ( ) B () B () B 4() B 6() B () B () B 4() B 6() I other words, the Beroulli umbers ca be obtaied by solvig (by forward substitutio) a lower triagular liear system (oe of the above two). For example, by forward solvig the first system, I have obtaied the first Beroulli umbers: 3 4 3/ 5/ 7/ B (), B () 6, B 4() 3, B 6() 4, B 8() 3, B () 5 66, B () 69 73, B 4() 7 6.6, B 6() ,... Beroulli umbers appear i the Euler-Maclauri summatio formula, ad, i particular, i the expressio of the error of the trapezoidal quadrature rule as sum of eve powers of the itegratio step h (the expressio that justifies the efficiecy of the Romberg-Trapezoidal quadrature method). Beroulli umbers are also ofte ivolved whe studyig the Riema-Zeta fuctio. For example, well kow is the followig Euler formula:,. ζ(s) + k k s, ζ() 4 B () π,,, 3,... ()! (see also [Riema s Zeta Fuctio, H. M. Edwards, 974]). The Ramauja s paper we refer i the followig is etitled Some properties of Beroulli s umbers (9).
3 The coefficiet matrices of the previous two lower triagular liear systems are submatrices of the matrix X displayed here below: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) X ( ) ( ) ( ) ( 3 ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) Oe ca easily observe that X ca be rewritte as a power series: Proof: [X] ij X + k (i j)! [Y i j ] ij 3 56 k! Y k, Y 3 4 j (i )(i ) (i j)! This remark is the startig poit i order to show that ( ) where φ ( ) i, j i. j ( ) ( ) ( ) ( ) ( ) (k + )! φk 6 6 6, k 4 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) (k + )! φk 5 5 5, k 4 ( ) ( ) ( ) ( ) ,, 3 4, 3 5 6,.... 3
4 It follows that the liear systems solved by the Beroulli umbers, ca be rewritte as follows, i terms of the matrix φ: / / B () + B () / /6 (k + )! φk B 4 () k B 6 () 3/3 4/56 / /4 : q e, (almosteve) 5/9 /8 + k (k + )! φk B () B () B 4 () B 6 () / (3/)/3 (5/)/5 (7/)/7 (9/)/9 / / / / Now we trasform φ ito a Toeplitz matrix. We have that d d DφD d 3 d d d d3 d xz, Z 3 d4 d 3 56 d5 d 4 : q o. (almostodd) d d d 3 d 4, iff d k xk d, k,, 3,..., iff (k )! D d D x, D x x! x 4! x ( )!. We are ready to itroduce the two eve ad odd lower triagular Toeplitz (l.t.t.) systems solved by the Beroulli umbers. Set B () b B () B 4 () where the B i (), i,,,..., are the Beroulli umbers. 4
5 The the (almosteve) system + k (k+)! φk b q e is equivalet to the system + k (k+)! (D xφd x ) k (D x b) D x q e, i.e. to the followig l.t.t. eve system: + k Idem, the (almostodd) system + k D x q o, i.e. to the followig l.t.t. odd system: x k (k + )! Zk (D x b) D x q e. (eve) (k+)! φk b q o is equivalet to the system + k + k (k+)! (D xφd x ) k (D x b) x k (k + )! Zk (D x b) D x q o. (odd) So, Beroulli umbers ca be computed by usig a l.t.t. liear system solver. Such solver yields the followig vector z: B () x z D x b! B () x 4! B 4() x s (s)! B s() x ( )! B () from which oe obtais the vector of the first Beroulli umbers: {b} {D x z}. Why x positive differet from may be useful A suitable choice of x ca make possible ad more stable the computatio via a l.t.t. solver of the etries z i of z for very large i. I fact, sice i i, x i (i)! B i() ( ) i+ p i, p i xi (i)! 4 πi (πe) i, p i+ x p i 4π, we have that xi (i)! B i() (+ ) if x < 4π (x > 4π ), both bad situatios. Istead, for x 4π the sequece xi (i)! B i(), i,,,..., should be lower ad upper bouded.... x (4)! B 4() iff x 6.84 x4 (8)! B 8() iff x 33. x8 (6)! B 6() iff x 36. x6 (3)! B 3() about iff x 6 (8.54) iff x 37.8 xs (s)! B s() about iff xs (s)! 4 πs ss (πe) iff x s (s)! (πe) s s s s 4... πs More geerally, the parameter x should be used to make more stable the l.t.t. solver. 5
6 Ramauja i a paper states that the Beroulli umbers B (), B 4 (), B 6 (),..., satisfy the followig lower triagular sparse system of liear equatios: B () B 4 () B 6 () B 8 () B () B () B 4 () B 6 () B 8 () B () B () /6 /3 /4 /45 /3 4/455 / /36 3/665 /3 /55 (actually, sice Ramauja-Beroulli umbers are the moduli of ours, his equatios, obtaied by a aalytical proof, are a bit differet). So, for example, from the above equatios I have easily computed the Beroulli umbers B 8 (), B (), B () (from the oes already computed): B 8 () , B () , B () Problem. Is it possible to obtai such Ramauja lower triagular sparse system of equatios from our odd ad eve l.t.t. liear systems Is it possible to obtai other sparse equatios (hopefully more sparse tha Ramauja oes) defiig the Beroulli umbers Note that the Ramauja matrix, say R, has ozero etries exactly i the places where the followig Toeplitz matrix + γ k (Z 3 ) k, Z k has ozero etries. But R it is ot a Toeplitz matrix!... break for some moths... Note that the vector of the idetermiates i the Ramauja system is Z T times our vector b: B () B () B 4 () B 6 () B () B 4 () ZT b. So, the Ramauja system ca be rewritte as follows R(Z T b) f, f [ f f f 3 ] T. 6
7 Remark The Ramauja matrix R satisfies the followig idetity ivolvig a sparse lower triagular Toeplitz matrix R: R! x 4! x 6! x 3! x 4! x 6! x 3 R, R + s x 3s (6s + )!(s + ) Z3s. RZ T b f, f [f f f 3 ] T iff R! x R! x x! 4! x 6! x 3 4! x 6! x 3 x 4! x 3 6! R x! x! x 4! x 4! B () B 4 () B 6 () x 3 6! x 3 6! x! ZT b f ZT b f So, the Ramauja system is is equivalet to the followig sparse l.t.t. system: R L(a R ) R (Z T D x b) x! x 4! x 3 6! f f f 3, x 4! where x 3 6! iff iff 8!3 x3 8!3 x3 8!3 x3 4!5 x6 8!3 x3 4!5 x6 8!3 x3 4!5 x6 8!3 x3!7 x9 4!5 x6 8!3 x3!7 x9 4!5 x6 8!3 x3 Note the ew otatio : a a a a L(a) f f f 3. a a a a a a, a R. 8!3 x3 4!5 x6!7 x9. 7
8 THEOREM Notatios: Z is the lower shift matrix Z L(a) is the lower triagular Toeplitz matrix with first colum a, i.e. a a a a + a, L(a) a a a i Z i a a a i a 3 a a a d(z) is the diagoal matrix with z i as diagoal etries. Set b B () B () B 4 (), where B i (), i,,,..., are the Beroulli umbers.,, D x diag ( xi, i,,,...), x R, (i)! The the vectors D x b ad Z T D x b solve the followig l.t.t. liear systems L(a) (D x b) D x q, L(a) (Z T D x b) d(z)z T D x q, where the vectors a (a i ) + i, q (q i) + i, ad z (z i) + i, ca assume respectively the values: a R x i i δ i mod 3 (i + )!( 3 i + ), qr i (i + )(i + ) ( δ 3 i mod ), i,,, 3,... 3 zi R δ i mod 3 i,, 3,..., 3i +, a o i a e i xi (i + )!, qe i z e i, i,,, 3,... i + i, i,, 3,..., i + x i (i + )!, i,,, 3,..., qo, qi o, i,, 3,... zi o i, i,, 3,.... i + 8
9 Problem (regardig the computatio of the Beroulli umbers). Ca the Ramauja l.t.t. sparse system L(a R )D x b D x q R, be obtaied as a cosequece of the eve ad odd l.t.t. system L(a e )D x b D x q e, L(a o )D x b D x q o fid â e, â o such that L(â e )L(a e ) L(a () ) L(â o )L(a o ), i.e. such that L(a e )â e a () L(a o )â o, with a () a R or a() more sparse tha a R. Importat: the computatio of such vectors â e, â o ad a () should be cheaper tha solvig the origial eve ad odd (dese) systems. A more geeral problem: is it possible to trasform efficietly a geeric l.t.t. matrix ito a more sparse l.t.t. matrix Questio: give a i, i,, 3,..., is it possible to obtai cheaply â i ad a () i a a a a 3 a a a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a â â â 3 â 4 â 5 â â â 3 â â â 3 â 4 â 5 a () a () a () ; a () ; γ, R b such that If the aswer is yes, the a dese l.t.t system ca be trasformed efficietly ito a sparse l.t.t. system: L(a)z c iff L(â)L(a)z L(â)c iff L(γ)z L(â)c. 9
10 DEFINITIONS: u E u u LEMMA: If u a a a, Eu u u PROBLEM: Give a, v u u C+, L(a), b, E s, b, L(Eu) v v, the a a a,, bs, I u T u I I u T u T L(Eu)Ev EL(u)L(v), L(E s u)e s v E s L(u)v, s N. a a, fid â L(a)â Ea () â â a () ad a() a () a () such that, b, L(a)L(â) L(Ea () )., b. Questios: Is it possible to obtai cheaply â i ad a () i There exist explicit formulas for the â i ad a () i At the momet, let us see i detail, with two examples, how the solutios of the above Problem ca lead to efficiet methods for solvig geeric l.t.t. liear systems.
11 EXAMPLE: 8 ( b k, b, k 3) a a () a a a 3 a 4 a 5 a 6 a 7, L(a) a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a, E. step : From a a () fid â â () such that L(a)â a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a â â â 3 â 4 â 5 â 6 â 7 : Ea () a () a () a () 3, L(a)L(â) L(Ea () ); step : From a () fid â () such that L(Ea () )Eâ () a () a () a () a () a () a () a () 3 a () a () a () 3 a () a () â () â () â () 3 : E a () a (), L(Ea () )L(Eâ () ) L(E a () ); step 3 log 8: From a () fid â () such that L(E a () )E â () a () a () a () a () â () : E 3 a (3), L(E a () )L(E â () ) L(E 3 a (3) ). L(E â () )L(Eâ () )L(â) [ L(a) ] L(E 3 a (3) ) [ I8 O ], so, oe realizes that we have performed a kid of Gaussia elimiatio.
12 How may arithmetic operatios (a.o.) Actually, give a i a () i, i,..., 7, we have to compute â i â () i, a () i â () i, a () a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a â () a () a () a () a () 3 a () a () a () â () â () â () â () 3 â â â 3 â 4 â 5 â 6 â 7 a (), ϕ a.o.. a () a () a () 3, ϕ 4 a.o.,, ϕ 8 a.o., The geeral case: k Give a i a () i, i,..., k, we have to compute â (j) i, a (j+) i a (j) a (j) a (j) a (j) j a (j) a (j) }{{} j j j,,..., k, k (j k : oly â (j) i ) â (j) â (j) â (j) j a (j+) a (j+) j+, ϕ j a.o., Total cost: k j ϕ j Remark. Note that, at step j, the a (j+) i are the j+ ozero etries of a matrix j j ( k j k j ) l.t.t. by vector product. So, if we assume such matrix by vector product computable i at most c k j (k j) a.o., for some costat c, the k ϕ j k c j j k k j (k j) + j j CostCompOf ( â (j) i ) k ( O( k k) + CostCompOf â (j) i, i,..., ) j
13 Computig the first colum of L(a) ( 3 8) Let v [v v v ] T be ay vector. From the idetity L(E â () )L(Eâ () )L(â) [ L(a) ] L(E 3 a (3) ) [ I8 O ] ad from the Lemma, it follows that L(a)z E v iff L(E 3 a (3) )z L(â)L(Eâ () )L(E â () )E v L(â)L(Eâ () )E L(â () )v L(â)EL(â () )EL(â () )v. So, the system L(a)z E v is equivalet to the system [ ] [ ] I8 O {z}8 L(E 3 a (3) )z L(â)EL(â () )EL(â () )v {z} 8 {L(â)} 8 {E} 8 {L(â () )} 8 {E} 8 {L(â () )} 8 {v} 8 {L(â)} 8 {E} 8,4 {L(â () )} 4 {E} 4, {L(â () )} {v}. Thus the vector z z â z â â z 3 z 4 â 3 â â â 4 â 3 â â z 5 â 5 â 4 â 3 â â z 6 â 6 â 5 â 4 â 3 â â z 7 â 7 â 6 â 5 â 4 â 3 â â â () â () â () â () 3 â () â () [ â () ] [ v v ] is such that How may arithmetic operatios (a.o.) Case k a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a z z z z 3 z 4 z 5 z 6 z 7 It is clear that the above procedure requires the computatio of matrix j j l.t.t. by vector products, with j,..., k (the vectors are sparse for j,..., k). So, if we assume such matrix by vector product computable i at most c j j a.o., for some costat c, the the above procedure requires at most arithmetic operatios. c k j j O( k k) j v v. 3
14 EXAMPLE: 9 ( b k, b 3, k ) a a () a a a 3 a 4 a 5 a 6 a 7 a 8, L(a) a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a a 8 a 7 a 6 a 5 a 4 a 3 a a, E. step : From a a () fid â â () such that L(a)â a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a a 8 a 7 a 6 a 5 a 4 a 3 a a â â â 3 â 4 â 5 â 6 â 7 â 8 : Ea () a () a (), L(a)L(â) L(Ea () ); step log 3 9: From a () fid â () such that L(Ea () )Eâ () a () a () a () a () a () a () a () a () a () â () â () : E a (), L(Ea () )L(Eâ () ) L(E a () ); L(Eâ () )L(â) [ L(a) ] L(E a () ) [ I9 O ], so, oe realizes that we have performed a kid of Gaussia elimiatio. 4
15 How may arithmetic operatios (a.o.) Actually, give a i a () i, i,..., 8, we have to compute â i â () i, a () i a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a a 8 a 7 a 6 a 5 a 4 a 3 a a â () i a () a () a () â () â (), â â â 3 â 4 â 5 â 6 â 7 â 7 â 8 ϕ 3 a.o.. a () a (), ϕ 9 a.o., The geeral case: 3 k Give a i a () i, i,..., 3 k, we have to compute â (j) i, a (j+) i a (j) a (j) a (j) a (j) 3 j a (j) a (j) }{{} 3 j 3 j j,,..., k, k (j k : oly â (j) i ) â (j) â (j) â (j) 3 j a (j+) a (j+) 3j+, ϕ 3 j a.o., Total cost: k j ϕ 3 j Remark. Note that, at step j, the a (j+) i are the 3 j+ ozero etries of a matrix 3 j 3 j (3 k j 3 k j ) l.t.t. by vector product. So, if we assume such matrix by vector product computable i at most c 3 k j (k j) a.o., for some costat c, the k ϕ j k c 3 j j k 3 k j (k j) + j j CostCompOf ( â (j) i ) k ( O(3 k k) + CostCompOf â (j) i, i,..., ) 3 j 5
16 Computig the first colum of L(a) (case 3 9): Let v [v v v ] T be ay vector. From the idetity L(Eâ () )L(â) [ L(a) ] L(E a () ) [ I9 O ] ad from the Lemma, it follows that L(a)z Ev iff L(E a () )z L(â)L(Eâ () )Ev L(â)EL(â () )v. So, the system L(a)z Ev is equivalet to the system [ ] [ ] I9 O {z}9 L(E a () )z L(â)EL(â () )v Thus the vector z z z z 3 z 4 z 5 z 6 z 7 z 8 {z} 9 {L(â)} 9 {E} 9 {L(â () )} 9 {v} 9 {L(â)} 9 {E} 9,3 {L(â () )} 3 {v} 3. â â â â 3 â â â 4 â 3 â â â 5 â 4 â 3 â â â 6 â 5 â 4 â 3 â â â 7 â 6 â 5 â 4 â 3 â â â 8 â 7 â 6 â 5 â 4 â 3 â â â () â () â () v v v is such that How may arithmetic operatios (a.o.) Case 3 k a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a a 8 a 7 a 6 a 5 a 4 a 3 a a z z z z 3 z 4 z 5 z 6 z 7 z 8 It is clear that the above procedure requires the computatio of matrix 3 j 3 j l.t.t. by vector products, with j,..., k (the vectors are sparse for j,..., k). So, if we assume such matrix by vector product computable i at most c 3 j j a.o., for some costat c, the the above procedure requires at most arithmetic operatios. For the geeral case b k see the Appedix. c k 3 j j O(3 k k) j v v v. 6
17 PROBLEM Aswer to the quotatio marks i the followig equality: L(a)â a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a Ea (), E. - There is ot a uique aswer to the. - There exists a aswer that allows to obtai from the eve ad odd systems, a system solved by Beroulli umbers where i the coefficiet matrix ull diagoals alterate with the o ull oes. Fid it... - If there exists a aswer such that the first j etries of â ca be computed i at most O( j j) arithmetic operatios, for all j k, the we have a algorithm of complexity O( k k) for solvig geeric k k lower triagular Toeplitz systems. Here below is a aswer such that the first j etries of â ca be computed with zero arithmetic operatios: L(a)â a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a a a 3 a 4 a 5 a a a 4 a a 3 + a a 6 a a 5 + a a 4 a 3 a 8 a a 7 + a a 6 a 3 a 5 + a 4 Ea (), L(a) ( e + + i ( ) i a i e i+ ) e + + i δ i mod ( ai + i j ( ) j a j a i j ) ei+. 7
18 PROBLEM Aswer to the quotatio marks i the followig equality: a a a L(a)â a 3 a a Ea (), E a 4 a 3 a a a 5 a 4 a 3 a a. - There is ot a uique aswer to the. - There exists a aswer that allows to obtai the Ramauja system solved by Beroulli umbers as a cosequece of the odd (eve) system. Fid it... - If there exists a aswer such that the first 3 j etries of â ca be computed i at most O(3 j j) arithmetic operatios, for all j k, the we have a algorithm of complexity O(3 k k) for solvig geeric 3 k 3 k lower triagular Toeplitz systems. Here below is a aswer: L(a)â a a a a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a a 8 a 7 a 6 a 5 a 4 a 3 a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a a + a a 3 a a a 4 a a 3 + a a 5 + a a 4 a a 3 a 6 a a 5 a a 4 + a 3 a 7 a a 6 + a a 5 a 3 a 4 a 8 + a a 7 a a 6 a 3 a 5 + a 4 a 9 a a 8 a a 7 + a 3 a 6 a 4 a 5 a a a 9 + a a 8 a 3 a 7 a 4 a 6 + a 5 a + a a a a 9 a 3 a 8 + a 4 a 7 a 5 a 6 a () a () a () 3 Ea (), a () 3a 3 3a a + a 3, a () 3a 6 3a a 5 3a a 4 + 3a 3 3a a a 3 + 3a a 4 + a 3, a () 3 3a 9 3a a 8 3a a 7 + 6a 3a 6 3a a a 6 3a a 3a 5 3a a 3a 4 + 3a a 7 + 3a a 4 3a 4a 5 + 3a 5a + a 3 3,... â i i r a ra i r + δ i mod a i + 3 { s 3 i 6 s 6 i 6 a i 3 i+3 +3sa 3s i odd a i 6 i+6 +3sa 3s i eve Ca such â i, i,,..., 3 j, be computed i at most O(3 j j) arithmetic operatios, i,,, 3, 4, 5,.... 8
19 ( a 3 a 3a a 3 3 3a a 4 3a 3a 3 3 a 5 3a 3a 4 3a a 6 3a 3 3 3a 5 3a a 7 3a 4 3a 3a 6 3a 3 3 a 8 3a 5 3a 3a 7 3a 4 3a a 9 3a 6 3a 3 3 3a 8 3a 5 3a a 3a 7 3a 4 3a 3a 9 3a 6 3a 3 3 a 3a 8 3a 5 3a 3a 3a 7 3a 4 3a a 3a 9 3a 6 3a 3 3 3a 3a 8 3a 5 3a a 3 3a 3a 7 3a 4 3a a a a a a a a 3 a a a 3 a a a 4 a 3 a a a 4 a 3 a a a 5 a 4 a 3 a a a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a a 7 a 6 a 5 a 4 a 3 a a a 8 a 7 a 6 a 5 a 4 a 3 a a a 8 a 7 a 6 a 5 a 4 a 3 a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a a a a a 9 a 8 a 7 a 6 a 5 a 4 a 3 a a ) a a a 3 a 4 a 5 a 6 a 7 a 8 a 9 a a a a 3 a 4 a 5 a 6 a 7 a 8 a 9 a a a a 3 a 4 a 5 a 6 Ca the above 3 j 3 j (7 7) matrix by vector product be computed i at most O(3 j j) arithmetic operatios If yes, the we would have a method which solves 3 k 3 k lower triagular Toeplitz liear systems i at most O(3 k k) arithmetic operatios. If o, the look for aother solutio â of the system L(a)â [ ] T such that { â } is computable 3 j from { a } i at most O(3 j j) arithmetic operatios... 3 j THE END 9
20 APPENDIX Itroduce low complexity l.t.t. liear system solvers: L(a) a a a, a() : a Fid â (), a () such that L(a () )â () Ea (), b, so that Fid â (), a () such that L(a () )L(â () ) L(Ea () ). L(a () )â () Ea (), b, so that L(Ea () )Eâ () E a (), b, L(Ea () )L(Eâ () ) L(E a () ). (use Lemma). Fid â (), a (3) such that... Fid â (k ), a (k ) such that L(a () )â () Ea (3) L(E a () )E â () E 3 a (3), b, L(E a () )L(E â () ) L(E 3 a (3) ). L(a (k ) )â (k ) Ea (k ) so that, b3,, b, L(E k a (k ) )E k â (k ) E k a (k ) L(E k a (k ) )L(E k â (k ) ) L(E k a (k ) ). so that, b k, Fid â (k ), a (k) such that L(a (k ) )â (k ) Ea (k), b, so that L(E k a (k ) )E k â (k ) E k a (k), bk, L(E k a (k ) )L(E k â (k ) ) L(E k a (k) ).
21 The b k {}}{ I O [ a (k) ] L(Ek a (k) ) L(E k â (k ) )L(E k â (k ) ) L(Eâ () )L(â () )L(a () ). This implies that Moreover, if L(a () )z c iff L(E k a (k) )z L(â () )L(Eâ () ) L(E k â (k ) )L(E k â (k ) )c. c E k v v v v, b, k where v (v i ) + i is ay vector (for example v e ), the by usig the Lemma, we obtai the followig result: L(a () )z c iff I b k O [ a (k) ] z L(Ek a (k) )z L(â () )EL(â () )E EL(â (k ) )EL(â (k ) )v. I other words, the vector { z }, bk, such that { } { } L(a) z (for example { L(a) } { } e { } { } { } z L(â () ) { } L(â () ) a a b k a { } z v v v b, b k, v, v i i ), ca be represeted as follows { } { } { } { } { } E L(â () ) E L(â (k ) ) E L(â (k ) ) {v} { } { } { } { } { } E L(â () ) E L(â (k ) ) E b b, b, b b k b k, b k { L(â (k ) ) } {v} b b k FIRST: Compute the first etries of â () ad the first b etries of a() (cost ϕ b k); compute the first b etries of â () ad the first b etries of a () (cost ϕ b k );... compute the first etries of â (k ) ad the b k first etries of a (k ) (cost ϕ b k b ); compute the first etries of â (k ) (cost ϕ b k b ). Total cost of this FIRST operatio: k j ϕ. b j SECOND: To such cost add k j cost( (b j b j l.t.t ) (b j vector) ) (the vector is sparse if j,..., k; the cost for j is zero if v e ). See also the ext page.
22 Amout of operatios. I the followig b k ad b : FIRST: For j,..., k compute, by performig ϕ b j I a(j+), i.e. scalars â (j) b j+ i ad a (j+) i such that a (j) (ote that there is o a (k) i a (j) a (j) a (j) b j a (j) a (j) }{{} b j b j Case b. I this case, sice â (j) i to be computed (the a (j+) i are the b j+ â (j) â (j) â (j) b j arithmetic operatios, the vectors I b j â(j) ad a (j+) a (j+) bj+ SECOND: Compute the b b l.t.t. by vector product { L(â (k ) ) } b k vector products of type COMMENTS, j,..., k to be computed). ( ) i a (j) i, oly b j b l.t.t. by vector products, j,..., k, eed j ozero etries of the resultig vectors). v v b, ad b j b j l.t.t. by { L(â (j) ) } b }{{ j } b j b j, j k,...,,. So, i case b, we have to perform j j l.t.t. by vector products, for j,..., k, twice. If we assume the cost of a j j l.t.t. by vector product bouded by c j j (c costat), the the total cost of the above operatios is smaller tha O( k k) O( log ). As a cosequece we have obtaied, i particular, a l.t.t. liear system solver of complexity O( log ) Aalogously, for b 3, if we assume both ϕ 3 j ad the cost of a 3 j 3 j l.t.t. by vector product bouded by c3 j j, the the total cost of the above operatios is smaller tha O(3 k k) O( log 3 ).... But is ϕ 3 j bouded by c3 j j... For me: rao/ramauja/collectedidex.html
Inverse Matrix. A meaning that matrix B is an inverse of matrix A.
Iverse Matrix Two square matrices A ad B of dimesios are called iverses to oe aother if the followig holds, AB BA I (11) The otio is dual but we ofte write 1 B A meaig that matrix B is a iverse of matrix
More informationTheorem: Let A n n. In this case that A does reduce to I, we search for A 1 as the solution matrix X to the matrix equation A X = I i.e.
Theorem: Let A be a square matrix The A has a iverse matrix if ad oly if its reduced row echelo form is the idetity I this case the algorithm illustrated o the previous page will always yield the iverse
More information(3) If you replace row i of A by its sum with a multiple of another row, then the determinant is unchanged! Expand across the i th row:
Math 5-4 Tue Feb 4 Cotiue with sectio 36 Determiats The effective way to compute determiats for larger-sized matrices without lots of zeroes is to ot use the defiitio, but rather to use the followig facts,
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 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 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 information(3) If you replace row i of A by its sum with a multiple of another row, then the determinant is unchanged! Expand across the i th row:
Math 50-004 Tue Feb 4 Cotiue with sectio 36 Determiats The effective way to compute determiats for larger-sized matrices without lots of zeroes is to ot use the defiitio, but rather to use the followig
More informationProperties and Tests of Zeros of Polynomial Functions
Properties ad Tests of Zeros of Polyomial Fuctios The Remaider ad Factor Theorems: Sythetic divisio ca be used to fid the values of polyomials i a sometimes easier way tha substitutio. This is show by
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 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 informationSeries III. Chapter Alternating Series
Chapter 9 Series III With the exceptio of the Null Sequece Test, all the tests for series covergece ad divergece that we have cosidered so far have dealt oly with series of oegative terms. Series with
More informationTMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
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 informationMatrix Algebra 2.2 THE INVERSE OF A MATRIX Pearson Education, Inc.
2 Matrix Algebra 2.2 THE INVERSE OF A MATRIX MATRIX OPERATIONS A matrix A is said to be ivertible if there is a matrix C such that CA = I ad AC = I where, the idetity matrix. I = I I this case, C is a
More informationA widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α
Nice plottig of proteis: I A widely used display of protei shapes is based o the coordiates of the alpha carbos - - C α -s. The coordiates of the C α -s are coected by a cotiuous curve that roughly follows
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 informationThe Discrete Fourier Transform
The Discrete Fourier Trasform Complex Fourier Series Represetatio Recall that a Fourier series has the form a 0 + a k cos(kt) + k=1 b k si(kt) This represetatio seems a bit awkward, sice it ivolves two
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationCHAPTER I: Vector Spaces
CHAPTER I: Vector Spaces Sectio 1: Itroductio ad Examples This first chapter is largely a review of topics you probably saw i your liear algebra course. So why cover it? (1) Not everyoe remembers everythig
More informationTHE ZETA FUNCTION AND THE RIEMANN HYPOTHESIS. Contents 1. History 1
THE ZETA FUNCTION AND THE RIEMANN HYPOTHESIS VIKTOR MOROS Abstract. The zeta fuctio has bee studied for ceturies but mathematicias are still learig about it. I this paper, I will discuss some of the zeta
More informationDeterminants of order 2 and 3 were defined in Chapter 2 by the formulae (5.1)
5. Determiats 5.. Itroductio 5.2. Motivatio for the Choice of Axioms for a Determiat Fuctios 5.3. A Set of Axioms for a Determiat Fuctio 5.4. The Determiat of a Diagoal Matrix 5.5. The Determiat of a Upper
More informationChapter 8. Euler s Gamma function
Chapter 8 Euler s Gamma fuctio The Gamma fuctio plays a importat role i the fuctioal equatio for ζ(s that we will derive i the ext chapter. I the preset chapter we have collected some properties of the
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 information(bilinearity), a(u, v) M u V v V (continuity), a(v, v) m v 2 V (coercivity).
Precoditioed fiite elemets method Let V be a Hilbert space, (, ) V a ier product o V ad V the correspodig iduced orm. Let a be a coercive, cotiuous, biliear form o V, that is, a : V V R ad there exist
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 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 informationMon Feb matrix inverses. Announcements: Warm-up Exercise:
Math 225-4 Week 6 otes We will ot ecessarily fiish the material from a give day's otes o that day We may also add or subtract some material as the week progresses, but these otes represet a i-depth outlie
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 informationThe Gamma function Michael Taylor. Abstract. This material is excerpted from 18 and Appendix J of [T].
The Gamma fuctio Michael Taylor Abstract. This material is excerpted from 8 ad Appedix J of [T]. The Gamma fuctio has bee previewed i 5.7 5.8, arisig i the computatio of a atural Laplace trasform: 8. ft
More informationCHAPTER 5. Theory and Solution Using Matrix Techniques
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
More information6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.
6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio
More informationStochastic Matrices in a Finite Field
Stochastic Matrices i a Fiite Field Abstract: I this project we will explore the properties of stochastic matrices i both the real ad the fiite fields. We first explore what properties 2 2 stochastic matrices
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 information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationChapter 8. Euler s Gamma function
Chapter 8 Euler s Gamma fuctio The Gamma fuctio plays a importat role i the fuctioal equatio for ζ(s) that we will derive i the ext chapter. I the preset chapter we have collected some properties of the
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 informationBernoulli, Ramanujan, Toeplitz and the triangular matrices
Bernoulli, Ramanujan, Toeplitz and the triangular matrices Carmine Di Fiore, Francesco Tudisco, Paolo Zellini Dipartimento di Matematica, Università di Roma Tor Vergata April 4, 3 Abstract By using one
More informationDiscrete-Time Systems, LTI Systems, and Discrete-Time Convolution
EEL5: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we begi our mathematical treatmet of discrete-time s. As show i Figure, a discrete-time operates or trasforms some iput sequece x [
More informationDefinitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients.
Defiitios ad Theorems Remember the scalar form of the liear programmig problem, Miimize, Subject to, f(x) = c i x i a 1i x i = b 1 a mi x i = b m x i 0 i = 1,2,, where x are the decisio variables. c, b,
More informationApply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.
Eigevalue-Eigevector Istructor: Nam Su Wag eigemcd Ay vector i real Euclidea space of dimesio ca be uiquely epressed as a liear combiatio of liearly idepedet vectors (ie, basis) g j, j,,, α g α g α g α
More informationLecture 8: October 20, Applications of SVD: least squares approximation
Mathematical Toolkit Autum 2016 Lecturer: Madhur Tulsiai Lecture 8: October 20, 2016 1 Applicatios of SVD: least squares approximatio We discuss aother applicatio of sigular value decompositio (SVD) of
More informationa for a 1 1 matrix. a b a b 2 2 matrix: We define det ad bc 3 3 matrix: We define a a a a a a a a a a a a a a a a a a
Math S-b Lecture # Notes This wee is all about determiats We ll discuss how to defie them, how to calculate them, lear the allimportat property ow as multiliearity, ad show that a square matrix A is ivertible
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 informationf(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction.
Theorem 3.6.4. [Liouville s Theorem] Every bouded etire fuctio is costat. Proof. Let f be a etire fuctio. Suppose that there is M R such that M for ay z C. The for ay z C ad R > 0 f (z) f(w) 2πi (w z)
More informationThe Choquet Integral with Respect to Fuzzy-Valued Set Functions
The Choquet Itegral with Respect to Fuzzy-Valued Set Fuctios Weiwei Zhag Abstract The Choquet itegral with respect to real-valued oadditive set fuctios, such as siged efficiecy measures, has bee used i
More informationarxiv: v1 [math.nt] 5 Jan 2017 IBRAHIM M. ALABDULMOHSIN
FRACTIONAL PARTS AND THEIR RELATIONS TO THE VALUES OF THE RIEMANN ZETA FUNCTION arxiv:70.04883v [math.nt 5 Ja 07 IBRAHIM M. ALABDULMOHSIN Kig Abdullah Uiversity of Sciece ad Techology (KAUST, Computer,
More informationLINEAR ALGEBRA. Paul Dawkins
LINEAR ALGEBRA Paul Dawkis Table of Cotets Preface... ii Outlie... iii Systems of Equatios ad Matrices... Itroductio... Systems of Equatios... Solvig Systems of Equatios... 5 Matrices... 7 Matrix Arithmetic
More informationSequences, Sums, and Products
CSCE 222 Discrete Structures for Computig Sequeces, Sums, ad Products Dr. Philip C. Ritchey Sequeces A sequece is a fuctio from a subset of the itegers to a set S. A discrete structure used to represet
More informationa for a 1 1 matrix. a b a b 2 2 matrix: We define det ad bc 3 3 matrix: We define a a a a a a a a a a a a a a a a a a
Math E-2b Lecture #8 Notes This week is all about determiats. We ll discuss how to defie them, how to calculate them, lear the allimportat property kow as multiliearity, ad show that a square matrix A
More informationNUMERICAL METHODS FOR SOLVING EQUATIONS
Mathematics Revisio Guides Numerical Methods for Solvig Equatios Page 1 of 11 M.K. HOME TUITION Mathematics Revisio Guides Level: GCSE Higher Tier NUMERICAL METHODS FOR SOLVING EQUATIONS Versio:. Date:
More informationMath 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 information, then cv V. Differential Equations Elements of Lineaer Algebra Name: Consider the differential equation. and y2 cos( kx)
Cosider the differetial equatio y '' k y 0 has particular solutios y1 si( kx) ad y cos( kx) I geeral, ay liear combiatio of y1 ad y, cy 1 1 cy where c1, c is also a solutio to the equatio above The reaso
More informationb i u x i U a i j u x i u x j
M ath 5 2 7 Fall 2 0 0 9 L ecture 1 9 N ov. 1 6, 2 0 0 9 ) S ecod- Order Elliptic Equatios: Weak S olutios 1. Defiitios. I this ad the followig two lectures we will study the boudary value problem Here
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationREVISION SHEET FP1 (MEI) ALGEBRA. Identities In mathematics, an identity is a statement which is true for all values of the variables it contains.
The mai ideas are: Idetities REVISION SHEET FP (MEI) ALGEBRA Before the exam you should kow: If a expressio is a idetity the it is true for all values of the variable it cotais The relatioships betwee
More informationREVIEW 1, MATH n=1 is convergent. (b) Determine whether a n is convergent.
REVIEW, MATH 00. Let a = +. a) Determie whether the sequece a ) is coverget. b) Determie whether a is coverget.. Determie whether the series is coverget or diverget. If it is coverget, fid its sum. a)
More informationRADICAL EXPRESSION. If a and x are real numbers and n is a positive integer, then x is an. n th root theorems: Example 1 Simplify
Example 1 Simplify 1.2A Radical Operatios a) 4 2 b) 16 1 2 c) 16 d) 2 e) 8 1 f) 8 What is the relatioship betwee a, b, c? What is the relatioship betwee d, e, f? If x = a, the x = = th root theorems: RADICAL
More informationEssential Question How can you recognize an arithmetic sequence from its graph?
. Aalyzig Arithmetic Sequeces ad Series COMMON CORE Learig Stadards HSF-IF.A.3 HSF-BF.A. HSF-LE.A. Essetial Questio How ca you recogize a arithmetic sequece from its graph? I a arithmetic sequece, the
More informationMATH 31B: MIDTERM 2 REVIEW
MATH 3B: MIDTERM REVIEW JOE HUGHES. Evaluate x (x ) (x 3).. Partial Fractios Solutio: The umerator has degree less tha the deomiator, so we ca use partial fractios. Write x (x ) (x 3) = A x + A (x ) +
More informationx x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula
NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS. If g( is cotiuous i [a,b], te uder wat coditio te iterative (or iteratio metod = g( as a uique solutio i [a,b]? '( i [a,b].. Wat
More informationn m CHAPTER 3 RATIONAL EXPONENTS AND RADICAL FUNCTIONS 3-1 Evaluate n th Roots and Use Rational Exponents Real nth Roots of a n th Root of a
CHAPTER RATIONAL EXPONENTS AND RADICAL FUNCTIONS Big IDEAS: 1) Usig ratioal expoets ) Performig fuctio operatios ad fidig iverse fuctios ) Graphig radical fuctios ad solvig radical equatios Sectio: Essetial
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 informationIn number theory we will generally be working with integers, though occasionally fractions and irrationals will come into play.
Number Theory Math 5840 otes. Sectio 1: Axioms. I umber theory we will geerally be workig with itegers, though occasioally fractios ad irratioals will come ito play. Notatio: Z deotes the set of all itegers
More informationDiscrete Mathematics for CS Spring 2007 Luca Trevisan Lecture 22
CS 70 Discrete Mathematics for CS Sprig 2007 Luca Trevisa Lecture 22 Aother Importat Distributio The Geometric Distributio Questio: A biased coi with Heads probability p is tossed repeatedly util the first
More informationMA131 - Analysis 1. Workbook 9 Series III
MA3 - Aalysis Workbook 9 Series III Autum 004 Cotets 4.4 Series with Positive ad Negative Terms.............. 4.5 Alteratig Series.......................... 4.6 Geeral Series.............................
More informationOnce we have a sequence of numbers, the next thing to do is to sum them up. Given a sequence (a n ) n=1
. Ifiite Series Oce we have a sequece of umbers, the ext thig to do is to sum them up. Give a sequece a be a sequece: ca we give a sesible meaig to the followig expressio? a = a a a a While summig ifiitely
More informationPolynomial Functions and Their Graphs
Polyomial Fuctios ad Their Graphs I this sectio we begi the study of fuctios defied by polyomial expressios. Polyomial ad ratioal fuctios are the most commo fuctios used to model data, ad are used extesively
More informationMATH301 Real Analysis (2008 Fall) Tutorial Note #7. k=1 f k (x) converges pointwise to S(x) on E if and
MATH01 Real Aalysis (2008 Fall) Tutorial Note #7 Sequece ad Series of fuctio 1: Poitwise Covergece ad Uiform Covergece Part I: Poitwise Covergece Defiitio of poitwise covergece: A sequece of fuctios f
More informationNAME: ALGEBRA 350 BLOCK 7. Simplifying Radicals Packet PART 1: ROOTS
NAME: ALGEBRA 50 BLOCK 7 DATE: Simplifyig Radicals Packet PART 1: ROOTS READ: A square root of a umber b is a solutio of the equatio x = b. Every positive umber b has two square roots, deoted b ad b or
More informationChimica Inorganica 3
himica Iorgaica Irreducible Represetatios ad haracter Tables Rather tha usig geometrical operatios, it is ofte much more coveiet to employ a ew set of group elemets which are matrices ad to make the rule
More informationMDIV. Multiple divisor functions
MDIV. Multiple divisor fuctios The fuctios τ k For k, defie τ k ( to be the umber of (ordered factorisatios of ito k factors, i other words, the umber of ordered k-tuples (j, j 2,..., j k with j j 2...
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 information2.1. The Algebraic and Order Properties of R Definition. A binary operation on a set F is a function B : F F! F.
CHAPTER 2 The Real Numbers 2.. The Algebraic ad Order Properties of R Defiitio. A biary operatio o a set F is a fuctio B : F F! F. For the biary operatios of + ad, we replace B(a, b) by a + b ad a b, respectively.
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 informationChapter 2. Periodic points of toral. automorphisms. 2.1 General introduction
Chapter 2 Periodic poits of toral automorphisms 2.1 Geeral itroductio The automorphisms of the two-dimesioal torus are rich mathematical objects possessig iterestig geometric, algebraic, topological ad
More informationMATH 304: MIDTERM EXAM SOLUTIONS
MATH 304: MIDTERM EXAM SOLUTIONS [The problems are each worth five poits, except for problem 8, which is worth 8 poits. Thus there are 43 possible poits.] 1. Use the Euclidea algorithm to fid the greatest
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 informationREVISION SHEET FP1 (MEI) ALGEBRA. Identities In mathematics, an identity is a statement which is true for all values of the variables it contains.
the Further Mathematics etwork wwwfmetworkorguk V 07 The mai ideas are: Idetities REVISION SHEET FP (MEI) ALGEBRA Before the exam you should kow: If a expressio is a idetity the it is true for all values
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 informationGenerating Functions. 1 Operations on generating functions
Geeratig Fuctios The geeratig fuctio for a sequece a 0, a,..., a,... is defied to be the power series fx a x. 0 We say that a 0, a,... is the sequece geerated by fx ad a is the coefficiet of x. Example
More informationIn algebra one spends much time finding common denominators and thus simplifying rational expressions. For example:
74 The Method of Partial Fractios I algebra oe speds much time fidig commo deomiators ad thus simplifyig ratioal epressios For eample: + + + 6 5 + = + = = + + + + + ( )( ) 5 It may the seem odd to be watig
More informationIP Reference guide for integer programming formulations.
IP Referece guide for iteger programmig formulatios. by James B. Orli for 15.053 ad 15.058 This documet is iteded as a compact (or relatively compact) guide to the formulatio of iteger programs. For more
More informationMath 475, Problem Set #12: Answers
Math 475, Problem Set #12: Aswers A. Chapter 8, problem 12, parts (b) ad (d). (b) S # (, 2) = 2 2, sice, from amog the 2 ways of puttig elemets ito 2 distiguishable boxes, exactly 2 of them result i oe
More informationDavid Vella, Skidmore College.
David Vella, Skidmore College dvella@skidmore.edu Geeratig Fuctios ad Expoetial Geeratig Fuctios Give a sequece {a } we ca associate to it two fuctios determied by power series: Its (ordiary) geeratig
More informationFrequency Response of FIR Filters
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we itroduce the idea of the frequecy respose of LTI systems, ad focus specifically o the frequecy respose of FIR filters.. Steady-state
More informationIf a subset E of R contains no open interval, is it of zero measure? For instance, is the set of irrationals in [0, 1] is of measure zero?
2 Lebesgue Measure I Chapter 1 we defied the cocept of a set of measure zero, ad we have observed that every coutable set is of measure zero. Here are some atural questios: If a subset E of R cotais a
More informationSigma notation. 2.1 Introduction
Sigma otatio. Itroductio We use sigma otatio to idicate the summatio process whe we have several (or ifiitely may) terms to add up. You may have see sigma otatio i earlier courses. It is used to idicate
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 information( 1) n (4x + 1) n. n=0
Problem 1 (10.6, #). Fid the radius of covergece for the series: ( 1) (4x + 1). For what values of x does the series coverge absolutely, ad for what values of x does the series coverge coditioally? Solutio.
More informationTHE ASYMPTOTIC COMPLEXITY OF MATRIX REDUCTION OVER FINITE FIELDS
THE ASYMPTOTIC COMPLEXITY OF MATRIX REDUCTION OVER FINITE FIELDS DEMETRES CHRISTOFIDES Abstract. Cosider a ivertible matrix over some field. The Gauss-Jorda elimiatio reduces this matrix to the idetity
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 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 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 informationEE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course
Sigal-EE Postal Correspodece Course 1 SAMPLE STUDY MATERIAL Electrical Egieerig EE / EEE Postal Correspodece Course GATE, IES & PSUs Sigal System Sigal-EE Postal Correspodece Course CONTENTS 1. SIGNAL
More informationMATH 10550, EXAM 3 SOLUTIONS
MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,
More information( ) ( ) ( ) notation: [ ]
Liear Algebra Vectors ad Matrices Fudametal Operatios with Vectors Vector: a directed lie segmets that has both magitude ad directio =,,,..., =,,,..., = where 1, 2,, are the otatio: [ ] 1 2 3 1 2 3 compoets
More informationSNAP Centre Workshop. Basic Algebraic Manipulation
SNAP Cetre Workshop Basic Algebraic Maipulatio 8 Simplifyig Algebraic Expressios Whe a expressio is writte i the most compact maer possible, it is cosidered to be simplified. Not Simplified: x(x + 4x)
More informationBernoulli numbers and the Euler-Maclaurin summation formula
Physics 6A Witer 006 Beroulli umbers ad the Euler-Maclauri summatio formula I this ote, I shall motivate the origi of the Euler-Maclauri summatio formula. I will also explai why the coefficiets o the right
More informationIt is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.
Taylor Polyomials ad Taylor Series It is ofte useful to approximate complicated fuctios usig simpler oes We cosider the task of approximatig a fuctio by a polyomial If f is at least -times differetiable
More informationOn the Inverse of a Certain Matrix Involving Binomial Coefficients
It. J. Cotemp. Math. Scieces, Vol. 3, 008, o. 3, 5-56 O the Iverse of a Certai Matrix Ivolvig Biomial Coefficiets Yoshiari Iaba Kitakuwada Seior High School Keihokushimoyuge, Ukyo-ku, Kyoto, 60-0534, Japa
More information