function [MU,R,B,T] = BoundingBox(D)
   [M,N] = size(D);
   MU=mean(D);
   for I=1:M 
      Q(I,:)=D(I,:)-MU;
   end
   [L,S,R] = svd(Q);
   for J=1:N
     P=Q*R(:,J);
     B(J) = min(P);
     T(J) = max(P);
   end
end

