<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On 26 September 2013 10:02, Daπid <span dir="ltr"><<a href="mailto:davidmenhur@gmail.com" target="_blank">davidmenhur@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">The simplest way is to do it in cartesian coordinates: take x, y, and z independently from N(0,1). If you want to generate only one normal number per step, consider the jacobian in the angles.</blockquote>

</div><br><br></div><div class="gmail_extra">Actually, this is wrong, as it would allow displacements (at 1 sigma) of 1 along the axis, but up to sqrt(3) along diagonals. What you actually want is a multivariate normal distribution with covariance proportional to the identity (uncorrelation between axis and isotropy).<br>

<br></div><div class="gmail_extra">See formulae here: <a href="http://en.wikipedia.org/wiki/Multivariate_normal_distribution">http://en.wikipedia.org/wiki/Multivariate_normal_distribution</a><br></div></div>