Population Dynamics in a Microfounded Predator-Prey Model. Thomas Christiaans. University of Siegen. Discussion Paper No
|
|
- Darren Leonard
- 5 years ago
- Views:
Transcription
1 Volkswirtschaftliche Diskussionsbeiträge V W L Population Dynamics in a Microfounded Predator-Prey Model Thomas Christiaans University of Siegen Discussion Paper No ISSN x UNIVERSITÄT SIEGEN FACHBEREICH WIRTSCHAFTSWISSENSCHAFTEN
2 Population Dynamics in a Microfounded Predator-Prey Model Thomas Christiaans Department of Economics, University of Siegen, Siegen, Germany christiaans@vwl.wiwi.uni-siegen.de Abstract. This paper analyzes the dynamics of a two-dimensional microfounded predatorprey model. It is shown that the dynamics closely resemble those of a model commonly used in mathematical biology if parameters of the latter are suitably restricted. The positive equilibrium of the microfounded model is globally asymptotically stable for positive initial values largely irrespective of the parameter values chosen. If a version of Allee s Law is included, however, species extinction becomes possible. Keywords: population dynamics; ratio dependence; species extinction JEL-Classification: Q0 Introduction Mathematical biologists have recently extensively studied an alternative to the classical Lotka-Volterra model and its variations. The classical approach has been losing ground since some of its predictions are not in line with many field observations. The alternative theory relies on the so-called ratio-dependent predator-prey models. Kuang and Beretta 998) have analyzed the global dynamic behavior of the following two-dimensional ratio-dependent type predator-prey model: ṅ = n a ṅ = n n ) bn K n mn f n d n mn, ). ) Here, n and n denote the prey and the predator population, respectively, and a, b, d, f, m, and K are positive parameters whose biological interpretation can be found in Kuang and Beretta 998), e.g. These authors have shown that the system possesses a unique and globally asymptotically stable positive equilibrium that is, an equilibrium where both n and n are positive) for positive initial values if f > d and am b. While f > d is a necessary condition for the existence of such an equilibrium, the condition am b is sufficient for its global stability but not necessary for existence. If am < b, the global dynamic behavior of the model changes substantially. E.g., extinction of both species becomes possible. A further analysis of these equations can be found in Hsu et al. 00), where it is shown that even limit cycles and heteroclinic cycles may exist if am< b. It should be noted that the debate on the appropriate type of models is not yet finished and that the ratio-dependent models do also have their critics, cf. e.g. Abrams and Ginzburg 000) and Deng et al. 003).
3 A common feature of conventional predator-prey models is that they are macro approaches in the sense that populations and their development in time are the basic units of analysis. As Eichner and Pethig 004a) put it, these macro approaches neglect the processes at the micro level of preying and being preyed upon which ultimately generate the growth functions describing the population dynamics. Such a microfoundation of predator-prey models based on economic methodology has been initiated by Hannon 976) and further developed by Tschirhart 000), Pethig and Tschirhart 00), and Eichner and Pethig 004b), e.g. The present paper focuses on the following microfounded two-dimensional system introduced by Eichner and Pethig 004b): ) α ṅ = n A r0 α )n e n e n ) z n α )n e n e α n e ṅ = n A n ) z α n α z n n e n e α )n e ) α γ, 3) ) α γ. 4) For i =,, e i, γ i, ñ i and z i as well as r 0 are positive constants, respectively, while 0<α i <. Cf. Eichner and Pethig 004b) for the ecological or economic interpretation of these parameters. The functions A i are defined by { A i n i ) := min, n } i, i =,. 5) ñ i The inclusion of the functions A i corresponds to the idea that the organism s generation of net offspring is the more hampered, the further n i drops below some critical population level ñ i > 0. According to Allee s Law, species i may be called an endangered species if n i < ñ i cf. Berryman, 003). Regarding the microfoundation, equations 3) and 4) are derived in Eichner and Pethig 004b) from the maximization of a net-offspring function B i x i, z i, n i )= A i n i )x α i i z i z i ) α i γ i, i =,, of the representative organism of species i, where 0<α i < and γ i > 0, subject to the resource constraint e i p i z i p i x i, i =,. x i is organism i s intake of biomass of its prey species i and z i is organism i s loss of own biomass to its predator, species i. For i =, x i = x 0 is the demand of organism for a basic resource whose total supply per period is r 0 > 0. e i > 0 is some exogenous lump-sum income of species i and p i is the market price of biomass of species i. Since species is the top predator its biomass price is set p 0. The equilibrium conditions r 0 = n x 0, n z = n x are used to determine the equilibrium values of p 0, p, x 0, x, and z, which are substituted into the net-offspring functions. Equations 3) and 4) then follow by substitution into the dynamic equations ṅ i = B i n i, i =,.
4 Eichner and Pethig 004b) have analyzed the resulting dynamics of a threedimensional food chain by means of numerical simulations. Among their results is that the positive equilibrium, which they have shown to be unique, is approached for all positive initial values if one sets A i n i ), i =,,3. Thus, extinction in the three-dimensional system seems to be impossible unless Allee s Law is explicitly considered. The purpose of the present paper is to give a complete characterization of the dynamics in case of a two-dimensional system with just two species. As a slight generalization of the two-dimensional case, the parameters α i are replaced by β i 0 and the restrictions 0<α i < by α i > 0. Repeating the derivations of Eichner and Pethig 004b) for the two species case yields the differential equations ṅ = n A r0 n ) n ṅ = n A n ) z β ) α β n e n e z β n e n e α n e n α z n n e β n e n e α n e ) β γ, ) α γ. As z β is just a constant in the two-dimensional setting, it will be assumed that β = 0. This assumption simplifies the analysis without substantially altering the results β > 0 would merely necessitate a modification of some parameter assumptions). In this case, the differential equations can be given a slightly simpler appearance: ṅ = n A r0 n ) n ) α β n e n e z β n e n e α n e ṅ = n A n α z e n ) β n e n e α n e ) β γ, 6) ) α γ. 7) A further simplification follows by setting A i n i ), i =,, for the time being. Letting equations 6) and 7) take the form m := α β β e e, f := α e β, c := z e, ) α ) r0 cn e n e β ṅ = n γ, 8) n n mn ) f cn α ṅ = n γ. 9) n mn While equation 9) is now qualitatively equivalent to ) if α =, equation 8) cannot be transformed into equation ) using the present setting. Eichner and Pethig 004a) have shown, however, that equation ) can be derived from a slightly modified micro model of the ecosystem. The remainder of this paper is devoted to an analysis of equations 6) and 7). An interesting question is whether the dynamic properties of these equations are similar to those of equations ) and ). To answer this question, Section starts with an analysis of the simplified equations 8) and 9). The results will then be 3
5 used to consider the original equations 6) and 7) in Section 3. Section 4 provides a concluding comparison of the microfounded and the conventional model. It will be shown that the dynamics of 8) and 9) if f c > γ /α and 0 < α < closely resemble those of ) and ) if f > d and am b. While according to ) and ) extinction is possible if f > d and am < b, extinction in the microfounded model is impossible under reasonable assumptions unless Allee s Law is included, that is, unless equations 6) and 7) instead of 8) and 9) are used. Analysis of Equations 8) and 9). Equilibria and Local Stability As the right-hand sides of 8) and 9) are not defined at n, n ) = 0,0), it will be assumed that ṅ = ṅ = 0 if n = n = 0 by definition. Calculating the limit for lim n,n ) 0,0) shows that under this assumption both ṅ and ṅ are continuous on R = {n, n ) R : n 0, n 0} if 0<α <, although not differentiable at 0,0), cf. Appendix A. A similar argument applies with respect to equations 6) and 7). All assumptions supposed to hold in the sequel are summarized as follows. They will not be repeated each time a proposition is stated. Assumptions. All parameters appearing in equations 8) and 9) as well as in 6) and 7) are positive. In addition, 0<α <, 0<α, and f c > γ /α. The initial values of n and n are non-negative. If n = n = 0, then ṅ = ṅ = 0. These assumptions serve the following purposes: 0<α < ensures continuity of 8) and 9) on R, f c > γ/α is necessary for the existence of a positive equilibrium, and 0<α will be needed when analyzing equation 7). There are three equilibria. To begin with, set n = 0 and n > 0, implying ṅ = 0 and ṅ =γ n. Thus, there is a first trivial) equilibrium E 0 = 0,0) and there can be no other equilibrium where n = 0. It is obvious that 0,0) is stable along the n -axis. Thus, if there is no prey, the predator will become extinct. Setting n = 0 and n > 0 implies ṅ = 0 and ṅ = r α 0 z β nα γ n. It is straightforward that there is a second equilibrium E at n, n )= z β /α ) r 0,0 γ /α that is stable along the n -axis. Thus, if there is no predator, the prey population reaches a steady state that is directly proportional to the size of the basic resource, r 0. 4
6 A third equilibrium, E, entails positive populations of both species. If n > 0 and n > 0, setting ṅ = 0 yields n = f c γ/α n mγ /α = qn. 0) }{{} =:q Hence, a positive equilibrium cannot exist unless the condition f c > γ /α ) holds, which has already been assumed. Since γ /α should usually be a small number in relation to f c, this assumption is rather natural. Upon subsitution of 0) into ṅ = 0 one gets n = r 0 γ /α ) ce qe ) β /α, mq which together with 0) describes the positive equilibrium E. Observe that in any of the three equilibria thus considered the population of both species is proportional to the size of the basic resource, r 0. As for the stability of the positive steady state, the Jacobian of system 8), 9) evaluated at E can be shown to have the following pattern of signs cf. Appendix B): ṅ / n ṅ / n sgnj)=sgn ṅ / n ṅ / n ) = ) Therefore, TrJ)<0 and J >0, implying that the equilibrium E is locally asymptotically stable by the Routh-Hurwitz criterion. These results are summarized in Proposition There are three equilibria. E 0 = 0,0) is locally asymptotically stable along the n -axis, E = r 0 z β /α /γ /α, 0) is locally asymptotically stable along the n -axis, and ) r 0 ce qe ) β /α E =, q r ) ) 0 ce qe ) β /α mq mq γ /α is locally asymptotically stable. γ /α ). Global Stability The global dynamic behavior of solutions can be determined using the phase diagram, whose derivation is as follows. From 0), the isocline ṅ = 0 is a positively sloped straight line through the origin unionized with the n -axis) in n, n )-space if conditon ) holds. Setting ṅ = 0 and solving for n yields the following expression unionized with the n -axis) for the isocline ṅ = 0 : β e n z γ /β n /r 0 ) α /β n = e α β )γ /β n /r 0 ) α /β β z e. 3) 5
7 Considering the denominator, it is straightforward that denominator 0 n r 0 β α β z ) β /α γ /α =: n. 4) Thus, there is a vertical asymptote at n. n > 0, numerator 0 n r 0 z β /α Similarly, under the precondition that γ /α =: n. 5) Therefore, the isocline cuts the n -axis at the equilibrium E where n = n > n and at the equilibrium E 0, where n = 0). Comparing the sign patterns of the numerator and the denominator shows that the isocline lies below the n -axis for 0<n < n and n > n, while it lies above the axis if n < n < n. n ṅ = 0 ṅ = 0 E E0 n E n Figure. Phase Diagram of Equations 8), 9) As has been shown in Appendix B, the cross partial derivatives of 8) and 9) are ṅ n < 0 and ṅ n > 0, respectively. These partial derivatives give rise to the and signs indicating the directions of motion off the zero isoclines. Putting all this information together yields the phase diagram. For the sake of completeness, it is shown in Appendix C that the slope of the isocline ṅ = 0 is negative if n < n < n. Due to the existence of the vertical asymptote, there is one specialty to be taken care of. Above but near the isocline ṅ = 0 in the region where n is negative, ṅ < 0. As ṅ > 0 below ṅ = 0 in the positive region, the question arises where the sign of ṅ changes. Inspection of equation 8) shows that a sign change off the isocline ṅ = 0 derived from 3) or from n = 0) is possible only if 8) has a vertical asymptote at n mn = 0. Thus, if n > 0 the sign change must occur in the irrelevant region where n < 0. 6
8 The phase diagram reveals that the equilibrium E is globally asymptotically stable for strictly positive initial values. The proof relies on the fact that it is always possible to draw a rectangular closed region from which the trajectories cannot escape cf. the dashed rectangle in Figur ). This proves that the differential equations 8) and 9) have a solution defined for all t 0 Hirsch and Smale, 974, p. 7). According to the Generalized Poincaré-Bendixson Theorem cf. Perko, 996, p. 43), any limit point of trajectories must be an equilibrium if there exists neither a closed orbit nor a separatrix cycle. As the considered region contains just one equilibrium which is locally asymptotically stable, there is no separatrix cycle. Closed orbits are ruled out by the direction of movements in the four regions separated by the zero-isoclines, or, more rigorously, by Dulac s criterion cf. Appendix D). It follows that there is just one possible limit point of trajectories if t, the equilibrium E. This proves global stability. Notice that Figure gives a qualitatively complete picture of the dynamics of system 8), 9). If both initial values are positive, the equilibrium E will be reached and both species will survive. If the intial value of the predator species is zero, a positive inital population of the prey species will reach the equilibrium E. Finally, if the initial value of the prey species is zero, E 0 will be reached and the predator species becomes extinct. The analysis is summarized in Proposition Each trajectory of system 8), 9) converges to an equilibrium. If n 0)=0< n 0), lim t n t ), n t ))=E 0. If n 0)>0= n 0), lim t n t ), n t ))= E. If n 0)>0 and n 0)>0, lim t n t ), n t ))=E. 3 Analysis of Equations 6) and 7) 3. Partitioning the Phase Diagram Using similar parameters as in equations 8) and 9), equations 6) and 7) read ) α ) ṅ = n A r0 cn e n e β n ) γ, 6 ) n n mn ) ṅ = n A f cn α n ) γ, 7 ) n mn where A i n i ) as defined in 5) is repeated here for convenience: { A i n i ) := min, n } i, i =,. ñ i As a first step, observe that the phase diagram can now be partitioned into three regions depending on the values of ñ i. As long as n i ñ i, the phase diagram applies as before. For the vector field generated by 8) and 9) always points inwards any closed rectangular region lying in the positive quadrant, cf. Figure. If n i < ñ i for i = and/or i =, the dynamics will change. The result will depend on the relative position of the equilibrium E. Figure indicates the situation for the case where the equilibrium values of n and n both exceed ñ and ñ, respectively. 7
9 n ṅ = 0 ṅ = 0 E ñ E0 ñ n E n Figure. Partitioning the Phase Diagram It remains to determine the shape of the isoclines in the regions to the left of ñ and below ñ. As there are a lot of feasible parameter configurations giving rise to various details of the implied dynamics, this section solely relies on the analysis of representative phase diagrams without analytically analyzing local stability or instability of equilibria, which is straightforward in most cases, however. Similarly, the exclusion of closed orbits is not explicitely considered as the phase diagrams reveal that trajectories are always trapped in either the basin of attraction of the positive equilibrium or in a region where at least one species eventually becomes extinct. 3. Small Predator Population Below ñ, the shape of the isocline ṅ = 0 is left unchanged, while ṅ = 0 can be solved for n to yield γ ñ ) /α mn n =, 6) f cn /α γ ñ ) /α which is positive at n = ñ due to ). At n a := γ ñ /f c) α < ñ, there is a horizontal asymptote. As n declines further, n in 6) becomes negative. In addition, it is shown in Appendix E that α is a sufficient condition for the slope to be negative. Thus, the relevant part of the isocline has the shape shown in Figure 3. There emerges a new equilibrium E 3, which is unstable, however. As the direction arrows indicate, the equilibrium E may even be approached if n < ñ, that is, if the predator is an endangered species. Extinction of the endangered species is certain, however, if n n a. The system then approaches the equilibrium E on the n -axis, where only the prey survives. Regarding the policy implications of the complete model of Eichner and Pethig 004b), it should be noted that humans may influence the ecosystem by variation of the size of the basic resource, r 0. While the isocline ṅ = 0 is independent of r 0, ṅ = 0 will shift to the left as r 0 decreases. It is therefore interesting to analyze the behavior of the system as E vanishes sliding along ṅ = 0 when r 0 decreases. Such 8
10 n ṅ = 0 ṅ = 0 E ñ E 3 n a E0 E n Figure 3. Dynamics in Case of n < ñ a case, where the equilibrium E does not exist as it would involve an equilibrium value of n which is smaller than ñ, is shown in Figure 4. As the direction arrows indicate, the equilibrium E will be approached. Notice, however, that this result presupposes that n ñ, as the modifications of the dynamic system arising if n < ñ have not been taken into account yet. n ṅ = 0 ṅ = 0 ñ n a E0 E n Figure 4. Dynamics in Case of n < ñ if E Vanishes 3.3 Small Prey Population To the left of ñ, the shape of the isocline ṅ = 0 is left unchanged, while ṅ = 0 from 6 ) can be solved for n to yield n = β e n z n /β γ ñ ) /β n /r 0 ) α /β e α β )γ ñ ) /β n /r 0 ) α /β β z e n /β. 7) Recall that 0<α < and β > 0. As to the denominator, it is straightforward that denominator 0 n r α /α ) β 0 z α β ) β /α ) γ ñ ) /α ) =: na. 8) 9
11 Thus, there is a vertical asymptote at n a. Similarly, numerator 0 n r α /α ) 0 z β /α ) γ ñ ) =: /α ) n0. 9) Thus, the isocline cuts the n -axis at n 0 < na. Similar considerations as in Section. lead to the graph of the isocline shown in Figure 5. As to the and signs beneath the isocline, the remarks given in Footnote apply analogously. n n 0 n a n Figure 5. The Isocline ṅ = 0 if n < ñ With respect to the dynamics, it is important whether a part of the isocline 7) in the positive region lies to the left of ñ. Setting n 0 ñ shows that this condition is equivalent to ñ n : n 0 ñ ñ n. Thus, if ñ < n, there is a part of 7) for n > n 0 lying to the left of ñ. Similarly, setting n a ñ yields n a ñ ñ n. Therefore, the vertical asymptote of 7) lies to the left of ñ if and only if ñ lies to the left of the asymptote at n of equation 3). Proceeding under the reasonable assumption that ñ < n, it follows that two cases must be considered. n ṅ = 0 ṅ = 0 ṅ = 0 ñ n a E0 n 0 n a ñ n n Figure 6. Dynamics in Case of n < ñ if ñ n 0
12 First, assume that ñ n. Figure 6 shows a possible configuration in this case. Other configurations are possible, depending on the relative position of ñ. In case of Figure 6, the predator will become extinct if the initial values of n and n are such that the system starts below the isocline ṅ = 0 or sufficiently far to the left of the ṅ = 0 isocline. Even the prey may become extinct as the equilibrium at n = n 0 is unstable. Second, let n < ñ < n. This implies that the right-hand part and the left-hand part of ṅ = 0 are both valid only up to the point where n = ñ. Moreover, as 6 ) is continuous at n = ñ, the isocline itself is continuous here. One possible configuration is shown in Figure 7, where the positive equilibrium E vanishes. Depending on parameter values, this equilibrium could as well persist. In case of Figure 7, the predator cannot survive in the long run. n ṅ = 0 ṅ = 0 ñ n a E0 ñ n Figure 7. Dynamics in Case of n < ñ if n < ñ < n Putting all information together yields the overall phase diagram, whose appearance depends on the specific values of ñ and ñ. Figure 8 provides one example. n ṅ = 0 ṅ = 0 ṅ = 0 E ñ E 3 E0 E 4 ñ E n Figure 8. A Complete Phase Diagram
13 Although not all possible configurations have been analyzed in detail in this section, it is straightforward to consider the principle possibilities in terms of further phase diagrams. As with the examples given here, it is fairly obvious that the only possible limit points for t are equilibria. Exact proofs could be given along the lines in Section.. The following proposition summarizes the main results. Proposition 3 Let 0,0) < ñ, ñ ) < E. Then each trajectory of system 6) and 7) converges to an equilibrium. The positive equilibrium E is asymptotically stable and its basin of attraction covers at least {n, n ) R n ñ, n ñ }. There exist initial values n 0, ñ ) and n 0, ñ ), respectively, such that the predator or both species eventually become extinct. 4 Discussion If the functions A i n i ) are set A i n i ), the dynamics implied by equations 8) and 9) under conditions f c > γ /α and 0 < α < widely resemble the dynamics of the conventional model ) and ) if f > d and am b, although the differential equations themselves are rather distinct. In both cases, there is a unique and globally stable equilibrium where both species survive if both initial values are positive. If there is no prey, the predator becomes extinct, and if there is no predator, the prey reaches a positive equilibrium. In fact, the phase diagram of ) and ) can be shown to look exactly like Figure if f > d and am > b. Empirically, it would be impossible to distinguish whether a given set of observations was generated by model ) and ) or by 8) and 9). If f > d and am= b, the isocline ṅ = 0 becomes a straight line, leaving the qualitative implications unchanged, however. The condition f c > γ /α is necessary and given other assumptions about parameters following from the micro approach) sufficient for the existence of a positive equilibrium in the microfounded model. In contrast, f > d is just a necessary condition in the conventional model. If am < b, such an equilibrium exists only if f is suitably bounded from above cf. Kuang and Beretta, 998, p. 39). Under such circumstances, equations ) and ) can generate entirely different dynamics. E.g., it is possible that the positive equilibrium E is locally but not globally asymptotically stable, and one or both species could become extinct. Moreover, even limit cycles or heteroclinic cycles are possible cf. Hsu et al., 00). All these cases are excluded in the microfounded model. If f d or f c γ /α, respectively), the positive equilibrium disappears in both models. However, this case is rather irrelevant considering the microfounded model. Notice that f c = α z e /β e ) and that e, e and α, β should reasonably be of comparable magnitude, respectively, while z, the maximum amount of biomass that the prey could use for transactions, should reasonably exceed the natural death rate of the predator, γ <, raised to the power of /α. Thus, the analysis of this case is merely of theoretical interest. The microfounded model thus leaves no room for the empirically relevant case of species extinction, which is possible in case of the conventional model for am< b even if f > d. This result shows that the economics approach to ecology, where species engage in a kind of maximization process, can resemble the dynamics of
14 settled biological models but adds more stability by excluding extinction. This phenomenon can be reintroduced, however, by taking Allee s Law into account. As the analysis of equations 6) and 7) has shown, adding the A i n i )-functions leaves the dynamics unaltered for a region around the equilibrium with positive populations of both species but adds the possibility of extinction if the respective initial values are sufficiently small. Appendix A Continuity of 8) and 9) at the Origin Using the definition of m, it is obvious that n e n e e n mn ) on R, from which Applying some transformations yields ) α r0 cn e n e γ n n n n mn 0 n e n e n mn e. ) β γ ) α r0 n ce ) β γ n. n As n 0, both interval boundaries converge to zero if 0<α <. Thus, ) α ) lim n r0 cn e n e β γ = 0, n,n ) 0,0) n n mn proving continuity of 8) on R if ṅ = 0 for n = n = 0 by definition. Continuity of equation 9) is proven similarly. B Derivation of ) The partial derivative of 8) with respect to n evaluated at ṅ = 0 is ṅ n ṅ =0= ) ) n n me e ) α ) β n mn )n e n e ) α r0 cn e n e β A) n n mn Using the definition of m, it follows that me e = α β )e /β e = α e /β > 0. Substituting into A) shows that the first term in parentheses and therefore the entire expression is negative. The partial derivative of 8) with respect to n, ṅ n = β e ) me ) α ) r0 cn e n e β, n n mn )n e n e ) n n mn is negative since e me < 0. Finally, it is straightforward that the partial derivatives of 9) are and ṅ α mn ) f cn α = > 0 n n n mn ) n mn ṅ n ṅ =0 = α ) mn f cn α < 0, A) n mn n mn proving ). Notice that the signs of the cross partials are determined even off the isoclines. 3
15 C The Slope of 3) Differentiation of n = β e n z γ /β n /r 0 ) α /β e α β )γ /β n /r 0 ) α /β β z e with respect to n, letting D be an abbreviation for the denominator, yields: n n ṅ =0= β e z γ /β n /r 0 ) α /β α e γ /β n /r 0 ) α /β )D D e e α α β )γ /β n /r 0 ) α /β z γ /β n /r 0 ) α /β D. Recall relations 4) and 5). As the denominator is positive if n > n, this expression is negative if β z γ /β n /r 0 ) α /β α γ /β n /r 0 ) α /β )D < e α α β )γ /β n /r 0 ) α /β z γ /β n /r 0 ) α /β. The right-hand side of this inequality is positive as z > γ /β n /r 0 ) α /β if n < n. The lefthand side is negative as n > n implies that D> 0 and β z α β )γ /β n /r 0 ) α /β < 0, respectively. This proves that the isocline ṅ = 0 is negatively sloped if n < n < n. D Exclusion of Closed Orbits Applying Dulac s criterion cf. Perko, 996, p. 6) to equations 8) and 9), there is no closed orbit lying entirely in R = {n, n ) R : n > 0, n > 0} if there exists a function B C R ) such that the trace of the Jacobian of Bṅ,Bṅ ) is not identically zero and does not change sign in R. Now consider the function B= /n n ). The partial derivative A) has been calculated under the assumption that ṅ = 0, which has had just the effect that the term in square brackets in 8) has been omitted in A). Thus, it is straightforward that Bṅ ) = ṅ n n n n ṅ =0 for all n, n ) R. An analogous argument shows that, using A), Bṅ ) = ṅ n n n n ṅ =0 for all n, n ) R. As it follows from Appendix B that both expressions are negative, the trace of the Jacobian of Bṅ,Bṅ ) is negative for all n, n ) R, proving that there are no closed orbits lying entirely in R. E The Slope of 6) Consider the region where n a < n < ñ. As the denominator of the derivative of 6) with respect to n is positive, it suffices to consider the numerator, which is γ ñ ) /α m f cn /α γ ñ ) /α f cn /α mγ ñ ) /α. α 4
16 A sufficient condition for this expression to be negative is that α : α γ ñ ) /α m f cn /α γ ñ ) /α f cn /α mγ ñ ) /α < 0 α )f cn /α α γ ñ ) /α < 0. References Abrams, P. A. and Ginzburg, L. R. 000): The Natur of Predation: Prey Dependent, Ratio Dependent or Neither, Trends in Ecology and Evolution, 5, Berryman, A. A. 003): On Principles, Laws and Theory in Population Ecology, Oikos, 03, Deng, B., Jessi, S., Ledder, G., Rand, A., and Srodulski, S. 003): Biological Control Does Not Imply Paradox A Case Against Ratio-Dependent Models, Discussion Paper, University of Nebraska-Lincoln. Eichner, T. and Pethig, R. 004a): An Analytical Foundation of the Ratio-Dependent Predator-Prey Model, Discussion Paper No. 7-04, University of Siegen. 004b): Economic Land Use, Ecosystem Services and Microfounded Species Dynamics, Discussion Paper No. 6-04, University of Siegen. Hannon, B. 976): Marginal Product Pricing in the Ecosystem, Journal of Theoretical Biology, 56, Hirsch, M. W. and Smale, S. 974): Differential Equations, Dynamical Systems, and Linear Algebra, New York: Academic Press. Hsu, S.-B., Hwang, T.-W., and Kuang, Y. 00): Global Analysis of the Michaelis-Menten- Type Ratio-Dependent Predator-Prey System, Journal of Mathematical Biology, 4, Kuang, Y. and Beretta, E. 998): Global Qualitative Analysis of a Ratio-Dependent Predator- Prey System, Journal of Mathematical Biology, 36, Perko, L. 996): Differential Equations and Dynamical Systems, nd edition, New York: Springer. Pethig, R. and Tschirhart, J. 00): Microfoundations of Population Dynamics, Journal of Bioeconomics, 3, Tschirhart, J. 000): General Equilibrium of an Ecosystem, Journal of Theoretical Biology, 03,
17 ! " # $ $ % & # " % ' % # ) *, - ). / # &. / %, " % 0 ) ' " % * ) 3 ) ' % % % 6 ) & 7 ' 0 # ) " % ) / 8 9: : ; " % < ' % ; : = 6-6 : # % :, " % 0 ) # ' ; / $ ; > 0 $ 3 ) # < 5 A, B % % % C % & ) $ # % % % % # ". / " % ) ) > ) / 8 9: : # ; ) 8. ; ) ' % ' / % - ) % ; > % & ) # ' % % #. /, " % 0 ) ' % % ), # % > " ) % D " ). / % % > " % # / $ & % # % E ) & ; ) ; F ; E / ' G % 6 ) H ' % I J H ' % ; K L M N O / / # 0 ) #. & #. " % 8 # 8 ) # % / & " & / 8 # 8 ) % # ) # 6 # # 0 % % " % ) / 8 9: : ; " % < ' % ; : = 6-6 : # % :, " % 0 ) # ' ; / $ ; H # ) % ' - / 8 # 8 < 5 A / % & ) $ # %. # % # 0 #.. # / 8 9: : # ; ) 8. ; ) ' ;. " % E # 8 ). # % 0 % P ) ) & ) $ / # " / ) ). P %. 8 % #. # & ) $ E ) & ; ) ; F ; E / ' G % 6 ) P & H ' % < I J H ' % Q ) $ # % P ; R J < 5 R S T H # < E % P % # $. % U % < > " % $ " V & V " < H. ) Q ) - / - / * # ) # 0 F " ) % H. # R 5 < 5 R S T H ) # ' % ) C % 6 % 8, % % ) F ' % 9 ) W # X #. / " $ $, % # % ) H, ) ) ", " ) < 5 R Z Z ; T V % # ) P E. P % \ " ) 8 9 O - # ) # \ " ) 8 # % % ) # 3 # %, # % ^ % \ " ) 8 # % " ) ) %. P Y 4 < 5 R Z Z ; T \. / # % ' F # O # ) '! % # % C % ) % # % # O ) # 9 O / C $ 8 ) # %. & / E. P V #, % ' W ) # $ - ), Y _ < 5 R ` a b " % c a L N ^ ' 8. / V # ), 6 ) / # % ) 3 # %, % Y R < 5 R Z M # % S T H 8. &. C % 8 " % $ 8 6 \ d " 0 ) # - /. ) # % ' F " ) % H. # Y Y < 5 R Z M # % S T > 6 ) F P 0. D P %, \ & &. Q % ) # & ) $ H. #. % $ Y I < 5 R S e M O * < V % 8 P > 6 ) % ' # % E ) ' ) # $ f " # P Y? < 5 R Z Z ; T O % ' % ' < F " % U " ) # P 9 \ $ 8 ). # \ 6 %. & ) Q A < " % ) - / H 8. # \ $ 8 / # % Q ) $ # % P Y A < 5 Y S e M \ & &. % V # % # ' $ % & X # ) f " # P Y J < 5 Y ` g U # " ) - %. / # &. / Q ) " % # ' % % # 7 ). / ) H 0 ) % ' " % ' 8 ) D % X # ) ) " ). % Y 5 < 5 Y S e M U %. 8 ) # 6 U # % # \ % 6 ) % $ % # E. # % C % ) % # % # # 8 # V 0 P < 5 Y ` g O / % % & ) \. '. # H " # % # $ % % / % & # " 0 < $ H & 8 " ) &. # % V & # % > ' ' ) ' # X # ) F. " ) I 4 < 5 I h i L * ) - % " % ' \ " ) < O / ) $ D " $ 3 - > % ' < " 8 < O / ) $ $ 8 D ) % % X ) < 8 ) ". / I _ < 5 I h i L j 0 ) B & & %. / Q 7 ) " % 7 0 ) Q 7 ) " $ ) % Q 0 ) # ". / $ # % %. / ) 6 # ) %, # % % I R < 5 I c ` L k Z C % ) % # % # \ % 6 ) % $ % # U ' # % - / $ 8 % # % ) F ) 0 " % I Y < 5 I T i # % l T % C % ' ) # C % > 8 8 X & # ) \. % $. 9 > H " $ $ # ) P ^ & 3 #. O / ) $ I I < 5 I S e M C % & ) $ # % # X ). / # & ' " I? < 5 I c ` L k Z > % \ 8 ) $ % # H " P % 3 # ) ' # % % ' 3 / # 6 ) % \. % $. # % \ /. # \ % 6 ) % $ % I A < 5? S e M \. '. # O # F & ) $ # % \ & &. %. P & O # # % 9 > E " 0. Q E ) 8. 6 I J < 5? T! " ) F # ) " % ' % ) 8 8 % 6 % % ) B, '. / % H " ) ) & ) $ I 5 < 5? Z T ) 8 ) # W % # %. # 8 # V # ), \ d " 0 ) " $ # % C % ) % # % # O # $ 8 % - / # 8 # C %. $ O < 5 A T > $ $ % % / C $ 8 #. & / C % # O # V % / 6 % & # % \ % 6 ) % $ % # O # F & ) $? 4 < 5 A S e M \ $ % O # F 6 % " % # Q ) - % ' \. % $ P? _ < 5 A Z T E # P < # < P " < ' E % % H P $ # C %. $ 8 H. # % ) #.? R < 5 A Z T H ) # '. 3 " % O # # % - / % W % # %. # % E ) &. % # ) \ % ' % "? Y < 5 A m L a E ) ". 6 % " $ 8 % # % Q ) - / % 6 8 % ' " % )? I < 5 J c S > ) % # 6 >, # % $. / # % $ % & 7 ) % F " % & " %, 8 ) ' ) # $ $ 0 % ' % ) E ) ' ) # $ $ d " #?? < 5 J l e K > $ 8 # ) % & H " % V / & ) F # 3 " % P. V? A < 5 J M a! " $ \ % & " n 6 % 6 % ) $ % '. / & % # " & # V # ), '. / ' -. / 0 # P $ $ ). / ) C % & ) $ # %? J < 5 J m L o > U % E " 0. Q 9 U % < \. " # 0 P C $ 8 p % % " $ # 0 P? 5 < 5 J c h H P D W #. # % C % ) % # % # 3 " % P. < O / Q ) $ # % # < 5 J m L K O / # #. # % " ) ) 9 # 0 " ) " 8 8 P # % % 6 $ % % / " $ # %. # 8 # A 4 < 5 J m L a > ' ' ) ' # \. % $. Q ) - / - / H " 0 %. % " $ 8 % A _ < 5 J Z T C $ 8 $ % % ' \ d " # 6 % ' % % % # W ) # % A R < 5 5 m L K # % S e M E ) ". ' % # % V # ), & ) F. P. % ' X # O ) # $ % # % 8 #
18 q r s t t u v w v x y v v w z { v x } ~ ƒ ˆ Š ˆ ƒ Œ Œ Ž ˆ ˆ Ž Œ ˆ ˆ Ž Œ ˆ ƒ ˆ Ž Œ ƒ ˆ ƒ q s t t z { š z x w { { œ } ˆ ž ƒ Œ Ÿ ž ˆ ž ƒ ž Ž ˆ ž Œ ž ƒ ž ˆ Œ ˆ q s t t z { z œ v x Œ ƒ ª «v x u v w z } Ž ž ˆ ƒ Œ ƒ ˆ ƒ Œ ƒ Œ ˆ ˆ ƒ ž ˆ ƒ Œ ƒ ± Œ ˆ ² Ž ˆ Œ ˆ ƒ q q s t t ª «v x u v w z } ³ ƒ ˆ Ž ˆ ž ƒ Ž ž ƒ ˆ ƒ Œ ž ƒ ž q s t t µ { x œ v } Ž ž Ž Œ } ž ± Œ ˆ Œ ƒ Œ ˆ ˆ ƒ } Œ ƒ Œ Ž ˆ Ž Ž ˆ q t s t t ¹ { v œ º» ¼ œ } Œ ˆ ƒ ˆ ƒ ž ƒ Œ ˆ Ž ˆ ˆ ƒ Œ ž ƒ ž Œ Ž ž ž ƒ ² ˆ ƒ ž ž ½ s t t ¹ { v œ º» ¼ œ } ˆ Œ ƒ ƒ ž ƒ ² ž ¾ ˆ ˆ ˆ ƒ ƒ ˆ Ž ž ƒ ž À s t t µ { x œ v } Ž s ˆ Œ ƒ ˆ ž ƒ s ˆ ˆ Œ ž ƒ ž ž ± Œ ˆ ƒ ˆ Ž Á ˆ Ž ˆ Ÿ ž ˆ ž ƒ ƒ Œ Ž Œ ˆ  ž ž Á ƒ Ž à s t t µ { x œ v } Ä ž ˆ ž ƒ Å ž ƒ ² Œ Æ Œ ž ƒ ƒ Ž Œ ˆ  ž ž ³ ž ž Å Ç s t t z { z œ v x Œ ƒ ª «v x u v w z } ˆ ƒ } Ž ž ˆ Ž ˆ ž ƒ Œ ƒ Ÿ ˆ ƒ Ž Œ ¾ ˆ ± Œ ˆ Œ ƒ Œ ˆ ˆ ƒ r s ½ ½ z { z œ v x ƒ ª «v x u v w z } Œ  ˆ È Ž ˆ ƒ ž ƒ ¾ ˆ ˆ Ž Œ ˆ ˆ ˆ Œ ƒ Ä Â É Ê Œ ˆ Ë ž ƒ ž s ˆ ƒ Ž È Œ ƒ s ½ ½ z { z œ v x ƒ ª «v x u v w z }  ˆ È Ž ˆ ƒ Ž Œ ˆ ˆ ƒ ƒ ˆ ƒ ˆ Œ Ž ˆ ˆ Ž ˆ ƒ ž ˆ ˆ Ž Ä Œ Ž Œ s ½ ½ ª «v x u v w z { œ «y { Ì Í v œ š z v œ } Ÿ Ž Œ  ž ž Ÿ ž ƒ ž ƒ Œ ƒ Ÿ Ž Œ Ÿ Œ Œ q s ½ ½ µ z { v Î { } ³ Œ ž ƒ ˆ Œ Ž ž ƒ Œ ƒ ³ ¾ ƒ ž Œ ƒ Œ ž ˆ Ž ˆ s ½ ½ z { z œ v x Œ ƒ µ { x œ v } ˆ ƒ Œ ƒ Œ ƒ Œ ˆ Œ ƒ Œ ˆ ˆ ƒ ž Ž ž Ž Œ Œ ƒ ˆ Œ t s ½ ½ Í { w v x º z x { œ «z { š z x w { { œ } ž ƒ ž ˆ ž ƒ Ä Ž ž ˆ Ž ƒ ˆ Ï ˆ ž ² Œ ˆ Œ ž Ž Œ  ˆ ƒ ˆ t ½ s ½ ½ Í { w v x º z x } Ä Œ Ž ž ˆ ž ƒ ž  Ž ž ž ˆ ž Ÿ ž ˆ ƒ ˆ ž ƒ t À s ½ ½ ¹ { v œ º» ¼ œ } Ÿ ž ˆ Ž Œ ž ƒ ž ˆ Œ ž Œ ƒ ˆ ž Ž ˆ Œ ž ƒ  Ž ž ƒ ˆ ž ƒ t à s ½ ½ y { Ì Í v œ š z v œ { œ «Ð œ «x v { Í { v œ v x } Ä Ž Œ ƒ ž ƒ Œ ž ƒ t Ç s ½ À Ñ ª x v œ z v œ }  ˆ ž ˆ Ž Œ ˆ ˆ ƒ Ò ˆ ˆ ˆ ƒ ˆ Œ ƒ ƒ ˆ Ž ž Ó ˆ ~ ˆ ƒ Ž Œ Œ ƒ Ô ˆ Ž ˆ t r s ½ À z { š z x w { { œ } ž ƒ ž  Ž ž } ˆ Œ ˆ Œ Œ ˆ ƒ } Œ ƒ Œ  ž ˆ ƒ ˆ ž ² t s ½ À z { š z x w { { œ } ž ƒ ž  Ž ž } Œ  ž ˆ ƒ ˆ ž ² } Œ ƒ ˆ Œ Ž ƒ ƒ ž ƒ t s ½ À µ z { v Î { } ˆ Ž ˆ ƒ Œ ƒ ž ƒ ˆ ˆ Œ ƒ ž Ž ž ƒ ƒ Œ ž Œ  ž ˆ Ž ˆ t q s ½ À ª «v x u v w z } Ä Ž Ž ˆ } ˆ ˆ Œ ƒ ˆ ˆ ž ˆ t s ½ À ¹ { v œ º» ¼ œ } Œ Œ ˆ Ž ž Ž Œ ˆ Ž Ž Ë Œ Æ ˆ Ž ƒ ƒ ˆ Ž Ä ˆ ƒ Œ ƒ ˆ ˆ ž Ž ˆ t t s ½ À z { z œ v x { œ «Ð œ «x v { Í { v œ v x } ž Ž ˆ ž ƒ Œ Ž Œ ˆ Ž Ÿ Œ Ž Œ ˆ Ž ¾ Œ ž ƒ ž Ä ˆ Ž ž ƒ Œ ƒ Ž ˆ ƒ ˆ À ½ ½ s ½ À ª «v x u v w z } Œ ˆ ƒ ˆ ˆ ƒ } ± ˆ Ž ƒ ƒ Ó Ž ˆ Ò ƒ ˆ Ž Œ ˆ ž Ž ˆ ˆ Ä ƒ Œ ˆ À ½ À s ½ Ã Õ { x Ì Ñ v Ö Õ z } ˆ ž ƒ Œ Œ ƒ ˆ  Ž ž Ï ³ ƒ ˆ Ä ƒ Œ ž  Ž ž ² Ž Œ ˆ ž Ž ˆ À ½ à s ½ à ª «v x u v w z } Š ž ž Á ƒ ˆ Ž ƒ Œ ¾ ˆ ž ž ƒ Æ ˆ Ž ƒ Œ ˆ ² Ž ž Å Æ ˆ Ž ˆ ƒ ƒ Å ² Œ Æ ˆ À ½ Ç s ½ à µ z { v Î { } ˆ Ž ˆ ƒ Œ ƒ ž ƒ ˆ ˆ Œ ƒ ž Ž ž ƒ ƒ Œ ž Œ  ž ˆ Ž ˆ Œ ƒ Ÿ Œ Œ À ½ r s ½ à ¹ { v œ º» ¼ œ } ƒ Œ ˆ ƒ Œ ž Ž ž ž ƒ ² ˆ ž Ž ƒ Á ƒ ˆ Ž ƒ Œ ž ƒ Œ ² Ž Œ ˆ Ä ž ˆ Ž ƒ Ä Ž ž Œ Œ ˆ ž ƒ ² ˆ ž Ž ž Œ À ½ s ½ Ç ª «v x u v w z } ² ˆ Å Œ ˆ Ž Œ Œ Œ ƒ ˆ Œ Ž ž Œ Å ž ž ž ƒ Ï ž Ž ƒ } Œ ž ƒ Œ ƒ Œ ˆ Œ ƒ ˆ À ½ s ½ Ç ª «v x u v w z { œ «Ð œ «x v { Í { v œ v x } Ž ž ² Œ Æ Ÿ ž ˆ ž ƒ Œ ƒ ž Ž Œ Ä ž Ž ž ƒ ˆ ƒ À ½ q s ½ Ç Í { w v x º z x } ± Œ ƒ Ž Œ Ž Ž ˆ Ø À ½ s ½ Ç z { z œ v x } Á ˆ Ž ˆ Ÿ ž ˆ ž ƒ ƒ ˆ ˆ ƒ Á ƒ Ž À ½ t s ½ Ç z { z œ v x { œ «ª «v x u v w z } ² ˆ Œ ž Œ Ž Œ ƒ Œ ƒ Œ ƒ ˆ ƒ ž ž Œ ƒ ž ƒ ˆ ˆ ž Œ ž ƒ ƒ Œ À À ½ s ½ Ç z { z œ v x { œ «ª «v x u v w z } Ä Ž ž ž ƒ Œ ž ƒ ž Ž ˆ Œ ž Ž s Ž ˆ ƒ Œ À À À s ½ Ç µ z { v Î { } Š Œ ˆ Ž ˆ ƒ ˆ ƒ Ÿ ž ƒ ž ƒ ƒ Œ Ž ˆ ž ˆ ž ˆ ƒ ˆ Ÿ ˆ À À à s ½ Ç z { š z x w { { œ } Ä ƒ ƒ Œ ˆ ž Œ Œ ² ˆ ž Ž ž Œ ž Ž ˆ Œ ƒ À À Ç s ½ Ç z { š z x w { { œ } ž ƒ s Œ ˆ  Ž ž } ƒ ž ˆ ƒ ž Ÿ ž Œ Ž Œ ˆ Ä Œ ƒ Œ ˆ } Œ ƒ Á ƒ Ž Œ ¾ Œ ž ƒ À À r s ½ r µ z { v Î { } ± Œ ˆ ƒ Œ ž Œ  ž ˆ ž ˆ ƒ ˆ Ÿ ˆ À À s ½ r z { z œ v x { œ «ª «v x u v w z } ˆ ƒ ƒ ž ƒ Œ ƒ Ž ž ž ˆ ƒ Ž ƒ Œ Ž ˆ Ž ž ˆ ž ƒ À À s ½ r z { z œ v x { œ «ª «v x u v w z } ž ƒ ž Œ ƒ ˆ } ˆ ž ˆ ˆ Ž ˆ Œ ƒ Ž ž ž ƒ ˆ ˆ ˆ ƒ Œ À À q s ½ r z { z œ v x { œ «ª «v x u v w z } Ä ƒ Œ ƒ Œ Œ ž ƒ Œ ž ƒ ž ˆ Ž Œ ž s ˆ ˆ ƒ ˆ ƒ Ž ˆ Œ ž Ž s Ž ˆ ž ˆ
Framework for functional tree simulation applied to 'golden delicious' apple trees
Purdue University Purdue e-pubs Open Access Theses Theses and Dissertations Spring 2015 Framework for functional tree simulation applied to 'golden delicious' apple trees Marek Fiser Purdue University
More information3.5 Competition Models: Principle of Competitive Exclusion
94 3. Models for Interacting Populations different dimensional parameter changes. For example, doubling the carrying capacity K is exactly equivalent to halving the predator response parameter D. The dimensionless
More informationThe construction and use of a phase diagram to investigate the properties of a dynamic system
File: Phase.doc/pdf The construction and use of a phase diagram to investigate the properties of a dynamic system 1. Introduction Many economic or biological-economic models can be represented as dynamic
More informationThe Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor
, pp. 35-46 http://dx.doi.org/10.14257/ijbsbt.2017.9.3.04 The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor Alemu Geleta Wedajo
More informationChecking Consistency. Chapter Introduction Support of a Consistent Family
Chapter 11 Checking Consistency 11.1 Introduction The conditions which define a consistent family of histories were stated in Ch. 10. The sample space must consist of a collection of mutually orthogonal
More informationGlobal Qualitative Analysis for a Ratio-Dependent Predator Prey Model with Delay 1
Journal of Mathematical Analysis and Applications 266, 401 419 (2002 doi:10.1006/jmaa.2001.7751, available online at http://www.idealibrary.com on Global Qualitative Analysis for a Ratio-Dependent Predator
More informationx 2 F 1 = 0 K 2 v 2 E 1 E 2 F 2 = 0 v 1 K 1 x 1
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 20, Number 4, Fall 1990 ON THE STABILITY OF ONE-PREDATOR TWO-PREY SYSTEMS M. FARKAS 1. Introduction. The MacArthur-Rosenzweig graphical criterion" of stability
More informationHARVESTING IN A TWO-PREY ONE-PREDATOR FISHERY: A BIOECONOMIC MODEL
ANZIAM J. 452004), 443 456 HARVESTING IN A TWO-PREY ONE-PREDATOR FISHERY: A BIOECONOMIC MODEL T. K. KAR 1 and K. S. CHAUDHURI 2 Received 22 June, 2001; revised 20 September, 2002) Abstract A multispecies
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 informationDYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD
Journal of Sustainability Science Management Volume 10 Number 2, December 2015: 16-23 ISSN: 1823-8556 Penerbit UMT DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD SARKER MD SOHEL
More informationConnection equations with stream variables are generated in a model when using the # $ % () operator or the & ' %
7 9 9 7 The two basic variable types in a connector potential (or across) variable and flow (or through) variable are not sufficient to describe in a numerically sound way the bi-directional flow of matter
More informationPredator-Prey Model with Ratio-dependent Food
University of Minnesota Duluth Department of Mathematics and Statistics Predator-Prey Model with Ratio-dependent Food Processing Response Advisor: Harlan Stech Jana Hurkova June 2013 Table of Contents
More information4.3 Laplace Transform in Linear System Analysis
4.3 Laplace Transform in Linear System Analysis The main goal in analysis of any dynamic system is to find its response to a given input. The system response in general has two components: zero-state response
More informationLA PRISE DE CALAIS. çoys, çoys, har - dis. çoys, dis. tons, mantz, tons, Gas. c est. à ce. C est à ce. coup, c est à ce
> ƒ? @ Z [ \ _ ' µ `. l 1 2 3 z Æ Ñ 6 = Ð l sl (~131 1606) rn % & +, l r s s, r 7 nr ss r r s s s, r s, r! " # $ s s ( ) r * s, / 0 s, r 4 r r 9;: < 10 r mnz, rz, r ns, 1 s ; j;k ns, q r s { } ~ l r mnz,
More informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
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 informationObservations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationApplications in Biology
11 Applications in Biology In this chapter we make use of the techniques developed in the previous few chapters to examine some nonlinear systems that have been used as mathematical models for a variety
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 informationDynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
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 informationAn Introduction to Optimal Control Applied to Disease Models
An Introduction to Optimal Control Applied to Disease Models Suzanne Lenhart University of Tennessee, Knoxville Departments of Mathematics Lecture1 p.1/37 Example Number of cancer cells at time (exponential
More information1. Allstate Group Critical Illness Claim Form: When filing Critical Illness claim, please be sure to include the following:
Dear Policyholder, Per your recent request, enclosed you will find the following forms: 1. Allstate Group Critical Illness Claim Form: When filing Critical Illness claim, please be sure to include the
More informationdv dt Predator-Prey Models
Predator-Prey Models This is a diverse area that includes general models of consumption: Granivores eating seeds Parasitoids Parasite-host interactions Lotka-Voterra model prey and predator: V = victim
More informationMax. Input Power (W) Input Current (Arms) Dimming. Enclosure
Product Overview XI025100V036NM1M Input Voltage (Vac) Output Power (W) Output Voltage Range (V) Output urrent (A) Efficiency@ Max Load and 70 ase Max ase Temp. ( ) Input urrent (Arms) Max. Input Power
More informationContinuous time population models
Continuous time population models Jaap van der Meer jaap.van.der.meer@nioz.nl Abstract Many simple theoretical population models in continuous time relate the rate of change of the size of two populations
More informationNa#onal Center for Theore#cal Sciences Mathema#cs Division, Taiwan
Na#onal Center for Theore#cal Sciences Mathema#cs Division, Taiwan - TAMKANG JOURNAL OF MATHEMATICS Volume 47, Number 1, 127-141, March 2016 doi:10.5556/j.tkjm.47.2016.1984 This paper is available online
More informationHarvesting in an integrated general equilibrium model
Harvesting in an integrated general equilibrium model Thomas Eichner and Rüdiger Pethig University of Siegen Department of Economics Hoelderlinstr. 3 D-57068 Siegen Germany phone: ++49 (0) 271 740 3164
More informationB œ c " " ã B œ c 8 8. such that substituting these values for the B 3 's will make all the equations true
System of Linear Equations variables Ð unknowns Ñ B" ß B# ß ÞÞÞ ß B8 Æ Æ Æ + B + B ÞÞÞ + B œ, "" " "# # "8 8 " + B + B ÞÞÞ + B œ, #" " ## # #8 8 # ã + B + B ÞÞÞ + B œ, 3" " 3# # 38 8 3 ã + 7" B" + 7# B#
More informationOptimal Control of PDEs
Optimal Control of PDEs Suzanne Lenhart University of Tennessee, Knoville Department of Mathematics Lecture1 p.1/36 Outline 1. Idea of diffusion PDE 2. Motivating Eample 3. Big picture of optimal control
More information! " # $! % & '! , ) ( + - (. ) ( ) * + / 0 1 2 3 0 / 4 5 / 6 0 ; 8 7 < = 7 > 8 7 8 9 : Œ Š ž P P h ˆ Š ˆ Œ ˆ Š ˆ Ž Ž Ý Ü Ý Ü Ý Ž Ý ê ç è ± ¹ ¼ ¹ ä ± ¹ w ç ¹ è ¼ è Œ ¹ ± ¹ è ¹ è ä ç w ¹ ã ¼ ¹ ä ¹ ¼ ¹ ±
More informationGrowth models for cells in the chemostat
Growth models for cells in the chemostat V. Lemesle, J-L. Gouzé COMORE Project, INRIA Sophia Antipolis BP93, 06902 Sophia Antipolis, FRANCE Valerie.Lemesle, Jean-Luc.Gouze@sophia.inria.fr Abstract A chemostat
More informationOC330C. Wiring Diagram. Recommended PKH- P35 / P50 GALH PKA- RP35 / RP50. Remarks (Drawing No.) No. Parts No. Parts Name Specifications
G G " # $ % & " ' ( ) $ * " # $ % & " ( + ) $ * " # C % " ' ( ) $ * C " # C % " ( + ) $ * C D ; E @ F @ 9 = H I J ; @ = : @ A > B ; : K 9 L 9 M N O D K P D N O Q P D R S > T ; U V > = : W X Y J > E ; Z
More informationOBSERVABILITY AND OBSERVERS IN A FOOD WEB
OBSERVABILITY AND OBSERVERS IN A FOOD WEB I. LÓPEZ, M. GÁMEZ, AND S. MOLNÁR Abstract. The problem of the possibility to recover the time-dependent state of a whole population system out of the observation
More information3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 REDCLIFF MUNICIPAL PLANNING COMMISSION FOR COMMENT/DISCUSSION DATE: TOPIC: April 27 th, 2018 Bylaw 1860/2018, proposed amendments to the Land Use Bylaw regarding cannabis
More informationWorkshop on Theoretical Ecology and Global Change March 2009
2022-3 Workshop on Theoretical Ecology and Global Change 2-18 March 2009 Stability Analysis of Food Webs: An Introduction to Local Stability of Dynamical Systems S. Allesina National Center for Ecological
More informationKey words and phrases. Bifurcation, Difference Equations, Fixed Points, Predator - Prey System, Stability.
ISO 9001:008 Certified Volume, Issue, March 013 Dynamical Behavior in a Discrete Prey- Predator Interactions M.ReniSagaya Raj 1, A.George Maria Selvam, R.Janagaraj 3.and D.Pushparajan 4 1,,3 Sacred Heart
More informationNUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD CHAIN WITH LOTKA-VOLTERRA LINEAR FUNCTIONAL RESPONSE
Journal of Sustainability Science and Management Volume 6 Number 1, June 2011: 44-50 ISSN: 1823-8556 Universiti Malaysia Terengganu Publisher NUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationA Discrete Model of Three Species Prey- Predator System
ISSN(Online): 39-8753 ISSN (Print): 347-670 (An ISO 397: 007 Certified Organization) Vol. 4, Issue, January 05 A Discrete Model of Three Species Prey- Predator System A.George Maria Selvam, R.Janagaraj
More information14 Periodic phenomena in nature and limit cycles
14 Periodic phenomena in nature and limit cycles 14.1 Periodic phenomena in nature As it was discussed while I talked about the Lotka Volterra model, a great deal of natural phenomena show periodic behavior.
More informationLEAD IN DRINKING WATER REPORT
181 US Hwy 46 Mine Hill, NJ 07803 (908) 654-8068 (800) 783-0567 Fax 908-654-8069 LEAD IN DRINKING WATER REPORT Performed At: Montrose School 356 Clark St., South Orange, NJ 07079 Performed For: South Orange
More informationTopic 6: Projected Dynamical Systems
Topic 6: Projected Dynamical Systems John F. Smith Memorial Professor and Director Virtual Center for Supernetworks Isenberg School of Management University of Massachusetts Amherst, Massachusetts 01003
More informationInternational Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: Vol.8, No.7, pp , 2015
International Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: 0974-4304 Vol.8, No.7, pp 99-, 05 Lotka-Volterra Two-Species Mutualistic Biology Models and Their Ecological Monitoring Sundarapandian
More informationA Producer-Consumer Model With Stoichiometry
A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More informationSlope Fields: Graphing Solutions Without the Solutions
8 Slope Fields: Graphing Solutions Without the Solutions Up to now, our efforts have been directed mainly towards finding formulas or equations describing solutions to given differential equations. Then,
More informationThe Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More informationRedoing the Foundations of Decision Theory
Redoing the Foundations of Decision Theory Joe Halpern Cornell University Joint work with Larry Blume and David Easley Economics Cornell Redoing the Foundations of Decision Theory p. 1/21 Decision Making:
More informationVectors. Teaching Learning Point. Ç, where OP. l m n
Vectors 9 Teaching Learning Point l A quantity that has magnitude as well as direction is called is called a vector. l A directed line segment represents a vector and is denoted y AB Å or a Æ. l Position
More informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More informationUsing Markov Chains To Model Human Migration in a Network Equilibrium Framework
Using Markov Chains To Model Human Migration in a Network Equilibrium Framework Jie Pan Department of Mathematics and Computer Science Saint Joseph s University Philadelphia, PA 19131 Anna Nagurney School
More informationScarf Instability and Production: A Simulation
18 th World IMACS / MODSIM Congress, Cairns, Australia 13-17 July 009 http://mssanz.org.au/modsim09 Scarf Instability and Production: A Simulation Fukiharu, T., K.Imamura Faculty of Economics, Hiroshima
More informationMATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation
MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation 1 Bifurcations in Multiple Dimensions When we were considering one-dimensional systems, we saw that subtle changes in parameter
More informationLab 5: Nonlinear Systems
Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting
More informationLecture 20/Lab 21: Systems of Nonlinear ODEs
Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species
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 informationAn Example file... log.txt
# ' ' Start of fie & %$ " 1 - : 5? ;., B - ( * * B - ( * * F I / 0. )- +, * ( ) 8 8 7 /. 6 )- +, 5 5 3 2( 7 7 +, 6 6 9( 3 5( ) 7-0 +, => - +< ( ) )- +, 7 / +, 5 9 (. 6 )- 0 * D>. C )- +, (A :, C 0 )- +,
More informationAN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS
AN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS First the X, then the AR, finally the MA Jan C. Willems, K.U. Leuven Workshop on Observation and Estimation Ben Gurion University, July 3, 2004 p./2 Joint
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 informationEconomics 201b Spring 2010 Solutions to Problem Set 1 John Zhu
Economics 201b Spring 2010 Solutions to Problem Set 1 John Zhu 1a The following is a Edgeworth box characterization of the Pareto optimal, and the individually rational Pareto optimal, along with some
More informationSOLAR MORTEC INDUSTRIES.
SOLAR MORTEC INDUSTRIES www.mortecindustries.com.au Wind Regions Region A Callytharra Springs Gascoyne Junction Green Head Kununurra Lord Howe Island Morawa Toowoomba Wittanoom Bourke Region B Adelaide
More informationM.5 Modeling the Effect of Functional Responses
M.5 Modeling the Effect of Functional Responses The functional response is referred to the predation rate as a function of the number of prey per predator. It is recognized that as the number of prey increases,
More information8 Ecosystem stability
8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 198-8 Construction of generalized pendulum equations with prescribed maximum
More informationComputational Aspects of Aggregation in Biological Systems
Computational Aspects of Aggregation in Biological Systems Vladik Kreinovich and Max Shpak University of Texas at El Paso, El Paso, TX 79968, USA vladik@utep.edu, mshpak@utep.edu Summary. Many biologically
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 informationTowards a theory of ecotone resilience: coastal vegetation on a salinity gradient
Towards a theory of ecotone resilience: coastal vegetation on a salinity gradient Jiang Jiang, Daozhou Gao and Donald L. DeAngelis Appendix A 1.1 The model without diffusion dn 1 = N 1 ρ 1 fs 11 N 1 1
More informationA Probabilistic Upper Bound on Differential Entropy
A Probabilistic Upper on ifferential Joseph estefano Member IEEE and Erik Learned-Miller Abstract A novel, non-trivial, probabilistic upper bound on the entropy of an unknown one-dimensional distribution,
More informationModeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool
Modeling the Immune System W9 Ordinary Differential Equations as Macroscopic Modeling Tool 1 Lecture Notes for ODE Models We use the lecture notes Theoretical Fysiology 2006 by Rob de Boer, U. Utrecht
More informationA Stability Analysis on Models of Cooperative and Competitive Species
Research Journal of Mathematical and Statistical Sciences ISSN 2320 6047 A Stability Analysis on Models of Cooperative and Competitive Species Abstract Gideon Kwadzo Gogovi 1, Justice Kwame Appati 1 and
More informationLotka 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 information2 One-dimensional models in discrete time
2 One-dimensional models in discrete time So far, we have assumed that demographic events happen continuously over time and can thus be written as rates. For many biological species with overlapping generations
More informationRelationships between phases of business cycles in two large open economies
Journal of Regional Development Studies2010 131 Relationships between phases of business cycles in two large open economies Ken-ichi ISHIYAMA 1. Introduction We have observed large increases in trade and
More informationLocal dynamics of holomorphic maps
5 Local dynamics of holomorphic maps In this chapter we study the local dynamics of a holomorphic map near a fixed point. The setting will be the following. We have a holomorphic map defined in a neighbourhood
More informationResearch Article Global Dynamics of a Competitive System of Rational Difference Equations in the Plane
Hindawi Publishing Corporation Advances in Difference Equations Volume 009 Article ID 1380 30 pages doi:101155/009/1380 Research Article Global Dynamics of a Competitive System of Rational Difference Equations
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 informationMATH 415, WEEKS 14 & 15: 1 Recurrence Relations / Difference Equations
MATH 415, WEEKS 14 & 15: Recurrence Relations / Difference Equations 1 Recurrence Relations / Difference Equations In many applications, the systems are updated in discrete jumps rather than continuous
More informationMAXIMA AND MINIMA CHAPTER 7.1 INTRODUCTION 7.2 CONCEPT OF LOCAL MAXIMA AND LOCAL MINIMA
CHAPTER 7 MAXIMA AND MINIMA 7.1 INTRODUCTION The notion of optimizing functions is one of the most important application of calculus used in almost every sphere of life including geometry, business, trade,
More informationPrincipal Secretary to Government Haryana, Town & Country Planning Department, Haryana, Chandigarh.
1 From To Principal Secretary to Government Haryana, Town & Country Planning Department, Haryana, Chandigarh. The Director General, Town & Country Planning Department, Haryana, Chandigarh. Memo No. Misc-2339
More informationLotka-Volterra Models Nizar Ezroura M53
Lotka-Volterra Models Nizar Ezroura M53 The Lotka-Volterra equations are a pair of coupled first-order ODEs that are used to describe the evolution of two elements under some mutual interaction pattern.
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationFeedback-mediated oscillatory coexistence in the chemostat
Feedback-mediated oscillatory coexistence in the chemostat Patrick De Leenheer and Sergei S. Pilyugin Department of Mathematics, University of Florida deleenhe,pilyugin@math.ufl.edu 1 Introduction We study
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 informationFIELDS INSTITUTE COMMUNICATIONS, VOLUME. 21, , 1999 RICH DYNAMICS OF GAUSE-TYPE RATIO-DEPENDENT PREDATOR-PREY SYSTEM
FIELDS INSTITUTE COMMUNICATIONS, VOLUME. 21, 325-337, 1999 RICH DYNAMICS OF GAUSE-TYPE RATIO-DEPENDENT PREDATOR-PREY SYSTEM YANG KUANG abstract Ratio-dependent predator-prey models are increasingly favored
More informationDiscrete Dynamical Systems
Discrete Dynamical Systems Justin Allman Department of Mathematics UNC Chapel Hill 18 June 2011 What is a discrete dynamical system? Definition A Discrete Dynamical System is a mathematical way to describe
More informationPharmacological and genomic profiling identifies NF-κB targeted treatment strategies for mantle cell lymphoma
CORRECTION NOTICE Nat. Med. 0, 87 9 (014) Pharmacoogica and genomic profiing identifies NF-κB targeted treatment strategies for mante ce ymphoma Rami Raha, Mareie Fric, Rodrigo Romero, Joshua M Korn, Robert
More informationOptimizing Stresses for Testing DRAM Cell Defects Using Electrical Simulation
Optimizing Stresses for Testing RAM Cell efects Using Electrical Simulation Zaid Al-Ars Ad J. van de Goor Faculty of Information Technology and Systems Section of Computer Engineering elft University of
More informationFeedback control for a chemostat with two organisms
Feedback control for a chemostat with two organisms Patrick De Leenheer and Hal Smith Arizona State University Department of Mathematics and Statistics Tempe, AZ 85287 email: leenheer@math.la.asu.edu,
More informationCOMPLEXITY OF SHORT RECTANGLES AND PERIODICITY
COMPLEXITY OF SHORT RECTANGLES AND PERIODICITY VAN CYR AND BRYNA KRA Abstract. The Morse-Hedlund Theorem states that a bi-infinite sequence η in a finite alphabet is periodic if and only if there exists
More informationExtinction and the Allee Effect in an Age Structured Population Model
AMS Special Session: Difference Equations and Applications Extinction and the Allee Effect in an Age Structured Population Model Nika Lazaryan and Hassan Sedaghat Department of Mathematics Virginia Commonwealth
More informationWe have two possible solutions (intersections of null-clines. dt = bv + muv = g(u, v). du = au nuv = f (u, v),
Let us apply the approach presented above to the analysis of population dynamics models. 9. Lotka-Volterra predator-prey model: phase plane analysis. Earlier we introduced the system of equations for prey
More informationVector analysis. 1 Scalars and vectors. Fields. Coordinate systems 1. 2 The operator The gradient, divergence, curl, and Laplacian...
Vector analysis Abstract These notes present some background material on vector analysis. Except for the material related to proving vector identities (including Einstein s summation convention and the
More informationQUESTIONS ON QUARKONIUM PRODUCTION IN NUCLEAR COLLISIONS
International Workshop Quarkonium Working Group QUESTIONS ON QUARKONIUM PRODUCTION IN NUCLEAR COLLISIONS ALBERTO POLLERI TU München and ECT* Trento CERN - November 2002 Outline What do we know for sure?
More informationETIKA V PROFESII PSYCHOLÓGA
P r a ž s k á v y s o k á š k o l a p s y c h o s o c i á l n í c h s t u d i í ETIKA V PROFESII PSYCHOLÓGA N a t á l i a S l o b o d n í k o v á v e d ú c i p r á c e : P h D r. M a r t i n S t r o u
More informationNarrowing confidence interval width of PAC learning risk function by algorithmic inference
Narrowing confidence interval width of PAC learning risk function by algorithmic inference Bruno Apolloni, Dario Malchiodi Dip. di Scienze dell Informazione, Università degli Studi di Milano Via Comelico
More informationPrice and Quantity Competition in Dynamic Oligopolies with Product Differentiation
Price and Quantity Competition in Dynamic Oligopolies with Product Differentiation Akio Matsumoto Chuo University Ferenc Szidarovszky University of Arizona Abstract This is a continuation of the earlier
More informationEXISTENCE OF POSITIVE PERIODIC SOLUTIONS OF DISCRETE MODEL FOR THE INTERACTION OF DEMAND AND SUPPLY. S. H. Saker
Nonlinear Funct. Anal. & Appl. Vol. 10 No. 005 pp. 311 34 EXISTENCE OF POSITIVE PERIODIC SOLUTIONS OF DISCRETE MODEL FOR THE INTERACTION OF DEMAND AND SUPPLY S. H. Saker Abstract. In this paper we derive
More informationMARKET AREA UPDATE Report as of: 1Q 2Q 3Q 4Q
Year: 2013 Market Area (Cit, State): Greater Portland, OR MARKET AREA UPDATE Report as of: 1Q 2Q 3Q 4Q Provided b (Compan / Companies): Coldwell Banker Seal, Archibald Relocation, Prudential Northwest
More informationPredator density-dependent prey dispersal in a patchy environment with a refuge for the prey
180 South African Journal of Science 104 May/June 2008 Biological Modelling Predator density-dependent prey dispersal in a patchy environment with a refuge for the prey K Dao Duc a* PAuger ab T Nguyen-Huu
More information