clc; clear;
xc=input('Xc= ');yc=input('yc= ');
r=input('r= ');
hold on;
i=0;
trace=zeros(1,8);
axis([-100 100 -100 100])
x=0; y=r;
while x<=y
plot(xc+x,yc+y,'g.');
plot(xc+x,yc-y,'b.');
plot(xc-x,yc+y,'k.');
plot(xc-x,yc-y,'r.');
plot(xc+y,yc+x,'b.');
plot(xc+y,yc-x,'k.');
plot(xc-y,yc+x,'g.');
plot(xc-y,yc-x,'b.');
da=(x+1)^2+y^2-r^2;
db=(x+1)^2+(y-1)^2-r^2;
i=i+1;
disp(' x y da x+1 y db x+1 y-1' )
trace(i,:)=[x y da x+1 y db x+1 y-1]
if abs(db)<abs(da)
y=y-1;
end
x=x+1;
end
{"html5":"htmlmixed","css":"css","javascript":"javascript","php":"php","python":"python","ruby":"ruby","lua":"text\/x-lua","bash":"text\/x-sh","go":"go","c":"text\/x-csrc","cpp":"text\/x-c++src","diff":"diff","latex":"stex","sql":"sql","xml":"xml","apl":"apl","asterisk":"asterisk","c_loadrunner":"text\/x-csrc","c_mac":"text\/x-csrc","coffeescript":"text\/x-coffeescript","csharp":"text\/x-csharp","d":"d","ecmascript":"javascript","erlang":"erlang","groovy":"text\/x-groovy","haskell":"text\/x-haskell","haxe":"text\/x-haxe","html4strict":"htmlmixed","java":"text\/x-java","java5":"text\/x-java","jquery":"javascript","mirc":"mirc","mysql":"sql","ocaml":"text\/x-ocaml","pascal":"text\/x-pascal","perl":"perl","perl6":"perl","plsql":"sql","properties":"text\/x-properties","q":"text\/x-q","scala":"scala","scheme":"text\/x-scheme","tcl":"text\/x-tcl","vb":"text\/x-vb","verilog":"text\/x-verilog","yaml":"text\/x-yaml","z80":"text\/x-z80"}