On modeling two immune effectors two strain antigen interaction

In this paper we consider the fractional order model with two immune effectors interacting with two strain antigen. The systems may explain the recurrence of some diseases e.g. tuberculosis (TB). The stability of equilibrium points are studied. Numerical solutions of this model are given. Using integer order system the system oscillates. Using fractional order system the system converges to a stable internal equilibrium. Ulam-Hyers stability of the system has been studied.


Introduction
Immune system (IS) is known to be multifunctional and multi-pathways i.e. a given function is performed by more than one effector. And each effector, typically perform more than one function [1]. This guarantees to a great extent the resilience of the immune network [2]. Also many antigens evolve with time hence they are multi strains. This explains why some diseases re-appear e.g. tuberculosis (Tb). Therefore modeling the interaction of two immune effectors with two strain antigen is an important problem.
The use of fractional-orders differential and integral operators in mathematical models has become increasingly widespread in recent years [3]. Several forms of fractional differential equations have been proposed in standard models.
Differential equations of fractional order have been the focus of many studies due to their frequent appearance in various applications in fluid mechanics, economic, viscoelasticity, biology, physics and engineering. Recently, a large amount of literatures developed concerning the application of fractional differential equations in nonlinear dynamics [3].
In this paper we study the fractional-order model with two immune effectors interacting with two strain antigen. In sec.2 we present the fractional-order model, study their equilibrium and their local stability and solve it numerically. In sec.3 the Ulam-Hyers stability is presented. In sec. 4 our conclusions are presented. Now we give the definition of fractional-order integration and fractional-order differentiation: Definition 1 The fractional integral of order b R + of the function f (t), t >0 is defined by and the fractional derivative of order a (n − 1, n) of f (t), t >0 is defined by The following properties are some of the main ones of the fractional derivatives and integrals (see [3][4][5][6][7][8][9]).
Let b, g R + and a (0, 1). Then The following lemma can be easily proved (see [7]).
The fractional-order model Let x 1 , x 2 be two strains of an antigen and y 1 , y 2 be two immune effectors then the fractionalorder IS model is given by: where 0 < a ≤ 1 and a 1 , a 2 , b 1 , b 2 , c 1 , c 2 , d 1 , d 2 are positive constants. The constants c 1 , c 2 are the mutation rates of the antigen strains.
There are several equilibria e.g. the zero equilibrium To study the stability we need the Jacobian matrix of (3) given by: It is direct to see that the zero solution is unstable (notice that by definition a 1 > > c 1 , a 2 > > c 2 The numerical simulations of (3) are given in Figures  12345678 for a 1 = 1.0, a 2  In Figure 1 we take a = 1.0. In Figure 2 we take a = 0.9. In Figure 3 we take a = 0.85. In Figure 4 we take a = 0.8. In Figure 5 we take a = 1.0. In Figure 6 we take a = 0.9. In Figure 7 we take a = 0.85. In Figure 8 we take a = 0.8. The relation between the two strains of an antigen x 1 (t) and x 2 (t) are given in Figures 1234 for different 0 < a ≤ 1. The relation between the two immune effectors y 1 (t) and y 2 (t) are given in Figures 5678 for different 0 < a ≤ 1.
The Figures 12345678 show that by using the integer order system (a = 1) the system oscillates (Figures 1  and 5) and by using the fractional order system (0 < a <1) the system converges to a stable internal equili-  . Mathematically the following system can be studied ( [10], [11]): the system (5) is Ulam-Hyers (UH) stable if it has an exact solution and if ∀ε >0 there is δ >0 such that if x a (t) is an approximation for the solution of (5) then there is an exact solution x(t) of (5) which is close to x a i.e., This definition has applicable significance since it means that if one is studying an UH stable system then one does not have to reach the exact solution (which usually is quite difficult or time consuming). All what is required is to get a function which satisfies (6). UH stability guarantees that there is a close exact solution. This is quite useful in many applications e.g. numerical analysis, optimization, biology and economics etc., where finding the exact solution is quite difficult. It also helps, if the stochastic effects are small, to use deterministic model to approximate a stochastic one.
We begin by realizing that UH stability is independent of the more familiar Lyapunov stability which states that the system (5) is Lyapunov stable if both x(t), y(t) are exact solutions of (5) and for all ε >0 there is δ >0 such that | x(0) − y(0) |< δ implies | x(t) − y(t) |< ε for all t >0.
A known counter-example proving the independence of the two concepts is the system: dx dt ax t a = > ( ), constant 0 (7) whose x = 0 solution is Lyapunov unstable while it is UH stable [11]. UH stability has been studied for functional equations [12], and linear differential equations [13]. Now we study local UH stability for nonlinear systems. Consider systems (5), (6), assume