Using numpy.polyfit to plot a line of best fit and extract uncertainty of the fitting/covariance matrix

Question

I have plotted some data with the x-axis variables as d_list and the y-axis variables as v_list. I would like to plot a line of best fit from this, extract the gradient of t