# Mathematical Analysis of the Prey-Predator System with Immigrant Prey Using the Soft Computing Technique

Authors: Naveed Ahmad Khan; Muhammad Sulaiman; Jamel Seidu; Fahad Sameer Alshammari
Journal: Discrete Dynamics in Nature and Society
 (2022)

---

## Abstract

In this paper, a mathematical model for the system of prey-predator with immigrant prey has been analyzed to find an approximate solution for immigrant prey population density, local prey population density, and predator population density. Furthermore, we present a novel soft computing technique named LeNN-WOA-NM algorithm for solving the mathematical model of the prey-predator system with immigrant prey. The proposed algorithm uses a function approximating ability of Legendre polynomials based on Legendre neural networks (LeNNs), global search ability of the whale optimization algorithm (WOA), and a local search mechanism of the Nelder–Mead algorithm. The LeNN-WOA-NM algorithm is applied to study the effect of variations on the growth rate, the force of interaction, and the catching rate of local prey and immigrant prey. The statistical data obtained by the proposed technique establish the effectiveness of the proposed algorithm when compared with techniques in the latest literature. The efficiency of solutions obtained by LeNN-WOA-NM is validated through performance measures including absolute errors, MAD, TIC, and ENSE.


---

## Body

## 1. Introduction

A prey-predator system is one of the prevalent phenomena in nature. The interaction between predators and prey in any environment was first introduced by Lotka and Volterra in 1926 [1]. Holling in 1966 [2] elaborated prey-predictor models with different kinds of functional responses for predation. All these models are inspired by biological phenomena and presented by nonlinear ordinary and partial differential equations. Danca et al. [3] study the detailed analysis of a nonlinear prey-predator model. A connection between predators and prey has a long history and will continue as a governing theme in biomathematics because of its universal significance [4]. Researchers have made many changes by introducing different facets to the predator-prey system such as delay in predator growth, harvesting prey and predators, and providing additional food to a predator for sustaining prey population and prey diseases [5–7]. Cai et al. [8] studied the dynamical system of prey-predator with Allee effects in prey growth. A stage-structured predator-prey model with gestation delay is investigated by [9]. In [10, 11], bifurcation, chaos and dynamic behavior of the nonlinear discrete-time predator-prey system is studied. Huang et al. [12] study the stability analysis of the model with consideration of the prey refuge. Studies by Din [13], Weide et al. [14], and Gong et al. [15, 16] are presented for the discrete-time nonlinear prey-predator types of model. Hadeler and Freedman [17] provide extensive references to real-world examples of three-species ecoepidemiological systems of sound prey, infected prey, and predators.In recent years, there has been a substantial increase in the amount of attention paid to population biology by scientists due to the significant applications it has in ecology. It bridges the gap between mathematics and biology. The dynamical systems in biology have been investigated to interpret various problems. The Lotka–Volterra population model is a well-known mathematical model that describes biological systems [18]. In real life, there is a strong correlation between size, age, and developmental stages of different populations of species. It is an important strategy to incorporate all these variables in the mathematical modeling of populations of different spices. To develop more realistic and accurate mathematical models, many scientists have suggested noise-induced models [19, 20] and spatial models [21]. In biomathematical problems, researchers utilize mathematical modeling and simulations along with biological structures to describe the phenomena by the system of nonlinear differential equations. These nonlinear models are considered stiff and unrealistic, and therefore, finding exact and semianalytical solutions for such problems is challenging because of nonlinearity. Analytical approaches for solving nonlinear problems are mostly based on Laplace or Fourier transformation, Laguerre’s integral formula, and the Grunwald–Letnikov concept [22, 23]. When tackling complex problems, these techniques may be challenging to use; in addition, the solution is provided in a closed form that necessitates the evaluation of special functions using complex expressions, such as the Mittag-Leffler function [24].In recent years, researchers have been working to develop new techniques for finding approximate solutions to nonlinear models; e.g., the Laplace Adomian decomposition method [25], the new coupled fractional reduced differential transform method [26], the Runge–Kutta–Fehlberg method [27], the finite element method [28], the Sumudu decomposition method [29], the implicit Adams methods [30], the confidence domain technique [31], and the homotopy analysis method [32] have become much more significant to get accurate solutions. All these deterministic approaches have their own advantages, applicability, and drawbacks. With great interest, it is noted that such techniques are gradient-based and call for information about the problem beforehand. The availability of several local optima, which leads to solutions where global optimality cannot be easily ensured, is one of the fundamental limitations of gradient-based approaches. Global optimality is sought in gradient-based approaches by randomly scanning the design space from various starting points. However, this causes the technique to become sluggish and computationally inefficient for complex nonlinear optimization problems [33, 34]. In addition, the Runge–Kutta methods are self-starting and stable techniques that can be easily implemented to calculate the solution to different problems. The main drawbacks of the Runge–Kutta methods are that they take longer time to calculate solutions than other multistep methods with equivalent precision, and it is difficult for them to get accurate global estimates of the truncation error [35]. To overcome these drawbacks, a stochastic metaheuristic approach with artificial neural networks is developed, which is free of a gradient and does not require any prior information about the problem. The majority of metaheuristic techniques are inspired by natural, physical, or biological processes and make use of a variety of operators to mimic the fundamental behavior. The harmony between exploration and exploitation is a recurring subject in all metaheuristics.In recent times, artificial neural network (ANN)-based stochastic algorithms with global and local search optimizers have been designed to solve differential equations representing physical phenomena including flow in a circular cylindrical conduit via electrohydrodynamics [36], a model of an immobilized enzyme system that follows the Michaelis–Menten (MM) kinetics for a microdisk biosensor [37], flow of Johnson–Segalman fluid on the surface of an infinitely long vertical cylinder [38], and beam-column designs [39]. The abovementioned techniques motivate authors to design a new soft computing algorithm, the LeNN-WOA-NM algorithm, to find approximate series solutions using Legendre polynomials for the model presenting the prey-predator system with immigrant prey. The salient features of the paper are summarized as follows:(i)
A mathematical model for a prey-predator system with immigrant prey is formulated and analyzed to study the influence of variations on the growth rate, force of interaction, and the catching rate of local and immigrant prey.(ii)
Artificial neural network-based weighted Legendre polynomials are used to construct the model of approximate solutions for the prey-predator model. A fitness function based on mean square errors is designed to assess unknown parameters with the help of global search ability of whale optimization and a local search mechanism of the Nelder–Mead algorithm.(iii)
The suggested technique can result in adequate solutions for nonlinear hard problems for which no exact algorithm exists that can solve them in a reasonable amount of time.(iv)
Performance of the proposed algorithm is validated in terms of absolute errors, mean absolute deviation, Theil’s inequality coefficient, Nash–Sutcliffe efficiency, and error in Nash–Sutcliffe efficiency.(v)
Approximate solutions, convergence of fitness, and performance measures obtained by the LeNN-WOA-NM algorithm for prey-predator systems with immigrant prey are shown through different graphs and tables, which shows the dominance and robustness of the proposed algorithm in solving real-world problems.(vi)
Unlike most classical methods, the LeNN-WOA-NM algorithm requires no gradient information and therefore can be used with nonanalytic, black-box, or simulation-based objective functions to approximate complex problems.


## 2. Problem Formulation

Goteti developed the mathematical model for a prey-predator system with immigrant prey to understand the interaction and communication between predators and immigrant prey. The model is assumed to follow mass action theory and consists of prey population density, which is defined as follows:(1)Nt=Xt+St,where St and Xt denote the population density of local and immigrant prey, respectively. Population density of the predator is denoted by Yt.The following assumptions are used in formulating the mathematical model for the prey-predator system with immigrant prey:(i)
In the presence of a predator, the population of prey is classified into two subcategories named local preySt and immigrant prey Xt.(ii)
In absence of the predator, the population growth of local prey logistically increases with an intrinsic growth rateα1 with environmental carrying capacity denoted by k1.(iii)
With availability of the local and immigrant prey population, the population of the predator grows logistically with growth ratesc1 and c2, while suffering loss of populations is denoted by μ1 and μ2.(iv)
Immigrant and local prey can reproduce, and therefore, it is assumed that birth rates should be positive. The growth rate of immigrant prey increases at the rateα2 with environmental carrying capacity denoted by k2.(v)
It is assumed that immigrant prey is a natural choice of predators.β1 and β2 denote the positive and negative force of interaction between local and immigrant prey.(vi)
It is assumed that the local prey and immigrant prey are caught by the predator at the rate ofγ1 and γ2, respectively.A mathematical model for a prey-predator system with immigrant prey is given by the following system of differential equations:(2)dSdt−Sα1−α1Sk1−β1SX+γ1SY=0,dXdt−Xα2−α2Xk2+β2SX+γ2XY=0,(3)dYdt−c1SY−c2XY+μ1Y+μ2Y2=0,with initial populations(4)S=S0,X=X0,Y=Y0att=0.


## 3. Approximate Solutions and Weighted Legendre Polynomials

The Legendre polynomials are denoted byLnt, where n denotes the order of Legendre polynomials. These polynomials constitute the set of orthogonal polynomials on −1,1. The first eleven Legendre polynomials are given in Table 1. High-order Legendre polynomials are generated by the following recursive formula:(5)Ln+1r=1n+12n+1rLnr−nLn−1r.Table 1
Eleven Legendre polynomials with an independent variablet.
nLnt011t21/23t2−131/25t3−3t41/835t4−30t2+351/863t5−70t3+15t61/16231t6−315t4+105t2−571/16429t7−693t5+315t3−35t81/1286435t8−12012t6+6930t4−1260t2+3591/12812155t9−25740t7+18018t5−4620t3+315t101/25646189t10−109395t8+90090t6−30030t4+3465t2−63We consider an approximate series solution for equations (2)–(4) representing the prey-predator model with immigrant prey as follows:(6)Sapproxt=∑n=134ζnLnψnt+θn,Xapproxt=∑n=3568ζnLnψnt+θn,Yapproxt=∑n=69102ζnLnψnt+θn,where ζn, ψn, and θn are unknown parameters.Sincenth order continuous derivatives of system (6) exist, we consider the first derivative of system (6) as follows:(7)ddtSapprox=∑n=134ζnLn′ψnt+θn,ddtXapprox=∑n=3568ζnLn′ψnt+θn,ddtYapprox=∑n=69102ζnLn′ψnt+θn,and we plug equations (6)–(7) in governing ordinary differential equations. Equations (2)–(4) will be transformed into an equivalent algebraic system of equations that can be solved for unknown parameters ζn, ψn, and θn using the LeNN-WOA-NM algorithm.


## 4. Fitness Function Formulation

The mean square error (MSE)-based fitness function for solving the prey-predator model for equations (2)-(4) can be written as follows:(8)Minimizeε=ε1+ε2+ε3+ε4+ε5+ε6,where ϵ1 to ϵ6 are the fitness-based errors for equations (2)–(4) along with the initial conditions and are given as follows:(9)ε1=1N∑n=134dSndt−Snα1−α1Snk1−β1SnXn+γ1SnYn2,ε2=1N∑n=3568dXndt−Xnα2−α2Xnk2+β2SnXn+γ2XnYn2,ε3=1N∑n=69102dYndt−c1SnYn−c2XnYn+μ1Yn+μ2Yn22,ε4=S0−0.52,ε5=X0−0.52,ε6=Y0−0.52.


## 5. Optimization Network

### 5.1. Whale Optimization Algorithm

The whale optimization algorithm is a nature-inspired metaheuristic algorithm designed by Mirjalili and Lewis [40]. The working strategy of WOA is inspired by foraging behavior of humpback whales. The humpback whales chase prey or krill by swimming around them in a molded way as shown in Figure 1. The mathematical model of each phase is explained below.Figure 1
Unique bubble-net feeding method and spiral moment for updating the position of humpback whales.


#### 5.1.1. Exploration Phase

Humpback whales encircle prey for hunting. Equations (10) and (11) mathematically model this behavior as follows:(10)E⟶=C⟶⋅Z⟶∗t−Z⟶t,(11)Z⟶∗t+1=Z⟶∗t−A⟶⋅D,where t denotes the current iteration, Z∗ has provided the best solution so far, and D gives the location of humpback whales to prey at each step. A and C are coefficient vectors which are defined as follows:(12)A⟶=2a⟶⋅r⟶−a⟶,(13)C⟶=2⋅r⟶.where r∈0,1 is an arbitrary nominated vector. The value of a reduces from two to zero during exploration as well as in the exploitation phase. During helix-shaped movement, the distance between prey and the humpback whale is given as follows: (14)Z⟶t+1=D′⋅ebl⋅cos2πl+Z⟶∗t,where D′ represents the location of the ith whale to the prey and is defined as follows:(15)D′=Z⟶∗t−Z⟶t.Furthermore,b denotes the state of the logarithmic helix and l∈−1,1 is any arbitrary number. The shrinking surrounding technique of humpback whales during contracting loop is summarized as follows:(16)Z⟶t+1=Z⟶t−A⟶⋅Difp≤0.5,D⟶′⋅ebl⋅cos2πl+Z⟶∗tifp≥0.5.


#### 5.1.2. Exploitation Phase

In the exploration phase (searching for prey), the heterogeneity of the vectorA is utilized. If A>1, the position of the search agent is updated, and the entire mechanism is modeled by the following equations:(17)D⟶=C⟶⋅Z⟶rand∗−Z⟶,(18)Z⟶t+1=Z⟶rand−A⟶⋅D⟶,where Zrand⟶ is an arbitrary position vector taken from the current population.Figure2 shows the flowchart of the WOA. It can be seen that the WOA creates a random initial population from candidate space and evaluates it using an error-based fitness function when the optimization process starts. After finding the best solution, the algorithm repeatedly executes the following steps until ending criterion is achieved.Figure 2
Graphical illustration of the working steps of the whale optimization algorithm and Nelder–Mead algorithm for the training of neurons in the LeNN architecture and the minimization of fitness functions.


### 5.2. Nelder–Mead Algorithm

A Nelder–Mead (NM) algorithm is a direct search method also known as a downhill simplex method developed by Nelder and Mead in 1965 to solve different problems without any information about the gradient [41]. NM is a single path following a local search optimizer that can find good results if initialized with a better initial solution. A simplex consisting of n+1 vertices is set up to minimize a function f with dimensions n [42]. The NM algorithm generates a sequence of simplices by following four basic procedures, namely, reflection, expansion, contraction, and shrink. Further details about the NM algorithm can be found in [43]. Figure 2 shows the working procedure of the NM algorithm.In recent times, the most successful and effective trend in optimization is the action of integrating components from different methods. The foremost motivation behind the hybridization of diverse algorithmic ideas is to acquire better performing systems, which exploit and coalesce benefits of different techniques. Therefore, in this study, we have combined the global and local search optimization algorithms to achieve more efficient and robust solutions. The hybridized working procedure of the designed LeNN-WOA-NM algorithm is illustrated in Figure2.


## 5.1. Whale Optimization Algorithm

The whale optimization algorithm is a nature-inspired metaheuristic algorithm designed by Mirjalili and Lewis [40]. The working strategy of WOA is inspired by foraging behavior of humpback whales. The humpback whales chase prey or krill by swimming around them in a molded way as shown in Figure 1. The mathematical model of each phase is explained below.Figure 1
Unique bubble-net feeding method and spiral moment for updating the position of humpback whales.


### 5.1.1. Exploration Phase

Humpback whales encircle prey for hunting. Equations (10) and (11) mathematically model this behavior as follows:(10)E⟶=C⟶⋅Z⟶∗t−Z⟶t,(11)Z⟶∗t+1=Z⟶∗t−A⟶⋅D,where t denotes the current iteration, Z∗ has provided the best solution so far, and D gives the location of humpback whales to prey at each step. A and C are coefficient vectors which are defined as follows:(12)A⟶=2a⟶⋅r⟶−a⟶,(13)C⟶=2⋅r⟶.where r∈0,1 is an arbitrary nominated vector. The value of a reduces from two to zero during exploration as well as in the exploitation phase. During helix-shaped movement, the distance between prey and the humpback whale is given as follows: (14)Z⟶t+1=D′⋅ebl⋅cos2πl+Z⟶∗t,where D′ represents the location of the ith whale to the prey and is defined as follows:(15)D′=Z⟶∗t−Z⟶t.Furthermore,b denotes the state of the logarithmic helix and l∈−1,1 is any arbitrary number. The shrinking surrounding technique of humpback whales during contracting loop is summarized as follows:(16)Z⟶t+1=Z⟶t−A⟶⋅Difp≤0.5,D⟶′⋅ebl⋅cos2πl+Z⟶∗tifp≥0.5.


### 5.1.2. Exploitation Phase

In the exploration phase (searching for prey), the heterogeneity of the vectorA is utilized. If A>1, the position of the search agent is updated, and the entire mechanism is modeled by the following equations:(17)D⟶=C⟶⋅Z⟶rand∗−Z⟶,(18)Z⟶t+1=Z⟶rand−A⟶⋅D⟶,where Zrand⟶ is an arbitrary position vector taken from the current population.Figure2 shows the flowchart of the WOA. It can be seen that the WOA creates a random initial population from candidate space and evaluates it using an error-based fitness function when the optimization process starts. After finding the best solution, the algorithm repeatedly executes the following steps until ending criterion is achieved.Figure 2
Graphical illustration of the working steps of the whale optimization algorithm and Nelder–Mead algorithm for the training of neurons in the LeNN architecture and the minimization of fitness functions.


## 5.1.1. Exploration Phase

Humpback whales encircle prey for hunting. Equations (10) and (11) mathematically model this behavior as follows:(10)E⟶=C⟶⋅Z⟶∗t−Z⟶t,(11)Z⟶∗t+1=Z⟶∗t−A⟶⋅D,where t denotes the current iteration, Z∗ has provided the best solution so far, and D gives the location of humpback whales to prey at each step. A and C are coefficient vectors which are defined as follows:(12)A⟶=2a⟶⋅r⟶−a⟶,(13)C⟶=2⋅r⟶.where r∈0,1 is an arbitrary nominated vector. The value of a reduces from two to zero during exploration as well as in the exploitation phase. During helix-shaped movement, the distance between prey and the humpback whale is given as follows: (14)Z⟶t+1=D′⋅ebl⋅cos2πl+Z⟶∗t,where D′ represents the location of the ith whale to the prey and is defined as follows:(15)D′=Z⟶∗t−Z⟶t.Furthermore,b denotes the state of the logarithmic helix and l∈−1,1 is any arbitrary number. The shrinking surrounding technique of humpback whales during contracting loop is summarized as follows:(16)Z⟶t+1=Z⟶t−A⟶⋅Difp≤0.5,D⟶′⋅ebl⋅cos2πl+Z⟶∗tifp≥0.5.


## 5.1.2. Exploitation Phase

In the exploration phase (searching for prey), the heterogeneity of the vectorA is utilized. If A>1, the position of the search agent is updated, and the entire mechanism is modeled by the following equations:(17)D⟶=C⟶⋅Z⟶rand∗−Z⟶,(18)Z⟶t+1=Z⟶rand−A⟶⋅D⟶,where Zrand⟶ is an arbitrary position vector taken from the current population.Figure2 shows the flowchart of the WOA. It can be seen that the WOA creates a random initial population from candidate space and evaluates it using an error-based fitness function when the optimization process starts. After finding the best solution, the algorithm repeatedly executes the following steps until ending criterion is achieved.Figure 2
Graphical illustration of the working steps of the whale optimization algorithm and Nelder–Mead algorithm for the training of neurons in the LeNN architecture and the minimization of fitness functions.


## 5.2. Nelder–Mead Algorithm

A Nelder–Mead (NM) algorithm is a direct search method also known as a downhill simplex method developed by Nelder and Mead in 1965 to solve different problems without any information about the gradient [41]. NM is a single path following a local search optimizer that can find good results if initialized with a better initial solution. A simplex consisting of n+1 vertices is set up to minimize a function f with dimensions n [42]. The NM algorithm generates a sequence of simplices by following four basic procedures, namely, reflection, expansion, contraction, and shrink. Further details about the NM algorithm can be found in [43]. Figure 2 shows the working procedure of the NM algorithm.In recent times, the most successful and effective trend in optimization is the action of integrating components from different methods. The foremost motivation behind the hybridization of diverse algorithmic ideas is to acquire better performing systems, which exploit and coalesce benefits of different techniques. Therefore, in this study, we have combined the global and local search optimization algorithms to achieve more efficient and robust solutions. The hybridized working procedure of the designed LeNN-WOA-NM algorithm is illustrated in Figure2.


## 6. Performance Measures

To examine the performance of the proposed algorithm, performance indices are defined in terms of mean absolute deviation (MAD), Theil’s inequality coefficient (TIC), and error in Nash–Sutcliffe efficiency (ENSE). Mathematical formulation forSapproxt, Xapproxt, and Yapproxt of the predator and prey model in the case of MAD, TIC, and ENSE is presented as follows:(19)MADS,MADX,MADY=1n∑i=1nSti−Sapproxti1n∑i=1nXti−Xapproxti1n∑i=1nYti−YapproxtiT,(20)TICS,TICX,TICY=1/n∑i=1nSti−Sapproxti21/n∑i=1nSti2+1/n∑i=1nSapproxti21/n∑i=1nXti−Xapproxti21/n∑i=1nXti2+1/n∑i=1nXapproxti21/n∑i=1nYti−Yapproxti21/n∑i=1nYti2+1/n∑i=1nYapproxti2T,(21)NSES,NSEX,NSEY=1−∑i=1nSti−Sapproxti2∑i=1nSti−S′ti2S′ti=1n∑i=1nSti1−∑i=1nXti−Xapproxti2∑i=1nXti−X′ti2X′ti=1n∑i=1nXti1−∑i=1nYti−Yapproxti2∑i=1nYti−Y′ti2Y′ti=1n∑i=1nYtiT,(22)ENSES,ENSEX,ENSEY=1−NSES,1−NSEX,1−NSEY,where n shows total input grid points. The values of performance measures of MAD, RMSE, and ENSE should be equal to zero for perfect modeling, while NSE should be equal to 1. The global versions of MAD, RMSE, and ENSE for the given mathematical model of the prey-predator system are formulated by given equations:(23)GMADS,GMADX,GMADY=1R∑r=1R1n∑i=1nSti−Sapproxti1R∑r=1R1n∑i=1nXti−Xapproxti1R∑r=1R1n∑i=1nYti−YapproxtiT,(24)GTICS,GTICX,GTICY=1R∑r=1R1/n∑i=1nSti−Sapproxti21/n∑i=1nSti2+1/n∑i=1nSapproxti21R∑r=1R1/n∑i=1nXti−Xapproxti21/n∑i=1nXti2+1/n∑i=1nXapproxti21R∑r=1R1/n∑i=1nYti−Yapproxti21/n∑i=1nYti2+1/n∑i=1nYapproxti2T,(25)GNSES,GNSEX,GNSEY=1R∑r=1R1−∑i=1nSti−Sapproxti2∑i=1nSti−S′ti2S′ti=1n∑i=1nSti1R∑r=1R1−∑i=1nXti−Xapproxti2∑i=1nXti−X′ti2X′ti=1n∑i=1nXti1R∑r=1R1−∑i=1nYti−Yapproxti2∑i=1nYti−Y′ti2Y′ti=1n∑i=1nYtiT,(26)GENSES,GENSEX,GENSEY=1−GNSES,1−GNSEX,1−GNSEY,where R denotes the number of independent runs. Global fitness (GFIT) is defined as the mean of fitness values attained in independent runs. Mathematically, GFIT is given as follows:(27)GFIT=1R∑r=1Rϵr.


## 7. Result and Discussion

The numerical solutions of the prey-predator system with immigrant prey under the influence of variations in various parameters are investigated with the proposed methodology. The exact solution for these nonlinear models is not known, so the comparative study is conducted with MATLAB solver ode45 and the homotopy perturbation method [44]. Six problems of the prey-predator model are considered for different cases depending on values of coefficients, i.e, α1,α2,β1,β2, γ1, and γ2. A graphical abstract of the paper is shown in Figure 3.Figure 3
Graphical overview of the paper.


### 7.1. Problem I: Influence of Variations inα1 on the Prey-Predator Model

In this problem, the effect of variations in the intrinsic growth rate of local preyα1 on population density is discussed. An error-based fitness function along with initial populations is given as follows:(28)Minimizeϵ=ϵ1+ϵ2+ϵ3+ϵ4+ϵ5+ϵ6,where ϵ1 to ϵ6 are defined by(29)ϵ1=1N∑n=134dSndt−Snα1−α1Sn50−0.2SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+0.1SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Five cases are considered, depending on the value ofα1.Case I:α1=0.06Case II:α1=0.08Case III:α1=0.10Case IV:α1=0.12Case V:α1=0.25The LeNN-WOA-NM algorithm is applied to prey-predator model equation (29) to study the influence of variations in the intrinsic growth rate of local prey. A comparison between approximate solutions obtained by the LeNN-WOA-NM algorithm for the prey-predator system with the homotopy perturbation method [44] for α1=0.12 is illustrated in Table 2. Approximate solutions for population densities of local prey, immigrant prey, and predator are given in Tables 3–5, respectively. Absolute errors are presented in Tables 6–8 and are graphically illustrated in Figure4. Table 9 represents the statistics of global values of performance indicators during 100 independent trails. Figure 5 shows the convergence of the fitness value. It can be seen that the fitness value for each case lies around 10−6 to 10−10. The percentage convergence of the fitness value and performance indicators during multiple runs is shown in Table 10. Trained neurons in the LeNN structure for obtaining best solutions are shown in Table 11. From Figure 6, the following conclusions are drawn:(i)
Population density of local prey has a direct relation with the intrinsic growth rate of local prey(ii)
Population density of immigrant prey has an inverse relation with the intrinsic growth rate of local prey(iii)
Population density of the predator varies directly with the intrinsic growth rate of local preyTable 2
Comparison between approximate solutions obtained by the LeNN-WOA-NM algorithm for the prey-predator system with the homotopy perturbation method [44] for α1=0.12.
tStXtYtHPOLeNN-WOA-NMHPOLeNN-WOA-NMHPOLeNN-WOA-NM0.00.5000000.500000000.5000000.500000000.5000000.500000000.20.5214580.521459560.4667620.466763210.5904400.591147650.40.5429420.542943660.4279690.427971170.6957770.699221870.60.5642440.564244220.3843290.384331040.8174680.826872840.80.5851550.585155970.3370080.337011390.9568440.977049171.00.6054780.605479600.2876060.287610081.1151271.15310441Table 3
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of local preySt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.00.500000290.500000240.500000010.500000000.500000000.20.515304610.517343940.519401660.521459560.535047430.40.530208040.534415820.538665660.542943660.571588100.60.544519360.551015320.557595650.564244220.609438840.80.558065660.566950840.575986240.585155970.648384251.00.570690560.582058780.593652210.605479600.68819372Table 4
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of immigrant preyXt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.00.499999790.499999970.499999990.500000000.500000000.20.466797890.466776600.466773630.466763210.466683200.40.428149410.428089980.428030110.427971170.427574260.60.384802960.384639710.384488970.384331040.383270750.80.337960500.337644890.337329770.337011390.334856361.00.289230600.288689110.288159860.287610080.28389875Table 5
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of the predatorYt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.00.499999910.499999190.499999990.500000000.500000000.20.590116940.590222920.590327560.590470070.591147650.40.694232910.694744100.695256280.695777350.699221870.60.813325170.814689840.816073830.817468430.826872840.80.948112170.950978270.953891980.956844910.977049171.01.099073441.104320821.109676021.115128371.15310441Table 6
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of local preySt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.04.25E − 109.71E − 114.27E − 119.01E − 113.05E − 130.23.87E − 103.87E − 111.80E − 098.37E − 119.38E − 140.45.99E − 111.16E − 101.85E − 093.57E − 113.30E − 130.61.46E − 101.33E − 102.21E − 093.05E − 111.05E − 130.81.37E − 102.15E − 113.80E − 097.04E − 111.82E − 121.09.93E − 112.11E − 106.82E − 116.49E − 114.90E − 14Table 7
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of immigrant preyXt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.03.33E − 111.68E − 099.24E − 123.85E − 123.05E − 130.25.66E − 127.98E − 114.55E − 103.72E − 129.38E − 140.44.40E − 114.29E − 104.03E − 101.09E − 113.30E − 130.66.22E − 112.53E − 095.15E − 103.52E − 121.05E − 130.82.46E − 123.20E − 091.24E − 091.07E − 121.82E − 121.02.09E − 111.17E − 105.55E − 114.44E − 124.90E − 14Table 8
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of the predatorYt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.02.17E − 124.68E − 115.32E − 101.16E − 103.05E − 130.23.04E − 111.21E − 103.32E − 102.52E − 119.38E − 140.47.91E − 117.16E − 125.91E − 101.83E − 103.30E − 130.67.12E − 111.43E − 105.80E − 121.36E − 101.05E − 130.82.17E − 182.46E − 122.41E − 103.16E − 121.82E − 121.01.60E − 104.16E − 121.40E − 107.66E − 114.90E − 14Figure 4
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations in the intrinsic growth rate of local prey on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Table 9
Statistics of global performance indices for variations in the intrinsic growth rate of local prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI1.21E – 072.47E − 051.30E − 053.82E – 051.21E − 072.26E – 051.61E − 053.21E − 061.55E − 073.70E − 051.30E − 052.00E − 06II5.15E – 081.64E − 058.52E − 068.13E – 065.25E − 081.69E – 051.19E − 051.15E − 065.66E − 082.67E − 059.42E − 063.87E − 07III8.12E – 072.74E − 051.39E − 052.61E – 051.45E − 072.14E – 051.51E − 052.90E − 061.61E − 073.89E − 051.40E − 051.15E − 06IV5.57E – 071.88E − 059.59E − 067.77E – 065.86E − 081.62E – 051.17E − 051.19E − 069.01E − 083.19E − 051.14E − 055.85E − 07V6.95E – 072.20E − 051.02E − 053.59E – 061.05E − 071.99E – 051.43E − 052.04E − 061.37E − 073.72E − 051.30E − 057.78E − 07Figure 5
Comparison between the box plots and the convergence graph of fitness evaluation over 100 independent runs for the prey-predator model with variations in the intrinsic growth rate of local prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inα1 (b) Convergence of the fitness value of the prey-predator model with variations in α1.
(a)(b)Table 10
Convergence analysis for variation inα1 on the prey-predator model.
CasesFITMADTICENSE≤10−7≤10−8≤10−9≤10−10≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−6≤10−7≤10−8StI98752551009935100996139223II803030100100491001005668368III7834201009441100996374333IV813740100994310010068834416V772750100974210010066876027XtI8137133964429854365236II85411101004311001005871203III744117197390100985264248IV86471641004641001005573379V804712198433100995768266YtI804116310095219957947147II864791100982810068546130III78368210091209957445100IV764293100942610060844106V702981100952110057741161Table 11
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations in the intrinsic growth rate of local preyα1.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.24298626−0.26686550.093285820.175246530.06571110.455184020.620505630.107119458−0.18453520.124780220.390269840.13404058−0.1205176−0.504835−0.1523680.29193615−0.062868140.04317448530.159463840.472050520.04884412−0.0142591−0.03468830.988866230.097152360.1815518110.86873273440.011206830.24610287−0.28379880.149467920.16097393−0.0744990.46391737−0.160730350.1133979795−0.0226690.09439050.09929740.341915940.582928480.009596280.06797262−0.112443580.0134443796−0.32161580.14101550.183250410.032286380.702878790.12845754−0.13618410.158063474−0.0583686170.487682090.02937720.17094612−0.00183540.148606320.651305640.676639640.3149841080.26112894980.083294060.04117897−0.0205859−0.1064367−0.1326281−0.0481534−0.0486244−0.333596420.56472758690.87211207−0.07438210.188828080.69221376−0.12135430.015845710.00730463−0.079596730.248981123100.746446580.015758840.170446530.131515920.136223070.26101806−0.1087744−0.09419355−32.2098458110.263135610.09923455−0.05061620.010091650.137956340.059686360.425062350.3027983990.094716386Case II10.33259121−0.36574490.210622760.259004420.159184370.410266290.29040796−1.70842597−0.2203902420.17889024−0.0431343−0.64111160.231101080.426257250.08496042−0.0683335−0.015372470.25253323230.072976370.1667513−0.0811545−0.2048892−0.1120065−0.0826385−0.6647359−0.16394699−0.009972364−0.13658610.47090897−0.06397773.43E − 050.1546024−0.1782350.46464940.1905217380.1954409665−0.0073580.131022910.24249896−0.1088503−0.16712490.16954851−0.7315874−0.03059029−0.0710649360.43362447−0.1556850.445635420.22079963−0.6459135−0.0878247−0.5471167−0.36068591−0.909664957−0.40725940.43108895−0.09730150.008594210.195191710.058313020.00294865−0.034408710.3723384358−0.015323−0.0127997−0.0797150.41948941−0.28784550.21735868−0.839503−0.048150650.02501122890.262173280.200284920.2222801−0.16303910.175203910.09747040.34891726−0.13212935−0.1419514100.092425250.19371262−0.15742670.341898510.284347720.084944610.46694404−0.1232987439.39674102110.13225195−0.0846055−0.02815680.142173720.318872310.429103770.267162070.0881731590.085561799Case III10.08792719−0.029981.21744620.691656860.06374932−0.04883220.114910280.090075909−0.3260432520.116100491.564065690.0795106−0.09029640.09097774−0.10179280.063716190.4922559960.40655735730.03262421−0.29831610.43612230.710230610.271307560.535740290.745466630.731579954−0.119253344−0.85342560.271025070.09970798−0.46308680.195882230.25655668−0.17334020.4924538−0.027388455−0.43735730.068925560.111137510.224249730.18678342−0.00955920.138609860.342147207−0.1281684660.85585046−0.00865720.708042980.165324530.463875910.19167055−0.41558450.3581745130.51958643670.4731327−0.01310840.538180720.00107541−0.23717140.215388680.270670760.007722887−0.116512018−0.1730569−0.05426561.004630150.333047160.421115070.54031591−0.0982948−0.309648890.22526495790.019505510.3195719−0.17817250.00222254−0.05279060.338791620.145361080.1850604020.02254770310−0.27583160.27171296−0.04286510.075659860.29894632−0.11823050.196849460.6264427160.062105445110.28574702−0.26285120.061022570.044035020.224146150.403643070.552061570.429077344−0.00256072Case IV10.25369644−0.4777306−0.8361847−0.06212640.27351167−0.60650990.68147388−1.43275825−0.1093238620.06991497−0.1139906−0.4735192−0.5854725−0.16611070.07770566−0.24603360.0524627850.1910779173−0.2351744−0.4870632−0.0716903−0.25560449.83E − 050.02156199−0.3140772−0.4208503−0.131827984−0.04898050.03121692−0.2705836−0.2176759−0.3388043−0.1570093−0.29472920.0008612380.09113584950.120660690.00528409−0.18766530.18129324−0.2476741−0.08439310.332233530.0378772150.5524868386−0.29176940.09885985−0.655587−0.58219360.02484774−0.2678441−0.12098250.019527337−0.204123270.21801060.10156035−0.0839076−0.04081370.202271790.08783601−1.14728740.088467804−0.10068238−0.3594382−0.0329465−0.0085748−0.0762450.071221−0.2677479−0.69126710.072673015−0.6012941490.12647778−0.0321146−0.1925977−0.1982251−0.1306954−0.4290151−0.4035842−0.01597822−0.0938865710−0.7218529−0.127281−0.253891−0.520109−0.148345−0.13018490.45904036−0.145556930.21087826111−0.0415557−0.0626802−0.0135924−0.2376142−0.26023340.31718881−0.0300309−0.077803080.119273654Case V10.115351690.093314110.064595980.137089550.025231310.144401870.536785940.2896713310.26549624421.20187280.66566882−0.01957540.553142020.209201250.535357010.30925340.735712480.32098752430.20319881−0.82568851.095402670.096363490.12623634−0.13614610.29987708−0.14730456−0.185248540.02729770.50122086−0.2264851−0.04105390.34541103−0.15725140.886125890.2192311830.29604182650.134366630.091248160.315029470.03670816−0.28421280.280475380.304938810.228728298−0.8384621860.590905350.22192090.358369330.132875320.072794670.57461590.052420130.220708110.05973650270.571477680.174590270.749635190.965947050.18063006−0.25056210.11584087−0.021910610.16070292580.140835490.51514220.061702260.181894870.194128650.221329980.9645296−0.155782240.4221313669−0.0987012−0.34473240.305476970.227528650.566319040.084389660.0392108−0.152703490.520225963100.053379540.389787040.324107670.040672550.406039520.21526779−0.0360576−0.00016971−1.80718751110.168828680.422474010.06696436−0.05341270.386903630.40861018−0.10325480.0398748980.21583529Figure 6
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations in the intrinsic growth rate of local prey on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


### 7.2. Problem II: Effect of Variations inβ1 on the Prey-Predator Model

In this problem, the effect of variations in the positive impact of force of interactionβ1 between local and immigrant prey on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(30)Minimizeϵ=ϵ1+ϵ2+ϵ3+ϵ4+ϵ5+ϵ6,where ϵ1 to ϵ6 are defined as follows:(31)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−β1SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+0.1SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Five cases are considered, depending on the value ofβ1.Case I:β1=0.05Case II:β1=0.10Case III:β1=0.15Case IV:β1=0.20Case V:β1=0.25Approximate solutions for the effect of variations inβ1 on population densities of the prey-predator model are given in Tables 12–14. Absolute errors in our solution for population densities are presented in Tables 15–17, respectively. The solution of the design scheme overlaps the exact solution with absolute errors that lie between 10−9 and 10−11 as illustrated in Figure 7. Figure 8 shows the behavior of fitness evaluation for the prey-predator model under the influence of β1. Convergence analysis of performance measures is given in Tables 18 and 19. Trained neurons obtained by the LeNN-WOA-NM algorithm for different cases of Eq (31) are shown in Table 20. From the solutions (see Figure 9), the following conclusions can be drawn:(i)
The negative force of interaction between local and immigrant prey has a direct relation with population density of local prey and the predator, while population density is inversely related toβ1Table 12
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of local preySt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.00.499998110.500005650.499999940.500000010.499999940.20.513941920.516445760.518937940.521454400.523991640.40.527961600.532937410.537904110.542936440.548018430.60.542034800.549368830.556749850.564242870.571838420.80.556109940.565646410.575311940.585152930.595162041.00.570094520.581677040.593450330.605473300.61773135Table 13
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of immigrant preyXt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.00.499999920.500005140.499999880.500000010.499999690.20.466800420.466793640.466776260.466762600.466734400.40.428179650.428125470.428041930.427970080.427895210.60.384888980.384706930.384519030.384329220.384132610.80.338106610.337755970.337374960.337012850.336631061.00.289435010.288834540.288230650.287607010.28697657Table 14
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of the predatorYt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.00.500000320.499992150.500000050.500000030.499999730.20.590041120.590162710.590306470.590440660.590577410.40.693908170.694515600.695164390.695780330.696409360.60.812601400.814195980.815842610.817467430.819123650.80.946901170.950165240.953485380.956842390.960258531.01.097427541.103224251.109124571.115130131.12125476Table 15
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of local preySt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.01.38E − 081.77E − 091.10E − 091.24E − 091.04E − 100.27.06E − 091.51E − 098.29E − 109.85E − 103.47E − 100.44.21E − 093.33E − 112.05E − 102.81E − 101.92E − 100.64.01E − 092.43E − 102.77E − 102.79E − 102.07E − 090.81.44E − 082.19E − 106.01E − 104.68E − 101.57E − 091.02.13E − 084.83E − 101.02E − 093.60E − 107.30E − 10Table 16
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of immigrant preyXt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.01.39E − 093.44E − 113.07E − 101.07E − 112.30E − 090.22.29E − 094.35E − 095.04E − 105.02E − 124.24E − 100.47.01E − 105.98E − 106.77E − 109.67E − 111.96E − 100.61.41E − 112.14E − 104.69E − 106.77E − 111.29E − 090.81.61E − 092.21E − 108.72E − 127.39E − 111.42E − 091.01.01E − 092.03E − 091.13E − 093.02E − 108.42E − 11Table 17
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of the predatorYt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.03.95E − 113.48E − 094.48E − 091.83E − 104.83E − 100.28.22E − 101.51E − 093.01E − 093.19E − 103.25E − 100.43.03E − 111.96E − 101.20E − 092.60E − 115.34E − 100.69.03E − 103.60E − 102.64E − 095.37E − 108.30E − 120.81.36E − 091.83E − 097.57E − 105.75E − 102.25E − 101.03.84E − 097.37E − 091.40E − 099.52E − 102.47E − 10Figure 7
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inβ1 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Figure 8
Comparison between the box plots and convergence graphs of fitness evaluation over 100 independent runs for the prey-predator model with variations in force of interaction between local and immigrant prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inβ1. (b) Convergence of the fitness value of the prey-predator model with variations in β1.
(a)(b)Table 18
Comparison through global performance indices for variations inβ1 on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI8.39E − 067.18E − 053.75E − 054.26E − 042.61E − 062.74E − 042.14E − 042.14E − 041.63E − 061.21E – 044.34E – 051.12E − 05II1.00E − 066.17E − 053.16E − 051.30E − 048.45E − 076.17E − 054.37E − 051.99E − 051.54E − 061.23E – 044.33E – 051.09E − 05III1.64E − 067.26E − 053.67E − 051.93E − 042.27E − 067.43E − 055.27E − 054.08E − 051.49E − 061.41E – 044.91E – 052.10E − 05IV1.44E − 066.46E − 053.23E − 051.17E − 041.73E − 066.86E − 054.88E − 053.59E − 051.91E − 061.23E – 044.33E – 051.55E − 05V5.98E – 075.97E − 052.96E − 057.46E − 051.07E − 066.29E − 054.45E − 052.23E − 051.63E − 061.04E – 043.72E – 059.00E − 06Table 19
Convergence analysis of population density of local, immigrant prey, and predator under the influence of variations inβ1.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI895516099836100913571384II87612561008410100942973397III988027210079121009425783712IV9877253100858100932286439V908440410081191009441836021XtI977536119977151008521946627II9981481410083111009017966927III98843312100827100909966519IV998439111008010100921471173V7727501008281008916936132YtI97753611996231009223998041II94783681006841009024977938III81374099655100872080499IV9774347996741009021814311V772750100974210010066876027Table 20
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations inβ1.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.430069150.19531180.585262780.20397331.09461748−0.30423830.451373750.6134726210.213003912−0.18170230.12967851−0.2238481−0.1299324−0.2205129−0.24340370.893839230.3475711850.1926919053−0.1639676−0.19069710.01654308−0.26832750.21586942−0.06508260.15826136−0.386370640.54484618540.6432707−0.2243964−0.0644709−0.14154830.00139286−0.12214450.21055624−0.284600750.3252831885−0.0633810.20479542−0.08377010.37097529−0.2551801−0.14574440.85822173−0.358983850.0958098676−0.05872150.09571579−0.21513340.238571310.306612510.34850165−0.0057450.6202316490.5993485570.31197992−0.0344083−0.05428880.01997741−0.1641653−0.05154490.000740050.000287462−0.3738123280.85365924−0.01697890.255722330.80492801−0.0023113−0.2852425−0.0522238−0.22752380.1463248649−0.2837938−0.069528−0.22218920.316929490.34731804−0.0409553−0.168709−0.125701740.02922959310−0.2720053−0.00832730.713032790.01306660.043472410.104905791.16491711−0.146201030.240782136110.02324570.044139090.57180799−0.39617790.13805184−0.02508270.1189892−0.099767430.676863876Case II10.200842570.018481050.011511120.188341850.077143851.001992210.330907110.5919545220.17208180120.777315120.20280783−0.1534665−0.0898678−0.02221640.050976590.39013928−0.16168208−0.874142883−0.89326870.01562423−0.8790331−0.8785872−0.24009780.011269390.62600687−0.074690290.2932000124−0.12540860.08769516−0.4580886−0.2608019−0.08707280.233163880.49916445−0.27253148−0.0084165550.05410689−0.0916165−0.4295043−0.5942252−0.03405130.03471072−0.0634013−0.11877784−0.7721202160.225648440.11645954−0.01944170.06716894−0.0067884−0.2770185−0.3074844−0.17361157−0.0294093470.257781660.028023310.328260150.0874406−0.92391470.57760720.12893563−0.074945180.1013149118−0.61522520.08467719−0.005507−0.00168850.269596760.000273920.01072972−0.059671890.7940112819−0.49165910.034230080.352213850.16735460.200979460.03970559−0.53290310.056857150.093143062100.093720710.45237283−0.42700230.122709540.01633451−0.02449120.32114085−0.085295710.896172613111.84E − 05−0.1783898−0.0967566−0.032303−0.3339006−0.05318620.052990490.002508360.0227459Case III10.328659270.73232241−0.04262260.36804480.12592558−0.17118130.68950930.3656759770.181511720.063523260.016153480.719430010.027656140.57709610.125083080.725793830.2342071480.27897832330.059020580.74414566−0.0387953−0.12858020.02597920.02360138−0.03907720.3343788710.37388573940.016729860.30780827−0.03051870.764788350.00922025−0.00031760.39416458−0.141058660.23126747550.168869060.16663601−0.15577560.14372478−0.00535840.398679360.538671370.083576797−0.09300856−0.1209806−0.37495880.216202920.046822170.16582226−0.18138350.31094709−0.118440160.0038172927−0.07560490.268685170.112320290.215547350.081722680.333028420.258253170.129740332−0.2709219380.17875086−0.01300430.7761683−0.2683637−0.05682850.08575422−0.34051820.078523738−0.103839490.521793020.034586650.064754270.069562310.250723130.027159480.259072610.038549325−0.02300298100.447951010.146474340.16811506−0.00085980.316891280.217170580.559240270.030536512−0.63332453110.07886187−0.12418510.139439470.01341213−0.62677380.34912288−0.1931707−0.012986310.579458203Case IV10.369260070.34021884−0.09357310.161692130.092468161.385538050.440585440.4999267590.12746079620.324814770.210042150.173602920.06611299−0.1524440.172302550.48717520.240390062−0.132472693−0.15700690.04398714−0.13721160.103939270.277123880.219293050.359345590.0724365430.59968018440.10042725−0.0405233−0.02761240.00768115−0.11527970.466710010.326976350.2188138030.1157854250.29815990.078954080.05929431−0.25757910.182107180.279497530.11208642−0.13550479−0.1930335560.017742780.07944380.343013220.3079287−0.02884290.352851730.08051491−0.03832890.30152455170.15448802−0.02900920.379211770.405773980.023177340.197600180.54536650.1901281440.13497225980.073032090.08335033−0.0272924−0.02684070.114950610.175658620.145480130.0504595460.99955956690.066380160.093651920.460467010.311571490.247886270.397456990.14745623−0.02038990.11010541100.145744110.270190120.08964107−0.0136530.69998870.147064260.143116930.0855874730.406653134110.00422723−0.02023260.097542644.84E − 050.011302540.45443510.20276643−0.021765530.041899718Case V10.416681550.551414170.317056470.53356928−2.0883662−0.26107040.24747335−0.495683480.1511660120.281111040.298731990.124158610.010351830.3912644−0.22508560.114638350.5725941980.6395297123−0.28800510.476324660.021482480.428317540.25753792−0.10679690.49570912−0.565745240.14021161640.283821540.074922320.25762031−0.00516160.449995−0.25836760.008623420.042826998−0.245326385−0.18289930.056100350.411578990.10823296−0.00221370.438408260.145317870.5040825650.0733579816−0.35365450.40455640.05410188−0.1039612−0.44871770.185451370.002075410.071838234−0.1819850970.23523025−0.11973780.227370420.27683837−0.45025420.53741643−0.19293150.0945797620.3349364458−0.01019840.607967930.06730913−0.07407840.15137890.63740825−0.3265244−0.272198550.2351746669−0.0038892−0.32395580.104053890.411949020.01303423−0.0585570.255794720.192432911−0.11924924100.209832440.246830530.091350390.04756230.44987648−0.0005977−0.29770780.0131149460.599953402110.176991850.19078854−0.4635923−0.0023057−0.20808410.138832930.21279433−0.199952010.280170884Figure 9
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations inβ1 on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


### 7.3. Problem III: Effect of Variations inα2 on the Prey-Predator Model

In this problem, the effect of variations in the increasing rate of immigrant preyα2 on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(32)Minimizeϵ=ϵ1+ϵ2+ϵ3+ϵ4+ϵ5+ϵ6,where ϵ1 to ϵ6 are defined as follows:(33)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−0.2SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xnα2−α2Xnk2+0.1SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Four cases are considered, depending on the value ofα2.Case I:α2=0.05Case II:α2=0.10Case III:α2=0.15Case IV:α2=0.25The LeNN-WOA-NM algorithm is used to optimize the population densities of equation (33). Table 21 represents the comparison between the Ranga–Kutta method and the proposed technique LeNN-WOA-NM. A comparison between population density of local prey St, immigrant prey Xt, and predator Yt with variations in α2 is shown in Table 22. Statistics of absolute errors are given in Table 23 and graphically presented in Figure 10. The behavior of the fitness function for different cases is shown on boxplots as demonstrated in Figure 11. The mean values of the fitness function lie around 10−7 as shown in Tables 24 and 25. Unknown parameters achieved by the LeNN-WOA-NM algorithm for the given problem are given in Table 26. From solutions (see Figure 12), the following conclusion can be drawn:(i)
Population density of local prey, immigrant prey, and predator varies directly with an increasing rate of immigrate preyα2Table 21
Comparison between solutions obtained by the LeNN-WOA-NM algorithm and ode45 for the prey-predator model withα2=0.15.
tStode45Xtode45Ytode450.00.500000600.50.499997130.50.500000850.50.20.521400430.52148360.462187510.469148750.590212060.590555110.40.542738690.543043030.419677190.432428940.694768210.696305340.60.563822410.564468040.373304840.390424020.814973650.818807100.80.584439230.585541390.324332830.344202130.952014790.959486471.00.604429480.606055270.274401980.295309561.107021861.11963769Table 22
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population densities of local, immigrant prey, and predator under the influence of the increasing rate of immigrant preyα2.
tStXtYtα2=0.05α2=0.10α2=0.15α2=0.25α2=0.05α2=0.10α2=0.15α2=0.25α2=0.05α2=0.10α2=0.15α2=0.250.00.499999960.499999700.500000600.500000180.499999970.500000230.499997130.499998700.4999999970.5000005620.5000008470.5000010330.20.521312910.521356890.521400430.521501120.453212050.457682610.462187510.471379240.5897703670.5899978480.5902120550.5906658560.40.542357970.542554470.542738690.543124360.403568430.411546230.419677190.436442150.6928196750.6937941060.6947682060.6967944650.60.562996940.563416230.563822410.564666160.352169170.362585310.373304840.395688610.8101597120.8125464260.8149736470.8200217540.80.583081610.583756220.584439230.585876620.300375870.312131320.324332830.350168880.9427785700.9473401200.9520147920.9618234251.00.602454640.603425500.604429480.606536820.249687630.261761580.274401980.301442451.0917048231.0992366981.1070218621.123532428Table 23
Absolute errors obtained by the LeNN-WOA-NM algorithm for population densities of local, immigrant prey, and predator under the influence ofα2.
tStXtYtα2=0.05α2=0.10α2=0.15α2=0.25α2=0.05α2=0.10α2=0.15α2=0.25α2=0.05α2=0.10α2=0.15α2=0.2501.38E − 081.54E − 092.74E − 091.77E – 091.18E − 091.88E − 104.68E − 096.17E − 093.17E − 106.73E − 118.97E – 102.19E − 100.26.09E − 091.40E − 093.92E − 093.57E – 096.24E − 101.22E − 106.51E − 109.47E − 091.15E − 114.00E − 119.21E – 107.22E − 110.44.30E − 091.63E − 101.22E − 093.29E – 101.59E − 098.95E − 115.03E − 091.78E − 117.33E − 101.29E − 101.48E – 104.81E − 100.62.00E − 096.47E − 102.54E − 101.89E – 095.33E − 101.49E − 112.95E − 098.36E − 101.06E − 099.40E − 125.36E – 105.34E − 100.88.76E − 097.02E − 101.16E − 093.77E – 099.98E − 119.28E − 117.06E − 103.10E − 111.31E − 102.41E − 112.99E – 104.31E − 1011.40E − 089.57E − 101.35E − 096.00E – 093.71E − 105.62E − 125.19E − 091.12E − 117.18E − 109.95E − 124.08E – 103.35E − 10Figure 10
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inα2 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Figure 11
Comparison between the box plots and convergence graphs of fitness evaluation over 100 independent runs for the prey-predator model with variations in the increasing rate of immigrant prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inα2. (b) Convergence of the fitness value of the prey-predator model with variations in α2.
(a)(b)Table 24
Statistics of global performance indices for variations in the intrinsic growth rate of immigrant prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI2.25E − 065.70E − 052.85E − 057.97E – 059.91E − 076.19E − 054.75E − 057.97E − 056.03E − 078.42E − 053.00E – 055.97E − 06II1.21E − 064.64E − 052.36E − 054.09E – 055.08E − 074.68E − 053.50E − 058.22E − 064.05E − 076.93E − 052.49E – 052.93E − 06III1.24E − 064.33E − 052.17E − 053.13E − 054.31E − 074.76E − 053.46E − 058.99E − 065.75E − 077.72E − 052.82E – 053.20E − 06IV5.22E − 074.55E − 052.29E − 054.13E − 053.65E − 074.91E − 053.39E − 051.15E − 058.39E − 078.64E − 053.06E – 054.58E − 06Table 25
Convergence analysis of population density of local prey, immigrant prey, and predator under the influence of variations in the increasing rate of immigrant preyα2.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI9462201100899100962989469II1007016110091141009826914214III996516110096111009927964711IV1008738610085201009842875321XtI9878451410085131009017956537II100895417100871110092191008638III999252111008981009717997932IV10089461610087121009621987127YtI998639121007651009527978972II10092461110079610097311009360III10083461510074810099261008756IV987843111006941009819999354Table 26
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations in the increasing rate of immigrant preyα2.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.29809866−0.03153660.580127350.472975030.168707710.254228270.044875360.4942418270.21904284220.07902509−0.142623−0.53265720.299332070.363340650.400516670.79930460.015001430.04649882930.37716273−0.0698453−0.5278093−0.59197780.06003798−0.6780052−0.6088520.6091109920.11634367240.064186860.242090730.119104830.29508708−0.00842750.0309784−0.0262412−0.34143893−0.536647465−0.0499254−0.02504030.051066860.063262910.384623390.077014170.121516760.1190981770.099603766−0.03959580.01489951−0.28834840.084806050.18102453−0.02914670.00187112−0.01144597−0.3412298770.156847550.086090640.207353430.518959580.40337040.210525920.498563420.2626537380.0345593378−0.03401640.048636210.222766390.044039820.213446940.318286370.242865030.0437745890.5295648479−0.25192170.0619910.136546150.18808650.144864010.376159850.468323640.1746203680.266799487100.54249116−0.03241190.225653230.107307140.217643420.300076760.077317380.2110310890.199472911110.129040850.04855899−0.46668450.25677349−0.0420107−0.1514589−0.17789850.104252082−0.3576282Case II10.22014881−0.09553890.019227330.440992690.169020980.210581670.566622760.214534099−0.125478432−0.5105560.025220920.565965080.2853307−0.01404990.1381917−0.00046460.4918297680.3823676313−0.16576350.01472461−0.11479830.393519120.21069293−0.26858750.395098870.072673362−0.1828391940.052823270.07711052−0.20514380.129927560.193529140.24990450.436458580.3821135730.05871400150.644992530.257626410.238082030.8495944−0.0043216−0.00171050.254571620.313158881−0.1355356160.09858966−0.02731450.487879870.093747910.21889407−0.0062504−0.13447480.044741355−0.214952670.23380371−0.12602150.32895940.05799592−0.3333640.15201786−0.4456943−0.069769270.00213826380.16035874−0.1136321−0.00365010.01897557−0.06503620.090919151.14914925−0.035859990.16591012790.016831660.030273560.013660120.167758170.200576770.14595760.538281720.1303285280.267951969100.267719360.122046520.2996044−0.11528120.188191760.17111444−0.38691720.095372589−1.0466345511−0.14301490.0042418−0.1108715−0.14027110.17676553−0.2193704−0.27470910.125441060.108988227Case III10.208334770.87977164−0.3298640.374971710.158035121.086854990.81716797−0.82852980.2076197242−0.009672−0.0834121−0.00629890.005063330.45907436−0.2830568−0.055188−1.296406970.3831862863−0.12604430.5466368−0.2547923−0.1978145−0.30147920.014824530.16026646−0.12624417−0.142116884−0.0459971−0.18508240.440100740.85262434−0.2293558−0.18033120.15934035−0.278759070.92589765350.61609669−0.0654918−0.22499821.06003476−0.06056080.0377206−0.24654680.095037929−0.003680466−0.42218770.19157710.29040894−0.1447329−0.1716272−0.28141120.102520450.1405726430.10066670870.12304762−0.0412171−0.2947437−0.33442590.15047463−0.0948965−0.3952389−0.186674990.2203322980.22471377−0.1369576−0.01359290.461183750.096751020.254991790.6880506−0.04146236−0.020095289−0.00285250.28384483−0.01379880.868556160.151195960.43141743−0.3108490.1770470290.639644264100.10012064−0.0365188−0.34958450.37524460.150820620.397374230.02610775−0.116321040.264941005110.6577301−0.17291860.28691512−0.277895−0.2620078−0.31836860.46666040.182729564−0.40213292Case IV10.275842730.144401320.05159123−0.00145010.029449330.438578970.415416940.3641066070.13763746721.03785820.22634893−0.13364490.374794840.323901380.116831911.106341910.7438702810.13372108130.305552340.141783720.098708630.1378461−0.01342310.252356050.158475340.5467479610.26755784840.49396717−0.00778540.128639560.3091184−0.04341350.984756940.07222518−0.00152740.0328421550.964611720.487610510.02679450.44318330.195001840.288898580.10389701−0.060481020.0127791256−0.0311005−0.24588980.013237040.392157780.075258550.0342729−0.02390090.3050476820.2408197467−0.00097850.278034190.434940140.08135220.165273610.10128603−0.02357030.1312367830.34259740880.032459480.055916940.9180906−0.0641160.192863010.46309077−0.02267510.2656931280.19181902290.04135314−0.0595634−0.018826−0.04538630.104730060.091887330.121393520.2816876520.039757676100.062064640.108719550.214806490.12965376−0.02754390.211693440.08952384−0.23209362−0.11227638110.41183634−0.01232830.31935826−0.06298370.007029690.084862440.225785890.1943565610.207691077Figure 12
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations in the increasing rate of immigrant prey on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


### 7.4. Problem IV

In this problem, the effect of variations in the negative force of interaction between local prey and immigrant prey on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(34)Minimizeϵ=ϵ1+ϵ2+ϵ3+ϵ4+ϵ5+ϵ6,where ϵ1 to ϵ6 are defined as follows:(35)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−0.2SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+β2SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Three cases are considered, depending on the value ofβ2.Case I:β2=0.05Case II:β2=0.10Case III:β2=0.20The LeNN-WOA-NM algorithm is used to optimize the population densities of prey-predator model equation (35). Population densities of local prey, immigrant prey, and the predator are given in Tables 27-28, while absolute errors are shown in Table 29. The absolute errors in the solutions of the proposed technique lie between 10−9 and 10−12 as shown in Figure 13. Table 30 shows global values of performance indices and convergence of the fitness value as shown in Figure 14. The statistics shown in Table 31 illustrate that the global values of different performance functions lie between 10−5 and 10−7, which highlights the robustness of the technique. The values of the weights in the LeNN structure for recreation of the approximate solutions are given in Table 32. From Figure 15, the following conclusion can be drawn:(i)
Population density of local prey varies directly with variations inβ2(ii)
Population density of immigrant prey and predator varies inversely with variations inβ2Table 27
Comparison between solutions obtained by the LeNN-WOA-NM algorithm and ode45 for the prey-predator model under the influence of variations inβ2.
tStode45Xtode45Ytode450.00.500002600.50.499996630.50.500004340.50.20.521405240.521408030.462026650.462027280.590204950.590209760.40.542730860.542742160.419188970.419193220.694726290.694733920.60.563798210.563803330.372432670.372431380.814839800.814839440.80.584394820.584398120.323085530.323078370.951682950.951689481.00.604344660.604353850.272810910.272808341.106366291.10638491Table 28
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence of the negative impact of the force of interaction on immigrant preyβ2.
tStXtYtβ2=0.05β2=0.10β2=0.20β2=0.05β2=0.10β2=0.20β2=0.05β2=0.10β2=0.200.00.499999970.499998390.500002600.499999980.500000380.499996630.499999980.499999710.5000043440.20.521474820.521466980.521405240.469148100.466766810.462026650.590550200.590446210.5902049530.40.543034740.542947310.542730860.432426210.427972540.419188970.696300310.695789680.6947262890.60.564467510.564244610.563798210.390424950.384325450.372432670.818802140.817473630.8148398000.80.585536730.585163100.584394820.344204810.337019100.323085530.959478280.956852740.9516829451.00.606046820.605480650.604344660.295308300.287602420.272810911.119627321.115148381.106366286Table 29
Absolute errors obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence of the negative impact of the force of interaction on immigrant preyβ2.
tStXtYtβ2=0.05β2=0.10β2=0.20β2=0.05β2=0.10β2=0.20β2=0.05β2=0.10β2=0.200.03.09E − 094.62E − 091.04E − 083.11E − 102.55E − 106.59E − 095.14E − 102.00E − 092.51E − 090.22.26E − 092.04E − 095.78E − 093.65E − 102.21E − 097.44E − 101.62E − 102.46E − 093.36E − 090.42.02E − 093.08E − 092.81E − 114.80E − 153.55E − 091.32E − 114.25E − 108.45E − 102.38E − 090.62.16E − 105.66E − 101.57E − 095.83E − 101.08E − 091.97E − 096.71E − 111.23E − 095.01E − 100.81.96E − 091.72E − 108.46E − 103.17E − 101.01E − 091.45E − 113.56E − 112.75E − 096.92E − 091.01.65E − 091.04E − 107.87E − 105.40E − 101.98E − 091.91E − 101.55E − 114.27E − 094.61E − 09Figure 13
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inβ2 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Table 30
Statistics of global performance indices for variations in the increasing rate of immigrant prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI5.58E − 074.70E − 052.36E − 053.69E − 054.72E − 075.09E − 053.59E − 051.10E − 058.00E − 079.48E − 053.40E – 054.71E − 06II1.48E − 065.02E − 052.49E − 055.05E − 054.76E − 075.11E − 053.63E − 051.24E − 057.09E − 079.42E − 053.35E – 055.21E − 06III2.13E − 065.09E − 052.58E − 056.78E − 058.11E − 075.28E − 053.85E − 051.18E − 057.99E − 077.69E − 052.79E – 053.93E − 06Figure 14
Comparison between the box plots and convergence graphs of fitness evaluation over 100 independent runs for the prey-predator model with variations inβ2. (a) Box plots for fitness evaluation of the prey-predator model with variations in β2. (b) Convergence of the fitness value of the prey-predator model with variations in β2.
(a)(b)Table 31
Convergence analysis of population density of local, immigrant prey, and predator under the influence of variations inβ2.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI10086323100921010010025924410II1006312010089101009926894810III966514210089151009737895215XtI10086451410088910094151006825II1009148121009071009419977330III9982421510088121009620986927YtI998039510067210094241008648II10076391410071310093231008842III9886481210077310096261008964Table 32
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations inβ2.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.275351440.31269649−0.2126230.355379880.472481030.108226330.521149510.185720875−0.0089027820.615168380.24930343−0.1133566−0.26988980.230757451.456621310.002854710.54326608−0.039481033−0.33317120.136830020.07415358−0.007558−0.14145220.00429373−0.2363309−0.53264189−0.126252194−0.1724685−0.1569281−0.1780917−0.02301890.06567477−0.1619849−0.074340.2804944630.72430884150.342352550.166015010.337796120.382142220.172063980.375675440.25394552−0.204242560.49782990260.00647555−0.06710070.418484980.4655630.178523990.259100010.62807088−0.25015769−0.0360779470.39917241−0.20864230.27916279−0.1982513−0.10032140.086122760.23112207−0.11265536−0.0463813180.37891578−0.14710920.255162720.47437651−0.04797260.39541230.233714660.1557852320.1115255199−0.15698220.21029864−0.280122−0.04986250.078988720.18498889−0.1404023−0.289746740.27339322810−0.1405838−0.0119767−0.16351250.11249313−0.06545650.69188379−0.0877791−0.1594545−0.77121428110.392969390.02457412−0.0200821−0.12014470.171618890.036300630.048113920.279716864−0.22688875Case II10.501052520.825363390.580454830.45518044−0.1160988−0.30930880.147352550.5765523590.19212727120.261418060.217962380.31815591−0.16766240.4234933−0.16140250.009758380.539177880.7758476183−0.1877622−0.1568166−0.38590220.34192011−0.11621830.787573110.40990792−0.291269520.09465586940.02599785−0.17645390.174915140.031499050.44911473−0.1695519−0.0583248−0.00447879−0.062893735−0.2923416−0.10657390.387648390.10732316−0.0035220.430096060.202684350.6275844490.023338846−0.38680090.44970016−0.0422409−0.2581683−0.41051990.016568080.000667930.075964394−0.1868564470.05739953−0.1193984−0.26071230.17181358−0.27136280.66311474−0.124996−0.079813080.5059381328−0.02163280.80670030.31701813−0.15115880.140927880.65770582−0.1918301−0.004073960.1210127269−8.30E − 07−0.0664668−0.11686640.63180611−0.1478398−0.15926770.270268430.1131427720.003053122100.09718671−0.0187463−0.0760148−0.04440810.299805990.00951493−0.1584802−0.011772520.821398692110.193967170.18738185−0.4407164−0.1400116−0.15198690.038500450.32211631−0.106169950.295182432Case III10.368776140.41870330.005517210.371063780.76553227−0.11350280.367662840.6565438770.07401287820.571231580.390292350.269251950.02897603−0.0177996−0.27249760.68341531−0.566598910.0352719053−0.2786739−0.1350866−0.40798240.18172611−0.1092217−0.65352410.19817369−0.181044410.12771873540.248960520.133223240.342458990.185354870.183366550.32712799−0.14222260.1528342880.02504536250.291759670.21119283−0.19578160.31132875−0.0532529−0.33059640.04517261−0.73533086−0.2105123960.386092490.275908470.58209501−0.535581−0.18238940.31125028−0.00330380.0826663750.32301675570.01164667−0.14217680.243283170.264800370.021735020.078484850.617754820.2318933740.07434687880.593813210.06375022−0.36897130.35998444−0.03686660.329072190.148773580.1434832470.33440679−0.2877027−0.08129760.585135370.114541340.04186864−0.38143560.28795352−0.00041389−0.1333593110−0.30739940.0111124−0.00688850.04196657−0.1275650.14629639−0.40819820.00302860.33443128110.06421859−0.09197220.34881249−0.0684058−0.14121790.4979470.198741510.0711107530.298703303Figure 15
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations inβ2 on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


### 7.5. Problem V: Effect of Variation inγ1 on the Prey-Predator Model

In this problem, the effect of variations in the catching rate of local preyγ1 on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(36)Minimizeε=ε1+ε2+ε3+ε4+ε5+ε6,where ϵ1 to ϵ6 are defined as follows:(37)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−0.2SnXn+γ1SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+0.1SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Four cases are considered, depending on variations inγ1.Case I:γ1=0.005Case II:γ1=0.015Case III:γ1=0.020Case IV:γ1=0.025The LeNN-WOA-NM algorithm is used to optimize the population densities of equation (37). Table 33 represents the comparison between the Ranga–Kutta method (ode45) and the proposed technique LeNN-WOA-NM algorithm. Approximate solutions for population density of local prey St, immigrant prey Xt, and the predator Yt with different values of γ1 are presented in Table 34. Statistics of absolute errors are given in Table 35 and graphically shown in Figure 16. The minimum and mean values of the fitness function for each case lie between 10−9 and 10−7, respectively as illustrated in Figure 17. Convergence analysis of the fitness value, MAD, TIC, and ENSE is given in Tables 36-37. Unknown neurons of LeNN are shown in Table 38. From Figure 18, the following conclusion can be drawn:(i)
Population density of local prey and population density of the predator vary inversely with variations inγ1(ii)
Population density of immigrant prey has a direct relation with variations inγ1Table 33
Comparison between solutions obtained by the LeNN-WOA-NM algorithm and ode45 for the prey-predator model withγ1=0.015.
tStode45Xtode45Ytode450.00.500000290.50.499999710.50.500000020.50.20.521179670.521174760.46675620.466764390.590433250.590425180.40.542302620.542298910.427974970.427979250.695720320.6957040.60.56315320.563151080.384347740.384354420.817264190.817258110.80.583510510.583506110.337055510.337060360.956375380.956370131.00.603152380.60315130.287692460.287699271.114209391.11419232Table 34
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence of variations inγ1.
tStXtYtγ1=0.005γ1=0.015γ1=0.020γ1=0.025γ1=0.005γ1=0.015γ1=0.020γ1=0.025γ1=0.005γ1=0.015γ1=0.020γ1=0.0250.00.499999820.500000290.500000000.499999970.499999820.499999710.499999980.500000170.5000005040.5000000180.4999999970.5000001310.20.521747000.521179670.520892690.520606610.466751740.466756200.466760070.466766150.5904495820.5904332540.5904113170.5903940500.40.543594720.542302620.541659740.541016310.427954660.427974970.427976180.427997140.6958466350.6957203220.6956306740.6955587970.60.565339290.563153200.562061770.560975590.384308290.384347740.384378340.384397480.8176786560.8172641910.8170462120.8168373420.80.586808480.583510510.581862690.580224420.336956440.337055510.337112360.337160850.9573164870.9563753780.9558944400.9554196061.00.607822510.603152380.600836190.598526840.287517770.287692460.287779660.287874201.1160620941.1142093921.1132594931.112332477Table 35
Absolute errors obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence of variations inγ1.
tStXtYtγ1=0.005γ1=0.015γ1=0.020γ1=0.025γ1=0.005γ1=0.015γ1=0.020γ1=0.025γ1=0.005γ1=0.015γ1=0.020γ1=0.0250.02.74E − 098.47E − 103.60E − 107.80E − 115.74E − 092.00E − 094.47E − 094.19E − 128.54E − 106.23E − 099.79E − 125.10E − 130.22.76E − 094.69E − 102.92E − 101.04E − 113.59E − 098.16E − 114.78E − 095.77E − 105.53E − 106.28E − 091.38E − 118.72E − 110.45.70E − 104.68E − 104.62E − 133.63E − 102.51E − 093.02E − 104.26E − 102.89E − 105.49E − 108.44E − 101.36E − 122.06E − 100.61.89E − 091.24E − 101.54E − 102.49E − 124.04E − 113.74E − 104.26E − 092.55E − 111.47E − 123.54E − 099.17E − 111.14E − 090.82.36E − 091.80E − 112.01E − 128.31E − 101.28E − 094.93E − 103.67E − 095.03E − 103.59E − 103.49E − 099.80E − 114.43E − 101.02.82E − 091.85E − 104.54E − 121.17E − 092.58E − 101.22E − 154.68E − 091.08E − 102.82E − 104.48E − 093.19E − 101.19E − 09Figure 16
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inγ1 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Figure 17
Comparison between the box plots and the convergence graph of fitness evaluation over 100 independent runs for the prey-predator model with variations in the catching rate of local prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inγ1. (b) Convergence of the fitness value of the prey-predator model with variations in γ1.
(a)(b)Table 36
Statistics of global performance indices for variations in the catching rate of local prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI2.50E − 065.21E − 052.69E − 052.50E − 069.20E − 074.63E − 053.33E − 051.51E − 059.19E − 078.88E − 053.10E − 058.98E − 06II7.32E − 066.20E − 053.14E − 055.63E − 052.43E − 066.86E − 054.92E − 057.58E − 051.44E − 061.06E − 043.80E − 051.70E − 05III2.86E − 064.87E − 052.48E − 059.34E − 057.60E − 075.39E − 053.89E − 051.86E − 051.39E − 068.53E − 053.09E − 055.88E − 06IV1.78E − 064.62E − 052.38E − 057.52E − 055.18E − 074.32E − 053.10E − 051.46E − 056.42E − 078.01E − 052.83E − 059.08E − 06Table 37
Convergence analysis of population density of local prey, immigrant prey, and predator under the influence of variations inγ1.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI957020510085211009437855621II94672529981231009443835721III957117210090191009642915519IV977829610091261009742905726XtI979059221009115100942797838II97865624998416999128967840III9985641710088151008921948133IV9993572210093201009631988146YtI99875115997771009525999159II96784815997381009330988762III9882491510075101009730998861IV9890601710093371009337999364Table 38
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations inγ1.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.10433411−0.9093816−0.05513840.10814164−0.0975817−0.40247590.110579010.030277215−0.076830932−0.1015821−0.23390630.11183609−0.00796490.12409457−0.00663660.031204470.0659841510.1698154793−0.17322190.05666544−0.3418627−0.15762190.205116340.23158339−0.3370419−0.154182380.19415431340.28706660.04587599−0.0750170.0317597−0.1567915−0.29606070.126614180.176984003−0.0830713250.18918051−0.26348620.092969590.177127140.066902930.055961240.06180474−0.14442377−0.9894472560.08011991−0.02418430.043841940.09323427−0.10295410.30404503−0.1508149−0.050515030.1910929767−0.4734878−0.18889670.139719480.249251060.002252150.00501846−0.9616715−0.197021870.01134063180.19753354−0.0151743−0.0358307−0.2646482−0.4318811−0.2034545−0.0870373−0.088720140.07936334790.073523210.043742590.17641082−0.004594−0.0709633−0.1885710.012284110.051586673−0.1707669100.05966701−0.0094660.11306597−1.0101745−0.11955560.30440396−0.165510.0727288860.05362364110.244192780.215236270.125015060.108878510.12112026−0.1899001−0.157371−0.159061210.128530508Case II10.425624360.224051350.025039560.35262391−0.22536640.262594550.271116860.3287966910.22892787220.416669880.5078151−0.18890930.121713910.209981050.38048580.303907470.3891165320.46070179630.105776950.060028460.21768786−0.53933980.390275590.139022260.473504190.4654198640.66197000740.465174070.09698879−0.16334780.110142970.197399070.66549310.05604469−0.003926970.27877739450.601000290.019320570.568304980.23543140.414040570.20437334−0.1267241−0.131288940.29384493960.081277360.106281460.386109420.235399620.171875290.385790850.352251370.0676494280.18848928170.354984110.044902120.049006640.064765870.352380740.502358060.362536390.203702979−0.079204880.384970770.018235750.34815044−0.07555940.494782790.380354950.16102938−0.124828140.3405400119−0.05617650.73291315−0.18902150.0053794−0.44388960.466410960.2889986−0.017214940.10349715710−8.42E − 05−0.01914780.330575130.011386710.210127450.280139280.417905120.0110224340.111373073110.03758884−0.00909620.298419490.061814350.282112330.10443310.76029441−0.053057410.400559476Case III10.34210927−0.1686215−0.06371520.234271170.252859750.392147380.269000940.5655202880.4852402582−0.0413096−0.00813010.106145370.60106519−0.07658960.421885660.341462510.5755837060.08123241630.29904262−0.075813−0.19675630.459391190.15935586−0.0737513−0.03714280.3957801950.52602507540.31966978−0.02968460.134928190.744499180.138230990.147330570.312965730.1703563410.44903041650.385111710.027429150.192458680.174809460.260677330.19256129−0.19667440.0266527150.64855673760.15405406−0.0855950.039002980.20316026−0.09393630.369996350.141163030.183493540.15728401470.011602360.19042424−0.0099165−0.1955586−0.01452490.51154810.092342830.0501830530.5295766948−0.0119424−0.00172220.19894821−0.04590190.169954560.099121220.442771620.2681979320.3604521989−0.0135090.028267380.13366696−0.09861210.125556210.29387726−0.03515110.15868406−0.02656018100.36233153−0.05044810.464221880.331541310.044883490.444237770.182334570.0026998190.25017744110.207808520.061875820.190128340.47309423−0.00014550.285446120.326557470.2322044930.038376682Case IV10.08743197−0.02745860.138384540.10216597−0.6205394−1.16254730.08890874−0.468214890.03190664120.85138918−0.25960890.27697631−0.0395679−0.31904840.47807164−1.07668990.1945224680.05147395530.32484395−0.03201280.049102590.59721285−0.7402654−0.5159275−0.54570820.5845145650.33420834241.265395290.0637656−0.32985420.007467950.12675579−0.21199380.06894660.0780600630.34124362250.27124475−0.0145706−0.8795409−0.2258574−1.04891440.18138843−0.7172783−0.446545880.2883523536−0.0615221−0.03571810.05167451−0.0093528−0.90763030.005299290.111646850.0841619050.0213049777−0.93304190.005870160.08238699−0.00384620.5318456−0.0029243−0.56249830.07126723−0.071945938−0.22228510.159644890.159293470.02677628−0.08595810.71257834−0.15642590.011627909−0.011415649−0.07467320.03595603−0.16882950.918538690.0432707−0.19964060.9708232−0.207186870.491048316100.111059850.03339826−0.0412419−0.9111904−0.07572990.383182960.140461090.0682333990.44790321911−0.9956577−0.03257510.958620420.01858614−0.03282860.2342197−0.8047726−0.035479630.018056695Figure 18
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations inγ1 on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


### 7.6. Problem VI: Effect of Variations inγ2 on the Prey-Predator Model

In this problem, the effect of variations in the catching rate of immigrant prey on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(38)Minimizeε=ε1+ε2+ε3+ε4+ε5+ε6.where ϵ1 to ϵ6 are defined as follows:(39)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−0.2SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+0.1SnXn+γ2XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Five cases are considered, depending on the value ofγ2.Case I:γ2=0.7Case II:γ2=0.8Case III:γ2=1.0Case IV:γ2=1.1The LeNN-WOA-NM algorithm is applied to prey-predator model equation (39) to study the influence of variations in the catching rate of immigrant prey. The results calculated by the designed algorithm are compared with those of the Runge–Kutta method using ode45 in MATLAB as shown in Table 39. Approximate solutions for population densities of local prey, immigrant prey, and the predator are given in Table 40. Absolute errors are presented in Table 41 and graphically shown in Figure 19. The minimum fitness values in Figure 20 reflect the accuracy of the solutions by the proposed algorithm. The convergence of our numerical approach is assessed by fitness values and statistical analysis of the performance indicators (MAD, TIC, and ENSE), see Table 42. Table 43 shows the correctness of the proposed algorithm to tackle real-world problems. Trained neurons in LeNN are shown in Table 44. From Figure 21, the following conclusions can be drawn:(i)
Population densities of local prey, immigrant prey, and the predator vary inversely with variation in the catching rate of immigrant preyγ2.Table 39
Comparison between solutions obtained by the LeNN-WOA-NM algorithm and ode45 for the prey-predator model withγ2=0.8.
tStode45Xtode45Ytode450.00.500000060.50.50000050.50.499999820.50.20.521510970.521174760.471866320.466764390.590685680.590425180.40.543161820.542298910.4381850.427979250.696940280.6957040.60.564756030.563151080.399355260.384354420.820536990.817258110.80.586083320.583506110.356171160.337060360.963191650.956370131.00.606933670.60315130.309865090.287699271.126496341.11419232Table 40
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence ofγ2.
tStXtYtγ2=0.7γ2=0.8γ2=1.0γ2=1.1γ2=0.7γ2=0.8γ2=1.0γ2=1.1γ2=0.7γ2=0.8γ2=1.0γ2=1.10.00.500000380.500000060.499999990.499999840.500000480.50000050.500000040.499999990.4999997790.4999998180.4999999840.4999996970.20.52156930.521510970.521410350.521350990.477021170.471866320.461725790.456737880.5909334430.5906856820.5902027520.5899718050.40.54339380.543161820.5427270.542506280.448656410.4381850.4180060.408282360.698110510.6969402820.6946498760.6935324810.60.565282850.564756030.563745210.563255510.415003050.399355260.36990490.356058780.8236969380.8205369860.8144934660.8115992330.80.587048620.586083320.5842640.583396270.376511730.356171160.318968250.301952150.9698202460.9631916450.9507666550.9449383661.00.608471360.606933670.604096540.602772190.33402430.309865090.267081790.248146241.1385322641.1264963371.1044111611.094260679Table 41
Absolute errors obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence ofγ2.
tStXtYtγ2=0.7γ2=0.8γ2=1.0γ2=1.1γ2=0.7γ2=0.8γ2=1.0γ2=1.1γ2=0.7γ2=0.8γ2=1.0γ2=1.10.01.26E − 091.32E − 104.79E − 105.02E − 101.48E − 111.87E − 107.79E − 105.53E − 109.67E − 101.58E − 095.89E − 101.90E − 090.21.31E − 094.54E − 111.43E − 105.38E − 101.04E − 104.83E − 128.26E − 102.09E − 103.34E − 111.58E − 096.20E − 102.77E − 100.47.22E − 103.13E − 124.65E − 107.47E − 117.86E − 103.53E − 102.65E − 097.30E − 101.16E − 093.29E − 107.51E − 111.87E − 090.69.27E − 102.41E − 121.09E − 101.51E − 103.37E − 126.92E − 101.12E − 092.57E − 101.29E − 091.11E − 095.35E − 101.29E − 090.82.37E − 095.90E − 112.34E − 112.37E − 109.11E − 101.53E − 101.69E − 103.50E − 113.98E − 101.50E − 095.10E − 101.19E − 101.02.91E − 091.44E − 101.20E − 122.88E − 109.82E − 101.61E − 101.58E − 094.63E − 115.13E − 101.68E − 097.98E − 104.24E − 10Figure 19
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inγ2 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Figure 20
Comparison between the box plots and convergence graphs of fitness evaluation over 100 independent runs for the prey-predator model with variations in the catching rate of immigrant prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inγ2. (b) Convergence of the fitness value of the prey-predator model with variations in γ2.
(a)(b)Table 42
Convergence analysis of population density of local prey, immigrant prey, and predator under the influence of variations inγ2.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI956524510088201009541895920II986424410092181009945936018III866312410090161009532905516IV986720410092211009744926221XtI9985522210085201009230946932II10087572210087231009127987333III9883461310082101009316987128IV9990541910094171009723998244YtI98844417100771110092281008757II998147189975131009731998755III998441151007561009426998956IV9985501410072910094341008755Table 43
Statistics of global performance indices for variations in the catching rate of local prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI1.98E − 064.80E − 052.43E − 055.78E − 055.99E − 075.04E − 053.39E − 052.68E − 059.28E − 079.13E − 053.24E − 055.44E − 06II2.17E − 063.83E − 051.95E − 052.91E − 055.21E − 074.70E − 053.30E − 051.44E − 051.36E − 069.17E − 053.25E − 051.01E − 05III2.16E − 064.95E − 052.52E − 055.55E − 051.01E − 065.58E − 054.16E − 051.53E − 057.10E − 078.79E − 053.17E − 056.81E − 06IV1.55E − 064.01E − 052.00E − 054.85E − 055.75E − 074.49E − 053.40E − 057.86E − 067.28E − 078.59E − 053.08E − 055.50E − 06Table 44
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations inγ1.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.27328231−0.0407082−0.00365680.190954050.196199130.253340580.154480320.3391994560.64626032920.002036930.007612690.246550490.67156041−0.12021850.410298230.354566580.167129547−0.0222010830.22854565−0.1965182−0.07997060.4234230.03898304−0.1160494−0.19322110.3272984570.58963410940.43068766−0.0047135−0.00677490.76144510.244368120.103860.45804558−0.105041050.09663247950.42381358−0.15026520.518274420.107579960.246969080.11302049−0.04908280.0544971930.54146435560.030643450.387909−0.09083380.20021796−0.10933340.498037490.085022350.2923121810.1689517997−0.02562720.181349440.04903339−0.14391810.016780890.015222870.122286450.0406107450.5042415958−0.00292340.053691650.161186290.092383710.35203171−0.07424580.582963550.3263946220.31123246990.01433551−0.00628260.114617840.000702810.081706550.27442616−0.06249750.046728012−0.00253833100.52132655−0.07458240.323731330.316258890.064695840.355187160.432659940.001675421−0.07508196110.04700892−0.14311380.190960340.610866830.058795620.071452470.316027320.2945685460.005492958Case II10.228880980.156627270.721183980.663963010.48400554−0.00149870.2930431−0.36246252−0.1982755620.320391690.213497160.26632024−0.1848203−0.0673308−0.06482460.131252990.0416570280.0299543743−0.08293880.030948320.003035220.62951247−0.19254380.24736655−0.19631210.4315314350.7080957674−0.1658907−0.1447463−0.6955688−0.232498−0.26160950.270252270.22782911−0.28937041−0.129944965−0.01747580.232151820.042900160.12334482−0.2128439−0.0360696−0.14592740.1433832820.1276542126−0.0712296−0.3626918−0.03486680.26004380.131684030.451310950.73463526−0.091098940.1000586577−0.0244551−0.0112202−0.1246274−0.34457320.31169855−0.06807980.148097380.0295262120.1760398998−0.29639480.276843650.02435043−0.1570692−0.10580240.01048308−0.02870220.0522047370.5002902529−0.0250215−0.10242610.056514850.61926963−0.2055780.43411127−0.37484230.146360873−0.4197865510−0.28473110.13336075−0.0903463−0.1464593−0.3554768−0.04750720.16657978−0.04184543−0.0441026211−0.00129140.02947063−0.2027891−0.01192550.472381010.094487350.42979806−0.05005385−0.00793862Case III10.17136911−0.6220119−0.06594510.160013380.06604026−0.62594620.417434420.347451460.12340272620.132480930.224071980.16846536−0.01740040.460482060.161721960.63635670.0278926630.3492171153−0.0296660.0085174−0.4187633−0.05734020.296738520.13113638−0.7934564−0.181700450.08725763940.509163990.04634821−0.14370190.08923057−0.2567021−0.29347460.26333479−0.21308607−0.0250535150.28783823−0.26621630.07660280.099144850.576658110.043000830.20662283−0.48985711−1.0974074360.08917337−0.10500090.238182560.02035607−0.10721920.2659477−0.0029764−0.131483110.3136653667−0.7503685−0.17044710.18719634−0.0233035−0.00889380.26343414−1.0737247−0.190186990.28010052480.342300890.00558946−0.1427919−0.2957063−0.64904810.105464590.273225470.041606539−0.1027815790.258797360.219539270.110529190.00206003−0.0680416−0.21899140.03786660.055450179−0.26693807100.114828690.042566210.24740165−1.02102880.108268520.10811710.289721620.107963625−0.01416421110.216458910.26405307−0.2682646−0.02512170.2673504−0.2407415−0.3201356−0.126449070.363928738Case IV10.163210441.001017390.446108510.19117524−0.18691310.032844490.487094060.4928390050.06881143120.050340140.274485020.055954270.34759862−0.0428559−0.46566320.412112910.5503501620.26125445730.023098730.00700133−0.3300406−0.032857−0.2059139−0.15721320.430692780.1197089380.1320454840.60643185−0.16770850.013515050.012693450.07414495−0.2533687−0.19908810.074189084−0.0087421950.06424168−0.03942680.040623410.103460710.24294981−0.14069910.15305910.007966010.1423245696−0.1607398−0.0067892−0.1142057−0.3542569−0.3179260.896468750.19029585−0.07772530.1382154397−0.0249375−0.1421790.211814490.00745531−0.0807929−0.0606121−0.3138917−0.04179489−0.0839599580.012718540.07904033−0.20324210.30678256−0.0970233−0.0973761−0.2681465−0.26023333−0.0574878590.157291660.020040370.084417750.50899351−0.01071120.16547740.05321702−0.55038421−0.01142195100.455268140.308482240.03885589−0.0670780.03709569−0.11034190.00132508−0.0601008−0.00809993110.00144682−0.0160469−0.14856550.07678165−0.0117808−0.039538−0.15537240.090185082−0.05537005Figure 21
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations inγ2 on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


## 7.1. Problem I: Influence of Variations inα1 on the Prey-Predator Model

In this problem, the effect of variations in the intrinsic growth rate of local preyα1 on population density is discussed. An error-based fitness function along with initial populations is given as follows:(28)Minimizeϵ=ϵ1+ϵ2+ϵ3+ϵ4+ϵ5+ϵ6,where ϵ1 to ϵ6 are defined by(29)ϵ1=1N∑n=134dSndt−Snα1−α1Sn50−0.2SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+0.1SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Five cases are considered, depending on the value ofα1.Case I:α1=0.06Case II:α1=0.08Case III:α1=0.10Case IV:α1=0.12Case V:α1=0.25The LeNN-WOA-NM algorithm is applied to prey-predator model equation (29) to study the influence of variations in the intrinsic growth rate of local prey. A comparison between approximate solutions obtained by the LeNN-WOA-NM algorithm for the prey-predator system with the homotopy perturbation method [44] for α1=0.12 is illustrated in Table 2. Approximate solutions for population densities of local prey, immigrant prey, and predator are given in Tables 3–5, respectively. Absolute errors are presented in Tables 6–8 and are graphically illustrated in Figure4. Table 9 represents the statistics of global values of performance indicators during 100 independent trails. Figure 5 shows the convergence of the fitness value. It can be seen that the fitness value for each case lies around 10−6 to 10−10. The percentage convergence of the fitness value and performance indicators during multiple runs is shown in Table 10. Trained neurons in the LeNN structure for obtaining best solutions are shown in Table 11. From Figure 6, the following conclusions are drawn:(i)
Population density of local prey has a direct relation with the intrinsic growth rate of local prey(ii)
Population density of immigrant prey has an inverse relation with the intrinsic growth rate of local prey(iii)
Population density of the predator varies directly with the intrinsic growth rate of local preyTable 2
Comparison between approximate solutions obtained by the LeNN-WOA-NM algorithm for the prey-predator system with the homotopy perturbation method [44] for α1=0.12.
tStXtYtHPOLeNN-WOA-NMHPOLeNN-WOA-NMHPOLeNN-WOA-NM0.00.5000000.500000000.5000000.500000000.5000000.500000000.20.5214580.521459560.4667620.466763210.5904400.591147650.40.5429420.542943660.4279690.427971170.6957770.699221870.60.5642440.564244220.3843290.384331040.8174680.826872840.80.5851550.585155970.3370080.337011390.9568440.977049171.00.6054780.605479600.2876060.287610081.1151271.15310441Table 3
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of local preySt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.00.500000290.500000240.500000010.500000000.500000000.20.515304610.517343940.519401660.521459560.535047430.40.530208040.534415820.538665660.542943660.571588100.60.544519360.551015320.557595650.564244220.609438840.80.558065660.566950840.575986240.585155970.648384251.00.570690560.582058780.593652210.605479600.68819372Table 4
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of immigrant preyXt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.00.499999790.499999970.499999990.500000000.500000000.20.466797890.466776600.466773630.466763210.466683200.40.428149410.428089980.428030110.427971170.427574260.60.384802960.384639710.384488970.384331040.383270750.80.337960500.337644890.337329770.337011390.334856361.00.289230600.288689110.288159860.287610080.28389875Table 5
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of the predatorYt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.00.499999910.499999190.499999990.500000000.500000000.20.590116940.590222920.590327560.590470070.591147650.40.694232910.694744100.695256280.695777350.699221870.60.813325170.814689840.816073830.817468430.826872840.80.948112170.950978270.953891980.956844910.977049171.01.099073441.104320821.109676021.115128371.15310441Table 6
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of local preySt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.04.25E − 109.71E − 114.27E − 119.01E − 113.05E − 130.23.87E − 103.87E − 111.80E − 098.37E − 119.38E − 140.45.99E − 111.16E − 101.85E − 093.57E − 113.30E − 130.61.46E − 101.33E − 102.21E − 093.05E − 111.05E − 130.81.37E − 102.15E − 113.80E − 097.04E − 111.82E − 121.09.93E − 112.11E − 106.82E − 116.49E − 114.90E − 14Table 7
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of immigrant preyXt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.03.33E − 111.68E − 099.24E − 123.85E − 123.05E − 130.25.66E − 127.98E − 114.55E − 103.72E − 129.38E − 140.44.40E − 114.29E − 104.03E − 101.09E − 113.30E − 130.66.22E − 112.53E − 095.15E − 103.52E − 121.05E − 130.82.46E − 123.20E − 091.24E − 091.07E − 121.82E − 121.02.09E − 111.17E − 105.55E − 114.44E − 124.90E − 14Table 8
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of the predatorYt under the influence of the intrinsic growth rate of local prey α1.
tα1=0.06α1=0.08α1=0.10α1=0.12α1=0.250.02.17E − 124.68E − 115.32E − 101.16E − 103.05E − 130.23.04E − 111.21E − 103.32E − 102.52E − 119.38E − 140.47.91E − 117.16E − 125.91E − 101.83E − 103.30E − 130.67.12E − 111.43E − 105.80E − 121.36E − 101.05E − 130.82.17E − 182.46E − 122.41E − 103.16E − 121.82E − 121.01.60E − 104.16E − 121.40E − 107.66E − 114.90E − 14Figure 4
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations in the intrinsic growth rate of local prey on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Table 9
Statistics of global performance indices for variations in the intrinsic growth rate of local prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI1.21E – 072.47E − 051.30E − 053.82E – 051.21E − 072.26E – 051.61E − 053.21E − 061.55E − 073.70E − 051.30E − 052.00E − 06II5.15E – 081.64E − 058.52E − 068.13E – 065.25E − 081.69E – 051.19E − 051.15E − 065.66E − 082.67E − 059.42E − 063.87E − 07III8.12E – 072.74E − 051.39E − 052.61E – 051.45E − 072.14E – 051.51E − 052.90E − 061.61E − 073.89E − 051.40E − 051.15E − 06IV5.57E – 071.88E − 059.59E − 067.77E – 065.86E − 081.62E – 051.17E − 051.19E − 069.01E − 083.19E − 051.14E − 055.85E − 07V6.95E – 072.20E − 051.02E − 053.59E – 061.05E − 071.99E – 051.43E − 052.04E − 061.37E − 073.72E − 051.30E − 057.78E − 07Figure 5
Comparison between the box plots and the convergence graph of fitness evaluation over 100 independent runs for the prey-predator model with variations in the intrinsic growth rate of local prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inα1 (b) Convergence of the fitness value of the prey-predator model with variations in α1.
(a)(b)Table 10
Convergence analysis for variation inα1 on the prey-predator model.
CasesFITMADTICENSE≤10−7≤10−8≤10−9≤10−10≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−6≤10−7≤10−8StI98752551009935100996139223II803030100100491001005668368III7834201009441100996374333IV813740100994310010068834416V772750100974210010066876027XtI8137133964429854365236II85411101004311001005871203III744117197390100985264248IV86471641004641001005573379V804712198433100995768266YtI804116310095219957947147II864791100982810068546130III78368210091209957445100IV764293100942610060844106V702981100952110057741161Table 11
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations in the intrinsic growth rate of local preyα1.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.24298626−0.26686550.093285820.175246530.06571110.455184020.620505630.107119458−0.18453520.124780220.390269840.13404058−0.1205176−0.504835−0.1523680.29193615−0.062868140.04317448530.159463840.472050520.04884412−0.0142591−0.03468830.988866230.097152360.1815518110.86873273440.011206830.24610287−0.28379880.149467920.16097393−0.0744990.46391737−0.160730350.1133979795−0.0226690.09439050.09929740.341915940.582928480.009596280.06797262−0.112443580.0134443796−0.32161580.14101550.183250410.032286380.702878790.12845754−0.13618410.158063474−0.0583686170.487682090.02937720.17094612−0.00183540.148606320.651305640.676639640.3149841080.26112894980.083294060.04117897−0.0205859−0.1064367−0.1326281−0.0481534−0.0486244−0.333596420.56472758690.87211207−0.07438210.188828080.69221376−0.12135430.015845710.00730463−0.079596730.248981123100.746446580.015758840.170446530.131515920.136223070.26101806−0.1087744−0.09419355−32.2098458110.263135610.09923455−0.05061620.010091650.137956340.059686360.425062350.3027983990.094716386Case II10.33259121−0.36574490.210622760.259004420.159184370.410266290.29040796−1.70842597−0.2203902420.17889024−0.0431343−0.64111160.231101080.426257250.08496042−0.0683335−0.015372470.25253323230.072976370.1667513−0.0811545−0.2048892−0.1120065−0.0826385−0.6647359−0.16394699−0.009972364−0.13658610.47090897−0.06397773.43E − 050.1546024−0.1782350.46464940.1905217380.1954409665−0.0073580.131022910.24249896−0.1088503−0.16712490.16954851−0.7315874−0.03059029−0.0710649360.43362447−0.1556850.445635420.22079963−0.6459135−0.0878247−0.5471167−0.36068591−0.909664957−0.40725940.43108895−0.09730150.008594210.195191710.058313020.00294865−0.034408710.3723384358−0.015323−0.0127997−0.0797150.41948941−0.28784550.21735868−0.839503−0.048150650.02501122890.262173280.200284920.2222801−0.16303910.175203910.09747040.34891726−0.13212935−0.1419514100.092425250.19371262−0.15742670.341898510.284347720.084944610.46694404−0.1232987439.39674102110.13225195−0.0846055−0.02815680.142173720.318872310.429103770.267162070.0881731590.085561799Case III10.08792719−0.029981.21744620.691656860.06374932−0.04883220.114910280.090075909−0.3260432520.116100491.564065690.0795106−0.09029640.09097774−0.10179280.063716190.4922559960.40655735730.03262421−0.29831610.43612230.710230610.271307560.535740290.745466630.731579954−0.119253344−0.85342560.271025070.09970798−0.46308680.195882230.25655668−0.17334020.4924538−0.027388455−0.43735730.068925560.111137510.224249730.18678342−0.00955920.138609860.342147207−0.1281684660.85585046−0.00865720.708042980.165324530.463875910.19167055−0.41558450.3581745130.51958643670.4731327−0.01310840.538180720.00107541−0.23717140.215388680.270670760.007722887−0.116512018−0.1730569−0.05426561.004630150.333047160.421115070.54031591−0.0982948−0.309648890.22526495790.019505510.3195719−0.17817250.00222254−0.05279060.338791620.145361080.1850604020.02254770310−0.27583160.27171296−0.04286510.075659860.29894632−0.11823050.196849460.6264427160.062105445110.28574702−0.26285120.061022570.044035020.224146150.403643070.552061570.429077344−0.00256072Case IV10.25369644−0.4777306−0.8361847−0.06212640.27351167−0.60650990.68147388−1.43275825−0.1093238620.06991497−0.1139906−0.4735192−0.5854725−0.16611070.07770566−0.24603360.0524627850.1910779173−0.2351744−0.4870632−0.0716903−0.25560449.83E − 050.02156199−0.3140772−0.4208503−0.131827984−0.04898050.03121692−0.2705836−0.2176759−0.3388043−0.1570093−0.29472920.0008612380.09113584950.120660690.00528409−0.18766530.18129324−0.2476741−0.08439310.332233530.0378772150.5524868386−0.29176940.09885985−0.655587−0.58219360.02484774−0.2678441−0.12098250.019527337−0.204123270.21801060.10156035−0.0839076−0.04081370.202271790.08783601−1.14728740.088467804−0.10068238−0.3594382−0.0329465−0.0085748−0.0762450.071221−0.2677479−0.69126710.072673015−0.6012941490.12647778−0.0321146−0.1925977−0.1982251−0.1306954−0.4290151−0.4035842−0.01597822−0.0938865710−0.7218529−0.127281−0.253891−0.520109−0.148345−0.13018490.45904036−0.145556930.21087826111−0.0415557−0.0626802−0.0135924−0.2376142−0.26023340.31718881−0.0300309−0.077803080.119273654Case V10.115351690.093314110.064595980.137089550.025231310.144401870.536785940.2896713310.26549624421.20187280.66566882−0.01957540.553142020.209201250.535357010.30925340.735712480.32098752430.20319881−0.82568851.095402670.096363490.12623634−0.13614610.29987708−0.14730456−0.185248540.02729770.50122086−0.2264851−0.04105390.34541103−0.15725140.886125890.2192311830.29604182650.134366630.091248160.315029470.03670816−0.28421280.280475380.304938810.228728298−0.8384621860.590905350.22192090.358369330.132875320.072794670.57461590.052420130.220708110.05973650270.571477680.174590270.749635190.965947050.18063006−0.25056210.11584087−0.021910610.16070292580.140835490.51514220.061702260.181894870.194128650.221329980.9645296−0.155782240.4221313669−0.0987012−0.34473240.305476970.227528650.566319040.084389660.0392108−0.152703490.520225963100.053379540.389787040.324107670.040672550.406039520.21526779−0.0360576−0.00016971−1.80718751110.168828680.422474010.06696436−0.05341270.386903630.40861018−0.10325480.0398748980.21583529Figure 6
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations in the intrinsic growth rate of local prey on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


## 7.2. Problem II: Effect of Variations inβ1 on the Prey-Predator Model

In this problem, the effect of variations in the positive impact of force of interactionβ1 between local and immigrant prey on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(30)Minimizeϵ=ϵ1+ϵ2+ϵ3+ϵ4+ϵ5+ϵ6,where ϵ1 to ϵ6 are defined as follows:(31)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−β1SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+0.1SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Five cases are considered, depending on the value ofβ1.Case I:β1=0.05Case II:β1=0.10Case III:β1=0.15Case IV:β1=0.20Case V:β1=0.25Approximate solutions for the effect of variations inβ1 on population densities of the prey-predator model are given in Tables 12–14. Absolute errors in our solution for population densities are presented in Tables 15–17, respectively. The solution of the design scheme overlaps the exact solution with absolute errors that lie between 10−9 and 10−11 as illustrated in Figure 7. Figure 8 shows the behavior of fitness evaluation for the prey-predator model under the influence of β1. Convergence analysis of performance measures is given in Tables 18 and 19. Trained neurons obtained by the LeNN-WOA-NM algorithm for different cases of Eq (31) are shown in Table 20. From the solutions (see Figure 9), the following conclusions can be drawn:(i)
The negative force of interaction between local and immigrant prey has a direct relation with population density of local prey and the predator, while population density is inversely related toβ1Table 12
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of local preySt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.00.499998110.500005650.499999940.500000010.499999940.20.513941920.516445760.518937940.521454400.523991640.40.527961600.532937410.537904110.542936440.548018430.60.542034800.549368830.556749850.564242870.571838420.80.556109940.565646410.575311940.585152930.595162041.00.570094520.581677040.593450330.605473300.61773135Table 13
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of immigrant preyXt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.00.499999920.500005140.499999880.500000010.499999690.20.466800420.466793640.466776260.466762600.466734400.40.428179650.428125470.428041930.427970080.427895210.60.384888980.384706930.384519030.384329220.384132610.80.338106610.337755970.337374960.337012850.336631061.00.289435010.288834540.288230650.287607010.28697657Table 14
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population density of the predatorYt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.00.500000320.499992150.500000050.500000030.499999730.20.590041120.590162710.590306470.590440660.590577410.40.693908170.694515600.695164390.695780330.696409360.60.812601400.814195980.815842610.817467430.819123650.80.946901170.950165240.953485380.956842390.960258531.01.097427541.103224251.109124571.115130131.12125476Table 15
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of local preySt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.01.38E − 081.77E − 091.10E − 091.24E − 091.04E − 100.27.06E − 091.51E − 098.29E − 109.85E − 103.47E − 100.44.21E − 093.33E − 112.05E − 102.81E − 101.92E − 100.64.01E − 092.43E − 102.77E − 102.79E − 102.07E − 090.81.44E − 082.19E − 106.01E − 104.68E − 101.57E − 091.02.13E − 084.83E − 101.02E − 093.60E − 107.30E − 10Table 16
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of immigrant preyXt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.01.39E − 093.44E − 113.07E − 101.07E − 112.30E − 090.22.29E − 094.35E − 095.04E − 105.02E − 124.24E − 100.47.01E − 105.98E − 106.77E − 109.67E − 111.96E − 100.61.41E − 112.14E − 104.69E − 106.77E − 111.29E − 090.81.61E − 092.21E − 108.72E − 127.39E − 111.42E − 091.01.01E − 092.03E − 091.13E − 093.02E − 108.42E − 11Table 17
Absolute errors obtained by the LeNN-WOA-NM algorithm for population density of the predatorYt under the influence of β1.
tβ1=0.05β1=0.10β1=0.15β1=0.20β1=0.250.03.95E − 113.48E − 094.48E − 091.83E − 104.83E − 100.28.22E − 101.51E − 093.01E − 093.19E − 103.25E − 100.43.03E − 111.96E − 101.20E − 092.60E − 115.34E − 100.69.03E − 103.60E − 102.64E − 095.37E − 108.30E − 120.81.36E − 091.83E − 097.57E − 105.75E − 102.25E − 101.03.84E − 097.37E − 091.40E − 099.52E − 102.47E − 10Figure 7
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inβ1 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Figure 8
Comparison between the box plots and convergence graphs of fitness evaluation over 100 independent runs for the prey-predator model with variations in force of interaction between local and immigrant prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inβ1. (b) Convergence of the fitness value of the prey-predator model with variations in β1.
(a)(b)Table 18
Comparison through global performance indices for variations inβ1 on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI8.39E − 067.18E − 053.75E − 054.26E − 042.61E − 062.74E − 042.14E − 042.14E − 041.63E − 061.21E – 044.34E – 051.12E − 05II1.00E − 066.17E − 053.16E − 051.30E − 048.45E − 076.17E − 054.37E − 051.99E − 051.54E − 061.23E – 044.33E – 051.09E − 05III1.64E − 067.26E − 053.67E − 051.93E − 042.27E − 067.43E − 055.27E − 054.08E − 051.49E − 061.41E – 044.91E – 052.10E − 05IV1.44E − 066.46E − 053.23E − 051.17E − 041.73E − 066.86E − 054.88E − 053.59E − 051.91E − 061.23E – 044.33E – 051.55E − 05V5.98E – 075.97E − 052.96E − 057.46E − 051.07E − 066.29E − 054.45E − 052.23E − 051.63E − 061.04E – 043.72E – 059.00E − 06Table 19
Convergence analysis of population density of local, immigrant prey, and predator under the influence of variations inβ1.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI895516099836100913571384II87612561008410100942973397III988027210079121009425783712IV9877253100858100932286439V908440410081191009441836021XtI977536119977151008521946627II9981481410083111009017966927III98843312100827100909966519IV998439111008010100921471173V7727501008281008916936132YtI97753611996231009223998041II94783681006841009024977938III81374099655100872080499IV9774347996741009021814311V772750100974210010066876027Table 20
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations inβ1.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.430069150.19531180.585262780.20397331.09461748−0.30423830.451373750.6134726210.213003912−0.18170230.12967851−0.2238481−0.1299324−0.2205129−0.24340370.893839230.3475711850.1926919053−0.1639676−0.19069710.01654308−0.26832750.21586942−0.06508260.15826136−0.386370640.54484618540.6432707−0.2243964−0.0644709−0.14154830.00139286−0.12214450.21055624−0.284600750.3252831885−0.0633810.20479542−0.08377010.37097529−0.2551801−0.14574440.85822173−0.358983850.0958098676−0.05872150.09571579−0.21513340.238571310.306612510.34850165−0.0057450.6202316490.5993485570.31197992−0.0344083−0.05428880.01997741−0.1641653−0.05154490.000740050.000287462−0.3738123280.85365924−0.01697890.255722330.80492801−0.0023113−0.2852425−0.0522238−0.22752380.1463248649−0.2837938−0.069528−0.22218920.316929490.34731804−0.0409553−0.168709−0.125701740.02922959310−0.2720053−0.00832730.713032790.01306660.043472410.104905791.16491711−0.146201030.240782136110.02324570.044139090.57180799−0.39617790.13805184−0.02508270.1189892−0.099767430.676863876Case II10.200842570.018481050.011511120.188341850.077143851.001992210.330907110.5919545220.17208180120.777315120.20280783−0.1534665−0.0898678−0.02221640.050976590.39013928−0.16168208−0.874142883−0.89326870.01562423−0.8790331−0.8785872−0.24009780.011269390.62600687−0.074690290.2932000124−0.12540860.08769516−0.4580886−0.2608019−0.08707280.233163880.49916445−0.27253148−0.0084165550.05410689−0.0916165−0.4295043−0.5942252−0.03405130.03471072−0.0634013−0.11877784−0.7721202160.225648440.11645954−0.01944170.06716894−0.0067884−0.2770185−0.3074844−0.17361157−0.0294093470.257781660.028023310.328260150.0874406−0.92391470.57760720.12893563−0.074945180.1013149118−0.61522520.08467719−0.005507−0.00168850.269596760.000273920.01072972−0.059671890.7940112819−0.49165910.034230080.352213850.16735460.200979460.03970559−0.53290310.056857150.093143062100.093720710.45237283−0.42700230.122709540.01633451−0.02449120.32114085−0.085295710.896172613111.84E − 05−0.1783898−0.0967566−0.032303−0.3339006−0.05318620.052990490.002508360.0227459Case III10.328659270.73232241−0.04262260.36804480.12592558−0.17118130.68950930.3656759770.181511720.063523260.016153480.719430010.027656140.57709610.125083080.725793830.2342071480.27897832330.059020580.74414566−0.0387953−0.12858020.02597920.02360138−0.03907720.3343788710.37388573940.016729860.30780827−0.03051870.764788350.00922025−0.00031760.39416458−0.141058660.23126747550.168869060.16663601−0.15577560.14372478−0.00535840.398679360.538671370.083576797−0.09300856−0.1209806−0.37495880.216202920.046822170.16582226−0.18138350.31094709−0.118440160.0038172927−0.07560490.268685170.112320290.215547350.081722680.333028420.258253170.129740332−0.2709219380.17875086−0.01300430.7761683−0.2683637−0.05682850.08575422−0.34051820.078523738−0.103839490.521793020.034586650.064754270.069562310.250723130.027159480.259072610.038549325−0.02300298100.447951010.146474340.16811506−0.00085980.316891280.217170580.559240270.030536512−0.63332453110.07886187−0.12418510.139439470.01341213−0.62677380.34912288−0.1931707−0.012986310.579458203Case IV10.369260070.34021884−0.09357310.161692130.092468161.385538050.440585440.4999267590.12746079620.324814770.210042150.173602920.06611299−0.1524440.172302550.48717520.240390062−0.132472693−0.15700690.04398714−0.13721160.103939270.277123880.219293050.359345590.0724365430.59968018440.10042725−0.0405233−0.02761240.00768115−0.11527970.466710010.326976350.2188138030.1157854250.29815990.078954080.05929431−0.25757910.182107180.279497530.11208642−0.13550479−0.1930335560.017742780.07944380.343013220.3079287−0.02884290.352851730.08051491−0.03832890.30152455170.15448802−0.02900920.379211770.405773980.023177340.197600180.54536650.1901281440.13497225980.073032090.08335033−0.0272924−0.02684070.114950610.175658620.145480130.0504595460.99955956690.066380160.093651920.460467010.311571490.247886270.397456990.14745623−0.02038990.11010541100.145744110.270190120.08964107−0.0136530.69998870.147064260.143116930.0855874730.406653134110.00422723−0.02023260.097542644.84E − 050.011302540.45443510.20276643−0.021765530.041899718Case V10.416681550.551414170.317056470.53356928−2.0883662−0.26107040.24747335−0.495683480.1511660120.281111040.298731990.124158610.010351830.3912644−0.22508560.114638350.5725941980.6395297123−0.28800510.476324660.021482480.428317540.25753792−0.10679690.49570912−0.565745240.14021161640.283821540.074922320.25762031−0.00516160.449995−0.25836760.008623420.042826998−0.245326385−0.18289930.056100350.411578990.10823296−0.00221370.438408260.145317870.5040825650.0733579816−0.35365450.40455640.05410188−0.1039612−0.44871770.185451370.002075410.071838234−0.1819850970.23523025−0.11973780.227370420.27683837−0.45025420.53741643−0.19293150.0945797620.3349364458−0.01019840.607967930.06730913−0.07407840.15137890.63740825−0.3265244−0.272198550.2351746669−0.0038892−0.32395580.104053890.411949020.01303423−0.0585570.255794720.192432911−0.11924924100.209832440.246830530.091350390.04756230.44987648−0.0005977−0.29770780.0131149460.599953402110.176991850.19078854−0.4635923−0.0023057−0.20808410.138832930.21279433−0.199952010.280170884Figure 9
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations inβ1 on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


## 7.3. Problem III: Effect of Variations inα2 on the Prey-Predator Model

In this problem, the effect of variations in the increasing rate of immigrant preyα2 on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(32)Minimizeϵ=ϵ1+ϵ2+ϵ3+ϵ4+ϵ5+ϵ6,where ϵ1 to ϵ6 are defined as follows:(33)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−0.2SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xnα2−α2Xnk2+0.1SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Four cases are considered, depending on the value ofα2.Case I:α2=0.05Case II:α2=0.10Case III:α2=0.15Case IV:α2=0.25The LeNN-WOA-NM algorithm is used to optimize the population densities of equation (33). Table 21 represents the comparison between the Ranga–Kutta method and the proposed technique LeNN-WOA-NM. A comparison between population density of local prey St, immigrant prey Xt, and predator Yt with variations in α2 is shown in Table 22. Statistics of absolute errors are given in Table 23 and graphically presented in Figure 10. The behavior of the fitness function for different cases is shown on boxplots as demonstrated in Figure 11. The mean values of the fitness function lie around 10−7 as shown in Tables 24 and 25. Unknown parameters achieved by the LeNN-WOA-NM algorithm for the given problem are given in Table 26. From solutions (see Figure 12), the following conclusion can be drawn:(i)
Population density of local prey, immigrant prey, and predator varies directly with an increasing rate of immigrate preyα2Table 21
Comparison between solutions obtained by the LeNN-WOA-NM algorithm and ode45 for the prey-predator model withα2=0.15.
tStode45Xtode45Ytode450.00.500000600.50.499997130.50.500000850.50.20.521400430.52148360.462187510.469148750.590212060.590555110.40.542738690.543043030.419677190.432428940.694768210.696305340.60.563822410.564468040.373304840.390424020.814973650.818807100.80.584439230.585541390.324332830.344202130.952014790.959486471.00.604429480.606055270.274401980.295309561.107021861.11963769Table 22
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population densities of local, immigrant prey, and predator under the influence of the increasing rate of immigrant preyα2.
tStXtYtα2=0.05α2=0.10α2=0.15α2=0.25α2=0.05α2=0.10α2=0.15α2=0.25α2=0.05α2=0.10α2=0.15α2=0.250.00.499999960.499999700.500000600.500000180.499999970.500000230.499997130.499998700.4999999970.5000005620.5000008470.5000010330.20.521312910.521356890.521400430.521501120.453212050.457682610.462187510.471379240.5897703670.5899978480.5902120550.5906658560.40.542357970.542554470.542738690.543124360.403568430.411546230.419677190.436442150.6928196750.6937941060.6947682060.6967944650.60.562996940.563416230.563822410.564666160.352169170.362585310.373304840.395688610.8101597120.8125464260.8149736470.8200217540.80.583081610.583756220.584439230.585876620.300375870.312131320.324332830.350168880.9427785700.9473401200.9520147920.9618234251.00.602454640.603425500.604429480.606536820.249687630.261761580.274401980.301442451.0917048231.0992366981.1070218621.123532428Table 23
Absolute errors obtained by the LeNN-WOA-NM algorithm for population densities of local, immigrant prey, and predator under the influence ofα2.
tStXtYtα2=0.05α2=0.10α2=0.15α2=0.25α2=0.05α2=0.10α2=0.15α2=0.25α2=0.05α2=0.10α2=0.15α2=0.2501.38E − 081.54E − 092.74E − 091.77E – 091.18E − 091.88E − 104.68E − 096.17E − 093.17E − 106.73E − 118.97E – 102.19E − 100.26.09E − 091.40E − 093.92E − 093.57E – 096.24E − 101.22E − 106.51E − 109.47E − 091.15E − 114.00E − 119.21E – 107.22E − 110.44.30E − 091.63E − 101.22E − 093.29E – 101.59E − 098.95E − 115.03E − 091.78E − 117.33E − 101.29E − 101.48E – 104.81E − 100.62.00E − 096.47E − 102.54E − 101.89E – 095.33E − 101.49E − 112.95E − 098.36E − 101.06E − 099.40E − 125.36E – 105.34E − 100.88.76E − 097.02E − 101.16E − 093.77E – 099.98E − 119.28E − 117.06E − 103.10E − 111.31E − 102.41E − 112.99E – 104.31E − 1011.40E − 089.57E − 101.35E − 096.00E – 093.71E − 105.62E − 125.19E − 091.12E − 117.18E − 109.95E − 124.08E – 103.35E − 10Figure 10
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inα2 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Figure 11
Comparison between the box plots and convergence graphs of fitness evaluation over 100 independent runs for the prey-predator model with variations in the increasing rate of immigrant prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inα2. (b) Convergence of the fitness value of the prey-predator model with variations in α2.
(a)(b)Table 24
Statistics of global performance indices for variations in the intrinsic growth rate of immigrant prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI2.25E − 065.70E − 052.85E − 057.97E – 059.91E − 076.19E − 054.75E − 057.97E − 056.03E − 078.42E − 053.00E – 055.97E − 06II1.21E − 064.64E − 052.36E − 054.09E – 055.08E − 074.68E − 053.50E − 058.22E − 064.05E − 076.93E − 052.49E – 052.93E − 06III1.24E − 064.33E − 052.17E − 053.13E − 054.31E − 074.76E − 053.46E − 058.99E − 065.75E − 077.72E − 052.82E – 053.20E − 06IV5.22E − 074.55E − 052.29E − 054.13E − 053.65E − 074.91E − 053.39E − 051.15E − 058.39E − 078.64E − 053.06E – 054.58E − 06Table 25
Convergence analysis of population density of local prey, immigrant prey, and predator under the influence of variations in the increasing rate of immigrant preyα2.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI9462201100899100962989469II1007016110091141009826914214III996516110096111009927964711IV1008738610085201009842875321XtI9878451410085131009017956537II100895417100871110092191008638III999252111008981009717997932IV10089461610087121009621987127YtI998639121007651009527978972II10092461110079610097311009360III10083461510074810099261008756IV987843111006941009819999354Table 26
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations in the increasing rate of immigrant preyα2.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.29809866−0.03153660.580127350.472975030.168707710.254228270.044875360.4942418270.21904284220.07902509−0.142623−0.53265720.299332070.363340650.400516670.79930460.015001430.04649882930.37716273−0.0698453−0.5278093−0.59197780.06003798−0.6780052−0.6088520.6091109920.11634367240.064186860.242090730.119104830.29508708−0.00842750.0309784−0.0262412−0.34143893−0.536647465−0.0499254−0.02504030.051066860.063262910.384623390.077014170.121516760.1190981770.099603766−0.03959580.01489951−0.28834840.084806050.18102453−0.02914670.00187112−0.01144597−0.3412298770.156847550.086090640.207353430.518959580.40337040.210525920.498563420.2626537380.0345593378−0.03401640.048636210.222766390.044039820.213446940.318286370.242865030.0437745890.5295648479−0.25192170.0619910.136546150.18808650.144864010.376159850.468323640.1746203680.266799487100.54249116−0.03241190.225653230.107307140.217643420.300076760.077317380.2110310890.199472911110.129040850.04855899−0.46668450.25677349−0.0420107−0.1514589−0.17789850.104252082−0.3576282Case II10.22014881−0.09553890.019227330.440992690.169020980.210581670.566622760.214534099−0.125478432−0.5105560.025220920.565965080.2853307−0.01404990.1381917−0.00046460.4918297680.3823676313−0.16576350.01472461−0.11479830.393519120.21069293−0.26858750.395098870.072673362−0.1828391940.052823270.07711052−0.20514380.129927560.193529140.24990450.436458580.3821135730.05871400150.644992530.257626410.238082030.8495944−0.0043216−0.00171050.254571620.313158881−0.1355356160.09858966−0.02731450.487879870.093747910.21889407−0.0062504−0.13447480.044741355−0.214952670.23380371−0.12602150.32895940.05799592−0.3333640.15201786−0.4456943−0.069769270.00213826380.16035874−0.1136321−0.00365010.01897557−0.06503620.090919151.14914925−0.035859990.16591012790.016831660.030273560.013660120.167758170.200576770.14595760.538281720.1303285280.267951969100.267719360.122046520.2996044−0.11528120.188191760.17111444−0.38691720.095372589−1.0466345511−0.14301490.0042418−0.1108715−0.14027110.17676553−0.2193704−0.27470910.125441060.108988227Case III10.208334770.87977164−0.3298640.374971710.158035121.086854990.81716797−0.82852980.2076197242−0.009672−0.0834121−0.00629890.005063330.45907436−0.2830568−0.055188−1.296406970.3831862863−0.12604430.5466368−0.2547923−0.1978145−0.30147920.014824530.16026646−0.12624417−0.142116884−0.0459971−0.18508240.440100740.85262434−0.2293558−0.18033120.15934035−0.278759070.92589765350.61609669−0.0654918−0.22499821.06003476−0.06056080.0377206−0.24654680.095037929−0.003680466−0.42218770.19157710.29040894−0.1447329−0.1716272−0.28141120.102520450.1405726430.10066670870.12304762−0.0412171−0.2947437−0.33442590.15047463−0.0948965−0.3952389−0.186674990.2203322980.22471377−0.1369576−0.01359290.461183750.096751020.254991790.6880506−0.04146236−0.020095289−0.00285250.28384483−0.01379880.868556160.151195960.43141743−0.3108490.1770470290.639644264100.10012064−0.0365188−0.34958450.37524460.150820620.397374230.02610775−0.116321040.264941005110.6577301−0.17291860.28691512−0.277895−0.2620078−0.31836860.46666040.182729564−0.40213292Case IV10.275842730.144401320.05159123−0.00145010.029449330.438578970.415416940.3641066070.13763746721.03785820.22634893−0.13364490.374794840.323901380.116831911.106341910.7438702810.13372108130.305552340.141783720.098708630.1378461−0.01342310.252356050.158475340.5467479610.26755784840.49396717−0.00778540.128639560.3091184−0.04341350.984756940.07222518−0.00152740.0328421550.964611720.487610510.02679450.44318330.195001840.288898580.10389701−0.060481020.0127791256−0.0311005−0.24588980.013237040.392157780.075258550.0342729−0.02390090.3050476820.2408197467−0.00097850.278034190.434940140.08135220.165273610.10128603−0.02357030.1312367830.34259740880.032459480.055916940.9180906−0.0641160.192863010.46309077−0.02267510.2656931280.19181902290.04135314−0.0595634−0.018826−0.04538630.104730060.091887330.121393520.2816876520.039757676100.062064640.108719550.214806490.12965376−0.02754390.211693440.08952384−0.23209362−0.11227638110.41183634−0.01232830.31935826−0.06298370.007029690.084862440.225785890.1943565610.207691077Figure 12
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations in the increasing rate of immigrant prey on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


## 7.4. Problem IV

In this problem, the effect of variations in the negative force of interaction between local prey and immigrant prey on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(34)Minimizeϵ=ϵ1+ϵ2+ϵ3+ϵ4+ϵ5+ϵ6,where ϵ1 to ϵ6 are defined as follows:(35)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−0.2SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+β2SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Three cases are considered, depending on the value ofβ2.Case I:β2=0.05Case II:β2=0.10Case III:β2=0.20The LeNN-WOA-NM algorithm is used to optimize the population densities of prey-predator model equation (35). Population densities of local prey, immigrant prey, and the predator are given in Tables 27-28, while absolute errors are shown in Table 29. The absolute errors in the solutions of the proposed technique lie between 10−9 and 10−12 as shown in Figure 13. Table 30 shows global values of performance indices and convergence of the fitness value as shown in Figure 14. The statistics shown in Table 31 illustrate that the global values of different performance functions lie between 10−5 and 10−7, which highlights the robustness of the technique. The values of the weights in the LeNN structure for recreation of the approximate solutions are given in Table 32. From Figure 15, the following conclusion can be drawn:(i)
Population density of local prey varies directly with variations inβ2(ii)
Population density of immigrant prey and predator varies inversely with variations inβ2Table 27
Comparison between solutions obtained by the LeNN-WOA-NM algorithm and ode45 for the prey-predator model under the influence of variations inβ2.
tStode45Xtode45Ytode450.00.500002600.50.499996630.50.500004340.50.20.521405240.521408030.462026650.462027280.590204950.590209760.40.542730860.542742160.419188970.419193220.694726290.694733920.60.563798210.563803330.372432670.372431380.814839800.814839440.80.584394820.584398120.323085530.323078370.951682950.951689481.00.604344660.604353850.272810910.272808341.106366291.10638491Table 28
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence of the negative impact of the force of interaction on immigrant preyβ2.
tStXtYtβ2=0.05β2=0.10β2=0.20β2=0.05β2=0.10β2=0.20β2=0.05β2=0.10β2=0.200.00.499999970.499998390.500002600.499999980.500000380.499996630.499999980.499999710.5000043440.20.521474820.521466980.521405240.469148100.466766810.462026650.590550200.590446210.5902049530.40.543034740.542947310.542730860.432426210.427972540.419188970.696300310.695789680.6947262890.60.564467510.564244610.563798210.390424950.384325450.372432670.818802140.817473630.8148398000.80.585536730.585163100.584394820.344204810.337019100.323085530.959478280.956852740.9516829451.00.606046820.605480650.604344660.295308300.287602420.272810911.119627321.115148381.106366286Table 29
Absolute errors obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence of the negative impact of the force of interaction on immigrant preyβ2.
tStXtYtβ2=0.05β2=0.10β2=0.20β2=0.05β2=0.10β2=0.20β2=0.05β2=0.10β2=0.200.03.09E − 094.62E − 091.04E − 083.11E − 102.55E − 106.59E − 095.14E − 102.00E − 092.51E − 090.22.26E − 092.04E − 095.78E − 093.65E − 102.21E − 097.44E − 101.62E − 102.46E − 093.36E − 090.42.02E − 093.08E − 092.81E − 114.80E − 153.55E − 091.32E − 114.25E − 108.45E − 102.38E − 090.62.16E − 105.66E − 101.57E − 095.83E − 101.08E − 091.97E − 096.71E − 111.23E − 095.01E − 100.81.96E − 091.72E − 108.46E − 103.17E − 101.01E − 091.45E − 113.56E − 112.75E − 096.92E − 091.01.65E − 091.04E − 107.87E − 105.40E − 101.98E − 091.91E − 101.55E − 114.27E − 094.61E − 09Figure 13
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inβ2 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Table 30
Statistics of global performance indices for variations in the increasing rate of immigrant prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI5.58E − 074.70E − 052.36E − 053.69E − 054.72E − 075.09E − 053.59E − 051.10E − 058.00E − 079.48E − 053.40E – 054.71E − 06II1.48E − 065.02E − 052.49E − 055.05E − 054.76E − 075.11E − 053.63E − 051.24E − 057.09E − 079.42E − 053.35E – 055.21E − 06III2.13E − 065.09E − 052.58E − 056.78E − 058.11E − 075.28E − 053.85E − 051.18E − 057.99E − 077.69E − 052.79E – 053.93E − 06Figure 14
Comparison between the box plots and convergence graphs of fitness evaluation over 100 independent runs for the prey-predator model with variations inβ2. (a) Box plots for fitness evaluation of the prey-predator model with variations in β2. (b) Convergence of the fitness value of the prey-predator model with variations in β2.
(a)(b)Table 31
Convergence analysis of population density of local, immigrant prey, and predator under the influence of variations inβ2.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI10086323100921010010025924410II1006312010089101009926894810III966514210089151009737895215XtI10086451410088910094151006825II1009148121009071009419977330III9982421510088121009620986927YtI998039510067210094241008648II10076391410071310093231008842III9886481210077310096261008964Table 32
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations inβ2.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.275351440.31269649−0.2126230.355379880.472481030.108226330.521149510.185720875−0.0089027820.615168380.24930343−0.1133566−0.26988980.230757451.456621310.002854710.54326608−0.039481033−0.33317120.136830020.07415358−0.007558−0.14145220.00429373−0.2363309−0.53264189−0.126252194−0.1724685−0.1569281−0.1780917−0.02301890.06567477−0.1619849−0.074340.2804944630.72430884150.342352550.166015010.337796120.382142220.172063980.375675440.25394552−0.204242560.49782990260.00647555−0.06710070.418484980.4655630.178523990.259100010.62807088−0.25015769−0.0360779470.39917241−0.20864230.27916279−0.1982513−0.10032140.086122760.23112207−0.11265536−0.0463813180.37891578−0.14710920.255162720.47437651−0.04797260.39541230.233714660.1557852320.1115255199−0.15698220.21029864−0.280122−0.04986250.078988720.18498889−0.1404023−0.289746740.27339322810−0.1405838−0.0119767−0.16351250.11249313−0.06545650.69188379−0.0877791−0.1594545−0.77121428110.392969390.02457412−0.0200821−0.12014470.171618890.036300630.048113920.279716864−0.22688875Case II10.501052520.825363390.580454830.45518044−0.1160988−0.30930880.147352550.5765523590.19212727120.261418060.217962380.31815591−0.16766240.4234933−0.16140250.009758380.539177880.7758476183−0.1877622−0.1568166−0.38590220.34192011−0.11621830.787573110.40990792−0.291269520.09465586940.02599785−0.17645390.174915140.031499050.44911473−0.1695519−0.0583248−0.00447879−0.062893735−0.2923416−0.10657390.387648390.10732316−0.0035220.430096060.202684350.6275844490.023338846−0.38680090.44970016−0.0422409−0.2581683−0.41051990.016568080.000667930.075964394−0.1868564470.05739953−0.1193984−0.26071230.17181358−0.27136280.66311474−0.124996−0.079813080.5059381328−0.02163280.80670030.31701813−0.15115880.140927880.65770582−0.1918301−0.004073960.1210127269−8.30E − 07−0.0664668−0.11686640.63180611−0.1478398−0.15926770.270268430.1131427720.003053122100.09718671−0.0187463−0.0760148−0.04440810.299805990.00951493−0.1584802−0.011772520.821398692110.193967170.18738185−0.4407164−0.1400116−0.15198690.038500450.32211631−0.106169950.295182432Case III10.368776140.41870330.005517210.371063780.76553227−0.11350280.367662840.6565438770.07401287820.571231580.390292350.269251950.02897603−0.0177996−0.27249760.68341531−0.566598910.0352719053−0.2786739−0.1350866−0.40798240.18172611−0.1092217−0.65352410.19817369−0.181044410.12771873540.248960520.133223240.342458990.185354870.183366550.32712799−0.14222260.1528342880.02504536250.291759670.21119283−0.19578160.31132875−0.0532529−0.33059640.04517261−0.73533086−0.2105123960.386092490.275908470.58209501−0.535581−0.18238940.31125028−0.00330380.0826663750.32301675570.01164667−0.14217680.243283170.264800370.021735020.078484850.617754820.2318933740.07434687880.593813210.06375022−0.36897130.35998444−0.03686660.329072190.148773580.1434832470.33440679−0.2877027−0.08129760.585135370.114541340.04186864−0.38143560.28795352−0.00041389−0.1333593110−0.30739940.0111124−0.00688850.04196657−0.1275650.14629639−0.40819820.00302860.33443128110.06421859−0.09197220.34881249−0.0684058−0.14121790.4979470.198741510.0711107530.298703303Figure 15
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations inβ2 on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


## 7.5. Problem V: Effect of Variation inγ1 on the Prey-Predator Model

In this problem, the effect of variations in the catching rate of local preyγ1 on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(36)Minimizeε=ε1+ε2+ε3+ε4+ε5+ε6,where ϵ1 to ϵ6 are defined as follows:(37)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−0.2SnXn+γ1SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+0.1SnXn+0.9XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Four cases are considered, depending on variations inγ1.Case I:γ1=0.005Case II:γ1=0.015Case III:γ1=0.020Case IV:γ1=0.025The LeNN-WOA-NM algorithm is used to optimize the population densities of equation (37). Table 33 represents the comparison between the Ranga–Kutta method (ode45) and the proposed technique LeNN-WOA-NM algorithm. Approximate solutions for population density of local prey St, immigrant prey Xt, and the predator Yt with different values of γ1 are presented in Table 34. Statistics of absolute errors are given in Table 35 and graphically shown in Figure 16. The minimum and mean values of the fitness function for each case lie between 10−9 and 10−7, respectively as illustrated in Figure 17. Convergence analysis of the fitness value, MAD, TIC, and ENSE is given in Tables 36-37. Unknown neurons of LeNN are shown in Table 38. From Figure 18, the following conclusion can be drawn:(i)
Population density of local prey and population density of the predator vary inversely with variations inγ1(ii)
Population density of immigrant prey has a direct relation with variations inγ1Table 33
Comparison between solutions obtained by the LeNN-WOA-NM algorithm and ode45 for the prey-predator model withγ1=0.015.
tStode45Xtode45Ytode450.00.500000290.50.499999710.50.500000020.50.20.521179670.521174760.46675620.466764390.590433250.590425180.40.542302620.542298910.427974970.427979250.695720320.6957040.60.56315320.563151080.384347740.384354420.817264190.817258110.80.583510510.583506110.337055510.337060360.956375380.956370131.00.603152380.60315130.287692460.287699271.114209391.11419232Table 34
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence of variations inγ1.
tStXtYtγ1=0.005γ1=0.015γ1=0.020γ1=0.025γ1=0.005γ1=0.015γ1=0.020γ1=0.025γ1=0.005γ1=0.015γ1=0.020γ1=0.0250.00.499999820.500000290.500000000.499999970.499999820.499999710.499999980.500000170.5000005040.5000000180.4999999970.5000001310.20.521747000.521179670.520892690.520606610.466751740.466756200.466760070.466766150.5904495820.5904332540.5904113170.5903940500.40.543594720.542302620.541659740.541016310.427954660.427974970.427976180.427997140.6958466350.6957203220.6956306740.6955587970.60.565339290.563153200.562061770.560975590.384308290.384347740.384378340.384397480.8176786560.8172641910.8170462120.8168373420.80.586808480.583510510.581862690.580224420.336956440.337055510.337112360.337160850.9573164870.9563753780.9558944400.9554196061.00.607822510.603152380.600836190.598526840.287517770.287692460.287779660.287874201.1160620941.1142093921.1132594931.112332477Table 35
Absolute errors obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence of variations inγ1.
tStXtYtγ1=0.005γ1=0.015γ1=0.020γ1=0.025γ1=0.005γ1=0.015γ1=0.020γ1=0.025γ1=0.005γ1=0.015γ1=0.020γ1=0.0250.02.74E − 098.47E − 103.60E − 107.80E − 115.74E − 092.00E − 094.47E − 094.19E − 128.54E − 106.23E − 099.79E − 125.10E − 130.22.76E − 094.69E − 102.92E − 101.04E − 113.59E − 098.16E − 114.78E − 095.77E − 105.53E − 106.28E − 091.38E − 118.72E − 110.45.70E − 104.68E − 104.62E − 133.63E − 102.51E − 093.02E − 104.26E − 102.89E − 105.49E − 108.44E − 101.36E − 122.06E − 100.61.89E − 091.24E − 101.54E − 102.49E − 124.04E − 113.74E − 104.26E − 092.55E − 111.47E − 123.54E − 099.17E − 111.14E − 090.82.36E − 091.80E − 112.01E − 128.31E − 101.28E − 094.93E − 103.67E − 095.03E − 103.59E − 103.49E − 099.80E − 114.43E − 101.02.82E − 091.85E − 104.54E − 121.17E − 092.58E − 101.22E − 154.68E − 091.08E − 102.82E − 104.48E − 093.19E − 101.19E − 09Figure 16
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inγ1 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Figure 17
Comparison between the box plots and the convergence graph of fitness evaluation over 100 independent runs for the prey-predator model with variations in the catching rate of local prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inγ1. (b) Convergence of the fitness value of the prey-predator model with variations in γ1.
(a)(b)Table 36
Statistics of global performance indices for variations in the catching rate of local prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI2.50E − 065.21E − 052.69E − 052.50E − 069.20E − 074.63E − 053.33E − 051.51E − 059.19E − 078.88E − 053.10E − 058.98E − 06II7.32E − 066.20E − 053.14E − 055.63E − 052.43E − 066.86E − 054.92E − 057.58E − 051.44E − 061.06E − 043.80E − 051.70E − 05III2.86E − 064.87E − 052.48E − 059.34E − 057.60E − 075.39E − 053.89E − 051.86E − 051.39E − 068.53E − 053.09E − 055.88E − 06IV1.78E − 064.62E − 052.38E − 057.52E − 055.18E − 074.32E − 053.10E − 051.46E − 056.42E − 078.01E − 052.83E − 059.08E − 06Table 37
Convergence analysis of population density of local prey, immigrant prey, and predator under the influence of variations inγ1.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI957020510085211009437855621II94672529981231009443835721III957117210090191009642915519IV977829610091261009742905726XtI979059221009115100942797838II97865624998416999128967840III9985641710088151008921948133IV9993572210093201009631988146YtI99875115997771009525999159II96784815997381009330988762III9882491510075101009730998861IV9890601710093371009337999364Table 38
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations inγ1.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.10433411−0.9093816−0.05513840.10814164−0.0975817−0.40247590.110579010.030277215−0.076830932−0.1015821−0.23390630.11183609−0.00796490.12409457−0.00663660.031204470.0659841510.1698154793−0.17322190.05666544−0.3418627−0.15762190.205116340.23158339−0.3370419−0.154182380.19415431340.28706660.04587599−0.0750170.0317597−0.1567915−0.29606070.126614180.176984003−0.0830713250.18918051−0.26348620.092969590.177127140.066902930.055961240.06180474−0.14442377−0.9894472560.08011991−0.02418430.043841940.09323427−0.10295410.30404503−0.1508149−0.050515030.1910929767−0.4734878−0.18889670.139719480.249251060.002252150.00501846−0.9616715−0.197021870.01134063180.19753354−0.0151743−0.0358307−0.2646482−0.4318811−0.2034545−0.0870373−0.088720140.07936334790.073523210.043742590.17641082−0.004594−0.0709633−0.1885710.012284110.051586673−0.1707669100.05966701−0.0094660.11306597−1.0101745−0.11955560.30440396−0.165510.0727288860.05362364110.244192780.215236270.125015060.108878510.12112026−0.1899001−0.157371−0.159061210.128530508Case II10.425624360.224051350.025039560.35262391−0.22536640.262594550.271116860.3287966910.22892787220.416669880.5078151−0.18890930.121713910.209981050.38048580.303907470.3891165320.46070179630.105776950.060028460.21768786−0.53933980.390275590.139022260.473504190.4654198640.66197000740.465174070.09698879−0.16334780.110142970.197399070.66549310.05604469−0.003926970.27877739450.601000290.019320570.568304980.23543140.414040570.20437334−0.1267241−0.131288940.29384493960.081277360.106281460.386109420.235399620.171875290.385790850.352251370.0676494280.18848928170.354984110.044902120.049006640.064765870.352380740.502358060.362536390.203702979−0.079204880.384970770.018235750.34815044−0.07555940.494782790.380354950.16102938−0.124828140.3405400119−0.05617650.73291315−0.18902150.0053794−0.44388960.466410960.2889986−0.017214940.10349715710−8.42E − 05−0.01914780.330575130.011386710.210127450.280139280.417905120.0110224340.111373073110.03758884−0.00909620.298419490.061814350.282112330.10443310.76029441−0.053057410.400559476Case III10.34210927−0.1686215−0.06371520.234271170.252859750.392147380.269000940.5655202880.4852402582−0.0413096−0.00813010.106145370.60106519−0.07658960.421885660.341462510.5755837060.08123241630.29904262−0.075813−0.19675630.459391190.15935586−0.0737513−0.03714280.3957801950.52602507540.31966978−0.02968460.134928190.744499180.138230990.147330570.312965730.1703563410.44903041650.385111710.027429150.192458680.174809460.260677330.19256129−0.19667440.0266527150.64855673760.15405406−0.0855950.039002980.20316026−0.09393630.369996350.141163030.183493540.15728401470.011602360.19042424−0.0099165−0.1955586−0.01452490.51154810.092342830.0501830530.5295766948−0.0119424−0.00172220.19894821−0.04590190.169954560.099121220.442771620.2681979320.3604521989−0.0135090.028267380.13366696−0.09861210.125556210.29387726−0.03515110.15868406−0.02656018100.36233153−0.05044810.464221880.331541310.044883490.444237770.182334570.0026998190.25017744110.207808520.061875820.190128340.47309423−0.00014550.285446120.326557470.2322044930.038376682Case IV10.08743197−0.02745860.138384540.10216597−0.6205394−1.16254730.08890874−0.468214890.03190664120.85138918−0.25960890.27697631−0.0395679−0.31904840.47807164−1.07668990.1945224680.05147395530.32484395−0.03201280.049102590.59721285−0.7402654−0.5159275−0.54570820.5845145650.33420834241.265395290.0637656−0.32985420.007467950.12675579−0.21199380.06894660.0780600630.34124362250.27124475−0.0145706−0.8795409−0.2258574−1.04891440.18138843−0.7172783−0.446545880.2883523536−0.0615221−0.03571810.05167451−0.0093528−0.90763030.005299290.111646850.0841619050.0213049777−0.93304190.005870160.08238699−0.00384620.5318456−0.0029243−0.56249830.07126723−0.071945938−0.22228510.159644890.159293470.02677628−0.08595810.71257834−0.15642590.011627909−0.011415649−0.07467320.03595603−0.16882950.918538690.0432707−0.19964060.9708232−0.207186870.491048316100.111059850.03339826−0.0412419−0.9111904−0.07572990.383182960.140461090.0682333990.44790321911−0.9956577−0.03257510.958620420.01858614−0.03282860.2342197−0.8047726−0.035479630.018056695Figure 18
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations inγ1 on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


## 7.6. Problem VI: Effect of Variations inγ2 on the Prey-Predator Model

In this problem, the effect of variations in the catching rate of immigrant prey on population densities of the prey-predator model is discussed. An error-based fitness function along with initial populations is given as follows:(38)Minimizeε=ε1+ε2+ε3+ε4+ε5+ε6.where ϵ1 to ϵ6 are defined as follows:(39)ϵ1=1N∑n=134dSndt−Sn0.12−0.12Sn50−0.2SnXn+0.01SnYn2,ϵ2=1N∑n=3568dXndt−Xn0.2−0.2Xnk2+0.1SnXn+γ2XnYn2,ϵ3=1N∑n=69102dYndt−0.9SnYn−0.8XnYn+0.01Yn+0.01Yn22,ϵ4=S0−0.52,ϵ5=X0−0.52,ϵ6=Y0−0.52.Five cases are considered, depending on the value ofγ2.Case I:γ2=0.7Case II:γ2=0.8Case III:γ2=1.0Case IV:γ2=1.1The LeNN-WOA-NM algorithm is applied to prey-predator model equation (39) to study the influence of variations in the catching rate of immigrant prey. The results calculated by the designed algorithm are compared with those of the Runge–Kutta method using ode45 in MATLAB as shown in Table 39. Approximate solutions for population densities of local prey, immigrant prey, and the predator are given in Table 40. Absolute errors are presented in Table 41 and graphically shown in Figure 19. The minimum fitness values in Figure 20 reflect the accuracy of the solutions by the proposed algorithm. The convergence of our numerical approach is assessed by fitness values and statistical analysis of the performance indicators (MAD, TIC, and ENSE), see Table 42. Table 43 shows the correctness of the proposed algorithm to tackle real-world problems. Trained neurons in LeNN are shown in Table 44. From Figure 21, the following conclusions can be drawn:(i)
Population densities of local prey, immigrant prey, and the predator vary inversely with variation in the catching rate of immigrant preyγ2.Table 39
Comparison between solutions obtained by the LeNN-WOA-NM algorithm and ode45 for the prey-predator model withγ2=0.8.
tStode45Xtode45Ytode450.00.500000060.50.50000050.50.499999820.50.20.521510970.521174760.471866320.466764390.590685680.590425180.40.543161820.542298910.4381850.427979250.696940280.6957040.60.564756030.563151080.399355260.384354420.820536990.817258110.80.586083320.583506110.356171160.337060360.963191650.956370131.00.606933670.60315130.309865090.287699271.126496341.11419232Table 40
Approximate solutions obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence ofγ2.
tStXtYtγ2=0.7γ2=0.8γ2=1.0γ2=1.1γ2=0.7γ2=0.8γ2=1.0γ2=1.1γ2=0.7γ2=0.8γ2=1.0γ2=1.10.00.500000380.500000060.499999990.499999840.500000480.50000050.500000040.499999990.4999997790.4999998180.4999999840.4999996970.20.52156930.521510970.521410350.521350990.477021170.471866320.461725790.456737880.5909334430.5906856820.5902027520.5899718050.40.54339380.543161820.5427270.542506280.448656410.4381850.4180060.408282360.698110510.6969402820.6946498760.6935324810.60.565282850.564756030.563745210.563255510.415003050.399355260.36990490.356058780.8236969380.8205369860.8144934660.8115992330.80.587048620.586083320.5842640.583396270.376511730.356171160.318968250.301952150.9698202460.9631916450.9507666550.9449383661.00.608471360.606933670.604096540.602772190.33402430.309865090.267081790.248146241.1385322641.1264963371.1044111611.094260679Table 41
Absolute errors obtained by the LeNN-WOA-NM algorithm for population densities of local prey, immigrant prey, and predator under the influence ofγ2.
tStXtYtγ2=0.7γ2=0.8γ2=1.0γ2=1.1γ2=0.7γ2=0.8γ2=1.0γ2=1.1γ2=0.7γ2=0.8γ2=1.0γ2=1.10.01.26E − 091.32E − 104.79E − 105.02E − 101.48E − 111.87E − 107.79E − 105.53E − 109.67E − 101.58E − 095.89E − 101.90E − 090.21.31E − 094.54E − 111.43E − 105.38E − 101.04E − 104.83E − 128.26E − 102.09E − 103.34E − 111.58E − 096.20E − 102.77E − 100.47.22E − 103.13E − 124.65E − 107.47E − 117.86E − 103.53E − 102.65E − 097.30E − 101.16E − 093.29E − 107.51E − 111.87E − 090.69.27E − 102.41E − 121.09E − 101.51E − 103.37E − 126.92E − 101.12E − 092.57E − 101.29E − 091.11E − 095.35E − 101.29E − 090.82.37E − 095.90E − 112.34E − 112.37E − 109.11E − 101.53E − 101.69E − 103.50E − 113.98E − 101.50E − 095.10E − 101.19E − 101.02.91E − 091.44E − 101.20E − 122.88E − 109.82E − 101.61E − 101.58E − 094.63E − 115.13E − 101.68E − 097.98E − 104.24E − 10Figure 19
Absolute errors obtained by the LeNN-WOA-NM algorithm under the influence of variations inγ2 on population densities of the prey-predator model. (a) Population density of local prey. (b) Population density of immigrant prey. (c) Population density of the predator.
(a)(b)(c)Figure 20
Comparison between the box plots and convergence graphs of fitness evaluation over 100 independent runs for the prey-predator model with variations in the catching rate of immigrant prey. (a) Box plots for fitness evaluation of the prey-predator model with variations inγ2. (b) Convergence of the fitness value of the prey-predator model with variations in γ2.
(a)(b)Table 42
Convergence analysis of population density of local prey, immigrant prey, and predator under the influence of variations inγ2.
CasesFITMADTICENSE≤10−6≤10−7≤10−8≤10−9≤10−4≤10−5≤10−6≤10−4≤10−5≤10−6≤10−5≤10−6≤10−7StI956524510088201009541895920II986424410092181009945936018III866312410090161009532905516IV986720410092211009744926221XtI9985522210085201009230946932II10087572210087231009127987333III9883461310082101009316987128IV9990541910094171009723998244YtI98844417100771110092281008757II998147189975131009731998755III998441151007561009426998956IV9985501410072910094341008755Table 43
Statistics of global performance indices for variations in the catching rate of local prey on the prey-predator model.
CasesStXtYtGFITGMADGTICGENSEGFITGMADGTICGENSEGFITGMADGTICGENSEI1.98E − 064.80E − 052.43E − 055.78E − 055.99E − 075.04E − 053.39E − 052.68E − 059.28E − 079.13E − 053.24E − 055.44E − 06II2.17E − 063.83E − 051.95E − 052.91E − 055.21E − 074.70E − 053.30E − 051.44E − 051.36E − 069.17E − 053.25E − 051.01E − 05III2.16E − 064.95E − 052.52E − 055.55E − 051.01E − 065.58E − 054.16E − 051.53E − 057.10E − 078.79E − 053.17E − 056.81E − 06IV1.55E − 064.01E − 052.00E − 054.85E − 055.75E − 074.49E − 053.40E − 057.86E − 067.28E − 078.59E − 053.08E − 055.50E − 06Table 44
Unknown parameters obtained by the LeNN-WOA-NM algorithm for the prey-predator model under the influence of variations inγ1.
IndexStXtYtζnψnθnζnψnθnζnψnθnCase I10.27328231−0.0407082−0.00365680.190954050.196199130.253340580.154480320.3391994560.64626032920.002036930.007612690.246550490.67156041−0.12021850.410298230.354566580.167129547−0.0222010830.22854565−0.1965182−0.07997060.4234230.03898304−0.1160494−0.19322110.3272984570.58963410940.43068766−0.0047135−0.00677490.76144510.244368120.103860.45804558−0.105041050.09663247950.42381358−0.15026520.518274420.107579960.246969080.11302049−0.04908280.0544971930.54146435560.030643450.387909−0.09083380.20021796−0.10933340.498037490.085022350.2923121810.1689517997−0.02562720.181349440.04903339−0.14391810.016780890.015222870.122286450.0406107450.5042415958−0.00292340.053691650.161186290.092383710.35203171−0.07424580.582963550.3263946220.31123246990.01433551−0.00628260.114617840.000702810.081706550.27442616−0.06249750.046728012−0.00253833100.52132655−0.07458240.323731330.316258890.064695840.355187160.432659940.001675421−0.07508196110.04700892−0.14311380.190960340.610866830.058795620.071452470.316027320.2945685460.005492958Case II10.228880980.156627270.721183980.663963010.48400554−0.00149870.2930431−0.36246252−0.1982755620.320391690.213497160.26632024−0.1848203−0.0673308−0.06482460.131252990.0416570280.0299543743−0.08293880.030948320.003035220.62951247−0.19254380.24736655−0.19631210.4315314350.7080957674−0.1658907−0.1447463−0.6955688−0.232498−0.26160950.270252270.22782911−0.28937041−0.129944965−0.01747580.232151820.042900160.12334482−0.2128439−0.0360696−0.14592740.1433832820.1276542126−0.0712296−0.3626918−0.03486680.26004380.131684030.451310950.73463526−0.091098940.1000586577−0.0244551−0.0112202−0.1246274−0.34457320.31169855−0.06807980.148097380.0295262120.1760398998−0.29639480.276843650.02435043−0.1570692−0.10580240.01048308−0.02870220.0522047370.5002902529−0.0250215−0.10242610.056514850.61926963−0.2055780.43411127−0.37484230.146360873−0.4197865510−0.28473110.13336075−0.0903463−0.1464593−0.3554768−0.04750720.16657978−0.04184543−0.0441026211−0.00129140.02947063−0.2027891−0.01192550.472381010.094487350.42979806−0.05005385−0.00793862Case III10.17136911−0.6220119−0.06594510.160013380.06604026−0.62594620.417434420.347451460.12340272620.132480930.224071980.16846536−0.01740040.460482060.161721960.63635670.0278926630.3492171153−0.0296660.0085174−0.4187633−0.05734020.296738520.13113638−0.7934564−0.181700450.08725763940.509163990.04634821−0.14370190.08923057−0.2567021−0.29347460.26333479−0.21308607−0.0250535150.28783823−0.26621630.07660280.099144850.576658110.043000830.20662283−0.48985711−1.0974074360.08917337−0.10500090.238182560.02035607−0.10721920.2659477−0.0029764−0.131483110.3136653667−0.7503685−0.17044710.18719634−0.0233035−0.00889380.26343414−1.0737247−0.190186990.28010052480.342300890.00558946−0.1427919−0.2957063−0.64904810.105464590.273225470.041606539−0.1027815790.258797360.219539270.110529190.00206003−0.0680416−0.21899140.03786660.055450179−0.26693807100.114828690.042566210.24740165−1.02102880.108268520.10811710.289721620.107963625−0.01416421110.216458910.26405307−0.2682646−0.02512170.2673504−0.2407415−0.3201356−0.126449070.363928738Case IV10.163210441.001017390.446108510.19117524−0.18691310.032844490.487094060.4928390050.06881143120.050340140.274485020.055954270.34759862−0.0428559−0.46566320.412112910.5503501620.26125445730.023098730.00700133−0.3300406−0.032857−0.2059139−0.15721320.430692780.1197089380.1320454840.60643185−0.16770850.013515050.012693450.07414495−0.2533687−0.19908810.074189084−0.0087421950.06424168−0.03942680.040623410.103460710.24294981−0.14069910.15305910.007966010.1423245696−0.1607398−0.0067892−0.1142057−0.3542569−0.3179260.896468750.19029585−0.07772530.1382154397−0.0249375−0.1421790.211814490.00745531−0.0807929−0.0606121−0.3138917−0.04179489−0.0839599580.012718540.07904033−0.20324210.30678256−0.0970233−0.0973761−0.2681465−0.26023333−0.0574878590.157291660.020040370.084417750.50899351−0.01071120.16547740.05321702−0.55038421−0.01142195100.455268140.308482240.03885589−0.0670780.03709569−0.11034190.00132508−0.0601008−0.00809993110.00144682−0.0160469−0.14856550.07678165−0.0117808−0.039538−0.15537240.090185082−0.05537005Figure 21
Comparison between solutions obtained by the LeNN-WOA-NM algorithm under the influence of variations inγ2 on population densities of the prey-predator model. (a) Approximate solutions for population density of local prey. (b) Approximate solutions for population density of immigrant prey. (c) Approximate solutions for population density of the predator.
(a)(b)(c)


## 8. Conclusion

In this paper, we have analyzed the system of the singular differential equation (SDE) representing the phenomena of the prey-predator model with immigrant prey. Generally, solving SDE is one of the challenging tasks, and therefore, we have designed a novel soft computing technique known as a LeNN-WOA-NM algorithm. Weighted Legendre polynomials are used to model approximate series solutions for the prey-predator model with variations in various coefficients, including the growth rate of local and immigrant preyα1,α2, force interaction between local and immigrant prey β1,β2, and the catching rate of local and immigrant prey γ1,γ2. We summarize our findings as follows:(i)
Variations inα1,α2,β1, and β2 has a direct impact on population density of local prey St, while population density of local prey is inversely affected by variations in γ1 and γ2.(ii)
Variations inα1,α2, β2, and γ2 has an inverse impact on population density of immigrant prey Xt, while population density of immigrant prey is directly affected by variations in β1 and γ2.(iii)
Variations inα1,α2 and β1 has a direct impact on population density of the predator St, while population density of the predator is inversely affected by variations in γ1, β2, and γ2.(iv)
Approximate solutions obtained by the LeNN-WOA-NM algorithm are compared with those obtained by the homotopy perturbation method and MATLAB solver ode45. The results show the dominance of the proposed technique.(v)
Lower absolute errors in our solutions and convergence analysis of fitness evaluation, MAD, TIC, and ENSE show the accuracy and robustness of the proposed algorithm for obtaining solutions to real-world problems.In future, this approach can be utilized to solve the complex nonlinear systems of fractional differential equations characterizing real-world problems, for instance, anomalous diffusion of contaminant from the fracture into the porous rock matrix, microbial survival and growth curves, smoking dynamics, parametric identification of Hammerstein systems with time delay, and asymmetric dead zones.



---

*Source: 1241761-2022-10-13.xml*