c1=1; c2= 1; c3=2; c4=0.1; c5=0.1; c6=1; b0=0.6;r0=0.3;l=150;abst=0.01; modell:={b'[t]==b[t] (c1 (1 - b[t]/c2) - (c3 r[t])/(b[t] + c4) ), r'[t]==r[t] c5 ( 1 - (c6 r[t])/b[t] ), b[0]==b0, r[0]==r0}; loesung=NDSolve[modell,{b,r},{t,0,l}]; abb=Plot[Evaluate[{b[t],r[t]}/.loesung],{t,l/4,l} ,PlotStyle->{{},{Dashing[{0.01}]}} ,AxesLabel->{"t","b,r"} ,ImageSize->{400,300}];