   m = 5;  ep = 1.e-10;
   q = randn(m,1);   q = q/norm(q);
   r1 = randn(m,1);
   r2 = randn(m,1);
   X = [q, q+ep*r1, q+ep*r2];
