If an incomplete block design is used to compare v treatments and later it is decided to fit a second order response surface using the fact that the v treatments are v combinations of levels of m different factors, how will one proceed to obtain this? This paper demonstrates how the usual analysis of an incomplete block design with or without recovery of interblock information can be augmented to yield estimates of response surface coefficients.