let root3 a =
  	let rec help x =
  		if abs_float(x ** 3. -. a) <= 1e-15 *. abs_float(a) then x
  		else help (x +. ( a /. (x ** 2.) -. x) /. 3.)
  	in help (if a >= 1. then a /. 3. else a)