I have a system of non-linear coupled PDEs, taken from a paper from the 1980s which I would like to numerically solve. I would prefer not to use a numerical Package like MatLab or Mathematica, though I will if I need to.(adsbygoogle = window.adsbygoogle || []).push({});

I would like to know if anyone knows how to solve non-linear coupled PDEs numerically or can point me to a text book/reference which can explain how to do so. I am most familiar with finite difference methods, so it would be preferable if I could get an algorithm which used a finite difference method, but I am flexible.

The system of equations are

∂b/∂t = dΔb + εnb - a(b,n)b

∂s/∂t = a(b,n)b

∂n/∂t = Δn - nb

where b(x,t), s(x,t), and n(x,t) are all functions of space and time; a(b,n) is some decreasing function of b and n (something simple, but not constant); ε and d are constants. The initial and boundary conditions are will be either dirichlet or neumann and the initial conditions are simple.

Any help would be appreciated. Thank you for your time.

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Finite Difference Numerical Solution to NL coupled PDEs

**Physics Forums | Science Articles, Homework Help, Discussion**