x = linspace(-pi,pi,7);
[xx,yy] = meshgrid(x);
u = cos(sqrt(xx.^2+yy.^2));
u(3,:) = NaN;

figure(1)
contourf(xx,yy,u)
xlabel('x'); ylabel('y');