Add compe.m.
Compe computes E which is higly related with R. In short, R is the transpose of the rotation matrix from polar(E).