TRIPLE FACTORIZATION OF SOME RIORDAN MATRICES. Paul Peart* Department of Mathematics, Howard University, Washington, D.C

Size: px
Start display at page:

Download "TRIPLE FACTORIZATION OF SOME RIORDAN MATRICES. Paul Peart* Department of Mathematics, Howard University, Washington, D.C"

Transcription

1 Pau Peart* Department of Mathematics, Howard University, Washington, D.C. 59 Leon Woodson Department of Mathematicss, Morgan State University, Batimore, MD 39 (Submitted June 99). INTRODUCTION When examining a combinatoria sequence, generating functions are often usefu. That is, if we are interested in anayzing the sequence a, a? a,..., we investigate the forma power series f(x) = a +a x + a x -\. In a recent paper [], techniques are discussed that assist in finding cosed-form expressions for the forma power series for a seect, but arge, set of combinatoria sequences. The methods invove using infinite matrices and the Riordan group. The Riordan group is defined in section of this paper. Each matrix, L, in the Riordan group is associated with a combinatoria sequence and with a matrix, S L, caed the Stietjes matrix. S L is defined in section 3. In this paper, we show that when S L is tridiagona, then L= PCF, where the first factor P is a Pasca-type matrix, the second factor C invoves the generating function for the Cataan numbers, and the third factor F invoves the Fibonacci generating function. The foowing is an exampe: I " o v r i - i o : - 5. r i The matrices in the Riordan group are infinite and ower trianguar. So the exampe shows ony the first seven rows. The first factor on the right is the Pasca matrix. The first coumn in the second factor has C(-x ) as generating function, where is the generating function for the Cataan numbers. The third factor has the Fibonacci numbers in each coumn. See section 6 for further exampes of this tripe factorization. * Partiay supported by NSF grant ]

2 In section, we define the Riordan group R and ist some properties that we use in the proofs of the propositions which are given in section 4. In section 3, we discuss the unique Stietjes matrix S L associated with each L in the Riordan group. In this paper, we concentrate on the subset of R given by R T - {L GR:S L is tridiagona}. In section 5, we derive a recurrence reation for the sequence associated with each member of R T, and we discuss the asymptotic behavior of these sequences. In section 6, we provide two exampes invoving we-known sequences. For each exampe, we give the tripe factorization, the Stietjes matrix, the recurrence reation and asymptotic behavior of the corresponding sequence.. THE RIORDAN GROUP A detaied description of this group is given in []. Here we provide a brief summary. Let M = (rnjj)jj>o be an infinite matrix with eements from C, the set of compex numbers. Let c f (x) be the generating function of the j * coumn ofm That is <*(*) = m n,i*"' We ca M a Riordan matrix if c t (x) - g(x)[/(x)j, where g(x) = + g x + g x +g 3 x 3 + -~, and f(x) = x + f x +f 3 x 3 + -~. In this case, we write M = (g(x), /(*)). We denote by R the set of Riordan matrices. R is a group under matrix mutipication with the foowing properties: ( (*<*), fix)) * (h(x), t(x)) = {g(x)h(f(x)), /(/(x))). (ii) / = (, x) is the identity eement. (iii) The inverse of Mis given by where / is the compositiona inverse of/ M~ = > / ( * ) jk/(*)) (iv) If (a,a y a,...) T is a coumn vector with generating function A(x), then mutipying M = (g(x\ /(*)) on the right by this coumn vector yieds a coumn vector with generating function B(x) = g(x)a(f(x)). 3. STIELTJES MATRI Let L be Riordan and et L be the matrix obtained from L by deeting the first row. For exampe, if I is the identity, we have / =. :. [MAY

3 Observe that L = IL. There exists a unique matrix, S L, such that LS L = L. We ca this matrix the Stietjes matrix of L. Exampe: If Z = x _ -x' -x) then s,= ". 4. PROPOSITIONS Proposition : If L = (g(x), f(x)) is Riordan and S L is tridiagona, then (a) S L=\ h, b b b b. ) / = jc(+ft/ + / ) and g = iff JS^ is as in (a). - Z> - yf ivoq/; Let S L = b y by A b 3 h. With c,(x) the generating function for the z* coumn of L, / >, we have c, = g/"'. By ooking at the first coumn of LS L and L, we obtain b xg + x xgf = g -, i.e., For / >, we obtain from LS L = L, *(*) = -h x- x xf 993] 3

4 c i =x{c i _ +b i c i + i+] c M ). o f = x(+b i f + A i+ f ).: = (ft, - bj)f + ( i+ - J+ )f for a i and; >.. -. ft, = bj and j+ = J+ for a / and; >..". we can take ft = ftj = ft = ft 3 = and A A = yj A^ =..-. / = x( + ft/ + ) - Remark: If 5^ is tridiagona, it has the form in (a) and then either (a) = and/ = -bx and# = -bx - (Z>i +b )x + (bb - Jx or * «A* -bx-j(b -4)x -bx +, (b) * and/ = and g =. x - b x - {xf Proposition : If L = (g, / ) is Riordan ; then S L = S L *+b if and ony ifl = P b L*, where - bx' - bx ft b ft ft 3 3ft 3b ft 4 4ft 3 6ft 4ft (cf. [3], p. 7) Proof: Note that v = ft ft ft ft. = b+. So, Z = P fr Z*^>/L = /P*Z* Conversey, suppose S L =b + $ L *. Then => Z = P*Z* = P'Sp* I* = P b (b+i)l* = bp b L*+P b IL* = bl + P b L* = bl + P b L* S L * = L(b + S L *). P b L* (ft/ + S L,) = MP*!* = P b L* = P b bil*+p b IL* = P b (b + T)L* = P b L* = I(P b L*) = P b L* 4 [MAY

5 Proposition 3: If L = {g, / ) is Riordan and S L = b + s + S b b \ b A A then Z, = P 6^, where Su = s + S. ' o A Proof: This foows immediatey from Proposition. Proposition 4 (PCFFactorization): In Proposition 3, Z^ = C x F es, where Q=(c(^),xc(;ix )) with C() = + [C(JC)] = -V-4x ^, i x ^ v " y ^^^ir^?-* Proof: Let J^ = (ft,/!). Then, from Proposition, we must have, when ^, x -V-4^x a n d /i= ^ & = Ifoc \-sx-{ + S)xf f x = xc(ax ) and g x = -ac-(z + S)x c(jbc )' Now, from section, property, we have But C A, = (c(^vc(^)).( T -^,*) = ( T c</x ) -8-(A ^ / ^ _ + S)x z c(ax z ) \ i -axc(ax z )-Sx z [c(ax z )] c(ax ).\T 6xc(/?x )-<5[xc(/fo: )],xc(ax ) o> - sxc(ax ) - 8x [c{hc )] - c(ax )-zxc(ax )-(A + S)x [c(ax )] o-c(x ) + Ax [c(ax )f=. 993] 5

6 5. RECURRENCE RELATIONS AND ASYMPTOTICS We have proved that when L = (/(x), g{x)) is Riordan and S L is tridiagona with the form s L = h A, b b b \ b. then and /(*) = Jf xn = ~ fa ~^" 4^x " fa + w= /bc =,-v-v/w Using the J.C. P. Mier formua (see Henrici [3]), we obtain for f n the three-term recurrence and for g n the five-term recurrence (#i + )/ + = (n + )Z>/ + ( -»)(ft - 4)/ _, nag = [{In-3)bA-nB]g^ + [(In-3)bB + (3-n)(b -A)A-nC]g n _ + [(n-3)bc + (3-n)(b -4)B]g n _ 3 +[(3-n)(b -A)C]g _ 4 where A-- x, B- x b + x \ -b, C = [ - x bb Q + b^. For the asymptotics, we use the methods described in Wif [4, Ch. 5]. For arge n, we obtain where b > A >. fn (n + \y V n+/ (b+4) I^^n Because there are too many cases to consider, we do not attempt to provide a genera formua for the asymptotic vaue of g n. However, the exampes in section 6 iustrate the techniques invoved. 6. EAMPLES Exampe Big Schroder Numbers: If we take -, b - 3, x = + 8-, and b - b + s -, then j. - 3 x - v x -6x +, - x - v x - 6 x + f= and g =. 4x x g is the generating function for the Big Schroder numbers []. 6 [MAY

7 L = (g,f)with S L and L--p%F W --[^, T^-y {^),^)),[ T L-, x ) [ o : i. 8 4 o ;] i Recurrence Reations: Here A =, B =, C =, (n + )/ + = 3(» +)/ + ( - «)/ _! for n >. / =,/ =. ng _ = 3(n-3)g n _ +(3-n)g n _ 3, (n + )g = 3(n-)g _ + (-n)g _, for Asymptotics: For arge «, fn = [*"]/'(*)' (n + \y3 ' (b+vi)" +/ _ (» +)" 3/ (3 + >/ ) n+/ A 3/4 V^ - 3 / 4 - ^. - x - V x - 6 x + x &-I t^k*'-«r + )»./.-fil + >'"P + j 5 ' «+/ Exampe Leendre Poynomias: We require We take *oo= 4x r - =, b = t, y = + S 4 -tx + t - 3/4 V^ and b n =.b + e = t. 993] 7

8 Tripe Factorization: L = (-ft-vx -ft + ) V* -ft + ' (f -\)y (f -)/4,(f -)/4' t'-i s L = o s-i = t ^ / Recurrence Reations: Here A = ±g-,b= '-^-, C = ± -. (n + )/ + = (#i +)// +( -»)/ _!, for «>. / =,/, =.» g - -(4»-3)rg- _ + (3-»)(+r )^_ +(4»-9)^_3+(3-»)^_ 4, for «>4. g =,g = t,g =^t -^,g 3 =^-^. Asymptotics: We assume that 7 >, so that the roots of x - tx + = are rea. Denote these roots by f and r with f \ < \ r\. We obtain v/ [*"]/(*) = -p ^ix -ft +) - ^ (H + )" 3 7 ^ _ ^ / / -J(r)" + -V^ I f Vr)" (r) - r - w n + V «#" [x"]g(x) = [x"](x i \-/ -ft + > <r)' W+ V^((r) -) REFERENCES. L. Comtet. Advanced Combinatorics. D. Reide Pubishing Company, S. Getu, L. Shapiro, W. Woan, & L. Woodson. "The Riordan Group." Discrete Appied Mathematics 34 (99): P. Henrici. Appied and Computationa Compex Anaysis. Vo. I. New York: Wiey & Sons, H. Wif. Generatingfunctionoogy. New York: Academic Press, 989. AMS Cassification numbers: 5A5, A4 8 [MAY

ALGORITHMIC SUMMATION OF RECIPROCALS OF PRODUCTS OF FIBONACCI NUMBERS. F. = I j. ^ = 1 ^ -, and K w = ^. 0) n=l r n «=1 -*/!

ALGORITHMIC SUMMATION OF RECIPROCALS OF PRODUCTS OF FIBONACCI NUMBERS. F. = I j. ^ = 1 ^ -, and K w = ^. 0) n=l r n «=1 -*/! ALGORITHMIC SUMMATIO OF RECIPROCALS OF PRODUCTS OF FIBOACCI UMBERS Staney Rabinowitz MathPro Press, 2 Vine Brook Road, Westford, MA 0886 staney@tiac.net (Submitted May 997). ITRODUCTIO There is no known

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

arxiv: v1 [math.nt] 12 Feb 2019

arxiv: v1 [math.nt] 12 Feb 2019 Degenerate centra factoria numbers of the second ind Taeyun Kim, Dae San Kim arxiv:90.04360v [math.nt] Feb 09 In this paper, we introduce the degenerate centra factoria poynomias and numbers of the second

More information

TRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS

TRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS Vo. 39 (008) ACTA PHYSICA POLONICA B No 8 TRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS Zbigniew Romanowski Interdiscipinary Centre for Materias Modeing Pawińskiego 5a, 0-106 Warsaw, Poand

More information

Theory of Generalized k-difference Operator and Its Application in Number Theory

Theory of Generalized k-difference Operator and Its Application in Number Theory Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication

More information

Discrete Bernoulli s Formula and its Applications Arising from Generalized Difference Operator

Discrete Bernoulli s Formula and its Applications Arising from Generalized Difference Operator Int. Journa of Math. Anaysis, Vo. 7, 2013, no. 5, 229-240 Discrete Bernoui s Formua and its Appications Arising from Generaized Difference Operator G. Britto Antony Xavier 1 Department of Mathematics,

More information

w = X ^ = ^ ^, (i*i < ix

w = X ^ = ^ ^, (i*i < ix A SUMMATION FORMULA FOR POWER SERIES USING EULERIAN FRACTIONS* Xinghua Wang Math. Department, Zhejiang University, Hangzhou 310028, China Leetsch C. Hsu Math. Institute, Dalian University of Technology,

More information

Assignment 7 Due Tuessday, March 29, 2016

Assignment 7 Due Tuessday, March 29, 2016 Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges

More information

Matrices and Determinants

Matrices and Determinants Matrices and Determinants Teaching-Learning Points A matri is an ordered rectanguar arra (arrangement) of numbers and encosed b capita bracket [ ]. These numbers are caed eements of the matri. Matri is

More information

Statistics for Applications. Chapter 7: Regression 1/43

Statistics for Applications. Chapter 7: Regression 1/43 Statistics for Appications Chapter 7: Regression 1/43 Heuristics of the inear regression (1) Consider a coud of i.i.d. random points (X i,y i ),i =1,...,n : 2/43 Heuristics of the inear regression (2)

More information

Pascal Eigenspaces and Invariant Sequences of the First or Second Kind

Pascal Eigenspaces and Invariant Sequences of the First or Second Kind Pascal Eigenspaces and Invariant Sequences of the First or Second Kind I-Pyo Kim a,, Michael J Tsatsomeros b a Department of Mathematics Education, Daegu University, Gyeongbu, 38453, Republic of Korea

More information

Partial permutation decoding for MacDonald codes

Partial permutation decoding for MacDonald codes Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics

More information

1 Linear Algebra Problems

1 Linear Algebra Problems Linear Algebra Problems. Let A be the conjugate transpose of the complex matrix A; i.e., A = A t : A is said to be Hermitian if A = A; real symmetric if A is real and A t = A; skew-hermitian if A = A and

More information

Course 2BA1, Section 11: Periodic Functions and Fourier Series

Course 2BA1, Section 11: Periodic Functions and Fourier Series Course BA, 8 9 Section : Periodic Functions and Fourier Series David R. Wikins Copyright c David R. Wikins 9 Contents Periodic Functions and Fourier Series 74. Fourier Series of Even and Odd Functions...........

More information

Math 220B - Summer 2003 Homework 1 Solutions

Math 220B - Summer 2003 Homework 1 Solutions Math 0B - Summer 003 Homework Soutions Consider the eigenvaue probem { X = λx 0 < x < X satisfies symmetric BCs x = 0, Suppose f(x)f (x) x=b x=a 0 for a rea-vaued functions f(x) which satisfy the boundary

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

Generalized Bell polynomials and the combinatorics of Poisson central moments

Generalized Bell polynomials and the combinatorics of Poisson central moments Generaized Be poynomias and the combinatorics of Poisson centra moments Nicoas Privaut Division of Mathematica Sciences Schoo of Physica and Mathematica Sciences Nanyang Technoogica University SPMS-MAS-05-43,

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

CALCULUS JIA-MING (FRANK) LIOU

CALCULUS JIA-MING (FRANK) LIOU CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion

More information

Consistent linguistic fuzzy preference relation with multi-granular uncertain linguistic information for solving decision making problems

Consistent linguistic fuzzy preference relation with multi-granular uncertain linguistic information for solving decision making problems Consistent inguistic fuzzy preference reation with muti-granuar uncertain inguistic information for soving decision making probems Siti mnah Binti Mohd Ridzuan, and Daud Mohamad Citation: IP Conference

More information

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,

More information

THINKING IN PYRAMIDS

THINKING IN PYRAMIDS ECS 178 Course Notes THINKING IN PYRAMIDS Kenneth I. Joy Institute for Data Anaysis and Visuaization Department of Computer Science University of Caifornia, Davis Overview It is frequenty usefu to think

More information

Parts Manual. EPIC II Critical Care Bed REF 2031

Parts Manual. EPIC II Critical Care Bed REF 2031 EPIC II Critical Care Bed REF 2031 Parts Manual For parts or technical assistance call: USA: 1-800-327-0770 2013/05 B.0 2031-109-006 REV B www.stryker.com Table of Contents English Product Labels... 4

More information

On the New q-extension of Frobenius-Euler Numbers and Polynomials Arising from Umbral Calculus

On the New q-extension of Frobenius-Euler Numbers and Polynomials Arising from Umbral Calculus Adv. Studies Theor. Phys., Vo. 7, 203, no. 20, 977-99 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/0.2988/astp.203.390 On the New -Extension of Frobenius-Euer Numbers and Poynomias Arising from Umbra

More information

Lecture Notes 4: Fourier Series and PDE s

Lecture Notes 4: Fourier Series and PDE s Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for

More information

Algorithms to solve massively under-defined systems of multivariate quadratic equations

Algorithms to solve massively under-defined systems of multivariate quadratic equations Agorithms to sove massivey under-defined systems of mutivariate quadratic equations Yasufumi Hashimoto Abstract It is we known that the probem to sove a set of randomy chosen mutivariate quadratic equations

More information

is equal to = 3 2 x, if x < 0 f (0) = lim h = 0. Therefore f exists and is continuous at 0.

is equal to = 3 2 x, if x < 0 f (0) = lim h = 0. Therefore f exists and is continuous at 0. Madhava Mathematics Competition January 6, 2013 Solutions and scheme of marking Part I N.B. Each question in Part I carries 2 marks. p(k + 1) 1. If p(x) is a non-constant polynomial, then lim k p(k) (a)

More information

The graded generalized Fibonacci sequence and Binet formula

The graded generalized Fibonacci sequence and Binet formula The graded generaized Fibonacci sequence and Binet formua Won Sang Chung,, Minji Han and Jae Yoon Kim Department of Physics and Research Institute of Natura Science, Coege of Natura Science, Gyeongsang

More information

Formulas for Angular-Momentum Barrier Factors Version II

Formulas for Angular-Momentum Barrier Factors Version II BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A

More information

Product Cosines of Angles between Subspaces

Product Cosines of Angles between Subspaces Product Cosines of Anges between Subspaces Jianming Miao and Adi Ben-Israe Juy, 993 Dedicated to Professor C.R. Rao on his 75th birthday Let Abstract cos{l,m} : i cos θ i, denote the product of the cosines

More information

On the Goal Value of a Boolean Function

On the Goal Value of a Boolean Function On the Goa Vaue of a Booean Function Eric Bach Dept. of CS University of Wisconsin 1210 W. Dayton St. Madison, WI 53706 Lisa Heerstein Dept of CSE NYU Schoo of Engineering 2 Metrotech Center, 10th Foor

More information

Applications of Riordan matrix functions to Bernoulli and Euler polynomials

Applications of Riordan matrix functions to Bernoulli and Euler polynomials Illinois Wesleyan University From the SelectedWors of Tian-Xiao He 06 Applications of Riordan matrix functions to Bernoulli and Euler polynomials Tian-Xiao He Available at: https://worsbepresscom/tian_xiao_he/78/

More information

March Algebra 2 Question 1. March Algebra 2 Question 1

March Algebra 2 Question 1. March Algebra 2 Question 1 March Algebra 2 Question 1 If the statement is always true for the domain, assign that part a 3. If it is sometimes true, assign it a 2. If it is never true, assign it a 1. Your answer for this question

More information

A natural differential calculus on Lie bialgebras with dual of triangular type

A natural differential calculus on Lie bialgebras with dual of triangular type Centrum voor Wiskunde en Informatica REPORTRAPPORT A natura differentia cacuus on Lie biagebras with dua of trianguar type N. van den Hijigenberg and R. Martini Department of Anaysis, Agebra and Geometry

More information

1. Introduction Definition 1.1. Let r 1 be an integer. The r-generalized Fibonacci sequence {G n } is defined as

1. Introduction Definition 1.1. Let r 1 be an integer. The r-generalized Fibonacci sequence {G n } is defined as SOME IDENTITIES FOR r-fibonacci NUMBERS F. T. HOWARD AND CURTIS COOPER Abstract. Let r 1 be an integer. The r-generalized Fibonacci sequence {G n} is defined as 8 >< 0, if 0 n < r 1; G n = 1, if n = r

More information

On nil-mccoy rings relative to a monoid

On nil-mccoy rings relative to a monoid PURE MATHEMATICS RESEARCH ARTICLE On ni-mccoy rings reative to a monoid Vahid Aghapouramin 1 * and Mohammad Javad Nikmehr 2 Received: 24 October 2017 Accepted: 29 December 2017 First Pubished: 25 January

More information

F. T. HOWARD AND CURTIS COOPER

F. T. HOWARD AND CURTIS COOPER SOME IDENTITIES FOR r-fibonacci NUMBERS F. T. HOWARD AND CURTIS COOPER Abstract. Let r 1 be an integer. The r-generalized Fibonacci sequence {G n} is defined as 0, if 0 n < r 1; G n = 1, if n = r 1; G

More information

CHAPTER I. Rings. Definition A ring R is a set with two binary operations, addition + and

CHAPTER I. Rings. Definition A ring R is a set with two binary operations, addition + and CHAPTER I Rings 1.1 Definitions and Examples Definition 1.1.1. A ring R is a set with two binary operations, addition + and multiplication satisfying the following conditions for all a, b, c in R : (i)

More information

Research Article Numerical Range of Two Operators in Semi-Inner Product Spaces

Research Article Numerical Range of Two Operators in Semi-Inner Product Spaces Abstract and Appied Anaysis Voume 01, Artice ID 846396, 13 pages doi:10.1155/01/846396 Research Artice Numerica Range of Two Operators in Semi-Inner Product Spaces N. K. Sahu, 1 C. Nahak, 1 and S. Nanda

More information

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia

More information

Extended central factorial polynomials of the second kind

Extended central factorial polynomials of the second kind Kim et a. Advances in Difference Equations 09 09:4 https://doi.org/0.86/s366-09-963- R E S E A R C H Open Access Extended centra factoria poynomias of the second ind Taeyun Kim,DaeSanKim,Gwan-WooJang and

More information

Multiplying two generating functions

Multiplying two generating functions Generating Functions 3.3 3.6 111 Multiplying two generating functions Let A(x) = a k x k and B(x) = b k x k. k 0 k 0 What is the coefficient of x k in A(x)B(x)? Generating Functions 3.3 3.6 111 Multiplying

More information

An explicit Jordan Decomposition of Companion matrices

An explicit Jordan Decomposition of Companion matrices An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

2M2. Fourier Series Prof Bill Lionheart

2M2. Fourier Series Prof Bill Lionheart M. Fourier Series Prof Bi Lionheart 1. The Fourier series of the periodic function f(x) with period has the form f(x) = a 0 + ( a n cos πnx + b n sin πnx ). Here the rea numbers a n, b n are caed the Fourier

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

A NOTE ON TRIGONOMETRIC MATRICES

A NOTE ON TRIGONOMETRIC MATRICES A NOTE ON TRIGONOMETRIC MATRICES garret j. etgen Introduction. Let Q(x) he an nxn symmetric matrix of continuous functions on X: 0^x

More information

Fitting affine and orthogonal transformations between two sets of points

Fitting affine and orthogonal transformations between two sets of points Mathematica Communications 9(2004), 27-34 27 Fitting affine and orthogona transformations between two sets of points Hemuth Späth Abstract. Let two point sets P and Q be given in R n. We determine a transation

More information

Odds find E n d s ENTRIES. and A L T O SOLO. iiik places for many species of liik JMtHSIDKVrS r ITCH I

Odds find E n d s ENTRIES. and A L T O SOLO. iiik places for many species of liik JMtHSIDKVrS r ITCH I LG Bx L G L L FY-F Y LLL G Y 593 JBL F GBB 3 P P P P P B P P P g J g -g ( J z x 000 P gg - F g L g g L YB J 933 g 95 Y LL g g g g- PY L FG LLL F LLX F q F PP06 BG P P P F g 000 g F L 93 PL - L g g G gg

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

1 The distributive law

1 The distributive law THINGS TO KNOW BEFORE GOING INTO DISCRETE MATHEMATICS The distributive law The distributive law is this: a(b + c) = ab + bc This can be generalized to any number of terms between parenthesis; for instance:

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

Legendre Polynomials - Lecture 8

Legendre Polynomials - Lecture 8 Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha

More information

On the Computation of (2-2) Three- Center Slater-Type Orbital Integrals of 1/r 12 Using Fourier-Transform-Based Analytical Formulas

On the Computation of (2-2) Three- Center Slater-Type Orbital Integrals of 1/r 12 Using Fourier-Transform-Based Analytical Formulas On the Computation of (2-2) Three- Center Sater-Type Orbita Integras of /r 2 Using Fourier-Transform-Based Anaytica Formuas DANKO ANTOLOVIC, HARRIS J. SILVERSTONE 2 Pervasive Technoogy Labs, Indiana University,

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

b n n=1 a n cos nx (3) n=1

b n n=1 a n cos nx (3) n=1 Fourier Anaysis The Fourier series First some terminoogy: a function f(x) is periodic if f(x ) = f(x) for a x for some, if is the smaest such number, it is caed the period of f(x). It is even if f( x)

More information

c 2000 Society for Industrial and Applied Mathematics

c 2000 Society for Industrial and Applied Mathematics SIAM J. SCI. COMPUT. Vo. 2, No. 5, pp. 909 926 c 2000 Society for Industria and Appied Mathematics A DEFLATED VERSION OF THE CONJUGATE GRADIENT ALGORITHM Y. SAAD, M. YEUNG, J. ERHEL, AND F. GUYOMARC H

More information

221B Lecture Notes Notes on Spherical Bessel Functions

221B Lecture Notes Notes on Spherical Bessel Functions Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,

More information

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain CORRECTIONS TO CLASSICAL PROCEDURES FOR ESTIMATING THURSTONE S CASE V MODEL FOR RANKING DATA Aberto Maydeu Oivares Instituto de Empresa Marketing Dept. C/Maria de Moina -5 28006 Madrid Spain Aberto.Maydeu@ie.edu

More information

Linear Network Coding for Multiple Groupcast Sessions: An Interference Alignment Approach

Linear Network Coding for Multiple Groupcast Sessions: An Interference Alignment Approach Linear Network Coding for Mutipe Groupcast Sessions: An Interference Aignment Approach Abhik Kumar Das, Siddhartha Banerjee and Sriram Vishwanath Dept. of ECE, The University of Texas at Austin, TX, USA

More information

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix

More information

Relations, Functions, Binary Relations (Chapter 1, Sections 1.2, 1.3)

Relations, Functions, Binary Relations (Chapter 1, Sections 1.2, 1.3) Relations, Functions, Binary Relations (Chapter 1, Sections 1.2, 1.3) CmSc 365 Theory of Computation 1. Relations Definition: Let A and B be two sets. A relation R from A to B is any set of ordered pairs

More information

Efficient Algorithms for Pairing-Based Cryptosystems

Efficient Algorithms for Pairing-Based Cryptosystems CS548 Advanced Information Security Efficient Agorithms for Pairing-Based Cryptosystems Pauo S. L. M. Barreto, HaeY. Kim, Ben Lynn, and Michae Scott Proceedings of Crypto, 2002 2010. 04. 22. Kanghoon Lee,

More information

Research Article Building Infinitely Many Solutions for Some Model of Sublinear Multipoint Boundary Value Problems

Research Article Building Infinitely Many Solutions for Some Model of Sublinear Multipoint Boundary Value Problems Abstract and Appied Anaysis Voume 2015, Artice ID 732761, 4 pages http://dx.doi.org/10.1155/2015/732761 Research Artice Buiding Infinitey Many Soutions for Some Mode of Subinear Mutipoint Boundary Vaue

More information

CONGRUENCES. 1. History

CONGRUENCES. 1. History CONGRUENCES HAO BILLY LEE Abstract. These are notes I created for a seminar tak, foowing the papers of On the -adic Representations and Congruences for Coefficients of Moduar Forms by Swinnerton-Dyer and

More information

Mixed Volume Computation, A Revisit

Mixed Volume Computation, A Revisit Mixed Voume Computation, A Revisit Tsung-Lin Lee, Tien-Yien Li October 31, 2007 Abstract The superiority of the dynamic enumeration of a mixed ces suggested by T Mizutani et a for the mixed voume computation

More information

Available online at ScienceDirect. Procedia Computer Science 96 (2016 )

Available online at  ScienceDirect. Procedia Computer Science 96 (2016 ) Avaiabe onine at www.sciencedirect.com ScienceDirect Procedia Computer Science 96 (206 92 99 20th Internationa Conference on Knowedge Based and Inteigent Information and Engineering Systems Connected categorica

More information

are left and right inverses of b, respectively, then: (b b 1 and b 1 = b 1 b 1 id T = b 1 b) b 1 so they are the same! r ) = (b 1 r = id S b 1 r = b 1

are left and right inverses of b, respectively, then: (b b 1 and b 1 = b 1 b 1 id T = b 1 b) b 1 so they are the same! r ) = (b 1 r = id S b 1 r = b 1 Lecture 1. The Category of Sets PCMI Summer 2015 Undergraduate Lectures on Fag Varieties Lecture 1. Some basic set theory, a moment of categorica zen, and some facts about the permutation groups on n etters.

More information

Two-Stage Least Squares as Minimum Distance

Two-Stage Least Squares as Minimum Distance Two-Stage Least Squares as Minimum Distance Frank Windmeijer Discussion Paper 17 / 683 7 June 2017 Department of Economics University of Bristo Priory Road Compex Bristo BS8 1TU United Kingdom Two-Stage

More information

HIRZEBRUCH χ y GENERA OF THE HILBERT SCHEMES OF SURFACES BY LOCALIZATION FORMULA

HIRZEBRUCH χ y GENERA OF THE HILBERT SCHEMES OF SURFACES BY LOCALIZATION FORMULA HIRZEBRUCH χ y GENERA OF THE HILBERT SCHEMES OF SURFACES BY LOCALIZATION FORMULA KEFENG LIU, CATHERINE YAN, AND JIAN ZHOU Abstract. We use the Atiyah-Bott-Berine-Vergne ocaization formua to cacuate the

More information

NIKOS FRANTZIKINAKIS. N n N where (Φ N) N N is any Følner sequence

NIKOS FRANTZIKINAKIS. N n N where (Φ N) N N is any Følner sequence SOME OPE PROBLEMS O MULTIPLE ERGODIC AVERAGES IKOS FRATZIKIAKIS. Probems reated to poynomia sequences In this section we give a ist of probems reated to the study of mutipe ergodic averages invoving iterates

More information

BALANCING REGULAR MATRIX PENCILS

BALANCING REGULAR MATRIX PENCILS BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity

More information

11.10a Taylor and Maclaurin Series

11.10a Taylor and Maclaurin Series 11.10a 1 11.10a Taylor and Maclaurin Series Let y = f(x) be a differentiable function at x = a. In first semester calculus we saw that (1) f(x) f(a)+f (a)(x a), for all x near a The right-hand side of

More information

Radiation Fields. Lecture 12

Radiation Fields. Lecture 12 Radiation Fieds Lecture 12 1 Mutipoe expansion Separate Maxwe s equations into two sets of equations, each set separatey invoving either the eectric or the magnetic fied. After remova of the time dependence

More information

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l Investigation on spectrum of the adjacency matrix and Lapacian matrix of graph G SHUHUA YIN Computer Science and Information Technoogy Coege Zhejiang Wani University Ningbo 3500 PEOPLE S REPUBLIC OF CHINA

More information

FOURIER SERIES ON ANY INTERVAL

FOURIER SERIES ON ANY INTERVAL FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series

More information

SYMMETRY AND SPECIALIZABILITY IN THE CONTINUED FRACTION EXPANSIONS OF SOME INFINITE PRODUCTS

SYMMETRY AND SPECIALIZABILITY IN THE CONTINUED FRACTION EXPANSIONS OF SOME INFINITE PRODUCTS SYMMETRY AND SPECIALIZABILITY IN THE CONTINUED FRACTION EXPANSIONS OF SOME INFINITE PRODUCTS J MC LAUGHLIN Abstract Let fx Z[x] Set f 0x = x and for n 1 define f nx = ff n 1x We describe several infinite

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Statistical Learning Theory: A Primer

Statistical Learning Theory: A Primer Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

School of Electrical Engineering, University of Bath, Claverton Down, Bath BA2 7AY

School of Electrical Engineering, University of Bath, Claverton Down, Bath BA2 7AY The ogic of Booean matrices C. R. Edwards Schoo of Eectrica Engineering, Universit of Bath, Caverton Down, Bath BA2 7AY A Booean matrix agebra is described which enabes man ogica functions to be manipuated

More information

Haar Decomposition and Reconstruction Algorithms

Haar Decomposition and Reconstruction Algorithms Jim Lambers MAT 773 Fa Semester 018-19 Lecture 15 and 16 Notes These notes correspond to Sections 4.3 and 4.4 in the text. Haar Decomposition and Reconstruction Agorithms Decomposition Suppose we approximate

More information

Laplace - Fibonacci transform by the solution of second order generalized difference equation

Laplace - Fibonacci transform by the solution of second order generalized difference equation Nonauton. Dyn. Syst. 017; 4: 30 Research Artice Open Access Sandra Pineas*, G.B.A Xavier, S.U. Vasantha Kumar, and M. Meganathan Lapace - Fibonacci transform by the soution of second order generaized difference

More information

Linear Algebra. Linear Algebra. Chih-Wei Yi. Dept. of Computer Science National Chiao Tung University. November 12, 2008

Linear Algebra. Linear Algebra. Chih-Wei Yi. Dept. of Computer Science National Chiao Tung University. November 12, 2008 Linear Algebra Chih-Wei Yi Dept. of Computer Science National Chiao Tung University November, 008 Section De nition and Examples Section De nition and Examples Section De nition and Examples De nition

More information

Unsolved Problems In Computational Science: I

Unsolved Problems In Computational Science: I Unsolved Problems In Computational Science: I Shanzhen Gao, Keh-Hsun Chen Department of Computer Science, College of Computing and Informatics University of North Carolina at Charlotte, Charlotte, NC 28223,

More information

Generating Functions and the Fibonacci Sequence

Generating Functions and the Fibonacci Sequence Department of Mathematics Nebraska Wesleyan University June 14, 01 Fibonacci Sequence Fun Fact: November 3rd is Fibonacci Day! (1, 1,, 3) Definition The Fibonacci sequence is defined by the recurrence

More information

School of Electrical Engineering, University of Bath, Claverton Down, Bath BA2 7AY

School of Electrical Engineering, University of Bath, Claverton Down, Bath BA2 7AY The ogic of Booean matrices C. R. Edwards Schoo of Eectrica Engineering, Universit of Bath, Caverton Down, Bath BA2 7AY A Booean matrix agebra is described which enabes man ogica functions to be manipuated

More information

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5]. PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform

More information

QUADRATIC FORMS AND FOUR PARTITION FUNCTIONS MODULO 3

QUADRATIC FORMS AND FOUR PARTITION FUNCTIONS MODULO 3 QUADRATIC FORMS AND FOUR PARTITION FUNCTIONS MODULO 3 JEREMY LOVEJOY AND ROBERT OSBURN Abstract. Recenty, Andrews, Hirschhorn Seers have proven congruences moduo 3 for four types of partitions using eementary

More information

MS 2001: Test 1 B Solutions

MS 2001: Test 1 B Solutions MS 2001: Test 1 B Solutions Name: Student Number: Answer all questions. Marks may be lost if necessary work is not clearly shown. Remarks by me in italics and would not be required in a test - J.P. Question

More information

Construction of Supersaturated Design with Large Number of Factors by the Complementary Design Method

Construction of Supersaturated Design with Large Number of Factors by the Complementary Design Method Acta Mathematicae Appicatae Sinica, Engish Series Vo. 29, No. 2 (2013) 253 262 DOI: 10.1007/s10255-013-0214-6 http://www.appmath.com.cn & www.springerlink.com Acta Mathema cae Appicatae Sinica, Engish

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

be a formal power series with h 0 = 1, v(h^) = \ and v(h n ) > 1 for n > 1. Define a n by

be a formal power series with h 0 = 1, v(h^) = \ and v(h n ) > 1 for n > 1. Define a n by THE POWER OF 2 DIVIDING THE COEFFICIENTS OF CERTAIN POWER SERIES F. T. Howard Wake Forest University, Box 7388 Reynolda Station, Winston-Salem, NC 27109 (Submitted August 1999-Final Revision January 2000)

More information

Inverses and Elementary Matrices

Inverses and Elementary Matrices Inverses and Elementary Matrices 1-12-2013 Matrix inversion gives a method for solving some systems of equations Suppose a 11 x 1 +a 12 x 2 + +a 1n x n = b 1 a 21 x 1 +a 22 x 2 + +a 2n x n = b 2 a n1 x

More information

MAS 315 Waves 1 of 8 Answers to Examples Sheet 1. To solve the three problems, we use the methods of 1.3 (with necessary changes in notation).

MAS 315 Waves 1 of 8 Answers to Examples Sheet 1. To solve the three problems, we use the methods of 1.3 (with necessary changes in notation). MAS 35 Waves of 8 Answers to Exampes Sheet. From.) and.5), the genera soution of φ xx = c φ yy is φ = fx cy) + gx + cy). Put c = : the genera soution of φ xx = φ yy is therefore φ = fx y) + gx + y) ) To

More information

Summation of p-adic Functional Series in Integer Points

Summation of p-adic Functional Series in Integer Points Fiomat 31:5 (2017), 1339 1347 DOI 10.2298/FIL1705339D Pubished by Facuty of Sciences and Mathematics, University of Niš, Serbia Avaiabe at: http://www.pmf.ni.ac.rs/fiomat Summation of p-adic Functiona

More information

Rapid and Stable Determination of Rotation Matrices between Spherical Harmonics by Direct Recursion

Rapid and Stable Determination of Rotation Matrices between Spherical Harmonics by Direct Recursion Chemistry Pubications Chemistry 11-1999 Rapid and Stabe Determination of Rotation Matrices between Spherica Harmonics by Direct Recursion Cheo Ho Choi Iowa State University Joseph Ivanic Iowa State University

More information

EXERCISES. a b = a + b l aq b = ab - (a + b) + 2. a b = a + b + 1 n0i) = oii + ii + fi. A. Examples of Rings. C. Ring of 2 x 2 Matrices

EXERCISES. a b = a + b l aq b = ab - (a + b) + 2. a b = a + b + 1 n0i) = oii + ii + fi. A. Examples of Rings. C. Ring of 2 x 2 Matrices / rings definitions and elementary properties 171 EXERCISES A. Examples of Rings In each of the following, a set A with operations of addition and multiplication is given. Prove that A satisfies all the

More information

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS

More information

Induction, sequences, limits and continuity

Induction, sequences, limits and continuity Induction, sequences, limits and continuity Material covered: eclass notes on induction, Chapter 11, Section 1 and Chapter 2, Sections 2.2-2.5 Induction Principle of mathematical induction: Let P(n) be

More information