• Submissions

# Epidemiological Nonlinear Modeling of Middle East Respiratory Syndrome Coronavirus in Saudi Arabia Using Homotopy Perturbation Method

1Department of Mathematics, AMET Deemed to be University, Chennai, Tamilnadu, India.

2Department of Science and Humanities, Hindusthan College of Engineering and Technology, Coimbatore, Tamilnadu, Chennai.

*Corresponding author: Rajendran.L, Department of Mathematics, AMET Deemed to be Univeristy, Chennai, India.

Submission: February 17, 2022; Published: March 15, 2022

DOI: 10.31031/COJRR.2022.03.000574

ISSN 2639-0590
Volum3 Issue5

#### Abstract

Middle East Respiratory Syndrome (MERS) is a viral respiratory disease caused by a coronavirus transmitted to humans by infected dromedary camels. It is a zoonotic virus transmitted between animals and people, and it is contactable through direct or indirect contact with infected animals. The very first case remained isolated in the Kingdom of Saudi Arabia (KSA) in 2012, and after that, MERS-covid exhibited outbreaks in several regions of the world. In this communication, we propose a mathematical model to analyze the transmission rate of this virus and remedial measures to reduce the transmission rate by applying the Homotopy Perturbation Method.

Keywords: Mathematical modeling; Middle east respiratory syndrome; Nonlinear differential equations; Homotopy perturbation method

#### Introduction

Since 2012, many countries have reported cases of MERS, including Algeria, Austria, Bahrain, China, Egypt, France, Germany, Greece, Islamic Republic of Iran, Italy, Jordan, Kuwait, Lebanon, Malaysia, the Netherlands, Oman, Philippines, Qatar, Republic of Korea, India, Kingdom of Saudi Arabia, Thailand, Tunisia, Turkey, United Arab Emirates, United Kingdom, United States, and Yemen [1-3]. This virus’s most common symptoms include fever, severe cough and shortness of breath. Some people also had diarrhea and nausea [4-8].

The virus uses its spike proteins as pro inflammatory adhesive factors for entering into the host through a specific receptor called Dipeptidyl Peptidase-4 (DPP4). This receptor is considered a critical factor in the motioning and activation of the developed and essential immune responses in infected patients. Using potent antigens in combination with solid adjuvant may effectively trigger the activation of specific MERS- covid cellular responses as well as the production of neutralizing antibodies. Regrettably, to date, there is no effective approved treatment or vaccine for MERS-covid [2,3]. The means of transmission to humans is not completely defined yet. Still, evidence proposes that MERS-covid can be transmitted to humans via both animals and humans, spreading among people in close contact. Transmission has also been observed from infected patients to healthcare personnel. Camels and bats are known as natural hosts, and evidence of camel to human transmission has been reported [4].

Corona virus’s virions are roughly spherical with a moderate degree of Pleomorphism [5]. The cube-like or petal-shaped spikes emerge from the virion surface as stalks with a bulb-like distal terminus. The mean particle diameter is 118- 136nm, including the contributions of the spikes, which projects 16-21nm from the virion envelope [6,7]. Coronaviruses have helically symmetric nucleocapsids, which consist of a canonical set of 4 major structural proteins: the Spike (S), Membrane (M) and Envelope (E), which are located in the membrane envelope and the Nucleocapsid (N) protein which is found in the ribonucleoprotein core.

A vast number of mathematical models have been proposed to understand the dynamics of the coronaviruses. Many researchers study the SIS, SIR and SEIR models to investigate the dynamics of the corona viruses [8-12]. This paper’s proposed a model that investigates the dynamics of the MERS-covid via the Homotopy Perturbation method [13-17] (Figure 1).

Figure 1: Schematic diagram of the MERS-Covid.

#### Mathematical Formulation of the Problem

Several models have been projected as tools to recognize epidemics. This section focuses on the SEIR model framework for camel-to-camel transmission [8]. We use the SEIR model to understand the spread of infectious diseases. The SEIR model is one of the epidemiological models. It groups the overall population into distinct sub-populations: a susceptible population S, an exposed population E, an infected population I, and a recovered population R. We make several further simplifying assumptions. We describe the epidemic model of nonlinear differential equations as follows [13]:

Initial conditions are [13]

S (0) = s0=0.79508, E (0) =e0=0.003552, I (0) =i0=0.003552, R (0) =r0=0.197815 (5)

where β is the transmission rate, μ is the natural death and birth rates, N is the population size of camels, σ and γ are rate at which infected camels change status from exposed to infectious and from infectious to recovered respectively. λ is the rate at which recovered camels lose immunity protection due to multiple factors, such as the evolution of the virus, natural loss of immunity [14,15], imports of susceptible juvenile camels, and exports or consumption of recovered adult camels. Susceptible camels could also be exported and consumed for meat. We assumed that the replacement of susceptible camels would be immediate, we also assumed no export or consumption of infectious camels, which is acceptable if the frequency, is less than 10%. The total population size of dromedary camels is assumed to be constant at 2, 70, 000 and the study period went from 1 January 2014 to 31 May 2016 [13,18].

#### Result and Discussion

So many authors have applied the stated method to various non-linear problems and demonstrated the efficiency of the HPM. Using the Homotopy Perturbation Method (refer to Appendix A), we obtained the approximate solutions of the eqns. (1)-(4) as follows:

#### Numerical Simulation

We execute numerical simulations of MERS model (1)-(4), for which we use a set of logical parameter values. The simulation is performed using graphmain3 solver written in MATLAB programming. Eqns. (6)-(9) is very closed expressions for the models Eqns. (1)-(4). We compared our analytical solutions obtained by HPM with the numerical simulation obtained by MATLAB software. To test the accuracy of our approximate analytical expressions, the system (1)-(4) is also solved numerically by using MATLAB for all possible values of parameters. The MATLAB program is also given in Appendix B.

The Homotopy Perturbation method (HPM) was applied to obtain the analytical solution of MERS-covid among dromedary camels in Saudi Arabia, given in the Mathematical model of the first-order nonlinear differential equation. (Figure 2-6) shows the comparison of analytical and numerical solutions. Figure 3 shows that the number of sensitive cases increases when (N) increases. In addition, we can study that HPM provides better accuracy and effective results for various values of parameters. Furthermore, we may wrap up a significant decrease in the spread of the virus when direct or indirect contact with infected animals decreases.

Figure 2: Comparison of our analytical expression of SEIR with the numerical result for the various values of the parameters using eqs.(6)-(9).

Figure 3: Plot of our analytical expression of susceptible for the various values of the parameters, when s0=0.79508, e0=0.003552, i0 0.003552, r0=0.197815, β=200, and γ=3 using equation (6).

Figure 4: Plot of our analytical expression of exposed for the various values of the parameter when s0=0.79508, e0=0.003552, i0=0.003552, r0=0.197815, β=200, and μ=3 using eq.(7).

Figure 5: Plot of our analytical expression of infected for the various values of the parameters, when s0=0.79508, e0=0.003552, i0=0.003552, r0=0.197815, β=200, γ=0.3 and μ=0.1 using eq.(8).

Figure 6: Plot of our analytical expression of infected for the various values of the parameters, when s0=0.79508, e0=0.003552, i0=0.003552, r0=0.197815, β=200, γ=0.25 and μ=0.09 using eq. (9).

#### Conclusion

This study examined the transmission rate of MERS-covid and preventive measures using the SEIR model with minimizing the direct and indirect contact of infected camels as a remedial action from the ODE system using the Homotopy Perturbation Method. Based on the information presented above, we can observe that the waves followed the waves in camel cases in humans, which implies that camel cases have led to human cases, and by taking adequate measures to identify, diagnose and isolate infected dromedary camels, the transmission rate of MERS-covid may be reduced.

Appendix A: Basic concept of Homotopy Perturbation method (HPM)

To illustrate the basic ideas of this method, we consider the following nonlinear functional equation:

With the following boundary condition:

where A is a general functional operator, B a boundary operator, f (r) is a known analytical function and Γ is the boundary of the domainΩ . The operator A can be decomposed into two operators L and N, where L is linear, and N is nonlinear operator.

Eqn. (19) can be, therefore, written as follows:

Using the Homotopy technique, we construct a Homotopy which U(r, p) :Ω×[0,1]→ R , satisfies:

where p∈[0,1] is an embedding parameter, 0 u is an initial approximation for the solution of Eqn. (SA2), which satisfies the boundary conditions. Obviously, from Eqns. (A4) and (A5) we will have:

The changing values of p from zero to unity are just that of U(r, p) from u0(r) to u(r) .In topology, this is called homotopy. According to HPM, we can first use the embedding parameter p as a small parameter and assume that the solution of Eqns. (A4) and (S5) as a power series in p:

Setting p = 1, results in the approximation to the solution of Eqn. (A8)

The combination of the Perturbation method and the Homotopy method is called the HPM; which has eliminated limitations of the traditional perturbation techniques. The series Eqn. (A9) is convergent for more cases (Table 1).

Appendix B: Analytical solution on nonlinear Eqns. (1)- (4) using Homotopy Perturbation Method (HPM).

Consider the differential equations

Initial conditions are

S (0) =0.79508, E (0)=0.003552, I(0)=0.003552, R(0)=0.197815 (B5)

The Homotopy form the Eqns. (B1) -(B4) can be constructed as follows:

Where p is the embedding parameter and p∈[0,1] ,The approximate solution of (B5) and (B6) are

Substituting Eqns. (B10) and (B11) into (B6) - (B9), Comparing the coefficients of like powers of p in Eqn. (B6) gives:

The initial conditions are,

The solution of the Eqns. (B12) and (B13) are given by

where

Comparing the coefficients of like powers of p in Eqn. (B7) gives:

The initial conditions are,

The solution of the Eqns. (B17) and (B18) are given by

Comparing the coefficients of like powers of p in Eqn. (B8) gives:

The initial conditions are,

The solution of the Eqns. (B22) - (B23) are given by

Comparing the coefficients of like powers of p in Eqn. (B9) gives:

The initial conditions are,

The solution of the Eqns. (B27) and (B28) are given by

Table 1: Nomenclature and Units.

Appendix C. Matlab program for the numerical solution of differential Eqns. (1) – (4)

function graphmain3
options= odeset(‘RelTol’,1e-6,’stats’,’on’);
X0=[1;1;1;1];
tspan=[0,3];
tic
[t,X]=ode45(@TestFunction,tspan,X0,options);
toc
figure
hold on
plot(t,X(:,1),’*’);
plot(t,X(:,2),’.’);
plot(t,X(:,3),’+’);
plot(t,X(:,4),’+’);
legend(‘x1’,’x2’,’x3’,’x4’,’x5’)
ylabel(‘x’)
xlabel(‘t’)
return
function [dx_dt] =TestFunction(t,x)
m=2; r=3; beta=0.01; N=2; y=0.1; e=1/2.
dx_dt(1) =m*N+y*x(4)-beta*x(1)*x(3)-m*x(1);
dx_dt(2) =beta*x(1)*x(3)-e*x(2)-m*x(2);
dx_dt(3) =e*x(2)-r*x(3)-m*x(3);
dx_dt(4) =r*x(3)-y*x(4)-m*x(4);
dx_dt = dx_dt’;
return

#### References

We at Crimson Publishing are a group of people with a combined passion for science and research, who wants to bring to the world a unified platform where all scientific know-how is available read more...