y = a x^3 + b x^2 + c x + d