input: from typing import List
  
  def linear_regression_prediction(X: List[float], Y: List[float]) -> float:
      """
      Perform linear regression on given dataset to predict the value of y for a new x.
      
      >>> linear_regression_prediction([1,2,3,4,5], [3,5,7,9,13])
      15.5
      """
output:     n = len(X)
      sum_x = sum(X)
      sum_y = sum(Y)
      sum_xy = sum(x * y for x, y in zip(X, Y))
      sum_x_squared = sum(x ** 2 for x in X)
      
      m = (n * sum_xy - sum_x * sum_y) / (n * sum_x_squared - sum_x ** 2)
      b = (sum_y - m * sum_x) / n
      
      x_new = X[-1] + 1
      y_pred = m * x_new + b
      
      return y_pred