How do you get the singular value decomposition of a matrix in Narray into a [u,x,v] triple of matrices? Then I want to grab the first column of u (u(:,1) in matlab). I am also interested in stacking a bunch of data matrices with various lags: c = column dimension of m stacked_m = [m(:, 1:c-2); m(:, 2:c-1); m(:,3:c)] Finally, just for grins, I want to assign 1 to all the elements of a matrix if which are less than 1: m(m<1) = 1 I think there is a *lot* of work, in various bits and pieces. (And I don't write Japanese, unfortunately.) However, it seems to be all over the place. (I was never very fond of the CPAN haphazard approach to "library" development.) I don't want to complain--I would be quite happy to help with a project. I just think we need a push to come up with a use-case based integration of all the tools, like scipy in Python (really like scipy + matplotlib), all in one coherent package.