# Combining 2-D plots for different values of a parameter into 3-D plot

So I have a system of non-linear differential equations that depend on a parameter $$\mu$$. What I'm trying to do is take the 2-D phase portraits for different values of $$\mu$$ and combine them into a 3-D plot with $$\mu$$ on the 3rd axis. The range of $$\mu$$ that I need is $$(-2,4)$$. The problem is that I have no idea how to go about it.

$\mu =-1$" />

$\mu =0$" />

$\mu =1$" />

Generation of phase portrait for a specific value of $$\mu$$

\[Mu] = 1;
eqn1 = x'[t] - y[t] + x[t] ((x[t])^2 + (y[t])^2 - \[Mu]) == 0;
eqn2 = y'[t] + x[t] + y[t] ((x[t])^2 + (y[t])^2 - \[Mu]) == 0;
eqn3 = x[0] == .01;
eqn4 = y[0] == .01;
sol = NDSolve[{eqn1, eqn2, eqn3, eqn4}, {x, y}, {t, 0, 100}];
ParametricPlot[Evaluate[{x[t], y[t]} /. sol], {t, 0, 100},
PlotRange -> All, AxesLabel -> {"X", "Y"}, PlotStyle -> Red,
PlotLabel -> "\[Mu] = 1"]


You can use ParametricNDSolveValue as follows:
ClearAll[x, y, μ, pndsv]

ParametricPlot3D[Evaluate[Table[pndsv[i], {i, -2, 4, 1}]], {t, 0, 100},
`