Facebook
From gg, 1 Month ago, written in Plain Text.
Embed
Download Paste or View Raw
Hits: 145
  1. clc; clear;
  2. xc=input('Xc= ');yc=input('yc= ');
  3. r=input('r= ');
  4. hold on;
  5. i=0;
  6. trace=zeros(1,8);
  7. axis([-100 100 -100 100])
  8. x=0; y=r;
  9. while x<=y
  10.  plot(xc+x,yc+y,'g.');
  11.  plot(xc+x,yc-y,'b.');
  12.  plot(xc-x,yc+y,'k.');
  13.  plot(xc-x,yc-y,'r.');
  14.  plot(xc+y,yc+x,'b.');
  15.  plot(xc+y,yc-x,'k.');
  16.  plot(xc-y,yc+x,'g.');
  17.  plot(xc-y,yc-x,'b.');
  18.   da=(x+1)^2+y^2-r^2;
  19.  db=(x+1)^2+(y-1)^2-r^2;
  20.   i=i+1;
  21.   disp('     x    y    da    x+1     y    db x+1    y-1'    )
  22.   trace(i,:)=[x y da x+1 y db x+1 y-1]
  23.  if abs(db)<abs(da)
  24.     y=y-1;
  25.  end  
  26.  x=x+1;
  27. end
  28.