How can I find 95% confidence interval for predicted value based on a multiple linear regression model?

Question

I have a model that I create in the following way:

X = df[["Height", "Age", "Gender"]]
Y = df["Weight"]
reg = linear_m