input: 
  import numpy as np
  from sklearn.linear_model import LinearRegression
  
  def linear_regression(X, Y):
      """Perform linear regression using the given dataset and predict the value of y for a new input.
      >>> linear_regression([1,2,3,4,5], [3,5,7,9,13])
      [15.]
      """
      X = np.array(X).reshape(-1, 1)
      Y = np.array(Y)
      
      model = LinearRegression()
      model.fit(X, Y)
      
      # Predicts the value of y for the new input
      return model.predict([[7]])
  
output:     X = np.array(X).reshape(-1, 1)
      Y = np.array(Y)
      
      model = LinearRegression()
      model.fit(X, Y)
      
      # Predicts the value of y for the new input
      return model.predict([[7]])