<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div class="">Thanks for all the help from the group and saving the animation works well.</div><div class="">However [had somebody noticed?], there was an error in the statistical formulation to generate a uniform random point distribution within the sphere. The more correct (but perhaps clumsy) version is below.</div><div class="">Heinz</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><pre style="font-family: Monospaced;" class=""><font size="4" class=""><span style="color:rgb(0,0,0);" class="">R</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(188,143,143);" class="">23</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">N</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(188,143,143);" class="">14500</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">nmax</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(188,143,143);" class="">1.1</span><span style="color:rgb(92,92,92);" class="">*</span><span style="color:rgb(188,143,143);" class="">6</span><span style="color:rgb(92,92,92);" class="">*</span><span style="color:rgb(0,0,0);" class="">N</span><span style="color:rgb(92,92,92);" class="">/</span><span style="color:rgb(218,112,214);" class="">%pi</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(100,174,100);font-style:italic;" class="">// Random points in cube 2R*2R*2R</span>
<span style="color:rgb(0,0,0);" class="">x</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(50,185,185);" class="">grand</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(0,0,0);" class="">nmax</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">1</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">'</span><span style="color:rgb(188,143,143);" class="">unf</span><span style="color:rgb(188,143,143);" class="">'</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(92,92,92);" class="">-</span><span style="color:rgb(0,0,0);" class="">R</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(0,0,0);" class="">R</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">y</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(50,185,185);" class="">grand</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(0,0,0);" class="">nmax</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">1</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">'</span><span style="color:rgb(188,143,143);" class="">unf</span><span style="color:rgb(188,143,143);" class="">'</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(92,92,92);" class="">-</span><span style="color:rgb(0,0,0);" class="">R</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(0,0,0);" class="">R</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">z</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(50,185,185);" class="">grand</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(0,0,0);" class="">nmax</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">1</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">'</span><span style="color:rgb(188,143,143);" class="">unf</span><span style="color:rgb(188,143,143);" class="">'</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(92,92,92);" class="">-</span><span style="color:rgb(0,0,0);" class="">R</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(0,0,0);" class="">R</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">R2</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(0,0,0);" class="">R</span><span style="color:rgb(92,92,92);" class="">*</span><span style="color:rgb(0,0,0);" class="">R</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(100,174,100);font-style:italic;" class="">// Exclude points outside of sphere with radius R</span>
<span style="color:rgb(0,0,0);" class="">T</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(0,0,0);" class="">x</span><span style="color:rgb(92,92,92);" class="">^</span><span style="color:rgb(188,143,143);" class="">2</span><span style="color:rgb(92,92,92);" class="">+</span><span style="color:rgb(0,0,0);" class="">y</span><span style="color:rgb(92,92,92);" class="">^</span><span style="color:rgb(188,143,143);" class="">2</span><span style="color:rgb(92,92,92);" class="">+</span><span style="color:rgb(0,0,0);" class="">z</span><span style="color:rgb(92,92,92);" class="">^</span><span style="color:rgb(188,143,143);" class="">2</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(92,92,92);" class="">></span><span style="color:rgb(0,0,0);" class="">R2</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">M</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(74,85,219);" class="">[</span><span style="color:rgb(0,0,0);" class="">x</span> <span style="color:rgb(0,0,0);" class="">y</span> <span style="color:rgb(0,0,0);" class="">z</span><span style="color:rgb(74,85,219);" class="">]</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">M</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(0,0,0);" class="">T</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(255,170,0);" class="">:</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(74,85,219);" class="">[</span><span style="color:rgb(74,85,219);" class="">]</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(74,85,219);" class="">[</span><span style="color:rgb(0,0,0);" class="">r</span> <span style="color:rgb(0,0,0);" class="">c</span><span style="color:rgb(74,85,219);" class="">]</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(50,185,185);" class="">size</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(0,0,0);" class="">M</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">M</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(0,0,0);" class="">N</span><span style="color:rgb(92,92,92);" class="">+</span><span style="color:rgb(188,143,143);" class="">1</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(255,170,0);" class="">:</span><span style="color:rgb(0,0,0);" class="">r</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">,</span> <span style="color:rgb(255,170,0);" class="">:</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(74,85,219);" class="">[</span><span style="color:rgb(74,85,219);" class="">]</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">x</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(0,0,0);" class="">M</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(255,170,0);" class="">:</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">1</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span> <span style="color:rgb(0,0,0);" class="">y</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(0,0,0);" class="">M</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(255,170,0);" class="">:</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">2</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span> <span style="color:rgb(0,0,0);" class="">z</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(0,0,0);" class="">M</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(255,170,0);" class="">:</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">3</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">f</span> <span style="color:rgb(92,92,92);" class="">=</span> <span style="color:rgb(50,185,185);" class="">figure</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">f</span><span style="color:rgb(92,92,92);" class="">.</span><span style="color:rgb(170,170,170);" class="">background</span> <span style="color:rgb(92,92,92);" class="">=</span> <span style="color:rgb(188,143,143);" class="">8</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">s</span> <span style="color:rgb(92,92,92);" class="">=</span> <span style="color:rgb(50,185,185);" class="">linspace</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(188,143,143);" class="">5</span><span style="color:rgb(0,0,0);" class="">,</span> <span style="color:rgb(188,143,143);" class="">1</span><span style="color:rgb(0,0,0);" class="">,</span> <span style="color:rgb(0,0,0);" class="">N</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline;" class="">scatter3d</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(0,0,0);" class="">x</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(0,0,0);" class="">y</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(0,0,0);" class="">z</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(0,0,0);" class="">msizes</span><span style="color:rgb(92,92,92);" class="">=</span><span style="color:rgb(188,143,143);" class="">1</span><span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(0,0,0);" class="">s</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">a</span> <span style="color:rgb(92,92,92);" class="">=</span> <span style="color:rgb(174,92,176);text-decoration:underline;" class="">gca</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">outgif</span> <span style="color:rgb(92,92,92);" class="">=</span> <span style="color:rgb(188,143,143);" class="">"</span><span style="color:rgb(188,143,143);" class="">Sphere01.gif</span><span style="color:rgb(188,143,143);" class="">"</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">idGif</span> <span style="color:rgb(92,92,92);" class="">=</span> <span style="color:rgb(174,92,176);text-decoration:underline;" class="">animaGIF</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(174,92,176);text-decoration:underline;" class="">gcf</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">,</span> <span style="color:rgb(0,0,0);" class="">outgif</span><span style="color:rgb(0,0,0);" class="">,</span> <span style="color:rgb(0,0,0);" class="">,</span><span style="color:rgb(188,143,143);" class="">2</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(160,32,240);" class="">for</span> <span style="color:rgb(0,0,0);" class="">i</span> <span style="color:rgb(92,92,92);" class="">=</span> <span style="color:rgb(188,143,143);" class="">20</span><span style="color:rgb(255,170,0);" class="">:</span><span style="color:rgb(188,143,143);" class="">60</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">a</span><span style="color:rgb(92,92,92);" class="">.</span><span style="color:rgb(170,170,170);" class="">rotation_angles</span> <span style="color:rgb(92,92,92);" class="">=</span> <span style="color:rgb(74,85,219);" class="">[</span><span style="color:rgb(188,143,143);" class="">45</span><span style="color:rgb(0,0,0);" class="">,</span> <span style="color:rgb(0,0,0);" class="">i</span><span style="color:rgb(74,85,219);" class="">]</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(0,0,0);" class="">idGif</span> <span style="color:rgb(92,92,92);" class="">=</span> <span style="color:rgb(174,92,176);text-decoration:underline;" class="">animaGIF</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(174,92,176);text-decoration:underline;" class="">gcf</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">,</span> <span style="color:rgb(0,0,0);" class="">idGif</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span> <span style="color:rgb(100,174,100);font-style:italic;" class="">// Adds the current figure to the GIF stream</span>
<span style="color:rgb(160,32,240);" class="">end</span><span style="color:rgb(0,0,0);" class="">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline;" class="">animaGIF</span><span style="color:rgb(74,85,219);" class="">(</span><span style="color:rgb(0,0,0);" class="">idGif</span><span style="color:rgb(74,85,219);" class="">)</span><span style="color:rgb(0,0,0);" class="">;</span> <span style="color:rgb(100,174,100);font-style:italic;" class="">// Closes the GIF stream</span></font></pre><div class=""><br class=""></div></div><div class=""><br class=""></div><br class=""><br class=""><blockquote type="cite" class="">On 02.07.2021, at 00:15, Heinz Nabielek <<a href="mailto:heinznabielek@me.com" class="">heinznabielek@me.com</a>> wrote:<br class=""><br class="">Colleagues:<br class=""><br class="">the code below generates 14,500 random points in a spherical volume.<br class="">How can I make a movie that shows how the sphere rotates around the z-axis? And how do I store the animation? animaGIF does not exist on my Scilab 6.1 in my macOS 10.15.7<br class=""><br class="">Heinz<br class=""><br class=""><br class="">R=23;<br class="">N=14500;<br class="">r=grand(N,1,'unf',0,R);<br class="">phi=grand(N,1,'unf',0,2*%pi);<br class="">theta=grand(N,1,'unf',0,%pi);<br class="">x=r.*cos(phi).*sin(theta);<br class="">y=r.*sin(phi).*sin(theta);<br class="">z=r.*cos(theta);<br class="">scatter3d(x,y,z,msizes=4);<br class="">gca().rotation_angles = [60, 60];<br class=""></blockquote><br class=""></body></html>