input: Write a python function to find the optimal solution of a given cost function using scipy.optimize.minimize.
  assert find_optimal_solution() == (2.0, 1.0)
output: def find_optimal_solution():
      from scipy.optimize import minimize
  
      def cost_function(x):
          return (x[0] - 2)**2 + (x[1] - 1)**2
  
      initial_guess = [0, 0]
      result = minimize(cost_function, initial_guess)
      return tuple(result.x)