This code:

```
eqs =
{Ca/u^2 +
r (6 a (c^2 - d^2) + 12 b c d + 6 (a^2 + b^2 + c^2 + d^2) Ca +
Ca^3) == 0,
a (1/u^2 - 1) +
r (3 a (a^2 + b^2) + 6 a (c^2 + d^2) + 3 c^2 Ca - 3 d^2 Ca +
3 a Ca^2) == 1/2,
b (1/u^2 - 1) +
r (3 b (a^2 + b^2) + 6 b (c^2 + d^2) + 6 c d Ca + 3 b Ca^2) == 0,
c (1/u^2 - 1/4) +
r (6 (a^2 + b^2) c + 3 c (c^2 + d^2) + 6 (a c + b d) Ca +
3 c Ca^2) == 0,
d (1/u^2 - 1/4) +
r (6 (a^2 + b^2) d + 3 d (c^2 + d^2) + 6 (-a d + b c) Ca +
3 d Ca^2) == 0, c > 0, d >= 0,
QQQ == c^2 + d^2, QQQ != 0
} // Rationalize(#, 0) &;
NSolve(eqs /. {u -> 5, r -> 0.04}, {a, b, c, d, Ca, QQQ}, Reals)
```

Gives me this answer:

```
{{a -> -0.732167, b -> 0, c -> 1.06332, d -> 0, Ca -> 0.443594,
QQQ -> 1.13066},
{a -> -0.732167, b -> 0, c -> 1.06332, d -> 0,
Ca -> 0.443594, QQQ -> 1.13066},
{a -> -0.732167, b -> 0,
c -> 1.06332, d -> 0, Ca -> 0.443594,
QQQ -> 1.13066}, {a -> -0.732167, b -> 0, c -> 1.06332, d -> 0,
Ca -> 0.443594, QQQ -> 1.13066}, {a -> -0.698614, b -> 0,
c -> 0.622043, d -> 0.622043, Ca -> 0,
QQQ -> 0.773876}, {a -> -0.698614, b -> 0, c -> 0.622043,
d -> 0.622043, Ca -> 0, QQQ -> 0.773876}, {a -> -0.698614, b -> 0,
c -> 0.622043, d -> 0.622043, Ca -> 0,
QQQ -> 0.773876}, {a -> -0.698614, b -> 0, c -> 0.622043,
d -> 0.622043, Ca -> 0, QQQ -> 0.773876}, {a -> -0.698614, b -> 0,
c -> 0.622043, d -> 0.622043, Ca -> 0,
QQQ -> 0.773876}, {a -> -0.698614, b -> 0, c -> 0.622043,
d -> 0.622043, Ca -> 0, QQQ -> 0.773876}, {a -> -0.698614, b -> 0,
c -> 0.622043, d -> 0.622043, Ca -> 0, QQQ -> 0.773876}}
```

But as you see, there are several identical answers.

And if a make WorkingPrecision -> 3, for example, this code is slow down.

Maybe there is other methods to make this code faster?