Calculating R^2 value for the data that fall within 95% prediction interval

Question

I have data as follows.

import numpy as np
import matplotlib.pyplot as plt
import scipy.stats as stats

x = np.array([50,52,53,54,58,60,62,64,66,67,68,70,72,74,76