grayplot of a surface defined by a function

deff('[z]=surf(x,y)','z=x**2+y**2');
fgrayplot(-1:0.1:1,-1:0.1:1,surf);
xtitle('fgrayplot ',' ',' ');