 = ( R + r \cos(\theta) )  \cos(\psi),\]</div>
<div class="math notranslate nohighlight">
\[y = ( R + r\cos(\theta) ) \sin(\psi),\]</div>
<div class="math notranslate nohighlight">
\[z = r \sin(\theta)\]</div>
<p>Need to draw <span class="math notranslate nohighlight">\(\theta\)</span> with distribution</p>
<div class="math notranslate nohighlight">
\[g(\theta) = (1+ r \cos(\theta)/R ) / (2\pi)\]</div>
<p>on <span class="math notranslate nohighlight">\(0 \leq \theta &lt; 2\pi\)</span>, and <span class="math