input: Write a python function to find the solutions of a quadratic equation.
  assert find_solutions(1, 5, 6) == (-3.0, -2.0)
  assert find_solutions(1, -7, 10) == (5.0, 2.0)
  assert find_solutions(1, -1, -6) == (3.0, -2.0)
output: def find_solutions(a, b, c):
      import math
      discriminant = b**2 - 4*a*c
      sol1 = (-b + math.sqrt(discriminant)) / (2*a)
      sol2 = (-b - math.sqrt(discriminant)) / (2*a)
      return (sol1, sol2)