function [Z,S]=scaled(X)
[U,S,V]=svd(X);
Z=X*pinv(S);
end