input: Write a function to perform linear regression for predicting the numbers in a given dataset. The function should take as input two lists, X and Y, representing the independent and dependent variables, respectively. The function should return a prediction for a given value of X.
  
  import numpy as np
  from sklearn.linear_model import LinearRegression
  
  def linear_regression(X, Y):
      X = np.array(X).reshape(-1,1)
      Y = np.array(Y)
      model = LinearRegression()
      model.fit(X, Y)
      return model
  
  model = linear_regression([1,2,3,4,5], [3,5,7,9,13])
  
  assert model.predict([[7]]) == [15.0]
output: def linear_regression(X, Y):
      X = np.array(X).reshape(-1,1)
      Y = np.array(Y)
      model = LinearRegression()
      model.fit(X, Y)
      return model