Jordan Journal of Mathematics and Statistics (JJMS) 9(2), 2016, pp
|
|
- Beryl Malone
- 5 years ago
- Views:
Transcription
1 Jordan Journal of Mathematics and Statistics JJMS) 92), 2016, pp HOPF BIFURCATION ANALYSIS IN A SYSTEM FOR CANCER VIROTHERAPY WITH EFFECT OF THE IMMUNE SYSTEM A. ASHYANI 1), H.M. MOHAMMADINEJAD 2) AND O. RABIEIMOTLAGH 3) Abstract. We consider a system of differential equations which is motivated biologically and simulates a cancer virotherapy. The existence of equilibrium points and their local stability are studied using the characteristic equation. We investigate Hopf bifurcation around the interior equilibrium point. 1. Introduction Cancer is usually nown by an unnatural growth in cell numbers which is called tumor and sometimes causes to death. Hence, Finding an effective method for controlling the growth rate of tumor always considered as a very important issue in medical science. Many efforts have been done to promote therapeutic ways for cancer. Surgery is the oldest one but it also contains collateral effects so that it is preferable that tumor is treated or controlled without surgery. Cancer virotherapy is the new way which is used recently and it is effective in treatment and it uses virus to damage tumor cells. This method of therapy has been nown in the early twentieth century. However, until 1949, serious experiments about it can not be found [6, 8, 4]. In the first experiments, virus found in the nature were used, but it failed because the immune system response damped the infection. Thus, prevented the virus from destroying the cancer [6, 3, 9]. The next step was done 2000 Mathematics Subject Classification. 34A34, 34D20, 34D23. Key words and phrases. Stability, Asymptotic Stability, Hopf bifurcation. Copyright c Deanship of Research and Graduate Studies, Yarmou University, Irbid, Jordan. Received: Aug. 10, 2015 Accepted: March 9,
2 94 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH by genetic researchers. They produced a new generation of viruses in the laboratory which were capable for cancer virotherapy. The most important fact in the cancer treatment is the fewer non-tumor cells destroyed, the better treatment it could be. The best choice to achieve this goal is cancer virotherapy. These days, such viruses are categorized as oncolytic viruses. There is a ind of iller cells in body which is nown as cytotoxic T lymphocytes CTL). It destroy foreign cells by nowing their special signals that have been sent to immune system. Recognizing cancer cells by immune system tumor-specific CTL starts to destroy them. Although, it may not damage all cancer cells, the cancer virotherapy after injection virus-specific CTL starts to now viruses and try to ill them. Therefore, both viruses and infected cells will be attaced by CTL. On the other hand, it helps to immune system for nowing more tumor cells which have not been recognized before. Then, tumor-specific CTL affects on both infected and non-infected cells to ill them [2, 5]. Experimental observations show that the interaction between the oncolytic virus and tumor might be too complicated. For this reason, it is almost impossible to figure out a complete theory, for finding all relations and analyzing them using the biological observations only. These convince researches to construct mathematical models simulating the virotherapy for studying unnown aspects of the phenomenon. Many of these effective models which are previously developed, based on system of differential equations [13, 12, 11, 7]. In 2001 Wodarz proposed the following system [13, 12]: 1.1) ẋ = rx1 x + y x + y ) dx βxy, ẏ = βxy + sy1 ) ay where yt) and xt) are the number of infected and non-infected cells, respectively. Parameters s and r are the growth rate of infected and non-infected cells in a logistic
3 HOPF BIFURCATION FOR CANCER VIROTHERAPY 95 model, respectively. Size of tumor of infected and non-infected tumor cells are considered to be limited by the carrying capacity. If the total number of infected and non-infected tumor cells becomes more than, the patient will die. Because of this fact we should have x + y <. Parameter a is the death rate of infected cells by virus, and d is the death rate of non-infected cells by the immune system. Finally, β is the rate of viral infection spread in tumor cells. In the same article, Wodarz presented the following system: 1.2) ẋ = rx1 x + y ) dx βxy P txz t, ẏ = βxy + sy1 x + y ) ay P vyz v P t yz t, Z v = C v yz v bz v, Z t = C t yz t x + y) bz t. This system contains four variables. Here, the number of virus-specific CTL and tumor-specific CTL are denoted by Z v and Z t, respectively. The rate of virus-specific CTL proliferate is C v yz v and the rate of illing the infected cells is P v yz v. Virusspecific CTL will die at a rate bz v in the absence of infection. Tumor-specific CTL propagates at a rate C t on infected and non-infected cells x + y). However, The warning signal, which is led to nowing tumor, has been recognized to be more useful in the presence of y virus. Finally, Tumor-specific CTL ills infected and non-infected cells at a rate P t Z t. Wodarz introduced the System 1.2), but he did not present further mathematical analysis or simulation for it. In 2014, we analyzed System 1.1) [1]. Although, we now immune system response is an important factor which effect on virus therapy. Thus, we want to improve System 1.1) by considering immune system response with separate variable, which helped us to provide results that is closer than biology happening. We consider both CTL responses in one variable because some structure
4 96 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH of responses for virus-specific CTL and tumor-specific CTL are similar. In addition, this assumption does not change the system and mae its analysis easier by restricting the system into a third dimensional model. In this paper, we introduce the system for interaction between non-infected cells, infected cells and CTL response as 1.3) 1.4) 1.5) ẋ = rx1 x + y ) βxy p 1xz, ẏ = sy1 x + y ) + βxy p 2yz ay, ż = c 1 xz + c 2 yz bz. The non-infected cells death with response of the immune system is presented by the term p 1 xz. Infected tumor cells die with two process. They die with infection at a rate a and they die with response immune system which is given by the term p 2 yz. On the other hand, CTL for tumor grows is presented by the term c 1 xz and CTL for virus grows by the term c 2 yz. Moreover, in the absence of any infection, the cells will die at a rate b. All the parameters of the system are non negative with initial population conditions x0) 0, y0) 0, z0) 0. This paper is organized as follows. Section 2 contains the positivity and the boundedness of the model. In Section 3, we investigate the conditions for existence at positive equilibrium points and stability of important equilibrium points. We demonstrate the occurrence of Hopf bifurcation in Section 4.In Section 5, we present the numerical results. Finally, we draw some conclusion in Section Positively and boundedness of the solutions The positivity and boundedness of the system ) are necessary to mae the system has biologically meaningful. Therefore, we study them first.
5 HOPF BIFURCATION FOR CANCER VIROTHERAPY 97 Theorem 2.1. All the solutions of system ) which starting in the R 3 + remain positive. Also, they remain bounded in the region Ω R 3 + defined by Ω = {x, y, z) R 3 + x, x + y, x + y + p r + s) z c m }. Proof. Equation 1.3) can be written as dx/x = φ 1 x, y, z)dt, where φ 1 x, y, z) = r1 x+y ) βy p 1z. Then xt) = x0)e t 0 φ1x,y,z)dt 0 for all t 0 as x0) 0. Equation 1.4) gives dy/y = φ 2 x, y, z)dt, where φ 2 x, y, z) = s1 x+y ) + βx p 2z a. Then yt) = y0)e t 0 φ2x,y,z)dt 0 for all t 0 as y0) 0. Similarly, zt) 0 for all t 0. On the other hand, dx rx1 x+y ) rx1 x ). Hence, we consider dt du/dt = ru1 U ), with solving this differential equation we have Ut) = which gives lim sup t Then, lim sup t d dt Ut) =. We now dx/dt du/dt, so lim sup t xt). Simple calculations yield to x + y) = rx1 x + y rx1 x + y 1 e t 0 rt ), xt) lim sup Ut). t ) p 1xz + sy1 x + y ) p 2yz ay ) + sy1 x + y x + y ) δx + y)1 ) where δ = maxr, s). Hence, similarly lim supxt) + yt)) for all t 0. t Let W t) = xt) + yt) + p zt) where p = minp c 1, p 2 ) and c = maxc 1, c 2 ). Then, dw dt = dx dt + dy dt + p dz c + p c c 1xz + p c c 2yz p c dt = rx1 x + y bz rx1 x + y ) p 1xz + sy1 x + y ) p 2yz ay ) + sy1 x + y ) p c bz rx r x2 + sy s y2 p c bz = r r x r) 2 + s s y s) 2 p bz rx sy r + s) mw c
6 98 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH where m = minr, s, b). Therefore, lim supw t)) r + s)/m. Thus, z is bounded t above and all the solutions of the system ) that started in R 3 + either approaches, enter or remain in the subset of Ω = {x, y, z) R 3 + x, x + y, x + y + p r + s) z c m }. Since the solutions of the system ) remain bounded in the positively invariant region Ω, the model is well posed in the epidemiologically and mathematically senses. 3. Existence of equilibrium points and stability analysis In this section, we study the stability of the equilibrium points Equilibrium points. System ) always has two equilibrium points E 0 = 0, 0, 0) and E 1 =, 0, 0). Other equilibrium points are given by E 2 = 0, s a), 0), E 3 = 0, b, bs + c 2s a) ), s c 2 c 2 p 2 E 4 = b, 0, c 1 b)r ar + βa s) ), E 5 = c 1 c 1 p 1 ββ + r s), rβ a) ββ + r s), 0), E 6 = ac 2p 1 b c 2 )p 2 r p 1 s) bp 2 β) c 1 c 2 )p 2 r p 1 s) + βc 1 p 2 c 2 p 1 ),, ac 1 p 1 b c 1 )p 2 r p 1 s) bp 1 β c 1 c 2 )p 2 r p 1 s) + βc 1 p 2 c 2 p 1 ) ac 2 r c 1 r + β)) + β c 2 r + c 1 s + br s + β)) ) := E = x, y, z ). c 1 c 2 )p 2 r p 1 s) + βc 1 p 2 c 2 p 1 ) Equilibrium point E 0 means that the immune system cells, non-infected and infected tumor cells do not exist. Hence, in biology it means that both non-infected and infected cells damaged in the absence of immune system. Equilibrium point E 1 means that only non-infected tumor cells remain and the tumor grows to its carrying capacity, which shows the inefficiency of virotherapy. This equilibrium has been nown as infection-free equilibrium.
7 HOPF BIFURCATION FOR CANCER VIROTHERAPY 99 The equilibrium points are biologically admissible if and only if all components are positive. If s > a, E 2 exists. Equilibrium point E 3 exists if c 2 s a) > bs. Furthermore, E 4 is meaningful if c > b. For existence of equilibrium point E 5 we have the following. Case1 If β + r s > 0, then the denominator of components of equilibrium point E 5 is positive. Thus, for positivity of components of E 5, the numerators should be positive. Therefore, ar + βa s) > 0 and β a > 0. Case2 If β + r s < 0, then the denominator of components of equilibrium point E 5 is negative. Thus, we want that numerators are negative. For this we have ar + βa s) < 0 and β a < 0. Equilibrium E 6 is called interior equilibrium point. It means that immune system cells, non-infected and infected cells exist. This is compatible with what exactly can be seen in reality. Hence, E 6 is the most interesting equilibrium point from the biological point of view. For existence of E 6 we have: Case1 If c 1 c 2 )p 2 r p 1 s)+βc 1 p 2 c 2 p 1 ) > 0, then the denominator of components of equilibrium point E 6 is positive. For existence of E 6, all components of E 6 should be positive. Therefore, all the numerators should be positive. Hence, we gain ac 2 p 1 b c 2 )p 2 r p 1 s) bp 2 β) > 0, ac 1 p 1 b c 1 )p 2 r p 1 s) bp 1 β > 0, ac 2 r c 1 r + β)) + β c 2 r + c 1 s + br s + β)) > 0. Case2 If c 1 c 2 )p 2 r p 1 s)+βc 1 p 2 c 2 p 1 ) < 0, then the denominator of components of equilibrium point E 6 is negative. Thus, all the numerators should be negative.
8 100 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH Hence, we get ac 2 p 1 b c 2 )p 2 r p 1 s) bp 2 β) < 0, ac 1 p 1 b c 1 )p 2 r p 1 s) bp 1 β < 0, ac 2 r c 1 r + β)) + β c 2 r + c 1 s + br s + β)) < Stability analysis. We compute the variational matrix of the system ) at any point x, y, z). Real part of eigenvalues of this matrix at any point determine the local stability of equilibrium points. Variational matrix is given by J = J ij ) 3 3, where J 11 = r1 x + y)/) rx)/ p 1 z yβ, J 12 = rx)/ xβ, J 13 = p 1 x, J 21 = sy)/ + yβ, J 22 = s1 x + y)/) sy)/ a p 2 z + xβ, J 23 = p 2 y, J 31 = c 1 z, J 32 = c 2 z, J 33 = c 1 x + c 2 y b. At the equilibrium point E 0, the eigenvalues of J are r, b and s a. Since r is positive, equilibrium point E 0 is always unstable. As a result, tumor will grow and therapy is failed in the absence of immune response and virus. In the following theorem, we study the local stability of the infection-free equilibrium. Theorem 3.1. The infection-free equilibrium of the system ) is unstable for > a/β) or > b/c 1 ) and locally asymptotically stable for < a/β) and < b/c 1 ). Proof. The variational matrix of the system ) at equilibrium point E 1 is given by r r β p 1 JE 1 ) = 0 β a c 1 b
9 HOPF BIFURCATION FOR CANCER VIROTHERAPY 101 Therefore, the eigenvalues are r, β a, and c 1 b. Because r is negative, stability of E 1 depends on the signs of β a and c 1 b. When > a/β) or > b/c 1 ), equilibrium point E 1 is unstable while it is asymptotically stable for < a/β) and < b/c 1 ). For local stability of interior equilibrium point we now that the characteristic equation of the variational matrix J at E 6 can be written as 3.1) DetJ λi 3 ) = λ 3 + a 2 λ 2 + a 1 λ + a 0 = 0 where a 0 = A 3 A 5 A 7 A 2 A 6 A 7 A 3 A 4 A 8 + A 1 A 6 A 8, 3.2) a 1 = A 2 A 4 A 1 A 5 + A 3 A 7 + A 6 A 8 ), a 2 = A 1 + A 5 ), where A 1 A 8 ) are introduced below. 3.3) A 1 = rac 2p 1 b c 2 )p 2 r p 1 s) bp 2 β) c 1 c 2 )p 2 r p 1 s) + c 2 p 1 + c 1 p 2 )β), A 2 = r + β)ac 2p 1 b c 2 )p 2 r p 1 s) bp 2 β), c 1 c 2 )p 2 r p 1 s) + c 2 p 1 + c 1 p 2 )β) A 3 = p 1 ac 2 p 1 + b c 2 )p 2 r p 1 s) + bp 2 β) c 1 + c 2 )p 2 r p 1 s) + c 2 p 1 c 1 p 2 )β, A 4 = s β) ac 1p 1 + b c 1 )p 2 r p 1 s) + bp 1 β), c 1 c 2 )p 2 r p 1 s) + c 2 p 1 + c 1 p 2 )β) A 5 = s ac 1p 1 + b c 1 )p 2 r p 1 s) + bp 1 β) c 1 c 2 )p 2 r p 1 s) + c 2 p 1 + c 1 p 2 )β), A 6 = p 2 ac 1 p 1 + b c 1 )p 2 r p 1 s) + bp 1 β) c 1 c 2 )p 2 r p 1 s) + c 2 p 1 + c 1 p 2 )β, A 7 = c 1ac 2 r c 1 r + β)) + β c 2 r + c 1 s + br s + β))), c 1 c 2 )p 2 r p 1 s) + c 2 p 1 + c 1 p 2 )β A 8 = c 2ac 2 r c 1 r + β)) + β c 2 r + c 1 s + br s + β))). c 1 c 2 )p 2 r p 1 s) + c 2 p 1 + c 1 p 2 )β
10 102 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH Hence, in the following theorem we present some conditions for local stability of interior equilibrium point E 6 = x, y, z ). Theorem 3.2. Let a 0, a 1 and a 2 be defined in 3.2). The interior equilibrium point of the system ) is locally asymptotically stable if a 0 > 0 and a 2 a 1 a 0 > 0. Proof. According to Routh-Herwitz criterion of stability, the necessary and the sufficient conditions to all roots of characteristic equation 3.1) to be negative real parts are a 2 > 0, a 0 > 0, and a 2 a 1 a 0 > 0. It is clear that a 2 > 0. Thus, the interior equilibrium point E 6 is locally asymptotically stable if a 0 > 0 and a 2 a 1 a 0 > Hopf bifurcation In this section, we find the conditions for existence of Hopf bifurcation at the interior equilibrium point E 6. First we linearized equation ) in E 6. Let u 1 t) = xt) x, u 2 t) = yt) y, and u 3 t) = zt) z. It leads us to E 6 to be the origin. Then ) can be written as 4.1) u 1 = A 1 u 1 + A 2 u 2 + A 3 u 3 u 1 ru 1 + u 2 ) + p 1 u 3 + u 2 β), u 2 = A 4 u 1 + A 5 u 2 + A 6 u 3 u 2 su 1 + u 2 ) + p 2 u 3 u 1 β), u 3 = A 7 u 1 + A 8 u 2 + c 1 u 1 + c 2 u 2 )u 3 where A 1 A 8 ) are introduced in 3.3). The linearization of 4.1) at 0, 0) is u 1 = A 1 u 1 + A 2 u 2 + A 3 u 3, u 2 = A 4 u 1 + A 5 u 2 + A 6 u 3, u 3 = A 7 u 1 + A 8 u 2 with the characteristic equation λ 3 + a 2 λ 2 + a 1 λ + a 0 = 0, where a 2 = A 1 + A 5 ), a 1 = A 2 A 4 A 1 A 5 +A 3 A 7 +A 6 A 8 ), and a 0 = A 3 A 5 A 7 A 2 A 6 A 7 A 3 A 4 A 8 +A 1 A 6 A 8. If a 0 = a 1 a 2, we have a pair of imaginary eigenvalues λ 1,2 = ±i a 1. In this case, λ 3 = a 2 is the real eigenvalue. Thus let ɛ = a 0 a 1 a 2 and using the Implicit Function Theorem, we have λ 1,2 = ɛ ± i 2a 1 +a 2 2 a ) 1 + ɛ a1 a 2 2a 1 a 1 +a 2 2 ) ) + oɛ2 ) and
11 HOPF BIFURCATION FOR CANCER VIROTHERAPY 103 λ 3 = a 2 ɛ a 1 +a oɛ 2 ). Furthermore, if ɛ = 0, we may have Hopf bifurcation. On the other hand, solving a 2, a 1 and a 0 with respect to A 1, A 2 and A 3 respectively, we get A 1 = A 5 a 2, A 2 = A A 3 A 7 + A 6 A 8 + a 1 + A 5 a 2, and A 4 A 3 = A 5 2 A 6 A 7 A 4 A 5 A 6 A 8 + A 2 6 A 7 A 8 A 4 a 0 + A 6 A 7 a 1 + A 5 A 6 A 7 a 2 A 4 A 5 A 7 + A 6 A 2 7 A 2 4 A 8 + A 4 A 6 A 8 a 2 A 4 A 5 A 7 + A 6 A 7 2 A 4 2 A 8. In this case, we have 2-dimensional center manifold and afterwards by projecting this system on the center manifold we gain 2-dimensional system. First, we use the below transformation to get the standard system; where U 1 U 2 U 3 = P 1 u 1 u 2 u 3 P 1 P 2 P 3, and P = P 4 P 5 P P 1 = A 5A 6 A 7 A 8 A 4 A 6 A 2 8 A 4 A 8 a 1, P A 5 A 7 + A 4 A 8 ) 2 + A 2 4 = A 5A 6 A 2 7 A 4 A 6 A 7 A 8 A 4 A 7 a 1, 7a 1 A 5 A 7 + A 4 A 8 ) 2 + A 2 7a 1 P 3 = A 6A 8 + A 5 a 2 + a 2 2 a1 A 2, P 2 = 5A 7 A 4 A 5 A 8 + A 6 A 7 A 8 + A 7 a 1 ), A 5 A 7 A 4 A 8 + A 7 a 2 A 5 A 7 + A 4 A 8 ) 2 + A 2 7a 1 a1 A 4 A 5 A 7 A 6 A A 2 P 5 = 4A 8 ) A 6 A 7 A 4 a 2, P A 5 A 7 + A 4 A 8 ) 2 + A 2 6 =. 7a 1 A 5 A 7 A 4 A 8 + A 7 a 2 Therefore, system 4.1) transforms into the following standard form; U 1 = a 1 U 2 + G 1 U 1, U 2, U 3 ), 4.2) U 2 = a 1 U 1 + G 2 U 1, U 2, U 3 ), U 3 = a 2 U 3 + G 3 U 1, U 2, U 3 )
12 104 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH where G 1 U 1, U 2, U 3 ) = 1 P 1 P 3 )P 5 + P 2 P 6 P 4 ) ) P 2 P 6 P 3 P 5 )U 1 + U 3 ) c 1 P 1 U 1 + P 2 U 2 + P 3 U 3 ) + c 2 P 4 U 1 + P 5 U 2 + P 6 U 3 )) ) + P 2 P 4 U 1 + P 5 U 2 + P 6 U 3 )p 2 U 1 + U 3 ) + sp 1 + P 4 )U 1 + P 2 + P 5 )U 2 + P 3 + P 6 )U 3 ) P 1 U 1 + P 2 U 2 + P 3 U 3 )β) P 5 P 1 U 1 + P 2 U 2 + P 3 U 3 )p 1 U 1 + U 3 ) + rp 1 + P 4 )U 1 ) + P 2 + P 5 )U 2 + P 3 + P 6 )U 3 ) + P 4 U 1 + P 5 U 2 + P 6 U 3 )β), G 2 U 1, U 2, U 3 ) = 1 P 1 P 3 )P 5 + P 2 P 6 P 4 ) ) P 3 P 4 P 1 P 6 )U 1 + U 3 ) c 1 P 1 U 1 + P 2 U 2 + P 3 U 3 ) + c 2 P 4 U 1 + P 5 U 2 + P 6 U 3 )) ) + P 1 P 3 )P 4 U 1 + P 5 U 2 + P 6 U 3 ) p 2 U 1 + U 3 ) sp 1 + P 4 )U 1 + P 2 + P 5 )U 2 + P 3 + P 6 )U 3 ) + P 1 U 1 + P 2 U 2 + P 3 U 3 )β) + P 4 P 6 )P 1 U 1 + P 2 U 2 + P 3 U 3 ) p 1 U 1 + U 3 ) + rp 1 + P 4 )U 1 ) + P 2 + P 5 )U 2 + P 3 + P 6 )U 3 ) + P 4 U 1 + P 5 U 2 + P 6 U 3 )β), G 3 U 1, U 2, U 3 ) = 1 P 1 P 3 )P 5 + P 2 P 6 P 4 ) ) P 1 P 5 P 2 P 4 )U 1 + U 3 ) c 1 P 1 U 1 + P 2 U 2 + P 3 U 3 ) + c 2 P 4 U 1 + P 5 U 2 + P 6 U 3 )) ) + P 5 P 1 U 1 + P 2 U 2 + P 3 U 3 )p 1 U 1 + U 3 ) + rp 1 + P 4 )U 1 + P 2 + P 5 )U 2 + P 3 + P 6 )U 3 ) + P 4 U 1 + P 5 U 2 + P 6 U 3 )β) P 2 P 4 U 1 + P 5 U 2 + P 6 U 3 )p 2 U 1 + U 3 ) + sp 1 + P 4 )U 1 ) + P 2 + P 5 )U 2 + P 3 + P 6 )U 3 ) P 1 U 1 + P 2 U 2 + P 3 U 3 )β).
13 HOPF BIFURCATION FOR CANCER VIROTHERAPY 105 By the existence theorem in the center manifold theory [10], there exists a center manifold which can be expressed locally as W c 0) = {U 1, U 2, U 3 ) R 3 U 3 = hu 1, U 2 ), U 3 < δ, h0, 0) = 0, Dh0, 0) = 0} for δ sufficiently small. We now compute W c 0). Assume that hu 1, U 2 ) have the following form: 4.3) U 3 = hu 1, U 2 ) = b 1 U b 2 U b 3 U 1 U 2 + ou 1, U 2 ). Based on the invariance of W c 0) under the dynamics of 4.2), the center manifold must satisfy 4.4) N U 1, U 2 ) = a 2 hu 1, U 2 ) + G 3 U 1, U 2, hu 1, U 2 )) h U 1, U 2 ) a 1 U 2 + G 1 U 1, U 2, hu 1, U 2 )) ) U 1 h U 2 U 1, U 2 ) a1 U 1 + G 2 U 1, U 2, hu 1, U 2 )) ).
14 106 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH Substituting 4.3) into 4.4) and equating terms of lie powers to zero, we obtain b 1 = 1 2a 2 2b 3 a1 + 2 A 8 c 1 + A 7 c 2 )A 2 5A 2 7 2A 4 A 5 A 7 A 8 + A 2 4A A 2 7a 1 ) A 5 A 6 A 7 A 4 A 6 A 8 + a 1 )) 2 + 1/)A 8 A 4 A 5 A 7 A 6 A A 2 4A 8 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a + 1))A 2 5A 2 7p 1 + A 2 4A 2 8p 1 A 5 A 7 2A 4 A 8 p 1 + A 6 A 7 r A 8 r + A 7 β)) + A 2 7p 1 a 1 + A 4 A 7 r A 8 r + A 7 β)a 6 A 8 + a 1 )) 1/)A 7 A 5 A 6 A 7 A 4 A 6 A 8 + a 1 ))A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 )) A 2 5A 2 7p 2 + A 2 4A 2 8p 2 A 5 A 7 2A 4 A 8 p 2 + A 6 A 7 s A 8 s + A 8 β)) + A 2 7p 2 a 1 + A 4 A 7 s A 8 s + A 8 β)a 6 A 8 + a 1 )))A 5 A 7 A 4 A 8 + A 7 a 2 ) )/ A 4 A 5 A 7 A 6 A A 2 4A 8 )A 2 5A 2 7 2A 4 A 5 A 7 A 8 + A 2 4A 2 8 ) + A 2 7a 1 ) 2 a 1 + a 2 2)), b 2 = 1/2 a 1 )) b 3 a 2 + a1 A 4 A 8 + A 7 A 5 + a 2 ))A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) A 2 5A 7 c 1 + A 2 4A 8 c 2 + A 6 A 7 A 8 c 1 A 7 c 2 ) A 4 A 5 A 8 c 1 + A 7 c 2 ) + A 7 c 1 a 1 ) A 5 A 6 A 7 A 4 A 6 A 8 + a 1 )) A 7 A 5 A 6 A 7 A 4 A 6 A 8 + a 1 ))A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 )) A 2 4A 8 A 4 A 5 A 7 + A 8 ) + A 7 A 2 5 A 6 A 7 + A 6 A 8 + a 1 ))s + A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))β))/ A 8 A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 ))A 2 4A 8 r + β) A 4 A 5 A 7 + A 8 )r + A 7 β) + A 7 A 2 5 A 6 A 7 + A 6 A 8 + a 1 )r A 6 A 7 β)))/ A 4 A 5 A 7 + A 6 A 2 7 A 2 4A 8 )A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))A 2 5A 2 7p 1 + A 2 4A 2 8p 1 + A 2 7a 1 p 1 + A 4 A 6 A 8 + a 1 )A 7 r A 8 r + A 7 β) A 5 A 7 2A 4 A 8 p 1 + A 6 A 7 r A 8 r + A 7 β))))/ + A 4 A 5 A 7 + A 6 A 2 7 A 2 4A 8 )A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8
15 HOPF BIFURCATION FOR CANCER VIROTHERAPY a 1 ))A 2 5A 2 7p 2 + A 2 4A 2 8p 2 + A 2 7a 1 p 2 + A 4 A 6 A 8 + a 1 )A 7 s A 8 s + A 8 β) A 5 A 7 2A 4 A 8 p 2 + A 6 A 7 s A 8 s + A 8 β))))/) )/ A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 )A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) 2 a 1 + a 2 2)) + 1/a 2 a1 2b3 a1 + 2 A 4 A 8 + A 7 A 5 + a 2 )) A 8 c 1 + A 7 c 2 )A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 )A 5 A 6 A 7 A 4 A 6 A 8 + a 1 )) 2 + A 8 A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 )) A 2 5A 2 7p 1 + A 2 4A 2 8p 1 + A 2 7a 1 p 1 + A 4 A 6 A 8 + a 1 )A 7 r A 8 r + A 7 β) A 5 A 7 2A 4 A 8 p 1 + A 6 A 7 r A 8 r + A 7 β))))/ 1/A7A 5 A 6 A 7 A 4 A 6 A 8 + a 1 ))A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))A 2 5A 2 7p 2 + A 2 4A 2 8p 2 + A 2 7a 1 p 2 + A 4 A 6 A 8 + a 1 )A 7 s A 8 s + A 8 β) A 5 A 7 2A 4 A 8 p 2 + A 6 A 7 s A 8 s / + A 8 β))))) A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 )A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) 2 a 1 + a 2 2)) )), b 3 = a1 A 5 A 7 A 4 A 8 + A 7 a 2 ) 2/)r s + β)a 1 A 3 4A 5 A A 2 7A A 6 A 8 + a 1 )A A 6 A 7 + A 8 ) + a 1 ) A 4 A 5 A 7 2A 6 A A 2 5A 7 + 2A 8 ) + A 7 + 2A 8 )a 1 ) + A 2 4A 8 A 2 52A 7 + A 8 ) + A 7 A 6 A 8 + a 1 ))) 2 A 4 A 5 A 7 A 6 A A 2 4A 8 A 8 c 1 + A 7 c 2 )A 2 5A 2 7 2A 4 A 5 A 7 A 8 + A 2 4A A 2 7a 1 )A 5 A 6 A 7 A 4 A 6 A 8 + a 1 )) 2 + A 8 A 4A 5 A 7 A 6 A A 2 4A 8 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 ))A 2 5A 2 7p 1 + A 2 4A 2 8p 1 A 5 A 7 2A 4 A 8 p 1 + A 6 A 7 r A 8 r + A 7 β)) + A 2 7p 1 a 1 + A 4 A 7 r A 8 r + A 7 β)a 6 A 8 + a 1 )) A 7 A 5A 6 A 7 A 4 A 6 A 8 + a 1 ))A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))A 2 5A 2 7p 2 + A 2 4A 2 8p 2 A 5 A 7 2A 4 A 8 p 2 + A 6 A 7 s A 8 s + A 8 β)) + A 2 7p 2 a 1 + A 4 A 7 s A 8 s + A 8 β)a 6 A 8 + a 1 ))) 1 A 4 A 5 A 7 A 6 A A 2 4A 8 A 2 5A 2 7 2A 4 A 5 A 7 A 8 + A 2 4A A 2 7a 1 )A 2 5A 7 c 1 + A 2 4A 8 c 2 + A 6 A 7 A 8 c 1 A 7 c 2 )
16 108 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH A 4 A 5 A 8 c 1 + A 7 c 2 ) + A 7 c 1 a 1 )A 5 A 6 A 7 A 4 A 6 A 8 + a 1 )) A 4A 5 A 7 A 6 A A 2 4A 8 A 2 5A 7 + A 4 A 5 A 8 A 7 A 6 A 8 + a 1 ))A 2 5A 2 7p 1 + A 2 4A 2 8p 1 A 5 A 7 2A 4 A 8 p 1 + A 6 A 7 r A 8 r + A 7 β)) + A 2 7p 1 a 1 + A 4 A 7 r A 8 r + A 7 β)a 6 A 8 + a 1 )) + A 4A 5 A 7 A 6 A A 2 4A 8 A 2 5A 7 + A 4 A 5 A 8 A 7 A 6 A 8 + a 1 ))A 2 5A 2 7p 2 + A 2 4A 2 8p 2 A 5 A 7 2A 4 A 8 p 2 + A 6 A 7 s A 8 s + A 8 β)) + A 2 7p 2 a 1 + A 4 A 7 s A 8 s + A 8 β)a 6 A 8 + a 1 )) A 8 A 4 A 5 A 7 A 6 A A 2 4A 8 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 ))A 2 4A 8 r + β) A 4 A 5 A 8 r + A 7 r + β)) + A 7 A 2 5r A 6 A 7 r + A 6 A 8 r A 6 A 7 β + ra 1 )) A 7 /)A 5 A 6 A 7 A 4 A 6 A 8 + a 1 ))A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 )) A 2 4A 8 s + A 4 A 5 A 7 s + A 8 s β)) + A 7 A 6 A 7 s A 6 A 8 s + A 6 A 8 β + A 2 5 s + β) sa 1 + βa 1 )))a 2 ) )/ A 2 5A 2 7 2A 4 A 5 A 7 A 8 + A 2 4A A 2 7a 1 ) 2 a 1 + a 2 2)4a 1 + a 2 2)) Finally, by substituting 4.3) into 4.2) we obtain the vector field reduced to the center manifold by U 1 = a 1 U 2 + η 1 U 1 U 2 + η 2 U η 3 U 2 2, 4.5) U 2 = a 1 U 1 + θ 1 U 1 U 2 + θ 2 U θ 3 U 2 2 where a1 1 η 1 = A 4 A 8 + A 7 A 5 + a 2 )) A 4 A 8 + A 7 A 5 + a 2 ) A 5A 7 A 4 A 8 ) 2 + A 7 c 1 a 1 ) + A 2 7a 1 )A 2 5A 7 c 1 + A 2 4A 8 c 2 + A 6 A 7 A 8 c 1 A 7 c 2 ) A 4 A 5 A 8 c 1 + A 7 c 2 )A 6 A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) + A 7 A 5 A 6 A 7 A 4 A 6 A 8 + a 1 ))a 2 + A 4 A 5 A 7 + A 6 A 2 7 A 2 4A 8 )a 2 2) + 1 A7 A 5 A 6 A 7 A 4 A 6 A 8 + a 1 ))A 2 5A 7
17 HOPF BIFURCATION FOR CANCER VIROTHERAPY 109 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 )) A 2 4A 8 A 4 A 5 A 7 + A 8 ) + A 7 A 2 5 A 6 A 7 + A 6 A 8 + a 1 ))s + A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))β) + A 8 A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 ))A 2 4A 8 r + β) A 4 A 5 A 7 + A 8 )r + A 7 β) + A 7 A 2 5 A 6 A 7 + A 6 A 8 + a 1 )r A 6 A 7 β)) + A 4 A 5 A 7 + A 6 A 2 7 A 2 4A 8 )A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))A 2 5A 2 7p 1 + A 2 4A 2 8p 1 + A 2 7a 1 p 1 + A 4 A 6 A 8 + a 1 )A 7 r A 8 r + A 7 β) A 5 A 7 2A 4 A 8 p 1 + A 6 A 7 r A 8 r + A 7 β))) A 4 A 5 A 7 + A 6 A 2 7 A 2 4A 8 )A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))A 2 5A 2 7p 2 + A 2 4A 2 8p 2 + A 2 7a 1 p 2 + A 4 A 6 A 8 + a 1 )A 7 s A 8 s + A 8 β) A 5 A 7 2A 4 A 8 p 2 + A 6 A 7 s A 8 s + A 8 β))) )))/ 2 A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 )A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) 2 a 1 + a 2 2)), η 2 = 1 A 4 A 8 + A 7 A 5 + a 2 )) A8 c 1 A 7 c 2 )A 5 A 7 A 4 A 8 A 7 A 5 + a 2 ) A 4 A 8 ) 2 + A 2 7a 1 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 ))A 6 A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) + A 7 A 5 A 6 A 7 A 4 A 6 A 8 + a 1 ))a 2 + A 4 A 5 A 7 + A 6 A 2 7 A 2 4A 8 )a 2 2) ) 1 A8 A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 ))A 2 5A 2 7p 1 + A 2 4A 2 8p 1 + A 2 7a 1 p 1 + A 4 A 6 A 8 + a 1 )A 7 r A 8 r + A 7 β) A 5 A 7 2A 4 A 8 p 1 + A 6 A 7 r A 8 r + A 7 β))) + A 7 A 5 A 6 A 7 A 4 A 6 A 8 + a 1 ))A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))A 2 5A 2 7p 2 + A 2 4A 2 8p 2 + A 2 7a 1 p 2 + A 4 A 6 A 8 + a 1 )A 7 s A 8 s + A 8 β) A 5 A 7 2A 4 A 8 p 2 + A 6 A 7 s A 8 s + A 8 β))) )))/ A 7 A 4 A 5 A 6 A 7 ) + A 2 4A 8 )A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) 2 a 1 + a 2 2)), η 3 = a 1 A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))A 2 4A 8 A 4 A 5 A 7 + A 8 ) )/ + A 7 A 2 5 A 6 A 7 + A 6 A 8 + a 1 )) A 4 A 8 + A 7 A 5 + a 2 ))r s + β) A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) 2 a 1 + a 2 2)),
18 110 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH θ 1 = A 4 A 8 + A 7 A 5 + a 2 )) A A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 )A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) 2 a 1 + a 2 5 A 7 2) A 4 A 8 ) 2 + A 2 7a 1 )A 2 5A 7 c 1 + A 2 4A 8 c 2 + A 6 A 7 A 8 c 1 A 7 c 2 ) A 4 A 5 A 8 c A 7 c 2 ) + A 7 c 1 a 1 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 ))a 2 + A 4 A 8 + A 7 A 5 + a 2 )) A 7 A 5 A 6 A 7 A 4 A 6 A 8 + a 1 )) A 4 A 7 A 8 A 5 a 1 + 2A A 6 A 8 + a 1 )a 2 + 2A 5 a 2 2) + A 2 4A 2 8a 1 + a 2 A 5 + a 2 )) + A 2 7A a 1 )a 2 A 5 + a 2 ) + A 6 A 8 a 1 + A 5 a 2 ))) A 2 4A 8 A 4 A 5 A 7 + A 8 ) + A 7 A 2 5 A 6 A 7 + A 6 A 8 + a 1 ))s + A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))β)) A 8 A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 ))A 4 A 8 a 2 + A 7 a 1 A 5 a 2 ))A 2 4A 8 r + β) A 4 A 5 A 7 + A 8 )r + A 7 β) + A 7 A 2 5 A 6 A 7 + A 6 A 8 + a 1 )r A 6 A 7 β))) + A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 )A 2 5A 7 A 4 A 5 A 8 + A 7 A 6 A 8 + a 1 ))A 4 A 8 a 2 + A 7 a 1 A 5 a 2 ))A 2 5A 2 7p 1 + A 2 4A 2 8p 1 + A 2 7a 1 p 1 + A 4 A 6 A 8 + a 1 )A 7 r A 8 r + A 7 β) A 5 A 7 2A 4 A 8 p 1 + A 6 A 7 r A 8 r + A 7 β)))) + A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 ) A 4 A 7 A 8 A 5 a 1 + 2A A 6 A 8 + a 1 )a 2 + 2A 5 a 2 2) + A 2 4A 2 8a 1 + a 2 A 5 + a 2 )) + A 2 7A a 1 )a 2 A 5 + a 2 ) + A 6 A 8 a 1 + A 5 a 2 )))A 2 5A 2 7p 2 + A 2 4A 2 8p 2 + A 2 7a 1 p 2 + A 4 A 6 A 8 + a 1 )A 7 s A 8 s + A 8 β) A 5 A 7 2A 4 A 8 p 2 + A 6 A 7 s )) A 8 s + A 8 β))), θ 2 = A 4 A 8 + A 7 A 5 + a 2 )) A 8 c 1 A 7 c 2 )A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 )A 5 A 6 A 7 A 4 A 6 A 8 + a 1 )) 2 1 a 2 + A 8 A 7 A 4 A 5 + A 6 A 7 ) A 4 A 8 + A 7 A 5 + a 2 )) + A 2 4A 8 ) A 5 A 6 A 7 + A 4 A 6 A 8 + a 1 ))A 4 A 8 a 2 + A 7 a 1 A 5 a 2 ))A 2 5A 2 7p 1 + A 2 4A 2 8p 1 + A 2 7a 1 p 1 + A 4 A 6 A 8 + a 1 )A 7 r A 8 r + A 7 β) A 5 A 7 2A 4 A 8 p 1 + A 6 A 7 r A 8 r + A 7 β))) + A 7 A 5 A 6 A 7 A 4 A 6 A 8 + a 1 )) A 4 A 7 A 8 A 5 a 1 + 2A A 6 A 8 + a 1 )a 2 + 2A 5 a 2 2) + A 2 4A 2 8a 1 + a 2 A 5 + a 2 )) + A 2 7A a 1 )a 2
19 HOPF BIFURCATION FOR CANCER VIROTHERAPY 111 θ 3 = A 5 + a 2 ) + A 6 A 8 a 1 + A 5 a 2 )))A 2 5A 2 7p 2 + A 2 4A 2 8p 2 + A 2 7a 1 p 2 + A 4 A 6 A 8 ))/ + a 1 )A 7 s A 8 s + A 8 β) A 5 A 7 2A 4 A 8 p 2 + A 6 A 7 s A 8 s + A 8 β))) A 7 A 4 A 5 + A 6 A 7 ) + A 2 4A 8 ) a 1 A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) 2 a 1 + a 2 2)), a 1 A 2 A 5 A 7 A 4 A 8 ) 2 + A 2 7a 1 ) 2 a 1 + a 5A A 4 A 5 A 8 A 7 A 6 A 8 2)) + a 1 ))A 4 A 8 a 2 + A 7 a 1 A 5 a 2 ))A 2 4A 8 r + β) A 4 A 5 A 7 + A 8 )r + A 7 β) + A 7 A 2 5 A 6 A 7 + A 6 A 8 + a 1 )r A 6 A 7 β)) A 4 A 7 A 8 A 5 a 1 + 2A A 6 A 8 + a 1 )a 2 + 2A 5 a 2 2) + A 2 4A 2 8a 1 + a 2 A 5 + a 2 )) + A 2 7A a 1 )a 2 A 5 + a 2 ) + A 6 A 8 a 1 + A 5 a 2 )))A 2 4A 8 s + A 7 A 2 5 A 6 A 7 + A 6 A 8 + a 1 )s A 7 A A 6 A 8 ) + a 1 )β + A 4 A 5 A 7 + A 8 )s + A 8 β)). At the bifurcation point, System 4.5) becomes where U 1 = 0 a 1 U 1 + f 1 U 1, U 2 ) U 2 a1 0 f 2 U 1, U 2 ) U 2 f 1 U 1, U 2 ) = η 1 U 1 U 2 + η 2 U η 3 U 2 2 and f 2 U 1, U 2 ) = θ 1 U 1 U 2 + θ 2 U θ 3 U 2 2. The coefficient a0) a is given by a = 1 16 f 1 U1 U 1 U 1 + f 1 U 1 U 2 U 2 + f 2 U 1 U 1 U 2 + f 2 U 2 U 2 U 2 ) a 1 f 1 U1 U 2 f 1 U 1 U 1 + f 1 U 2 U 2 ) f 2 U 1 U 2 f 2 U 1 U 1 + f 2 U 2 U 2 ) f 1 U 1 U 1 f 2 U 1 U 1 + f 1 U 2 U 2 f 2 U 2 U 2 ). Therefore, a = 1 8 a 1 η1 η 2 + η 3 ) θ 1 θ 2 + θ 3 ) 2η 2 θ 2 + 2η 3 θ 3 ) 0. Hence, from Hopf bifurcation theory in [10] we get Hopf bifurcation. Also, the periodic orbit is asymptotically stable for a < 0 and unstable for a > 0.
20 112 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH 5. Numerical simulation In this section, we carry out some numerical simulations to verify the previous results. Suppose that r = 1.1, β = 1, s = 1, = 0.7, a = , 5.1) b = 0.1, p 1 = 1, p 2 = 2, c 1 = 2, c 2 = 0.1. In this case, E 6 = , , ) and System ) is given as ẋ = 1.1x1 x + y ) xy xz, 0.7 ẏ = y1 x + y ) + xy 2yz y, 0.7 ż = 2xz + 0.1yz 0.1z. Furthermore, restricted system to the center manifold is approximately given by U 1 = U U 1 U U U 2 2, U 2 = U U 1 U U U 2 2. Thus a 0 = a 1 a , λ 1,2 ±i and λ 3 = In this case a = Therefore, we have Hopf bifurcation that bifurcating periodic orbit is stable. Thus, E 6 is unstable. These results are presented in Figures 1 and 2).
21 HOPF BIFURCATION FOR CANCER VIROTHERAPY 113 Figure 1. Components of the solution U 1 t) and U 2 t) with parameters given in 5.1). Figure 2. The bifurcating periodic solution. 6. Conclusion In this paper, we consider a virus therapy for cancer and analyze stability of equilibrium points of this system. Equilibrium E 0 is unstable. It means that with this therapy, tumor does not completely destroy. Thus, we loo for on other equilibrium points. Infection-free equilibrium point E 1 exists and may be stable. Stability of E 1 is not useful because it means that non-infected tumor cells exist then therapy fails. We note that the interior equilibrium point E 6 is the most interesting equilibrium point from the biological point of view since its existence means that both of the
22 114 A. ASHYANI, H.M. MOHAMMADINEJAD AND O. RABIEIMOTLAGH non-infected and infected tumor cells and CTL exist. Its stability means that the tumor growth is controlled in a way that it can not reach to the carrying capacity. Hence, under the conditions of the parameters in theorem 3.2, it means that with this therapy we could control the size of the tumor which is x + y, but tumor exists and not completely wasted. In addition, we noted that when a 0 = a 1 a 2, we have Hopf bifurcation in E 6. Moreover, theorem 3.2 means this equilibrium point is unstable. Hence, existence of Hopf bifurcation is not useful because Hopf cycle means that the size of tumor cells decreases. Exactly when we suppose that therapy is useful, tumor size increases, and this behavior repeated See Figure 2). Although, this behavior means when we suppose that tumor treated, after some times we see that tumor grows. Thus, patient will die. Therefore, we attempted that value of parameters must be controlled in the way that we have not Hopf cycle. References [1] A. Ashyani, H. M. Mohammadinejad, O. RabieiMotlagh Stability analysis of mathematical model of virus therapy for cancer.to be appear in Iranian Journal of Mathematical Sciences and Informatics 2016). [2] E. J. Fuchs, P. Matzinger. Is cancer dangerous to the immune system? Semin Immuno. 18, 1996), , [3] R. J. Huebner, W. P. Rowe, W. E. Schatten, R. R. Smith, L. B. Thomas. Studies on the use of viruses in the treatment of carcinoma of the cervix. Cancer. 9, 1956), [4] C. M. Kunin. Cellular susceptibility to enteroviruses. Bacteriological Reviews. 28, No. 4, 1964), [5] P. Matzinger, An innate sense of danger. Semin Immunol. 10, 1998), , [6] A.E. Moore. The destructive effect of the virus of Russian far east encephalitis on the transplantable mouse sarcoma 180. Cancer. 2, 1949), [7] A. S. Novozhilov, F. S. Berezovsaya, E. V. Koonin, G. P. Karev. Mathematical modeling of tumor therapy with oncolytic viruses: regimes with complete tumor elimination within the framewor of deterministic models. Biology Direct. 1, 2006), 6 24.
23 HOPF BIFURCATION FOR CANCER VIROTHERAPY 115 [8] L. Pelner, GA. Fowler, HC. Nauts. Effects of concurrent infections and their toxins on the course of leuemia. Acta. Med. Scand. Suppl. 338, 1958), [9] A. R. Pond, E. E. Manuelidis. Oncolytic effect of poliomyelitis virus on human epidermoid carcinoma Hela tumor) heterologously transplanted to guinea pigs. The American J. of Pathology. 45, No. 2, 1964), [10] S. Wiggins: Introduction to applied nonlinear dynamical systems and chaos, Berlin, Springer, [11] D. Wodarz, N. L. Komarova. Computational biology of cancer: lecture notes and mathematical modeling. World Scientific Publishing Company, Singapore, [12] D. Wodarz. Computational approaches to study oncolytic virus therapy: insights and challenges. Gene. Ther. Mol. Biol. 8, 2004), [13] D. Wodarz. Viruses as antitumor weapons: defining conditions for tumor remission. Cancer Research. 61, 2001), ) Department of Mathematics, University of Birjand, Birjand, Iran. address: a.ashyani@birjand.ac.ir 2) Department of Mathematics, University of Birjand, Birjand, Iran. address: hmohammadin@birjand.ac.ir 3) Department of Mathematics, University of Birjand, Birjand, Iran. address: orabieimotlagh@birjand.ac.ir
Lotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationSTUDY OF THE DYNAMICAL MODEL OF HIV
STUDY OF THE DYNAMICAL MODEL OF HIV M.A. Lapshova, E.A. Shchepakina Samara National Research University, Samara, Russia Abstract. The paper is devoted to the study of the dynamical model of HIV. An application
More informationBifurcation and Stability Analysis of a Prey-predator System with a Reserved Area
ISSN 746-733, England, UK World Journal of Modelling and Simulation Vol. 8 ( No. 4, pp. 85-9 Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area Debasis Mukherjee Department
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
More informationDynamics of a Hepatitis B Viral Infection Model with Logistic Hepatocyte Growth and Cytotoxic T-Lymphocyte Response
Nonlinear Analysis and Differential Equations, Vol. 4, 16, no. 3, 19-1 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/nade.16.5164 Dynamics of a Hepatitis B Viral Infection Model with Logistic Hepatocyte
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationSensitivity and Stability Analysis of Hepatitis B Virus Model with Non-Cytolytic Cure Process and Logistic Hepatocyte Growth
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 12, Number 3 2016), pp. 2297 2312 Research India Publications http://www.ripublication.com/gjpam.htm Sensitivity and Stability Analysis
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
More informationCalculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y
Calculus for the Life Sciences II Assignment 6 solutions Find the tangent plane to the graph of the function at the point (0, π f(x, y = 3π 3 cos 2x + 2 sin 3y Solution: The tangent plane of f at a point
More informationGlobal Stability Analysis on a Predator-Prey Model with Omnivores
Applied Mathematical Sciences, Vol. 9, 215, no. 36, 1771-1782 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.512 Global Stability Analysis on a Predator-Prey Model with Omnivores Puji Andayani
More informationResearch Article On the Stability Property of the Infection-Free Equilibrium of a Viral Infection Model
Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume, Article ID 644, 9 pages doi:.55//644 Research Article On the Stability Property of the Infection-Free Equilibrium of a Viral
More informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More information6.2 Brief review of fundamental concepts about chaotic systems
6.2 Brief review of fundamental concepts about chaotic systems Lorenz (1963) introduced a 3-variable model that is a prototypical example of chaos theory. These equations were derived as a simplification
More informationMath 232, Final Test, 20 March 2007
Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.
More informationZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS
This is a preprint of: Zero-Hopf bifurcation for a class of Lorenz-type systems, Jaume Llibre, Ernesto Pérez-Chavela, Discrete Contin. Dyn. Syst. Ser. B, vol. 19(6), 1731 1736, 214. DOI: [doi:1.3934/dcdsb.214.19.1731]
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationMODELING AND ANALYSIS OF THE SPREAD OF CARRIER DEPENDENT INFECTIOUS DISEASES WITH ENVIRONMENTAL EFFECTS
Journal of Biological Systems, Vol. 11, No. 3 2003 325 335 c World Scientific Publishing Company MODELING AND ANALYSIS OF THE SPREAD OF CARRIER DEPENDENT INFECTIOUS DISEASES WITH ENVIRONMENTAL EFFECTS
More informationProject 1 Modeling of Epidemics
532 Chapter 7 Nonlinear Differential Equations and tability ection 7.5 Nonlinear systems, unlike linear systems, sometimes have periodic solutions, or limit cycles, that attract other nearby solutions.
More informationON THE EXISTENCE OF HOPF BIFURCATION IN AN OPEN ECONOMY MODEL
Proceedings of Equadiff-11 005, pp. 9 36 ISBN 978-80-7-64-5 ON THE EXISTENCE OF HOPF BIFURCATION IN AN OPEN ECONOMY MODEL KATARÍNA MAKOVÍNYIOVÁ AND RUDOLF ZIMKA Abstract. In the paper a four dimensional
More informationViral evolution model with several time scales
Viral evolution model with several time scales AA Archibasov Samara National Research University, 4 Moskovskoe Shosse, 4486, Samara, Russia Abstract In this paper a viral evolution model with specific
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationCONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS
International Journal of Bifurcation and Chaos, Vol. 12, No. 6 (22) 1417 1422 c World Scientific Publishing Company CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS JINHU LÜ Institute of Systems
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationA Mathematical Model for Transmission of Dengue
Applied Mathematical Sciences, Vol. 10, 2016, no. 7, 345-355 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.510662 A Mathematical Model for Transmission of Dengue Luis Eduardo López Departamento
More informationStochastic Viral Dynamics with Beddington-DeAngelis Functional Response
Stochastic Viral Dynamics with Beddington-DeAngelis Functional Response Junyi Tu, Yuncheng You University of South Florida, USA you@mail.usf.edu IMA Workshop in Memory of George R. Sell June 016 Outline
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationGlobal Stability of a Computer Virus Model with Cure and Vertical Transmission
International Journal of Research Studies in Computer Science and Engineering (IJRSCSE) Volume 3, Issue 1, January 016, PP 16-4 ISSN 349-4840 (Print) & ISSN 349-4859 (Online) www.arcjournals.org Global
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 0
MATH 56A: STOCHASTIC PROCESSES CHAPTER 0 0. Chapter 0 I reviewed basic properties of linear differential equations in one variable. I still need to do the theory for several variables. 0.1. linear differential
More informationAn Application of Perturbation Methods in Evolutionary Ecology
Dynamics at the Horsetooth Volume 2A, 2010. Focused Issue: Asymptotics and Perturbations An Application of Perturbation Methods in Evolutionary Ecology Department of Mathematics Colorado State University
More informationNONSTANDARD NUMERICAL METHODS FOR A CLASS OF PREDATOR-PREY MODELS WITH PREDATOR INTERFERENCE
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 67 75. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationHopf Bifurcation and Limit Cycle Analysis of the Rikitake System
ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Science Vol.4(0) No.,pp.-5 Hopf Bifurcation and Limit Cycle Analysis of the Rikitake System Xuedi Wang, Tianyu Yang, Wei Xu Nonlinear
More informationGlobal Analysis of a HCV Model with CTL, Antibody Responses and Therapy
Applied Mathematical Sciences Vol 9 205 no 8 3997-4008 HIKARI Ltd wwwm-hikaricom http://dxdoiorg/02988/ams20554334 Global Analysis of a HCV Model with CTL Antibody Responses and Therapy Adil Meskaf Department
More informationScenarios for the transition to chaos
Scenarios for the transition to chaos Alessandro Torcini alessandro.torcini@cnr.it Istituto dei Sistemi Complessi - CNR - Firenze Istituto Nazionale di Fisica Nucleare - Sezione di Firenze Centro interdipartimentale
More information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More informationChaos. Lendert Gelens. KU Leuven - Vrije Universiteit Brussel Nonlinear dynamics course - VUB
Chaos Lendert Gelens KU Leuven - Vrije Universiteit Brussel www.gelenslab.org Nonlinear dynamics course - VUB Examples of chaotic systems: the double pendulum? θ 1 θ θ 2 Examples of chaotic systems: the
More informationEE 380. Linear Control Systems. Lecture 10
EE 380 Linear Control Systems Lecture 10 Professor Jeffrey Schiano Department of Electrical Engineering Lecture 10. 1 Lecture 10 Topics Stability Definitions Methods for Determining Stability Lecture 10.
More informationDynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs
Dynamics of a opulation Model Controlling the Spread of lague in rairie Dogs Catalin Georgescu The University of South Dakota Department of Mathematical Sciences 414 East Clark Street, Vermillion, SD USA
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationDynamical Systems & Scientic Computing: Homework Assignments
Fakultäten für Informatik & Mathematik Technische Universität München Dr. Tobias Neckel Summer Term Dr. Florian Rupp Exercise Sheet 3 Dynamical Systems & Scientic Computing: Homework Assignments 3. [ ]
More informationHopf Bifurcation of a Nonlinear System Derived from Lorenz System Using Centre Manifold Approach ABSTRACT. 1. Introduction
Malaysian Journal of Mathematical Sciences 10(S) March : 1-13 (2016) Special Issue: The 10th IMT-GT International Conference on Mathematics, Statistics and its Applications 2014 (ICMSA 2014) MALAYSIAN
More informationSystem Control Engineering 0
System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More informationStability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations
Stability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations Sreelakshmi Manjunath Department of Electrical Engineering Indian Institute of Technology Madras (IITM), India JTG Summer
More informationLecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps
Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit
More informationA Study of Nonlinear Dynamics in Mathematical Biology
UNF Digital Commons UNF Theses and Dissertations Student Scholarship 2013 A Study of Nonlinear Dynamics in Mathematical Biology Joseph Ferrara University of North Florida Suggested Citation Ferrara, Joseph,
More informationIntroduction: What one must do to analyze any model Prove the positivity and boundedness of the solutions Determine the disease free equilibrium
Introduction: What one must do to analyze any model Prove the positivity and boundedness of the solutions Determine the disease free equilibrium point and the model reproduction number Prove the stability
More informationON FRACTIONAL ORDER CANCER MODEL
Journal of Fractional Calculus and Applications, Vol.. July, No., pp. 6. ISSN: 9-5858. http://www.fcaj.webs.com/ ON FRACTIONAL ORDER CANCER MODEL E. AHMED, A.H. HASHIS, F.A. RIHAN Abstract. In this work
More informationEen vlinder in de wiskunde: over chaos en structuur
Een vlinder in de wiskunde: over chaos en structuur Bernard J. Geurts Enschede, November 10, 2016 Tuin der Lusten (Garden of Earthly Delights) In all chaos there is a cosmos, in all disorder a secret
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationModels Involving Interactions between Predator and Prey Populations
Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationDynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting
Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting August 2012 Overview Last Model ẋ = αx(1 x) a(1 m)xy 1+c(1 m)x H(x) ẏ = dy + b(1 m)xy 1+c(1 m)x (1) where
More informationBifurcation, stability, and cluster formation of multi-strain infection models
J. Math. Biol. (23) 67:57 532 DOI.7/s285-2-6-3 Mathematical Biology Bifurcation, stability, and cluster formation of multi-strain infection models Bernard S. Chan Pei Yu Received: January 22 / Revised:
More informationState will have dimension 5. One possible choice is given by y and its derivatives up to y (4)
A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To
More informationLecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:
Lecture 6 Chaos Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Chaos, Attractors and strange attractors Transient chaos Lorenz Equations
More informationQualitative Analysis of Tumor-Immune ODE System
of Tumor-Immune ODE System L.G. de Pillis and A.E. Radunskaya August 15, 2002 This work was supported in part by a grant from the W.M. Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1. Simplified System
More informationCHEM 515: Chemical Kinetics and Dynamics
Alejandro J. Garza S01163018 Department of Chemistry, Rice University, Houston, TX email: ajg7@rice.edu, ext. 2657 Submitted December 12, 2011 Abstract Spontaneous antispiral wave formation was observed
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationStability Implications of Bendixson s Criterion
Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,
More informationGerardo Zavala. Math 388. Predator-Prey Models
Gerardo Zavala Math 388 Predator-Prey Models Spring 2013 1 History In the 1920s A. J. Lotka developed a mathematical model for the interaction between two species. The mathematician Vito Volterra worked
More informationM469, Fall 2010, Practice Problems for the Final
M469 Fall 00 Practice Problems for the Final The final exam for M469 will be Friday December 0 3:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
More informationTHE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL
THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationGlobal Analysis for an HIV Infection Model with CTL Immune Response and Infected Cells in Eclipse Phase
Article Global Analysis for an HIV Infection Model with CTL Immune Response and Infected Cells in Eclipse Phase Karam Allali 1 ID Jaouad Danane 1 ID and Yang Kuang 2 * ID 1 Laboratory of Mathematics and
More informationGLOBAL STABILITY OF SIR MODELS WITH NONLINEAR INCIDENCE AND DISCONTINUOUS TREATMENT
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 304, pp. 1 8. SSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GLOBAL STABLTY
More informationPERIODIC SOLUTIONS OF EL NIÑO MODEL THROUGH THE VALLIS DIFFERENTIAL SYSTEM
This is a preprint of: Periodic solutios of El niño model thorugh the Vallis differential system, Rodrigo D. Euzébio, Jaume Llibre, Discrete Contin. Dyn. Syst., vol. 34(9), 3455 3469, 214. DOI: [1.3934/dcds.214.34.3455]
More informationHepatitis C Mathematical Model
Hepatitis C Mathematical Model Syed Ali Raza May 18, 2012 1 Introduction Hepatitis C is an infectious disease that really harms the liver. It is caused by the hepatitis C virus. The infection leads to
More informationNew ideas in the non-equilibrium statistical physics and the micro approach to transportation flows
New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows Plenary talk on the conference Stochastic and Analytic Methods in Mathematical Physics, Yerevan, Armenia,
More informationBackward bifurcation underlies rich dynamics in simple disease models
Backward bifurcation underlies rich dynamics in simple disease models Wenjing Zhang Pei Yu, Lindi M. Wahl arxiv:1504.05260v1 [math.ds] 20 Apr 2015 October 29, 2018 Abstract In this paper, dynamical systems
More informationThe Dynamic Properties of a Deterministic SIR Epidemic Model in Discrete-Time
Applied Mathematics, 05, 6, 665-675 Published Online September 05 in SciRes http://wwwscirporg/journal/am http://dxdoiorg/046/am056048 The Dynamic Properties of a Deterministic SIR Epidemic Model in Discrete-Time
More informationStability analysis of a prey-predator model with a reserved area
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 4, 5(3):93-3 ISSN: 976-86 CODEN (USA): AASRFC Stability analysis of a prey-predator model with a reserved area Neelima
More informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
More informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationEpidemics in Two Competing Species
Epidemics in Two Competing Species Litao Han 1 School of Information, Renmin University of China, Beijing, 100872 P. R. China Andrea Pugliese 2 Department of Mathematics, University of Trento, Trento,
More information1 What is a differential equation
Math 10B - Calculus by Hughes-Hallett, et al. Chapter 11 - Differential Equations Prepared by Jason Gaddis 1 What is a differential equation Remark 1.1. We have seen basic differential equations already
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationIntermediate Differential Equations. John A. Burns
Intermediate Differential Equations Delay Differential Equations John A. Burns jaburns@vt.edu Interdisciplinary Center for Applied Mathematics Virginia Polytechnic Institute and State University Blacksburg,
More informationHopf bifurcations analysis of a three-dimensional nonlinear system
BULETINUL ACADEMIEI DE ŞTIINŢE A REPUBLICII MOLDOVA. MATEMATICA Number 358), 28, Pages 57 66 ISSN 124 7696 Hopf bifurcations analysis of a three-dimensional nonlinear system Mircea Craioveanu, Gheorghe
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationA model for transfer phenomena in structured populations
A model for transfer phenomena in structured populations Peter Hinow 1, Pierre Magal 2, Glenn F. Webb 3 1 Institute for Mathematics and its Applications, University of Minnesota, Minneapolis, MN 55455,
More informationarxiv: v1 [math.ds] 20 Sep 2016
THE PITCHFORK BIFURCATION arxiv:1609.05996v1 [math.ds] 20 Sep 2016 Contents INDIKA RAJAPAKSE AND STEVE SMALE Abstract. We give development of a new theory of the Pitchfork bifurcation, which removes the
More informationHarvesting Model for Fishery Resource with Reserve Area and Modified Effort Function
Malaya J. Mat. 4(2)(2016) 255 262 Harvesting Model for Fishery Resource with Reserve Area and Modified Effort Function Bhanu Gupta and Amit Sharma P.G. Department of Mathematics, JC DAV College, Dasuya
More informationThe Big, Big Picture (Bifurcations II)
The Big, Big Picture (Bifurcations II) Reading for this lecture: NDAC, Chapter 8 and Sec. 10.0-10.4. 1 Beyond fixed points: Bifurcation: Qualitative change in behavior as a control parameter is (slowly)
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationLecture 3 : Bifurcation Analysis
Lecture 3 : Bifurcation Analysis D. Sumpter & S.C. Nicolis October - December 2008 D. Sumpter & S.C. Nicolis General settings 4 basic bifurcations (as long as there is only one unstable mode!) steady state
More informationCan multiple species of Malaria co-persist in a region? Dynamics of multiple malaria species
Can multiple species of Malaria co-persist in a region? Dynamics of multiple malaria species Xingfu Zou Department of Applied Mathematics University of Western Ontario London, Ontario, Canada (Joint work
More informationSMA 208: Ordinary differential equations I
SMA 208: Ordinary differential equations I Modeling with First Order differential equations Lecturer: Dr. Philip Ngare (Contacts: pngare@uonbi.ac.ke, Tue 12-2 PM) School of Mathematics, University of Nairobi
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationDynamical systems tutorial. Gregor Schöner, INI, RUB
Dynamical systems tutorial Gregor Schöner, INI, RUB Dynamical systems: Tutorial the word dynamics time-varying measures range of a quantity forces causing/accounting for movement => dynamical systems dynamical
More informationPreprint núm May Stability of a stochastic model for HIV-1 dynamics within a host. L. Shaikhet, A. Korobeinikov
Preprint núm. 91 May 2014 Stability of a stochastic model for HIV-1 dynamics within a host L. Shaikhet, A. Korobeinikov STABILITY OF A STOCHASTIC MODEL FOR HIV-1 DYNAMICS WITHIN A HOST LEONID SHAIKHET
More informationDeterministic Changes - Chapter 5 of Heinz
Deterministic Changes - Chapter 5 of Heinz Mathematical Modeling, Spring 2019 Dr Doreen De Leon 1 Introduction - Section 51 of Heinz Our plan now is to use mathematics to explain changes in variables There
More information