I am trying to use a explicit FDM for transient 1d conditions with linear elements for specific diffusion equation:
ds/dt=D(x)*d2s/dx2+R(s)
the problem is that I am using different, nonlinear functions describing diffusion constant D(x) and reaction rate R(s).
The diffusion parameter is...