soln[μ_, 
          t_] := NDSolve[{y''[x] - μ (1 - y[x]^2) y'[x] + 
                  y[x] == 0, y[0] ==  .01, y'[0] == 0}, 
            y[x], {x, 0, t}][[1, 1, 2, 0]];

plot[t_, m_] := 
      Block[{$DisplayFunction = Identity}, 
{ParametricPlot[
            Evaluate[{x, soln[m, t][x]}], {x, 0, t}, 
            PlotStyle -> Red, 
AxesLabel -> 
              TraditionalForm /@ {x, y[x]}, 
            PlotLabel -> 
              TraditionalForm[μ == 
                  m]], 
ParametricPlot[
            Evaluate[{soln[m, t][x], D[soln[m, t][x], x]}], {x, 0, 
              t}, AspectRatio -> Automatic, PlotStyle -> Red, 
            PlotPoints -> 100, 
            AxesLabel -> 
              TraditionalForm /@ {y[x], y'[
                    x]}, 
PlotRange -> 
              All]
}]

Show[
      GraphicsArray[
plot[100, #] & /@ { .2, 1, 
            5}
], GraphicsSpacing -> {- .15, 0}, 
      ImageSize -> 600]

solnb[ϵ_, f_, ν_, 
          t_] := NDSolve[{y''[
                    x] - ϵ (1 - y[x]^2) y'[x] + 
                  y[x] == ϵ f Cos[(1 + ϵ ν) x], y[0] ==  .01, y'[0] == 0}, 
            y[x], {x, 0, t}][[1, 1, 2, 0]];

plotb[t_, ϵ_, ν_, ff_] := 
      Block[{$DisplayFunction = Identity}, 
{ParametricPlot[
            Evaluate[{x, solnb[ϵ, ff, ν, t][x]}], {x, 0, t}, 
            PlotStyle -> Red, 
AxesLabel -> 
              TraditionalForm /@ {x, y[x]}, 
            PlotLabel -> 
              TraditionalForm[
                f == ff]], 
ParametricPlot[
            Evaluate[{solnb[ϵ, ff, ν, t][x], 
                D[solnb[ϵ, ff, ν, t][x], x]}], {x, 0, t}, 
            AspectRatio -> Automatic, PlotStyle -> Red, 
            PlotPoints -> 500, 
            AxesLabel -> 
              TraditionalForm /@ {y[x], y'[
                    x]}, 
PlotRange -> 
              All]
}]

plotC[t_, ϵ_, ν_, ff_] := 
      Block[{$DisplayFunction = Identity}, 
{ParametricPlot[
            Evaluate[{x, solnb[ϵ, ff, ν, t][x]}], {x, 0, t}, 
            PlotStyle -> Red, 
            PlotPoints -> 500, 
AxesLabel -> 
              TraditionalForm /@ {x, y[x]}, 
            PlotLabel -> 
              TraditionalForm[f == ff]]
}]

Show[
      GraphicsArray[
plotC[500, 0.1, 0.4, #] & /@ {1.0, 
            1.3,  1.5}
], 
      GraphicsSpacing -> {- .15, 0}, 
      ImageSize -> 500]