x=-.5:.005:.5;
y=x;
[X,Y]=meshgrid(x,y);
Z=sin(1./(X.^2+Y.^2+eps));
pcolor(X,Y,Z);
shading interp
axis('equal','square','off')
