function SincF = sinc123(x,Type) % Voir Matlab #10
 
N=length(x);
 
if Type==1 % 1D
    SincF=sinc(x);
    
elseif Type==2 % 2D
    sinc_2D=zeros(N,N);
    y=x;
    for i=1:N
        for j=1:N
            sinc_x=sinc(x(i));
            sinc_y=sinc(y(j));
            sinc_2D(i,j)=sinc_x*sinc_y;
        end;
    end;
    SincF=sinc_2D;
else % 3D
    sinc_3D=zeros(N,N,N);
    y=x; z=x;
    for i=1:N
        for j=1:N
            for k=1:N
                sinc_x=sinc(x(i));
                sinc_y=sinc(y(j));
                sinc_z=sinc(z(k));
                sinc_3D(i,j,k)=sinc_x*sinc_y*sinc_z;
            end;
        end;
    end;
    SincF=sinc_3D;
end;
 
 
end