[continued from previous message]   
      
   This make it much easier to program our dynamical model, because we can   
   now easily check for (a) by having one of the components of G be dL/dt,   
   and we can check for (b) by having one of the components of G be the   
   distance between an adjacent pair of particles.  And, with this type of   
   programming, we can return to specifying only the goal times when we want   
   to print the solution, and let the ODE integration routine automatically   
   choose efficient internal step sizes.   
      
   So, at the cost of having to use a somewhat more complicated set of   
   ODE integration software, we get both easier programming *and* a much   
   more efficient program.   
   <>   
      
   If I find the time, I'll write a program implementing the above and   
   post some of the output.   
      
   ciao,   
   --   
   -- "Jonathan Thornburg [remove -color to reply]"    
      (he/him; currently on the west coast of Canada)   
      "The ideal subject of totalitarian rule is not the convinced Nazi or   
       the convinced Communist, but people for whom the distinction between   
       fact and fiction (i.e., the reality of experience) and the distinction   
       between true and false (i.e., the standards of thought) no longer exist."   
                                -- Hannah Arendt, The Origins of Totalitarianism   
      
   --- SoupGate-DOS v1.05   
    * Origin: you cannot sedate... all the things you hate (1:229/2)