let epsilon = 1e-9

let sin_taylor x =
  let rec accumulate s y k =
    if abs_float y < epsilon then s else