How to incorporate the uncertainty of the model coefficients in the prediction interval of a multiple linear regression
One possible solution is Bayesian linear regression. Bayesian linear regression estimates a posterior distribution for each coefficient. From that posterior distribution, a credible interval can be calculated.