a=0.04; b=0.003; g=0.0015; i0=0.06; modell={v'[t]==v[t] (a-v[t]) (v[t]-1) - w[t] + i0, w'[t]==b v[t] - g w[t], v[0]==v0, w[0]==w0}; v0=0.1; w0=0; tmax=1000; p:={loesung=NDSolve[modell,{v,w},{t,0,tmax},MaxSteps->5000]; bild=ParametricPlot[Evaluate[{v[t],w[t]}/.loesung],{t,0,tmax}, PlotRange->All, DisplayFunction->Identity]}; p; bild1=bild; w0=0.1; v0=0.4; p; bild2=bild; w0=0.15; p; bild3=bild; w0=0.25; v0=1.2; p; bild4=bild; w0=0.025; v0=-0.25; p; bild5=bild; bild=Show[bild1,bild2,bild3,bild4,bild5,DisplayFunction->$DisplayFunction, AxesLabel->{"v","w"},ImageSize->{400,300}];