How to compute (and plot) the confidence intervals of linear regression estimators?

Question

I have some data x_data, y_data and I estimate the linear regression line using the package stats from skipy. My code is this