Should it be R=a(1-e^2)/(1-ecos(theta))?
Otherwise, great article.

CAPTCHA
This question is for testing whether you are a human visitor and to prevent automated spam submissions.