// Factoriser un entier donné n>=2.
// Chager la fonction Erato.sci.
function F=factorisation(n)
Facteurs=Erato(n);
Puissances=[];
for i=1:taille(Facteurs)
    j=1;
    while reste(n,Facteurs(i)^j)==0
        j=j+1
    end
    Puissances=[Puissances,j-1];
end
K=find(Puissances==0);
Facteurs(K)=[];
Puissances(K)=[];
F=[Facteurs;Puissances];
endfunction