exec('init_tp_image.sce'); exec('contours.sci');
n = 200;
// image de taille 200x200 // avec un fond noir
   im = zeros(n,n);
   for i = 1:n
        for j = 1:n
              I = i - n/2;
              J = j - n/2;
              r = sqrt(I^2 + J^2);
              if r < n/4 then
            // le disque au centre contient un degrade
                im(i,j) = int(255*(j-n/5)/n*1.75);
              end
        end
    end
// on lisse l'image pour avoir un bord moins abrupt
G = W_gauss_2D(1.5);
im = int(conv2(im, G, "same"));


afficher_image(im);