We propose an accurate non numerical solution of the Fisher Equation (FE), capable of reproducing the known analytical solutions and those obtained from a numerical analysis. The form we propose is based on educated guesses concerning the possibility of merging diffusive and logistic behavior into a single formula.
In this paper, we propose a “quasi-exact” solution of the Fisher equation and compare it with solutions available in the literature. We also develop an extensive numerical analysis along with an accurate comparison with the solution presented here.
The Fisher Equation (FE) introduced to describe the spreading of genes [
In the above equation F is a function of the spatial coordinate and time and the initial condition is given by g(x), supposed to be a continuous infinitely differentiable function. From the dynamical point of view the equation describes a spatial diffusion in a homogeneous medium (the coefficient D is assumed to be independent of the spatial coordinate) embedded with a growth of logistic nature [
with g(x) such that the integral in Equation (1.2) converges.
In the case of reduces to
In absence of the diffusion, (D = 0), Equation (1.1) describes a purely logistic process2 and its solution reads [
Where g0 is the initial condition associated with Equation (1.1) and is characterized by an initial exponential growth successively undergoing a saturation, when the F values approaching the carrying capacity. We can provide a tentative solution for the general case, by merging the characteristics implicit in Equations (1.3) and (1.4) through the following expression
possessing both features of spreading and logistic saturation. The above solution and the relevant derivation may sound as the consequences of a naive procedure without any rigor. We will provide a more sounded justification in the concluding section using an appropriate Hopf-Cole transform, here we will follow a pragmatic approach and accept it on the basis of its agreement with already exact known solutions and with full numerical solutions. Examples of evolution at different times for an initially Gaussian distribution are given in
The distribution have been normalized at the value taken at x = 0 at different times and the growth of the amplitude at the origin and for x = 7 is shown in
As we will see in the forthcoming section, the process described by Equation (1.1) can be viewed as a kind of travelling wave front at a given velocity, within the frame-
work of the previous solution the front velocity can be evaluated as
The behavior of the velocity vs. the time is reported in
while for short “times”
We have so far shown that our guessed solution is compatible with the phenomenology of the FE, in the forthcoming sections we will compare it with known analytical solutions and with numerical solutions obtained with an ad hoc developed code.
Various forms of exact solutions for the Fisher equation have been provided in the past. They are limited to specific cases of the initial function and/or of the constant entering the equation itself. Albeit limited to specific cases these solution provides an important test for comparison with our quasi exact form. One of these solutions, proposed in [7-9], is given by
and C is an arbitrary constant. To make a comparison we use Equation (1.5) with
The results of the integration are given in
The logistic behaviour of the solution (2.1) is shown in
In a more recent paper [
The comparison with the quasi exact form is shown in
Other examples could be discussed, but will not be reported here, the general consistency of our semi-analytical approach is confirmed at least for the simple form endowed in Equation (1.1). In the forthcoming section we will consider a comparison with a full numerical procedure along with some concluding remarks.
As already stressed, the exact solutions of the Fisher equation, used in the previous section as benchmarking of the validity of Equation (1.5), refer to specific cases of initial conditions. They are of limited usefulness in the application because they refer to conditions hardly applicable to problems concerning evolution of populations, genes or saturation mechanisms in laser Physics.
We will therefore further benchmark our solution with a series of numerical test using an ad hoc developed code.
The solution of Equation (1.1) will be afforded numerically by means of an iterative procedure employing a symmetric split decomposition [11-13]. We reconsider the Equation (1.1) in the form
where Λ is assumed to be a function of the spatial coordinate (as well as of time) but independent of the function F. The formal solution of our problem can accordingly be written as
Since the function Λ is depending on the spatial coordinate we cannot use a naïve disentanglement of the exponential, a convenient solution based on a symmetric decomposition scheme is eventually provided by
Where o(δ3) specifies the error of the procedure, which is associated (among other things discussed in the following) with the integration step δ.
If we divide the integration interval in N interval each one of size δ, the solution of our problem (3.2) can be cast in the form
which yields the following recursion
We can now go back to the Fisher equation and make the following substitution in Equation (3.5)
Thus finally ending up with the “linearized solution”
From the dynamical point of view we are left with a diffusive evolution problem counteracted by a spatially dependent growth rate. An example of comparison between numerical and solution with the previously outlined method and Equation (1.5) is given in
The error is accordingly specified by the discretization step and by the commutator involving the operators A, B, namely
Abrupt variations of the function itself at different evolution steps may provide a large error even for small δ, furthermore a large number of iterations may determine an accumulation of error which can be computed from
The integrals regarding the Weierstrass transform are explicitly derived using a Simpson scheme. The homemade code even though accurate is not optimized and we have used repeatedly the combinations of two numerical errors due to the split operator and to the Simpson integration method. To make sure of the comparison we have attempted a second numerical procedure, based on FFT algorithm implemented in Mathematica.
A 3-D perspective of the numerical integration is given in
In this case too the comparison between analytical and numerical solutions yields a satisfactory agreement.
The solution we have proposed, albeit naïve, seems to be rather efficient in solving the problem of providing the evolution of a population governed by the Fisher equation. Our solution is based on the requirement that it should provide the logistic and diffusive behavior, implicit in the assumption underlying its original derivation in [
It should be not surprising that this family of equations can be solved in a simple form since the logistic equation, even though non-linear can be reduced to a linear form by an appropriate change of variable and more in general of Burger type equations, which can be reduced to a heat
equation by means of a Hopf-Cole transformation [
The following example can be useful to support the previous argument, by considering the differential equation
In this case we are not dealing with a diffusive Fisher equation but with a simpler form involving a logistic saturation and a coordinate translation3. According to the same arguments invoked for the FE, we “guess” that the solution of Equation (3.11) writes4
namely the composition of a translation and of a logistic. The above solution can be obtained in rigorous terms after setting5
The transformation (3.13) eliminates the non-linearity in Equation (3.11) and this allows the solution in exact form (see [
Regarding the Fisher equation, the presence of the second derivative does not allow the elimination of the nonlinearity, the corresponding equation reads indeed
which contains a non-linearity which can be neglected whenever
which is not verified for large values of F(x, t), namely in the region in which the saturation occurs.
In a forthcoming investigation we will extend the analysis to multi-dimensional (2-D, 3-D) configurations and discuss specific examples of applications.
The authors express their sincere appreciation to Prof. F. Pacella for the kind interest and for a critical reading of the manuscript.