f = '1+0.3*sin(4*pi*(x+y)/((1+sqrt(3))/2))' ;
ys = funcdefn(f,5) ;
is = (1:3^5)*3 - 2 ;
g = defnstep(ys(is,3),7,5) ;
