p=0; for i=1:1000 for j=1:1000 if (i-10)^2+(j-10)^2<35 or (i-10)^2+((j-10)^2/4)<50 p=p+1; end end end