how to generate a function from a linear model

Hi All,

I am trying to generate a function from a linear model. I think there 
should be build-in function that perform this action but I've had no 
luck finding it.

For example, I have a model created using lm().

	model <- lm(sat.d~1+sat.n+I(sat.n^2))

What I would like to have is a function (similar to the one generated 
by splinefun()) so that I can use it on different data-sets.

Thanks in advance for the help.

Tony Han Bao
tonybao at mac.com