l=10;
r=0.55;q=8;c=0.2;
f[x_]:=r x (1-x/q)-x^2/(1+x^2);
modell={u'[s]==v[s],
v'[s]==-c v[s] - f[u[s]],u[0]==0,v[0]==0.7};
loesung=NDSolve[modell,{u,v},{s,0,l}];
abb=Plot3D[u[x-c t]/.loesung[[1]],{x,5,10},{t,0,5},
PlotPoints->30,AxesLabel->{"x","t","u"},ImageSize->{400,300}];