Some applications of the Hermite matrix polynomials series expansions 1
|
|
- Clara Martha Miles
- 5 years ago
- Views:
Transcription
1 Journal of Computational and Applied Mathematics 99 (1998) Some applications of the Hermite matrix polynomials series expansions 1 Emilio Defez, Lucas Jodar Departamento de Matematica Aplicada, Universidad Politecnica de Valencia, Spain Received 0 October 1997; received in revised form 5 May 1998 Abstract This paper deals with Hermite matrix polynomials expansions of some relevant matrix functions appearing in the solution of dierential systems. Properties of Hermite matrix polynomials such as the three terms recurrence formula permit an ecient computation of matrix functions avoiding important computational drawbacks of other well-known methods. Results are applied to compute accurate approximations of certain dierential systems in terms of Hermite matrix polynomials. c 1998 Elsevier Science B.V. All rights reserved. AMS classication: 15A60; 33C5; 34A50; 41A10 Keywords: Hermite matrix polynomials; Matrix functions; Dierential systems; Approximate solutions; Error bounds; Sylvester equation 1. Introduction and preliminaries The evaluation of matrix functions is frequent in the solution of dierential systems. So, the system Y = AY; Y (0) = y 0 ; (1) where A is matrix and y 0 is a vector, arises in the semidiscretization of the heat equation [17]. The matrix dierential problem Y + A Y =0; Y(0) = P; Y (0) = Q; () Corresponding author. edefez@mat.upv.es. 1 This work has been partially supported by the Generalitat Valenciana grant GV-C-CN and the Spanish D.G.I.C.Y.T grant PB CO /98/$ see front matter c 1998 Elsevier Science B.V. All rights reserved. PII: S (98)
2 106 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) where A is a matrix and P and Q are vectors, arises from semidiscretization of the wave equation [15]. The Sylvester matrix dierential equation X = AX + XB; X(0) = C; (3) where A, B and C are matrices, appears in systems stability and control [4, p. 6]. The solutions of problems (1) (3) can be expressed in terms of exp(at); cos(at); sin(at) and exp(bt) and the computation of these matrix functions has motivated many and varied approaches. An excelent survey about the matrix exponential is [11] and the study of cos(at) is treated in [1, 14, 15]. Some of the drawbacks of the existing methods are (i) The computation of eigenvalues or eigenvectors [11, 16]; (ii) The computation of the inverses of matrices (see Pade methods, [5, p. 557]) [11, 14]; (iii) Storage problems and expensive computational time [5, Chap. 1]; (iv) Round-o accumulation errors [5, p. 551]; (v) Diculties for computing approximations of matrix functions with a prexed accuracy [15, 16]. In this paper we propose a new method for computing the above matrix functions using Hermite matrix polynomials which avoids the quoted computational diculties. Results are applied to construct approximations of problems (1) (3) with a prexed accuracy in a bounded domain. If D 0 is the complex plane cut along the negative real axis and log(z) denotes the principal logarithm of z, [13, p. 7], then z 1= represents exp( 1log(z)). If A is a matrix in Cr r ; its -norm denoted A is dened by A = Ax = x, where for a vector y in C r, y denotes the usual euclidean norm of y, y =(y T y) 1=. The set of all the eigenvalues of A is denoted by (A). If f(z) and g(z) are holomorphic functions of the complex variable z, which are dened in an open set of the complex plane, and A is a matrix in C r r such that (A), then from the properties of the matrix functional calculus, [3, p. 558], it follows that f(a)g(a)=g(a)f(a). If A is a matrix with (A) D 0, then A 1= = A denotes the image by z 1= of the matrix functional calculus acting on the matrix A. We say that A is a positive stable matrix if Re(z) 0 for all z (A). For the sake of clarity in the presentation we recall some properties of the Hermite matrix polynomials which will be used below and that have been established in [7], see also [8]. If A is a positive stable matrix in C r r, the nth Hermite matrix polynomial is dened by [n=] H n (x; A)=n! ( 1) k ( A) n k x n k ; (4) k!(n k)! and satisfy the three terms recurrence relationship H n (x; A)=Ix (A)H n 1 (x; A) (n 1)H n (x; A); n 1 H 1 (x; A)=0; H 0 (x; A)=I; (5) where I is the identity matrix in C r r. By [7] we also have e xt A t I = 1 n! H n(x; A)t n ; t + ; (6)
3 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) and if H n (x) denotes the classical nth scalar Hermite polynomial, then one gets H n (x; A) 6 H n( ia A ) ; x a; n N; (7) i n H n ( A ia ) n!i n t n = exp( t a A + t ); t : (8) If A(k; n) and B(k; n) are matrices on C r r Lemma 11 of [1, p. 57] it follows that for ; k 0, in an analogous way to the proof of [n=] A(k; n k); k 0 A(k; n)= B(k; n)= k 0 B(k; n k): (9) If B is a matrix in C r r and n 0 is a positive integer we denote by A(B; n 0 ), D(B; n 0 ) and E(B; n 0 ) the real numbers: n 0 [n=] A(B; n 0 )= E(B; n 0 )= n=0 n 0 n=0 ( B ) n k k!(n k)! ; D(B; n 0)= n 0 n=0 ( B ) (n k) k!((n k))! ; (10) ( B ) (n k)+1 k!((n k) + 1)! : (11) This paper is organized as follows. In Section some new properties of Hermite matrix polynomials are established. Section 3 deals with the Hermite matrix polynomial series expansions of e At, sin(at) and cos(at) of an arbitrary matrix as well as with theirs nite series truncation with a prexed accuracy in a bounded domain. Finally, in Section 4 analytic-numerical approximations of problems (1) (3) are contructed in terms of Hermite matrix polynomial series expansions. Given an admissible error 0 and a bounded domain D, an approximation in terms of Hermite matrix polynomials is contructed so that the error with respect to the exact solution is uniformly upper bounded by in D.. On Hermite matrix polynomials Let A be a positive stable matrix. By (6) it follows that e xt A+t I = H n (x; A) t n : n! Hence ( ( Ax) n t n = n! )( t n n! ) H n (x; A) t n n! (1)
4 108 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) and by (9) one gets ( )( ) t n H n (x; A) t n = [n=] H n k (x; A) n! n! k!(n k)! t n : (13) By identication of the coecient of t n in (1) and (13) one gets x n I =( [n=] A) n n! k!(n k)! H n k(x; A); x + : (14) Lemma.1. Let be A a positive stable matrix; K and integer. Then H n (x; A) 6 n! n K n e x ; x K A : (15) Proof. It is clear that for n = 0 the inequality (15) holds true. Suppose that (15) is true for k =0; 1;:::;n. Taking norms in the three terms formula (5) and using the induction hypothesis one gets H n+1 (x; A) 6 x (A) H n (x; A) +n H n 1 (x; A) 6 K n! n K n e x +n (n 1)! n 1 K n 1 e x = n! n K n 1 e x (K + n ) = (n + 1)! n K n 1 e x K n + : n +1 n +1 (16) Since K ( n + K 6 + ) = K + 6 K ; n +1 n +1 by (16) one gets (15) for n +1. The following result that may be regarded as a matrix version of the algorithm given in [10], provides an ecient procedure for computing nite matrix polynomial series expansions in terms of matrix polynomials satisfying a three terms formula. In particular, it is applicable to the Hermite matrix polynomials. Note the remarkable fact from a computational point of view that the evaluation of a matrix polynomial at x is only expressed in terms of P 0 (x). Theorem.1. Let {P n (x)} be a sequence of matrix polynomials such that A n P n (x)=(xi B n )P n 1 (x) C n P n (x); n 1;
5 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) where A n is an invertible matrix in C r r and the degree of P n (x) is n. Let Q(x) be a matrix polynomial dened by Q(x)= E j P j (x); j=0 where E j is a matrix in C r r. Let x be any real number and consider the sequence of matrices dened by D n = E n ; D n 1 = E n 1 + D n A 1 n (xi B n ); and for j = n ;:::;0; D j = E j + D j+1 A 1 j+1(xi B j+1 ) D j+ A 1 j+c j+ : Then Q(x)=D 0 P 0 (x). Proof. See [10]. 3. Hermite matrix polynomials series expansions We begin this section with Hermite matrix polynomial series expansion of exp(bt); sin(bt) and cos(bt) for matrices satisfying the spectral property Re(z) Im(z) for all z (B): (17) Theorem 3.1. Let B be a matrix in C r r satisfying (17). Then e Bx = e cos(bx)= 1 e sin(bx)= 1 e H n (x; 1 B ) ; x + ; (18) n! ( 1) n (n)! H n ( 1) n (n + 1)! H n+1 (x; 1 B ) ; x + ; (19) (x; 1 B ) ; x + : (0) Let E(B; x; n); C(B; x; n) and S(B; x; n) be respectively the nth partial sum of series (18) (0) respectively. Given 0 and c 0; let a = c + and let n 0 ;n 1 and n be the rst positive integers satisfying A(Ba; n 0 ) exp(a B +1) e ; D(Ba; n 1 ) e(cosh(a B ) ); E(Ba; n ) e(sinh(a B ) ); (1) () (3)
6 110 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) where A; D and E are dened by (10) and (11). Then for x c and n max{n 0 ;n 1 ;n } one gets e Bx E(B; x; n) cos(bx) C(B; x; n) ; sin(bx) S(B; x; n) : (4) (5) (6) Proof. Let A = 1 B. By the spectral mapping theorem [3, p. 569] and (17) it follows that { } ( ) 1 1 (A)= b ; b (B) ; Re b = 1 { (Re(b)) (Im(b)) } 0; b (B): Thus A is a positive stable matrix and taking t = 1 in (6), B = A, one gets e Bx = e n! H n B ; x + : If n 0 1 one gets n 0 ( e ebx k H k! x; 1 ) B 6 ( e k n 0 k! H k x; 1 ) B = e H k (x; 1 B ) : (7) k n 0 k! By (7) for 0, c 0, a = c + and x a, it follows that e H k (x; 1 B ) 6e H k ( B ai ) n 0 H k ( B ai ) e : (8) k n 0 k! k!i k k!i k k 0 By (8), (10), (7) and (8) one gets n 0 ( e ebx k! H k x; 1 ) n 0 H k ( B ai ) B 6 exp(a B +1) e k!i k By [1, p. 57] one gets lim n [k=] l=0 ( B a) k l l!(k l)! = exp(a B +1): = exp(a B +1) A(Ba; n 0 ): (9) Hence, taking the rst positive integer n 0 satisfying (1), then by (9) one gets n 0 ( e ebx k! H k x; 1 ) B 6; x c; n n 0 : Considering (14) for the positive stable matrix A = 1 B, it follows that x n I = B n (n)! k!((n k))! H (n k) B :
7 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) Taking into account the series expansion of cos(bx) and (9) we can write cos(bx)= ( 1) n B n x n = ( 1) n (n)! k! ((n k))! H (n k) B = ( 1) n+k k! (n)! H n B = ( ) ( 1) k ( 1) n k! (n)! H n B : k 0 Hence cos(bx)= 1 e By [1, p. 57] one gets lim n k l=0 k 0 ( 1) n (n)! H n B ; x + : (30) ( B a) (k l) l!((k l))! = e cosh(a B ): Let us consider the n 1 th partial sum of the series (30). Then in an analogous way to the previous computations, for x c one gets cos(bx) C(B; x; n 1 ) 6 cosh(a B ) 1 e D(aB; n 1); a= c + ; where D(aB; n 1 ) is given by (10). Taking the rst integer n 1 1 satisfying () one gets (5). By (14) we also have x n+1 I = B (n+1) (n + 1)! k!((n k) + 1)! H (n k)+1 B : and in an analogous way to the series expansion of cos(bx), one gets (0). By [1, p. 57] one gets lim n k l=0 ( B a) (k l)+1 l!((k l) + 1)! = e sinh(a B ): Then in an analogous way to the previous computations, one gets (6). Hence the result is established. Now we show that condition (17) can be removed in the construction of accurate Hermite matrix polynomial expansions of exp(at); cos(at) and sin(at) of an arbitrary matrix A in C r r. Lemma 3.1. Let A be a matrix in C r r and let any positive number such that max{ Re(z) + Im(z) ; z (A)}: Then the matrix B = A + I satises (17). Proof. By the spectral mapping theorem (A + I)= {z + ; z (A)} and for z (A) one gets Re(z + ) = Re(z) + Re(z) Im(z) = Im(z + ). Hence the result is established.
8 11 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) Corollary 3.1. Let A be a matrix in C r r and let be a number satisfying the condition of Lemma 3.1. Let 0; c 0; a= c +. With the notacion of Theorem 3.1 it follows that (i) If n 0 satises (1) for B = A + I; then e Ax e x E(A + I;x;n) ; n n 0 ; x c: (31) (ii) Let n 1 and n be positive integers satisfying () and (3) respectively for == ; and B = A + I; then for x c and n max{n 1 ;n }; sin(ax) {cos(x)s(a + I;x;n) sin(x)c(a + I;x;n)} ; (3) cos(ax) {cos(x)c(a + I;x;n)+ sin(x)s(a + I;x;n)} : (33) Proof. (i) The result is a consequence of Theorem 3.1 and the formula e Ax =e x e (A+I)x. The proof of part (ii) follows from the expressions sin(ax) = sin[(a + I)x] cos(x) cos[(a + I)x] sin(x); cos(ax) = cos[(a + I)x] cos(x)+ sin[(a + I)x] sin(x); and Theorem 3.1. In the following example we illustrate the use of Hermite matrix polynomials for computing e A, sin A and cos A for a matrix A satisfying (17). As it has been proved in Corollary 3.1, condition (17) is not necessary taking an appropiate value of in accordance with Lemma 3.1. It is interesting to point out that algorithm of Theorem.1 has been adapted for the Hermite matrix polynomials using (5). Computations have been perfomed using Mathematica version ::1. Example 3.1. Consider the matrix A = 0 1 ; 1 1 where (A)={1; ; }. Using the minimal theorem [3, p. 571] one gets the exact value of e A : e e e e A = e +e e e e e + e e e e 14: : : = 1: : : : 4: : :
9 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) Given an admissible error =10 5, from (1), taking n 0 = 30 the approximation E(A; 1; 30) provides the required accuracy 30 1 E (A; 1; 30) = e n! H n (1; 1 ) A n=0 14: : : = 1: : : ; 4: : : e A E (A; 1; 30) =3: : Of course, in practice the number of terms required to obtain a prexed accuracy uses to be smaller than the one provided by (1), because the error bounds given by Theorem 3.1 is valid for any matrix satisfying (17). So for instance taking n 0 = 19 one gets 14: : : E (A; 1; 19) = 1: : : ; 4: : : e A E (A; 1; 19) =6: : In an analogous way we have sin() + cos() cos() cos() sin A = sin(1) + sin() + cos() sin(1) cos() cos() sin(1) + sin() sin(1) sin() sin() 0: : : = 0: : : : 0: : : Taking n 0 = 8 one gets the Hermite matrix polynomial approximation with S (A; 1; 8) = 1 8 ( 1) n e (n + 1)! H n+1 (1; 1 ) A n=0 0: : : = 0: : : ; 0: : : sin A S(A; 1; 8) =4: :
10 114 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) The number of terms provided by (3) for the accuracy =10 5 is n 0 = 19. Finally approximating cos A by Hermite matrix polynomial series we have cos() sin() sin() sin() cos A = cos(1) + cos() sin() cos(1) + sin() sin() cos(1) + cos() cos(1) cos() cos() 1: : : = 1: : : : 0: : : Taking n 0 = 8 one gets the approximation C (A; 1; 8) = 1 8 ( 1) n e (n)! H n (1; 1 ) A n=0 1: : : = 1: : : ; 0: : : with error cos A C (A; 1; 8) =9: : For a prexed accuracy =10 5, the expression () gives n 0 = Applications In this section we construct matrix polynomial approximations of problems (1) (3) expressed in terms of Hermite matrix polynomials. It is well known that the solution of problem (1) is Y (x)=e Ax y 0 : Let be as in Lemma 3.1, 0, c 0, a = c +. Let n 0 be the rst positive integer such that A((A + I)a; n 0 ) exp(a A + I +1) e(1 + y 0 ) : Then by Corollary 3.1 it follows that Thus e Ax y 0 e x E(A + I;x;n)y 0 ; x c; n n 0 : Ỹ n (x)=e (1 x) [ H k (x; 1(A + ] I) ) y 0 ; (34) k!
11 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) is an approximate solution of problem (1) such that if Y (x) is the exact solution one gets Y (x) Ỹn(x) ; x c; n n 0 : (35) Problem () can be solved considering the extended system [ ] [ ] [ ] Y 0 I P Z = Y ; Z = A Z; Z(0) = ; 0 Q but such an approach increases the computational cost [1], and involves a lack of explicitness in terms of two vector parameters that is very interesting to study boundary value problems associated to () using the shooting method [9]. Consider problem () where A is an invertible matrix in C r r. By [6] the pair {cos(ax); sin(ax)} is a fundamental set of solutions of the equation Y +A Y =0; x + ; because Y 1 (x) = cos(ax), Y (x) = sin(ax) satisfy [ ] [ ] Y1 (0) Y W (0) = (0) I 0 Y 1(0) Y (0) = is invertible in C 0 A r r : Hence, the unique solution of problem () is given by Y (x) = cos(ax)p+ sin(ax)a 1 Q; x + : (36) If is given by Lemma 3.1, the expresion (36) and Corollary 3.1 suggest the approximation Ỹ n (x)={cos(x)c (A + I;x;n) + sin(x) (A + I;x;n)} P + {cos(x)s (A + I;x;n) sin(x)c (A + I;x;n)} A 1 Q = C (A + I;x;n) ( cos(x)p sin(x)a 1 Q ) +S (A + I;x;n) ( sin(x)p+ cos(x)a 1 Q ) ; [ ]( ( 1) k H k (x; 1 Ỹ n (x)= (A+I) ) cos(x) P ) (k)! e Q sin(x)a 1 e [ ]( ( 1) k H k+1 (x; 1 + (A+I) ) sin(x) P ) (k + 1)! e + Q cos(x)a 1 : (37) e By Corollary 3.1, for n max{n 1 ;n }, where n 1 and n are given by the Corollary 3.1, one gets Y (x) Ỹn(x) 6( P + A 1 Q ); x c: (38) We conclude this section with the construction of Hermite matrix polynomials approximations of the solution of problem (3). By [, p. 195] the solution of (3) is given by X (t)=e At Ce Bt : (39)
12 116 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) Let max{ A ; B } where and A max{ Re(z)+ Im(z) ; z (A)}; B max{ Re(z)+ Im(z) ; z (B)}: Corollary 3.1 suggests the approximation X n (t)=e t E(A + I;t;n)Ce t E(B + I;t;n); (40) and note that X (t) X n (t) =(e At e t E(A + I;t;n))Ce Bt +e t E(B + I;t;n)C(e Bt e t E(B + I;t;n)): (41) Consider the domain t c and take K = max{c( B +); }+1: By (15) and (18) it follows that t K=( B +) and E(B + I;t;n) 6e c +1 (K) j = L: (43) j! j 0 Taking norms in (41) and using (43) one gets X (t) X n (t) 6 e At e t E(A + I;t;n) C e c B + L C e Bt Taking the rst positive integer n 0 such that (4) e t E(B + I;t;n) : (44) e c B A((A+I)a; n 0) exp(w A+I a+1) e( C +1) ; and the rst positive integer n 1 such that A((B + I)a; n 1) exp( B + I a +1) L( C +1)e ; then for n max{n 0;n 1} and t c, by Corollary 3.1 and (39), (40), (43) it follows that X (t) X n (t) ; t c: References [1] T. Apostol, Explicit formulas for solutions of the second-order matrix dierential equation Y = AY, Amer. Math. Monthly 8 (1975) [] R. Bellman, Introduction to Matrix Analysis, McGraw-Hill, New York, 1965.
13 E. Defez, L. Jodar / Journal of Computational and Applied Mathematics 99 (1998) [3] N. Dunford, J. Schwartz, Linear Operators, Part I, Interscience, New York, [4] Z. Gajic, M. Qureshi, Lyapunov Matrix Equationin Systems Stability and Control, Academic Press, New York, [5] G. Golub, C.F. Van Loan, Matrix Computations, Johns Hopkins Univ. Press, Baltimore, M.A, [6] L. Jodar, Explicit solutions for second order operator dierential equations with two boundary value conditions, Linear Algebra Appl. 103 (1988) [7] L. Jodar, R. Company, Hermite matrix polynomials and second order matrix dierential equations, J. Approx. Theory Appl. 1 () (1996) [8] L. Jodar, E. Defez, On Hermite matrix polynomials and Hermite matrix functions, J. Approx. Theory Appl. 14 (1) (1998) [9] H.B. Keller, Numerical Solution of Two Point Boundary Value Problems, SIAM Regional Conf. Series in Appl. Maths., vol. 4, Philadelphia, [10] A.G. Law, C.N. Zhang, A. Rezazadeh, L. Jodar, Evaluation of a rational function, Numer. Algorithms 3 (199) [11] C.B. Moler, C.F. Van Loan, Nineteen dubious ways to compute the exponential of a matrix, SIAM Rev. 0 (1978) [1] E.D. Rainville, Special Functions, Chelsea, New York, [13] S. Saks, A. Zygmund, Analytic Functions, Elsevier, Amsterdam, [14] S.M. Serbin, Rational approximations of trigonometric matrices with applications to second-order systems of dierential equations, Appl. Math. Comput. 5 (1979) [15] S.M. Serbin, S.A. Blalock, An algorithm for computing the matrix cosine, SIAM J. Sci. Statist. Comput. 1 () (1980) [16] S. Serbin, C. Serbin, A time-stepping procedure for X = A 1X + XA + D; X (0) = C, IEEE Trans. Autom. Control. 5 (1980) [17] R.S. Varga, On higher order stable implicit methods for solving parabolic partial dierential equations, J. Math. Phys. 40 (1961) [18] R.C. Ward, Numerical computation of the matrix exponential with accuracy estimate, SIAM J. Numer. Anal. 14 (1977)
On Tricomi and Hermite-Tricomi Matrix Functions of Complex Variable
Communications in Mathematics and Applications Volume (0), Numbers -3, pp. 97 09 RGN Publications http://www.rgnpublications.com On Tricomi and Hermite-Tricomi Matrix Functions of Complex Variable A. Shehata
More informationBESSEL MATRIX DIFFERENTIAL EQUATIONS: EXPLICIT SOLUTIONS OF INITIAL AND TWO-POINT BOUNDARY VALUE PROBLEMS
APPLICATIONES MATHEMATICAE 22,1 (1993), pp. 11 23 E. NAVARRO, R. COMPANY and L. JÓDAR (Valencia) BESSEL MATRIX DIFFERENTIAL EQUATIONS: EXPLICIT SOLUTIONS OF INITIAL AND TWO-POINT BOUNDARY VALUE PROBLEMS
More informationSome properties of Hermite matrix polynomials
JOURNAL OF THE INTERNATIONAL MATHEMATICAL VIRTUAL INSTITUTE ISSN p) 33-4866, ISSN o) 33-4947 www.imvibl.org /JOURNALS / JOURNAL Vol. 515), 1-17 Former BULLETIN OF THE SOCIETY OF MATHEMATICIANS BANJA LUKA
More informationMatrix Analogues of Some Properties for Bessel Matrix Functions
J. Math. Sci. Univ. Tokyo (15, 519 53. Matrix Analogues of Some Properties for Bessel Matrix Functions By Bayram Çek ım and Abdullah Altın Abstract. In this study, we obtain some recurrence relations for
More informationComputing the matrix cosine
Numerical Algorithms 34: 13 26, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. Computing the matrix cosine Nicholas J. Higham and Matthew I. Smith Department of Mathematics, University
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationLecture 2: Computing functions of dense matrices
Lecture 2: Computing functions of dense matrices Paola Boito and Federico Poloni Università di Pisa Pisa - Hokkaido - Roma2 Summer School Pisa, August 27 - September 8, 2018 Introduction In this lecture
More informationContents. 6 Systems of First-Order Linear Dierential Equations. 6.1 General Theory of (First-Order) Linear Systems
Dierential Equations (part 3): Systems of First-Order Dierential Equations (by Evan Dummit, 26, v 2) Contents 6 Systems of First-Order Linear Dierential Equations 6 General Theory of (First-Order) Linear
More informationConstrained Leja points and the numerical solution of the constrained energy problem
Journal of Computational and Applied Mathematics 131 (2001) 427 444 www.elsevier.nl/locate/cam Constrained Leja points and the numerical solution of the constrained energy problem Dan I. Coroian, Peter
More informationImproved Newton s method with exact line searches to solve quadratic matrix equation
Journal of Computational and Applied Mathematics 222 (2008) 645 654 wwwelseviercom/locate/cam Improved Newton s method with exact line searches to solve quadratic matrix equation Jian-hui Long, Xi-yan
More informationExponentials of Symmetric Matrices through Tridiagonal Reductions
Exponentials of Symmetric Matrices through Tridiagonal Reductions Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A simple and efficient numerical algorithm
More informationElsevier Editorial System(tm) for Applied. Title: On the polynomial approximation of matrix functions
Mathematics and Computation Elsevier Editorial System(tm) for Applied Manuscript Draft Manuscript Number: Title: On the polynomial approximation of matrix functions Article Type: Full Length Article Keywords:
More informationOn the convergence of interpolatory-type quadrature rules for evaluating Cauchy integrals
Journal of Computational and Applied Mathematics 49 () 38 395 www.elsevier.com/locate/cam On the convergence of interpolatory-type quadrature rules for evaluating Cauchy integrals Philsu Kim a;, Beong
More informationGegenbauer Matrix Polynomials and Second Order Matrix. differential equations.
Gegenbauer Matrix Polynomials and Second Order Matrix Differential Equations Polinomios Matriciales de Gegenbauer y Ecuaciones Diferenciales Matriciales de Segundo Orden K. A. M. Sayyed, M. S. Metwally,
More information290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f
Numer. Math. 67: 289{301 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Least supported bases and local linear independence J.M. Carnicer, J.M. Pe~na? Departamento de Matematica
More informationA quadrature rule of interpolatory type for Cauchy integrals
Journal of Computational and Applied Mathematics 126 (2000) 207 220 www.elsevier.nl/locate/cam A quadrature rule of interpolatory type for Cauchy integrals Philsu Kim, U. Jin Choi 1 Department of Mathematics,
More informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More informationPolyGamma Functions of Negative Order
Carnegie Mellon University Research Showcase @ CMU Computer Science Department School of Computer Science -998 PolyGamma Functions of Negative Order Victor S. Adamchik Carnegie Mellon University Follow
More informationA NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS. Masaru Nagisa. Received May 19, 2014 ; revised April 10, (Ax, x) 0 for all x C n.
Scientiae Mathematicae Japonicae Online, e-014, 145 15 145 A NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS Masaru Nagisa Received May 19, 014 ; revised April 10, 014 Abstract. Let f be oeprator monotone
More informationTaylor series based nite dierence approximations of higher-degree derivatives
Journal of Computational and Applied Mathematics 54 (3) 5 4 www.elsevier.com/locate/cam Taylor series based nite dierence approximations of higher-degree derivatives Ishtiaq Rasool Khan a;b;, Ryoji Ohba
More informationNP-hardness of the stable matrix in unit interval family problem in discrete time
Systems & Control Letters 42 21 261 265 www.elsevier.com/locate/sysconle NP-hardness of the stable matrix in unit interval family problem in discrete time Alejandra Mercado, K.J. Ray Liu Electrical and
More informationMoore Penrose inverses and commuting elements of C -algebras
Moore Penrose inverses and commuting elements of C -algebras Julio Benítez Abstract Let a be an element of a C -algebra A satisfying aa = a a, where a is the Moore Penrose inverse of a and let b A. We
More informationOn some properties of elementary derivations in dimension six
Journal of Pure and Applied Algebra 56 (200) 69 79 www.elsevier.com/locate/jpaa On some properties of elementary derivations in dimension six Joseph Khoury Department of Mathematics, University of Ottawa,
More informationTwo Results About The Matrix Exponential
Two Results About The Matrix Exponential Hongguo Xu Abstract Two results about the matrix exponential are given. One is to characterize the matrices A which satisfy e A e AH = e AH e A, another is about
More informationBOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION
K Y BERNETIKA VOLUM E 46 ( 2010), NUMBER 4, P AGES 655 664 BOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION Guang-Da Hu and Qiao Zhu This paper is concerned with bounds of eigenvalues of a complex
More informationSTABILITY OF INVARIANT SUBSPACES OF COMMUTING MATRICES We obtain some further results for pairs of commuting matrices. We show that a pair of commutin
On the stability of invariant subspaces of commuting matrices Tomaz Kosir and Bor Plestenjak September 18, 001 Abstract We study the stability of (joint) invariant subspaces of a nite set of commuting
More informationLinearized methods for ordinary di erential equations
Applied Mathematics and Computation 104 (1999) 109±19 www.elsevier.nl/locate/amc Linearized methods for ordinary di erential equations J.I. Ramos 1 Departamento de Lenguajes y Ciencias de la Computacion,
More informationMATHEMATICAL FORMULAS AND INTEGRALS
MATHEMATICAL FORMULAS AND INTEGRALS ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom Academic Press San Diego New York Boston London
More informationPositive Denite Matrix. Ya Yan Lu 1. Department of Mathematics. City University of Hong Kong. Kowloon, Hong Kong. Abstract
Computing the Logarithm of a Symmetric Positive Denite Matrix Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A numerical method for computing the logarithm
More informationHEINEMANN HIGHER CHECKLIST
St Ninian s High School HEINEMANN HIGHER CHECKLIST I understand this part of the course = I am unsure of this part of the course = Name Class Teacher I do not understand this part of the course = Topic
More informationNonsingularity and group invertibility of linear combinations of two k-potent matrices
Nonsingularity and group invertibility of linear combinations of two k-potent matrices Julio Benítez a Xiaoji Liu b Tongping Zhu c a Departamento de Matemática Aplicada, Instituto de Matemática Multidisciplinar,
More informationRobust linear optimization under general norms
Operations Research Letters 3 (004) 50 56 Operations Research Letters www.elsevier.com/locate/dsw Robust linear optimization under general norms Dimitris Bertsimas a; ;, Dessislava Pachamanova b, Melvyn
More informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More informationDiscrete Orthogonal Polynomials on Equidistant Nodes
International Mathematical Forum, 2, 2007, no. 21, 1007-1020 Discrete Orthogonal Polynomials on Equidistant Nodes Alfredo Eisinberg and Giuseppe Fedele Dip. di Elettronica, Informatica e Sistemistica Università
More informationPROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS
PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS Abstract. We present elementary proofs of the Cauchy-Binet Theorem on determinants and of the fact that the eigenvalues of a matrix
More informationAbsolutely indecomposable symmetric matrices
Journal of Pure and Applied Algebra 174 (2002) 83 93 wwwelseviercom/locate/jpaa Absolutely indecomposable symmetric matrices Hans A Keller a; ;1, A Herminia Ochsenius b;1 a Hochschule Technik+Architektur
More informationMATHEMATICAL FORMULAS AND INTEGRALS
HANDBOOK OF MATHEMATICAL FORMULAS AND INTEGRALS Second Edition ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom ACADEMIC PRESS A Harcourt
More informationSpurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics
More informationA note on continuous behavior homomorphisms
Available online at www.sciencedirect.com Systems & Control Letters 49 (2003) 359 363 www.elsevier.com/locate/sysconle A note on continuous behavior homomorphisms P.A. Fuhrmann 1 Department of Mathematics,
More informationChapter 7: Techniques of Integration
Chapter 7: Techniques of Integration MATH 206-01: Calculus II Department of Mathematics University of Louisville last corrected September 14, 2013 1 / 43 Chapter 7: Techniques of Integration 7.1. Integration
More informationInterval solutions for interval algebraic equations
Mathematics and Computers in Simulation 66 (2004) 207 217 Interval solutions for interval algebraic equations B.T. Polyak, S.A. Nazin Institute of Control Sciences, Russian Academy of Sciences, 65 Profsoyuznaya
More informationMATH 425-Spring 2010 HOMEWORK ASSIGNMENTS
MATH 425-Spring 2010 HOMEWORK ASSIGNMENTS Instructor: Shmuel Friedland Department of Mathematics, Statistics and Computer Science email: friedlan@uic.edu Last update April 18, 2010 1 HOMEWORK ASSIGNMENT
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationNonhomogeneous Linear Differential Equations with Constant Coefficients - (3.4) Method of Undetermined Coefficients
Nonhomogeneous Linear Differential Equations with Constant Coefficients - (3.4) Method of Undetermined Coefficients Consider an nth-order nonhomogeneous linear differential equation with constant coefficients:
More informationSample Quiz 8, Problem 1. Solving Higher Order Constant-Coefficient Equations
Sample Quiz 8 Sample Quiz 8, Problem. Solving Higher Order Constant-Coefficient Equations The Algorithm applies to constant-coefficient homogeneous linear differential equations of order N, for example
More informationApproximation theory
Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2
More informationA scale-free goodness-of-t statistic for the exponential distribution based on maximum correlations
Journal of Statistical Planning and Inference 18 (22) 85 97 www.elsevier.com/locate/jspi A scale-free goodness-of-t statistic for the exponential distribution based on maximum correlations J. Fortiana,
More informationBounds for Variable Degree Rational L Approximations to the Matrix Cosine
Bounds for Variable Degree Rational L Approximations to the Matrix Cosine Ch. Tsitouras a, V.N. Katsikis b a TEI of Sterea Hellas, GR34400, Psahna, Greece b TEI of Piraeus, GR12244, Athens, Greece Abstract
More informationKrylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations
Krylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations Peter Benner and Tobias Breiten Abstract We discuss Krylov-subspace based model reduction
More informationError formulas for divided difference expansions and numerical differentiation
Error formulas for divided difference expansions and numerical differentiation Michael S. Floater Abstract: We derive an expression for the remainder in divided difference expansions and use it to give
More informationFunctions: A Fourier Approach. Universitat Rostock. Germany. Dedicated to Prof. L. Berg on the occasion of his 65th birthday.
Approximation Properties of Multi{Scaling Functions: A Fourier Approach Gerlind Plona Fachbereich Mathemati Universitat Rostoc 1851 Rostoc Germany Dedicated to Prof. L. Berg on the occasion of his 65th
More informationFree download from not for resale. Apps 1.1 : Applying algebraic skills to rectilinear shapes.
Free download from, not for resale. Apps 1.1 : Applying algebraic skills to rectilinear shapes. Gradients m = tanθ Distance Formula Midpoint Formula Parallel lines Perpendicular lines y = mx + c y - b
More informationChapter 1 Numerical approximation of data : interpolation, least squares method
Chapter 1 Numerical approximation of data : interpolation, least squares method I. Motivation 1 Approximation of functions Evaluation of a function Which functions (f : R R) can be effectively evaluated
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationConstrained controllability of semilinear systems with delayed controls
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 56, No. 4, 28 Constrained controllability of semilinear systems with delayed controls J. KLAMKA Institute of Control Engineering, Silesian
More information1 Holomorphic functions
Robert Oeckl CA NOTES 1 15/09/2009 1 1 Holomorphic functions 11 The complex derivative The basic objects of complex analysis are the holomorphic functions These are functions that posses a complex derivative
More informationarxiv: v1 [math.ca] 31 Dec 2018
arxiv:181.1173v1 [math.ca] 31 Dec 18 Some trigonometric integrals and the Fourier transform of a spherically symmetric exponential function Hideshi YAMANE Department of Mathematical Sciences, Kwansei Gakuin
More informationDiscrete Halanay-type inequalities and applications
Nonlinear Analysis 55 (2003) 669 678 www.elsevier.com/locate/na Discrete Halanay-type inequalities and applications Eduardo Liz a;, Anatoli Ivanov b, Juan Bosco Ferreiro c a Departamento de Matematica
More informationThe generalized Schur complement in group inverses and (k + 1)-potent matrices
The generalized Schur complement in group inverses and (k + 1)-potent matrices Julio Benítez Néstor Thome Abstract In this paper, two facts related to the generalized Schur complement are studied. The
More informationCOMMUTATIVITY OF OPERATORS
COMMUTATIVITY OF OPERATORS MASARU NAGISA, MAKOTO UEDA, AND SHUHEI WADA Abstract. For two bounded positive linear operators a, b on a Hilbert space, we give conditions which imply the commutativity of a,
More informationOn matrix equations X ± A X 2 A = I
Linear Algebra and its Applications 326 21 27 44 www.elsevier.com/locate/laa On matrix equations X ± A X 2 A = I I.G. Ivanov,V.I.Hasanov,B.V.Minchev Faculty of Mathematics and Informatics, Shoumen University,
More informationMcGill University Math 325A: Differential Equations LECTURE 12: SOLUTIONS FOR EQUATIONS WITH CONSTANTS COEFFICIENTS (II)
McGill University Math 325A: Differential Equations LECTURE 12: SOLUTIONS FOR EQUATIONS WITH CONSTANTS COEFFICIENTS (II) HIGHER ORDER DIFFERENTIAL EQUATIONS (IV) 1 Introduction (Text: pp. 338-367, Chap.
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K. R. MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 1998 Comments to the author at krm@maths.uq.edu.au Contents 1 LINEAR EQUATIONS
More informationLegendre s Equation. PHYS Southern Illinois University. October 18, 2016
Legendre s Equation PHYS 500 - Southern Illinois University October 18, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 18, 2016 1 / 11 Legendre s Equation Recall We are trying
More informationThe reflexive and anti-reflexive solutions of the matrix equation A H XB =C
Journal of Computational and Applied Mathematics 200 (2007) 749 760 www.elsevier.com/locate/cam The reflexive and anti-reflexive solutions of the matrix equation A H XB =C Xiang-yang Peng a,b,, Xi-yan
More informationApplications of CAS to analyze the step response of a system with parameters
Applications of CAS to analyze the step response of a system with parameters Takuya Kitamoto kitamoto@yamaguchi-u.ac.jp Faculty of Education Yamaguchi University 753-8513 Japan Abstract Recently, Computer
More informationMatrix Exponential Formulas
Matrix Exponential Formulas Linear Analysis May 3, 016 Abstract We present some tricks for quickly calculating the matrix exponential of certain special classes of matrices, such as matrices. 1 matrix
More information2 Computing complex square roots of a real matrix
On computing complex square roots of real matrices Zhongyun Liu a,, Yulin Zhang b, Jorge Santos c and Rui Ralha b a School of Math., Changsha University of Science & Technology, Hunan, 410076, China b
More informationALGEBRA AND TRIGONOMETRY
ALGEBRA AND TRIGONOMETRY correlated to the Alabama Course of Study for Algebra 2 with Trigonometry CC2 6/2003 2004 Algebra and Trigonometry 2004 correlated to the Number and Operations Students will: 1.
More informationIntroductory Analysis I Fall 2014 Homework #9 Due: Wednesday, November 19
Introductory Analysis I Fall 204 Homework #9 Due: Wednesday, November 9 Here is an easy one, to serve as warmup Assume M is a compact metric space and N is a metric space Assume that f n : M N for each
More informationChapter 8 Integration Techniques and Improper Integrals
Chapter 8 Integration Techniques and Improper Integrals 8.1 Basic Integration Rules 8.2 Integration by Parts 8.4 Trigonometric Substitutions 8.5 Partial Fractions 8.6 Numerical Integration 8.7 Integration
More informationAlgorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations
dmitry.lyakhov@kaust.edu.sa Symbolic Computations May 4, 2018 1 / 25 Algorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations Dmitry A. Lyakhov 1 1 Computational Sciences
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationON THE HÖLDER CONTINUITY OF MATRIX FUNCTIONS FOR NORMAL MATRICES
Volume 10 (2009), Issue 4, Article 91, 5 pp. ON THE HÖLDER CONTINUITY O MATRIX UNCTIONS OR NORMAL MATRICES THOMAS P. WIHLER MATHEMATICS INSTITUTE UNIVERSITY O BERN SIDLERSTRASSE 5, CH-3012 BERN SWITZERLAND.
More informationMATH NEW HOMEWORK AND SOLUTIONS TO PREVIOUS HOMEWORKS AND EXAMS
MATH. 4433. NEW HOMEWORK AND SOLUTIONS TO PREVIOUS HOMEWORKS AND EXAMS TOMASZ PRZEBINDA. Final project, due 0:00 am, /0/208 via e-mail.. State the Fundamental Theorem of Algebra. Recall that a subset K
More informationav 1 x 2 + 4y 2 + xy + 4z 2 = 16.
74 85 Eigenanalysis The subject of eigenanalysis seeks to find a coordinate system, in which the solution to an applied problem has a simple expression Therefore, eigenanalysis might be called the method
More informationAn explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary
Journal of Computational and Applied Mathematics 156 (2003) 319 343 www.elsevier.com/locate/cam An explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary
More information8.5 Taylor Polynomials and Taylor Series
8.5. TAYLOR POLYNOMIALS AND TAYLOR SERIES 50 8.5 Taylor Polynomials and Taylor Series Motivating Questions In this section, we strive to understand the ideas generated by the following important questions:
More informationProblem Set 9 Due: In class Tuesday, Nov. 27 Late papers will be accepted until 12:00 on Thursday (at the beginning of class).
Math 3, Fall Jerry L. Kazdan Problem Set 9 Due In class Tuesday, Nov. 7 Late papers will be accepted until on Thursday (at the beginning of class).. Suppose that is an eigenvalue of an n n matrix A and
More informationRational Trigonometry. Rational Trigonometry
There are an infinite number of points on the unit circle whose coordinates are each rational numbers. Indeed every Pythagorean triple gives one! 3 65 64 65 7 5 4 5 03 445 396 445 3 4 5 5 75 373 5 373
More informationSome suggested repetition for the course MAA508
Some suggested repetition for the course MAA58 Linus Carlsson, Karl Lundengård, Johan Richter July, 14 Contents Introduction 1 1 Basic algebra and trigonometry Univariate calculus 5 3 Linear algebra 8
More informationMS 3011 Exercises. December 11, 2013
MS 3011 Exercises December 11, 2013 The exercises are divided into (A) easy (B) medium and (C) hard. If you are particularly interested I also have some projects at the end which will deepen your understanding
More informationADDITIONAL MATHEMATICS
ADDITIONAL MATHEMATICS GCE NORMAL ACADEMIC LEVEL (016) (Syllabus 4044) CONTENTS Page INTRODUCTION AIMS ASSESSMENT OBJECTIVES SCHEME OF ASSESSMENT 3 USE OF CALCULATORS 3 SUBJECT CONTENT 4 MATHEMATICAL FORMULAE
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 998 Comments to the author at krm@mathsuqeduau All contents copyright c 99 Keith
More informationThe restarted QR-algorithm for eigenvalue computation of structured matrices
Journal of Computational and Applied Mathematics 149 (2002) 415 422 www.elsevier.com/locate/cam The restarted QR-algorithm for eigenvalue computation of structured matrices Daniela Calvetti a; 1, Sun-Mi
More informationNew series expansions of the Gauss hypergeometric function
New series expansions of the Gauss hypergeometric function José L. López and Nico. M. Temme 2 Departamento de Matemática e Informática, Universidad Pública de Navarra, 36-Pamplona, Spain. e-mail: jl.lopez@unavarra.es.
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationPolynomial functions over nite commutative rings
Polynomial functions over nite commutative rings Balázs Bulyovszky a, Gábor Horváth a, a Institute of Mathematics, University of Debrecen, Pf. 400, Debrecen, 4002, Hungary Abstract We prove a necessary
More informationThe symmetric linear matrix equation
Electronic Journal of Linear Algebra Volume 9 Volume 9 (00) Article 8 00 The symmetric linear matrix equation Andre CM Ran Martine CB Reurings mcreurin@csvunl Follow this and additional works at: http://repositoryuwyoedu/ela
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationApplied Mathematics Letters. Comparison theorems for a subclass of proper splittings of matrices
Applied Mathematics Letters 25 (202) 2339 2343 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Comparison theorems for a subclass
More informationMath Exam 2, October 14, 2008
Math 96 - Exam 2, October 4, 28 Name: Problem (5 points Find all solutions to the following system of linear equations, check your work: x + x 2 x 3 2x 2 2x 3 2 x x 2 + x 3 2 Solution Let s perform Gaussian
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationCOMPLEX ANALYSIS HW 3
COMPLEX ANALYSIS HW 3 CLAY SHONKWILER 6 Justify the swap of limit and integral carefully to prove that for each function f : D C continuous on an open set D C, and for each differentiable curve : [0, D,
More informationEXAMPLES OF PROOFS BY INDUCTION
EXAMPLES OF PROOFS BY INDUCTION KEITH CONRAD 1. Introduction In this handout we illustrate proofs by induction from several areas of mathematics: linear algebra, polynomial algebra, and calculus. Becoming
More informationSection 8.7. Taylor and MacLaurin Series. (1) Definitions, (2) Common Maclaurin Series, (3) Taylor Polynomials, (4) Applications.
Section 8.7 Taylor and MacLaurin Series (1) Definitions, (2) Common Maclaurin Series, (3) Taylor Polynomials, (4) Applications. MATH 126 (Section 8.7) Taylor and MacLaurin Series The University of Kansas
More informationA trigonometric orthogonality with respect to a nonnegative Borel measure
Filomat 6:4 01), 689 696 DOI 10.98/FIL104689M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat A trigonometric orthogonality with
More informationAn FKG equality with applications to random environments
Statistics & Probability Letters 46 (2000) 203 209 An FKG equality with applications to random environments Wei-Shih Yang a, David Klein b; a Department of Mathematics, Temple University, Philadelphia,
More informationPositive definite solutions of some matrix equations
Available online at www.sciencedirect.com Linear Algebra and its Applications 49 (008) 40 44 www.elsevier.com/locate/laa Positive definite solutions of some matrix equations Aleksandar S. Cvetković, Gradimir
More informationChapter 2: Differentiation
Chapter 2: Differentiation Winter 2016 Department of Mathematics Hong Kong Baptist University 1 / 75 2.1 Tangent Lines and Their Slopes This section deals with the problem of finding a straight line L
More information