surf.gls function

Fits a Trend Surface by Generalized Least-squares