// inefficient cuberoot function float cbrt(float x) { return pow(x, 1.0/3.0); }