function x = homosolve(u)

utu = u'*u;
[evec ev] = eig(utu);
[minev minidx] = min(diag(ev));
x = evec(:,minidx);
