Numerical solution of multiple subdomain diffusion-reaction model with fully-implicit and linear implicit approach is implemented in the coupled dune-copasi framework.