Quick and simple cubic polynomial fit to surface - no checks
where:
z = c[0] + c[1].x + c[2].y + c[3].x^2 + c[4].y^2 + c[5].xy
quadfit()