MISCELLANEOUS DYNAMIC EQUATIONS. Elvan Akın Bohner and Martin Bohner. 1. Introduction

Size: px
Start display at page:

Download "MISCELLANEOUS DYNAMIC EQUATIONS. Elvan Akın Bohner and Martin Bohner. 1. Introduction"

Transcription

1 MISCELLANEOUS DYNAMIC EQUATIONS Elvan Akın Bohner and Marin Bohner We consider several dynamic equaions and resen mehods on how o solve hese equaions. Among hem are linear equaions of higher order, Euler Cauchy equaions of higher order, logisic equaions or Verhuls equaions), Bernoulli equaions, Riccai equaions, and Clairau equaions.. Inroducion The calculus on ime scales has been inroduced in order o unify he heories of differenial equaions and of difference equaions and in order o exend hose heories o a more general class of so-called dynamic equaions. For an inroducion ino he calculus of ime scales we refer o [2, 3, 4]. Here we jus would like o menion ha any nonemy and closed subse T of R is called a ime scale, and ha he forward jum oeraor σ : T T is defined by σ) inf{s T : s }. Similarly, a backward jum oeraor ρ is defined, and he graininess µ of T is given by µ) σ). In his aer we consider cerain classes of simle dynamic equaions and resen mehods on how o solve hem. Discree versions of some of hose equaions are sudied in [, Chaer 3]. Firs, in Secion 2, we consider higher order linear dynamic equaions wih consan coefficiens, i.e., equaions of he form.) α k y k. For a reader no familiar wih he ime scales calculus, i is helful before sudying [2]) o hink of he oeraor as he usual derivaive if T R and he usual forward difference oeraor if T Z. Of course, as in he heory of ordinary differenial equaions, he so-called characerisic olynomial lays a cenral rôle, and in aricular he case of mulile roos of he characerisic olynomial will yield some ineresing resuls. Using a similar mehod, in Secion 3 we carry he sudy of Euler Cauchy equaions from he second order case as resened in [2, Secion 3.7] o he higher order case. However, in he higher order case we are forced due o he fac ha forward jum oeraors are no necessarily differeniable o facor hose Euler Cauchy equaions aroriaely, namely as.2) α k M k y, where he oeraors M k are defined recursively by.3) M y y and M k+ y M k y) for k N.

2 2 ELVAN AKIN BOHNER AND MARTIN BOHNER As in he higher order linear case, an associaed characerisic olynomial is imoran, and as for linear equaions, he case of mulile roos of his characerisic olynomial urns ou o be of aricular ineres. Moivaed by he fac ha for ordinary differenial equaions) he recirocal of a soluion of a linear equaion is a soluion of a logisic equaion, we inroduce our wo versions wo, as here are wo versions of linear equaions denoe v σ v σ).4) u )u + f) and v )v σ + f) of logisic equaions. These logisic equaions will aear in he forms.5) y [ ) + f)y)] y and x [) f)x)] x, where is defined as follows: The se R of all regressive i.e., + µ)) for all T) and rd-coninuous i.e., coninuous a oins wih σ) and lef-sided limis exis a oins wih ρ) ) funcions is an Abelian grou he so-called regressive grou) under he addiion defined by q + q + µq, and is he addiive inverse of R, i.e., / + µ) see [2, Exercise 2.26]). Now, knowledge of soluions of he linear equaions.4) hen enables us o exhibi soluions of he logisic equaions.5). Nex, we are aiming o inroduce Bernoulli equaions in such a way ha heir soluions raised o he αh ower are soluions of logisic equaions.5). However, in order o do so, we firs need o use he chain rule [2, Theorem.9] for finding a simle form of he derivaive of y α. This leads o he inroducion of a circle do mulilicaion defined by.6) α )) : α) + µ))h) α dh, which hen may be used o sudy Bernoulli equaions of he form [ )].7) x ) α f)xα ) x, and soluions of Bernoulli equaions see Secion 6) hen have our above described desired roery, so ha hey can readily be found as we know soluions of he logisic equaions.5). Furhermore, he circle do roduc urns ou o be ineresing in iself, and we devoe Secion 5 of his aer o he sudy of his roduc. Ineresing roeries of his roduc are he formulas.8) x α ) x α α x x, + µα ) + µ)α, and e α e α. Noe ha for R, he unique soluion of y )y saisfying y ) y see [2, Theorem 2.33]) is denoed as e, ), and he hird formula in.8) nicely sulemens he known relaions see [2, Theorem 2.36], and his is essenially a consequence of he roduc rule q) q + σ q and he quoien rule /q) q q )/qq σ )).9) e e q e q and e e q. e q Moreover, he se of osiively regressive funcions R +, i.e., he se of all rdconinuous funcions ha saisfy + µ)) > for all T see [2, Definiion

3 MISCELLANEOUS DYNAMIC EQUATIONS ]), ogeher wih addiion and real) scalar mulilicaion, urns ou o be a vecor sace, which we call he regressive vecor sace. In Secion 7 we also consider Riccai equaions of he form ) 2 z.) z + q) + q)z σ + ), ) where he circle square of x R is defined o be x 2 x) x). Proeries of he circle square are given in [2, Theorem 2.9], and we use hem o relae soluions of he Riccai equaion o soluions of he logisic equaions.5), rovided one aricular soluion of he Riccai equaion is known. Finally, in Secion 8 we sudy Clairau equaions aearing as.) y y + fy ), where f : R R may be any coninuously differeniable funcion. 2. Linear Equaions We consider homogeneous linear dynamic equaions of higher order wih consan coefficiens.), where y k are defined recursively by y y and y k+ y k) for k N and α k R for all k n wih α n. We call 2.) ϕλ) α k λ k he associaed characerisic olynomial. Theorem 2.. If λ R is a zero of 2.), hen e λ, ) solves.). Proof. According o our assumions, λ R, and hence y) e λ, ) is well defined. We find y k λ k y for all k N, and herefore α k y k α k λ k y ϕλ)y. Hence, if ϕλ), hen y solves.). Definiion 2.2. Equaion.) is called regressive if λ R for all zeros λ of ϕ. Theorem 2.3. Equaion.) is regressive if and only if 2.2) α k µ)) n k for all T. Proof. To show his, le λ, λ 2,..., λ n be he zeros of ϕ couning muliliciies. Then n ϕλ) α k λ k α n λ λ ).

4 4 ELVAN AKIN BOHNER AND MARTIN BOHNER The condiion λ R for all n is equivalen o n + µ)λ ). If n µ), hen α n + µ)λ ) α n n α k µ)) n k. If µ), hen n n [ α n + µ)λ ) α n µ)) n ) ] λ µ) µ)) n ϕ ) µ) n µ)) n α k ) k µ) α k µ)) n k. This roves our resul. Remark 2.4. Le λ, λ 2,..., λ n be he zeros of he characerisic olynomial 2.). Then i is easy o see ha.) is regressive iff n n λ i R, where λ i : λ λ 2... λ n. i Examle 2.5. Suose T 2 N and consider he equaion i 2.3) y 3y + 2y. The characerisic olynomial for 2.3) is λ 2 3λ + 2 and hence has zeros and 2. The regressiviy condiion 2.2) for 2.3) reads 2µ 2 ) + 3µ) + for all T, and i is clearly saisfied. Hence wo soluions of 2.3) are y ) e, ) log k ) and y 2 ) e 2, ) log k ). Now we consider linear equaions in he mulile roo case. The following lemma collecs some ideniies, ha are useful when rying o solve his roblem by he variaion of arameers mehod. Lemma 2.6. Suose λ R, le be differeniable, and assume ha here exiss a sequence { k } k N saisfying 2.4) and k k+ + λµ) for k N. Le y) )e λ, ). Then we have { k } k 2.5) y k ) λ k e λ for all k N and 2.6) { n } ϕ ) λ) α k y k e λ.! k

5 MISCELLANEOUS DYNAMIC EQUATIONS 5 Proof. We show 2.5) by inducion. Firs y y e λ e λ, and so 2.5) holds for k. If 2.5) holds for some k N, hen use he roduc rule) { k } ) k y k+ ) λ k e λ k ) k λ k { e λ + e σ } λ k ) k λ k { λe λ + } + µλ) e λ { k { ) k λ k+ + k k [ ) k k λ k+ + + } k + ) λ k+ { k+ } k ) + λ k e λ e λ )] λ k+ + k+ } so ha 2.5) holds for k +. This roves ha 2.5) holds for all k N. Now { k } k α k y k α k ) λ k e λ { n } k α k ) λ k e λ k { n } k! α k! k )! λk e λ k { n }! ϕ) λ) e λ, and his roves 2.6). Theorem 2.7. Suose λ R, and assume ha λ is a zero of 2.) wih muliliciy a leas m N. Pu m, k for all k m, and recursively k ) k+τ) for m 2 k. + λµτ) Then y) )e λ, ) is a soluion of.), where. Proof. Noe ha he sequence { k } k N defined as above is well defined only rdconinuous funcions are inegraed) and clearly saisfies 2.4). Since ϕ ) λ) for all < m, 2.6) from Lemma 2.6 imlies ha y) )e λ, ) solves.). e λ

6 6 ELVAN AKIN BOHNER AND MARTIN BOHNER Examle 2.8. If λ is a rile zero of ϕ, hen e λ, ), e λ, ) + λµτ), and e λ, ) τ s +λµs) + λµτ). are hree soluions of.). Noe ha for T R we have indeed e.g., + λµτ) dτ. 3. Euler Equaions Throughou his secion we assume ha T is a ime scale wih T, ). Moivaed by he fac ha soluions of Euler differenial equaions i.e., for T R) are of he form y) λ wih y ) λ y), we wan o look for soluions of he form y) e λ, ) also in he case of an Euler dynamic equaion on an arbirary ime scale. Noe ha for such y we have y ) λ y) and y ) λλ ) y), σ) and unforunaely y may no be differeniable see [2, Examle.56]) again for a general ime scale. Bu noe ha y ) λy), y )) λ2 y), y )) λ 2 y), y )) ) λ3 y),, and hence we wish o consider equaions of he form.2), where he oeraors M k are defined recursively by.3). We call.2) an Euler equaion of nh order, and as before 2.) is called he associaed characerisic olynomial. Theorems 3. and 3.3 below may be shown in a way comleely analogous o he roofs of Theorems 2. and 2.3, and hence we omi heir roofs. Theorem 3.. If λ/ R and λ is a zero of 2.), hen e λ, ) solves.2). Definiion 3.2. Equaion.2) is called regressive if λ Theorem 3.3. Equaion.2) is regressive if and only if 3.) α k k µ)) n k for all T. R for all zeros λ of ϕ. Examle 3.4. If n 2, hen.2) can be rewrien as ) α 2 y ) + α y + α y α 2 [σ)y α + + y + α ] y. α 2 α 2 Hence.2) is equivalen o 3.2) σ)y + ay + by wih a α α 2 +, b α α 2.

7 MISCELLANEOUS DYNAMIC EQUATIONS 7 The characerisic olynomial of 3.2) has he same zeros as λ 2 + a )λ + b. The regressiviy condiion 3.) for 3.2) is equivalen o σ) aµ) + bµ)) 2 for all T. Examle 3.5. If n 3, hen.2) can be exanded as in Examle 3.4 if T is a ime scale wih differeniable forward jum σ. E.g., if T R, hen 3.3) α 3 y ) ) + α 2 y ) + α y + α y, which is equivalen o 3.4) 3 y +a 2 y +by +cy wih a 3+ α 2 α 3, b + α α 3 + α 2 α 3, c α α 3. Noe ha he characerisic equaion for 3.4) reads λ 3 + a 3)λ 2 + b a + 2)λ + c while he characerisic equaion for 3.3) reads in a sense more naural and suggesive simly α 3 λ 3 + α 2 λ 2 + α λ + α. This remark alies o Euler equaions of any order and on any ime scale: Characerisic equaions for Euler equaions in facored form are easier o remember han hose for he equaions in exanded form. Examle 3.6. If ϕ) or ϕ), hen we have on any ime scale) e, ) and e, ). Examle 3.7. I is easy o verify e λ, ) ) λ if T R, Γ + λ) Γ ) e λ, ) Γ) Γ + λ) where Γ is he gamma funcion, and if T Z, e λ, ) log q [+q )λ] if T q N wih q >. Now we consider Euler equaions in he mulile roo case. The following lemma collecs some ideniies ha are useful when rying o solve his roblem by he variaion of arameers mehod. Is roof is comleely analogous o he roof of Lemma 2.6, and herefore we omi i. Lemma 3.8. Suose λ R, le be differeniable, and assume ha here exiss a sequence { k } k N saisfying 3.5) and k k+ + λµ) for k N. Le y) )e λ, ). Then we have { k 3.6) M k y } k ) λ k e λ

8 8 ELVAN AKIN BOHNER AND MARTIN BOHNER for all k N and 3.7) { n α k M k y } ϕ ) λ)! As in Theorem 2.7, we now may use Lemma 3.8 o obain he following resul. Theorem 3.9. Suose λ R, and assume ha λ is a zero of 2.) wih muliliciy a leas m N. Pu m, k for all k m, and recursively k ) e λ. k+τ) for m 2 k. τ + λµτ) Then y) )e λ, ) is a soluion of.2), where. Examle 3.. If λ is a rile zero of ϕ, hen e λ, ), e λ, ) τ + λµτ), and e, λ ) τ are hree soluions of.2). Noe ha for T R we have indeed e.g., τ + λµτ) dτ τ log. s s+λµs) τ + λµτ) Examle 3.. Suose λ is a double zero of a characerisic olynomial 2.) for T q N where q > ), and le. Then wo soluions of he corresonding Euler q-difference equaion.2) are given by y ) log q [+q )λ] see Examle 3.7) and y 2 ) [ + q )λ] y ) τ + λµτ) y ) τ. Noe ha q m τ m σq k ) q k τ m µq k m ) q )q k q k q k mq ), and so y 2 ) q ) log q [+q )λ] log q. In aricular, if q 2, hen y ) log 2 +λ) and y 2 ) log 2 +λ) log 2. Examle 3.2. Le T N 2 and consider he equaion + ) 2 y y + y. Here, σ) + ) 2 and µ) + 2. The characerisic olynomial is λ 2 2λ +. Hence wo soluions are given by see Examle 3.6) y ) e, ) and y 2) y ) Noe ha m 2 m στ) k σk 2 ) k 2 m στ) k τ + µτ) µk 2 m ) σk 2 ) m k στ). ) 2 k, k +

9 MISCELLANEOUS DYNAMIC EQUATIONS 9 and hence y 2 ) k ) 2 k. k + 4. Logisic Equaions Le us firs recall he following resuls from [2, Theorems 2.77 and 2.74]. Theorem 4. Variaion of Consans). Suose R and f C rd. Le T. i) Le u R. The unique soluion u saisfying he lef equaion in.4) and u ) u is given by u) e, )u + e, στ))fτ). ii) Le v R. The unique soluion v saisfying he righ equaion in.4) and v ) v is given by v) e, )v + e, τ)fτ). Moivaed by he fac ha for differenial equaions) u /y solves a linear equaion if y is a soluion of he logisic equaion, we assume ha u is a soluion of he firs linear equaion in.4), and hen y /u saisfies ) y u u uu σ u + f uu σ + fy)y σ [ + fy)] y. Therefore we call he firs equaion in.5) wih + fy R) a logisic dynamic equaion or Verhuls equaion). Similarly, we could sar wih a soluion v of he second linear equaion in.4) and find ha x /v saisfies ) x v v vv σ vσ f vv σ fx σ )x [ fx] x. We also call he second equaion in.5) wih fx R) a logisic dynamic equaion. Throughou we assume R and f C rd. Noe ha if u, hen y /u solves he firs equaion in.5) and auomaically saisfies + fy R, as can be seen from he calculaion + µ + fy) u + µu + f) u u + µu u uσ u. Similarly, if v, hen x /v solves he second equaion in.5) and saisfies + µfx v + µf v vσ + µf v ) v vσ + µv σ v vσ + µ) v so ha fx R follows. Using Theorem 4., i is now easy o solve logisic equaions as follows. Theorem 4.2. Suose R and f C rd.

10 ELVAN AKIN BOHNER AND MARTIN BOHNER i) Le y. If u) e, ) + e, στ))fτ) for all T, y hen y) u) e, ) y + e στ), )fτ) solves he firs equaion in.5) and saisfies y ) y. ii) Le x. If hen v) e, ) x + e, τ)fτ) for all T, x) v) e, ) x + e τ, )fτ) solves he second equaion in.5) and saisfies x ) x. In alicaions e.g., oulaion dynamics) one ofen assumes ha here exiss a consan N such ha 4.) ) Nf) for all T. For he remainder of his secion we assume 4.). Then we can evaluae he inegral in he denominaor of he soluion given in Theorem 4.2 ii) exlicily as e τ, )fτ) e τ, ) τ) N N [e, ) ]. Hence we obain he following corollary from Theorem 4.2 ii). Corollary 4.3. Suose R and le N be a consan. Le x. If hen 4.2) x) saisfies x ) x and solves x N + e, ) for all T, N e, ) x N + e,) N 4.3) x x x N + µ N x ) N + wih x N x N R. ) e, ) Remark 4.4. As is common in oulaion dynamics, we call he funcion from Corollary 4.3 he inrinsic growh funcion, while we refer o N as he sauraion level or environmenal carrying caaciy. Noe also ha he funcions x ) and x 2 ) N are soluions of 4.3), he so-called equilibrium soluions or criical soluions). If he saring value of x is beween hese wo soluions, i.e., x, N), hen /x /N > and hence x N + e, ) > e, ) N N

11 MISCELLANEOUS DYNAMIC EQUATIONS so ha < x) < e, ) N e, ) N, rovided R+. Noe ha R + ensures ha e, ) is osiive see [2, Theorem 2.48]). x > N, hen /x /N < and hence x) > N, rovided R +. Assume now ha ) > for all T. Then R +. If x, N), hen < x) < N and ) x x x N + µ N x > so ha x is sricly increasing. Similarly, if x > N, hen x is sricly decreasing. In any case, if x and if T is unbounded above, hen lim x) N, rovided lim e, ). Examle 4.5. Le T R, le α be a consan, and consider he coninuous) logisic equaion 4.4) x [α αx)] x, i.e., x αx x). Everyhing in Remark 4.4 alies o 4.4) here ) α and N ), and by Corollary 4.3 he soluion of 4.4) wih xs) z where s R and z R \ {}) is x) e α s) z + eα s) + z ) e. αs ) Examle 4.6. Le T Z and consider he discree) logisic equaion 4.5) x [ x)] x, i.e., x x x) + x. Everyhing in Remark 4.4 alies o 4.5) here ) and N ), and by Corollary 4.3 he soluion of 4.5) wih xs) z where s Z and z R \ {}) is x)! s! z +! s! + z ). s!! Remark 4.7. One ofen finds an equaion of he form 4.6) y n+ αy n y ) n N referred o as a logisic difference equaion. Observe ha his equaion and our logisic difference equaion inroduced in his secion are no he same. Our equaion has he advanage ha i is easily accessible wih he mehods used o sudy he normal logisic differenial equaion. Even hough 4.6) aears as a naural analogue of he logisic differenial equaion, he resuls given in his secion indicae ha our equaion deserves he name logisic equaion raher han 4.6) which is of course also worhy of sudy for is richness in dynamics). If

12 2 ELVAN AKIN BOHNER AND MARTIN BOHNER 5. The Regressive Vecor Sace For α R, we wan o define he Bernoulli equaion in such a way, ha he subsiuion x x α ransforms a soluion x of he Bernoulli equaion ino a soluion x of he logisic equaion. To do so, we now find a formula for he derivaive of x α. This ask leads us in his secion o he inroducion of a circle do mulilicaion, which urns ou o make R +,, ) ino a vecor sace, he so-called regressive vecor sace. To begin wih, we use he chain rule [2, Theorem.9] o calculae If x), hen x ) x α ) ) x ) 5.) x ) x) x ) x) α x) + µ)x )h ) α dh. ) α α + µ) x ) x) h dh. In order o have everyhing well defined, we wan o assume ha, if α R \ N, 5.2) + µ) x ) h > for all T and all h [, ]. x) I is easy o see ha x /x R + is sufficien for 5.2). Le us inroduce he noaion { R if α N Rα) : R + if α R \ N. Noe ha R + imlies ha + µ))h > for all T and all h [, ]. Definiion 5.. For α R and Rα) we define α by.6). Remark 5.2. Noe ha α α rovided T R. Theorem 5.3. Le α R. If α N, suose ha x) for all T. If α N, suose ha x)x σ ) > for all T. Then he firs formula in.8) holds. Proof. Firs noe ha +µx /x x σ /x imlies x /x Rα). Then he saemen follows from 5.). Examle 5.4. I is easy o see ha 2, Remark 5.5. The hree formulas µ, and 2 ) 2. xy) x xy x y y, x/y) x x/y x y y, x α ) x α α x x could be used o define a logarihm as log f, ) f τ) fτ),

13 MISCELLANEOUS DYNAMIC EQUATIONS 3 bu hen he usual logarihm rules would no hold as we had formulas involving an addiional inegral, e.g., log xy, ) log x, ) + log y, ) + µτ)x τ)y τ). xτ)yτ) I is sill an oen roblem o define a logarihm ha saisfies somehow smooher logarihm rules. Theorem 5.6. Le α R. If Rα), hen α R. More recisely, he second formula in.8) holds. Proof. Le α R and Rα). Then + µα ) + µα + µh) α dh + +µ where we used he subsiuion s + µh. Therefore, α R. αs α ds + µ) α, The following resul emhasizes again he imorance of our circle do mulilicaion inroduced in Definiion 5.. Theorem 5.7. If α R and Rα), hen he hird formula in.8) holds. Proof. Firs noe ha Rα) imlies α R by Theorem 5.6. Now we le T and u y e α, ). Then y ) and by Theorem 5.3 we ski he second argumen ) ) y e α ) Hence y solves he iniial value roblem α e e e α α )y. y α ))y, y ). Bu his iniial value roblem has exacly one soluion, namely e α, ). We will use Theorem 5.7 o show ha R +,, ) saisfies he axioms of a vecor sace. To do so, he following easy auxiliary resul is helful. Lemma 5.8. Suose, q R. If e e q, hen q. Proof. Differeniae boh sides of e e q o obain his resul. Now we are ready o rove he main resul of his secion. Theorem 5.9. R +,, ) is a real vecor sace. Proof. As is known see e.g., [2, Lemma 2.47]), R +, ) is an Abelian grou. We now rove he remaining four vecor sace roeries 5.3) α β ) αβ) for all α, β R, R + ; 5.4) for all R + ; 5.5) α q) α ) α q) for all, q R +, α R; 5.6) α + β) α ) β ) for all α, β R, R +.

14 4 ELVAN AKIN BOHNER AND MARTIN BOHNER Clearly, 5.4) follows direcly from Definiion 5.. For 5.3), 5.5), and 5.6) we use Theorem 5.7, Lemma 5.8, he rules.9) for he exonenial funcion, and he normal rules of exonens as follows: imlies 5.3), e α β ) e α β e α ) β e αβ e αβ) imlies 5.5), and imlies 5.6). e α q) e α q e e q ) α e α e α q e α e α q e α ) α q) e α+β) e α+β e α e β e α e β e α ) β ) Remark 5.. Noe ha 5.4) of course holds for each R. Furhermore, if α, β N, hen 5.3), 5.5), and 5.6) also hold for, q R. 6. Bernoulli Equaions Now we le α R \ {} and consider he Bernoulli equaion.7). Noe ha.7) is in he form of he second equaion in.5) if α. Inroducing x x α, we find if x never vanishes and solves.7), hen x x xα ) x α α x x α [ )] α fxα α ) f x, where we used Theorem 5.3 and he vecor sace roeries 5.3) and 5.5). Hence x solves he logisic equaion 6.) x [α ) f x)] x. Equaion 6.) is of he form.5), so ha we may aly Theorem 4.2 ii) use also Theorem 5.7) o find x as e α, ) x) x + e α τ, )fτ). I is hence easy o show he following resul. Theorem 6.. Suose α R \ {}, Rα), and f C rd. Le x. If hen x α + e α τ, )fτ) > for all T, x) [ x α solves he Bernoulli equaion.7). e, ) + e α τ, )fτ) Examle 6.2. Consider he iniial value roblem 6.2) x x + x 4, x). ] /α

15 MISCELLANEOUS DYNAMIC EQUATIONS 5 Noe ha he differenial equaion in 6.2) can be rewrien as a dynamic wih T R) Bernoulli equaion [ )] x α fxα x wih ), f) 3, α 3. According o Theorem 6. wih and x ), he soluion of 6.2) is x) [ x α e, ) + e α τ, )fτ) ] /α Examle 6.3. Consider he iniial value roblem e [ ) e 3 ] /3. 6.3) x 5x x3 + 5x + x 2 + x x 2, x). Noe ha 5 x x 2 + x x 5 x 2 ) x x2 imlies ha he difference equaion in 6.3) can be rewrien as a dynamic wih T Z) Bernoulli equaion [ )] x α fxα x wih ) 5, f), α 2. According o Theorem 6. wih and x ), he soluion of 6.3) is x) [ x α e, ) + ] /α e α τ, )fτ) 7. Riccai Equaions In his secion we consider he Riccai equaion 7.) z + q) + r)z σ + z 2 ) + µ)z, where, q, r are assumed o be rd-coninuous such ha ) for all T. Noe ha 7.) can be wrien in he form.). We now assume ha z is a known soluion of.). We ick any oher soluion z of.) and consider he difference y z z. I follows ha z y q rz σ [ z ) 2 ry σ + ) 2 + q + rz σ + ) ] 2 z. ) 2 z To coninue his calculaion, he following lemma is useful and can be roved by a direc comuaion. Lemma 7.. For, q R we have 2 q 2 [ q)] q ).

16 6 ELVAN AKIN BOHNER AND MARTIN BOHNER We may use Lemma 7. o coninue our calculaion from above as follows: [ z ) 2 ) ] 2 z y ry σ + [ ry σ z + z )] z z ) [ ry σ z + y z )]. ) Denoing s z z, we coninue o find y ry σ + ys ry rµy + ys ys r) rµy and herefore solve his las equaion for y ) 7.2) y y s r ys r). + µr Now, aly he formula o find where a b + c) a + b + c + µab + µac a b) + c + µa) s r) z z r z [ 2 z ) ] r 7.3) f ) + µ r z r z + y and [ + µ ) + y z r Hence, using 7.2), we find ha y saisfies he equaion 7.4) y [ g) + f)y)] y. )] g r 2 z ). g + fy, Equaion 7.4) is a logisic equaion as in.5). Hence he following resul holds. Theorem 7.2. Suose z is a soluion of he Riccai equaion.). Define f and g by 7.3) and le y be a soluion of he logisic equaion 7.4). Then z z +y is also a soluion of he Riccai equaion.). Corollary 7.3. Suose z is a soluion of he Riccai equaion.). Define f and g by 7.3). Then a soluion z of.) saisfying z ) z is given by e g, ) z) z ) + z + z ) e g, στ))fτ), rovided none of he denominaors is zero. Proof. This follows from Theorem 4.2, where soluions of.5) are given.

17 MISCELLANEOUS DYNAMIC EQUATIONS 7 Examle 7.4. We consider he equaion 7.5) z + [ )] z σ + µ)z. Equaion 7.5) is a Riccai equaion of he form 7.) wih coefficiens ), q), and r) ). By examinaion, z ) solves 7.5). We have z ) ) and r z ) ) [ )] ) so ha f and g from 7.3) are given by f) g). Hence, by Corollary 7.3, we find ha anoher soluion of 7.5) is given by z 2 z) ), z e ), ) rovided none of he demoninaors is zero. This soluion z saisfies z ) z. 8. Clairau Equaions The Clairau dynamic equaion aears as in.), where f : R R is some coninuously differeniable funcion. We roceed o rea.) as is coninuous analogue y y + fy ) and erform he subsiuion v y. Then 8.) y) v) + fv)). Differeniaing 8.) and using [2, Theorem.9], we conclude so ha v) y ) σ)v ) + v) + v ) { v ) σ) + Hence we obain he following resul. Theorem 8.. For any consan c R, 8.2) y) c + fc) f v) + hµ)v ) ) dh f v) + hµ)v ) ) } dh. is a soluion of.). Furher soluions of.) may be obained by solving 8.3) Examle 8.2. Consider he equaion f v) + hµ)v ) ) dh σ), y v. 8.4) y y + y ) 2, which is of he form.) wih fx) x 2. By 8.2) of Theorem 8., we find ha y) c + c 2 is a soluion of 8.4), where c R is an arbirary consan. Wih v y, he lef hand side of 8.3) becomes f v) + hµ)v ) ) dh 2 v) + hµ)v ) ) dh 2v) + µ)v ).

18 8 ELVAN AKIN BOHNER AND MARTIN BOHNER Hence we mus solve 8.5) 2v + µ)v σ), i.e., v + v σ σ). If T R, hen 8.5) becomes 2v), and herefore y) 2 4 of 8.4). Now we assume ha 8.6) T { k : k N } wih < < 2 <.... Then we may divide 8.5) by µ) o arrive a v 2 µ) vσ + σ) µ), is also a soluion which is an equaion as in.4) wih 2/µ and f σ/µ. By Theorem 4., a soluion is given by noe ha 2/µ) 2/µ) We find v) e 2/µ, )v + e 2/µ, τ) στ) µτ). m v m ) ) m v + ) m k k+ ) m c ) [α m + ] m ) k k, where we relaced v by c α. The value of he above sum deends on he ime scale, and we will do he calculaion for a few ime scales ha saisfy 8.6). Firs, m m T N imlies ) k k ) m 2 ) 4 4, k so we choose α /4 and find ha y) ) ) c ) c2 k c ) ) is also a soluion of 8.4) for any c R see also [, Examle 3..]). Nex, m T q N imlies ) k k ) m qm+ q + q q +, k so we choose α q/q + ) and derive ha y) q + q c )log q q2 q + ) 2 + c2 c ) log q is also a soluion of 8.4) for any c R. Finally, we noe ha m T N 2 imlies ) k m mm + ) k ), 2 k and hence we choose α and find ha y) c ) c2 is also a soluion of 8.4) for any c R. c ) ) 2 q ) 2 2q + ) 4 2 ) 2 2 4

19 MISCELLANEOUS DYNAMIC EQUATIONS 9 References [] R. P. Agarwal. Difference Equaions and Inequaliies. Marcel Dekker, Inc., New York, 992. [2] M. Bohner and A. Peerson. Dynamic Equaions on Time Scales: An Inroducion wih Alicaions. Birkhäuser Boson Inc, Boson, MA, 2. [3] S. Hilger. Analysis on measure chains a unified aroach o coninuous and discree calculus. Resuls Mah., 8:8 56, 99. [4] V. Lakshmikanham, B. Kaymakçalan, and S. Sivasundaram. Dynamic Sysems on Measure Chains, volume 37 of Mahemaics and is Alicaions. Kluwer Academic Publishers Grou, Dordrech, 996. Received Ocober 8, 22 Universiy of Missouri Rolla, Dearmen of Mahemaics and Saisics, Rolla, MO 654 address: akine@umr.edu, bohner@umr.edu h:// akine, h:// bohner

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

THREE POSITIVE SOLUTIONS OF A THREE-POINT BOUNDARY VALUE PROBLEM FOR THE p-laplacian DYNAMIC EQUATION ON TIME SCALES 1.

THREE POSITIVE SOLUTIONS OF A THREE-POINT BOUNDARY VALUE PROBLEM FOR THE p-laplacian DYNAMIC EQUATION ON TIME SCALES 1. Commun. Oim. Theory 218 (218, Aricle ID 13 hs://doi.org/1.23952/co.218.13 THREE POSITIVE SOLUTIONS OF A THREE-POINT BOUNDARY VALUE PROBLEM FOR THE -LAPLACIAN DYNAMIC EQUATION ON TIME SCALES ABDULKADIR

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

Stopping Brownian Motion without Anticipation as Close as Possible to its Ultimate Maximum

Stopping Brownian Motion without Anticipation as Close as Possible to its Ultimate Maximum Theory Probab. Al. Vol. 45, No.,, (5-36) Research Reor No. 45, 999, De. Theore. Sais. Aarhus Soing Brownian Moion wihou Aniciaion as Close as Possible o is Ulimae Maximum S. E. GRAVERSEN 3, G. PESKIR 3,

More information

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL: Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

Undetermined coefficients for local fractional differential equations

Undetermined coefficients for local fractional differential equations Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani

More information

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Novi Sad J. Mah. Vol. 32, No. 2, 2002, 95-108 95 POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Hajnalka Péics 1, János Karsai 2 Absrac. We consider he scalar nonauonomous neural delay differenial

More information

Chapter 3 Common Families of Distributions

Chapter 3 Common Families of Distributions Chaer 3 Common Families of Disribuions Secion 31 - Inroducion Purose of his Chaer: Caalog many of common saisical disribuions (families of disribuions ha are indeed by one or more arameers) Wha we should

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

Optimality Conditions for Unconstrained Problems

Optimality Conditions for Unconstrained Problems 62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

arxiv: v1 [math.ca] 15 Nov 2016

arxiv: v1 [math.ca] 15 Nov 2016 arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy

More information

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at BOX-JENKINS MODEL NOAION he Box-Jenkins ARMA(,q) model is denoed b he equaion + + L+ + a θ a L θ a 0 q q. () he auoregressive (AR) ar of he model is + L+ while he moving average (MA) ar of he model is

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Hamilton Jacobi equations

Hamilton Jacobi equations Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion

More information

Ordinary dierential equations

Ordinary dierential equations Chaper 5 Ordinary dierenial equaions Conens 5.1 Iniial value problem........................... 31 5. Forward Euler's mehod......................... 3 5.3 Runge-Kua mehods.......................... 36

More information

On the Existence, Uniqueness and Stability Behavior of a Random Solution to a Non local Perturbed Stochastic Fractional Integro-Differential Equation

On the Existence, Uniqueness and Stability Behavior of a Random Solution to a Non local Perturbed Stochastic Fractional Integro-Differential Equation On he Exisence, Uniqueness and Sabiliy ehavior of a Random Soluion o a Non local Perurbed Sochasic Fracional Inegro-Differenial Equaion Mahmoud M. El-orai,*, M.A.Abdou, Mohamed Ibrahim M. Youssef Dearmen

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

f(t) dt, x > 0, is the best value and it is the norm of the

f(t) dt, x > 0, is the best value and it is the norm of the MATEMATIQKI VESNIK 66, 1 (214), 19 32 March 214 originalni nauqni rad research aer GENERALIZED HAUSDORFF OPERATORS ON WEIGHTED HERZ SPACES Kuang Jichang Absrac. In his aer, we inroduce new generalized

More information

Math Spring 2015 PRACTICE FINAL EXAM (modified from Math 2280 final exam, April 29, 2011)

Math Spring 2015 PRACTICE FINAL EXAM (modified from Math 2280 final exam, April 29, 2011) ame ID number Mah 8- Sring 5 PRACTICE FIAL EXAM (modified from Mah 8 final exam, Aril 9, ) This exam is closed-book and closed-noe You may use a scienific calculaor, bu no one which is caable of grahing

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

Description of the MS-Regress R package (Rmetrics)

Description of the MS-Regress R package (Rmetrics) Descriion of he MS-Regress R ackage (Rmerics) Auhor: Marcelo Perlin PhD Suden / ICMA Reading Universiy Email: marceloerlin@gmail.com / m.erlin@icmacenre.ac.uk The urose of his documen is o show he general

More information

Perturbations of the half-linear Euler Weber type differential equation

Perturbations of the half-linear Euler Weber type differential equation J. Mah. Anal. Al. 323 2006 426 440 www.elsevier.com/locae/jmaa Perurbaions of he half-linear Euler Weber ye differenial equaion Ondřej Došlý Dearmen of Mahemaics, Masaryk Universiy, Janáčkovo nám. 2a,

More information

Dirac s hole theory and the Pauli principle: clearing up the confusion.

Dirac s hole theory and the Pauli principle: clearing up the confusion. Dirac s hole heory and he Pauli rincile: clearing u he conusion. Dan Solomon Rauland-Borg Cororaion 8 W. Cenral Road Moun Prosec IL 656 USA Email: dan.solomon@rauland.com Absrac. In Dirac s hole heory

More information

A note to the convergence rates in precise asymptotics

A note to the convergence rates in precise asymptotics He Journal of Inequaliies and Alicaions 203, 203:378 h://www.journalofinequaliiesandalicaions.com/conen/203//378 R E S E A R C H Oen Access A noe o he convergence raes in recise asymoics Jianjun He * *

More information

t 2 B F x,t n dsdt t u x,t dxdt

t 2 B F x,t n dsdt t u x,t dxdt Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.

More information

Rise-Time Distortion of Signal without Carrying Signal

Rise-Time Distortion of Signal without Carrying Signal Journal of Physics: Conference Series PAPER OPEN ACCESS Rise-Time Disorion of Signal wihou Carrying Signal To cie his aricle: N S Bukhman 6 J. Phys.: Conf. Ser. 738 8 View he aricle online for udaes and

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

The Miki-type identity for the Apostol-Bernoulli numbers

The Miki-type identity for the Apostol-Bernoulli numbers Annales Mahemaicae e Informaicae 46 6 pp. 97 4 hp://ami.ef.hu The Mii-ype ideniy for he Aposol-Bernoulli numbers Orli Herscovici, Toufi Mansour Deparmen of Mahemaics, Universiy of Haifa, 3498838 Haifa,

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

EIGENVALUE PROBLEMS FOR SINGULAR MULTI-POINT DYNAMIC EQUATIONS ON TIME SCALES

EIGENVALUE PROBLEMS FOR SINGULAR MULTI-POINT DYNAMIC EQUATIONS ON TIME SCALES Elecronic Journal of Differenial Equaions, Vol. 27 (27, No. 37, pp. 3. ISSN: 72-669. URL: hp://ejde.mah.xsae.edu or hp://ejde.mah.un.edu EIGENVALUE PROBLEMS FOR SINGULAR MULTI-POINT DYNAMIC EQUATIONS ON

More information

APPLIED ECONOMETRIC TIME SERIES (2nd edition)

APPLIED ECONOMETRIC TIME SERIES (2nd edition) INSTRUCTOR S RESOURCE GUIDE TO ACCOMPANY APPLIED ECONOMETRIC TIME SERIES (2nd ediion) Waler Enders Universiy of Alabama Preared by Pin Chung American Reinsurance Comany and Iowa Sae Universiy Waler Enders

More information

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013 IMPLICI AND INVERSE FUNCION HEOREMS PAUL SCHRIMPF 1 OCOBER 25, 213 UNIVERSIY OF BRIISH COLUMBIA ECONOMICS 526 We have exensively sudied how o solve sysems of linear equaions. We know how o check wheher

More information

First-Order Recurrence Relations on Isolated Time Scales

First-Order Recurrence Relations on Isolated Time Scales Firs-Order Recurrence Relaions on Isolaed Time Scales Evan Merrell Truman Sae Universiy d2476@ruman.edu Rachel Ruger Linfield College rruger@linfield.edu Jannae Severs Creighon Universiy jsevers@creighon.edu.

More information

Math 2214 Sol Test 2B Spring 2015

Math 2214 Sol Test 2B Spring 2015 Mah 14 Sol Tes B Sring 015 roblem 1: An objec weighing ounds sreches a verical sring 8 fee beond i naural lengh before coming o res a equilibrium The objec is ushed u 6 fee from i s equilibrium osiion

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

Sections 2.2 & 2.3 Limit of a Function and Limit Laws

Sections 2.2 & 2.3 Limit of a Function and Limit Laws Mah 80 www.imeodare.com Secions. &. Limi of a Funcion and Limi Laws In secion. we saw how is arise when we wan o find he angen o a curve or he velociy of an objec. Now we urn our aenion o is in general

More information

On Oscillation of a Generalized Logistic Equation with Several Delays

On Oscillation of a Generalized Logistic Equation with Several Delays Journal of Mahemaical Analysis and Applicaions 253, 389 45 (21) doi:1.16/jmaa.2.714, available online a hp://www.idealibrary.com on On Oscillaion of a Generalized Logisic Equaion wih Several Delays Leonid

More information

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant). THE WAVE EQUATION 43. (S) Le u(x, ) be a soluion of he wave equaion u u xx = 0. Show ha Q43(a) (c) is a. Any ranslaion v(x, ) = u(x + x 0, + 0 ) of u(x, ) is also a soluion (where x 0, 0 are consans).

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0. Time-Domain Sysem Analysis Coninuous Time. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 1. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 2 Le a sysem be described by a 2 y ( ) + a 1

More information

Lecture 6 - Testing Restrictions on the Disturbance Process (References Sections 2.7 and 2.10, Hayashi)

Lecture 6 - Testing Restrictions on the Disturbance Process (References Sections 2.7 and 2.10, Hayashi) Lecure 6 - esing Resricions on he Disurbance Process (References Secions 2.7 an 2.0, Hayashi) We have eveloe sufficien coniions for he consisency an asymoic normaliy of he OLS esimaor ha allow for coniionally

More information

Essential Maps and Coincidence Principles for General Classes of Maps

Essential Maps and Coincidence Principles for General Classes of Maps Filoma 31:11 (2017), 3553 3558 hps://doi.org/10.2298/fil1711553o Published by Faculy of Sciences Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma Essenial Maps Coincidence

More information

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function Elecronic Journal of Qualiaive Theory of Differenial Equaions 13, No. 3, 1-11; hp://www.mah.u-szeged.hu/ejqde/ Exisence of posiive soluion for a hird-order hree-poin BVP wih sign-changing Green s funcion

More information

CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS

CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS SARAJEVO JOURNAL OF MATHEMATICS Vol.10 (22 (2014, 67 76 DOI: 10.5644/SJM.10.1.09 CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS ALMA OMERSPAHIĆ AND VAHIDIN HADŽIABDIĆ Absrac. This paper presens sufficien

More information

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15. SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a

More information

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange

More information

TO our knowledge, most exciting results on the existence

TO our knowledge, most exciting results on the existence IAENG Inernaional Journal of Applied Mahemaics, 42:, IJAM_42 2 Exisence and Uniqueness of a Periodic Soluion for hird-order Delay Differenial Equaion wih wo Deviaing Argumens A. M. A. Abou-El-Ela, A. I.

More information

Generalized Chebyshev polynomials

Generalized Chebyshev polynomials Generalized Chebyshev polynomials Clemene Cesarano Faculy of Engineering, Inernaional Telemaic Universiy UNINETTUNO Corso Viorio Emanuele II, 39 86 Roma, Ialy email: c.cesarano@unineunouniversiy.ne ABSTRACT

More information

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

On Carlsson type orthogonality and characterization of inner product spaces

On Carlsson type orthogonality and characterization of inner product spaces Filoma 26:4 (212), 859 87 DOI 1.2298/FIL124859K Published by Faculy of Sciences and Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma On Carlsson ype orhogonaliy and characerizaion

More information

arxiv: v1 [math.gm] 7 Nov 2017

arxiv: v1 [math.gm] 7 Nov 2017 A TOUR ON THE MASTER FUNCTION THEOPHILUS AGAMA arxiv:7.0665v [mah.gm] 7 Nov 07 Absrac. In his aer we sudy a funcion defined on naural numbers having eacly wo rime facors. Using his funcion, we esablish

More information

Morning Time: 1 hour 30 minutes Additional materials (enclosed):

Morning Time: 1 hour 30 minutes Additional materials (enclosed): ADVANCED GCE 78/0 MATHEMATICS (MEI) Differenial Equaions THURSDAY JANUARY 008 Morning Time: hour 30 minues Addiional maerials (enclosed): None Addiional maerials (required): Answer Bookle (8 pages) Graph

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

Math Final Exam Solutions

Math Final Exam Solutions Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2.

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2. THE BERNOULLI NUMBERS The Bernoulli numbers are defined here by he exponenial generaing funcion ( e The firs one is easy o compue: (2 and (3 B 0 lim 0 e lim, 0 e ( d B lim 0 d e +e e lim 0 (e 2 lim 0 2(e

More information

Solution of Integro-Differential Equations by Using ELzaki Transform

Solution of Integro-Differential Equations by Using ELzaki Transform Global Journal of Mahemaical Sciences: Theory and Pracical. Volume, Number (), pp. - Inernaional Research Publicaion House hp://www.irphouse.com Soluion of Inegro-Differenial Equaions by Using ELzaki Transform

More information

Short Introduction to Fractional Calculus

Short Introduction to Fractional Calculus . Shor Inroducion o Fracional Calculus Mauro Bologna Deparameno de Física, Faculad de Ciencias Universidad de Tarapacá, Arica, Chile email: mbologna@ua.cl Absrac In he pas few years fracional calculus

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population

1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population Scienific Comuing I Module : Poulaion Modelling Coninuous Models Michael Bader Par I ODE Models Lehrsuhl Informaik V Winer 7/ Discree vs. Coniuous Models d d = F,,...) ) =? discree model: coninuous model:

More information

THE TRACE OF FROBENIUS OF ELLIPTIC CURVES AND THE p-adic GAMMA FUNCTION

THE TRACE OF FROBENIUS OF ELLIPTIC CURVES AND THE p-adic GAMMA FUNCTION THE TRACE OF FROBENIUS OF ELLIPTIC CURVES AND THE -ADIC GAMMA FUNCTION DERMOT McCARTHY Absrac. We define a funcion in erms of quoiens of he -adic gamma funcion which generalizes earlier work of he auhor

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

The Linear Quadratic Regulator on Time Scales

The Linear Quadratic Regulator on Time Scales Inernaional Journal of Difference Equaions ISSN 973-669, Volume 5, Number 2, pp 149 174 (21) hp://campusmsedu/ijde The Linear Quadraic Regulaor on Time Scales Marin Bohner and Nick Winz Missouri Universiy

More information

Ordinary Differential Equations

Ordinary Differential Equations Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

OSCILLATION OF SECOND-ORDER DELAY AND NEUTRAL DELAY DYNAMIC EQUATIONS ON TIME SCALES

OSCILLATION OF SECOND-ORDER DELAY AND NEUTRAL DELAY DYNAMIC EQUATIONS ON TIME SCALES Dynamic Sysems and Applicaions 6 (2007) 345-360 OSCILLATION OF SECOND-ORDER DELAY AND NEUTRAL DELAY DYNAMIC EQUATIONS ON TIME SCALES S. H. SAKER Deparmen of Mahemaics and Saisics, Universiy of Calgary,

More information

Positive continuous solution of a quadratic integral equation of fractional orders

Positive continuous solution of a quadratic integral equation of fractional orders Mah. Sci. Le., No., 9-7 (3) 9 Mahemaical Sciences Leers An Inernaional Journal @ 3 NSP Naural Sciences Publishing Cor. Posiive coninuous soluion of a quadraic inegral equaion of fracional orders A. M.

More information

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation Commun Theor Phys Beijing, China 43 2005 pp 591 596 c Inernaional Academic Publishers Vol 43, No 4, April 15, 2005 An Invariance for 2+1-Eension of Burgers Equaion Formulae o Obain Soluions of KP Equaion

More information

Orthogonal Rational Functions, Associated Rational Functions And Functions Of The Second Kind

Orthogonal Rational Functions, Associated Rational Functions And Functions Of The Second Kind Proceedings of he World Congress on Engineering 2008 Vol II Orhogonal Raional Funcions, Associaed Raional Funcions And Funcions Of The Second Kind Karl Deckers and Adhemar Bulheel Absrac Consider he sequence

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information