This function is a helper function for  roots_quartic, so it is not needed to make checks on the domain
of the coefficients; I renamed it  `_roots_quartic_euler ` to  make this clear; I also changed its arguments to avoid
useless computations. The input cannot be multivariate.