I am trying to solve a coupled reaction-diffusion equations, using Crank Nicolson (implicit Finite Differences Method). I know how to solve them separately, but not simultaneously.
I transform it to Crank Nicholson form and I obtain:
Is the equation for Cb(x,t) correct? How can I couple them? I know how to solve them separately applying boundary conditions (Dirichlet both sides for Cf), and Neumann both sides for Cb.
Thank you very much!