这是由于方程的解有好多不连续的点造成的,
将其排除即可
sol = Solve[(y (x^2 + 1) - mu)^2 == (x^2 - 2 x + 1) (1 +
1/2 (1 + mu^2) y^2 - 2 mu y) &&
y == (1 - x) mu/((x - 3) x^2 + (x + 1)), {y, x}, Reals];
Plot[Evaluate[y /. sol], {mu, -1, 1}, PlotRange -> {0.1, 0.5},
Frame -> True, AxesOrigin -> {-1, 0.1},
Exclusions -> {-0.845, -0.61, 0, 0.61, 0.845}]
将其排除即可
sol = Solve[(y (x^2 + 1) - mu)^2 == (x^2 - 2 x + 1) (1 +
1/2 (1 + mu^2) y^2 - 2 mu y) &&
y == (1 - x) mu/((x - 3) x^2 + (x + 1)), {y, x}, Reals];
Plot[Evaluate[y /. sol], {mu, -1, 1}, PlotRange -> {0.1, 0.5},
Frame -> True, AxesOrigin -> {-1, 0.1},
Exclusions -> {-0.845, -0.61, 0, 0.61, 0.845}]