Convergence of a linear recursive sequence

Size: px
Start display at page:

Download "Convergence of a linear recursive sequence"

Transcription

1 int. j. math. educ. sci. technol., 2004 vol. 35, no. 1, Convergence of a linear recursive sequence E. G. TAY*, T. L. TOH, F. M. DONG and T. Y. LEE Mathematics and Mathematics Education, National Institute of Education, Nanyang Technological University, Singapore egtay@nie.edu.sg. (Received 4 November 2002) A necessary and sufficient condition is found for a linear recursive sequence to be convergent, no matter what initial values are given. Its limit is also obtained when the sequence is convergent. Methods from various areas of mathematics are used to obtain the results. 1. Introduction Let k be a positive integer and p 1, p 2,, p k be real numbers. Define a sequence fa n g 1 n¼1 as follows: a nþkþ1 ¼ Xk p i a nþi, n ¼ 0, 1, 2, ð1þ In this paper, we shall study the convergence of fa n g 1 n¼1 under the condition: p i 5 0 for i ¼ 1, 2,, k. Hence in this paper, we always assume that p i 5 0 for i ¼ 1, 2,, k. The convergence of the above sequence was first studied by Zhao et al. [3]. A partial result for k ¼ 3 was obtained. The convergence of the sequence fa n g 1 n¼1 sometimes depends not only on p 1, p 2,, p k, but also on the initial values a 1, a 2,, a k. For example, if k >1, p 1 ¼ 1 and p i ¼ 0 for 2 4 i 4 k, then {a n } is divergent, unless a 1 ¼ a 2 ¼¼a k. But if p k ¼ 1 and p i ¼ 0fori ¼ 1, 2,, k 1, then fa n g 1 n¼1 is always convergent, no matter what the initial values a 1, a 2,, a k are. In this paper, we do not consider special initial values a 1, a 2,, a k. Our purpose is to find a necessary and sufficient condition on p 1, p 2,, p k such that fa n g 1 n¼1 is always convergent, no matter what the initial values a 1, a 2,, a k are. Integers a 1, a 2,, a n, where n 5 2, are called coprime if there is no integer q with q 5 2 such that a i is divisible by q for all i ¼ 1, 2,, n. Our goal in this paper is to prove that this sequence is convergent for all initial values a 1, a 2,, a k if and only if either (i) p 1 þ p 2 þþp k < 1, or (ii) p 1 þ p 2 þþp k ¼ 1 and the integers in the following set are coprime: ¼fkþ1 i : p i 6¼ 0, i ¼ 1, 2,, kg * The author to whom correspondence should be addressed. International Journal of Mathematical Education in Science and Technology ISSN X print/issn online # 2004 Taylor & Francis Ltd DOI: /

2 52 E. G. Tay et al. We also find its limit if the sequence is convergent, P k P i lim a p j a i n ¼ P n!1 k ðk þ 1 iþp i If k ¼ 1, then a n ¼ a 1 p n 1 1 for n 5 1. So the above result is obvious when k ¼ 1. Thus in the following, we always assume that k 5 2. In section 2, we use an analysis method to prove that if p 1 þ p 2 þþp k ¼ 1 and the integers in are coprime, then fa n g 1 n¼1 is always convergent, and its limit is also obtained. In section 3, we prove that fa n g 1 n¼1 is not always convergent if p 1 þ p 2 þþ p k ¼ 1 and the integers in are not coprime. In section 4, we use a matrix method to prove fa n g 1 n¼1 is convergent under the condition that p 1 þ p 2 þþp k ¼ 1 and the integers in are coprime. Its limit is found by another method. In section 5, we finally settle the last two cases: p 1 þ p 2 þþp k < 1 and p 1 þ p 2 þþp k > 1. Thus we more than extend the results of Zhao et al.: we are also motivated by the pedagogical implications of using methods from different areas of mathematics, as stated in [3]. In the teaching of mathematics, an example demonstrating the application of methods or techniques in one branch of mathematics to solve a problem formulated in another field of mathematics can not only help students to see the intrinsic links between different areas of mathematics but also motivate them to develop a passion and interest in their learning of mathematics. 2. The integers in are coprime and p 1 þ p 2 þþp k ¼ 1 Recall that k 5 2 and p i 5 0 for i ¼ 1, 2,, k. In this section, we always assume that p 1 þ p 2 þþp k ¼ 1. We shall show that the sequence fa n g 1 n¼1 is always convergent if all integers in are coprime. The limit of this sequence is also obtained if it is convergent. Define another two sequences fm n g 1 n¼1 and fm ng 1 n¼1 : and M n ¼ maxfa nþi : 0 4 i 4 k 1g n ¼ 1, 2, ð2þ m n ¼ minfa nþi : 0 4 i 4 k 1g n ¼ 1, 2, ð3þ Lemma 2.1. fm n g 1 n¼1 is a descending sequence. Proof: For any n 5 1, by definition, a nþk ¼ Xk p i a n 1þi 4 Xk p i M n ¼ M n Hence M nþ1 ¼ maxfa nþi : 1 4 i 4 kg 4 maxfm n, a nþk g¼m n

3 Convergence of a linear recursive sequence 53 Similarly we can prove the following result. Lemma 2.2. fm n g 1 n¼1 is an ascending sequence. Since M 1 5 M n 5 m n 5 m 1 for all n 5 1, by Lemmas 2.1 and 2.2, the two sequences fm n g 1 n¼1 and fm ng 1 n¼1 are convergent. Thus if the following sequence tends to zero: M 1 m 1, M 2 m 2, M 3 m 3, then fa n g 1 n¼1 is convergent. Let ¼fkþ1 i : p i > 0, 1 4 i 4 kg Let t ¼jj (i.e. the cardinality of ) and ¼fr 1, r 2,, r t g where 1 4 r 1 < r 2 < < r t 4 k. Define ( ) ¼ Xt w i r i : w i is a non-negative integer for i ¼ 1, 2,, t Lemma 2.3. The integers in are coprime if and only if r, r þ 1 2 for some integer r. Proof. (Necessity). Assume that r 1, r 2,, r t are coprime. Then there exist integers w 1, w 2,, w t such that X t w i r i ¼ 1 Let Then r ¼ X ð w i Þ r i 1 4 i 4 t w i <0 r þ 1 ¼ X w i r i 1 4 i 4 t w i >0 Hence r, r þ 1 2. (Sufficiency) Assume that r, r þ 1 2 for some r. Let r ¼ Xt w i r i and r þ 1 ¼ Xt w 0 i r i

4 54 E. G. Tay et al. Then 1 ¼ Xt ðw 0 i w iþr i implying that r 1, r 2,, r t are coprime. Lemma 2.4. The integers in are coprime if and only if fg, g þ 1, g þ 2, g for some positive integer g. Proof. By Lemma 2.3, we just need to prove the necessity of this lemma. Assume that the integers in are coprime. By Lemma 2.3, r, r þ 1 2 for some r. Thus fðr 1 1Þr þ i : i ¼ 0, 1,, r 1 1g since ðr 1 1Þr þ i ¼ iðr þ 1Þþðr 1 i 1Þr 2. Let g ¼ðr 1 1Þr. For any integer n > g, we have n ¼ g þ wr 1 þ i for some non-negative integers w and i with 0 4 i 4 r 1 1. Observe that g þ i 2. Hence n 2. This proves the necessity. By definition, for any integer n, there exist f i (n) for i ¼ 1, 2,, k such that a n ¼ f 1 ðnþa 1 þ f 2 ðnþa 2 þþf k ðnþa k ð4þ where f i (n) is a polynomial of p 1, p 2,, p k with non-negative integer coefficients, and f i (n) is independent of a 1, a 2,, a k. Lemma 2.5. In (4), f i ðnþ 5 0 and f i (n) is unique for all i ¼ 1, 2,, k and n 5 1. Further, if p 1 þ p 2 þþp k ¼ 1, then f 1 ðnþþf 2 ðnþþþf k ðnþ ¼1 Proof. As p j 5 0 for all j with 1 4 j 4 k and f i (n) is a polynomial of p 1, p 2,, p k with non-negative integer coefficients, we have f i ðnþ 5 0. Assume that there exists another expression for a n : a n ¼ f1 0 ðnþa 1 þ f2 0 ðnþa 2 þþfk 0 ðnþa k Then ð f i ðnþ f 0 i ðnþþa i ¼ 0 Letting a i ¼ 1 and a j ¼ 0 for j 6¼ i, we have f i ðnþ ¼fi 0ðnÞ. Hence f iðnþ ¼fi 0 ðnþ for i ¼ 1, 2,, k. Thus the lemma holds. Assume that p 1 þ p 2 þþp k ¼ 1. If a 1 ¼ a 2 ¼¼a k ¼ 1, then we have a n ¼ 1 for all n 5 1 by definition. Hence f 1 ðnþþf 2 ðnþþþf k ðnþ ¼1. Lemma 2.5 shows that f i (n) is well defined for any integers n and i with n 5 1 and 1 4 i 4 k. In the following, we show some properties on f i (n), and these results will be used in the proof of the main result.

5 Convergence of a linear recursive sequence 55 Lemma 2.6. For any integers n and i: n 5 k þ 1 and 1 4 i 4 k, f i ðnþ > 0 if and only if f i ðn r j Þ > 0 for some j with 1 4 j 4 t. Proof: By definition, a n ¼ Xk p s a n kþs 1 s¼1 ¼ Xk s¼1 p s f i ðn k þ s 1Þa i So ¼ Xk s¼1 p s f i ðn k þ s 1Þa i f i ðnþ ¼ Xk s¼1 p s f i ðn k þ s 1Þ ¼ Xt p kþ1 rj f i ðn r j Þ Since p kþ1 rj > 0 and f i ðn r j Þ 5 0 for j ¼ 1, 2,, t, the result follows immediately. Corollary 2.1. If n 5 k and f i ðnþ > 0, then f i ðn þ uþ > 0 for any u 2. Proof. Assume that u > 0. Since n þ r j > k for j ¼ 1, 2,, t, the result follows immediately by Lemma 2.6 repeatedly. Corollary 2.2. Let i and j be any integers with 1 4 i 4 k and 1 4 j 4 t. If i þ r j 5 k þ 1, then f i ði þ u þ r j Þ > 0 for any u 2. Proof. Note that f i ðiþ ¼1 > 0. Since i þ r j > k, we have f i ði þ r j Þ > 0 by Lemma 2.6. Thus the result follows from Corollary 2.1. If the integers in are coprime, then by Lemma 2.4, there exists g 2 such that n 2 if n 5 g. Lemma 2.7. Assume that p 1 > 0 and the integers in are coprime. Let g be the minimum integer such that u 2 for every integer u 5 g. Then f i ðnþ > 0 for i ¼ 1, 2,, k and n 5 g þ 2k. Proof. Since p 1 > 0, we have r t ¼ k. Thus i þ r t ¼ i þ k > k. Let n 0 ¼ n i k 5 g. Then by Corollary 2.2, we have f i ðnþ ¼f i ðn 0 þ i þ r t Þ > 0

6 56 E. G. Tay et al. Define a matrix B ¼ðb i, j Þ kk : (i) b i, iþ1 ¼ 1for, 2,, k 1; (ii) b k, j ¼ p j for j ¼ 1, 2,, k; (iii) b i, j ¼ 0ifi6¼ k or j 6¼ i þ 1. Let X n ¼ða n, a nþ1,, a nþk 1 Þ be a column vector. Then X nþ1 ¼ BX n, n ¼ 1, 2, Thus X nþ1 ¼ B n X 1, n ¼ 1, 2, Lemma 2.8. For any positive integers n and j, Proof. a nþjk ¼ f 1 ðnþa jkþ1 þ f 2 ðnþa jkþ2 þþf k ðnþa ðjþ1þk We just need to show that it holds when j ¼ 1, i.e. a nþk ¼ f 1 ðnþa kþ1 þ f 2 ðnþa kþ2 þþf k ðnþa 2k Observe that X n ¼ B n 1 X 1 and X nþk ¼ B n 1 þ1. By Lemma 2.5, the matrix B n 1 is unique in the expression X n ¼ B n 1 X 1. Hence the result holds. Assume that p 1 > 0 and the integers in are coprime. Let g be the minimum integer such that u 2 for every integer u 5 g. Let f ¼ minf f i ðuþ : 1 4 i 4 k, g þ 2k 4 u 4 g þ 3k 1g ð5þ By Lemma 2.7 and the result P k f iðnþ ¼1, we have 0 < f < 1. Lemma 2.9. Assume that p 1 > 0 and the integers in are coprime. If n 5 g þ 2k, then a n 4 fm jkþ1 þð1 fþm jkþ1 where j ¼bðn gþ=kc 2. Proof. Observe that g þðj þ 2Þk 4 n < g þðj þ 3Þk. Let n 0 ¼ n jk. Then g þ 2k 4 n 0 4 g þ 3k 1. By Lemma 2.8, a n ¼ a n0 þjk ¼ Xk f i ðn 0 Þa jkþi By definition, m jkþ1 ¼ a jkþs for some s with 1 4 s 4 k. Thus a n ¼ f s ðn 0 Þa jkþs þ X f i ðn 0 Þa jkþi 1 4 i 4 k i6¼s 4 f s ðn 0 Þm jkþ1 þ X 1 4 i 4 k i6¼s f i ðn 0 ÞM jkþ1 ¼ f s ðn 0 Þm jkþ1 þð1 f s ðn 0 ÞÞM jkþ1 4 fm jkþ1 þð1 fþm jkþ1 where the last inequality follows from the fact that f 4 f s ðn 0 Þ and m jkþ1 4 M jkþ1. The proof is thus completed.

7 Now we are going to establish the main result in this section. Theorem 2.1. If p 1 þ p 2 þþp k ¼ 1 and the integers in are coprime, then the sequence fa n g 1 n¼1 is always convergent. Proof. There are two cases: p 1 >0 or p 1 ¼ 0. Case 1. p 1 >0. We first prove that for any non-negative integer j 5 0, M gþðjþ2þk m gþðjþ2þk 4 ð1 fþðm jkþ1 m jkþ1 Þ By definition, M gþðjþ2þk ¼ a n for some n with g þðj þ 2Þk 4 n < g þðj þ 3Þk. Then by Lemma 2.9, M gþð jþ2þk ¼ a n 4 fm jkþ1 þð1 fþm jkþ1 So M gþð jþ2þk m gþð jþ2þk 4 fm jkþ1 þð1 fþm jkþ1 m jkþ1 ¼ð1 fþðm jkþ1 m jkþ1 Þ Let c ¼dðg 1Þ=ke. Then g þðj þ 2Þk 4 ðc þ j þ 2Þk þ 1 and M ðcþjþ2þkþ1 m ðcþjþ2þkþ1 4 M gþðjþ2þk m gþðjþ2þk Hence Convergence of a linear recursive sequence 57 4 ð1 fþðm jkþ1 m jkþ1 Þ M n m n 4 ð1 fþ bðn 1Þ=ðcþ2Þkc ðm 1 m 1 Þ As 0 < f < 1, lim ðm n m n Þ¼0 n!1 Therefore the sequence fa n g 1 n¼1 is convergent by the definition of M n and m n. Case 2. p 1 ¼ 0. Let j ¼ minfi : p i 6¼ 0, i ¼ 1, 2,, kg. So24j4k. Then a nþkþ1 ¼ Xk p i a n kþi for any integer n 5 0. Define another sequence {b n }, where b n ¼ a nþj 1 n ¼ 1, 2, 3,. Let k 0 ¼ k j þ 1 and p 0 i ¼ p iþj 1 for i ¼ 1, 2,, k 0. Then i¼j for b nþk0 þ1 ¼ a nþk0 þ1þj 1 ¼ a nþkþ1 ¼ Xk p i a n kþi ¼ Xk0 p 0 i b n k 0 þi i¼j for any integer n 5 0. Observe that p 0 i 5 0 for i ¼ 1, 2,, k0, p 0 1 > 0 and p 0 1 þ p0 2 þþp0 k ¼ 1. It is easy to verify that 0 0 ¼fk 0 i þ 1 : p 0 i > 0, i ¼ 1, 2,, k0 g¼ since p 0 i ¼ p iþj 1 and k 0 ¼ k j þ 1. Thus the integers in 0 are coprime. By the result in case 1, the sequence b 1, b 2, is convergent and thus fa n g 1 n¼1 is convergent.

8 58 E. G. Tay et al. Now we assume that the sequence fa n g 1 n¼1 is convergent and we are going to determine d ¼ lim n!1 a n. Theorem 2.2. to d, then If p 1 þ p 2 þþp k ¼ 1 and the sequence fa n g 1 n¼1 is convergent d ¼ lim n!1 a n ¼ P k P i p j a i P k ðk i þ 1Þp i Proof. It is easy to verify that for any n 5 0, X i p j!ða nþiþ1 a nþi Þ¼0 Consider the sum of the above expression for n ¼ 0, 1, 2,, t. Using this result repeatedly, we have X i p j!ða tþiþ1 a i Þ¼0 When t tends to infinity,! X i 0 ¼ lim p j ða tþiþ1 a i Þ¼ Xk t!1 Hence! d Xk X i p j ¼ Xk X i X i p j!ð lim t!1 a tþiþ1 a i Þ p j!a i Observe that So X i p j! ¼ Xk ðk i þ 1Þ p i P k P i p j a i d ¼ P k ðk i þ 1Þp i 3. The integers in are not coprime and p 1 þ p 2 þþp k ¼ 1 For any u 2 and 1 4 i 4 k, define 8 >< 1, if u ¼ 0 Jði, uþ ¼ 1, if i þ r j > k for some j with w j > 0 >: 0, otherwise where u ¼ P t w jr j.

9 Convergence of a linear recursive sequence 59 Lemma 3.1. For any integers n and i: n 5 1 and 1 4 i 4 k, f i ðnþ > 0 if and only if n ¼ u þ i for some u 2 with Jði, uþ ¼1. Proof: (Necessity) Assume that f i ðnþ > 0. If 1 4 n 4 k, then f i ðnþ > 0 if and only if n ¼ i and so the necessity holds. Now let n > k. By Lemma 2.6, f i ðn r j Þ > 0 for some j. By induction, n r j ¼ i þ u 0 for some u 0 2 with Jði, u 0 Þ¼1. Let u ¼ u 0 þ r j. It is easy to verify that Jði, uþ ¼1. Hence the necessity holds. (Sufficiency) Assume that n ¼ u þ i for some u 2 with Jði, uþ ¼1. It is obvious that fðnþ > 0ifu¼0. If u 6¼ 0, then u ¼ P t w jr j and w j > 0 for some j such that i þ r j > k. Let u 0 ¼ u r j. Then n ¼ u 0 þ r j þ i and u 0 2. By Corollary 2.1, f i ðr j þ iþ > 0. By Lemma 2.6 repeatedly, we have f i ðr j þ i þ u 0 Þ > 0, i.e. f i ðnþ > 0. The sufficiency holds. Let q be the largest common factor of integers in. Lemma 3.2. For 1 4 i 4 k, ifn i is not divisible by q, then f i ðnþ ¼0. Proof. Assume that n i is not divisible by q. By Lemma 3.1, we just need to prove that n i 62. Otherwise, n i ¼ Xt for some non-negative integers w j, j ¼ 1, 2,, t. As q is a common factor of r 1, r 2,, r t, n i is divisible by q, a contradiction. By Lemma 3.2, for any integer l 5 0 and s ¼ 1, 2,, q, w j r j a sþlq ¼ bðk sþ=qc X c¼0 f sþcq ðs þ lqþa sþcq ð6þ So the following result is obtained immediately. Theorem 3.1. If p 1 þ p 2 þþp k ¼ 1 and the integers in are not coprime, then the sequence fa n g 1 n¼1 is not always convergent. Proof. Choose special initial values: for 1 4 i 4 k, ifi 1 is divisible by q, then a i ¼ 1, and otherwise, a i ¼ 0. Then by (6), for any n 5 1, if ðn 1Þ is divisible by q, then a n ¼ 1 and otherwise, a n ¼ 0. This shows that fa n g 1 n¼1 is divergent for some initial values. Hence the theorem holds. 4. By algebraic approach In the following, another method is used to prove the convergence of the sequence fa n g 1 n¼1 if p 1 þ p 2 þþp k ¼ 1 and the integers in are coprime. Recall the column vector X n and a matrix B defined in section 2. We have X nþ1 ¼ BX n, n ¼ 1, 2,

10 60 E. G. Tay et al. and thus X nþ1 ¼ B n X 1, n ¼ 1, 2, Observe that B is a stochastic matrix, since each entry is non-negative and the sum of entries in each row of B is 1. We shall prove that if p 1 þ p 2 þþp k ¼ 1 and the integers in are coprime, then all entries in B s are positive for some s, and thus B is a regular stochastic matrix. Construct another k k matrix A ¼ða i, j Þ kk : a i, j ¼!ðb i, j Þ, i ¼ 1, 2,, k, j ¼ 1, 2,, k where!ðxþ ¼1ifx > 0 and!ðxþ ¼0 otherwise. Lemma 4.1. For any integer s 5 1, the entry of B s at (i, j ) is positive if and only if the entry of A s at (i, j ) is positive. Proof. Observe that all entries in both A s and B s are non-negative for all s. By the definition of A, the lemma can be proved by induction on s. Let D be the digraph with vertex set VðDÞ ¼fx 1, x 2,, x k g, and arc set fx i x j : a i, j ¼ 1, i, j ¼ 1, 2,, kg The following is the digraph D when k ¼ 6 and p 1 > 0, p 3 > 0 and p 4 >0. Then A is an adjacency matrix of D. Let A s ¼ða i, j ðsþþ kk. We refer to the following theorem in Graph Theory (see, for example, [1]). Theorem 4.1. For any s 5 1, a i, j ðsþ is the number of walks in D from x i to x j of length s. Lemma 4.2. If p 1 > 0 and all integers in are coprime, then there exists an integer s > 0 such that all entries in B s are positive. Proof. By Lemma 2.4, n 2 for all n 5 g. We shall show that for any i, j with 1 4 i 4 k and 1 4 j 4 k, D has a walk from x i to x j of length g þ 2k. Then the lemma follows from Lemma 4.1 and Theorem 4.1. We first prove that D has a walk from x k to x j of length g þ k for j ¼ 1, 2,, k. For any i ¼ 1, 2,, t, D has a basic circle of length r i, denoted by C i : x k! x k ri þ1! x k ri þ2!!x k For any j with 1 4 j 4 k, we have g þ k j 5 g. Thus g k þ j ¼ Xt v¼1 w v r v

11 for some non-negative integers w v. Then the following walk from x k to x j is of length g þ k: start from x k ; go along circle C v, w v times, for v ¼ 1, 2,, t; finally take the walk x k! x 1! x 2!!x j : We then prove that for any s > g þ k, D has a walk from x k to x j of length s for j ¼ 1, 2,, k. This can be done by induction. By induction, D has a walk P s 1 ðx k, x j 1 Þ from x k to x j 1 of length s 1, where x 0 ¼ x k. So P s 1 ðx k, x j 1 Þ followed by x j 1 x j is a walk from x k to x j of length s. Now let s ¼ g þ 2k and i, j be integers with 1 4 i 4 k and 1 4 j 4 k. By the above result, D has a walk P from x k to x j of length s k þ i, since s k þ i 5 g þ k. Thus the walk x i! x iþ1!!x k followed by P is a walk from x i to x j of length s. By Theorem 4.1, all entries of A gþ2k and hence all entries of B gþ2k are positive. By Lemmas 4.1 and 4.2, B is a regular stochastic matrix. By a theorem on regular stochastic matrices [2], we have (a) B has a unique fixed probability vector t, i.e., t is a row vector satisfying tb ¼ t (b) the sequence B, B 2, B 3, approaches the matrix T whose rows are each the fixed probability vector t. By (b), the sequence fa n g 1 n¼1 approaches tx 1, where t is the unique probability vector satisfying tb ¼ t. Now we are going to determine the probability vector t satisfying tb ¼ t. Let t ¼ðt 1, t 2,, t k Þ. Then this system of equations is t 1 ¼ p 1 t k t iþ1 ¼ t i þ p iþ1 t k, i ¼ 1, 2,, k 1 It is easy to show that Convergence of a linear recursive sequence 61 t i ¼ðp 1 þ p 2 þþp i Þt k for i ¼ 1, 2,, k 1. As t is a probability vector, we have P k t i ¼ 1. Also notice that p 1 þ p 2 þþp k ¼ 1. Thus ð p 1 þþp i Þt k ¼ 1 Hence t k t k ¼ ðk i þ 1Þp i ¼ 1 1 P k ðk i þ 1Þp i

12 62 E. G. Tay et al. implying that for i ¼ 1, 2,, k. Therefore t i ¼ p 1 þ p 2 þþp i P k ðk i þ 1Þp i lim n!1 a n ¼ðt 1, t 2,, t k ÞX 1 ¼ P k P i p j a i P k ðk i þ 1Þp i 5. p 1 þ p 2 þþp k 6¼ 1 Let U ¼ p 1 þ p 2 þþp k. In sections 2, 3 and 4, we study the convergence of fa n g 1 n¼1 under the condition that U ¼ 1. In this section, we consider the other two cases: U > 1 and U<1. Lemma 5.1. For any n 5 1, M n 4 U bðn 1Þ=kc M 1 and m n 5 U bðn 1Þ=kc m 1 Proof. By definition, for any n 5 1, a nþk ¼ Xk p i a nþi 1 4 Xk p i M n ¼ UM n Hence a nþkþi 4 UM nþi 4 UM n for any i 5 0. Thus M nþk ¼ maxfa nþkþi 1 : i ¼ 1, 2,, kg 4 UM n Therefore M n 4 U bðn 1Þ=kc M 1. Similarly, we have m n 5 U bðn 1Þ=kc m 1. Theorem 5.1. If U<1, then the sequence fa n g 1 n¼1 always tends to zero. If U > 1 and a i > 0 for i ¼ 1, 2,, k, then the sequence fa n g 1 n¼1 is divergent. Proof. (1) Assume that U<1. By Lemma 5.1, for any integer n 5 1, U bðn 1Þ=kc m 1 4 m n 4 M n 4 U bðn 1Þ=kc M 1 By Lemmas 2.1 and 2.2, the two sequences {M n } and {m n } are convergent. Hence lim M n ¼ lim m n ¼ 0 n!1 n!1 implying that lim a n ¼ 0. n!1 (2) Assume that U > 1 and a i > 0 for i ¼ 1, 2,, k. By Lemma 5.1, m n 5 U bðn 1Þ=kc m 1 Since m 1 > 0 and U >1, {m n } tends to infinity. Hence the sequence fa n g 1 n¼1 is divergent.

13 Convergence of a linear recursive sequence 63 Remarks (i) It is shown in this paper that fa n g 1 n¼1 is not always convergent if p 1 þ p 2 þþp k ¼ 1 and the integers in are not coprime. A problem arises naturally: what initial values a 1, a 2,, a k can be chosen such that fa n g 1 n¼1 is convergent if fa ng 1 n¼1 is not always convergent? (ii) In this paper, we study the convergence of the sequence fa n g 1 n¼1 under the condition that a i 5 0 for all i ¼ 1, 2,, k. Ifp i <0 for some i, the problem appears quite difficult. We have verified a few cases by Excel, and no patterns were found. We think it is worth studying the problem: find conditions on p 1, p 2,, p k, where p 1 þ p 2 þþp k 5 0 and p i <0 for some i, such that fa n g 1 n¼1 is always convergent. Acknowledgment We thank the referees for their valuable comments. References [1] HARARY, F., 1969, Graph Theory (Reading, MA: Addison-Welsey). [2] LIPSCHUTZ, S., 1983, Schaum s Outline of Theory and Problems of Finite Mathematics (New York: McGraw-Hill). [3] DONGSHENG, Z., YEONG, L. T., SENG, L. C., and FWE, Y. S., 2002, Int. J. Math. Educ. Sci. Technol., 33, 123.

Possible numbers of ones in 0 1 matrices with a given rank

Possible numbers of ones in 0 1 matrices with a given rank Linear and Multilinear Algebra, Vol, No, 00, Possible numbers of ones in 0 1 matrices with a given rank QI HU, YAQIN LI and XINGZHI ZHAN* Department of Mathematics, East China Normal University, Shanghai

More information

Combinatorial proofs of Honsberger-type identities

Combinatorial proofs of Honsberger-type identities International Journal of Mathematical Education in Science and Technology, Vol. 39, No. 6, 15 September 2008, 785 792 Combinatorial proofs of Honsberger-type identities A. Plaza* and S. Falco n Department

More information

Sequential dynamical systems over words

Sequential dynamical systems over words Applied Mathematics and Computation 174 (2006) 500 510 www.elsevier.com/locate/amc Sequential dynamical systems over words Luis David Garcia a, Abdul Salam Jarrah b, *, Reinhard Laubenbacher b a Department

More information

Determinants and polynomial root structure

Determinants and polynomial root structure International Journal of Mathematical Education in Science and Technology, Vol 36, No 5, 2005, 469 481 Determinants and polynomial root structure L G DE PILLIS Department of Mathematics, Harvey Mudd College,

More information

THE ALGORITHM TO CALCULATE THE PERIOD MATRIX OF THE CURVE x m þ y n ¼ 1

THE ALGORITHM TO CALCULATE THE PERIOD MATRIX OF THE CURVE x m þ y n ¼ 1 TSUKUA J MATH Vol 26 No (22), 5 37 THE ALGORITHM TO ALULATE THE PERIOD MATRIX OF THE URVE x m þ y n ¼ y Abstract We show how to take a canonical homology basis and a basis of the space of holomorphic -forms

More information

Differentiation matrices in polynomial bases

Differentiation matrices in polynomial bases Math Sci () 5 DOI /s9---x ORIGINAL RESEARCH Differentiation matrices in polynomial bases A Amiraslani Received January 5 / Accepted April / Published online April The Author(s) This article is published

More information

Preservation of local dynamics when applying central difference methods: application to SIR model

Preservation of local dynamics when applying central difference methods: application to SIR model Journal of Difference Equations and Applications, Vol., No. 4, April 2007, 40 Preservation of local dynamics when applying central difference methods application to SIR model LIH-ING W. ROEGER* and ROGER

More information

On Moore Bipartite Digraphs

On Moore Bipartite Digraphs On Moore Bipartite Digraphs M. A. Fiol, 1 * J. Gimbert, 2 J. Gómez, 1 and Y. Wu 3 1 DEPARTMENT DE MATEMÀTICA APLICADA IV TELEMÀTICA UNIVERSITAT POLITÈCNICA DE CATALUNYA JORDI GIRONA 1-3, MÒDUL C3 CAMPUS

More information

THE DOMINATION NUMBER OF GRIDS *

THE DOMINATION NUMBER OF GRIDS * SIAM J. DISCRETE MATH. Vol. 2, No. 3, pp. 1443 143 2011 Society for Industrial and Applied Mathematics THE DOMINATION NUMBER OF GRIDS * DANIEL GONÇALVES, ALEXANDRE PINLOU, MICHAËL RAO, AND STÉPHAN THOMASSÉ

More information

Number of Complete N-ary Subtrees on Galton-Watson Family Trees

Number of Complete N-ary Subtrees on Galton-Watson Family Trees Methodol Comput Appl Probab (2006) 8: 223 233 DOI: 10.1007/s11009-006-8549-6 Number of Complete N-ary Subtrees on Galton-Watson Family Trees George P. Yanev & Ljuben Mutafchiev Received: 5 May 2005 / Revised:

More information

Universally bad integers and the 2-adics

Universally bad integers and the 2-adics Journal of Number Theory 17 (24) 322 334 http://www.elsevier.com/locate/jnt Universally bad integers and the 2-adics S.J. Eigen, a Y. Ito, b and V.S. Prasad c, a Northeastern University, Boston, MA 2115,

More information

Semiconjugate factorizations of higher order linear difference equations in rings

Semiconjugate factorizations of higher order linear difference equations in rings Journal of Difference Equations and Applications ISSN: 1023-6198 (Print) 1563-5120 (Online) Journal homepage: http://www.tandfonline.com/loi/gdea20 Semiconjugate factorizations of higher order linear difference

More information

Meromorphic functions sharing three values

Meromorphic functions sharing three values J. Math. Soc. Japan Vol. 56, No., 2004 Meromorphic functions sharing three values By Xiao-Min Li and Hong-Xun Yi* (Received Feb. 7, 2002) (Revised Aug. 5, 2002) Abstract. In this paper, we prove a result

More information

ON EXCHANGE RINGS WITH BOUNDED INDEX OF NILPOTENCE

ON EXCHANGE RINGS WITH BOUNDED INDEX OF NILPOTENCE COMMUNICATIONS IN ALGEBRA, 29(7), 3089 3098 (2001) ON EXCHANGE RINGS WITH BOUNDED INDEX OF NILPOTENCE Tongso W Department of Applied Mathematics, Shanghai, Jiaotong University, Shanghai 200030, P. R. China

More information

Stability analysis for a class of Takagi Sugeno fuzzy control systems with PID controllers

Stability analysis for a class of Takagi Sugeno fuzzy control systems with PID controllers International Journal of Approximate Reasoning 46 (2007) 109 119 www.elsevier.com/locate/ijar Stability analysis for a class of Takagi Sugeno fuzzy control systems with PID controllers Le Hung Lan * Institute

More information

On memory gradient method with trust region for unconstrained optimization*

On memory gradient method with trust region for unconstrained optimization* Numerical Algorithms (2006) 41: 173 196 DOI: 10.1007/s11075-005-9008-0 * Springer 2006 On memory gradient method with trust region for unconstrained optimization* Zhen-Jun Shi a,b and Jie Shen b a College

More information

Appendix A Conventions

Appendix A Conventions Appendix A Conventions We use natural units h ¼ 1; c ¼ 1; 0 ¼ 1; ða:1þ where h denotes Planck s constant, c the vacuum speed of light and 0 the permittivity of vacuum. The electromagnetic fine-structure

More information

An interior point type QP-free algorithm with superlinear convergence for inequality constrained optimization

An interior point type QP-free algorithm with superlinear convergence for inequality constrained optimization Applied Mathematical Modelling 31 (2007) 1201 1212 www.elsevier.com/locate/apm An interior point type QP-free algorithm with superlinear convergence for inequality constrained optimization Zhibin Zhu *

More information

Semiregular Modules with Respect to a Fully Invariant Submodule #

Semiregular Modules with Respect to a Fully Invariant Submodule # COMMUNICATIONS IN ALGEBRA Õ Vol. 32, No. 11, pp. 4285 4301, 2004 Semiregular Modules with Respect to a Fully Invariant Submodule # Mustafa Alkan 1, * and A. Çiğdem Özcan 2 1 Department of Mathematics,

More information

Asymptotic stability for difference equations with decreasing arguments

Asymptotic stability for difference equations with decreasing arguments Journal of Difference Equations and Applications ISSN: 1023-6198 (Print) 1563-5120 (Online) Journal homepage: http://www.tandfonline.com/loi/gdea20 Asymptotic stability for difference equations with decreasing

More information

Optimal algorithm for minimizing production cycle time of a printed circuit board assembly line

Optimal algorithm for minimizing production cycle time of a printed circuit board assembly line int. j. prod. res., 1 december 2004, vol. 42, no. 23, 5031 504 Optimal algorithm for minimizing production cycle time of a printed circuit board assembly line D. M. KODEK* and M. KRISPER The problem of

More information

Linear Algebra and its Applications

Linear Algebra and its Applications Linear Algebra and its Applications 432 2010 661 669 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: wwwelseviercom/locate/laa On the characteristic and

More information

Constructions of difference covering arrays $

Constructions of difference covering arrays $ Journal of Combinatorial Theory, Series A 104 (00) 7 9 Constructions of difference covering arrays $ Jianxing Yin Department of Mathematics, Suzhou University, Suzhou 15006, China Received 19 February

More information

8. Limit Laws. lim(f g)(x) = lim f(x) lim g(x), (x) = lim x a f(x) g lim x a g(x)

8. Limit Laws. lim(f g)(x) = lim f(x) lim g(x), (x) = lim x a f(x) g lim x a g(x) 8. Limit Laws 8.1. Basic Limit Laws. If f and g are two functions and we know the it of each of them at a given point a, then we can easily compute the it at a of their sum, difference, product, constant

More information

Convergence, Periodicity and Bifurcations for the Two-parameter Absolute-Difference Equation

Convergence, Periodicity and Bifurcations for the Two-parameter Absolute-Difference Equation Journal of Difference Equations and pplications ISSN: 123-6198 (Print) 1563-512 (Online) Journal homepage: http://www.tandfonline.com/loi/gdea2 Convergence, Periodicity and Bifurcations for the Two-parameter

More information

Commun Nonlinear Sci Numer Simulat

Commun Nonlinear Sci Numer Simulat Commun Nonlinear Sci Numer Simulat 14 (2009) 4041 4056 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Symbolic computation

More information

CONSISTENT ESTIMATION THROUGH WEIGHTED HARMONIC MEAN OF INCONSISTENT ESTIMATORS IN REPLICATED MEASUREMENT ERROR MODELS

CONSISTENT ESTIMATION THROUGH WEIGHTED HARMONIC MEAN OF INCONSISTENT ESTIMATORS IN REPLICATED MEASUREMENT ERROR MODELS ECONOMETRIC REVIEWS, 20(4), 507 50 (200) CONSISTENT ESTIMATION THROUGH WEIGHTED HARMONIC MEAN OF INCONSISTENT ESTIMATORS IN RELICATED MEASUREMENT ERROR MODELS Shalabh Department of Statistics, anjab University,

More information

Signal Processing 93 (2013) Contents lists available at SciVerse ScienceDirect. Signal Processing

Signal Processing 93 (2013) Contents lists available at SciVerse ScienceDirect. Signal Processing Signal Processing 93 (03) 7 730 Contents lists available at SciVerse ScienceDirect Signal Processing journal homepage: www.elsevier.com/locate/sigpro Lainiotis filter, golden section and Fibonacci sequence

More information

Link invariant and G 2

Link invariant and G 2 Hiroshima Math. J. 47 (2017), 19 41 Link invariant and G 2 web space Takuro Sakamoto and Yasuyoshi Yonezawa (Received September 1, 2015) (Revised October 24, 2016) Abstract. In this paper, we reconstruct

More information

MATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS

MATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS MATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS 1. We have one theorem whose conclusion says an alternating series converges. We have another theorem whose conclusion says an alternating series diverges.

More information

NONSTATIONARY DISCRETE CHOICE: A CORRIGENDUM AND ADDENDUM. PETER C. B. PHILIPS, SAINAN JIN and LING HU COWLES FOUNDATION PAPER NO.

NONSTATIONARY DISCRETE CHOICE: A CORRIGENDUM AND ADDENDUM. PETER C. B. PHILIPS, SAINAN JIN and LING HU COWLES FOUNDATION PAPER NO. NONSTATIONARY DISCRETE CHOICE: A CORRIGENDUM AND ADDENDUM BY PETER C. B. PHILIPS, SAINAN JIN and LING HU COWLES FOUNDATION PAPER NO. 24 COWLES FOUNDATION FOR RESEARCH IN ECONOMICS YALE UNIVERSITY Box 2828

More information

Numerical solution of hyperbolic heat conduction in thin surface layers

Numerical solution of hyperbolic heat conduction in thin surface layers International Journal of Heat and Mass Transfer 50 (007) 9 www.elsevier.com/locate/ijhmt Numerical solution of hyperbolic heat conduction in thin surface layers Tzer-Ming Chen * Department of Vehicle Engineering,

More information

Ch. 7.6 Squares, Squaring & Parabolas

Ch. 7.6 Squares, Squaring & Parabolas Ch. 7.6 Squares, Squaring & Parabolas Learning Intentions: Learn about the squaring & square root function. Graph parabolas. Compare the squaring function with other functions. Relate the squaring function

More information

Algorithmic analysis of the Geo/Geo/c retrial queue

Algorithmic analysis of the Geo/Geo/c retrial queue European Journal of Operational Research xxx (2007) xxx xxx www.elsevier.com/locate/ejor Algorithmic analysis of the Geo/Geo/c retrial queue Jesus R. Artalejo a, *, Antonis Economou b, Antonio Gómez-Corral

More information

Higher-order sliding modes, differentiation and output-feedback control

Higher-order sliding modes, differentiation and output-feedback control INT. J. CONTROL, 2003, VOL. 76, NOS 9/10, 924 941 Higher-order sliding modes, differentiation and output-feedback control ARIE LEVANT{ Being a motion on a discontinuity set of a dynamic system, sliding

More information

MATH FINAL EXAM REVIEW HINTS

MATH FINAL EXAM REVIEW HINTS MATH 109 - FINAL EXAM REVIEW HINTS Answer: Answer: 1. Cardinality (1) Let a < b be two real numbers and define f : (0, 1) (a, b) by f(t) = (1 t)a + tb. (a) Prove that f is a bijection. (b) Prove that any

More information

Infinite Continued Fractions

Infinite Continued Fractions Infinite Continued Fractions 8-5-200 The value of an infinite continued fraction [a 0 ; a, a 2, ] is lim c k, where c k is the k-th convergent k If [a 0 ; a, a 2, ] is an infinite continued fraction with

More information

Dimension of the Mesh Algebra of a Finite Auslander Reiten Quiver

Dimension of the Mesh Algebra of a Finite Auslander Reiten Quiver COMMUNICATIONS IN ALGEBRA Õ Vol. 31, No. 5, pp. 2207 2217, 2003 Dimension of the Mesh Algebra of a Finite Auslander Reiten Quiver Ragnar-Olaf Buchweitz 1, * and Shiping Liu 2 1 Department of Mathematics,

More information

Chapter 2 Linear Systems

Chapter 2 Linear Systems Chapter 2 Linear Systems This chapter deals with linear systems of ordinary differential equations (ODEs, both homogeneous and nonhomogeneous equations Linear systems are etremely useful for analyzing

More information

A constructive converse Lyapunov theorem on asymptotic stability for nonlinear autonomous ordinary differential equations

A constructive converse Lyapunov theorem on asymptotic stability for nonlinear autonomous ordinary differential equations Dynamical Systems, Vol. 20, No. 3, September 2005, 281 299 A constructive converse Lyapunov theorem on asymptotic stability for nonlinear autonomous ordinary differential equations S. F. HAFSTEIN* University

More information

Edge colored complete bipartite graphs with trivial automorphism groups

Edge colored complete bipartite graphs with trivial automorphism groups Edge colored complete bipartite graphs with trivial automorphism groups Michael J. Fisher Garth Isaak Abstract We determine the values of s and t for which there is a coloring of the edges of the complete

More information

Minimal almost convexity

Minimal almost convexity J. Group Theory 8 (2005), 239 266 Journal of Group Theory ( de Gruyter 2005 Minimal almost convexity Murray Elder and Susan Hermiller* (Communicated by D. J. S. Robinson) Abstract. In this article we show

More information

Applied Mathematics and Computation

Applied Mathematics and Computation Applied Mathematics and Computation 245 (2014) 86 107 Contents lists available at ScienceDirect Applied Mathematics and Computation journal homepage: www.elsevier.com/locate/amc An analysis of a new family

More information

RESEARCH ARTICLE. An extension of the polytope of doubly stochastic matrices

RESEARCH ARTICLE. An extension of the polytope of doubly stochastic matrices Linear and Multilinear Algebra Vol. 00, No. 00, Month 200x, 1 15 RESEARCH ARTICLE An extension of the polytope of doubly stochastic matrices Richard A. Brualdi a and Geir Dahl b a Department of Mathematics,

More information

Fuzzy age-dependent replacement policy and SPSA algorithm based-on fuzzy simulation

Fuzzy age-dependent replacement policy and SPSA algorithm based-on fuzzy simulation Available online at wwwsciencedirectcom Information Sciences 178 (2008) 573 583 wwwelseviercom/locate/ins Fuzzy age-dependent replacement policy and SPSA algorithm based-on fuzzy simulation Jiashun Zhang,

More information

A peak factor for non-gaussian response analysis of wind turbine tower

A peak factor for non-gaussian response analysis of wind turbine tower Journal of Wind Engineering and Industrial Aerodynamics 96 (28) 227 2227 www.elsevier.com/locate/jweia A peak factor for non-gaussian response analysis of wind turbine tower Luong Van Binh a,, Takeshi

More information

Daniel Onofrei a b & Bogdan Vernescu a c a Department of Mathematics, University of Utah, JWB, Building,

Daniel Onofrei a b & Bogdan Vernescu a c a Department of Mathematics, University of Utah, JWB, Building, This article was downloaded by: [University of Houston], [Daniel Onofrei] On: 08 March 2012, At: 10:29 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954

More information

The Product of Like-Indexed Terms in Binary Recurrences

The Product of Like-Indexed Terms in Binary Recurrences Journal of Number Theory 96, 152 173 (2002) doi:10.1006/jnth.2002.2794 The Product of Like-Indexed Terms in Binary Recurrences F. Luca 1 Instituto de Matemáticas UNAM, Campus Morelia, Ap. Postal 61-3 (Xangari),

More information

SAMPLE SIZE AND OPTIMAL DESIGNS IN STRATIFIED COMPARATIVE TRIALS TO ESTABLISH THE EQUIVALENCE OF TREATMENT EFFECTS AMONG TWO ETHNIC GROUPS

SAMPLE SIZE AND OPTIMAL DESIGNS IN STRATIFIED COMPARATIVE TRIALS TO ESTABLISH THE EQUIVALENCE OF TREATMENT EFFECTS AMONG TWO ETHNIC GROUPS MARCEL DEKKER, INC. 70 MADISON AVENUE NEW YORK, NY 006 JOURNAL OF BIOPHARMACEUTICAL STATISTICS Vol., No. 4, pp. 553 566, 00 SAMPLE SIZE AND OPTIMAL DESIGNS IN STRATIFIED COMPARATIVE TRIALS TO ESTABLISH

More information

Weak (a,d)-skew Armendariz ideals

Weak (a,d)-skew Armendariz ideals Matematical Sciences (28) 2:235 242 tts://doiorg/7/s46-8-263-3 (2345678(),-volV)(2345678(),-volV) ORIGINL PPER Weak (a,d)-skew rmendariz ideals R Majdabadi Faraani F Heydari H Tavallaee M Magasedi Received:

More information

Primitive Matrices with Combinatorial Properties

Primitive Matrices with Combinatorial Properties Southern Illinois University Carbondale OpenSIUC Research Papers Graduate School Fall 11-6-2012 Primitive Matrices with Combinatorial Properties Abdulkarem Alhuraiji al_rqai@yahoo.com Follow this and additional

More information

Finite difference/spectral approximations for the time-fractional diffusion equation q

Finite difference/spectral approximations for the time-fractional diffusion equation q Journal of Computational Physics 225 (2007) 533 552 www.elsevier.com/locate/jcp Finite difference/spectral approximations for the time-fractional diffusion equation q Yumin Lin, Chuanju Xu * School of

More information

Mathematics MFP2 (JUN15MFP201) General Certificate of Education Advanced Level Examination June Unit Further Pure TOTAL

Mathematics MFP2 (JUN15MFP201) General Certificate of Education Advanced Level Examination June Unit Further Pure TOTAL Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Further Pure 2 Tuesday 16 June 2015 General Certificate of Education Advanced

More information

A QBD approach to evolutionary game theory

A QBD approach to evolutionary game theory Applied Mathematical Modelling (00) 91 9 wwwelseviercom/locate/apm A QBD approach to evolutionary game theory Lotfi Tadj a, *, Abderezak Touzene b a Department of Statistics and Operations Research, College

More information

THE RIEMANN INTEGRAL USING ORDERED OPEN COVERINGS

THE RIEMANN INTEGRAL USING ORDERED OPEN COVERINGS ROCKY MOUNTAIN JOURNAL OF MATHEMATICS THE RIEMANN INTEGRAL USING ORDERED OPEN COVERINGS ZHAO DONGSHENG AND LEE PENG YEE ABSTRACT. We define the Riemann integral for bounded functions defined on a general

More information

UNDERSTANDING RULER AND COMPASS CONSTRUCTIONS WITH FIELD THEORY

UNDERSTANDING RULER AND COMPASS CONSTRUCTIONS WITH FIELD THEORY UNDERSTANDING RULER AND COMPASS CONSTRUCTIONS WITH FIELD THEORY ISAAC M. DAVIS Abstract. By associating a subfield of R to a set of points P 0 R 2, geometric properties of ruler and compass constructions

More information

A novel variable structure control scheme for chaotic synchronization

A novel variable structure control scheme for chaotic synchronization Chaos, Solitons and Fractals 18 (2003) 275 287 www.elsevier.com/locate/chaos A novel variable structure control scheme for chaotic synchronization Chun-Chieh Wang b,c, Juhng-Perng Su a, * a Department

More information

FACTOR GRAPH OF NON-COMMUTATIVE RING

FACTOR GRAPH OF NON-COMMUTATIVE RING International Journal of Advanced Research in Engineering and Technology (IJARET) Volume 9, Issue 6, November - December 2018, pp. 178 183, Article ID: IJARET_09_06_019 Available online at http://www.iaeme.com/ijaret/issues.asp?jtype=ijaret&vtype=9&itype=6

More information

Sequences, Series, and the Binomial Formula

Sequences, Series, and the Binomial Formula CHAPTER Sequences, Series, nd the Binomil Formul. SEQUENCES. ; ; ; ; 6 ; 6 6. ðþ ; ðþ ; ð Þ 6; ðþ ; ðþ 6; 6 ð6þ. ðþ ; ðþ : ðþ ; ðþ ; ðþ ; 6 ðþ 6 6 6. ; ; ; ; ; 6 6 6. 0 ; ; ; 8 ; 6 8 ; 6. 0; ; 6 ; ; 6

More information

A null space method for solving system of equations q

A null space method for solving system of equations q Applied Mathematics and Computation 149 (004) 15 6 www.elsevier.com/locate/amc A null space method for solving system of equations q Pu-yan Nie 1 Department of Mathematics, Jinan University, Guangzhou

More information

Optimization of machining parameters of Wire-EDM based on Grey relational and statistical analyses

Optimization of machining parameters of Wire-EDM based on Grey relational and statistical analyses int. j. prod. res., 2003, vol. 41, no. 8, 1707 1720 Optimization of machining parameters of Wire-EDM based on Grey relational and statistical analyses J. T. HUANG{* and Y. S. LIAO{ Grey relational analyses

More information

An Explicit Formula for the Discrete Power Function Associated with Circle Patterns of Schramm Type

An Explicit Formula for the Discrete Power Function Associated with Circle Patterns of Schramm Type Funkcialaj Ekvacioj, 57 (2014) 1 41 An Explicit Formula for the Discrete Power Function Associated with Circle Patterns of Schramm Type By Hisashi Ando1, Mike Hay2, Kenji Kajiwara1 and Tetsu Masuda3 (Kyushu

More information

A note on solution sensitivity for Karush Kuhn Tucker systems*

A note on solution sensitivity for Karush Kuhn Tucker systems* Math Meth Oper Res (25) 6: 347 363 DOI.7/s86449 A note on solution sensitivity for Karush Kuhn ucker systems* A. F. Izmailov, M. V. Solodov 2 Faculty of Computational Mathematics and Cybernetics, Department

More information

How do our representations change if we select another basis?

How do our representations change if we select another basis? CHAPTER 6 Linear Mappings and Matrices 99 THEOREM 6.: For any linear operators F; G AðV Þ, 6. Change of Basis mðg FÞ ¼ mðgþmðfþ or ½G FŠ ¼ ½GŠ½FŠ (Here G F denotes the composition of the maps G and F.)

More information

Some spectral inequalities for triangle-free regular graphs

Some spectral inequalities for triangle-free regular graphs Filomat 7:8 (13), 1561 1567 DOI 198/FIL138561K Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://wwwpmfniacrs/filomat Some spectral inequalities for triangle-free

More information

The Singapore Copyright Act applies to the use of this document.

The Singapore Copyright Act applies to the use of this document. Title On graphs whose low polynomials have real roots only Author(s) Fengming Dong Source Electronic Journal of Combinatorics, 25(3): P3.26 Published by Electronic Journal of Combinatorics This document

More information

Mathematics MPC1 (JUN15MPC101) General Certificate of Education Advanced Subsidiary Examination June Unit Pure Core TOTAL

Mathematics MPC1 (JUN15MPC101) General Certificate of Education Advanced Subsidiary Examination June Unit Pure Core TOTAL Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Pure Core 1 Wednesday 13 May 2015 General Certificate of Education Advanced

More information

New design of estimators using covariance information with uncertain observations in linear discrete-time systems

New design of estimators using covariance information with uncertain observations in linear discrete-time systems Applied Mathematics and Computation 135 (2003) 429 441 www.elsevier.com/locate/amc New design of estimators using covariance information with uncertain observations in linear discrete-time systems Seiichi

More information

The Number of Independent Sets in a Regular Graph

The Number of Independent Sets in a Regular Graph Combinatorics, Probability and Computing (2010) 19, 315 320. c Cambridge University Press 2009 doi:10.1017/s0963548309990538 The Number of Independent Sets in a Regular Graph YUFEI ZHAO Department of Mathematics,

More information

THE NUMBER OF INDEPENDENT DOMINATING SETS OF LABELED TREES. Changwoo Lee. 1. Introduction

THE NUMBER OF INDEPENDENT DOMINATING SETS OF LABELED TREES. Changwoo Lee. 1. Introduction Commun. Korean Math. Soc. 18 (2003), No. 1, pp. 181 192 THE NUMBER OF INDEPENDENT DOMINATING SETS OF LABELED TREES Changwoo Lee Abstract. We count the numbers of independent dominating sets of rooted labeled

More information

Commun Nonlinear Sci Numer Simulat

Commun Nonlinear Sci Numer Simulat Commun Nonlinear Sci Numer Simulat xxx ) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: wwwelseviercom/locate/cnsns A note on the use of Adomian

More information

Classification of a subclass of low-dimensional complex filiform Leibniz algebras

Classification of a subclass of low-dimensional complex filiform Leibniz algebras Linear Multilinear lgebra ISSN: 008-087 (Print) 56-59 (Online) Journal homepage: http://www.tfonline.com/loi/glma20 Classification of a subclass of low-dimensional complex filiform Leibniz algebras I.S.

More information

Global attractivity in a rational delay difference equation with quadratic terms

Global attractivity in a rational delay difference equation with quadratic terms Journal of Difference Equations and Applications ISSN: 1023-6198 (Print) 1563-5120 (Online) Journal homepage: http://www.tandfonline.com/loi/gdea20 Global attractivity in a rational delay difference equation

More information

Elastic solutions for stresses in a transversely isotropic half-space subjected to three-dimensional buried parabolic rectangular loads

Elastic solutions for stresses in a transversely isotropic half-space subjected to three-dimensional buried parabolic rectangular loads INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal. Meth. Geomech., 2002; 26:1449 1476 (DOI: 10.1002/nag.253) Elastic solutions for stresses in a transversely

More information

To get horizontal and slant asymptotes algebraically we need to know about end behaviour for rational functions.

To get horizontal and slant asymptotes algebraically we need to know about end behaviour for rational functions. Concepts: Horizontal Asymptotes, Vertical Asymptotes, Slant (Oblique) Asymptotes, Transforming Reciprocal Function, Sketching Rational Functions, Solving Inequalities using Sign Charts. Rational Function

More information

A four-class association scheme derived from a hyperbolic quadric in PG(3, q)

A four-class association scheme derived from a hyperbolic quadric in PG(3, q) Adv. Geom. 4 (24), 5 7 Advances in Geometry ( de Gruyter 24 A four-class association scheme derived from a hyperbolic quadric in PG(3, q) (Communicated by E. Bannai) Abstract. We prove the existence of

More information

Progression-free sets in finite abelian groups

Progression-free sets in finite abelian groups Journal of Number Theory 104 (2004) 162 169 http://www.elsevier.com/locate/jnt Progression-free sets in finite abelian groups Vsevolod F. Lev Department of Mathematics, University of Haifa at Oranim, Tivon

More information

1 Lecture 8: Interpolating polynomials.

1 Lecture 8: Interpolating polynomials. 1 Lecture 8: Interpolating polynomials. 1.1 Horner s method Before turning to the main idea of this part of the course, we consider how to evaluate a polynomial. Recall that a polynomial is an expression

More information

On (k, d)-multiplicatively indexable graphs

On (k, d)-multiplicatively indexable graphs Chapter 3 On (k, d)-multiplicatively indexable graphs A (p, q)-graph G is said to be a (k,d)-multiplicatively indexable graph if there exist an injection f : V (G) N such that the induced function f :

More information

Math 113 (Calculus 2) Exam 4

Math 113 (Calculus 2) Exam 4 Math 3 (Calculus ) Exam 4 November 0 November, 009 Sections 0, 3 7 Name Student ID Section Instructor In some cases a series may be seen to converge or diverge for more than one reason. For such problems

More information

Multicriteria decision-making method using the correlation coefficient under single-valued neutrosophic environment

Multicriteria decision-making method using the correlation coefficient under single-valued neutrosophic environment International Journal of General Systems, 2013 Vol. 42, No. 4, 386 394, http://dx.doi.org/10.1080/03081079.2012.761609 Multicriteria decision-making method using the correlation coefficient under single-valued

More information

NAG Library Routine Document D02JAF.1

NAG Library Routine Document D02JAF.1 D02 Ordinary Differential Equations NAG Library Routine Document Note: before using this routine, please read the Users Note for your implementation to check the interpretation of bold italicised terms

More information

10. Smooth Varieties. 82 Andreas Gathmann

10. Smooth Varieties. 82 Andreas Gathmann 82 Andreas Gathmann 10. Smooth Varieties Let a be a point on a variety X. In the last chapter we have introduced the tangent cone C a X as a way to study X locally around a (see Construction 9.20). It

More information

Primitive Digraphs with Smallest Large Exponent

Primitive Digraphs with Smallest Large Exponent Primitive Digraphs with Smallest Large Exponent by Shahla Nasserasr B.Sc., University of Tabriz, Iran 1999 A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of MASTER OF SCIENCE

More information

Analysis-3 lecture schemes

Analysis-3 lecture schemes Analysis-3 lecture schemes (with Homeworks) 1 Csörgő István November, 2015 1 A jegyzet az ELTE Informatikai Kar 2015. évi Jegyzetpályázatának támogatásával készült Contents 1. Lesson 1 4 1.1. The Space

More information

Parameter Estimation for Partially Complete Time and Type of Failure Data

Parameter Estimation for Partially Complete Time and Type of Failure Data Biometrical Journal 46 (004), 65 79 DOI 0.00/bimj.00004 arameter Estimation for artially Complete Time and Type of Failure Data Debasis Kundu Department of Mathematics, Indian Institute of Technology Kanpur,

More information

Graphs and Their Applications (7)

Graphs and Their Applications (7) 11111111 I I I Graphs and Their Applications (7) by K.M. Koh* Department of Mathematics National University of Singapore, Singapore 117543 F.M. Dong and E.G. Tay Mathematics and Mathematics Education National

More information

Proximal-Based Pre-correction Decomposition Methods for Structured Convex Minimization Problems

Proximal-Based Pre-correction Decomposition Methods for Structured Convex Minimization Problems J. Oper. Res. Soc. China (2014) 2:223 235 DOI 10.1007/s40305-014-0042-2 Proximal-Based Pre-correction Decomposition Methods for Structured Convex Minimization Problems Yuan-Yuan Huang San-Yang Liu Received:

More information

CHAPTER 1. Relations. 1. Relations and Their Properties. Discussion

CHAPTER 1. Relations. 1. Relations and Their Properties. Discussion CHAPTER 1 Relations 1. Relations and Their Properties 1.1. Definition of a Relation. Definition 1.1.1. A binary relation from a set A to a set B is a subset R A B. If (a, b) R we say a is Related to b

More information

Chapter 2 Analysis Methods

Chapter 2 Analysis Methods Chapter Analysis Methods. Nodal Analysis Problem.. Two current sources with equal internal resistances feed a load as shown in Fig... I a ¼ 00 A; I b ¼ 00 A; R ¼ 00 X; R L ¼ 00 X: (a) Find the current

More information

Asymptotic efficiency and small sample power of a locally most powerful linear rank test for the log-logistic distribution

Asymptotic efficiency and small sample power of a locally most powerful linear rank test for the log-logistic distribution Math Sci (2014) 8:109 115 DOI 10.1007/s40096-014-0135-4 ORIGINAL RESEARCH Asymptotic efficiency and small sample power of a locally most powerful linear rank test for the log-logistic distribution Hidetoshi

More information

CHARACTERIZING MANIFOLDS MODELED ON CERTAIN DENSE SUBSPACES OF NON-SEPARABLE HILBERT SPACES

CHARACTERIZING MANIFOLDS MODELED ON CERTAIN DENSE SUBSPACES OF NON-SEPARABLE HILBERT SPACES TSUKUBA J. MATH. Vol. 7 No. 1 (003), 143 159 CHARACTERIZING MANIFOLDS MODELED ON CERTAIN DENSE SUBSPACES OF NON-SEPARABLE HILBERT SPACES By Abstract. For an infinite set G, let l f ðgþ be the linear span

More information

Solutions to generalized Sylvester matrix equation by Schur decomposition

Solutions to generalized Sylvester matrix equation by Schur decomposition International Journal of Systems Science Vol 8, No, May 007, 9 7 Solutions to generalized Sylvester matrix equation by Schur decomposition BIN ZHOU* and GUANG-REN DUAN Center for Control Systems and Guidance

More information

5 Flows and cuts in digraphs

5 Flows and cuts in digraphs 5 Flows and cuts in digraphs Recall that a digraph or network is a pair G = (V, E) where V is a set and E is a multiset of ordered pairs of elements of V, which we refer to as arcs. Note that two vertices

More information

MATH 205 HOMEWORK #3 OFFICIAL SOLUTION. Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. (a) F = R, V = R 3,

MATH 205 HOMEWORK #3 OFFICIAL SOLUTION. Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. (a) F = R, V = R 3, MATH 205 HOMEWORK #3 OFFICIAL SOLUTION Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. a F = R, V = R 3, b F = R or C, V = F 2, T = T = 9 4 4 8 3 4 16 8 7 0 1

More information

A Note on the 2 F 1 Hypergeometric Function

A Note on the 2 F 1 Hypergeometric Function A Note on the F 1 Hypergeometric Function Armen Bagdasaryan Institution of the Russian Academy of Sciences, V.A. Trapeznikov Institute for Control Sciences 65 Profsoyuznaya, 117997, Moscow, Russia E-mail:

More information

Get started [Hawkes Learning] with this system. Common final exam, independently administered, group graded, grades reported.

Get started [Hawkes Learning] with this system. Common final exam, independently administered, group graded, grades reported. Course Information Math 095 Elementary Algebra Placement No placement necessary Course Description Learning Outcomes Elementary algebraic topics for students whose mathematical background or placement

More information

On the Hyers-Ulam Stability Problem for Quadratic Multi-dimensional Mappings on the Gaussian Plane

On the Hyers-Ulam Stability Problem for Quadratic Multi-dimensional Mappings on the Gaussian Plane Southeast Asian Bulletin of Mathematics (00) 6: 483 50 Southeast Asian Bulletin of Mathematics : Springer-Verlag 00 On the Hyers-Ulam Stability Problem f Quadratic Multi-dimensional Mappings on the Gaussian

More information

Communicated by Alireza Abdollahi. 1. Introduction. For a set S V, the open neighborhood of S is N(S) = v S

Communicated by Alireza Abdollahi. 1. Introduction. For a set S V, the open neighborhood of S is N(S) = v S Transactions on Combinatorics ISSN print: 2251-8657, ISSN on-line: 2251-8665 Vol. 01 No. 2 2012, pp. 49-57. c 2012 University of Isfahan www.combinatorics.ir www.ui.ac.ir ON THE VALUES OF INDEPENDENCE

More information

Applied Mathematics Letters

Applied Mathematics Letters Applied Mathematics Letters 25 (2012) 974 979 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml On dual vector equilibrium problems

More information