• 検索結果がありません。

Dynamics of a Higher-Order Nonlinear Difference Equation

N/A
N/A
Protected

Academic year: 2022

シェア "Dynamics of a Higher-Order Nonlinear Difference Equation"

Copied!
15
0
0

読み込み中.... (全文を見る)

全文

(1)

Volume 2010, Article ID 534947,15pages doi:10.1155/2010/534947

Research Article

Dynamics of a Higher-Order Nonlinear Difference Equation

Guo-Mei Tang,

1

Lin-Xia Hu,

2

and Xiu-Mei Jia

3

1School of Mathematics and Computer Science, Northwest University for Nationalities, Lanzhou, Gansu 730030, China

2Department of Mathematics, Tianshui Normal University, Tianshui, Gansu 741001, China

3Department of Mathematics, Hexi University, Zhangye, Gansu 734000, China

Correspondence should be addressed to Guo-Mei Tang,[email protected] Received 20 April 2010; Accepted 18 July 2010

Academic Editor: Leonid Berezansky

Copyrightq2010 Guo-Mei Tang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We consider the higher-order nonlinear difference equationxn1 αxn/ABxnxn−k, n 0,1, . . ., where parameters are positive real numbers and initial conditionsx−k, . . . , x0 are nonnegative real numbers,k ≥ 2. We investigate the periodic character, the invariant intervals, and the global asymptotic stability of all positive solutions of the abovementioned equation. We show that the unique equilibrium of the equation is globally asymptotically stable under certain conditions.

1. Introduction and Preliminaries

In this paper, we will investigate the global behavior of solutions of the following nonlinear difference equation:

xn1 αxn

ABxnxn−k, n0,1, . . ., 1.1 where parameters are positive real numbers and initial conditionsx−k, . . . , x0are nonnegative real numbers,k≥2.

In 2003, the authors in1considered the difference equation

xn1 αβxn

ABxnCxn−1, n0,1, . . ., 1.2

(2)

with nonnegative parameters α, β, A, B, C and nonnegative initial conditions x−1, x0. They obtained some global asymptotic stability results for the solutions of1.2. For1.2, we can also see2–4.

For the global behavior of solutions of some related equations, see 5–13. Other related results can be found in14–20. For the sake of convenience, we recall some definitions and theorems which will be useful in the sequel.

Definition 1.1. LetIbe some interval of real numbers, and let

f:Im1−→I 1.3

be a continuously differential function. Then for every set of initial conditionsy−k,. . .,y−1,y0I, the difference equation

yn1f

yn, yn−1, . . . , yn−k

, n0,1,2, . . . , 1.4

has a unique solution{yn}n−k.

A pointyis called an equilibrium point of1.4if yf

y, y, . . . , y

. 1.5

That is,

yny, forn≥0, 1.6

is a solution of1.4, or equivalentlyyis a fixed point off.

Definition 1.2. Letybe an equilibrium point of1.4. Then the following are considered.

iThe equilibriumyis called locally stableor stableif, for everyε > 0,there exists δ >0 such that, for ally−k, . . . , y−1, y0Iwithi0

i−k|yiy|< δ,we have|yny|< ε for allnk.

iiThe equilibriumyof1.4is called locally asymptotically stableasymptotic stable if it is locally stable and if there existsγ >0 such that, for ally−k, . . . , y−1, y0Iwith i0

i−k|yiy|< γ,we have limn→ ∞yny.

iiiThe equilibriumyof1.4is called a global attractor if, for everyy−k, . . . , y−1, y0I, we have limn→ ∞yny.

ivThe equilibriumyof1.4is globally asymptotically stable if it is locally stable and is a global attractor.

vThe equilibriumyof1.4is called unstable if it is not stable.

viThe equilibriumyof1.4is called a source, or a repeller, if there existsr >0 such that, for ally−k, . . . , y−1, y0Iwithi0

i−k|yiy|< γ,there existsN ≥1 such that

|yNy| ≥r.

(3)

An intervalJIis called an invariant interval for1.4if

y−k, . . . , y0JynJ ∀n >0. 1.7

That is, every solution of1.4with initial conditions inJremains inJ.

The linearized equation associated with1.4about the equilibriumyis

yn1k

io

∂f

∂ui

y, . . . , y

yn−i, n0,1,2, . . .. 1.8

Its characteristic equation is

λk1k

io

∂f

∂ui

y, . . . , y

λk−i. 1.9

Theorem 1.3see10. Assume thatfis aC1function, and letybe an equilibrium of 1.4. Then the following statements are true.

iIf all the roots of 1.9lie in the open unit disk|λ|<1, then the equilibriumyof 1.4is asymptotically stable.

iiIf at least one root of 1.9has absolute value greater than one, then the equilibriumyof 1.4is unstable.

Theorem 1.4see10. Assume thatP, QRandk∈ {1,2, . . .}. Then

|P||Q|<1 1.10

is a sufficient condition for the asymptotic stability of the difference equation

yn1PynQyn−k, n0,1, . . . . 1.11 Lemma 1.5see8. Consider the difference equation

yn1f

yn, yn−k

, n0,1, . . ., 1.12

wherek∈ {1,2, . . .}.LetI a, bbe some interval of real numbers and assume that

f :a, b×a, b−→a, b 1.13

is a continuous function satisfying the following properties.

afx, yis nondecreasing inx∈a, bfor eachy∈a, b, andfx, yis nonincreasing in y∈a, bfor eachx∈a, b.

(4)

bIfm, M∈a, b×a, bis a solution of the system

fm, M m, fM, m M, 1.14

thenmM.

Then1.12has a unique equilibriumx∈a, band every solution of 1.12converges tox.

Lemma 1.6see8. Consider the difference equation yn1f

yn, yn−k

, n0,1, . . ., 1.15

wherek∈ {1,2, . . .}.LetI a, bbe some interval of real numbers and assume that

f :a, b×a, b−→a, b 1.16

is a continuous function satisfying the following properties.

afx, yis nonincreasing in each of its arguments.

bIfm,M∈a, b×a, bis a solution of the system

mfM, M, Mfm, m, 1.17

thenmM.

Then1.15has a unique equilibriumy∈a, band every solution of 1.15converges toy.

2. Local Stability and Period-Two Solutions

The equilibria of1.1are the solutions of the equation

x αx

ABxx. 2.1

So1.1possesses the unique positive equilibrium

x 1−A

1−A24αB1

2B1 . 2.2

The linearized equation associated with1.1about the positive equilibrium is zn1 Bx−1

ABxxzn x

ABxxzn−k0. 2.3

The next result follows fromTheorem 1.4.

(5)

Theorem 2.1. Assume that

eitherA1 orA <1, B−11−A24B2α >0. 2.4

Then the positive equilibriumxof 1.1is locally asymptotically stable.

Theorem 2.2. Equation1.1has no nonnegative prime period-two solution.

Proof. Assume for the sake of contradiction that there exist distinct nonnegative real numbers φandψsuch that

. . . , φ, ψ, φ, ψ, . . . 2.5

is a prime period-two solution of1.1.

aAssume thatkis odd. Thenxn1xn−kandφ, ψsatisfy the following system:

φ αψ

ABψφ, ψ αφ

ABφψ. 2.6

Subtracting both sides of the above two equations, we obtain φψ

φψ A1

0. 2.7

Ifφ /ψ, thenφψ−A1; this contradicts the hypothesis thatφ, ψ≥0.

bAssume thatkis even. Thenxnxn−kandφ, ψsatisfy the following system:

φ αψ

ABψψ, ψ αφ

ABφφ. 2.8

Subtracting both sides of the above two equations, we obtain φψ

A1 0. 2.9

Ifφ /ψ, thenA−1; this contradicts the hypothesis thatA≥0.

The proof is complete.

3. Boundedness and Invariant Interval

In this section, we will investigate the boundedness and invariant interval of1.1.

(6)

Theorem 3.1. Every solution of 1.1is bounded from above and from below by positive constants.

Proof. Let {xn}n−k be a positive solution of1.1. Clearly, if the solution is bounded from above by a constantM, then

xn1α

A B1M, forn≥ −k, 3.1

and so it is also bounded from below. Now for the sake of contradiction assume that the solution is not bounded from above. Then there exists a subsequence{xnm1}m0such that

mlim→ ∞nm∞, lim

m→ ∞xnm1∞ 3.2

and also

xnm1max{xn:nnm} form≥0. 3.3

From1.1we see that

xn1< α A 1

Axn forn≥0, 3.4

and so

mlim→ ∞xnm1 lim

m→ ∞xnm ∞. 3.5

Hence, for sufficiently largem,

0≤xnm1xnm αxnm

ABxnmxnm−kxnm α 1−ABxnmxnm−kxnm

ABxnmxnm−k <0, 3.6 which is a contradiction.

The proof is complete.

Let

f x, y

αx

ABxy. 3.7

Then the following statements are true.

Lemma 3.2. aAssume thatABα.Thenfx, yis increasing inxfor eachyand decreasing in yfor eachx.

bAssume that A < Bα.Thenfx, yis decreasing in y for eachx, decreasing in xfor y∈0, Bα−A, and increasing inxfory∈Bα−A,∞.

Proof. The proofs ofaandbare simple and will be omitted.

(7)

Theorem 3.3. Equation1.1possesses the following invariant intervals:

a 0,1/B,whenBαA;

b Bα−A,1/B,whenA < Bα < A1/B;

c 0, α/A,whenBαA1/B;

d 1/B, Bα−A,whenA1/B < Bα < Aα/A;

e 1/B, α/A,whenBαAα/A.

Proof. aSetgx αx/ABx, sogxis nondecreasing forxandg1/B≤ 1/Bif A, whenx−k, . . . , x0∈0,1/B; then we have

x1 αx0

ABx0x−kαx0

ABx0g 1 B

≤ 1

B. 3.8

The proof follows by induction.

bIn view ofLemma 3.2b, by using the monotonic character of the functionfx, y and the conditionA < Bα < A 1/B,whenx−k, . . . , x0∈Bα−A,1/B,we can get

x1 αx0

ABx0x−k fxo, x−kf A,1 B

> BαA, x1 αx0

ABx0x−k fxo, x−kf 1

B, Bα−A

1 B.

3.9

The proof follows by induction.

cSethx αx/ABxα/Aandgx αx/ABx, sohxis increasing andgxis decreasing forxifBαA1/B. In view ofLemma 3.2b, by using the monotonic character of the functionfx, y, whenx−k, . . . , x0 ∈0, α/A,we have

x1 αx0

ABx0x−kαx0

ABx0 α/Ah0>0, x1 αx0

ABx0x−kαx0

ABx0g0 α A.

3.10

The proof follows by induction.

dIn view ofLemma 3.2b, by using the monotonic character of the functionfx, y and the conditionA1/B < Bα < Aα/A, whenx−k, . . . , x0∈1/B, Bα−A,we obtain

x1 αx0

ABx0x−k fxo, x−kfBαA, BαA> 1 B, x1 αx0

ABx0x−k fxo, x−kf 1 B, 1

B

1

ABB1 < BαA.

3.11

The proof follows by induction.

(8)

e In view of the conditionAα/A, we can get Aα/A; by using the monotonic character of the function fx, y and the condition Aα/A, when x−k, . . . , x0∈1/B, α/A, we have

x1 αx0

ABx0x−k fxo, x−k A,α

A

Aαα A2Bαα≥ 1

B, x1 αx0

ABx0x−k fxo, x−kf 1 B,1

B

1 ABB1 < α

A.

3.12

The proof follows by induction.

The proof is complete.

4. Semicycles Analysis

We now give the definitions of positive and negative semicycles of a solution of1.4relative to an equilibrium pointx.

A positive semicycle of a solution{xn} of 1.4 consists of a string of terms {xl, xl1, . . . , xm}, all greater than or equal to the equilibriumx, withl≥ −kandm≤ ∞and such that

eitherl−k orl >−k, xl−1< x, either m∞orm <∞, xm1< x.

4.1

A negative semicycle of a solution {xn} of 1.4 consists of a string of terms {xl, xl1, . . . , xm}, all less than the equilibriumx, withl≥ −kandm≤ ∞and such that

eitherl−k orl >−k, xl−1x, eitherm∞ orm <∞, xm1x.

4.2

Theorem 4.1 see 12. Assume that fC0,∞ ×0,∞,0,∞ is such that fx, y is increasing inxfor each fixedyand is decreasing inyfor each fixedx. Letxbe a positive equilibrium of 1.12. Then the following are considered.

aIfk1, then every solution of 1.12has semicycles of length at least two.

bIfk2, then every solution of 1.12has semicycles that are either of length at leastk1 or of length at mostk1.

(9)

Let{xn}be a positive solution of1.1. Then one has the following identities:

xn1− 1 B 1

B

Axn−k

ABxnxn−k , n∈N0, 4.3

xn1α A −1

A

Bα−Axnαxn−k

ABxnxn−k , n∈N0, 4.4

xn1−Bα−A B1/B−Bα−AxnBα1/Bxn−k Axn−k−Bα−A

ABxnxn−k , n∈N0,

4.5 xn1x xxxn−k Bx−1/Bx−xn

ABxnxn−k , n∈N0, 4.6

x2k1n1x2k1n

Bx2k1n1−1

1/B−x2k1n

Ax2k1nBx2k1nk

ABx2k1n1−1

ABxn1x2k1n

αx2k1nk

1ABx2k1nk

αB/1ABx2k1n A

αAx2k1nx22k1n

ABx2k1n1−1

ABx2k1nkx2k1n

αx2k1nk , n∈N0.

4.7

If 1/BA, thenx1/Band4.3,4.7change into

xn1− 1 B 1

B

1/B−xn−k

ABxnxn−k, n∈N0, 4.8

x2k1n1x2k1n

1/B−x2k1n

ABx2k1n1−1Bx2k1nx2k1n1−1

ABx2k1n1−1

ABx2k1nkx2k1n

αx2k1nk

1/B−x2k1n

B2x2k1nkx2k1n1−1B2αx2k1nkAx2k1nABα

ABx2k1n1−1

ABx2k1nkx2k1n

αx2k1nk ,

n∈N0. 4.9

The following lemmas are straightforward consequences of identities4.3−4.9.

Lemma 4.2. Assume that

A, 4.10

(10)

and let{xn}be a solution of 1.1. Then the following statements are true.

ixn≤1/Bfor alln1.

iiIf for someN0,xN−kxandxNx, thenxN1x.

iiiIf for someN0,xN−k> xandxN< x, thenxN1< x.

iv0< x <1/B.

Lemma 4.3. Assume that

A < Bα < A 1

B, 4.11

and let{xn}be a solution of 1.1. Then the following statements are true.

iIf for someN0,xN< BαA, thenxNk1>1/B.

iiIf for someN0,xNA, thenxNk11/B.

iiiIf for someN0,xN> BαA, thenxNk1<1/B.

ivIf for someN0,BαA < xN<1/B, thenA < xNk1<1/B.

vIf for someN0,xN−kxandxNx, thenxN1x.

viIf for someN0,xN−k> xandxN< x, thenxN1< x.

viiIf for someN0,x2k1N< BαA, thenx2k1N1> x2k1N. viiiIf for someN0,x2k1N>1/B, thenx2k1N1< x2k1N.

ixA < x <1/B.

Lemma 4.4. Assume that

BαA 1

B, 4.12

and let{xn}be a solution of 1.1. Then the following statements are true.

iIf for someN0,xN>1/B, thenxNk1 <1/B.

iiIf for someN0,xN1/B, thenxNk1 1/B.

iiiIf for someN0,xN<1/B, thenxNk1 >1/B.

ivIf for someN0,x2k1N>1/B, thenx2k1N1< x2k1N. vIf for someN0,x2k1N<1/B, thenx2k1N1> x2k1N. vix1/B.

Lemma 4.5. Assume that

A 1

B < Bα < A α

A, 4.13

(11)

and let{xn}be a solution of 1.1. Then the following statements are true.

iIf for someN0,xN< BαA, thenxNk1>1/B.

iiIf for someN0,xNA, thenxNk11/B.

iiiIf for someN0,xN> BαA, thenxNk1<1/B.

ivIf for someN0, 1/B < xN< BαA, then 1/B < xNk1< BαA.

vIf for someN0,xN−kxandxNx, thenxN1x.

viIf for someN0,xN−k> xandxN> x, thenxN1< x.

viiIf for someN0,x2k1N<1/B, thenx2k1N1> x2k1N. viiiIf for someN0,x2k1N> BαA, thenx2k1N1< x2k1N.

ix1/B < x < Bα−A.

Lemma 4.6. Assume that

A α

A, 4.14

and let{xn}be a solution of 1.1. Then the following statements are true.

ixn< α/Afor alln1.

iiIf for someN0,xN< α/A, thenxNk1>1/B.

iiiIf for someN0, 1/B < xN< α/A, then 1/B < xNk1< α/A.

ivIf for someN0,xN−kxandxNx, thenxN1x.

vIf for someN0,xN−k> xandxN> x, thenxN1< x.

viIf for someN0,x2k1N<1/B, thenx2k1N1> x2k1N. viiIf for someN0,x2k1N> α/A, thenx2k1N1< x2k1N. viii1/B < x < α/A.

The following result is a consequence ofTheorem 4.1and Lemmas4.2−4.6.

Theorem 4.7. Let{xn}n−kbe a nontrivial solution of 1.1. Then the following statements are true.

aAssume thatBαA.Then, except possibly for the first semicycle, every oscillatory solution of 1.1has semicycles that are either of length at leastk1, or of length at mostk1.

bAssume that A < Bα < A1/B. Then, except possibly for the first semicycle, every oscillatory solution of 1.1 which lies in the invariant interval Bα − A,1/B has semicycles that are either of length at leastk1, or of length at mostk1.

cAssume that A1/B. Then, except possibly for the first semicycle,{xn}n−k is oscillatory and the sum of the lengths of two consecutive semicycles is equal to 2k1.

dAssume thatA1/B < Bα < Aα/A.Then, except possibly for the first semicycle, every oscillatory solution of 1.1which lies in the invariant interval1/B, Bα−Ahas semicycles at mostk1.

eAssume thatBαAα/A.Then, except possibly for the first semicycle, every oscillatory solution of 1.1which lies in the invariant interval1/B, α/Ahas semicycles at most k1.

(12)

5. Global Stability Proof

In this section, we will investigate the global stability of all positive solutions of1.1.

Theorem 5.1. Let{xn}n−kbe a positive solution of1.1. Then the following statements are true.

aAssume thatBαA.Then every solution of 1.1eventually enters the interval0,1/B.

bAssume thatA < Bα < A1/B.Then every solution of 1.1eventually enters the interval Bα−A,1/B.

cAssume thatA1/B < Bα < Aα/A.Then every solution of 1.1eventually enters the interval1/B, Bα−A.

dAssume thatBαAα/A.Then every solution of 1.1eventually enters the interval 1/B, α/A.

Proof. aIn view ofLemma 4.2, we know thatxn≤1/Bfor alln≥1 andx∈0,1/B; that is, all solutions of1.1eventually enter the interval0,1/B.

b If x−k, . . . , x0 ∈ Bα− A,1/B, by Theorem 3.3b, then we have xn ∈ Bα − A,1/B,for alln ≥ 0. If the initial conditions are not in the interval Bα− A,1/B, then we consider the 2k1th subsequences {x2k1nj}2k1j0 of the solution {xn}. We will give the proof for the subsequence{x2k1n}. The proof for all the other subsequences is similar and will be omitted. Without loss of generality, we assume that there existsN sufficiently large such thatx2k1N < BαAif x2k1N > 1/B, then the proof is similar and will be omitted; then in view of Lemmas4.3iiandiv, we know thatxk12N1 >1/B > Bα−A

and x2k1N1 < 1/B. If x2k1N1A, then, by induction, we know that the

former assertion implies that the result is true. Ifx2k1N1 < BαA, byLemma 4.3viii, then we can get x2k1N1 > x2k1N. It follows by induction that the subsequence {x2k1Nm}m0is increasing, and becausex2k1Nm< Bα−A, so limm→ ∞x2k1Nmexists and limm→ ∞x2k1NmA. However, taking limits by4.7, we get a contradiction.

cThe proof is similar tob, so will be omitted.

dIn view ofLemma 4.6, we know thatxn < α/Afor alln≥ 1; that is, all solutions of1.1eventually enter the interval 0, α/A. Furthermore, byTheorem 3.3,1/B, α/Ais an invariant interval of 1.1. Now, assume for the sake of contradiction that all solutions never enter the interval1/B, α/A, then the subsequence{x2k1Nm}m0enters the interval 0,1/B. Because x2k1N ≤ 1/Band Aα/A, then, by Lemma 4.6, we know that x2k1N1 > x2k1N; it follows by induction that the subsequence {x2k1Nm}m0 is increasing in the interval0,1/B. So limm→ ∞x2k1Nm exists and limm→ ∞x2k1Nm ≤ 1/B, which is a contradiction because1.1has no equilibrium point in the interval0,1/B.

The proof is complete.

Theorem 5.2. Assume that2.4holds. Then the positive equilibriumxof 1.1is a global attractor of all positive solutions of1.1.

We consider the following five cases.

Case 1. Assume thatBαA. By Theorems3.3aand5.1a, we know that1.1possesses an invariant interval0,1/Band every solution of1.1eventually enters the interval0,1/B.

Further, it is easy to see thatfx, yincreases inxand decreases inyin0,1/B.

(13)

Letm, M∈0,1/Bbe a solution of the system αm

ABmM m, αM

ABMmM, 5.1

which is equivalent to

αmAmBm2Mm, αMAMBM2Mm. 5.2

Hence

m−M1ABMm 0. 5.3

Now ifmM / 1−A/B, thenmM. For instance, this is the case ifA≥1 is satisfied.

IfmM 1−A/B, thenmandMsatisfy the system mM 1−A

B , mM α

1−B 5.4

and the equation

BB−1m2 B−1A−1m−0, 5.5

whose discriminant is

Δ B−1

B−1A−124B2α

. 5.6

Clearly, in this case,B < 1, and in view of condition2.4, we have Δ < 0, from which it follows thatmM. In view ofLemma 1.5,1.1has a unique equilibriumx∈0,1/Band every solution of1.1converges tox.

Case 2. Assume thatA < Bα < A1/B. By Theorems3.3band5.1b, we know that1.1 possesses an invariant intervalBα−A,1/Band every solution of1.1eventually enters the intervalBα−A,1/B. Further, it is easy to see thatfx, yincreases inxand decreases in yinBα−A,1/B. Then using the same argument in Case1,1.1has a unique equilibrium x∈Bα−A,1/Band every solution of1.1converges tox.

Case 3. Assume thatBαA1/B. Considering the 2k1th subsequences{x2k1nj}n0j∈ {0,1, . . . ,2k 1}, n ≥ 0, then by Lemma 4.4, we know that each one of the 2k 1th subsequences is above 1/B, below 1/B, or identically equal to 1/B. Furthermore, by the identity4.9, we can get that all 2k1th subsequences converge monotonically to limits, and for alln∈N,

x2k1n1x2k1n iffx2k1n 1

B. 5.7

(14)

So all the 2k1th subsequences{x2k1n1j}n0j ∈ {0,1, . . . ,2k1}converge to 1/B.

That is,x1/Bis a global attractor of1.1.

Case 4. Assume thatA1/B < Bα < Aα/A.By Theorems3.3dand5.1c, we know that 1.1possesses an invariant interval 1/B, Bα−A and every solution of1.1 eventually enters the interval 1/B, Bα− A. Furthermore, it is easy to see that the function fx, y decreases in each of its arguments in the interval1/B, Bα−A. Letm, M∈1/B, Bα−Abe a solution of the system

αm

ABmm M, αM

ABMM m, 5.8

that is, the solution of the system

αmAM B1mM, αMAm B1mM. 5.9

Thenm−MA1 0, which implies thatmM. EmployingLemma 1.6, we see that1.1 has a unique equilibriumx∈1/B, Bα−Aand every solution of1.1converges tox.

Case 5. Assume that Aα/A. By Theorems 3.3e and 5.1d, we know that 1.1 possesses an invariant interval1/B, α/Aand every solution of1.1eventually enters the interval1/B, α/A. Further, it is clear to see that the functionfx, ydecreases in each of its arguments in the interval1/B, α/A. Then, using the same argument as in Case4,1.1has a unique equilibriumx∈1/B, α/Aand every solution of1.1converges tox.

The proof is complete.

In view of Theorems2.1and5.2, we have the following result.

Theorem 5.3. Assume that2.4holds. Then the unique positive equilibriumxof 1.1is globally asymptotically stable.

References

1 M. R. S. Kulenovi´c, G. Ladas, L. F. Martins, and I. W. Rodrigues, “The dynamics ofxn1 α βxn/ABxnCxn−1: facts and conjectures,” Computers & Mathematics with Applications, vol. 45, no. 6–9, pp. 1087–1099, 2003.

2 M. R. S. Kulenovi´c and G. Ladas, Dynamics of Second Order Rational Difference Equations, Chapman Hall/CRC, Boca Raton, Fla, USA, 2002.

3 E. Camouzis and G. Ladas, Dynamics of Third-Order Rational Difference Equations with Open Problems and Conjectures, vol. 5 of Advances in Discrete Mathematics and Applications, Chapman Hall/CRC, Boca Raton, Fla, USA, 2008.

4 Y.-H. Su and W.-T. Li, “Global asymptotic stability of a second-order nonlinear difference equation,”

Applied Mathematics and Computation, vol. 168, no. 2, pp. 981–989, 2005.

5 R. P. Agarwal, W.-T. Li, and P. Y. H. Pang, “Asymptotic behavior of a class of nonlinear delay difference equations,” Journal of Difference Equations and Applications, vol. 8, no. 8, pp. 719–728, 2002.

6 A. M. Amleh, E. Camouzis, and G. Ladas, “On second-order rational difference equation. I,” Journal of Difference Equations and Applications, vol. 13, no. 11, pp. 969–1004, 2007.

7 A. M. Amleh, E. Camouzis, and G. Ladas, “On second-order rational difference equations. II,” Journal of Difference Equations and Applications, vol. 14, no. 2, pp. 215–228, 2008.

8 R. DeVault, W. Kosmala, G. Ladas, and S. W. Schultz, “Global behavior ofyn1 pyn−k/qyn yn−k,” Nonlinear Analysis: Theory, Methods & Applications, vol. 47, no. 7, pp. 4743–4751, 2001.

(15)

9 L.-X. Hu, W.-T. Li, and H.-W. Xu, “Global asymptotical stability of a second order rational difference equation,” Computers & Mathematics with Applications, vol. 54, no. 9-10, pp. 1260–1266, 2007.

10 V. L. Koci´c and G. Ladas, Global Behavior of Nonlinear Difference Equations of Higher Order with Applications, vol. 256 of Mathematics and Its Applications, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1993.

11 W. A. Kosmala, M. R. S. Kulenovi´c, G. Ladas, and C. T. Teixeira, “On the recursive sequenceyn1 pyn−1/qynyn−1,” Journal of Mathematical Analysis and Applications, vol. 251, no. 2, pp. 571–586, 2000.

12 M. Dehghan and N. Rastegar, “On the global behavior of a high-order rational difference equation,”

Computer Physics Communications, vol. 180, no. 6, pp. 873–878, 2009.

13 W.-T. Li and H.-R. Sun, “Dynamics of a rational difference equation,” Applied Mathematics and Computation, vol. 163, no. 2, pp. 577–591, 2005.

14 M. R. S. Kulenovi´c, G. Ladas, and N. R. Prokup, “A rational difference equation,” Computers &

Mathematics with Applications, vol. 41, no. 5-6, pp. 671–678, 2001.

15 M. R. S. Kulenovi´c, G. Ladas, and N. R. Prokup, “On the recursive sequencexn1 αxnβxn−1/1 xn,” Journal of Difference Equations and Applications, vol. 6, no. 5, pp. 563–576, 2000.

16 L.-X. Hu, W.-T. Li, and S. Stevi´c, “Global asymptotic stability of a second order rational difference equation,” Journal of Difference Equations and Applications, vol. 14, no. 8, pp. 779–797, 2008.

17 M. R. S. Kulenovi´c and O. Merino, “Global attractivity of the equilibrium ofxn1 pxnxn−1/qxn xn−1forq < p,” Journal of Difference Equations and Applications, vol. 12, no. 1, pp. 101–108, 2006.

18 R. Mazrooei-Sebdani and M. Dehghan, “Global stability ofyn1 pqynryn−k/1yn,” Applied Mathematics and Computation, vol. 182, no. 1, pp. 621–630, 2006.

19 S. Stevi´c, “On positive solutions of ak1th order difference equation,” Applied Mathematics Letters, vol. 19, no. 5, pp. 427–431, 2006.

20 Y.-H. Su and W.-T. Li, “Global attractivity of a higher order nonlinear difference equation,” Journal of Difference Equations and Applications, vol. 11, no. 10, pp. 947–958, 2005.

参照

関連したドキュメント

Luckhaus, Existence and regularity of weak solutions to the Dirichlet problem for semilinear elliptic systems of higher order, J.. ˆ Otani, Existence and nonexistence of

We also give a more concrete representation of the solutions to the nonlinear difference equation by calculating the special solution to the third-order homogeneous linear dif-

Ladas, Global Behavior of Nonlinear Difference Equations of Higher Order with Applications, Mathematics and Its Applications, volH. Ladas, Dynamics of Second Order Rational

Using a Razumikhin type theorem, we deduce sufficient conditions that guarantee the uniform asymptotic stability and boundedness of solutions of a scalar real fourth-order

Migda, “On the asymptotic behavior of solutions of higher order nonlinear difference equations,” Nonlinear Analysis: Theory, Methods &amp; Applications, vol.. Thandapani, “On

Global asymptotic stability of a second-order nonlinear difference equation, Appl. Sun, Dynamics of a rational difference

Xi, “Global asymptotic stability of a higher order rational di ff erence equation,” Journal of Mathematical Analysis and Applications, vol. Xi, “Global attractivity for a family

Ak¬n, Oscillatory behaviour of a higher-order nonlinear neutral type functional di¤erence equation with oscillating coe¢ cients, Appl.. Yildirim, Oscillation criteria for a certain