function franges_circ (r,I,rgb)
% r, I sont des vecteurs de meme longueur n
% rgb doit avoir dimensions (n,3) ou (1,n,3)

a = linspace(0, 2*pi, 361);
m = 0;

for n = 1: length(r)
    for i = 1:length(a)
        z(i,n) = I(n)+0.5;
        [x(i,n), y(i,n)] = pol2cart (a(i), r(n));
        for k = 1:3
            if (ndims(rgb) == 3)
                C(i,n,k) = rgb(1,n,k);
            else
                C(i,n,k) = rgb(n,k);
            end
        end
    end
end
figure; surf (x,y,z,C)
xlabel('x (m)'); ylabel('y (m)')
title ('intensité des franges')
view(2)
shading interp
axis('square')