function cercle(a,b,cote,n)
%CERCLE(A,B,COTE,N) représente la parité des cercles

% On prend X entre A et A+COTE et Y entre B et B+COTE; COTE indique
% donc le pas. On calcule X^2+Y^2, puis sa partie entière et enfin
% on prend le reste de la division par N.
% Lire Pour la Science, novembre 1986.
% Exemple: prendre CERCLE(-15,-20,87,7)

[X,Y]=meshdom(a:cote/400:a+cote,b:cote/400:b+cote);
Z=fix(X.^2+Y.^2);

colormap(hot(n));
image(rem(Z,n));
axis equal
axis off

% Gadget: faites évoluer la palette des couleurs de façon cyclique.
%figure(gcf)
%spinmap(inf,1)
