I am looking to solve some coupled pdes with neumann zero-flux bc and random ic. Ideally I would like to have a 2d or 3d solution but 1d will suffice for now. The two species involve a functional diffusion model which considers the difference in diffusion between the two species. This difference is then considered by a sigmoidal function to define positive or negative flow. There is also exponential decay of each of the species.
du/dt = - a1*u + 1/(1+exp(a2*d2u/dx2-a3*d2v/dx2)) - 1/2
dv/dt = - a4*v + 1/(1+exp(a5*d2u/dx2-a6*d2v/dx2)) - 1/2
Does anyone know of any off the shelf tools in matlab that can do this? I have the code running in mathematica but it is painfully slow.
Thanks for your help in advance :)
No products are associated with this question.