ISSN:
0748-8025
Keywords:
Engineering
;
Engineering General
Source:
Wiley InterScience Backfile Collection 1832-2000
Topics:
Mathematics
,
Technology
Notes:
A second-order explicit method is developed for the numerical solution of the Ricatti (logistic) initial-value problem u′ ≡ du/dt = αu(1 - u), t 〉 0, u(0) = U0, in which α ≠ 0 is a real parameter. The method is based on two first-order methods which appeared in an earlier paper by the authors (Twizell et al.1). In addition to being chaos-free and of higher order, the novel method is seen to converge to the correct, stable, steady-state solution for any value of the parameter α, provided the denominator of the method does not vanish. Convergence is monotonic or oscillatory depending on the magnitude of the product αl, where l is the parameter in the discretization of the independent variable t. This dependence of the type of convergence on αl is likened to the behaviour of the well known Crank-Nicolson method for solving the simple heat equation. Conversion of the numerical method to a reliable, empirical model for predicting the limited growth of successive generations of a population is given. When extended to the numerical solution of Fisher's equation, in which the quadratic polynomial αu(1 - u) appears as the reaction term, the numerical solution is found by solving a linear algebraic system at each time step, as opposed to solving a non-linear system, which often happens when solving non-linear partial differential equations.
Additional Material:
1 Ill.
Type of Medium:
Electronic Resource
URL:
http://dx.doi.org/10.1002/cnm.1630080805
Permalink
