Calculate polar radius from flattening

Only the semi-major axis and the flattening (or inverse flattening in this case) are defining constants, the semi-minor axis is derived from them. This is reflected in this change.
Merge pull request #2 from benelsen/master

Calculate polar radius from flattening