From: Morgan on
Does MATLAB have the capability to find numerical approximations for systems of partial differential equations in two spatial variables as well as time?
Specifically I'm interested in a pair of reaction-diffusion equations each involving a simple laplace operator:

dF/dt = (d^2F)/(dx^2)+(d^2F)/(dy^2) + Z(F,G)
dG/dt = (d^2F)/(dx^2)+(d^2F)/(dy^2) + W(F,G)

Additionally, if I do need to use only one spatial variable, is it possible to use pdepe defining a value of F(0,0)=c (a constant), or can pdepe only be used with initial and boundary conditions (ex F(t,0)=1 and F(0,x)=0)?

Is there any way to find numerical/graphical approximations for this system given the following constraints?
G(0,x,y) = 2
F(0,0,0) = c (some positive constant)