function B = disque(Rayon,pas);
%   B=disque(R,pas) 
%
%Génère la fonction caractéristique d'un disque de rayon R, avec un pas
%donné (pas). 

x = -Rayon:pas:Rayon-pas;
L = length(x);
S = x'*x;
f = ones(L);
for i = 1:L 
   for j = 1:L 
      a = x(i)^2+ x(j)^2; 
      if a > (Rayon)^2 
         f(i,j) = 0; 
      end;
   end;
end;
B = f; 