fd_predator_prey_test, a Fortran77 code which calls fd_predator_prey(), which applies the finite difference method (FDM) to estimate solutions of a pair of ordinary differential equations that model the behavior of a pair of predator and prey populations.
The information on this web page is distributed under the MIT license.
fd_predator_prey, a Fortran77 code which applies the finite difference method (FDM) to estimate solutions of a pair of ordinary differential equations that model the behavior of a pair of predator and prey populations.
TRF_100 uses 100 timesteps, which are not enough.
TRF_1000 uses 1000 timesteps; the solution does not explode, and seems to show periodicity, except that it is clearly growing.
TRF_10000 uses 10000 timesteps. The cyclic nature of the solution is clear.