function [f]=factors(number)
    x=factor(number);
    u=(unique(x))';
    v=zeros(size(u,1),1);
    u=[u v];
    for i=1:size(u,1);
        u(i,2)=sum(u(i,1)==x);
    end
    f=u
endfunction