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}];