Finds the real roots of ax³ + bx² + cx + d = 0. Shifting x by b/(3a) removes the squared term, and the remaining t³ + pt + q = 0 is solved according to the sign of the discriminant. There is always at least one real root.