let solve_deg_2 (a,b,c) =
  let delta = ((b ** 2.0) -. 4.0 *. a *. c) in
  !print_check (Printf.sprintf "solve_deg_2\na=%f" a);
  !print_check (Printf.sprintf "b=%f" b);
  !print_check (Printf.sprintf "c=%f" c);
  !print_check (Printf.sprintf "delta=%f" delta);
  let x1 = ((-. b) -. (sqrt delta)) /. (2.0 *. a) in
  let x2 = ((-. b) +. (sqrt delta)) /. (2.0 *. a) in
  let check x =
    !print_check (Printf.sprintf "a %f^2 + b %f + = %f"
                     x x
                     (a*. (x**2.0) +. b *. x +. c))
  in
  check x1;
  check x2;
  (`Real x1, `Real x2)