let atan(angle:float) (n:int) : float =
  let result = ref angle;
  for i=1 to n do
        result := !result + (1 - ((i mod 2) * 2 )) *. gen_atan_term(angle, n);
  done;
  !result;;