input: Write a function to optimize a given multivariate function using the scipy.optimize.minimize method. The function should take the multivariate function and an initial guess as inputs and return the optimized parameters.
  assert optimize_function(lambda x: x[0]**2 + x[1]**2, [1, 1])==[0.0, 0.0]
  assert optimize_function(lambda x: (x[0]-1)**2 + (x[1]-2)**2, [0, 0])==[1.0, 2.0]
  assert optimize_function(lambda x: x[0]**2 + x[1]**2 + x[2]**2, [3, 4, 5])==[0.0, 0.0, 0.0]
output: from scipy.optimize import minimize
  
  def optimize_function(func, initial_guess):
      result = minimize(func, initial_guess)
      return result.x