c = 0; for x = 2:18 for y = 4:16 if (x-10)^2+(y-10)^2<35 c = c + 1; elseif (x-10)^2+((y-10)^2)/4<50 c = c+1; end end end