somcub(n)

               function x = somcub(n)
               m = (1:n).^3;                  % vecteur ligne des cubes
               m = m(ones(1,n),:);            % rptition des lignes
               m = triu(m+m');                % matrice triangulaire i^3+j^3
               m = sort(m(:)');               % tris des lments de m
               m(m==0)=[];                    % suppression des zros
               x = m(logical([0,~diff(m)]));  % deux conscutifs gaux -> rsultat