[求助]这种偏微分方程组的数值模拟如何着手

[wptom]

d *y(x,z)'' + a* y(x,z)'=0 (一般的扩散方程)
dy/dt=y*(1-y)/(b-y)-c*y
其中b>1
a c d是常数，0<y<1
这个数值模拟如何着手？