n=0; for i=-50:1:50 for j=-50:1:50 if or(((i-10)^2+(j-10)^2)<35,((i-10)^2+(j-10)^2/4)<50)==1 n=n+1; end end end fprintf('n=%d\n',n)