k=0.6; l=0.3; eps=0.001; end=0.004;
modell={u'[t]==-u[t](u[t]+k-l)v[t],eps v'[t]==u[t]-(u[t]+k)v[t],
u[0]==1,v[0]==0};
loesung=NDSolve[modell,{u,v},{t,0,end}];
abb=Plot[{u[t]/.loesung,v[t]/.loesung},{t,0,end},
          PlotRange->{0,1},
          PlotStyle->{{},Dashing[{0.01}]},
          AxesLabel->{"s","u,v"},
          ImageSize->{400,300}];