input: def predict_with_linear_regression(X: List[int], Y: List[int]) -> np.ndarray:
      """Perform linear regression on given dataset to predict numbers.
      >>> predict_with_linear_regression([1,2,3,4,5], [3,5,7,9,13])
      array([15.])
      """
output:     import numpy as np
      X = np.array(X)
      Y = np.array(Y)
      X_mean = np.mean(X)
      Y_mean = np.mean(Y)
      numerator = np.sum((X - X_mean) * (Y - Y_mean))
      denominator = np.sum((X - X_mean) ** 2)
      slope = numerator / denominator
      intercept = Y_mean - slope * X_mean
      prediction = slope * (X[-1] + 1) + intercept
      return np.array([prediction])