Welcome to MathOverflow

Question

0

1

Can you use the Crank-Nicolson method to get a numerical approximation to the fisher-kolmogorov equation? If not what would be the easiest way to model the equation using matlab?

Thanks and sorry its so basic but i'm new to these sorts of numerical approximations

flag	asked Jan 12 2012 at 18:20 Tony Flowers 5●1

magnetohydrodynamic · Accepted Answer · 2012-01-19 17:17:43Z UTC

Crank-Nicolson will work fine, though a fully implicit scheme may actually perform a bit better especially if you have discontinuous initial conditions. Since Crank-Nicolson is effectively half explicit, discontinuities in the initial conditions $u(x,0)$ propagate into a lot of noise echoing around your solution grid.

Make sure you choose $\Delta t$ sufficiently small that your propagation matrix remains diagonally dominant. That is to say, choose it small in comparison to the initial size $|u(x,0)|$.

S. Mabuza · Answer 2 · 2012-01-19 17:32:12Z UTC

Another approach especially in higher dimensions would be to use a Fourier - Spectral method for spatial discretization followed by a fourth order Runge - Kutta algorithm to the algebraic differential system which is in "Fourier Space". The Inverse Fast Fourier Transform can then be used to recover the solution in the original domain. This technique is useful especially with coupled reaction diffusion equations.

Numerical solution to Fisher-Kolmogorov equation

Remember to vote up questions/answers you find interesting or helpful (requires 15 reputation points)

2 Answers

You can accept an answer to one of your own questions by clicking the check mark next to it. This awards 15 reputation points to the person who answered and 2 reputation points to you.