In[]:=
n=750;range=Range[n];cf[x_]:=Blend[{Darker@Blue,White},x]munchingSquare=Outer[BitXor,range,range];(*plotstandardMunchingsquare*)ArrayPlot[munchingSquare,ColorFunction->cf,PlotLabel->"750 x 750 Munching square"](*Createabinarymatrixwhereblackpixel=primevalueandwhitepixel=nonprimevalue*)plot=(Map[Boole@*PrimeQ,munchingSquare,{2}]//ArrayPlot[#,PlotLabel->"750 x 750 Munching square, \n where Black = prime and White = nonprime"]&)
Out[]=
Out[]=
