I am trying to solve the following system of PDE's
dn := diff(n(x,t),t) = d*diff(n(x,t),[x,x])+n(x,t)*(r(n(x,t)+m(x,t)))*s[1](C[1](x,t))-n(x,t);
dm := diff(m(x,t),t) = d*diff(m(x,t),[x,x])+m(x,t)*(r(n(x,t)+m(x,t)))*s[1](C[1](x,t))-m(x,t);
dc1 := diff(C[1](x,t),t) = d[1]*diff(C[1](x,t),[x,x])+P[1]+(n(x,t)+m(x,t))*p[1](n(x,t)+m(x,t))+H*m(x,t)*p[1](n(x,t)+
m(x,t))-(P[1]+p[1](1))*C[1](x,t);

with initial and boundary conditions
IBC := n(x,0)=piecewise(x^2<1,0,x^2>1,1),m(x,0)=piecewise(x^2>1,0,x^2<1,1),C[1](x,0)=1,
C[1](infinity,t)=1,n(infinity,t)=1,m(infinity,t)=0,n(-infinity,t)=1,m(-infinity,t)=0,