a=0.04; b=0.003; g=0.0015; i0=0; 
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->1000];
bild=ParametricPlot[Evaluate[{v[t],w[t]}/.loesung],{t,0,tmax},
PlotRange->All, DisplayFunction->Identity]};
p; bild1=bild;
v0=0.06;  p; bild2=bild;
v0=0.08; p; bild3=bild;
v0=0.1; p; bild4=bild;
v0=0.5;  p; bild5=bild;
bild=Show[bild1,bild2,bild3,bild5,DisplayFunction->$DisplayFunction,
AxesLabel->{"v","w"},ImageSize->{400,300}];