julia cheat sheet... svd - singular value decomposition a = readdlm("/tmp/analysis") b = ones(size(a)[1])' * a c = ones(size(a)[1]) * b / size(a)[1] d = a - c