Matlab - consulta

function  [n,factores,factores1,ndiv]= divisores1(n)

if n==1

    factores=1;

    factores1=1;

    ndiv=1;

else

    factores=factor(n);

    factores1=unique(factores);

    ndiv=1;

    for i=1:length(factores1)

        cantidad=length(find(factores1(i)==factores));

        ndiv=ndiv*(cantidad+1);

    end

end

>> [n,factores,factores1,ndiv]= divisores1(15)

n =

    15

factores =

     3     5

factores1 =

     3     5

ndiv =

     4

>> [n,factores,factores1,ndiv]= divisores1(20)

n =

    20

factores =

     2     2     5

factores1 =

     2     5

ndiv =

     6