function [ell] = schm(u)
v = normalize(u);
s = size(v);
n = ss(1);

[U,D,V]=svd(v);
disp('U''');
disp(U');
disp('V');
disp(V);

for i = 1:n
	disp(i)
end

ell = diag(D);
disp('Schmidt coefficients');
disp(ell);
