eq:a*x^3+b*x^2+c*x+d = 0;
sol:solve(eq,x);
cubic(a,b,c,d) := ev([rhs(sol[1]),rhs(sol[2]),rhs(sol[3])]);