for j=2:2000 
c=sum(diviseurs(j))-j
d=sum(diviseurs(c))-c
if (d==j) then disp ([j,c])
end
end