Stability of block LDL T factorization of a symmetric tridiagonal matrix
|
|
- Bennett Allison
- 5 years ago
- Views:
Transcription
1 Liner Algebr nd its Applictions 287 (1999) 181±189 Stbility of block LDL T fctoriztion of symmetric tridigonl mtrix Nichols J. Highm 1 Deprtment of Mthemtics, University of Mnchester, Mnchester M13 9PL, UK Received 30 September 1997; ccepted 20 April 1998 Submitted by D.D. Olesky Dedicted to Ludwig Elsner on the occsion of his 60th birthdy Abstrct For symmetric inde nite tridigonl mtrices, block LDL T fctoriztion without interchnges is shown to hve excellent numericl stbility when pivoting strtegy of Bunch is used to choose the dimension (1 or 2) of the pivots. Ó 1999 Elsevier Science Inc. All rights reserved. AMS clssi ction: 65F05; 65G05 Keywords: Tridigonl mtrix; Symmetric inde nite mtrix; Digonl pivoting method; LDL T fctoriztion; Growth fctor; Numericl stbility; Rounding error nlysis; LAPACK; LINPACK 1. Introduction Liner systems involving symmetric inde nite tridigonl mtrices rise in number of situtions. For exmple, Asen's method with prtil pivoting [1] produces fctoriztion PAP T ˆ LTL T of symmetric mtrix A, where P is permuttion mtrix, L is unit lower tringulr, nd T is tridigonl. To solve liner system Ax ˆ b using Asen's method it is necessry to solve system with coe cient mtrix T. A recent ppliction tht produces liner systems with symmetric tridigonl coe cient mtrices is Lnczos-bsed trust region method for unconstrined optimiztion of Gould et l. [8]. 1 E-mil: highm@m.mn.c.uk /99/$ ± see front mtter Ó 1999 Elsevier Science Inc. All rights reserved. PII: S ( 9 8 )
2 182 N.J. Highm / Liner Algebr nd its Applictions 287 (1999) 181±189 Symmetric tridigonl liner systems re most commonly solved by Gussin elimintion with prtil pivoting (GEPP) or by LDL T fctoriztion without pivoting. Neither method is completely stisfctory. GEPP destroys the symmetry, nd therefore cnnot be used to determine the inerti, while n LDL T fctoriztion yields the inerti of A directly from the digonl of D, but cn fil to exist nd its computtion cn be numericlly unstble when it does exist. A method tht promises to combine the bene ts of GEPP nd LDL T fctoriztion ws proposed by Bunch [3], but hs received little ttention in the literture. Bunch's ide is to compute block LDL T fctoriztion without interchnges, with prticulr strtegy for choosing the pivot size (1 or 2) t ech stge of the fctoriztion. Bunch's method requires less storge but slightly more computtion thn GEPP (see [3] for the detils). The purpose of this work is to exmine the numericl stbility of block LDL T fctoriztion with Bunch's pivoting strtegy. In Section 2 we de ne the pivoting strtegy nd explin how Bunch's derivtion of it yields bound of order 1 for the growth fctor. In Section 3 we show tht klk=kak cn be rbitrrily lrge nd explin why numericl stbility is therefore not consequence of error nlysis for generl block LU fctoriztion. We prove normwise bckwrd stbility of the method in Section 3, mking use of results of Highm [10] on the stbility of generl block LDL T fctoriztion. 2. Block LDL T fctoriztion nd the choice of pivot Consider the computtion of block LDL T fctoriztion without interchnges of symmetric tridigonl mtrix A 2 R nn. In the rst stge of the fctoriztion we choose n integer s ˆ 1 or 2 nd prtition A ˆ s s n s E C T : 2:1 n s C B If E is singulr for both choices of s then 11 ˆ 21 ˆ 0, but 21 ˆ 0 mens tht the rst row nd column is lredy in digonl form nd we cn skip to the next stge of the fctoriztion. Therefore, we cn ssume tht E is nonsingulr. Then we cn fctorize I s 0 E 0 Is E 1 C T A ˆ : 2:2 CE 1 I n s 0 B CE 1 C T 0 I n s This process cn be repeted recursively on the n s n s Schur complement S ˆ B CE 1 C T :
3 N.J. Highm / Liner Algebr nd its Applictions 287 (1999) 181± The result is fctoriztion A ˆ LDL T ; 2:3 where L is unit lower tringulr nd D is block digonl with ech digonl block hving dimension 1 or 2. While the fctoriztion lwys exists, whether it cn be computed in numericlly stble wy depends on the choice of pivots. Bunch's strtegy [3] for choosing the pivot size s t ech stge of the fctoriztion is fully de ned by describing the choice of the rst pivot. Algorithm 1 (Bunch's pivoting strtegy). This lgorithm determines the pivot size, s, for the rst stge of block LDL T fctoriztion pplied to symmetric tridigonl mtrix A 2 R nn. r :ˆ mxfj p ij j: i; j ˆ 1: ng (compute once, t the strt of the fctoriztion) :ˆ 5 1 =2 0:62 if rj 11 j P 2 21 s ˆ 1 else s ˆ 2 end Bunch excludes 11 from the mximiztion de ning r; we nd it more nturl to include it, becuse it increses the probbility tht 1 1 pivot will be chosen, while hving no e ect on the nlysis below. Bunch's choice of pivot cn be explined by considering element growth in the fctoriztion [3]. Since A is tridigonl, the mtrix C in (2.1) hs the form C ˆ s 1;s e 1 e T s, for unit vectors e 1 2 R n s, e s 2 R s. Hence the Schur complement S ˆ B 2 s 1;s et s E 1 e s e 1 e T 1 ; 2:4 which shows tht only the (1,1) element of the Schur complement di ers from the corresponding element of A. We now exmine the possible element growth in this position. Consider rst the cse s ˆ 1. We hve s 11 ˆ = 11: Hence, from the conditions in Algorithm 1, js 11 j 6 r r : The choice s ˆ 2 is mde when rj 11 j < 2 21 : 2:5
4 184 N.J. Highm / Liner Algebr nd its Applictions 287 (1999) 181±189 For s ˆ 2 we therefore hve det E ˆ j j j 22j=r < 0; 2:6 since < 1. Hence E is inde nite. We hve E 1 ˆ :7 det E nd so, from (2.4), s 11 ˆ =det E. Hence, using (2.5), 2 32 j 11j js 11 j 6 j 33 j r r2 1 r ˆ r 1 : We hve obtined bounds depending only on nd r for the size of the (1,1) element of the Schur complement. This element is not subsequently modi ed nd becomes digonl element of D. It follows tht growth in ny prticulr element tkes plce over single stge of the fctoriztion nd is not cumultive. The vlue of cn therefore be determined by equting the mximl element growth for n s ˆ 1 step with tht for n s ˆ 2 step. Hence we set r r ˆ r 1 ; p which is qudrtic in hving the positive root :ˆ 5 1 =2. With so chosen, the growth fctor q n for the fctoriztion stis es q n :ˆ mx i;j jd ij j mx i;j j ij j p 5 3 2:62: 3. Error nlysis Tht the growth fctor is nicely bounded does not, by itself, imply tht computtion of the block LDL T fctoriztion is numericlly stble process; see [10] for discussion in the cse of block LDL T fctoriztion of generl symmetric mtrices. From results on block LU fctoriztion [6], numericl stbility could be deduced if we could show tht klk=kak is suitbly bounded. We therefore exmine the size of the block CE 1 of L in (2.2). For s ˆ 1 we hve kce 1 k 1 ˆ j 21j j 11 j 6 r j 21 j ; nd the bound is shrp. It follows tht klk=kak cn be rbitrrily lrge. A prmetrized exmple is given by
5 N.J. Highm / Liner Algebr nd its Applictions 287 (1999) 181± A ˆ 1=2 ; 0 6 1; 1=2 2 for which the rst pivot is 1 1 nd 1 0 L ˆ 1=2 1 For s ˆ 2, ; D ˆ kce 1 k 1 ˆ j 32 jke T 2 E 1 k 1 6 j 32 j j 21j j 11 j j 32j 1 1 j 21 j ; r ; klk 1 =kak 1 1=2 =2: using (2.5) gin. This bound is shrp nd gin it esy to construct prmetrized exmple in which kce 1 k 1 =kak 1 cn be rbitrrily lrge. We conclude tht numericl stbility does not follow from results on generl block LU fctoriztion. Highm [10] proves the following generl result. We employ the usul model of oting point rithmetic fl x op y ˆ x op y 1 d ; jdj 6 u; op ˆ ; ; =; where u is the unit roundo. Absolute vlues of mtrices nd inequlities between mtrices re to be interpreted componentwise. Theorem 3.1. Let block LDL T fctoriztion with ny pivoting strtegy be pplied to symmetric mtrix A 2 R nn to yield the computed fctoriztion PAP T ^L ^D^L T, where P is permuttion mtrix nd D hs digonl blocks of dimension 1 or 2. Let ^x be the computed solution to Ax ˆ b obtined using the fctoriztion. Assume tht for ll liner systems Ey ˆ f involving 2 2 pivots E the computed solution ^x stis es E DE ^y ˆ f ; jdej 6 cu O u 2 jej; 3:1 where c is constnt. Then, P A DA 1 P T ˆ ^L ^D^L T ; A DA 2 ^x ˆ b; where jda i j 6 p n u jaj P T j^ljj ^Djj^L T jp O u 2 ; i ˆ 1: 2; 3:2 with p liner polynomil.
6 186 N.J. Highm / Liner Algebr nd its Applictions 287 (1999) 181±189 For our tridigonl A we cn set the polynomil p in Theorem 3.1 to be of zero degree, nd we hve P ˆ I. However, to verify tht the theorem is pplicble, we hve to check condition (3.1). It su ces to consider the rst stge of the fctoriztion. Suppose, rst, tht GEPP is used to solve Ey ˆ f. For 2 2 pivot E to be selected we must hve rj 11 j < j 21jr; which implies j 11 j < j 21 j < j 21 j: Hence GEPP interchnges rows 1 nd 2 of E nd fctorizes " PE ˆ ˆ 1 0 # ˆ LU: 21 From ([9], Theorem. 9.4), we hve the bckwrd error result PE DE ^y ˆ Pf ; jdej 6 6u O u 2 j^ljj ^Uj: Now, using (2.5), j 21 j j 22 j jljjuj 6 j 11 j j 21j j 22 j : 21 j 11 j 2 1 j 21 j p Hence jljjuj 6 p 5 PjEj. It follows tht (3.1) holds with c ˆ 6 5. Another wy to solve the liner systems Ey ˆ f is by the use of the explicit inverse, s is done in LINPACK [7] nd LAPACK [2] in their implementtions of block LDL T fctoriztion with the pivoting strtegyof Bunch nd Kufmn [4] for generl symmetric mtrices. The formul used in LINPACK nd LA- PACK is suitble here too 1 22 y ˆ f : 3:3 21 It is not hrd to show tht condition (3.1) holds when the formul (3.3) is used; the proof is very similr to tht in [10] for the pivoting strtegy of Bunch nd Kufmn. We hve now estblished tht Theorem 3.1 is pplicble. To deduce stbility of the fctoriztion we hve to show tht jljjdjjl T j is suitbly bounded in norm (we hve replced the computed L nd D by their exct counterprts, which ffects only the second order term of (3.2)). We write
7 N.J. Highm / Liner Algebr nd its Applictions 287 (1999) 181± I jej I jl T jljjdjjl T j ˆ 21 j jl 21 j jl S j jd S j jl T S " j jej jejjl T 21 ˆ j # jl 21 jjej jl 21 jjejjl T 21 j jl SjjD S jjl T S j ; 3:4 where L 21 nd E re from the rst stge of the fctoriztion. We rst bound F :ˆ jl 21 jjej jcjje 1 jjej: For s ˆ 1 we hve, trivilly, kf k 1 ˆ j 21 j 6 r. For s ˆ 2, using (2.5)±(2.7), je 1 1 j 22 j j 21 j j11 j j 21 j jjej j 21 j j 11 j j 21 j j 22 j 1 j 22 jj 11 j j 22 jj 21 j ˆ j 21 jj 11 j 2 21 j 11jj 22 j j 22 j 1 2 j 22j r j 21 j j 11j j 21 1 j 22j j r j 22j j 21 j 4 5: 3: j 11j j 21 j Hence kf k 1 6 j 32 jke 1 e T 2 je 1 jjejk 1 6 j 32j 1 2 j 11j j 21 j 1 6 r 1 2 j 21j r r 1 Now we bound G :ˆ jl 21 jjejjl T 21j. For s ˆ 1, < 8r: 3:6 kgk 1 ˆ 2 21 j 11 j 6 r < 2r: For s ˆ 2, G 6 jcjje 1 jjejje 1 jc T j ˆ 2 32 et 2 je 1 jjejje 1 j e 2 e 1 e T 1 : We bound the (2,2) element of je 1 jjejje 1 j strting with (3.5) nd nd tht 2 32 kgk j 11 j r 3 6 r < 16r: 3:7 2 1
8 188 N.J. Highm / Liner Algebr nd its Applictions 287 (1999) 181±189 We hve now bounded ll the terms in (3.4) except the term jl S jjd S jjl T S j. But L S nd D S re block LDL T fctors of the Schur complement of D in A, nd every Schur complement S stis es ksk M 6 q n kak M 6 2:62kAk M ; 3:8 where kak M ˆ mxj ij j: i;j From the bounds (3.6)±(3.8) nd the structure of the (2,2) block in (3.4) we deduce tht kjljjdjjl T jk M :62kAk M < 42kAk M : The following result summrizes the stbility of block LDL T fctoriztion with Bunch's pivoting strtegy. Theorem 3.2. Let block LDL T fctoriztion with the pivoting strtegy of Algorithm 1 be pplied to symmetric tridigonl mtrix A 2 R nn to yield the computed fctoriztion A ^L ^D^L T, nd let ^x be the computed solution to Ax ˆ b obtined using the fctoriztion. Assume tht ll liner systems Ey ˆ f involving 2 2 pivots E re solved by GEPP or by using the explicit inverse formul (3.3). Then A DA 1 ˆ ^L ^D^L T ; A DA 2 ^x ˆ b; where kda i k M 6 cukak M O u 2 ; i ˆ 1: 2; 3:9 with c constnt. 4. Conclusions Theorem 3.2 shows tht block LDL T fctoriztion with the pivoting strtegy of Algorithm 1 is normwise bckwrd stble wy to fctorize symmetric tridigonl mtrix A nd to solve liner system Ax ˆ b. Block LDL T fctoriztion therefore provides n ttrctive lterntive to GEPP for solving such liner systems. Since the inerti of A is the sme s tht of the block digonl fctor D, the fctoriztion lso provides normwise bckwrd stble wy to compute the inerti. However, for computing inertis of symmetric tridigonl mtrices stndrd LDL T fctoriztion without pivoting hs the stronger componentwise reltive form of bckwrd stbility ([5], Lemm 5.3), nd so is preferble in the bisection method for computing eigenvlues, for exmple.
9 Acknowledgements N.J. Highm / Liner Algebr nd its Applictions 287 (1999) 181± I thnk Nick Gould for pointing out the open question of the numericl stbility of block LDL T fctoriztion with Bunch's pivoting strtegy. References [1] J.O. Asen, On the reduction of symmetric mtrix to tridigonl form, BIT 11 (1971) 233± 242. [2] E. Anderson, Z. Bi, C.H. Bischof, J.W. Demmel, J.J. Dongrr, J.J. Du Croz, A. Greenbum, S.J. Hmmrling, A. McKenney, S. Ostrouchov, D.C. Sorensen, LAPACK Users' Guide, Relese nd ed., Society for Industril nd Applied Mthemtics, Phildelphi, PA, USA, 1995, xix pp. ISBN [3] J.R. Bunch, Prtil pivoting strtegies for symmetric mtrices, SIAM J. Numer. Anl. 11 (3) (1974) 521±528. [4] J.R. Bunch, L. Kufmn, Some stble methods for clculting inerti nd solving symmetric liner systems, Mth. Comp. 31 (137) (1977) 163±179. [5] J.W. Demmel, Applied numericl liner lgebr, Society for Industril nd Applied Mthemtics, Phildelphi, PA, USA, xi+419 pp. ISBN [6] J.W. Demmel, N.J. Highm, R.S. Schreiber, Stbility of block LU fctoriztion, Numericl Liner Algebr with Applictions 2 (2) (1995) 173±190. [7] J.J. Dongrr, J.R. Bunch, C.B. Moler, G.W. Stewrt, LINPACK users' guide, Society for Industril nd Applied Mthemtics, Phildelphi, PA, USA, ISBN X. [8] N.I.M. Gould, S. Lucidi, M. Rom, P.L. Toint, Solving the trust-region subproblem using the Lnczos method, Report RAL , Atls Centre, Rutherford Appleton Lbortory, Didcot, Oxon, UK, 1997, p. 28. [9] N.J. Highm, Accurcy nd stbility of numericl lgorithms, Society for Industril nd Applied Mthemtics, Phildelphi, PA, USA, 1996, xxviii+688 pp. ISBN [10] N.J. Highm, Stbility of the digonl pivoting method with prtil pivoting, SIAM J. Mtrix Anl. Appl. 18 (1) (1997) 52±65.
Contents. Outline. Structured Rank Matrices Lecture 2: The theorem Proofs Examples related to structured ranks References. Structure Transport
Contents Structured Rnk Mtrices Lecture 2: Mrc Vn Brel nd Rf Vndebril Dept. of Computer Science, K.U.Leuven, Belgium Chemnitz, Germny, 26-30 September 2011 1 Exmples relted to structured rnks 2 2 / 26
More informationMatrices, Moments and Quadrature, cont d
Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk bout solving systems of liner equtions. These re problems tht give couple of equtions with couple of unknowns, like: 6 2 3 7 4
More informationHW3, Math 307. CSUF. Spring 2007.
HW, Mth 7. CSUF. Spring 7. Nsser M. Abbsi Spring 7 Compiled on November 5, 8 t 8:8m public Contents Section.6, problem Section.6, problem Section.6, problem 5 Section.6, problem 7 6 5 Section.6, problem
More informationNumerical Linear Algebra Assignment 008
Numericl Liner Algebr Assignment 008 Nguyen Qun B Hong Students t Fculty of Mth nd Computer Science, Ho Chi Minh University of Science, Vietnm emil. nguyenqunbhong@gmil.com blog. http://hongnguyenqunb.wordpress.com
More informationHere we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system.
Section 24 Nonsingulr Liner Systems Here we study squre liner systems nd properties of their coefficient mtrices s they relte to the solution set of the liner system Let A be n n Then we know from previous
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk out solving systems of liner equtions. These re prolems tht give couple of equtions with couple of unknowns, like: 6= x + x 7=
More informationElements of Matrix Algebra
Elements of Mtrix Algebr Klus Neusser Kurt Schmidheiny September 30, 2015 Contents 1 Definitions 2 2 Mtrix opertions 3 3 Rnk of Mtrix 5 4 Specil Functions of Qudrtic Mtrices 6 4.1 Trce of Mtrix.........................
More informationMath 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of
More informationN 0 completions on partial matrices
N 0 completions on prtil mtrices C. Jordán C. Mendes Arújo Jun R. Torregros Instituto de Mtemátic Multidisciplinr / Centro de Mtemátic Universidd Politécnic de Vlenci / Universidde do Minho Cmino de Ver
More informationEngineering Analysis ENG 3420 Fall Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00
Engineering Anlysis ENG 3420 Fll 2009 Dn C. Mrinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00 Lecture 13 Lst time: Problem solving in preprtion for the quiz Liner Algebr Concepts Vector Spces,
More information1 Linear Least Squares
Lest Squres Pge 1 1 Liner Lest Squres I will try to be consistent in nottion, with n being the number of dt points, nd m < n being the number of prmeters in model function. We re interested in solving
More informationIntroduction to Determinants. Remarks. Remarks. The determinant applies in the case of square matrices
Introduction to Determinnts Remrks The determinnt pplies in the cse of squre mtrices squre mtrix is nonsingulr if nd only if its determinnt not zero, hence the term determinnt Nonsingulr mtrices re sometimes
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationThe Algebra (al-jabr) of Matrices
Section : Mtri lgebr nd Clculus Wshkewicz College of Engineering he lgebr (l-jbr) of Mtrices lgebr s brnch of mthemtics is much broder thn elementry lgebr ll of us studied in our high school dys. In sense
More informationPart IB Numerical Analysis
Prt IB Numericl Anlysis Theorems with proof Bsed on lectures by G. Moore Notes tken by Dexter Chu Lent 2016 These notes re not endorsed by the lecturers, nd I hve modified them (often significntly) fter
More informationChapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY
Chpter 3 MTRIX In this chpter: Definition nd terms Specil Mtrices Mtrix Opertion: Trnspose, Equlity, Sum, Difference, Sclr Multipliction, Mtrix Multipliction, Determinnt, Inverse ppliction of Mtrix in
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More informationMath Lecture 23
Mth 8 - Lecture 3 Dyln Zwick Fll 3 In our lst lecture we delt with solutions to the system: x = Ax where A is n n n mtrix with n distinct eigenvlues. As promised, tody we will del with the question of
More informationS. S. Dragomir. 2, we have the inequality. b a
Bull Koren Mth Soc 005 No pp 3 30 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Abstrct Compnions of Ostrowski s integrl ineulity for bsolutely
More informationLecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.
Lecture 3 3 Solving liner equtions In this lecture we will discuss lgorithms for solving systems of liner equtions Multiplictive identity Let us restrict ourselves to considering squre mtrices since one
More informationLecture 4: Piecewise Cubic Interpolation
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 4: Piecewise Cubic Interpoltion Compiled 5 September In this lecture we consider piecewise cubic interpoltion
More informationLecture Note 9: Orthogonal Reduction
MATH : Computtionl Methods of Liner Algebr 1 The Row Echelon Form Lecture Note 9: Orthogonl Reduction Our trget is to solve the norml eution: Xinyi Zeng Deprtment of Mthemticl Sciences, UTEP A t Ax = A
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationQuantum Physics II (8.05) Fall 2013 Assignment 2
Quntum Physics II (8.05) Fll 2013 Assignment 2 Msschusetts Institute of Technology Physics Deprtment Due Fridy September 20, 2013 September 13, 2013 3:00 pm Suggested Reding Continued from lst week: 1.
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationJim Lambers MAT 169 Fall Semester Lecture 4 Notes
Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of
More informationThe Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 11
The Islmic University of Gz Fculty of Engineering Civil Engineering Deprtment Numericl Anlysis ECIV 6 Chpter Specil Mtrices nd Guss-Siedel Associte Prof Mzen Abultyef Civil Engineering Deprtment, The Islmic
More informationS. S. Dragomir. 1. Introduction. In [1], Guessab and Schmeisser have proved among others, the following companion of Ostrowski s inequality:
FACTA UNIVERSITATIS NIŠ) Ser Mth Inform 9 00) 6 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Dedicted to Prof G Mstroinni for his 65th birthdy
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationTANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION
Indin Journl of Mthemtics nd Mthemticl Sciences Vol. 7, No., (June ) : 9-38 TANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More informationMAC-solutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
More informationENGI 3424 Engineering Mathematics Five Tutorial Examples of Partial Fractions
ENGI 44 Engineering Mthemtics Five Tutoril Exmples o Prtil Frctions 1. Express x in prtil rctions: x 4 x 4 x 4 b x x x x Both denomintors re liner non-repeted ctors. The cover-up rule my be used: 4 4 4
More informationUNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3
UNIFORM CONVERGENCE Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 Suppose f n : Ω R or f n : Ω C is sequence of rel or complex functions, nd f n f s n in some sense. Furthermore,
More informationModule 6: LINEAR TRANSFORMATIONS
Module 6: LINEAR TRANSFORMATIONS. Trnsformtions nd mtrices Trnsformtions re generliztions of functions. A vector x in some set S n is mpped into m nother vector y T( x). A trnsformtion is liner if, for
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationHomework 5 solutions
Section.: E ; AP, Section.: E,,; AP,,, Section. Homework solutions. Consider the two upper-tringulr mtrices: b b b A, B b b. b () Show tht their product C = AB is lso upper-tringulr. The product is b b
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More informationCHAPTER 4a. ROOTS OF EQUATIONS
CHAPTER 4. ROOTS OF EQUATIONS A. J. Clrk School o Engineering Deprtment o Civil nd Environmentl Engineering by Dr. Ibrhim A. Asskk Spring 00 ENCE 03 - Computtion Methods in Civil Engineering II Deprtment
More informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry dierentil eqution (ODE) du f(t) dt with initil condition u() : Just
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationNumerical Integration
Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl
More informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
More informationAlgebra Of Matrices & Determinants
lgebr Of Mtrices & Determinnts Importnt erms Definitions & Formule 0 Mtrix - bsic introduction: mtrix hving m rows nd n columns is clled mtrix of order m n (red s m b n mtrix) nd mtrix of order lso in
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More informationSemigroup of generalized inverses of matrices
Semigroup of generlized inverses of mtrices Hnif Zekroui nd Sid Guedjib Abstrct. The pper is divided into two principl prts. In the first one, we give the set of generlized inverses of mtrix A structure
More informationAQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system
Complex Numbers Section 1: Introduction to Complex Numbers Notes nd Exmples These notes contin subsections on The number system Adding nd subtrcting complex numbers Multiplying complex numbers Complex
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationChapter 3 Solving Nonlinear Equations
Chpter 3 Solving Nonliner Equtions 3.1 Introduction The nonliner function of unknown vrible x is in the form of where n could be non-integer. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,
More informationState space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies
Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response
More informationThe use of a so called graphing calculator or programmable calculator is not permitted. Simple scientific calculators are allowed.
ERASMUS UNIVERSITY ROTTERDAM Informtion concerning the Entrnce exmintion Mthemtics level 1 for Interntionl Bchelor in Communiction nd Medi Generl informtion Avilble time: 2 hours 30 minutes. The exmintion
More informationLecture 3. Limits of Functions and Continuity
Lecture 3 Limits of Functions nd Continuity Audrey Terrs April 26, 21 1 Limits of Functions Notes I m skipping the lst section of Chpter 6 of Lng; the section bout open nd closed sets We cn probbly live
More informationNumerical quadrature based on interpolating functions: A MATLAB implementation
SEMINAR REPORT Numericl qudrture bsed on interpolting functions: A MATLAB implementtion by Venkt Ayylsomyjul A seminr report submitted in prtil fulfillment for the degree of Mster of Science (M.Sc) in
More informationIntegral equations, eigenvalue, function interpolation
Integrl equtions, eigenvlue, function interpoltion Mrcin Chrząszcz mchrzsz@cernch Monte Crlo methods, 26 My, 2016 1 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion
More informationCBE 291b - Computation And Optimization For Engineers
The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn
More informationA recursive construction of efficiently decodable list-disjunct matrices
CSE 709: Compressed Sensing nd Group Testing. Prt I Lecturers: Hung Q. Ngo nd Atri Rudr SUNY t Bufflo, Fll 2011 Lst updte: October 13, 2011 A recursive construction of efficiently decodble list-disjunct
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More informationChapter 3 Polynomials
Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling
More information4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for.
4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX Some reliminries: Let A be rel symmetric mtrix. Let Cos θ ; (where we choose θ π for Cos θ 4 purposes of convergence of the scheme)
More informationMath 554 Integration
Mth 554 Integrtion Hndout #9 4/12/96 Defn. A collection of n + 1 distinct points of the intervl [, b] P := {x 0 = < x 1 < < x i 1 < x i < < b =: x n } is clled prtition of the intervl. In this cse, we
More informationSOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 2014
SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 014 Mrk Scheme: Ech prt of Question 1 is worth four mrks which re wrded solely for the correct nswer.
More informationMath 360: A primitive integral and elementary functions
Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:
More informationMath 4310 Solutions to homework 1 Due 9/1/16
Mth 4310 Solutions to homework 1 Due 9/1/16 1. Use the Eucliden lgorithm to find the following gretest common divisors. () gcd(252, 180) = 36 (b) gcd(513, 187) = 1 (c) gcd(7684, 4148) = 68 252 = 180 1
More informationMath 61CM - Solutions to homework 9
Mth 61CM - Solutions to homework 9 Cédric De Groote November 30 th, 2018 Problem 1: Recll tht the left limit of function f t point c is defined s follows: lim f(x) = l x c if for ny > 0 there exists δ
More informationIs there an easy way to find examples of such triples? Why yes! Just look at an ordinary multiplication table to find them!
PUSHING PYTHAGORAS 009 Jmes Tnton A triple of integers ( bc,, ) is clled Pythgoren triple if exmple, some clssic triples re ( 3,4,5 ), ( 5,1,13 ), ( ) fond of ( 0,1,9 ) nd ( 119,10,169 ). + b = c. For
More informationDeterminants Chapter 3
Determinnts hpter Specil se : x Mtrix Definition : the determinnt is sclr quntity defined for ny squre n x n mtrix nd denoted y or det(). x se ecll : this expression ppers in the formul for x mtrix inverse!
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
More informationChapter 2. Determinants
Chpter Determinnts The Determinnt Function Recll tht the X mtrix A c b d is invertible if d-bc0. The expression d-bc occurs so frequently tht it hs nme; it is clled the determinnt of the mtrix A nd is
More information1.9 C 2 inner variations
46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationMatrix Algebra. Matrix Addition, Scalar Multiplication and Transposition. Linear Algebra I 24
Mtrix lger Mtrix ddition, Sclr Multipliction nd rnsposition Mtrix lger Section.. Mtrix ddition, Sclr Multipliction nd rnsposition rectngulr rry of numers is clled mtrix ( the plurl is mtrices ) nd the
More informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationDiscrete Least-squares Approximations
Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve
More information1 Online Learning and Regret Minimization
2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in
More informationQuadratic Forms. Quadratic Forms
Qudrtic Forms Recll the Simon & Blume excerpt from n erlier lecture which sid tht the min tsk of clculus is to pproximte nonliner functions with liner functions. It s ctully more ccurte to sy tht we pproximte
More informationJournal of Inequalities in Pure and Applied Mathematics
Journl o Inequlities in Pure nd Applied Mthemtics http://jipm.vu.edu.u/ Volume 6, Issue 4, Article 6, 2005 MROMORPHIC UNCTION THAT SHARS ON SMALL UNCTION WITH ITS DRIVATIV QINCAI ZHAN SCHOOL O INORMATION
More informationTHE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS. Circa 1870, G. Zolotarev observed that the Legendre symbol ( a p
THE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS PETE L CLARK Circ 1870, Zolotrev observed tht the Legendre symbol ( p ) cn be interpreted s the sign of multipliction by viewed s permuttion of the set Z/pZ
More informationAn approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95
An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech
More informationarxiv: v2 [math.nt] 2 Feb 2015
rxiv:407666v [mthnt] Fe 05 Integer Powers of Complex Tridigonl Anti-Tridigonl Mtrices Htice Kür Duru &Durmuş Bozkurt Deprtment of Mthemtics, Science Fculty of Selçuk University Jnury, 08 Astrct In this
More informationMonte Carlo method in solving numerical integration and differential equation
Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationQuadrature Rules for Evaluation of Hyper Singular Integrals
Applied Mthemticl Sciences, Vol., 01, no. 117, 539-55 HIKARI Ltd, www.m-hikri.com http://d.doi.org/10.19/ms.01.75 Qudrture Rules or Evlution o Hyper Singulr Integrls Prsnt Kumr Mohnty Deprtment o Mthemtics
More informationThe Riemann-Lebesgue Lemma
Physics 215 Winter 218 The Riemnn-Lebesgue Lemm The Riemnn Lebesgue Lemm is one of the most importnt results of Fourier nlysis nd symptotic nlysis. It hs mny physics pplictions, especilly in studies of
More informationMatching with Multiple Applications: The Limiting Case
Mtching with Multiple Applictions: The Limiting Cse Jmes W. Albrecht y Georgetown Uniersity Susn B. Vromn Georgetown Uniersity August 003 Pieter A. Gutier Ersmus Uniersity Tinbergen Institute Abstrct We
More information1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.
1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt
More informationMatrix Solution to Linear Equations and Markov Chains
Trding Systems nd Methods, Fifth Edition By Perry J. Kufmn Copyright 2005, 2013 by Perry J. Kufmn APPENDIX 2 Mtrix Solution to Liner Equtions nd Mrkov Chins DIRECT SOLUTION AND CONVERGENCE METHOD Before
More informationAN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
Applied Mthemtics E-Notes, 5(005), 53-60 c ISSN 1607-510 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ AN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties
More informationPARTIAL FRACTION DECOMPOSITION
PARTIAL FRACTION DECOMPOSITION LARRY SUSANKA 1. Fcts bout Polynomils nd Nottion We must ssemble some tools nd nottion to prove the existence of the stndrd prtil frction decomposition, used s n integrtion
More informationA Bernstein polynomial approach for solution of nonlinear integral equations
Avilble online t wwwisr-publictionscom/jns J Nonliner Sci Appl, 10 (2017), 4638 4647 Reserch Article Journl Homepge: wwwtjnscom - wwwisr-publictionscom/jns A Bernstein polynomil pproch for solution of
More informationMTH 5102 Linear Algebra Practice Exam 1 - Solutions Feb. 9, 2016
Nme (Lst nme, First nme): MTH 502 Liner Algebr Prctice Exm - Solutions Feb 9, 206 Exm Instructions: You hve hour & 0 minutes to complete the exm There re totl of 6 problems You must show your work Prtil
More information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
More informationCzechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction
Czechoslovk Mthemticl Journl, 55 (130) (2005), 933 940 ESTIMATES OF THE REMAINDER IN TAYLOR S THEOREM USING THE HENSTOCK-KURZWEIL INTEGRAL, Abbotsford (Received Jnury 22, 2003) Abstrct. When rel-vlued
More informationarxiv:math/ v2 [math.ho] 16 Dec 2003
rxiv:mth/0312293v2 [mth.ho] 16 Dec 2003 Clssicl Lebesgue Integrtion Theorems for the Riemnn Integrl Josh Isrlowitz 244 Ridge Rd. Rutherford, NJ 07070 jbi2@njit.edu Februry 1, 2008 Abstrct In this pper,
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationThe Velocity Factor of an Insulated Two-Wire Transmission Line
The Velocity Fctor of n Insulted Two-Wire Trnsmission Line Problem Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 Mrch 7, 008 Estimte the velocity fctor F = v/c nd the
More information