Números amigos en un rango octave o Matlab

% |----> Este script me da los números amigos entre 1 y n <----|

function z = amigo(n);
tic
z = [ ];
n1 = [1 : n];

for i = 1:n
v(i) = sum( divisores( i ) ) ;
end

for i=1:n
for j=1:n

if ( n1( i ) == v( j ) && n1( j ) == v( i ) )
z=[z ; n1( i ) n1( j )] ;
end

end
end

toc
end

%Para desarrollar este algoritmo debemos tener guardado en la misma carpeta el siguiente algoritmo: