- #1

- 367

- 13

## \frac {\partial A}{\partial t} = aAC + bBD - cE ##

## \frac {\partial B}{\partial t} = - c(E+F) ##

## \frac {\partial C}{\partial t} = dAD - eE ##

## \frac {\partial D}{\partial t} = dBC + - eF ##

## \frac {\partial E}{\partial z} = cCD - fE ##

## \frac {\partial F}{\partial z} = cCD - fF ##

(Upper case letters stand for the solutions sought; lower case letters represent constants)

I've tried existing packages such as FiPy and tried my own basic solver using FEM, but neither has been successful aside from a few nontrivial cases. With the two independent variables (z and t) especially, this seems to complicate the solver. Any advice you have on existing libraries or numerical methods I should implement would be of great interest. (I am most familiar with Python, but am open to responses for any programming language.)