Using fec to smooth a grayplot

deff('[z]=surf(x,y)','z=x**3+y');
Sfgrayplot(-1:0.1:1,-1:0.1:1,surf);