From,
http://en.wikipedia.org/wiki/Geodetic_system
The transformation matrix from ENU-coordinates to ECEF is:
Figure 1 - Transformation Matrix
Working with this, I have found that when (x,y,z) and (Xr, Yr, Zr) represent positions vectors, the transformation works fine.
However...