I took the equation from: http://www.aanda.org/articles/aa/full/2002/45/aah3860/aah3860.right.html#eq:CROTA2PC