Solve a system of coupled PDEs

I have a set of equations,

Unless specified as functions, the other terms are constant, f is a noise term.

These have boundary conditions given by

Answer 1

[Alan Moses]

You may refer to the pdepe function and Solve System of PDEs.