|
Playing a bit with the function.
The very strange thing is that I get a numerical instability if I set
nu=0. This means that the problem has nothing to do with the advection
term. It is odd since a pure diffusion equation really shouldn't be
unstable under the Crank-Nichoson solver.
|