ELA

Size: px
Start display at page:

Download "ELA"

Transcription

1 Volume 18, pp , August GENERALIZED PASCAL TRIANGLES AND TOEPLITZ MATRICES A R MOGHADDAMFAR AND S M H POOYA Abstract The purpose of this article is to study determinants of matrices which are known as generalized Pascal triangles (see R Bacher Determinants of matrices related to the Pascal triangle J Théor Nombres Bordeaux, 14:19 41, 2002 This article presents a factorization by expressing such a matrix as a product of a unipotent lower triangular matrix, a Toeplitz matrix, and a unipotent upper triangular matrix The determinant of a generalized Pascal matrix equals thus the determinant of a Toeplitz matrix This equality allows for the evaluation of a few determinants of generalized Pascal matrices associated with certain sequences In particular, families of quasi-pascal matrices are obtained whose leading principal minors generate any arbitrary linear subsequences (F nr+s n 1 or (L nr+s n 1 of the Fibonacci or Lucas sequence New matrices are constructed whose entries are given by certain linear non-homogeneous recurrence relations, and the leading principal minors of which form the Fibonacci sequence Key words Determinant, Matrix factorization, Generalized Pascal triangle, Generalized symmetric (skymmetric Pascal triangle, Toeplitz matrix, Recursive relation, Fibonacci (Lucas, Catalan sequence, Golden ratio AMS subject classifications 15A15, 11C20 1 Introduction Let P ( be the infinite symmetric matrix with entries P i,j = ( i+j i for i, j 0 The matrix P ( is hence the famous Pascal triangle yielding the binomial coefficients The entries of P ( satisfy the recurrence relation P i,j = P i 1,j + P i,j 1 Indeed, this matrix has the following form: P ( = Received by the editors January 19, 2009 Accepted for publication August 1, 2009 Handling Editor: Michael J Tsatsomeros Department of Mathematics, Faculty of Science, K N Toosi University of Technology, P O Box , Tehran, Iran (moghadam@kntuacir School of Mathematics, Institute for Studies in Theoretical Physics and Mathematics (IPM, PO Box , Tehran, Iran (moghadam@mailipmir This research was in part supported by a grant from IPM (No Department of Mathematics, Faculty of Science, K N Toosi University of Technology, P O Box , Tehran, Iran 564

2 Volume 18, pp , August One can easily verify that (see [2, 6]: Generalized Pascal Triangles and Toeplitz Matrices 565 (11 P ( =L( L( t, where L( is the infinite unipotent lower triangular matrix (12 L( =(L i,j i,j 0 = with entries L i,j = ( i j WedenotebyL(n the finite submatrix of L( withentries L i,j,0 i, j n 1 To introduce our result, we first present some notation and definitions We recall that a matrix T ( =(t i,j i,j 0 is said to be Toeplitz if t i,j = t k,l whenever i j = k l Let α =(α i i 0 and β =(β i i 0 be two sequences with α 0 = β 0 We shall denote by T α,β ( =(t i,j i,j 0 the Toeplitz matrix with t i,0 = α i and t 0,j = β j We also denote by T α,β (n the submatrix of T α,β ( consisting of the entries in its first n rows and columns We come now back to the (11 In fact, one can rewrite it as follows: P ( =L( I L( t, where matrix I (identity matrix is a particular case of a Toeplitz matrix In [1], Bacher considers determinants of matrices generalizing the Pascal triangle P ( He introduces generalized Pascal triangles as follows Let α =(α i i 0 and β = (β i i 0 be two sequences starting with a common first term α 0 = β 0 = γ Then, the generalized Pascal triangle associated with α and β, is the infinite matrix P α,β ( =(P i,j i,j 0 with entries P i,0 = α i, P 0,j = β j (i, j 0 and P i,j = P i 1,j + P i,j 1, for i, j 1 We denote by P α,β (n the finite submatrix of P α,β ( withentriesp i,j,0 i, j n 1 An explicit formula for entry P i,j of P α,β (n is also given by the following formula (see [1]: ( i + j ( i ( i + j s ( j ( i + j t P i,j = γ + (α s α s 1 + (β t β t 1 j j i s=1 t=1

3 Volume 18, pp , August ARMoghaddamfarandSMHPooya For an arbitrary sequence α =(α i i 0, we define the sequences ˆα =(ˆα i i 0 and ˇα =(ˇα i i 0 as follows: (13 ˆα i = i ( i ( 1 i+k α k and ˇα i = k k=0 i k=0 ( i α k k With these definitions we can now state our main result Indeed, the purpose of this article is to obtain a factorization of the generalized Pascal triangle P α,β (n associated with the arbitrary sequences α and β, as a product of a unipotent lower triangular matrix L(n, a Toeplitz matrix T ˆα, ˆβ(n and a unipotent upper triangular matrix U(n (see Theorem 31, that is Similarly, we show that P α,β (n =L(n T ˆα, ˆβ(n U(n T α,β (n =L(n 1 Pˇα, ˇβ(n U(n 1 In fact, we obtain a connection between generalized Pascal triangles and Toeplitz matrices In view of these factorizations, we can easily see that det(p α,β (n = det(t ˆα, ˆβ(n Finally, we present several applications of Theorem 31 to some other determinant evaluations We conclude the introduction with notation and terminology to be used throughout the article For convenience, we will let 1 denote the complex number i C By x we denote the integer part of x, ie, the greatest integer that is less than or equal to x Wealsodenoteby x the smallest integer greater than or equal to x Given a matrix A, wedenotebyr i (A andc j (A therowi and the column j of A, respectively We use the notation A t for the transpose of A In general, an n n matrix of the following form: [ ] [ A B ( A B resp C P α,β (n k C T α,β (n k ] where A, B and C are arbitrary matrices of order k k, k (n k and(n k k, respectively, is called a quasi-pascal (resp quasi-toeplitz matrix Throughout this article we assume that:

4 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 567 F =(F i i 0 =(0, 1, 1, 2, 3, 5, 8,,F i, F =(F i i 1 =(1, 1, 2, 3, 5, 8,,F i, L =(L i i 0 =(2, 1, 3, 4, 7, 11, 18,,L i, C =(C i i 0 =(1, 1, 2, 5, 14, 42, 132,,C i, I =(I i i 0 =(0!, 1!, 2!, 3!, 4!, 5!, 6!,,I i = i!, I =(I i i 1 =(1!, 2!, 3!, 4!, 5!, 6!,,I i = i!, (Fibonacci numbers, (Fibonacci numbers 0, (Lucas numbers, (Catalan numbers, The rest of this article is organized as follows In Section 2, we derive some preparatory results In Section 3, we prove the main result (Theorem 31 Section 4 deals with applications of Theorem 31 In Section 5, we present two matrices whose entries are recursively defined, and we show that the leading principal minor sequence of these matrices is the sequence F =(F i i 1 2 Preliminary Results As we mentioned in the Introduction, if α =(α i i 0 is an arbitrary sequence, then we define the sequences ˆα =(ˆα i i 0 and ˇα =(ˇα i i 0 as in (13 For certain sequences α the associated sequences ˆα and ˇα seem also to be of interest since they have appeared elsewhere In Tables 1 and 2, we have presented some sequences α and the associated sequences ˆα and ˇα Table 1 Some sequences α and associated sequences ˇα α Reference ˇα (0, 1, 1, 2, 3, 5, 8, A in [8] F (1, 0, 1, 1, 2, 3, 5, 8, A in [8] F (2, 1, 3, 4, 7, 11, 18, A in [8] L (1, 0, 1, 1, 3, 6, 15, A in [8] C (1, 0, 1, 2, 9, 44, 265, A in [8] I (1, 1, 3, 11, 53, 309, 2119, A in [8] I Table 2 Some sequences α and associated sequences ˆα α Reference ˆα (0, 1, 3, 8, 21, 55, 144, A in [8] F (1, 2, 5, 13, 34, 89, 233, A in [8] F (2, 3, 7, 18, 47, 123, 322, A in [8] L (1, 2, 5, 15, 51, 188, 731, A in [8] C (1, 2, 5, 16, 65, 326, 1957, A in [8] I (1, 3, 11, 49, 261, 1631, A in [8] I

5 Volume 18, pp , August ARMoghaddamfarandSMHPooya As another nice example, consider α =R i (L( where L is introduced as in (12 Then, we have ˇα =R i (P ( Lemma 21 Let α be an arbitrary sequence Then, we have ˆˇα = ˇˆα = α Proof Suppose α =(α i i 0 and ˆˇα =(ˆˇα i i 0 Then, we have ˆˇα i = i k=0 ( 1i+k( i kˇαk = i k=0 ( 1i+k( i k ( k k s=0 s αs = i k k=0 s=0 ( 1i+k( i k k( s αs = ( 1 i i l=0 α i l h=l ( 1h( i h h( l = ( 1 i i l=0 α i l h=l ( 1h( ( i i l l h l = ( 1 i i l=0 α i i l( l h=l ( 1h( i l h l But, if l<i,thenwehave i h=l ( 1h( i l h l = 0 Therefore, we obtain and hence ˆˇα = α ˆˇα i =( 1 i i ( i ( i i l α l ( 1 h l h l l=i h=l The proof of second part is similar to the previous case = α i, Lemma 22 Let i, j be positive integers Then, we have i j ( ( { i k + j 0 if i j, ( 1 k = k + j j 1 if i = j k=0 The proof follows from the easy identity ( ( i k + j = k + j j ( i j ( i j k The following Lemma is a special case of a general result due to Krattenthaler (see Theorem 1 in [9] Lemma 23 Let α =(α i i 0 and β =(β j j 0 be two geometric sequences with α i = ρ i and β j = σ j Then, we have det(p α,β (n = (ρ + σ ρσ

6 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices Main Result Now, we are in the position to state and prove the main result of this article Theorem 31 Let α =(α i i 0 and β =(β i i 0 be two sequences starting with a common first term α 0 = β 0 = γ Then, we have (31 P α,β (n =L(n T ˆα, ˆβ(n U(n, and (32 T α,β (n =L(n 1 Pˇα, ˇβ(n U(n 1, where L(n =(L i,j 0 i,j<n is a lower triangular matrix with L i,j = { 0 if i<j ( i j if i j, and U(n =L(n t In particular, we have det(p α,β (n = det(t ˆα, ˆβ(n Proof First, we claim that P α,β (n =L(n Q(n, where L(n =(L i,j 0 i,j<n is a lower triangular matrix with L i,j = { 0 if i<j ( i j if i j, and Q(n =(Q i,j 0 i,j<n with Q i,0 =ˆα i, Q 0,i = β i and (33 Q i,j = Q i 1,j 1 + Q i,j 1, 1 i, j < n For instance, when n = 4 the matrices L(4 and Q(4 are given by: and Q(4 = L(4 = , γ β 1 β 2 β 3 γ + α 1 α 1 β 1 + α 1 β 1 + β 2 + α 2 γ 2α 1 + α 2 α 1 + α 2 α 2 β 1 + α 1 + α 2 γ +3α 1 3α 2 + α 3 α 1 2α 2 + α 3 α 2 + α 3 α 3

7 Volume 18, pp , August ARMoghaddamfarandSMHPooya Note that the entries of L(n satisfying in the following recurrence (34 L i,j = L i 1,j 1 + L i 1,j, 1 i, j < n For the proof of the claimed factorization we compute the (i, j-th entry of L(n Q(n, that is n (L(n Q(n i,j = L i,k Q k,j k=1 In fact, it suffices to show that R 0 (L(n Q(n = R 0 (P α,β (n, C 0 (L(n Q(n = C 0 (P α,β (n and (35 (L(n Q(n i,j =(L(n Q(n i,j 1 +(L(n Q(n i 1,j, for 1 i, j < n First, suppose that i = 0 Then, we obtain (L(n Q(n 0,j = L 0,k Q k,j = L 0,0 Q 0,j = β j, k=0 and so R 0 (L(n Q(n = R 0 (P α,β (n = (β 0, β 1,, β Next, suppose that i 1andj = 0 In this case, we have (L(n Q(n i,0 = = i k=0 k=0 L i,kq k,0 { (i k = i j=0 β j k l=0 ( 1l+k( } k l βl { i j t=0 ( 1t( ( i t+j } t+j j = β i, (by Lemma 22 and so C 0 (L(n Q(n = C 0 (P α,β (n = (α 0, α 1,, α t

8 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 571 Finally, we must establish (35 Therefore, we assume that 1 i, j < n Inthis case, we have (L(n Q(n i,j = k=0 L i,kq k,j which is (35 Next, we claim that = L i,0 Q 0,j + k=1 L i,kq k,j = L i,0 Q 0,j + k=1 L i,k(q k 1,j 1 + Q k,j 1 (by (33 = L i,0 Q 0,j + k=1 L i,kq k 1,j 1 + k=1 L i,kq k,j 1 = L i,0 Q 0,j + k=1 (L i 1,k 1 + L i 1,k Q k 1,j 1 + k=0 L i,kq k,j 1 L i,0 Q 0,j 1 (by (34 = L i,0 Q 0,j + k=1 L i 1,k 1Q k 1,j 1 + k=1 L i 1,k(Q k,j Q k,j 1 +(L(n Q(n i,j 1 L i,0 Q 0,j 1 (by (33 = L i,0 Q 0,j + n 2 k=0 L i 1,kQ k,j 1 + k=1 L i 1,kQ k,j k=0 L i 1,kQ k,j 1 + L i 1,0 Q 0,j 1 +(L(n Q(n i,j 1 L i,0 Q 0,j 1 = L i,0 Q 0,j + k=0 L i 1,kQ k,j 1 + k=0 L i 1,kQ k,j L i 1,0 Q 0,j k=0 L i 1,kQ k,j 1 + L i 1,0 Q 0,j 1 +(L(n Q(n i,j 1 L i,0 Q 0,j 1 (note that L i 1, =0 = (L i,0 L i 1,0 Q 0,j +(L(n Q(n i 1,j +(L i 1,0 L i,0 Q 0,j 1 +(L(n Q(n i,j 1 = (L(n Q(n i 1,j +(L(n Q(n i,j 1, (note that L i,0 = L i 1,0 =1 Q(n =T ˆα, ˆβ(n U(n, where U(n =L(n t and T α,β (n =(T i,j 0 i,j<n with T i,0 =ˆα i, T 0,j = ˆβ j,and Note that, we have T i,j = T i 1,j 1, 1 i, j < n (36 U i,j = U i 1,j 1 + U i,j 1, 1 i, j < n

9 Volume 18, pp , August ARMoghaddamfarandSMHPooya For instance, when n = 4 the matrices T ˆα, ˆβ(4 and U(4 are given by: 2 γ γ + β 1 γ 2β 1 + β 2 γ +3β 1 3β 2 + β 3 Tˆα, ˆβ(4 γ + α 1 γ γ + β 1 γ 2β 1 + β 2 = 6 4 γ 2α 1 + α 2 γ + α 1 γ γ + β 1 γ +3α 1 3α 2 + α 3 γ 2α 1 + α 2 γ + α 1 γ and U(4 = As before, the proof of the claim requires some calculations If we have i =0, then (T ˆα, ˆβ(n U(n 0,j = = j k=0 = j i=0 β i k=0 T 0,kU k,j {( k (j l=0 ( 1l+k( k l } βl k { j i t=0 ( 1t( ( j t+i } t+i i = β j, (by Lemma 22 which implies that R 0 (T ˆα, ˆβ(n U(n = R 0 (Q(n = (β 0, β 1,, β If j =0, then we obtain (T ˆα, ˆβ(n U(n i,0 = T i,k U k,1 = T i,0 U 0,0 = ˆα i, k=0 and so C 0 (T ˆα, ˆβ(n U(n = C 0 (Q(n = (ˆα 0, ˆα 1,, ˆα t Finally, we assume that 1 i, j < n 1 and establish (33 Indeed, by calculations we observe that (T ˆα, ˆβ(n U(n i,j = k=0 T i,ku k,j = T i,0 U 0,j + k=1 T i,ku k,j = T i,0 U 0,j + k=1 T i,k(u k 1,j 1 + U k,j 1 (by (36 = T i,0 U 0,j + k=1 T i,ku k 1,j 1 + k=1 T i,ku k,j 1 = T i,0 U 0,j + k=1 T i 1,k 1U k 1,j 1 + k=0 T i,ku k,j 1 T i,0 U 0,j 1 (note that, T i,k = T i 1,k 1 = T i,0 (U 0,j U 0,j 1 + k=0 T i 1,kU k,j 1 +(T ˆα, ˆβ(n U(n i,j 1 (note that U,j 1 =0 = (T ˆα, ˆβ(n U(n i 1,j 1 +(T ˆα, ˆβ(n U(n i,j 1,

10 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 573 which is (33 The proof of (31 is now complete To prove of (32, we observe that L(n 1 Pˇα, ˇβ(n U(n 1 = L(n 1 L(n Tˆˇα,ˆˇβ(n U(n U(n 1 (by (31 = Tˆˇα,ˆˇβ(n = T α,β (n (by Lemma 21 The proof is complete 4 Some Applications Let A =(a i,j i,j 0 be an infinite matrix and let D n be the nth leading principal minor of A consisting of the entries in its first n rows and columns We will mainly be interested in the sequence of leading principal minors (D 1,D 2,D 3,, especially, in the case that it forms a Fibonacci or Lucas (subsequence 41 Generalized Pascal Triangle Associated With an Arithmetic or Geometric Sequence It is of interest to evaluate the determinant of generalized Pascal triangle P α,β (n, where one of the sequences α or β is an arithmetic or geometric sequence Corollary 41 Let α =(α i i 0 be an arithmetic sequence with α i = a + id, and let β =(β i i 0 be an arbitrary sequence with β 0 = a WesetD n =det(p α,β (n Then, we have with D 0 =1 D n = ( d k ˆβk D n k 1, k=0 Proof By Theorem 31, we deduce that D n =det(t ˆα, ˆβ(n, where ˆα =(ˆα i i 0 with ˆα i = i k=0 ( 1i+k( i k αk ( β0 + kd = i k=0 ( 1i+k( i k = β i 0 k=0 ( 1i+k( i k + d i k=0 ( 1i+k( i k k β 0 if i =0, = d if i =1, 0 if i>1 Now, expanding through the first row of T ˆα, ˆβ(n, we obtain the result

11 Volume 18, pp , August ARMoghaddamfarandSMHPooya Corollary 42 Let α =(α i i 0 be an arithmetic sequence with α i = a + id, and let β =(β i i 0 be an alternating sequence with β i =( 1 i a Then, we have det(p α,β (n = a(2d + a Proof LetD n =det(p α,β (n By Corollary 41, we have D n = ( d k ˆβk D n k 1 (n 1, k=0 with D 0 = 1 An easy calculation shows that ˆβ k = k l=0 ( 1k+l( k l βl Hence, we have = k l=0 ( 1k+l( k l ( 1 l a = a( 1 k k ( k l=0 l = a( 2 k D n = a (2d k D n k 1 k=0 Now, replacing n by n 1, we obtain n 2 D = a (2d k 1 D n k 2, k=0 and by calculation it follows that ( D n 2dD = a k=0 (2dk D n k 1 n 2 k=0 (2dk 1 D n k 2 ( = a k=0 (2dk D n k 1 k=1 (2dk D n k 1 = ad, or equivalently But, this implies that D n = a(2d + a D n =(2d + ad Corollary 43 Let α =(α i i 0 be an arithmetic sequence with α i = id, and let β =(β i i 0 be the square sequence, ie, β i = i 2 If D n =det(p α,β (n, thenwe have (41 D n = dd n 2 +2d 2 D n 3

12 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 575 Proof By Corollary 41, we get D n = ( d k ˆβk D n k 1 k=0 We claim that ˆβ 0 =0, ˆβ 1 =1, ˆβ 2 =2and ˆβ k =0fork 3 Now, it is obvious that our claim implies the validity of (41 Clearly ˆβ 0 =0, ˆβ 1 =1and ˆβ 2 = 2 Now, we assume that k 3 In this case we have k ( k k ( k ˆβ k = ( 1 k+j β j = ( 1 k+j j 2 j j j=0 We define the functions f and g as follows: k ( k f(x =(1 x k = ( x j j and g(x = xf (x j=0 Now, an easy calculation shows that and putting x =1,weget g (x =(1 x k 2 ( = k j=0 j=0 j=0 k j=0 ( k j 2 ( 1 j 1 =0 j Now, by multiplying both sides by ( 1 k+1,weobtain k ( k j 2 ( 1 k+j =0, j as claimed ( k j 2 ( x j 1, j Corollary 44 Let α =(α i i 0 and β =(β j j 0 be two geometric sequences with α i = ρ i and β j = σ j Then, we have det(t α,β (n = (1 ρσ Proof By Theorem 31, we have det(t α,β (n = det(pˇα, ˇβ(n On the other hand, straightforward computations show that ˇα =(ˇα i i 0 with ˇα i =(1+ρ i and similarly ˇβ =(ˇβ j j 0 with ˇβ j =(1+σ j By applying Lemma 23, we conclude the assertion

13 Volume 18, pp , August ARMoghaddamfarandSMHPooya 42 Certain Generalized Pascal Triangles Proposition 45 Let a, b, c C and let n be a positive integer Let α =(α i i 0 and β =(β j j 0 be two sequences with α i =(2 i 1a + c and β j =(2 j 1b + c Then, we have det(p α,β (n = 1 n c if a = b = c, [ ] c + a(n 1 (c a if a = b c, b b a (c an + a b a (c bn if a b Proof By Theorem 31, we have det(p α,β (n = det(t ˆα, ˆβ(n A straightforward computation shows that ˆα =(c,a,a,a, and ˆβ =(c,b,b,b, Therefore, in the notation of [7], we have T ˆα, ˆβ(n =M n (b, a, c, and since M n (a, b, c = M n (b, a, c t we have det(t ˆα, ˆβ(n = det(m n (a, b, c Now, the proof follows the lines in the proof of Theorem 2 in [7] Proposition 46 Let a, b, c C and let n be a positive integer Let α =(α i i 0 and β =(β j j 0 be two sequences with α i =2 i 1 (ia +2c and β j =2 j 1 (jb +2c Then, we have det(p α,β (n = ( 1 n+1 (a + b n 2[ c(a + b+(n 1ab ] Proof Again from Theorem 31, we have det(p α,β (n = det(t ˆα, ˆβ(n, where ˆα and ˆβ are two arithmetic sequences as ˆα =(ˆα i i 0 =(c, c + a, c +2a,,c+ ia, and ˆβ =(ˆβ j j 0 =(c, c + b, c +2b,,c+ jb, Now we compute the determinant of T ˆα, ˆβ(n elementary column operations: To do this, we apply the following C j C j C j 1, j = n 1,n 2,,2;

14 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 577 and we obtain the following quasi-toeplitz matrix: c b b b c + a c +2a T λ,µ (n 1 c +(a, where λ =( a, a, a,andµ =( a,b,b, Again, we subtract column j from column j +1, j = n 2,n 3,,2 It is easy to see that, step by step, the rows and columns are emptied until finally the determinant c b 0 0 c + a a det(t ˆα, ˆβ(n c +2a a (a + bi (n 2 (n 2 = det, c +(n 2a a c +(a a 0 0 is obtained The proposition follows now immediately, by expanding the determinant along the last row 43 Fibonacci and Lucas Numbers as Leading Principal Minors of a Quasi-Pascal Matrix There are several infinite matrices for which the leading principal minors form a Fibonacci or Lucas (subsequence For instance, in [10], we have presented a family of tridiagonal matrices with the following form: 1 λ λ λ λ F λ ( = λ 2 0 ( λ where λ =(λ i i 0 with λ i C = C\{0} Indeed, the leading principal minors of these matrices for every λ form the sequence ( F n+1 (Theorem 1 in [10] Also, n 1 for the special cases λ 0 = λ 1 = {1, 1}, see [3, 4] and [11] In ([12], pp , Strang presents the infinite tridiagonal (Toeplitz matrices: (43 P = T (3,t,0,0,,(3,t,0,0, (, where t = ±1, and it is easy to show that the leading principal minors of T form the subsequence ( F 2n+2 from the Fibonacci sequence As another example, the n 1

15 Volume 18, pp , August ARMoghaddamfarandSMHPooya leading principal minors of the Toeplitz matrices: (44 Q = T (2,1,1,1,,(2,t,0,0, (, where t = ±1, form the sequence ( F n+2 n 1 for t = 1 and the sequence ( F 2n+1 n 1 for t = 1 ([3], Examples 1, 2 We can summarize the above results in the following proposition Proposition 47 ([3, 4, 11, 12] Let n be a natural number, α = (α i i 0 and β =(β i i 0 be two sequences, and let D n be the nth leading principal minor of T α,β ( Then, the following hold (1 If α = β =(1, 1, 0, 0,, thend n = F n+1 (2 If α = β =(3,t,0, 0, 0, where t = ±1, thend n = F 2n+2 (3 If α =(1, 1, 0, 0, and β =(1, 1, 0, 0,, thend n = F n+1 (4 If α =(2, 1, 1, 1, and β =(2, 1, 0, 0,, thend n = F 2n+1 (5 If α =(2, 1, 1, 1, and β =(2, 1, 0, 0,, thend n = F n+2 Let φ = , the golden ratio, and Φ = 1 5 2, the golden ratio conjugate The recent article of Griffin, Stuart and Tsatsomeros [7] gives the following result: Proposition 48 ([7], Lemma 7 For each positive integer n, let P (n =T (1,Φ,Φ,,(1,φ,φ, (n, and Q(n =T (0, Φ, Φ,,(0, φ, φ, (n Then, we have det(p (n = F n+1, and det(q(n = F Using Propositions 47, 48 and Theorem 31, we immediately deduce the following corollary Corollary 49 Let n be a natural number, α =(α i i 0 and β =(β i i 0 be two sequences, and let D n be the nth leading principal minor of P α,β ( Then, the following hold (1 If α i = β i =1+i 1, thend n = F n+1 (2 If α i = β i =3 i, thend n = F 2n+2 (3 If α i = β i =3+i, thend n = F 2n+2 (4 If α i =1 i and β i =1+i, thend n = F n+1

16 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 579 (5 If α i =2 i +1 and β i =2 i, thend n = F 2n+1 (6 If α i =2 i +1 and β i =2+i, thend n = F n+2 (7 If α i =(2 i 1Φ + 1 and β i =(2 i 1φ +1,thenD n = F n+1 (8 If α i =(1 2 i Φ and β i =(1 2 i φ, thend n = F In the sequel, we study together the sequences F and L, and, in order to unify our treatment, we introduce the following useful notations For ε {+, } we let F ε n = F n if ε =+;andf ε n = L n if ε = Theorem 410 Let r be a non-negative integer and s be a positive integer Suppose that F ε φ r,s = 2r+s and ψ r,s = φ r,s Fr+s ε Fε 2r+s F ε r+s Then, the leading principal minors of the following infinite quasi-pascal matrix: F ε r+s ψ r,s 0 0 ψ r,s φ r,s ( 1 r ( 1 r P [r,s] ( = 0 ( 1 r 0 ( 1 r P α,α ( where α =(α i i 0 is an arithmetic sequence with α i = Fr subsequence {Fnr+s ε } n=1 from Fibonacci or Lucas sequences + i ( 1 r, form the Proof Cahill and Narayan in [5] introduce the following quasi-toeplitz matrices: T [r,s] ( = Fr+s ε ψ r,s 0 0 ψ r,s φ r,s ( 1 r 0 0 ( 1 r 0 0 T β,β ( where β =(F r, ( 1 r, 0, 0, Moreover, they show that (45 det(t [r,s] (n = F ε nr+s Now, we decompose the matrix T [r,s] (n as follows: (46 T [r,s] (n = L(nP [r,s] (n L(n t,

17 Volume 18, pp , August ARMoghaddamfarandSMHPooya where L(n =I 2 2 L 1 (n 2 The proof of (46 is similar to the proof of Theorem 31 and we omit it here Now, using (45 and (46, we easily see that det(p [r,s] (n = F ε nr+s, and the proof of theorem is complete 44 Generalized Pascal Triangle Associated With a Constant Sequence In this subsection, we study the generalized Pascal triangle P α,β (n, where one of the sequences α or β is a constant sequence The following result is another consequence of Theorem 31 Corollary 411 Let α = (α i i 0 and β = (β i i 0 be two sequences with α 0 = β 0 = γ Ifα or β is a constant sequence, then we have det(p α,β (n = γ n Proof By Theorem 31, we have det(p α,β (n = det(t ˆα, ˆβ(n But in both cases, thetoeplitzmatrixt ˆα, ˆβ(n is a lower triangular matrix or an upper triangular one with γ on its diagonal This implies the corollary The generalized Pascal triangle P α,α ( associated with the pair of identical sequences α and α, is called the generalized symmetric Pascal triangle associated with α and yields symmetric matrices P α,α (n by considering principal submatrices consisting of the first n rows and columns of P α,α ( For an arbitrary sequence α =(α i i 0 with α 0 = 0, we define α =( α i i 0 where α i =( 1 i α i for all i Then, the generalized Pascal triangle P α, α ( associated with the sequences α and α, is called the generalized skymmetric Pascal triangle associated with α and α, and yields skymmetric matrices P α, α (n by considering principal submatrices consisting of the first n rows and columns of P α, α ( Example 412 Let n 2 be a natural number Then, we have: (i The generalized symmetric Pascal triangle P F,F (n has determinant 2 n 2 (ii The generalized skymmetric Pascal triangle P F, F(n has determinant 2 n 2 All assertions in this example follow from Theorem 31 of [1] However, we will reprove them independently

18 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 581 Proof (i Consider the generalized symmetric Pascal triangle P F,F (n = F F n F F n+1 Now, we apply the following elementary row operations: R i R i R i 1 R i 2, i = n 1,n 2,,2 It is easy to see that F F n+1 det ( P F,F (n (F n 1 = det = det F F n P λ,µ (n 2 where λ =(2, 2, 2,andµ =(2(F 3 1, 2(F 4 1, 2(F 5 1, 2(F 6 1, Now, by Corollary 411, we get ( 0 1 det(p F,F (n = det 1 2 det ( P λ,µ (n 2 = 2 n 2, as desired

19 Volume 18, pp , August ARMoghaddamfarandSMHPooya (ii Here, we consider the generalized skymmetric Pascal triangle P F, F(n = ( 1 F F Similarly, we apply the following elementary column operations: C j C j +C j 1 C j 2, j = n 1,n 2,,2; and we obtain det(p F, F (n = det = det F F n P ν,λ (n 2 F F n+1 2 where λ =(2, 2, 2, Again, by Corollary 411, we get as desired ( 0 1 det(p F, F(n = det 1 0 det ( P ν,λ (n 2 =2 n 2, Example 413 Let n 2 be a natural number Then, the generalized Pascal triangle P F,I (n has determinant ( 1n

20 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 583 Proof Consider the following generalized Pascal triangle: P F,I (n = n! F F n+1 Again, we use the similar elementary row operations as Example 1(i: R i R i R i 1 R i 2, i = n 1,n 2,,2 Therefore, we deduce that n! det ( P F,I (n = det n! = det 0 0 P λ,µ (n where λ =( 1, 1, 1,andµ =( 1, 10, 73, Now, by Corollary 411, we get as desired det ( ( 1 2 P F,I (n =det 1 3 det ( P λ,µ (n 2 =( 1 n 2 =( 1 n, 5 New Matrices Whose Leading Principal Minors Form the Fibonacci Sequence In this section, we present two matrices whose entries are recursively

21 Volume 18, pp , August ARMoghaddamfarandSMHPooya defined, and we show that the leading principal minor sequence of these matrices is the sequence F =(F n n 1 It is worth mentioning that to construct these matrices we use non-homogeneous recurrence relations Theorem 51 Let a, b C and let n be a natural number Let (r i,j i,j 0 be the doubly indexed sequence given by the recurrence r i,j = r i,j 1 + r i 1,j + ai + bj, for i, j 1, and the initial conditions r i,0 = r 0,j =1, i, j 0 Let D n =det(r i,j 0 i,j Then, D n satisfies the following recursion: D 1 =1, (51 D 2 =1+a + b, D n =(1+a + bd abd n 2 (n 3 In particular, we have D n = F n if and only if (a, b {(1, 1, ( 1, 1} Proof LetR(n denotethematrix(r i,j 0 i,j First, we claim that where U(n =L(n t and 1 T (n = R(n =L(n T (n U(n, ω b a ω b a ω b a ω where ω =1+a + b Now it is easy to see that D n =det(t (n and by expansion through the last row of T (n, we obtain (51 In what follows, for convenience, we will let R = R(n, L = L(n, T = T (n and U = U(n Now, for the proof of the claimed factorization we compute the (i, j-entry of L T U, thatis (52 (L T U i,j = L i,r T r,s U s,j s=0 In fact, so as to prove the theorem, we should establish n n R 0 (L T U =R 0 (R =(1, 1,,1,,

22 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 585 C 0 (L T U =C 0 (A =(1, 1,,1, and (53 (L T U i,j =(L T U i 1,j 1 +(L T U i 1,j + ai + bj, for 1 i, j n 1 Let us do the required calculations First, suppose that i = 0 Then, we have (L T U 0,j = L 0,r T r,s U s,j = T 0,s U s,j = U 0,j =1, s=0 and so R 0 (L T U =R 0 (R =(1, 1,,1 Next, assume that j = 0 In this case, we obtain (L T U i,0 = L i,r T r,s U s,0 = L i,r T r,0 = L i,0 =1, s=0 s=0 and hence we have C 0 (L T U =C 0 (R =(1, 1,,1 Finally, we must establish (53 At the moment, let us assume that 1 i, j n 1 In this case we have (L T U i,j = s=0 = = = L i,r T r,s U s,j L i,r T r,0 U 0,j + L i,r T r,0 U 0,j + ( by (36 + = L i,r T r,0 U 0,j + s=1 + s=1 s=1 s=1 L i,r T r,s U s,j 1 L i,r T r,0 U 0,j + s=0 + s=1 r=1 s=1 L i,r T r,s U s,j 1 L i,0 T 0,s U s 1,j 1 L i,r T r,s U s,j L i,r T r,s ( Us 1,j 1 + U s,j 1 L i,r T r,s U s 1,j 1 L i,r T r,s U s 1,j 1 L i,r T r,0 U 0,j 1

23 Volume 18, pp , August ARMoghaddamfarandSMHPooya = + = = L i,r T r,0 U 0,j + s=1 r=1 s=1 L i,0 T 0,s U s 1,j 1 L i,r T r,0 U 0,j + r=1 s=1 +(L T U i,j s=1 L i,r T r,0 U 0,j + s=2 + = s=1 + = r=2 s=2 ( Li 1,r 1 + L i 1,r Tr,s U s 1,j 1 +(L T U i,j 1 L i,r T r,0 U 0,j 1 ( by (34 and (52 L i 1,r 1 T r,s U s 1,j 1 + L i,0 T 0,s U s 1,j 1 L i 1,0 T 1,s U s 1,j 1 + L i,0 T 0,s U s 1,j 1 L i,r T r,0 U 0,j + s=2 r=2 s=2 r=1 s=1 L i,r T r,0 U 0,j 1 L i 1,r 1 T r,s U s 1,j 1 + L i 1,0 T 1,s U s 1,j 1 + +(L T U i,j s=1 L i,r T r,0 U 0,j + r=1 + = r=1 s=1 r=2 s=2 r=1 s=1 r=1 L i 1,r T r,s U s 1,j 1 L i 1,r 1 T r,1 U 0,j 1 L i 1,r T r,s U s 1,j 1 +(L T U i,j 1 L i,r T r,0 U 0,j 1 L i 1,r 1 T r,s U s 1,j 1 + r=1 s=1 L i,0 T 0,s U s 1,j 1 L i 1,r 1 T r,1 U 0,j 1 + r=1 L i 1,r T r,s ( Us,j U s,j 1 L i 1,r 1 T r 1,s 1 U s 1,j 1 s=2 L i 1,r T r,s U s,j +(L T U i,j 1 + s=1 ( by the structure of T + L i,r T r,0 U 0,j + s=2 r=1 s=1 r=1 s=1 r=1 s=1 L i 1,0 T 1,s U s 1,j 1 L i 1,r T r,s U s,j 1 L i,0 T 0,s U s 1,j 1 L i 1,0 T 1,s U s 1,j 1 + L i 1,r T r,s U s,j 1 + r=1 s=0 L i 1,r 1 T r,1 U 0,j 1 L i,r T r,0 U 0,j 1 ( by (36 L i,r T r,0 U 0,j 1 r=1 L i 1,r T r,s U s,j L i 1,r T r,s U s,j 1 +(L T U i,j 1 + s=1 L i 1,r 1 T r,1 U 0,j 1 r=1 L i,r T r,0 U 0,j 1 ( note that Li 1, = U,j 1 =0 L i 1,r T r,0 U 0,j L i,0 T 0,s U s 1,j 1

24 Volume 18, pp , August Generalized Pascal Triangles and Toeplitz Matrices 587 = + = L i,r T r,0 U 0,j + s=0 r=1 L i 1,r T r,s U s,j +(L T U i,j 1 + s=1 L i,r T r,0 U 0,j + r=1 +(L T U i 1,j s=0 +(L T U i,j 1 + s=1 L i 1,r 1 T r,1 U 0,j 1 + s=0 s=2 L i 1,0 T 0,s U s,j L i,0 T 0,s U s 1,j 1 L i 1,r 1 T r,1 U 0,j 1 + L i 1,0 T 0,s U s,j L i,0 T 0,s U s 1,j 1 L i,r T r,0 U 0,j 1 ( by (52 = (L T U i 1,j +(L T U i,j 1 +Ψ i,j, r=1 s=2 L i 1,0 T 1,s U s 1,j 1 r=1 L i 1,r T r,0 U 0,j L i,r T r,0 U 0,j 1 L i 1,0 T 1,s U s 1,j 1 L i 1,r T r,0 U 0,j where Ψ i,j = L i,r T r,0 U 0,j + s=0 r=1 L i 1,0 T 0,s U s,j L i,r T r,0 U 0,j 1 L i 1,r 1 T r,1 U 0,j 1 + r=1 s=2 L i 1,r T r,0 U 0,j + s=1 L i 1,0 T 1,s U s 1,j 1 L i,0 T 0,s U s 1,j 1 But by an easy calculation one can show that Ψ i,j = ai + bj, which implies the first part of theorem The second part of theorem is now obvious Acknowledgments We would like to thank the referee for the careful reading of an earlier version of this article and for suggestions and comments The first author would like to thank IPM for their financial support

25 Volume 18, pp , August ARMoghaddamfarandSMHPooya REFERENCES [1] R Bacher Determinants of matrices related to the Pascal triangle J Théor Nombres Bordeaux, 14:19 41, 2002 [2] R Bacher and R Chapman Symmetric Pascal matrices modulo p European J Combin, 25: , 2004 [3] N D Cahill, J R D Errico, D A Narayan, and J Y Narayan Fibonacci determinants The College Math J, 33: , 2002 [4] N D Cahill, J R D Errico, and J P Spence Complex factorizations of the Fibonacci and Lucas numbers Fibonacci Quart, 41:13 19, 2003 [5] N D Cahill and D A Narayan Fibonacci and Lucas numbers as tridiagonal matrix determinants Fibonacci Quart, 42: , 2004 [6] A Edelman and G Strang Pascal matrices Amer Math Monthly, 111: , 2004 [7] KGriffin,JLStuart,andMJTsatsomeros NoncirculantToeplitzmatricesallofwhose powers are Toeplitz Czechoslovak Math J, 58: , 2008 [8] Integer-sequences, njas/sequences/indexhtml [9] C Krattenthaler Evaluations of some determinants of matrices related to the Pascal triangle Sém Lothar Combin, Article B47g, 19 pp, 2001/2002 [10] A R Moghaddamfar, S M H Pooya, S Navid Salehy, and S Nima Salehy Fibonacci and Lucas sequences as the principal minors of some infinite matrices J Algebra Appl, (to appear [11] G Strang Introduction to Linear Algebra Third Edition Wellesley-Cambridge Press, 2003 [12] G Strang and K Borre Linear Algebra Geodesy, and GPS Wellesley-Cambridge Press, 1997

More calculations on determinant evaluations

More calculations on determinant evaluations Electronic Journal of Linear Algebra Volume 16 Article 007 More calculations on determinant evaluations A. R. Moghaddamfar moghadam@kntu.ac.ir S. M. H. Pooya S. Navid Salehy S. Nima Salehy Follow this

More information

Certain matrices related to Fibonacci sequence having recursive entries

Certain matrices related to Fibonacci sequence having recursive entries Electronic Journal of Linear Algebra Volume 17 Volume 17 (2008) Article 40 2008 Certain matrices related to Fibonacci sequence having recursive entries Ali Reza Moghaddamfar moghadam@mail.ipm.ir S. Navid

More information

The k-fibonacci matrix and the Pascal matrix

The k-fibonacci matrix and the Pascal matrix Cent Eur J Math 9(6 0 403-40 DOI: 0478/s533-0-0089-9 Central European Journal of Mathematics The -Fibonacci matrix and the Pascal matrix Research Article Sergio Falcon Department of Mathematics and Institute

More information

Transformations Preserving the Hankel Transform

Transformations Preserving the Hankel Transform 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 10 (2007), Article 0773 Transformations Preserving the Hankel Transform Christopher French Department of Mathematics and Statistics Grinnell College Grinnell,

More information

Upper triangular matrices and Billiard Arrays

Upper triangular matrices and Billiard Arrays Linear Algebra and its Applications 493 (2016) 508 536 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Upper triangular matrices and Billiard Arrays

More information

A matrix over a field F is a rectangular array of elements from F. The symbol

A matrix over a field F is a rectangular array of elements from F. The symbol Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted

More information

A Note on the Determinant of Five-Diagonal Matrices with Fibonacci Numbers

A Note on the Determinant of Five-Diagonal Matrices with Fibonacci Numbers Int. J. Contemp. Math. Sciences, Vol. 3, 2008, no. 9, 419-424 A Note on the Determinant of Five-Diagonal Matrices with Fibonacci Numbers Hacı Civciv Department of Mathematics Faculty of Art and Science

More information

MATRICES. a m,1 a m,n A =

MATRICES. a m,1 a m,n A = MATRICES Matrices are rectangular arrays of real or complex numbers With them, we define arithmetic operations that are generalizations of those for real and complex numbers The general form a matrix of

More information

Linear Systems and Matrices

Linear Systems and Matrices Department of Mathematics The Chinese University of Hong Kong 1 System of m linear equations in n unknowns (linear system) 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.......

More information

MCR3U Unit 7 Lesson Notes

MCR3U Unit 7 Lesson Notes 7.1 Arithmetic Sequences Sequence: An ordered list of numbers identified by a pattern or rule that may stop at some number or continue indefinitely. Ex. 1, 2, 4, 8,... Ex. 3, 7, 11, 15 Term (of a sequence):

More information

The Determinant: a Means to Calculate Volume

The Determinant: a Means to Calculate Volume The Determinant: a Means to Calculate Volume Bo Peng August 16, 2007 Abstract This paper gives a definition of the determinant and lists many of its well-known properties Volumes of parallelepipeds are

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY LINEAR ALGEBRA ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,

More information

Lectures on Linear Algebra for IT

Lectures on Linear Algebra for IT Lectures on Linear Algebra for IT by Mgr. Tereza Kovářová, Ph.D. following content of lectures by Ing. Petr Beremlijski, Ph.D. Department of Applied Mathematics, VSB - TU Ostrava Czech Republic 11. Determinants

More information

ACI-matrices all of whose completions have the same rank

ACI-matrices all of whose completions have the same rank ACI-matrices all of whose completions have the same rank Zejun Huang, Xingzhi Zhan Department of Mathematics East China Normal University Shanghai 200241, China Abstract We characterize the ACI-matrices

More information

Elementary maths for GMT

Elementary maths for GMT Elementary maths for GMT Linear Algebra Part 2: Matrices, Elimination and Determinant m n matrices The system of m linear equations in n variables x 1, x 2,, x n a 11 x 1 + a 12 x 2 + + a 1n x n = b 1

More information

Section Summary. Sequences. Recurrence Relations. Summations Special Integer Sequences (optional)

Section Summary. Sequences. Recurrence Relations. Summations Special Integer Sequences (optional) Section 2.4 Section Summary Sequences. o Examples: Geometric Progression, Arithmetic Progression Recurrence Relations o Example: Fibonacci Sequence Summations Special Integer Sequences (optional) Sequences

More information

a 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2.

a 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2. Chapter 1 LINEAR EQUATIONS 11 Introduction to linear equations A linear equation in n unknowns x 1, x,, x n is an equation of the form a 1 x 1 + a x + + a n x n = b, where a 1, a,, a n, b are given real

More information

EVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS

EVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS EVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS CHARLES HELOU AND JAMES A SELLERS Abstract Motivated by a recent work about finite sequences where the n-th term is bounded by n, we evaluate some classes

More information

Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI?

Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI? Section 5. The Characteristic Polynomial Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI? Property The eigenvalues

More information

Linear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.

Linear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0. Matrices Operations Linear Algebra Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0 The rectangular array 1 2 1 4 3 4 2 6 1 3 2 1 in which the

More information

Determinant evaluations for binary circulant matrices

Determinant evaluations for binary circulant matrices Spec Matrices 014; :187 199 Research Article Open Access Christos Kravvaritis* Determinant evaluations for binary circulant matrices Abstract: Determinant formulas for special binary circulant matrices

More information

Chapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations

Chapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations Chapter 1: Systems of linear equations and matrices Section 1.1: Introduction to systems of linear equations Definition: A linear equation in n variables can be expressed in the form a 1 x 1 + a 2 x 2

More information

Math/CS 466/666: Homework Solutions for Chapter 3

Math/CS 466/666: Homework Solutions for Chapter 3 Math/CS 466/666: Homework Solutions for Chapter 3 31 Can all matrices A R n n be factored A LU? Why or why not? Consider the matrix A ] 0 1 1 0 Claim that this matrix can not be factored A LU For contradiction,

More information

Intrinsic products and factorizations of matrices

Intrinsic products and factorizations of matrices Available online at www.sciencedirect.com Linear Algebra and its Applications 428 (2008) 5 3 www.elsevier.com/locate/laa Intrinsic products and factorizations of matrices Miroslav Fiedler Academy of Sciences

More information

1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )

1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i ) Direct Methods for Linear Systems Chapter Direct Methods for Solving Linear Systems Per-Olof Persson persson@berkeleyedu Department of Mathematics University of California, Berkeley Math 18A Numerical

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY 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 information

The 4-periodic spiral determinant

The 4-periodic spiral determinant The 4-periodic spiral determinant Darij Grinberg rough draft, October 3, 2018 Contents 001 Acknowledgments 1 1 The determinant 1 2 The proof 4 *** The purpose of this note is to generalize the determinant

More information

Some constructions of integral graphs

Some constructions of integral graphs Some constructions of integral graphs A. Mohammadian B. Tayfeh-Rezaie School of Mathematics, Institute for Research in Fundamental Sciences (IPM), P.O. Box 19395-5746, Tehran, Iran ali m@ipm.ir tayfeh-r@ipm.ir

More information

Matrices. Chapter Definitions and Notations

Matrices. Chapter Definitions and Notations Chapter 3 Matrices 3. Definitions and Notations Matrices are yet another mathematical object. Learning about matrices means learning what they are, how they are represented, the types of operations which

More information

Math 240 Calculus III

Math 240 Calculus III The Calculus III Summer 2015, Session II Wednesday, July 8, 2015 Agenda 1. of the determinant 2. determinants 3. of determinants What is the determinant? Yesterday: Ax = b has a unique solution when A

More information

SOME RESULTS ON q-analogue OF THE BERNOULLI, EULER AND FIBONACCI MATRICES

SOME RESULTS ON q-analogue OF THE BERNOULLI, EULER AND FIBONACCI MATRICES SOME RESULTS ON -ANALOGUE OF THE BERNOULLI, EULER AND FIBONACCI MATRICES GERALDINE M. INFANTE, JOSÉ L. RAMÍREZ and ADEM ŞAHİN Communicated by Alexandru Zaharescu In this article, we study -analogues of

More information

arxiv: v1 [math.nt] 6 Apr 2018

arxiv: v1 [math.nt] 6 Apr 2018 THE GEOMETRY OF SOME FIBONACCI IDENTITIES IN THE HOSOYA TRIANGLE RIGOBERTO FLÓREZ, ROBINSON A. HIGUITA, AND ANTARA MUKHERJEE arxiv:1804.02481v1 [math.nt] 6 Apr 2018 Abstract. In this paper we explore some

More information

A q-analogue OF THE GENERALIZED FACTORIAL NUMBERS

A q-analogue OF THE GENERALIZED FACTORIAL NUMBERS J. Korean Math. Soc. 47 (2010), No. 3, pp. 645 657 DOI 10.4134/JKMS.2010.47.3.645 A q-analogue OF THE GENERALIZED FACTORIAL NUMBERS Seok-Zun Song, Gi-Sang Cheon, Young-Bae Jun, and LeRoy B. Beasley Abstract.

More information

MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2

MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2 MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a

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

MATH39001 Generating functions. 1 Ordinary power series generating functions

MATH39001 Generating functions. 1 Ordinary power series generating functions MATH3900 Generating functions The reference for this part of the course is generatingfunctionology by Herbert Wilf. The 2nd edition is downloadable free from http://www.math.upenn. edu/~wilf/downldgf.html,

More information

Determinants Chapter 3 of Lay

Determinants Chapter 3 of Lay Determinants Chapter of Lay Dr. Doreen De Leon Math 152, Fall 201 1 Introduction to Determinants Section.1 of Lay Given a square matrix A = [a ij, the determinant of A is denoted by det A or a 11 a 1j

More information

Math 324 Summer 2012 Elementary Number Theory Notes on Mathematical Induction

Math 324 Summer 2012 Elementary Number Theory Notes on Mathematical Induction Math 4 Summer 01 Elementary Number Theory Notes on Mathematical Induction Principle of Mathematical Induction Recall the following axiom for the set of integers. Well-Ordering Axiom for the Integers If

More information

MAT 2037 LINEAR ALGEBRA I web:

MAT 2037 LINEAR ALGEBRA I web: MAT 237 LINEAR ALGEBRA I 2625 Dokuz Eylül University, Faculty of Science, Department of Mathematics web: Instructor: Engin Mermut http://kisideuedutr/enginmermut/ HOMEWORK 2 MATRIX ALGEBRA Textbook: Linear

More information

MATH Topics in Applied Mathematics Lecture 12: Evaluation of determinants. Cross product.

MATH Topics in Applied Mathematics Lecture 12: Evaluation of determinants. Cross product. MATH 311-504 Topics in Applied Mathematics Lecture 12: Evaluation of determinants. Cross product. Determinant is a scalar assigned to each square matrix. Notation. The determinant of a matrix A = (a ij

More information

Refined Inertia of Matrix Patterns

Refined Inertia of Matrix Patterns Electronic Journal of Linear Algebra Volume 32 Volume 32 (2017) Article 24 2017 Refined Inertia of Matrix Patterns Kevin N. Vander Meulen Redeemer University College, kvanderm@redeemer.ca Jonathan Earl

More information

Cullen Numbers in Binary Recurrent Sequences

Cullen Numbers in Binary Recurrent Sequences Cullen Numbers in Binary Recurrent Sequences Florian Luca 1 and Pantelimon Stănică 2 1 IMATE-UNAM, Ap. Postal 61-3 (Xangari), CP 58 089 Morelia, Michoacán, Mexico; e-mail: fluca@matmor.unam.mx 2 Auburn

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY 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 information

Sign Patterns with a Nest of Positive Principal Minors

Sign Patterns with a Nest of Positive Principal Minors Sign Patterns with a Nest of Positive Principal Minors D. D. Olesky 1 M. J. Tsatsomeros 2 P. van den Driessche 3 March 11, 2011 Abstract A matrix A M n (R) has a nest of positive principal minors if P

More information

Fibonacci and Lucas numbers via the determinants of tridiagonal matrix

Fibonacci and Lucas numbers via the determinants of tridiagonal matrix Notes on Number Theory and Discrete Mathematics Print ISSN 30 532, Online ISSN 2367 8275 Vol 24, 208, No, 03 08 DOI: 07546/nntdm2082403-08 Fibonacci and Lucas numbers via the determinants of tridiagonal

More information

Impulse Response Sequences and Construction of Number Sequence Identities

Impulse Response Sequences and Construction of Number Sequence Identities Impulse Response Sequences and Construction of Number Sequence Identities Tian-Xiao He Department of Mathematics Illinois Wesleyan University Bloomington, IL 6170-900, USA Abstract As an extension of Lucas

More information

Computing the Determinant and Inverse of the Complex Fibonacci Hermitian Toeplitz Matrix

Computing the Determinant and Inverse of the Complex Fibonacci Hermitian Toeplitz Matrix British Journal of Mathematics & Computer Science 9(6: -6 206; Article nobjmcs30398 ISSN: 223-085 SCIENCEDOMAIN international wwwsciencedomainorg Computing the Determinant and Inverse of the Complex Fibonacci

More information

GENERALIZATION OF AN IDENTITY OF ANDREWS

GENERALIZATION OF AN IDENTITY OF ANDREWS GENERALIZATION OF AN IDENTITY OF ANDREWS arxiv:math/060480v1 [math.co 1 Apr 006 Eduardo H. M. Brietze Instituto de Matemática UFRGS Caixa Postal 15080 91509 900 Porto Alegre, RS, Brazil email: brietze@mat.ufrgs.br

More information

Fundamentals of Engineering Analysis (650163)

Fundamentals of Engineering Analysis (650163) Philadelphia University Faculty of Engineering Communications and Electronics Engineering Fundamentals of Engineering Analysis (6563) Part Dr. Omar R Daoud Matrices: Introduction DEFINITION A matrix is

More information

MTH 309 Supplemental Lecture Notes Based on Robert Messer, Linear Algebra Gateway to Mathematics

MTH 309 Supplemental Lecture Notes Based on Robert Messer, Linear Algebra Gateway to Mathematics MTH 309 Supplemental Lecture Notes Based on Robert Messer, Linear Algebra Gateway to Mathematics Ulrich Meierfrankenfeld Department of Mathematics Michigan State University East Lansing MI 48824 meier@math.msu.edu

More information

The generalized order-k Fibonacci Pell sequence by matrix methods

The generalized order-k Fibonacci Pell sequence by matrix methods Journal of Computational and Applied Mathematics 09 (007) 33 45 wwwelseviercom/locate/cam The generalized order- Fibonacci Pell sequence by matrix methods Emrah Kilic Mathematics Department, TOBB University

More information

On the Shifted Product of Binary Recurrences

On the Shifted Product of Binary Recurrences 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 13 (2010), rticle 10.6.1 On the Shifted Product of Binary Recurrences Omar Khadir epartment of Mathematics University of Hassan II Mohammedia, Morocco

More information

On the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination

On the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination On the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination J.M. Peña 1 Introduction Gaussian elimination (GE) with a given pivoting strategy, for nonsingular matrices

More information

The Fibonacci Identities of Orthogonality

The Fibonacci Identities of Orthogonality The Fibonacci Identities of Orthogonality Kyle Hawins, Ursula Hebert-Johnson and Ben Mathes January 14, 015 Abstract In even dimensions, the orthogonal projection onto the two dimensional space of second

More information

Math Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88

Math Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88 Math Camp 2010 Lecture 4: Linear Algebra Xiao Yu Wang MIT Aug 2010 Xiao Yu Wang (MIT) Math Camp 2010 08/10 1 / 88 Linear Algebra Game Plan Vector Spaces Linear Transformations and Matrices Determinant

More information

TOPOLOGICAL COMPLEXITY OF 2-TORSION LENS SPACES AND ku-(co)homology

TOPOLOGICAL COMPLEXITY OF 2-TORSION LENS SPACES AND ku-(co)homology TOPOLOGICAL COMPLEXITY OF 2-TORSION LENS SPACES AND ku-(co)homology DONALD M. DAVIS Abstract. We use ku-cohomology to determine lower bounds for the topological complexity of mod-2 e lens spaces. In the

More information

Symmetric functions and the Vandermonde matrix

Symmetric functions and the Vandermonde matrix J ournal of Computational and Applied Mathematics 172 (2004) 49-64 Symmetric functions and the Vandermonde matrix Halil Oruç, Hakan K. Akmaz Department of Mathematics, Dokuz Eylül University Fen Edebiyat

More information

A hyperfactorial divisibility Darij Grinberg *long version*

A hyperfactorial divisibility Darij Grinberg *long version* A hyperfactorial divisibility Darij Grinberg *long version* Let us define a function H : N N by n H n k! for every n N Our goal is to prove the following theorem: Theorem 0 MacMahon We have H b + c H c

More information

Positive entries of stable matrices

Positive entries of stable matrices Positive entries of stable matrices Shmuel Friedland Department of Mathematics, Statistics and Computer Science, University of Illinois at Chicago Chicago, Illinois 60607-7045, USA Daniel Hershkowitz,

More information

Logic and Discrete Mathematics. Section 6.7 Recurrence Relations and Their Solution

Logic and Discrete Mathematics. Section 6.7 Recurrence Relations and Their Solution Logic and Discrete Mathematics Section 6.7 Recurrence Relations and Their Solution Slides version: January 2015 Definition A recurrence relation for a sequence a 0, a 1, a 2,... is a formula giving a n

More information

Week 9-10: Recurrence Relations and Generating Functions

Week 9-10: Recurrence Relations and Generating Functions Week 9-10: Recurrence Relations and Generating Functions April 3, 2017 1 Some number sequences An infinite sequence (or just a sequence for short is an ordered array a 0, a 1, a 2,..., a n,... of countably

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY LINEAR ALGEBRA ELEMENTARY LINEAR ALGEBRA K. R. MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Corrected Version, 7th April 013 Comments to the author at keithmatt@gmail.com Chapter 1 LINEAR EQUATIONS 1.1

More information

RECOUNTING DETERMINANTS FOR A CLASS OF HESSENBERG MATRICES. Arthur T. Benjamin Department of Mathematics, Harvey Mudd College, Claremont, CA 91711

RECOUNTING DETERMINANTS FOR A CLASS OF HESSENBERG MATRICES. Arthur T. Benjamin Department of Mathematics, Harvey Mudd College, Claremont, CA 91711 INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 7 (2007), #A55 RECOUNTING DETERMINANTS FOR A CLASS OF HESSENBERG MATRICES Arthur T. Benjamin Department of Mathematics, Harvey Mudd College,

More information

On Systems of Diagonal Forms II

On Systems of Diagonal Forms II On Systems of Diagonal Forms II Michael P Knapp 1 Introduction In a recent paper [8], we considered the system F of homogeneous additive forms F 1 (x) = a 11 x k 1 1 + + a 1s x k 1 s F R (x) = a R1 x k

More information

Spectrally arbitrary star sign patterns

Spectrally arbitrary star sign patterns Linear Algebra and its Applications 400 (2005) 99 119 wwwelseviercom/locate/laa Spectrally arbitrary star sign patterns G MacGillivray, RM Tifenbach, P van den Driessche Department of Mathematics and Statistics,

More information

Discrete Applied Mathematics

Discrete Applied Mathematics Discrete Applied Mathematics 157 (2009 1696 1701 Contents lists available at ScienceDirect Discrete Applied Mathematics journal homepage: www.elsevier.com/locate/dam Riordan group involutions and the -sequence

More information

Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems

Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September

More information

DISTRIBUTION OF FIBONACCI AND LUCAS NUMBERS MODULO 3 k

DISTRIBUTION OF FIBONACCI AND LUCAS NUMBERS MODULO 3 k DISTRIBUTION OF FIBONACCI AND LUCAS NUMBERS MODULO 3 k RALF BUNDSCHUH AND PETER BUNDSCHUH Dedicated to Peter Shiue on the occasion of his 70th birthday Abstract. Let F 0 = 0,F 1 = 1, and F n = F n 1 +F

More information

Background on Linear Algebra - Lecture 2

Background on Linear Algebra - Lecture 2 Background on Linear Algebra - Lecture September 6, 01 1 Introduction Recall from your math classes the notion of vector spaces and fields of scalars. We shall be interested in finite dimensional vector

More information

Math 3108: Linear Algebra

Math 3108: Linear Algebra Math 3108: Linear Algebra Instructor: Jason Murphy Department of Mathematics and Statistics Missouri University of Science and Technology 1 / 323 Contents. Chapter 1. Slides 3 70 Chapter 2. Slides 71 118

More information

Application of Theorems of Schur and Albert

Application of Theorems of Schur and Albert University of South Carolina Scholar Commons Faculty Publications Mathematics, Department of 9-1-1976 Application of Theorems of Schur and Albert Thomas L. Markham University of South Carolina - Columbia,

More information

Equality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.

Equality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same. Introduction Matrix Operations Matrix: An m n matrix A is an m-by-n array of scalars from a field (for example real numbers) of the form a a a n a a a n A a m a m a mn The order (or size) of A is m n (read

More information

Counting Matrices Over a Finite Field With All Eigenvalues in the Field

Counting Matrices Over a Finite Field With All Eigenvalues in the Field Counting Matrices Over a Finite Field With All Eigenvalues in the Field Lisa Kaylor David Offner Department of Mathematics and Computer Science Westminster College, Pennsylvania, USA kaylorlm@wclive.westminster.edu

More information

Section Summary. Definition of a Function.

Section Summary. Definition of a Function. Section 2.3 Section Summary Definition of a Function. Domain, Codomain Image, Preimage Injection, Surjection, Bijection Inverse Function Function Composition Graphing Functions Floor, Ceiling, Factorial

More information

1300 Linear Algebra and Vector Geometry

1300 Linear Algebra and Vector Geometry 1300 Linear Algebra and Vector Geometry R. Craigen Office: MH 523 Email: craigenr@umanitoba.ca May-June 2017 Matrix Inversion Algorithm One payoff from this theorem: It gives us a way to invert matrices.

More information

Leapfrog Constructions: From Continuant Polynomials to Permanents of Matrices

Leapfrog Constructions: From Continuant Polynomials to Permanents of Matrices Leapfrog Constructions: From Continuant Polynomials to Permanents of Matrices Alberto Facchini Dipartimento di Matematica Università di Padova Padova, Italy facchini@mathunipdit André Leroy Faculté Jean

More information

Spectral inequalities and equalities involving products of matrices

Spectral inequalities and equalities involving products of matrices Spectral inequalities and equalities involving products of matrices Chi-Kwong Li 1 Department of Mathematics, College of William & Mary, Williamsburg, Virginia 23187 (ckli@math.wm.edu) Yiu-Tung Poon Department

More information

MATRICES ARE SIMILAR TO TRIANGULAR MATRICES

MATRICES ARE SIMILAR TO TRIANGULAR MATRICES MATRICES ARE SIMILAR TO TRIANGULAR MATRICES 1 Complex matrices Recall that the complex numbers are given by a + ib where a and b are real and i is the imaginary unity, ie, i 2 = 1 In what we describe below,

More information

Iyad T. Abu-Jeib and Thomas S. Shores

Iyad T. Abu-Jeib and Thomas S. Shores NEW ZEALAND JOURNAL OF MATHEMATICS Volume 3 (003, 9 ON PROPERTIES OF MATRIX I ( OF SINC METHODS Iyad T. Abu-Jeib and Thomas S. Shores (Received February 00 Abstract. In this paper, we study the determinant

More information

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES

AN 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 information

A PERIODIC APPROACH TO PLANE PARTITION CONGRUENCES

A PERIODIC APPROACH TO PLANE PARTITION CONGRUENCES A PERIODIC APPROACH TO PLANE PARTITION CONGRUENCES MATTHEW S. MIZUHARA, JAMES A. SELLERS, AND HOLLY SWISHER Abstract. Ramanujan s celebrated congruences of the partition function p(n have inspired a vast

More information

ON QUADRAPELL NUMBERS AND QUADRAPELL POLYNOMIALS

ON QUADRAPELL NUMBERS AND QUADRAPELL POLYNOMIALS Hacettepe Journal of Mathematics and Statistics Volume 8() (009), 65 75 ON QUADRAPELL NUMBERS AND QUADRAPELL POLYNOMIALS Dursun Tascı Received 09:0 :009 : Accepted 04 :05 :009 Abstract In this paper we

More information

Lecture 12 (Tue, Mar 5) Gaussian elimination and LU factorization (II)

Lecture 12 (Tue, Mar 5) Gaussian elimination and LU factorization (II) Math 59 Lecture 2 (Tue Mar 5) Gaussian elimination and LU factorization (II) 2 Gaussian elimination - LU factorization For a general n n matrix A the Gaussian elimination produces an LU factorization if

More information

GOLDEN SEQUENCES OF MATRICES WITH APPLICATIONS TO FIBONACCI ALGEBRA

GOLDEN SEQUENCES OF MATRICES WITH APPLICATIONS TO FIBONACCI ALGEBRA GOLDEN SEQUENCES OF MATRICES WITH APPLICATIONS TO FIBONACCI ALGEBRA JOSEPH ERCOLANO Baruch College, CUNY, New York, New York 10010 1. INTRODUCTION As is well known, the problem of finding a sequence of

More information

b jσ(j), Keywords: Decomposable numerical range, principal character AMS Subject Classification: 15A60

b jσ(j), Keywords: Decomposable numerical range, principal character AMS Subject Classification: 15A60 On the Hu-Hurley-Tam Conjecture Concerning The Generalized Numerical Range Che-Man Cheng Faculty of Science and Technology, University of Macau, Macau. E-mail: fstcmc@umac.mo and Chi-Kwong Li Department

More information

Lecture Notes 1: Matrix Algebra Part C: Pivoting and Matrix Decomposition

Lecture Notes 1: Matrix Algebra Part C: Pivoting and Matrix Decomposition University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 46 Lecture Notes 1: Matrix Algebra Part C: Pivoting and Matrix Decomposition Peter J. Hammond Autumn 2012, revised Autumn 2014 University

More information

. a m1 a mn. a 1 a 2 a = a n

. a m1 a mn. a 1 a 2 a = a n Biostat 140655, 2008: Matrix Algebra Review 1 Definition: An m n matrix, A m n, is a rectangular array of real numbers with m rows and n columns Element in the i th row and the j th column is denoted by

More information

SCHUR IDEALS AND HOMOMORPHISMS OF THE SEMIDEFINITE CONE

SCHUR IDEALS AND HOMOMORPHISMS OF THE SEMIDEFINITE CONE SCHUR IDEALS AND HOMOMORPHISMS OF THE SEMIDEFINITE CONE BABHRU JOSHI AND M. SEETHARAMA GOWDA Abstract. We consider the semidefinite cone K n consisting of all n n real symmetric positive semidefinite matrices.

More information

22m:033 Notes: 3.1 Introduction to Determinants

22m:033 Notes: 3.1 Introduction to Determinants 22m:033 Notes: 3. Introduction to Determinants Dennis Roseman University of Iowa Iowa City, IA http://www.math.uiowa.edu/ roseman October 27, 2009 When does a 2 2 matrix have an inverse? ( ) a a If A =

More information

DYNAMICS OF THE ZEROS OF FIBONACCI POLYNOMIALS M. X. He Nova Southeastern University, Fort Lauderdale, FL. D, Simon

DYNAMICS OF THE ZEROS OF FIBONACCI POLYNOMIALS M. X. He Nova Southeastern University, Fort Lauderdale, FL. D, Simon M. X. He Nova Southeastern University, Fort Lauderdale, FL D, Simon Nova Southeastern University, Fort Lauderdale, FL P. E. Ricci Universita degli Studi di Roma "La Sapienza," Rome, Italy (Submitted December

More information

CALCULATING EXACT CYCLE LENGTHS IN THE GENERALIZED FIBONACCI SEQUENCE MODULO p

CALCULATING EXACT CYCLE LENGTHS IN THE GENERALIZED FIBONACCI SEQUENCE MODULO p CALCULATING EXACT CYCLE LENGTHS IN THE GENERALIZED FIBONACCI SEQUENCE MODULO p DOMINIC VELLA AND ALFRED VELLA. Introduction The cycles that occur in the Fibonacci sequence {F n } n=0 when it is reduced

More information

Series of Error Terms for Rational Approximations of Irrational Numbers

Series of Error Terms for Rational Approximations of Irrational Numbers 2 3 47 6 23 Journal of Integer Sequences, Vol. 4 20, Article..4 Series of Error Terms for Rational Approximations of Irrational Numbers Carsten Elsner Fachhochschule für die Wirtschaft Hannover Freundallee

More information

MULTISECTION METHOD AND FURTHER FORMULAE FOR π. De-Yin Zheng

MULTISECTION METHOD AND FURTHER FORMULAE FOR π. De-Yin Zheng Indian J. pure appl. Math., 39(: 37-55, April 008 c Printed in India. MULTISECTION METHOD AND FURTHER FORMULAE FOR π De-Yin Zheng Department of Mathematics, Hangzhou Normal University, Hangzhou 30036,

More information

A New Shuffle Convolution for Multiple Zeta Values

A New Shuffle Convolution for Multiple Zeta Values January 19, 2004 A New Shuffle Convolution for Multiple Zeta Values Ae Ja Yee 1 yee@math.psu.edu The Pennsylvania State University, Department of Mathematics, University Park, PA 16802 1 Introduction As

More information

CHAPTER 4 Knot insertion

CHAPTER 4 Knot insertion CHAPTER 4 Knot insertion In Chapter we were led to B-splines, defined via the recurrence relation, as a convenient way to represent spline functions In Chapters 2 and 3 we then established some of the

More information

SPRING OF 2008 D. DETERMINANTS

SPRING OF 2008 D. DETERMINANTS 18024 SPRING OF 2008 D DETERMINANTS In many applications of linear algebra to calculus and geometry, the concept of a determinant plays an important role This chapter studies the basic properties of determinants

More information

EIGENVALUES AND EIGENVECTORS 3

EIGENVALUES AND EIGENVECTORS 3 EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices

More information

q-counting hypercubes in Lucas cubes

q-counting hypercubes in Lucas cubes Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math (2018) 42: 190 203 c TÜBİTAK doi:10.3906/mat-1605-2 q-counting hypercubes in Lucas cubes Elif SAYGI

More information

A SUFFICIENT CONDITION FOR STRICT TOTAL POSITIVITY OF A MATRIX. Thomas Craven and George Csordas

A SUFFICIENT CONDITION FOR STRICT TOTAL POSITIVITY OF A MATRIX. Thomas Craven and George Csordas A SUFFICIENT CONDITION FOR STRICT TOTAL POSITIVITY OF A MATRIX Thomas Craven and George Csordas Abstract. We establish a sufficient condition for strict total positivity of a matrix. In particular, we

More information

#A87 INTEGERS 18 (2018) A NOTE ON FIBONACCI NUMBERS OF EVEN INDEX

#A87 INTEGERS 18 (2018) A NOTE ON FIBONACCI NUMBERS OF EVEN INDEX #A87 INTEGERS 8 (208) A NOTE ON FIBONACCI NUMBERS OF EVEN INDEX Achille Frigeri Dipartimento di Matematica, Politecnico di Milano, Milan, Italy achille.frigeri@polimi.it Received: 3/2/8, Accepted: 0/8/8,

More information