circlegen.m

有趣的可视的数值方法 出自网站http://www.mathworks.com/moler

function circlegen(h)
%CIRCLEGEN  Generate approximate circles.
%   CIRCLEGEN(h) uses step size h.
%   CIRCLE with no arguments uses h = 0.20906.

if nargin < 1
   h = sqrt(2*(1-cos(2*pi/30)));
end

shg
clf
set(gcf,'doublebuffer','on')
z = exp(2*pi*i*(0:256)/256);
grey = [.8 .8 .8];
p = plot(real(z),imag(z),'-',0,0,'o',1,0,'.');
set(p(1:2),'color',grey);
p = p(3);
set(p,'color','black','erasemode','none')
axis([-2 2 -2 2])
axis square
title(['h = ' num2str(h)])
stop = uicontrol('style','toggle','string','stop');

x = 1;
y = 0;
while ~get(stop,'value')
   x = x + h*y;
   y = y - h*x;
   set(p,'xdata',x,'ydata',y);
   drawnow
end

set(stop,'string','close','value',0,'callback','close')