Expressing the roots of a cubic as polynomials in one root All roots of $8x^3-6x+1$ are real. (*)
The discriminant of $8x^3-6x+1$ is $5184=72^2$ and so the splitting field of $8x^3-6x+1$ has degree $3$.
Therefore, all three roots can be expressed as polynomials in any one given root.
Indeed, if $a$ is a root, then the others are $2a^2-1$ and $-2a^2-a+1$.

This can be easily checked. But how can we find these expressions in the first place?

I've tried this: let $b,c$ be the other roots. Then from Vieta's formulas we get $$b+c=-a, \qquad bc=-\dfrac{1}{8a}$$
The corresponding quadratic has discriminant $3-3a^2$ but it is not at all clear that this is the square of a polynomial in $a$. (It is $(4 a^2+a-2
)^2=(b-c)^2$, as it should be.) I'm stuck here.

More generally, is there an algorithm that
  given a cubic with integer coefficients and
  having real splitting field of degree $3$,
  expresses all three roots as polynomials in any one given root?

(*) From the triple-angle formula $\sin (3\theta) = - 4\sin^3\theta + 3\sin\theta$ when $\sin (3\theta) = 1/2$, these roots are $\sin(10^\circ)$, $\sin(50^\circ)$, $\sin(-70^\circ)$, but perhaps this is immaterial here.

Based on the answers in this question, I'm able to finish the general case.
Let $x^3+px+q$ be a cubic with rational coefficients and having real splitting field of degree $3$. Let its roots be $a,b,c$. Then
$$b+c=-a, \qquad bc=-\dfrac{q}{a}$$
and $b,c$ are the roots of the quadratic
$$
h(x)=(x-b)(x-c)=x^2+ax-\dfrac{q}{a}
$$
The discriminant of $h$ is $b-c$. Now comes the nice idea from those answers:
$$
d=(a-b)(a-c)(b-c)=h(a)(b-c)
$$
where $d^2=-4p^3-27q^2$ is the discriminant of the original cubic. The hypothesis on the cubic implies that $d$ is rational. Therefore,
$$
b-c = \dfrac{d}{h(a)}
$$
To write this as a polynomial in $a$, we solve the linear system on $A,B,C$ implied by
$$
d=(Aa^2+Ba+C)h(a)=(Aa^2+Ba+C)(2a^2-\dfrac{q}{a})
$$
or
$$
ad=(Aa^2+Ba+C)(2a^3-q)=(Aa^2+Ba+C)(-2pa-3q)
$$
The solution is
$$
A=-\dfrac{6p}{d},
\quad
B=\dfrac{9q}{d},
\quad
C=-\dfrac{4p^2}{d}
$$
which gives
$$
b,c = -\dfrac{a}{2} \pm \dfrac{1}{2}(Aa^2+Ba+C)
$$
