Get all perfect numbers on n

    n=10
    a=zeros(n);
    for i=1:n
        for j=1:(i/2)
        if (rem(i,j)==0)
                a(j,i)=j
        endif
    endfor
    endfor
    b=sum(a);
    perfectos=0;
    for z=1:n
        if (b(z)-z==z)
            perfectos=z
    endif
    endfor