fn cubic_bezier(f: f64, x1: f64, y1: f64, x2: f64, y2: f64) -> f64
y for the x = f on the cubic bezier (0,0) (x1,y1) (x2,y2) (1,1) — Newton iterations on the x polynomial.
f