{"id":1585,"date":"2005-02-19T18:32:25","date_gmt":"2005-02-20T02:32:25","guid":{"rendered":"http:\/\/www.ogre.nu\/wp\/?p=1585"},"modified":"2010-05-26T15:12:43","modified_gmt":"2010-05-26T23:12:43","slug":"the-stubborn-spiral","status":"publish","type":"post","link":"https:\/\/bendwavy.org\/wp\/?p=1585","title":{"rendered":"the stubborn spiral"},"content":{"rendered":"<table align=center>\n<tr>\n<td> Rusin&#8217;s disco ball <br \/> <img decoding=\"async\" src=\"\/pack\/480Rusin19.png\"\/> <\/td>\n<td> golden angles <br \/> <img decoding=\"async\" src=\"\/pack\/480gold.png\"\/> <\/td>\n<td> Saff &#038; Kuijlaars <br \/> <img decoding=\"async\" src=\"\/pack\/480SK.png\"\/> <\/td>\n<\/tr>\n<\/table>\n<p>Examples of <a href=\"\/pack\/pack.htm\">three algorithms<\/a> for distributing nodes fairly evenly over a sphere.  Those on the middle and right slice the sphere into parallel bands of equal area (much narrower than the white discs), and put one node (center of a disc) somewhere in each band.  Saff &#038; Kuijlaars place the nodes along a spiral path across the bands, keeping the distance between turns of the spiral roughly constant.  Failing to grok how their rule does that, I approach it from another angle.<br \/>\n<!--more--><br \/>\nIn the ideal limit (N approaching infinity), the unit sphere is covered exactly once by a strip of width 2&radic;(&pi;\/N) and length 2&radic;(N&pi;); the strip is divided into N nearly square pieces, with a node in the center of each.  What parametric functions of -1 &lt; t &lt; 1 describe the midline of the strip?<\/p>\n<p>(In what follows, &theta; and &phi; have their conventional meanings of longitude and co-latitude, i.e. angular distance from a pole.)<\/p>\n<p>The equal slices rule gives cos&phi; = t.  I seek a function &theta;(t) such that the derivative along the path has constant magnitude:<\/p>\n<blockquote><p>\nds&sup2; = d&phi;&sup2; + sin&sup2;&phi; d&theta;&sup2; = N&pi; dt&sup2;\n<\/p><\/blockquote>\n<p>I proceed:<\/p>\n<blockquote><p>\n&ndash;sin&phi; d&phi; = dt; d&phi; = &ndash;dt \/ sin&phi; = &ndash;dt \/ &radic;(1&ndash;t&sup2;)<\/p>\n<p>ds&sup2; = N&pi; dt&sup2; = ( dt&sup2; \/ (1&ndash;t&sup2;) ) + (1&ndash;t&sup2;) d&theta;&sup2;<\/p>\n<p>(1&ndash;t&sup2;) d&theta;&sup2; = dt&sup2; (N&pi; &ndash; 1 \/ (1&ndash;t&sup2;))<br \/>\n&nbsp; = dt&sup2; (N&pi; &ndash; N&pi;t&sup2; &ndash; 1) \/ (1&ndash;t&sup2;)<\/p>\n<p>d&theta;&sup2; = dt&sup2; (N&pi; &ndash; N&pi;t&sup2; &ndash; 1) \/ (1&ndash;t&sup2;)&sup2;\n<\/p><\/blockquote>\n<p>So I ask <a href=\"http:\/\/integrals.wolfram.com\/\">the Wolfram Integrator<\/a> to integrate<\/p>\n<blockquote><p><code> Sqrt[ n*Pi - n*Pi*x^2 - 1 ] \/ (1-x^2) <\/code><\/p><\/blockquote>\n<p> and it says <center><img decoding=\"async\" align=center src=\"\/doodle\/integr1.gif\"\/><\/center><\/p>\n<p>Now, I do make mistakes; I&#8217;ve worked the problem several times now and got subtly different answers.  Last time around, in fact, there were no imaginary coefficients; so I made this in <a href=\"http:\/\/www.povray.org\/\">PoV-Ray<\/a>:<br \/>\n<center><img decoding=\"async\" src=\"\/doodle\/badhelix.png\"\/><\/center><\/p>\n<p>Not quite what I&#8217;m looking for.  Can you spot where I went wrong?<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Rusin&#8217;s disco ball golden angles Saff &#038; Kuijlaars Examples of three algorithms for distributing nodes fairly evenly over a sphere. Those on the middle and right slice the sphere into parallel bands of equal area (much narrower than the white &hellip; <a href=\"https:\/\/bendwavy.org\/wp\/?p=1585\">Continue reading <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[12],"tags":[],"class_list":["post-1585","post","type-post","status-publish","format-standard","hentry","category-mathematics"],"_links":{"self":[{"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=\/wp\/v2\/posts\/1585","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=1585"}],"version-history":[{"count":1,"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=\/wp\/v2\/posts\/1585\/revisions"}],"predecessor-version":[{"id":2531,"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=\/wp\/v2\/posts\/1585\/revisions\/2531"}],"wp:attachment":[{"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1585"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1585"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/bendwavy.org\/wp\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1585"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}