Solutions to the recurrence relation u n+1 = v n+1 + u n v n in terms of Bell polynomials
|
|
- Marilynn Reynolds
- 5 years ago
- Views:
Transcription
1 Volume 31, N. 2, pp , 2012 Copyright 2012 SBMAC ISSN / ISSN (Online) Solutions to the recurrence relation u n+1 = v n+1 + u n v n in terms of Bell polynomials CHRISTOPHER S. WITHERS 1 and SARALEES NADARAJAH 2 1 Applied Mathematics Group, Industrial Research Limited, Lower Hutt, New Zealand 2 School of Mathematics, University of Manchester, Manchester M13 9PL, UK s: c.withers@irl.cri.nz / Saralees.Nadarajah@manchester.ac.uk Abstract. Motivated by time series analysis, we consider the problem of solving the recurrence relation u n+1 = v n+1 + u n v n for n = 0 and u n, given the sequence v n. A solution is given as a Bell polynomial. When v n can be written as a weighted sum of nth powers, then the solution u n also takes this form. Mathematical subject classification: 33E99. Key words: autoregressive processes, Bell polynomials, convolution, maximum. 1 Introduction and summary We define the convolution of sequences {a n, b n : n m} as a n b n = n a j b n j j=0 for n 0. We consider the recurrence equation for u n, u n+1 = v n+1 + u n v n (1.1) for n 0, where u 0 = v 0, and v n is a given sequence. The need for solutions to (1.1) arises with respect to the distribution of the maximum of first order autoregressive processes and more generally to that of #CAM-249/10. Received: 21/VIII/10. Accepted: 02/V/11. Corresponding author.
2 246 SOLUTIONS TO THE RECURRENCE RELATION u n+1 = v n+1 + u n v n autoregressive processes of any order, see Withers and Nadarajah (2010). There are many papers studying the distribution of the maximum of autoregressive processes, see, for example, Chernick and Davis (1982), McCormick and Mathew (1989), McCormick and Park (1992), Borkovec (2000), Ol shanskii (2004) and Elek and Zempléni (2008). However, the results either give the limiting extreme value distributions or assume that the errors come from a specific class (for example, uniform distributed errors, negative binomial errors, A RC H(1) errors, etc). We are aware of no work giving the exact distribution of the maximum of autoregressive processes. As explained in Withers and Nadarajah (2010), solutions to (1.1) will lead to the exact distribution of the maximum of autoregressive processes of any order. The aim of this note to provide solutions for (1.1) accessible for all scientists, not just mathematicians. These solutions are given in terms of Bell polynomials. In-built routines for Bell polynomials are available in most computer algebra packages. For example, see BellB in Mathematica and IncompleteBellPoly in Matlab. So, the solutions given will be accessible to most practitioners. In Section 2, a solution for (1.1) is given as a Bell polynomial. In order to investigate the behavior of u n for large n, an assumption needs to made on the behavior of v n for large n. In Section 3, we show that when v n can be written as a weighted sum of nth powers, then the solution u n also has this form. This assumption is extended in Section 4 to the case when the weights are not constants, but polynomials in n. Some conclusions and future work are noted in Section 5. 2 The solution as a Bell polynomial Theorem 2.1 provides an explicit solution of the recurrence relation (1.1) in terms of the complete ordinary Bell polynomial, B n (w), a function of (w 1,..., w n ), defined for any sequence w = (w 1, w 2,...) by the formal generating function (1 W (t)) 1 = B n (w)t n, where W (t) = n=0 w n t n. n=1
3 CHRISTOPHER S. WITHERS and SARALEES NADARAJAH 247 The solution given by Theorem 2.1 can be computed by a single call to the inbuilt routine, BellB, in Mathematica or some other equivalent computer package. We believe that this is the most direct and the most efficient way to calculate a solution for the recurrence relation (1.1). However, in the absence of computer packages, the complete ordinary Bell polynomial can be calculated using the recurrence relation derived by Theorem 2.2. Theorem 2.1. The recurrence relation (1.1) has the solution: for n 0, where w n = v n 1. u n = B n+1 (w) (2.1) Proof. Multiply (1.1) by t n and sum from n = 0 to obtain (U(t) V (t)) /t = U(t)V (t), where U(t) = u n t n, V (t) = v n t n, n=0 n=0 U(t) = (1 W (t)) 1 V (t) since W (t) = tv (t), and 1 + tu(t) = 1 + (1 W (t)) 1 W (t) = (1 tv (t)) 1. (2.2) Taking the coefficient of t n in the last line gives the explicit solution (2.1). Theorem 2.2. A recurrence relation for b n = B n (w) is given by b n = w n b n for n 1, where w 0 = 0 and b 0 = 1. For example, b 1 = w 0 b 1 + w 1 b 0 = w 1, b 2 = w 0 b 2 + w 1 b 1 + w 2 b 0 = w 1 b 1 + w 2 = w1 2 + w 2, b 3 = w 0 b 3 + w 1 b 2 + w 2 b 1 + w 3 b 0 = w 1 b 2 + w 2 b 1 + w 3 = w w 1w 2 + w 3, and so on.
4 248 SOLUTIONS TO THE RECURRENCE RELATION u n+1 = v n+1 + u n v n Proof. Follows by taking the coefficient of t n in (1 x) 1 1 = x(1 x) 1, where x = W (t). An alternative to the complete ordinary Bell polynomial is the partial ordinary Bell polynomial, B nj (w), defined by for 0 j n. For example, W (t) j = B nj (w)t n n= j B n0 (w) = δ n0, B n1 (w) = w n, B nn (w) = w n 1, where δ 00 = 1 and δ n0 = 0 for n 0. These polynomials are tabled on page 309 of Comtet (1974) for 1 n 10. Recurrence formulas for them are also given by Comtet (1974). They can be computed by a single call to the in-built routine, IncompleteBellPoly, in Matlab. Theorem 2.3 states the relationship between the complete ordinary Bell polynomial and the partial ordinary Bell polynomial. It also provides a recurrence relation for the latter. Corollary 2.1 derives the relationship between u n and v n. Corollary 2.2 derives the reciprocal relationship between v n and u n. Theorem 2.3. We have B n (w) = n B nj (w). (2.3) j=0 A recurrence relation for b nj = B nj (w) is for j 1, j 2 0. For example, for j 1. b n, j+1 = b n, j1 + j 2 = n 1 +n 2 =n n n 1 +n 2 =n b n1 j 1 b n2 j 2 (2.4) n 1 b n1 jw n2 = b kj w n k
5 CHRISTOPHER S. WITHERS and SARALEES NADARAJAH 249 Proof. Take the coefficient of t n in the expansion for (1 W (t)) 1 to obtain (2.3). The recurrence relation (2.4) follows by taking the coefficient of t n in W (t) j 1+ j 2 = W (t) j 1W(t) j 2. The proof is complete. Corollary 2.1. We have u 0 = v 0, u 1 = v 1 + v 2 0, u 2 = v 2 + 2v 0 v 1 + v 3 0, u 3 = v 3 + ( 2v 0 v 2 + v 2 1) + 3v 2 0 v 1 + v 4 0, u 4 = v 4 + (2v 0 v 3 + 2v 1 v 2 ) + ( 3v0 2 v 2 + 3v 0 v1) 2 + 4v 3 0 v 1 + v0 5, u 5 = v 5 + ( ( ) 2v 0 v 4 + 2v 1 v 3 + v2) 2 + 3v 2 0 v 3 + 6v 0 v 1 v 2 + v1 3 + ( 4v 3 0 v 2 + 6v 2 0 v2 1) + 5v 4 0 v 1 + v 6 0, u 6 = v 6 + (2v 0 v 5 + 2v 1 v 4 + 2v 2 v 3 ) + ( 3v0 2 v 4 + 6v 0 v 1 v 3 + 3v 0 v v2 1 v ) 2 + ( 4v0 3 v v0 2 v ( 1v 2 + 4v 0 v1) 3 + 5v 4 0 v v0 1) 3 v2 + 6v 5 0 v 1 + v0 7, u 7 = v 7 + ( ) 2v 0 v 6 + 2v 1 v 5 + 2v 2 v 4 + v3 2 + ( 3v0 2 v 5 + 6v 0 v 1 v 4 + 6v 0 v 2 v 3 + 3v1 2 v ) 3 + 3v 1 v2 2 + ( 4v0 3 v v0 2 v 1v 3 + 6v0 2 v v 0v1 2 v ) 2 + v1 4 + ( 5v0 4 v v0 3 v ) 1v v0 2 v3 1 + ( 6v 5 0 v v 4 0 v2 1) + 7v 6 0 v 1 + v 8 0, u 8 = v 8 + (2v 0 v 7 + 2v 1 v 6 + 2v 2 v 5 + v 3 v 4 ) + ( 3v0 2 v 6 + 6v 0 v 1 v 5 + 6v 0 v 2 v 4 + 3v 0 v v2 1 v ) 4 + 6v 1 v 2 v 3 + v2 3 + ( 4v0 3 v v0 2 v 1v v0 2 v 2v v 0 v1 2 v v 0 v 1 v v3 1 v ) 2 + ( 5v0 4 v v0 3 v 1v v0 3 v v2 0 v2 1 v ) 2 + 5v 0 v1 4 + ( 6v0 5 v v0 4 v ( ) 1v v0 1) 3 v3 + 7v 6 0 v v0 5 v v 7 0 v 1 + v 9 0,
6 250 SOLUTIONS TO THE RECURRENCE RELATION u n+1 = v n+1 + u n v n u 9 = v 9 + ( ) 2v 0 v 8 + 2v 1 v 7 + 2v 2 v 6 + 2v 3 v 5 + v4 2 + ( 3v 2 0 v 7 + 6v 0 v 1 v 6 + 6v 0 v 2 v 5 + 6v 0 v 3 v 4 + 3v 2 1 v 5 + 6v 1 v 2 v 4 + 3v 1 v v2 2 v ( 3) + 4v 3 0 v v0 2 v 1v v0 2 v 2v 4 + 6v0 2 v v 0v1 2 v v 0 v 1 v 2 v 3 + 4v 0 v v3 1 v ) 3 + 6v1 2 v2 2 + ( 5v 4 0 v v 3 0 v 1v v 3 0 v 2v v 2 0 v2 1 v v 2 0 v 1v v 0 v1 3 v ( 2 + v1) 5 + 6v 5 0 v v0 4 v 1v v0 4 v v0 3 v2 1 v ( v0 1) 2 v4 + 7v 6 0 v v0 5 v ) 1v v0 4 v3 1 + ( 8v0 7 v ) v0 6 v v 8 0 v 1 + v0 10. (2.5) Proof. Applying (2.3) to (2.1), and reading the partial polynomials from Comtet s table, we obtain the results. Corollary 2.2. We have v n = B n+1 (x) for n 0, where x n = u n 1. Proof. This follows from W (t) = 1 (1 X (t)) 1, where X (t) = tu(t) = n=1 x nt n. Alternatively, it follows by writing (1.1) as u n+1 = v n+1 + u n v n, n 0, u 0 = v 0, where u n = v n, v n = u n and applying (2.1). 3 The weighted sum of powers solution The solution (2.1) gives no indication of the behavior of u n for large n. To obtain this we need to make an assumption on the behavior of v n for large n. Then if V (t) is tractable, we can obtain u n 1, n 1, as the coefficient of t n in (2.2). Theorems 3.1 and 3.2 illustrate this by two methods. Theorem 3.1. Suppose v n is a weighted sum of powers, at least for large enough n, say r v n = b j ν n j (3.1) for n n 0, where 1 r, n 0 0. Suppose also {δ j } are the roots of r b k ν n 0 k / (δ ν k) = p n0 (δ), (3.2)
7 CHRISTOPHER S. WITHERS and SARALEES NADARAJAH 251 where p n+1 (δ) = δ n+1 v n δ n. Assume that {b j } are all non-zero and that {ν j } are all non-zero and distinct. Assume also that the r roots of (3.2) are all distinct. Then u n has the form u n = J γ j δ n j (3.3) for n m 0, where J I = r + n 0, m 0 = 2n 0 1 if n 0 1 and m 0 = 0 if n 0 = 0. Proof. Having found {δ j }, {γ j } are the roots of J A jn0 (ν)γ j = q n0 (ν) (3.4) for ν = ν 1,..., ν I, where A jn (ν) = δ n j /(δ j ν) and q n+1 (ν) = ν n+1 + u n δ n. Note (3.4) can be written A n0 γ = Q n0, (3.5) where (A n ) kj = A jn (ν k ), Q n = (Q n1,..., Q ni ) and Q nk = q n (ν k ). So, if J = r, a solution is γ = A 1 n 0 Q n0. If r =, numerical solutions can be found by truncating the infinite matrix A n and infinite vectors (Q n, γ ) to N N matrix and N-vectors, then increasing N until the desired precision is reached. The proof, which is by substitution, assumes that {δ j, ν j } are all distinct. The proof relies on the fact that if J a jr n j = 0 for ( 1 n J and r 1 ),..., r J are distinct, then a 1 = = a J = 0, since det r n j : 1 n, j J = 0. If J < r, a solution is given by dropping r J rows of (3.5). If J > r, there are not enough equations for a solution by this method. If n 0 2, the values of u n for n < 2n 0 2 can be found from (1.1) or (2.1) or the extension of (2.5). Now suppose that n 0 1 and n = 2n 0 1. Then s 2 = 0 so that the result remains true. Set f (δ) equal to the left and right hand sides of (3.2). Then a sufficient condition that its roots are distinct, is that f (δ) = 0 implies f (δ) = 0.
8 252 SOLUTIONS TO THE RECURRENCE RELATION u n+1 = v n+1 + u n v n The second and more general method is to compute u n from its generating function via the generating function of v n and (2.2). The advantages of this method are: (i) it always applies, provided that the generating functions exist near t = 0; (ii) we shall see that it becomes clear how to extend the method when multiple roots exist. Theorem 3.2. Suppose (3.1) holds. Then u n has the form u n 1 = R Tj n p j (n) (3.6) for some R, T j and p j (n), a polynomial of some degree n j. Proof. Again we begin with (3.1). The generating function is V (t) = V 0 + V 1, where V 0 = n 0 1 n=0 v nt n, V 1 = r b js n 0 j /(1 s j ), s j = ν j t. Set D = r (1 s j), N = DV 1. Then D(1 tv (t)) = L, where L = D(1 tv 0 ) t N is a polynomial of degree J = r + n 0 > r, assuming that n 0 > 0. So, we can write L = J (1 tt k) say, and 1 + tu(t) = (1 tv (t)) 1 = D/L. Suppose first that {t k } are distinct. Then by the usual partial fractions expansion, D/L = J q k (t k ) / (1 tt k ), where q k (t) = D/L k and L k = L/(1 tt k ). So, u n 1 = J q k (t k ) tk n for n 1. If n 0 = 0 then V 0 = 0, U(t) = N/L = r m k(t k )/(1 tt k ), where m k (t) = M/L k so that u n = r m k (t k ) tk n for n 0. Now suppose first that {t k } are not distinct. Then we can write L = R (1 tt j) n j, where R n j = J and {T k } are distinct. By the general
9 CHRISTOPHER S. WITHERS and SARALEES NADARAJAH 253 partial fraction expansion, (compare Section 2.10 of Gradshteyn and Ryzhik (2007)), n R j D/L = c jk / ( ) k 1 tt j, where c j,n j k+1 = Q (k 1) j (T 1 j )/(k 1)! and Q j (t) = (1 tt j ) n j D/L. So, (3.6) holds, where n j ( ) k + n 1 p j (n) = c jk k 1 for n 1. So, p j (n) is a polynomial of degree n j. If n 0 = 0 just replace D/L by N/L as for the case of distinct roots; then u n is equal to the right hand side of (3.6) for n 0. Corollaries 3.1 to 3.3 deduce the asymptotics of u n and v n as n when (3.1), (3.3) and (3.6) are satisfied. Further particular cases are considered by Corollaries 3.4 to 3.7. Corollary 3.1. If (3.1) holds then as n, where v n B n r n 1 B n = j { b j exp ( iθ j n ) : ν j = r1 } = O(1), r 1 = max I ν j, ν j = r j exp ( ) iθ j. Suppose in addition r = 2 and n 0 = 1. By (3.1), V (t) = 2 b j(1 ν j t) 1 = N(t)/D when max 2 ν jt < 1, where N(t) = b 1 (1 ν 2 t) + b 2 (1 ν 1 t), D = 2 i=0 d it i for d 0 = 1, d 1 = ν 1 ν 2 and d 2 = ν 1 ν 2. So, 1 tv (t) = M(t)/D, where M(t) = D t N(t) = 1 c 1 + c 2 t 2, c 1 = ν 1 ν 2 + b 1 + b 2 and c 2 = ν 1 ν 2 + b 1 ν 2 + b 2 ν 1. Corollary 3.2. If (3.3) holds then u n C n r n 2
10 254 SOLUTIONS TO THE RECURRENCE RELATION u n+1 = v n+1 + u n v n as n, where C n = j { b j exp ( iψ j n ) : δ j = R 1 } = O(1), R 1 = max J δ j, δ j = R j exp ( ) iψ j. Corollary 3.3. If (3.6) holds, set T j = r j exp(iψ j ), r = max r r j and N = max{n j : r j = r}. Then u n 1 r n exp ( ( ) ) N + n 1 inψ j c j N N 1 r j =r,n j =n r n n N 1 D n /(N 1)! for D n = exp ( ) inψ j c j N = O(1). r j =r,n j =n Corollary 3.4. If c 2 = 0 and M(t) has two distinct roots then the solution (3.3) holds with J = 2. Corollary 3.5. Suppose c 2 = 0 = c 1. Then M(t) has one root and (3.3) holds with J = 1. Alternatively, since the right hand side of (2.2) is equal to D(1 c 1 t) 1 1, we obtain u n 1 = c1 n 2 i=0 d ic i 1 for n 1. Corollary 3.6. Suppose c 1 = c 2 = 0. Then M(t) = 1 and the right hand side of (2.2) is equal to D, giving u 0 = d 1, u 1 = d 2 and u n = 0 for n 2. Corollary 3.7. Suppose c 2 = 0 and M(t) has two equal roots, say t = t 1. The root satisfies M(t) = Ṁ(t) = 0, giving t 1 = c 1 /(2c 2 ) = 2/c 1 and 4c 2 = c 2 1. So, M(t) = (1 t/t 1 ) 2 and the right hand side of (2.2) is equal to D ( j + 1)(t/t 1 ) j 1, j=0 giving u n 1 = t1 n{n nd 1t 1 + (n 1)d 2 t1 2 } for n 1. So, in this case, the solution is a weighted power with the weight linear in n.
11 CHRISTOPHER S. WITHERS and SARALEES NADARAJAH An extension to polynomial weights The weighted sum of powers assumption for v n arose naturally in Withers and Nadarajah (2010) assuming that a certain matrix had diagonal Jordan form, or at least if the eigenvalues of the non-diagonal Jordan blocks are zero. When this is not the case, we showed in Withers and Nadarajah (2010) that v n = r i=1 min(n,m i ) 1 k=0 ( n 1 k ) w ik ν n k 1 i (4.1) for n 1, so that v 1 = r i=1 w i0. For this more general case, the method of obtaining u n from v n via its generating function, (2.2), still holds, as shown by Theorem 4.1. Theorem 4.1. Suppose (4.1) holds. Then u n has the form u n 1 = J c k tk n (4.2) for n 1 for some J, c k and t k. Note (4.2) is of the form (3.3) with n 0 = 1. Proof. Setting s = νt and D = d/ds, ( ) t n n 1 ν n k 1 = t k+1 A k, k n=1 where ( ) n 1 A k = s n k 1 = k n=1 D k s n 1 /k! = D k (1 s) 1 /k! = (1 s) k 1 n=1 for s < 1. So, setting s i = ν i t for max r i=1 s i < 1, (4.1) gives the generating function for {v n } V (t) = 1 + V 0 (t) + V 1 (t), where V 0 (t) = ν i =0 m i 1 k=0 w ik t k+1, V 1 (t) = ν i =0 m i 1 k=0 w ik t k+1 (1 s i ) k 1.
12 256 SOLUTIONS TO THE RECURRENCE RELATION u n+1 = v n+1 + u n v n Let {θ j, j = 1,..., R} be the distinct non-zero values of {ν i, i = 1,..., r}. Set n 0 = max {m i : ν i = 0}, J = M n 0, M j = max { } m i : ν i = θ j, R m = max M j = max m i, M = M j, c jk = { } wik : ν i = θ j, j i i N j (t) = M j 1 k=0 t k+1 ( 1 tθ j ) M j k 1 c jk, D(t) = R ( ) M 1 tθ j j. Then V 0 (t) is a polynomial of degree n 0, N j (t) is a polynomial of degree M j, D(t) is a polynomial of degree M, and where V 1 (t) = R N j (t)/ ( ) M 1 tθ j j, D(t)V 1 (t) = R N j (t)q j (t) = N(t) say, Q j (t) is a polynomial of degree M M j, and N(t) is a polynomial of degree M. So, D(t) (1 tv (t)) = D(t) (1 v 0 t tv 0 (t)) t N(t) = L(t) = say, is a polynomial of degree J, giving 1 + tu(t) = (1 tv (t)) 1 = D(t)/ J (1 tt k ). J (1 tt k ) Expanding in partial fractions (see Section 2.10 of Gradshteyn and Ryzhik (2007)) and taking the coefficient of t n gives u n 1 as a weighted sum nth powers of {t k } with the weights polynomials in n. If {t k } are all distinct, then the partial fraction expansion has the form 1 + tu(t) = (1 tv (t)) 1 = J c k / (1 tt k ) so that (4.2) follows. If {t k } are not distinct, then proceed as for the case given in Section 3.
13 CHRISTOPHER S. WITHERS and SARALEES NADARAJAH 257 Suppose that m 1 = 2, ν 1 = 0, (m i, ν i ) = (1, 0) for i > 1. Then R = 1, θ 1 = ν 1, M = M 1 = m 1 = 2, J = 4, M 2 = m 2 = 1 and v n = w 10 ν n (n 1)w 11 ν n 2 1 for n 2. So, setting s = ν 1 t and c 0 = r i=2 w i0, V 0 (t) = c 0 t, V 1 (t) = w 10 t(1 s) 1 + w 11 t 2 (1 s) 2, V (t) = v 0 + c 0 t + N 1 (t)(1 s) 2, N(t) = N 1 (t) = w 10 t(1 s) + w 11 t 2, Q 1 (t) = 1, D(t) (1 tv (t)) = (1 s) 2 ( 1 tv 0 t 2 c 0 ) t N(t) = L(t) = So, if {t i } are distinct, then 4 (1 t i t). i=1 1 + tu(t) = (1 tv (t)) 1 = (1 s) 2 /L(t) = 4 c k / (1 tt k ) say, giving (4.2) with J = 4. For analytic expressions for the roots of a quartic, see Section 3.8.3, page 17 of Abramowitz and Stegun (1964). 5 Conclusions We have solved the recurrence relation u n+1 = v n+1 + u n v n, n = 0 for u n, given the sequence v n. The solution for u n is in terms of Bell polynomials. This form is convenient since in-built routines for computing Bell polynomials are widely available. So, the solutions can be directly applied to derive the distribution of the maximum for autoregressive processes. We have also established the behavior of u n for large n by assuming some form for the behavior of v n for large n. The assumed forms are (3.1), a weighted sum of powers, and (4.1). As shown in Withers and Nadarajah (2010), one of these assumptions always holds. So, the assumptions are not at all restrictive. The work presented can be extended in several ways: 1) consider solving u n+1 = v n+1 + ω n+1 + u n v n + u n ω n, n = 0 for u n, given the sequences v n and ω n ; 2) consider solving u n+1 = v n+1 + ω n+1 + μ n+1 + u n v n + u n ω n + u n μ n, n = 0 for u n, given the sequences v n, ω n and μ n ; and, 3) consider solving multivariate forms of (1.1) taking the form u n+1 = v n+1 + u n v n,
14 258 SOLUTIONS TO THE RECURRENCE RELATION u n+1 = v n+1 + u n v n where the equality and the convolution operator are interpreted element wise. We hope to address some of these problems in a future paper. Acknowledgments. The authors would like to thank the Editor, the Associate Editor and the referee for careful reading and for their comments which greatly improved the paper. REFERENCES [1] M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions. U.S. Department of Commerce, National Bureau of Standards, Applied Mathematics Series, 55 (1964). [2] M. Borkovec, Extremal behavior of the autoregressive process with ARC H(1) errors. Stochastic Processes and Their Applications, 85 (2000), [3] M.R. Chernick and R.A. Davis, Extremes in autoregressive processes with uniform marginal distributions. Statistics and Probability Letters, 1 (1982), [4] L. Comtet, Advanced Combinatorics. Reidel, Dordrecht (1974). [5] P. Elek and A. Zempléni, Tail behaviour and extremes of two-state Markovswitching autoregressive models. Computers and Mathematics with Applications, 55 (2008), [6] I.S. Gradshteyn and I.M. Ryzhik, Tables of Integrals, Series and Products, seventh edition. Academic Press, New York (2007). [7] W.P. McCormick and G. Mathew, Asymptotic results for an extreme value estimator of the autocorrelation coefficient for a first order autoregressive sequence. In: Extreme Value Theory (Oberwolfach, 1987), pp , Lecture Notes in Statistics, 51 (1989), Springer, New York. [8] W.P. McCormick and Y.S. Park, Asymptotic analysis of extremes from autoregressive negative binomial processes. Journal of Applied Probability, 29 (1992), [9] K.A. Ol shanskii, On the extremal index of a thinned autoregression process. (Russian) Vestnik Moskovskogo Universiteta. Seriya I. Matematika, Mekhanika, 70 (2004), Translation in Moscow University Mathematics Bulletin, 59 (2004), [10] C.S. Withers and S. Nadarajah, The distribution of the maximum of a first order autoregressive process: the continuous case. Metrika, doi: /s (2010).
THE nth POWER OF A MATRIX AND APPROXIMATIONS FOR LARGE n. Christopher S. Withers and Saralees Nadarajah (Received July 2008)
NEW ZEALAND JOURNAL OF MATHEMATICS Volume 38 (2008), 171 178 THE nth POWER OF A MATRIX AND APPROXIMATIONS FOR LARGE n Christopher S. Withers and Saralees Nadarajah (Received July 2008) Abstract. When a
More informationEXPANSIONS FOR FUNCTIONS OF DETERMINANTS OF POWER SERIES
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 18, Number 1, Spring 010 EXPANSIONS FOR FUNCTIONS OF DETERMINANTS OF POWER SERIES CHRISTOPHER S. WITHERS AND SARALEES NADARAJAH ABSTRACT. Let Aɛ) be an analytic
More informationEXPANSIONS FOR QUANTILES AND MULTIVARIATE MOMENTS OF EXTREMES FOR HEAVY TAILED DISTRIBUTIONS
REVSTAT Statistical Journal Volume 15, Number 1, January 2017, 25 43 EXPANSIONS FOR QUANTILES AND MULTIVARIATE MOMENTS OF EXTREMES FOR HEAVY TAILED DISTRIBUTIONS Authors: Christopher Withers Industrial
More informationGENERALIZATION OF UNIVERSAL PARTITION AND BIPARTITION THEOREMS. Hacène Belbachir 1
#A59 INTEGERS 3 (23) GENERALIZATION OF UNIVERSAL PARTITION AND BIPARTITION THEOREMS Hacène Belbachir USTHB, Faculty of Mathematics, RECITS Laboratory, Algiers, Algeria hbelbachir@usthb.dz, hacenebelbachir@gmail.com
More informationSimplifying Coefficients in a Family of Ordinary Differential Equations Related to the Generating Function of the Laguerre Polynomials
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Applications and Applied Mathematics: An International Journal (AAM Vol. 13, Issue 2 (December 2018, pp. 750 755 Simplifying Coefficients
More informationEXPLICIT EXPRESSIONS FOR MOMENTS OF χ 2 ORDER STATISTICS
Bulletin of the Institute of Mathematics Academia Sinica (New Series) Vol. 3 (28), No. 3, pp. 433-444 EXPLICIT EXPRESSIONS FOR MOMENTS OF χ 2 ORDER STATISTICS BY SARALEES NADARAJAH Abstract Explicit closed
More informationA Symbolic Operator Approach to Several Summation Formulas for Power Series
A Symbolic Operator Approach to Several Summation Formulas for Power Series T. X. He, L. C. Hsu 2, P. J.-S. Shiue 3, and D. C. Torney 4 Department of Mathematics and Computer Science Illinois Wesleyan
More informationApplications. More Counting Problems. Complexity of Algorithms
Recurrences Applications More Counting Problems Complexity of Algorithms Part I Recurrences and Binomial Coefficients Paths in a Triangle P(0, 0) P(1, 0) P(1,1) P(2, 0) P(2,1) P(2, 2) P(3, 0) P(3,1) P(3,
More informationResearch Article The Approximate Solution of Fredholm Integral Equations with Oscillatory Trigonometric Kernels
Applied Mathematics, Article ID 72327, 7 pages http://ddoiorg/055/204/72327 Research Article The Approimate Solution of Fredholm Integral Equations with Oscillatory Trigonometric Kernels Qinghua Wu School
More informationA BIVARIATE GENERALISATION OF GAMMA DISTRIBUTION
A BIVARIATE GENERALISATION OF GAMMA DISTRIBUTION J. VAN DEN BERG, J.J.J. ROUX and A. BEKKER Department of Statistics, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria,, South
More informationA 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 informationDeterminants of Partition Matrices
journal of number theory 56, 283297 (1996) article no. 0018 Determinants of Partition Matrices Georg Martin Reinhart Wellesley College Communicated by A. Hildebrand Received February 14, 1994; revised
More informationSome Congruences for the Partial Bell Polynomials
3 47 6 3 Journal of Integer Seuences, Vol. 009), Article 09.4. Some Congruences for the Partial Bell Polynomials Miloud Mihoubi University of Science and Technology Houari Boumediene Faculty of Mathematics
More informationArithmetic properties of lacunary sums of binomial coefficients
Arithmetic properties of lacunary sums of binomial coefficients Tamás Mathematics Department Occidental College 29th Journées Arithmétiques JA2015, July 6-10, 2015 Arithmetic properties of lacunary sums
More informationHIGHER-ORDER DIFFERENCES AND HIGHER-ORDER PARTIAL SUMS OF EULER S PARTITION FUNCTION
ISSN 2066-6594 Ann Acad Rom Sci Ser Math Appl Vol 10, No 1/2018 HIGHER-ORDER DIFFERENCES AND HIGHER-ORDER PARTIAL SUMS OF EULER S PARTITION FUNCTION Mircea Merca Dedicated to Professor Mihail Megan on
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationResearch Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method
Applied Mathematics Volume 2012, Article ID 605741, 10 pages doi:10.1155/2012/605741 Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method M. Zarebnia
More informationIntroduction to Matrices
214 Analysis and Design of Feedback Control Systems Introduction to Matrices Derek Rowell October 2002 Modern system dynamics is based upon a matrix representation of the dynamic equations governing the
More informationResearch Article The Adjacency Matrix of One Type of Directed Graph and the Jacobsthal Numbers and Their Determinantal Representation
Applied Mathematics Volume 20, Article ID 423163, 14 pages doi:101155/20/423163 Research Article The Adjacency Matrix of One Type of Directed Graph and the Jacobsthal Numbers and Their Determinantal Representation
More informationInvestigating Geometric and Exponential Polynomials with Euler-Seidel Matrices
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 14 (2011), Article 11.4.6 Investigating Geometric and Exponential Polynomials with Euler-Seidel Matrices Ayhan Dil and Veli Kurt Department of Mathematics
More informationACI-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 informationON SIMILARITIES BETWEEN EXPONENTIAL POLYNOMIALS AND HERMITE POLYNOMIALS
Journal of Applied Mathematics and Computational Mechanics 2013, 12(3), 93-104 ON SIMILARITIES BETWEEN EXPONENTIAL POLYNOMIALS AND HERMITE POLYNOMIALS Edyta Hetmaniok, Mariusz Pleszczyński, Damian Słota,
More informationSums of Products of Bernoulli Numbers*
journal of number theory 60, 234 (996) article no. 00 Sums of Products of Bernoulli Numbers* Karl Dilcher Department of Mathematics, Statistics, and Computing Science, Dalhousie University, Halifax, Nova
More informationCharacteristic Polynomial
Linear Algebra Massoud Malek Characteristic Polynomial Preleminary Results Let A = (a ij ) be an n n matrix If Au = λu, then λ and u are called the eigenvalue and eigenvector of A, respectively The eigenvalues
More informationarxiv: v1 [math.co] 3 Nov 2014
SPARSE MATRICES DESCRIBING ITERATIONS OF INTEGER-VALUED FUNCTIONS BERND C. KELLNER arxiv:1411.0590v1 [math.co] 3 Nov 014 Abstract. We consider iterations of integer-valued functions φ, which have no fixed
More informationBESSEL FUNCTIONS APPENDIX D
APPENDIX D BESSEL FUNCTIONS D.1 INTRODUCTION Bessel functions are not classified as one of the elementary functions in mathematics; however, Bessel functions appear in the solution of many physical problems
More informationFast Polynomial Multiplication
Fast Polynomial Multiplication Marc Moreno Maza CS 9652, October 4, 2017 Plan Primitive roots of unity The discrete Fourier transform Convolution of polynomials The fast Fourier transform Fast convolution
More informationReceived: 2/7/07, Revised: 5/25/07, Accepted: 6/25/07, Published: 7/20/07 Abstract
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 7 2007, #A34 AN INVERSE OF THE FAÀ DI BRUNO FORMULA Gottlieb Pirsic 1 Johann Radon Institute of Computational and Applied Mathematics RICAM,
More informationAssignment 10. Arfken Show that Stirling s formula is an asymptotic expansion. The remainder term is. B 2n 2n(2n 1) x1 2n.
Assignment Arfken 5.. Show that Stirling s formula is an asymptotic expansion. The remainder term is R N (x nn+ for some N. The condition for an asymptotic series, lim x xn R N lim x nn+ B n n(n x n B
More informationMulti Stage Queuing Model in Level Dependent Quasi Birth Death Process
International Journal of Statistics and Systems ISSN 973-2675 Volume 12, Number 2 (217, pp. 293-31 Research India Publications http://www.ripublication.com Multi Stage Queuing Model in Level Dependent
More informationResearch Article A New Class of Meromorphic Functions Associated with Spirallike Functions
Applied Mathematics Volume 202, Article ID 49497, 2 pages doi:0.55/202/49497 Research Article A New Class of Meromorphic Functions Associated with Spirallike Functions Lei Shi, Zhi-Gang Wang, and Jing-Ping
More informationCONVOLUTION TREES AND PASCAL-T TRIANGLES. JOHN C. TURNER University of Waikato, Hamilton, New Zealand (Submitted December 1986) 1.
JOHN C. TURNER University of Waikato, Hamilton, New Zealand (Submitted December 986). INTRODUCTION Pascal (6-66) made extensive use of the famous arithmetical triangle which now bears his name. He wrote
More informationAn Application of Catalan Numbers on Cayley Tree of Order 2: Single Polygon Counting
BULLETIN of the Malaysian Mathematical Sciences Society http://math.usm.my/bulletin Bull. Malays. Math. Sci. Soc. (2) 31(2) (2008), 175 183 An Application of Catalan Numbers on Cayley Tree of Order 2:
More informationMade available courtesy of Elsevier:
A note on processes with random stationary increments By: Haimeng Zhang, Chunfeng Huang Zhang, H. and Huang, C. (2014). A Note on Processes with Random Stationary Increments. Statistics and Probability
More informationThe Truncated Pentagonal Number Theorem
The Truncated Pentagonal Number Theorem George E. Andrews Department of Mathematics The Pennsylvania State University University Park, PA 16802 USA Mircea Merca Doctoral School in Applied Mathematics University
More informationOn the eigenvalues of Euclidean distance matrices
Volume 27, N. 3, pp. 237 250, 2008 Copyright 2008 SBMAC ISSN 00-8205 www.scielo.br/cam On the eigenvalues of Euclidean distance matrices A.Y. ALFAKIH Department of Mathematics and Statistics University
More informationEVALUATION 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 informationMULTISECTION 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 informationThe Generating Functions for Pochhammer
The Generating Functions for Pochhammer Symbol { }, n N Aleksandar Petoević University of Novi Sad Teacher Training Faculty, Department of Mathematics Podgorička 4, 25000 Sombor SERBIA and MONTENEGRO Email
More informationRecurrence Relations between Symmetric Polynomials of n-th Order
Applied Mathematical Sciences, Vol. 8, 214, no. 15, 5195-522 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.214.47525 Recurrence Relations between Symmetric Polynomials of n-th Order Yuriy
More information1. Introduction The incomplete beta function I(a, b, x) is defined by [1, p.269, Eq and p.944, Eq ]
MATHEMATICS OF COMPUTATION Volume 65, Number 215 July 1996, Pages 1283 1288 AN ASYMPTOTIC EXPANSION FOR THE INCOMPLETE BETA FUNCTION B.G.S. DOMAN Abstract. A new asymptotic expansion is derived for the
More informationarxiv:quant-ph/ v1 20 Apr 1995
Combinatorial Computation of Clebsch-Gordan Coefficients Klaus Schertler and Markus H. Thoma Institut für Theoretische Physik, Universität Giessen, 3539 Giessen, Germany (February, 008 The addition of
More informationCOMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE
BIT 6-85//41-11 $16., Vol. 4, No. 1, pp. 11 118 c Swets & Zeitlinger COMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE WALTER GAUTSCHI Department of Computer Sciences,
More informationMULTIPLIERS OF THE TERMS IN THE LOWER CENTRAL SERIES OF THE LIE ALGEBRA OF STRICTLY UPPER TRIANGULAR MATRICES. Louis A. Levy
International Electronic Journal of Algebra Volume 1 (01 75-88 MULTIPLIERS OF THE TERMS IN THE LOWER CENTRAL SERIES OF THE LIE ALGEBRA OF STRICTLY UPPER TRIANGULAR MATRICES Louis A. Levy Received: 1 November
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationResearch Article Localization and Perturbations of Roots to Systems of Polynomial Equations
International Mathematics and Mathematical Sciences Volume 2012, Article ID 653914, 9 pages doi:10.1155/2012/653914 Research Article Localization and Perturbations of Roots to Systems of Polynomial Equations
More informationLetting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc.
1 Polynomial Matrices 1.1 Polynomials Letting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc., n ws ( ) as a
More informationCounting 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 informationMULTI-RESTRAINED STIRLING NUMBERS
MULTI-RESTRAINED STIRLING NUMBERS JI YOUNG CHOI DEPARTMENT OF MATHEMATICS SHIPPENSBURG UNIVERSITY SHIPPENSBURG, PA 17257, U.S.A. Abstract. Given positive integers n, k, and m, the (n, k)-th m- restrained
More informationOn the stirling expansion into negative powers of a triangular number
MATHEMATICAL COMMUNICATIONS 359 Math. Commun., Vol. 5, No. 2, pp. 359-364 200) On the stirling expansion into negative powers of a triangular number Cristinel Mortici, Department of Mathematics, Valahia
More informationNumerical Analysis: Solving Systems of Linear Equations
Numerical Analysis: Solving Systems of Linear Equations Mirko Navara http://cmpfelkcvutcz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office
More informationCayley-Hamilton Theorem
Cayley-Hamilton Theorem Massoud Malek In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n Let A be an n n matrix Although det (λ I n A
More informationADVANCED PROBLEMS AND SOLUTIONS
Edited by RAYMOND E. WHITNEY Please send all communications concerning to RAYMOND E. WHITNEY, MATHEMATICS DEPARTMENT, LOCK HAVEN UNIVERSITY, LOCK HAVEN, PA 17745. This department especially welcomes problems
More informationMath Exam 3 Solutions
Math 6 - Exam 3 Solutions Thursday, July 3rd, 0 Recast the following higher-order differential equations into first order systems If the equation is linear, be sure to give the coefficient matrix At and
More informationMatrices and Determinants
Chapter1 Matrices and Determinants 11 INTRODUCTION Matrix means an arrangement or array Matrices (plural of matrix) were introduced by Cayley in 1860 A matrix A is rectangular array of m n numbers (or
More informationPAijpam.eu ADAPTIVE K-S TESTS FOR WHITE NOISE IN THE FREQUENCY DOMAIN Hossein Arsham University of Baltimore Baltimore, MD 21201, USA
International Journal of Pure and Applied Mathematics Volume 82 No. 4 2013, 521-529 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v82i4.2
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationAsymptotics of Integrals of. Hermite Polynomials
Applied Mathematical Sciences, Vol. 4, 010, no. 61, 04-056 Asymptotics of Integrals of Hermite Polynomials R. B. Paris Division of Complex Systems University of Abertay Dundee Dundee DD1 1HG, UK R.Paris@abertay.ac.uk
More informationResearch Article The Laplace Likelihood Ratio Test for Heteroscedasticity
International Mathematics and Mathematical Sciences Volume 2011, Article ID 249564, 7 pages doi:10.1155/2011/249564 Research Article The Laplace Likelihood Ratio Test for Heteroscedasticity J. Martin van
More informationOn prediction and density estimation Peter McCullagh University of Chicago December 2004
On prediction and density estimation Peter McCullagh University of Chicago December 2004 Summary Having observed the initial segment of a random sequence, subsequent values may be predicted by calculating
More informationMATH39001 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 informationResearch Article Parametric Evaluations of the Rogers-Ramanujan Continued Fraction
International Mathematics and Mathematical Sciences Volume 011, Article ID 940839, 11 pages doi:10.1155/011/940839 Research Article Parametric Evaluations of the Rogers-Ramanujan Continued Fraction Nikos
More informationON THE CONTINUED FRACTION EXPANSION OF THE UNIQUE ROOT IN F(p) OF THE EQUATION x 4 + x 2 T x 1/12 = 0 AND OTHER RELATED HYPERQUADRATIC EXPANSIONS
ON THE CONTINUED FRACTION EXPANSION OF THE UNIQUE ROOT IN F(p) OF THE EQUATION x 4 + x 2 T x 1/12 = 0 AND OTHER RELATED HYPERQUADRATIC EXPANSIONS by A. Lasjaunias Abstract.In 1986, Mills and Robbins observed
More informationBias Reduction when Data are Rounded
Bias Reduction when Data are Rounded Christopher S. Withers & Saralees Nadarajah First version: 10 December 2013 Research Report No. 12, 2013, Probability and Statistics Group School of Mathematics, The
More informationON THE MOMENTS OF ITERATED TAIL
ON THE MOMENTS OF ITERATED TAIL RADU PĂLTĂNEA and GHEORGHIŢĂ ZBĂGANU The classical distribution in ruins theory has the property that the sequence of the first moment of the iterated tails is convergent
More information#A69 INTEGERS 13 (2013) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES
#A69 INTEGERS 3 (203) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES William D. Banks Department of Mathematics, University of Missouri, Columbia, Missouri bankswd@missouri.edu Greg Martin Department of
More informationFE 5204 Stochastic Differential Equations
Instructor: Jim Zhu e-mail:zhu@wmich.edu http://homepages.wmich.edu/ zhu/ January 20, 2009 Preliminaries for dealing with continuous random processes. Brownian motions. Our main reference for this lecture
More informationResearch Article New Partition Theoretic Interpretations of Rogers-Ramanujan Identities
International Combinatorics Volume 2012, Article ID 409505, 6 pages doi:10.1155/2012/409505 Research Article New Partition Theoretic Interpretations of Rogers-Ramanujan Identities A. K. Agarwal and M.
More informationELEMENTARY 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 informationBivariate Weibull-power series class of distributions
Bivariate Weibull-power series class of distributions Saralees Nadarajah and Rasool Roozegar EM algorithm, Maximum likelihood estimation, Power series distri- Keywords: bution. Abstract We point out that
More informationDr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review
Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics
More informationStein s Method and Characteristic Functions
Stein s Method and Characteristic Functions Alexander Tikhomirov Komi Science Center of Ural Division of RAS, Syktyvkar, Russia; Singapore, NUS, 18-29 May 2015 Workshop New Directions in Stein s method
More information6.4 Incomplete Beta Function, Student s Distribution, F-Distribution, Cumulative Binomial Distribution
6.4 Incomplete Beta Function, Student s Distribution, F-Distribution, Cumulative Binomial Distribution 29 CITED REFERENCES AND FURTHER READING: Stegun, I.A., and Zucker, R. 974, Journal of Research of
More informationLinear 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 informationChapter 7. Linear Algebra: Matrices, Vectors,
Chapter 7. Linear Algebra: Matrices, Vectors, Determinants. Linear Systems Linear algebra includes the theory and application of linear systems of equations, linear transformations, and eigenvalue problems.
More informationLinks Between Sums Over Paths in Bernoulli s Triangles and the Fibonacci Numbers
Links Between Sums Over Paths in Bernoulli s Triangles and the Fibonacci Numbers arxiv:1611.09181v1 [math.co] 28 Nov 2016 Denis Neiter and Amsha Proag Ecole Polytechnique Route de Saclay 91128 Palaiseau
More informationA 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 informationCHAPTER 3 Further properties of splines and B-splines
CHAPTER 3 Further properties of splines and B-splines In Chapter 2 we established some of the most elementary properties of B-splines. In this chapter our focus is on the question What kind of functions
More informationMulticomponent diffusion in gases and plasma mixtures
High Temperatures ^ High Pressures, 2002, volume 34, pages 109 ^ 116 15 ECTP Proceedings pages 1337 ^ 1344 DOI:10.1068/htwu73 Multicomponent diffusion in gases and plasma mixtures Irina A Sokolova Institute
More informationALTERNATIVE DERIVATION OF SOME REGULAR CONTINUED FRACTIONS
ALTERNATIVE DERIVATION OF SOME REGULAR CONTINUED FRACTIONS R. F. C. WALTERS (Received 21 July 1966, revised 20 February 1967) In this paper we find an expression for g* as the limit of quotients associated
More informationSpherically Symmetric Logistic Distribution
Journal of Multivariate Analysis 7, 2226 (999) Article ID jmva.999.826, available online at httpwww.idealibrary.com on Spherically Symmetric Logistic Distribution Nikolai A. Volodin The Australian Council
More informationOn Certain Sums of Stirling Numbers with Binomial Coefficients
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 18 (2015, Article 15.9.6 On Certain Sums of Stirling Numbers with Binomial Coefficients H. W. Gould Department of Mathematics West Virginia University
More informationThe Expansion of the Confluent Hypergeometric Function on the Positive Real Axis
Applied Mathematical Sciences, Vol. 12, 2018, no. 1, 19-26 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.712351 The Expansion of the Confluent Hypergeometric Function on the Positive Real
More informationUDC Anderson-Darling and New Weighted Cramér-von Mises Statistics. G. V. Martynov
UDC 519.2 Anderson-Darling and New Weighted Cramér-von Mises Statistics G. V. Martynov Institute for Information Transmission Problems of the RAS, Bolshoy Karetny per., 19, build.1, Moscow, 12751, Russia
More informationSOME IDENTITIES INVOLVING GENERALIZED GENOCCHI POLYNOMIALS AND GENERALIZED FIBONACCI-LUCAS SEQUENCES
SOME IDENTITIES INVOLVING GENERALIZED GENOCCHI POLYNOMIALS AND GENERALIZED FIBONACCI-LUCAS SEQUENCES Zhizheng Zhang and Jingyu Jin Department of mathematics, Luoyang Teachers College, Luoyang, Henan, 471022,
More informationAlmost Product Evaluation of Hankel Determinants
Almost Product Evaluation of Hankel Determinants Ömer Eğecioğlu Department of Computer Science University of California, Santa Barbara CA 93106 omer@csucsbedu Timothy Redmond Stanford Medical Informatics,
More informationOn the Convergence of the Summation Formulas Constructed by Using a Symbolic Operator Approach
On the Convergence of the Summation Formulas Constructed by Using a Symbolic Operator Approach Tian-Xiao He 1, Leetsch C. Hsu 2, and Peter J.-S. Shiue 3 1 Department of Mathematics and Computer Science
More informationHIGHER ORDER MULTI-TERM TIME-FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS INVOLVING CAPUTO-FABRIZIO DERIVATIVE
Electronic Journal of Differential Equations, Vol. 27 27), No. 243, pp.. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu HIGHER ORDER MULTI-TERM TIME-FRACTIONAL PARTIAL DIFFERENTIAL
More informationand the compositional inverse when it exists is A.
Lecture B jacques@ucsd.edu Notation: R denotes a ring, N denotes the set of sequences of natural numbers with finite support, is a generic element of N, is the infinite zero sequence, n 0 R[[ X]] denotes
More informationExponential tail behavior for solutions to the homogeneous Boltzmann equation
Exponential tail behavior for solutions to the homogeneous Boltzmann equation Maja Tasković The University of Texas at Austin Young Researchers Workshop: Kinetic theory with applications in physical sciences
More informationLaplace Transform of Spherical Bessel Functions
Laplace Transform of Spherical Bessel Functions arxiv:math-ph/01000v 18 Jan 00 A. Ludu Department of Chemistry and Physics, Northwestern State University, Natchitoches, LA 71497 R. F. O Connell Department
More information1 Determinants. 1.1 Determinant
1 Determinants [SB], Chapter 9, p.188-196. [SB], Chapter 26, p.719-739. Bellow w ll study the central question: which additional conditions must satisfy a quadratic matrix A to be invertible, that is to
More informationOn Kesten s counterexample to the Cramér-Wold device for regular variation
On Kesten s counterexample to the Cramér-Wold device for regular variation Henrik Hult School of ORIE Cornell University Ithaca NY 4853 USA hult@orie.cornell.edu Filip Lindskog Department of Mathematics
More informationSOLUTION OF GENERALIZED LINEAR VECTOR EQUATIONS IN IDEMPOTENT ALGEBRA
, pp. 23 36, 2006 Vestnik S.-Peterburgskogo Universiteta. Matematika UDC 519.63 SOLUTION OF GENERALIZED LINEAR VECTOR EQUATIONS IN IDEMPOTENT ALGEBRA N. K. Krivulin The problem on the solutions of homogeneous
More informationResearch Article Fourier Series of the Periodic Bernoulli and Euler Functions
Abstract and Applied Analysis, Article ID 85649, 4 pages http://dx.doi.org/.55/24/85649 Research Article Fourier Series of the Periodic Bernoulli and Euler Functions Cheon Seoung Ryoo, Hyuck In Kwon, 2
More informationDeterminants. Chia-Ping Chen. Linear Algebra. Professor Department of Computer Science and Engineering National Sun Yat-sen University 1/40
1/40 Determinants Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra About Determinant A scalar function on the set of square matrices
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 1982 NOTES ON MATRIX METHODS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 198 NOTES ON MATRIX METHODS 1. Matrix Algebra Margenau and Murphy, The Mathematics of Physics and Chemistry, Chapter 10, give almost
More informationOn the asymptotics of Bessel functions in the Fresnel regime
On the asymptotics of Bessel functions in the Fresnel regime Zhu Heitman b, James Bremer c, Vladimir Rokhlin b, Bogdan Vioreanu a a Department of Mathematics, University of Michigan b Department of Computer
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationAlgebraic aspects of some Riordan arrays related to binary words avoiding a pattern
Algebraic aspects of some Riordan arrays related to binary words avoiding a pattern D. Merlini and R. Sprugnoli Dipartimento di Sistemi e Informatica viale Morgagni 65, 5034, Firenze, Italia merlini@dsi.unifi.it,
More information