Licson Lee says to Ow< [MSFW]
double smoothstep(double x) { return x * x * (3 - 2 * x); }