Findn(triu()) in Julia

[digital_carver]

For future readers:

Findn() function not found in Julia

It looks like findn(A) was a function in quite old versions of Julia, which was deprecated in Julia 0.7 in favor of findall(!iszero, A). So, you should update your code to aij = findall(!iszero, triu(NA,1)) where aij is now an array of CartesianIndex, essentially a kind of multi-component index for indexing multidimensional arrays.