x=-50:0.5:.50;
y=x;
[X,Y]=meshgrid(x,y);
Z=(cos(X))^2+(cos(Y))^2;
pcolor(X,Y,Z);
shading interp
axis('equal','square','off')
