input: Write a python function to solve a quadratic equation given in the form of ax^2 + bx + c.
  assert solve_equation(2, -3, 1) == (1.0, 0.5)
  assert solve_equation(1, -2, 1) == (1.0,)
  assert solve_equation(1, 0, -1) == (1.0, -1.0)
output: import math
  
  def solve_equation(a, b, c):
      discriminant = b**2 - 4*a*c
      if discriminant > 0:
          root1 = (-b + math.sqrt(discriminant)) / (2*a)
          root2 = (-b - math.sqrt(discriminant)) / (2*a)
          return (root1, root2)
      elif discriminant == 0:
          root = -b / (2*a)
          return (root,)
      else:
          return ()