Fit a polynomial shape, imposing f(x,y)>0

Hi, I would like to fit 3 series of experimental data (x,y, and z) using a polynomial surface.

It works easily with the command "fit". f = fit( [x, y], z, 'poly32' ) How can I force the fitted surface to be always positive ("z(x,y) > 0")

Thanks