When a linear relationship has been fitted by least squares, the methods for securing a prediction interval for the response at some fixed value of the independent variable are explained in many statistical text books. This paper describes the somewhat more complex problem of determining the joint prediction interval for the responses at each ofKseparate settings of the independent variables when allKpredictions must be based upon the original fitted model.