Biology, Medicine, & Natural Product Chemistry ISSN: 2089-6514 Volume 5, Number 1, 2016 | Pages: 9-14 | DOI: 10.14421/biomedich.2016.51.9-14 Local Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid) Dewi Anggreini Lecturer of Study Program Mathematics Education, STKIP PGRI Tulungagung Jl. Mayor Sujadi Timur No. 7 Tulungagung, 66221 Tel/Fax: (0355) 321426 Author correspondency: anggreini_004@yahoo.com Abstract Mathematical model has many benefits in life, especially the development of science and application to other fields. The mathematical model seeks to represent real-life problems formulated mathematically to get the right solution. This research is the application of mathematical models in the field of biology that examines the interaction of the two populations that host populations and parasitoid populations. This study differs from previous studies that examine the interaction of two more species that prey and predators where predators kill prey quickly. In this study the parasitoid population slowly killing the host population by living aboard and take food from the host population it occupies. In this study of differential equations are used to construct a mathematical model was particularly focused on the stability of the local mathematical model of interaction of two differential equations that host and parasitoid populations. Stability discussed in this study are stable equilibrium points are obtained from the characteristic equation systems of differential equations host and parasitoid interactions. Type the stability of the equilibrium point is determined on the eigenvalues of the Jacobian matrix. Analysis of stability is obtained by determining the eigenvalues of the Jacobian matrix around equilibrium points. Having obtained the stable equilibrium points are then given in the form of charts and portraits simulation phase to determine the behavior of the system in the future. Keywords: Mathematical Model, The interaction of two Populations Model, Stability Introduction Mathematics has many benefits in the areas of science and everyday life, especially in the field of Biology. According Widowati and Sutimin (2007:1). The mathematical model is a mathematical field that seeks to represent and describe physical systems or a problem in the real world in a mathematical statement, in order to obtain an understanding of these real world problems to be more precise. Representation of this process is known as a mathematical model. Construction analysis and use of mathematical models is seen as one of the most important applications of mathematics. Have been many studies that discuss the interaction of two species or two populations such as predator-prey or predator – prey. Unlike the predator prey models that describe the interaction where a predator as prey and prey as prey, predators will kill its prey directly in the process of natural selection. While in this study established the relationship between host and parasitoid where the host does not kill quickly, but slowly with the way of life of a ride or take in nutrients from the host they occupy. This research is exciting to do because there are still few studies that discuss mathematical models of host and parasitoid interactions in particular discuss the stability of the system. In this study, given how differential equations can be used to construct a mathematical model was particularly focused on the stability of differential equations mathematical model of the interaction of the two populations that host and parasitoid. Stability discussed in this study are stable equilibrium points are obtained from the characteristic equation systems of differential equations host and parasitoid interactions. Having obtained the stable equilibrium points are then given a simulation to determine the type of stable equilibrium points using matlab program. Type the stability of the equilibrium point is determined on the eigenvalues of the Jacobian matrix. Some types of stability based on the roots of the characteristic equation is as follows. If the value of the roots of the characteristic equation real, distinct, and equally a sign of the system asymptotically stable equilibrium point. If the value of the roots of the characteristic equation complex conjugate purely imaginary but not the asymptotically stable equilibrium point. If the value of the roots of the characteristic equation real, distinct, and not as a sign of the equilibrium point of the system is unstable. If the value of the roots of the characteristic equation purely imaginary ekuilibrim point system is stable but not asymptotically stable. Based on the interaction of host and parasitoid, a mathematical model in the form of differential equations and systems of differential equations is a suitable way to describe the interaction between host and parasitoid. Parasitoid is an organism that spends most of his life history to rely on a single host organism which eventually kills (and often took food) in the process. Parasitoids laid eggs in the body of the target animal after hatching larvae then sucking the bodily fluids of the target animal to death When the host-parasitoid find a place to live, it must be appropriate host physiology 10 Biology, Medicine, & Natural Product Chemistry 5 (1), 2016: 9-14 and nutrition for the successful development of the parasitoid. The goal in conducting this study was to examine how the stability of differential equations of mathematical models the interaction of two populations (host and parasitoids) and examines how the stability of mathematical models simulating the interaction of the two populations of differential equations (host and parasitoid) with matlab program. This research is the development of a mathematical model, especially the model of differential equations formed from the interaction between host and parasitoid. Applications using matlab help to find graph and phase portraits of systems of differential equations stability so they can know the behavior of the system in the future. As study materials for researchers who like or are interested in the development of mathematical models especially differential equations. The interaction of two populations models (Predator- Prey) The basic model Predator-Prey or predator prey interaction of the two populations is given by the equation Lotka-Volterra (1926) *dh mh a hp dt   *dp qp ra hp dt   In this equation ℎ is the prey populations and is a predator population, with t is time. The equation states the prey population changes over time and the equation stated predator population changes over time. So the equation: * b sN a hT (1) With captured prey populations of prey each unit time.   * 1 a h f h bh   (2) With b states the time to deal with a host, f(h) stating the number of hosts attacked parasitoid per unit time or the number of hosts diparasit, A stated figures hosted by parasitoid attack . Equation (2) is Holling type II functional response to the host parasitoid interactions. Dynamic systems In general, the dynamic system is defined as a real problem that is modeled mathematically using differential equations in the equation which contains parameters that are interconnected , as well as changes in the parameters of the equation will lead to changes in the stability of the equilibrium point . The equilibrium point is one of the key concepts in a dynamic system. More general system can be expressed in the following form.̇1 = 1( 1, 2,…, , )2̇ = 2( 1, 2,…, , )⋮̇ = ( 1, 2,…, , ) (3) With =( , ,…, , ), =1,2,…, n is a general function of , =1,2,…, and time is t. The system can be further simplified into a system that is not dependent on the function of time (autonomous system) as followṡ1 = 1( 1, 2,…, )2̇ = 2( 1, 2,…, )⋮̇ = ( 1, 2,…, ) (3) With =( , ,…, , ), =1,2,…, n is a function that does not depend explicitly on the time t . The definition of equilibrium point (Perko, 1991) Point =( , ,…, ) ∈ called the equilibrium point of (3) if ( ̅)̇ = 0, =1,2,…, Stability of Equilibrium Point Definition of Local Stability the equilibrium point ∈ on the system (3) said: 1) Local Stable if for every å>0 there ä>0 such that for every solution ( ) that satisfies ‖ ( )− ̅‖<ä apply ‖ ( )− ̅‖< å for each ≥ . 2) Local asymptotically stable if the equilibrium point∈ stable and there >0 such that for every solution ( ) that satisfies x ‖ ( )− ̅‖< applylim→ ( )= ̅. 3) Do not be stable if the equilibrium point ∈ does not meet 1. Below is given a definition of linear and nonlinear systems granted system(3), With ⊆ and : → continuous function on E. System (3) is said to be linear if , ,…, each linear with respect , ,…, . so the system (3) can be written in the form ̇1 = 11 1+ 12 2+⋯+ 1 ,̇2 = 21 1+ 22 2+⋯+ 2 ,̇ = 1 1+ 2 2+⋯+ , (4) Dewi Aggreini. – Local Stability Analysis of a Mathematical Model of the Interaction of … 11 With ̇ = , continuous on ≤ ≤ , , ,=1,2,…, . Then the system (3) can be expressed in the form ̇ =̇ , with ∈ , ̇ = ( ̇ , ̇ ,…, ̇ ) and A matrix size nxn . Granted system (3) with ⊆ and : → continuous function on E. System (3) is said to be nonlinear if there is i such that nonlinear. The nature of the solution around the equilibrium point of nonlinear systems (3) can be determined through linearization around the equilibrium point of the system. Definition 2.1 Jacobian matrix If the function belongs =( , ,…, ) on the system(3) with ∈ ( ), =1,2,3,…, . Matriks ̅ ̅ … ̅̅ ̅ … ̅⋮ ⋮ ⋱ ⋮̅ ̅ … ̅ Named Jacobian matrix of f in point ̅. (Kocak, 1991). By using Jacobian matrix ( ̅), point stability properties ekuilibrim ̅ it can be seen as long as the point is hyperbolic Eigenvalues and Eigenvectors Formally the definition of eigenvalues and eigenvectors is as follows. Definitions (Anton and Rorres, 2014) Let A nxn matrix and ∈ , ≠0. Vector x is called an eigenvector / characteristic vector of A if=ë for a ë∈ . Numbers ë that satisfies the above equation is called the eigenvalues / value characteristics. Vector x is called an eigenvector corresponding to ë. Stability Analysis of Fixed Point (Boyce & Diprima, 2001) Consider the linear system = with a, b , c and d constants (5). Linear systems (5) commonly called "data model species" in population dynamics. If the matrix A = and such λ are the eigenvalues of matrix A, then: A =A ↔ −ë −ë = 00 Phase portrait and linear systems (Boyce & DiPrima, 2001) Belongs to the following linear systeṁ = (6) Suppose det (A) ≠ 0, so X are the eigenvalues of matrix A , ie ë is the root of the characteristic equationë − =0 (7) Images phase and the system (6) is almost entirely dependent on its eigenvaluesFigure. 2 shows the absorbance graphic of plot toward thin film wavelength of TiO2 that has been soaked in dye with variation of temperature extract. 1. If the eigenvalues different real negative is called a node, all trajectories toward zero, which means no points of zero is stable. 2. If the eigenvalues of positive reals different then called nodal source, all trajectories out of the points become unstable 3. If the eigenvalues of different real opposite in sign, with this so-called saddle point, all trajectories will be away to infinity along the eigenvectors, this resulted critical point will always be unstable. 4. If the eigenvalues of different real opposite in sign, with this so-called saddle point, all trajectories will be away to infinity along the eigenvectors, this resulted critical point will always be unstable. 5. If the eigenvalues equal, with two linearly independent eigenvectors, you will get what is called a star point or propernode, when λ < 0 then the point of criticism would be stable and unstable for λ > 0. 6. If the eigenvalues equal to the eigenvectors will be obtained by so-called improper node, when λ < 0 then the point of criticism would be stable and directions trajectory will go to zero. Whereas for λ > 0 directions trajectory going out leaving the zero point and the point of criticism would be unstable. If the eigenvalues are complex numbers ë± = ± with <0, it will produce a stable behavior called Spiral, 7. If the eigenvalues are complex numbers ë± = ± dengan >0, it will generate a behavior called unstable spiral, all trajectories going out leaving the zero point and the point of criticism would be stable. 12 Biology, Medicine, & Natural Product Chemistry 5 (1), 2016: 9-14 Method This research is a study of Library Studies, which is used in literature studies of relevant literature sources used to gather the necessary information in the study. A literature study by collecting literature sources which may include books, texts, papers and so on. After the collected literature sources ollowed by a review of the literature sources. At the end of the library resources as a basis for analyzing the problem. Mathematical modeling the interaction of two populations of host and parasitoid in this research is needed to determine how the shape of the equation, and then to look for form completion hosted model and the parasitoid. Steps to resolve the problem that is used to analyze the shape and model of the study are as follows: a. Lowering the mathematical model the interaction of two populations of host and parasitoid with Holling response function of type II. b. Analysis of the models begins with finding the equilibrium point of the model systems of differential equations. c. Determining the characteristic equation and the eigenvalues of the Jacobian matrix system which is calculated at each point of equilibrium. d. Then check the stability of the equilibrium point. e. Creating a numerical simulation of the results of the analysis of mathematical models the interaction of two populations of host and parasitoid with Holling response function of type II. Research Result In the discussion this time will be discussed on a mathematical model of interaction between host and parasitoid host. In this model there are two species that interact namely host and parasitoid. The population of the host is assumed to have enough food. Prior to their predators, the intrinsic growth rate (birth rate minus death rate) host proportional to its population. In this case the host becomes the main food for parasitoid so before their host, the parasitoid population growth will decline due to the mortality rate of the parasitoid. Parasitoid attack based on the functional response of the host. In contrast with these interactions, the parasitoid population will be converted to its growth. The ratio of parasitoid immature to adult parasitoids is constant and the population has the maximum growth rate. The parameters given to model the interaction between the host and parasitoid is ℎ( ) which states the host population at time t, ( ) states parasitoid population at the time, (ℎ) stated functional response, stating the effect of oviposition on adult parasitoid population, r states the intrinsic growth rate of the host ,d declared the death rate of the parasitoid . All of these parameters is positive. Functional response (ℎ) assumed to be a function : → , , (ℎ) function with a continuous rise (0)=0, and (ℎ) diferensiabel kontinu. Based on the assumptions above were obtained host and parasitoid interaction model as follows: phfrh dt dh )( ( ( ) ) . dp af h d p dt   (8) With (ℎ)= ∗ a Holling type II functional response to the interaction of host and parasitoid. Here is a system of differential equations after the interaction of the two populations substituted into the equation (8), namely: ℎ= ℎ− ∗ℎ1+ ℎ= ∗ℎ1+ ℎ− (9) From the point Equilibrium Model In this section we will look for a solution to the 2R , because the population must not be negative and the solutions are beyond 2R cannot be interpreted biologically. Equilibrium point equation (9) is obtained if: 0 dt dh (10) 0. dp dt  (11) From equation (10) and (11) was obtained 0 1 *         bh pa rh (12) * 0. 1 a h a d p bh       (13) From equation (12) and (13) was obtained 3 obtained equilibrium point           dbaa ar ph *11 ,0,           0,, *22 dbaa d ph           dbaa ar dbaa d ph **33 ,, . Dewi Aggreini. – Local Stability Analysis of a Mathematical Model of the Interaction of … 13 On the condition that the equilibrium point ∗ −>0 or ∗ > Based on these results it can be concluded that: if∗ > then between host and parasitoid populations can coexis. ∗ > means that the effects of oviposition on adult parasitoid population multiplied by the time available for the parasitoids to chase and finds a host is greater than the mortality rate of parasitoid multiplied by time to handle one host . If ∗ > then between host and parasitoid populations cannot coexist. Stability of Equilibrium Point One way to determine the stability of equilibrium points of the system (9) is to use linearization. Based Systems (9), for example ℎ= ( , ) and = ( , ) with,   p bh ha rhyxg         1 , * 1 and  yxg ,2 pdbh ha a        1 * . obtainable       * * 1 1 2 * * 2 2 2 11 , . 11 a p a hg g r bhbhh p J h p g g a ap a ah d h p bhbh                           Local Stability Analysis In the following theorems will discuss the local stability of equilibrium points:(ℎ , )= ∗ − , ∗ −(ℎ , )= 0, ∗ and (ℎ , )= ∗ ,0 Theorem If ∗ > , ∗ > ( + ), ∗∗ > − ,∗∗ > and ∗∗ + − < ∗∗ − then the equilibrium point (ℎ , )= 0, ∗ local asymptotically stable. Theorem If ∗ > , ∗ > ( + ) and( ∗ )∗ − < ( ∗ )∗ , and ( ∗ )∗ > then the equilibrium point           dbaa ar dbaa d ph **33 ,, Local asymptotically stable. Theorem If 0 d d r a         and 2 4 d rd d r rd a a                then the equilibrium point  2 2 *, , 0 d h p a a db      unstable. Conclusion The following will be given some conclusions related to this research. 1. From the system of differential equations interaction of two species , namely: = ℎ− ∗ and = ∗ − obtained equilibrium points(ℎ , )= 0, ∗ , (ℎ , )= ∗ , ∗ and (ℎ , )= ∗ ,0 . 2. If ∗ > , ∗ > ( + ), ∗∗ > − ,∗∗ > and ∗∗ + − < 4 ∗∗ , , then the equilibrium point (ℎ , )= 0, ∗ local asymptotically stable. 3. If ∗∗ + − <4 ∗∗ , ,∗∗ > − , and ∗∗ > thenë , = ∗∗ ± ∗∗ ∗∗ have a negative real part . So based on the equilibrium point theorem (ℎ , )= 0, ∗ local asymptotically stable. 4. By Theorem if: ∗ > , ∗ > ( + ),∗∗ > − , ∗∗ > and ∗∗ + −< 4 ∗∗ − , then at the beginning of time very small number of host population , while close to the parasitoid population ∗ . So for a long time, the parasitoid population will approach ∗ while the host population will become extinct. 14 Biology, Medicine, & Natural Product Chemistry 5 (1), 2016: 9-14 5. If ∗ > , ∗ > ( + ), and ( ∗ )∗ −< ( ∗ )∗ , and ( ∗ )∗ > then the equilibrium point (ℎ , )= ∗ , ∗ local asymptotically stable. Because (ℎ , )= ∗ , ∗ have negative real part, the equilibrium point (ℎ , )= ∗ , ∗ local asymptotically stable. By Theorem if: ( ∗ )∗ − − ( ∗ )∗ , ∗> , ∗ > ( + ), and ( ∗ )∗ > , then at the beginning of time close to the host population ∗ while close to the parasitoid population∗ . So for a long time, will host population approaching ∗ and the parasitoid population will approach ∗ . If + + >0 and+ + >4 + then the equilibrium point (ℎ , )= ∗ ,0 unstable. References Anton, H. and Rorres, C., 2014, Elementary Linear Algebra, John Wiley & Sons, Inc., New York. Boyce, W.E. & DiPrima R.C. 2001. Elementary Differential Equations and Boundary Value Problems. Seventh Edition, John Wiley & Sons: New York. Edelstein, L., Keshet, 2005, Mathematical Models in Biology, Siam, New York. Haberman, R. 1977. Mathematical Models in Mechanical Vibrations, Population Dynamics, and Traffic Flow. Prentice- Hall: New Jersey. Hanh, W., 1967, Stability of Motion, Springer-Verlag, New York. Hale, J. K. & H. Kocak. 1991. Dynamics and bifurcations. New York: Springer Verlag. Kuznetsov, Y.A., 1998, Elements of Applied bifurcation Theory, Springer-Verlag, New York. Ni'mah, Khoirun, 2015, Analysis of Model S-I-P Interaction of Two Species Predator-Prey With Response Function Holling Type II, Semarang State University. Perko, L., 2001. Differential Equations and Dynamical Systems, Springer-Verlag, New York. Sulistyowati, E. et al., 2001, Functional Response Parasitoid Cephalonomia Stephanoderis Betr Against Coffee Fruit powder, Hypothenemus hampei Ferr, Pest Plant Science Program, University of Gadjah Mada. Tarumingkeng, R.C., Population Dynamics (Quantitative Ecological Assessment), 1994, Pustaka Sinar Harapan and Activities of Discourse Christian University, Jakarta.