look into eigen vectors, covariance matrices, to get the orientation of the ellipsoid. As an experiment, I used a hacked version of the Rapid Collision Detection code to find bounding boxes around sets of points. then I derived my own

another code sample

http://www.cs.unc.edu/~hoff/research/vfculler/vfculler/obbtree.cpp