input: 
  import numpy as np
  from sklearn.linear_model import LinearRegression
  
  X = np.array([1,2,3,4,5]).reshape(-1,1)
  Y = np.array([3,5,7,9,13])
  
  model = LinearRegression()
  model.fit(X, Y)
  
  # Predicts the value of y
  print(model.predict([[7]]))
  
output:     
  # Calculates the coefficient of determination (R^2)
  print(model.score(X, Y))
  