p(x,y)=cos(x*sin(y)) n(x,y)=floor(p(x,y)/sin(x*cos(y))) f(x,y)=0.5*(sin(256*x)+sin(256*y))-p(x,y)+n(x,y)*sin(x*cos(y)) g(x,y)=f(x,y) > 0 ? 1: 0 set pm3d map unset colorbox set palette gray set size square unset key unset xtics unset ytics set samples 512 set isosamples 512,512 set term png tiny crop R=2*pi set xrange[0:R] set yrange[0:-R] # set out "g011.png" splot g(x,y) # set out "f011.png" splot f(x,y) # set palette color set out "f011c.png" splot f(x,y) |