clear;
n=100; a=1; b=2;
A=repmat([-n:n].^2,2*n+1,1);
B=(sqrt(a*A+b*A')<n-1)+(sqrt(a*A+b*A')>n+1);
figure(1); clf; imagesc(B);