ICA algorithm

[W,s,v]  = svd((repmat(sum(x.*x,1),szie(x,1),1).*x)*x');