Spatially distributed Model I (Related to the Model 3 in Holmes et al (2017) R= Rac, P = Rho, E = ECM Rac-Rho equations d R/dt= {b_R }/{1+\rho_a^3}R_I - \delta R +D_R Laplacian( R), d \rho/dt ={b_\rho }{1+R_a^3}\rho_I - \delta \rho +D_\rho Laplacian( Rho) ECM equation dE/dt= epsilon*((K+GammaR*f_R1())-E*(kP+GammaP*f_P1())) where the f's are Hill functions of Rac or Rho.

1 dx*space.x (bR/(1+P^n))*RI-1*delta*R (kE+GammaE*f_E1())*PI/(1+R^n)-1*P E^3/(E0^3+E^3) epsilon*((K+GammaR*f_R1())-E*(kP+GammaP*f_P1())) R^n/(R0^n+R^n) P^n/(P0^n+P^n) (-bR/(1+P^n))*RI+1*delta*R -(kE+GammaE*f_E1())*PI/(1+R^n)+1*P E min(R,2)