b=floor((50)^(1/2)+10); c=floor((200)^(1/2)+10); e=floor(10-(50)^(1/2))+1; f=floor(10-(200)^(1/2))+1; d=0; for x=e:b for y=f:c if ((((x-10)^2)+((y-10)^2))<35)||((((x-10)^2)+((y-10)^2)/4)<50) d=d+1; end end end disp(d)