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.Model1RacRhoECMPDEsIn1D1dx*space.x(bR/(1+P^n))*RI-1*delta*R(kE+GammaE*f_E1())*PI/(1+R^n)-1*PE^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*PEmin(R,2)