FORM 4.0 (Jun 30 2012) 64-bits                   Run: Sat Aug  3 17:45:52 2013
    *   Title:  SMALL CHARGED SPHERE:
    *           RdotU comparison to Schoonschip
    *    File:  rdotu-ck.frm
    *  Author:  David N. Williams
    * License:  Creative Commons Attribution-Share Alike
    * Started:  July 8, 2012
    * Revised:  July 8-10,28, 2012
    *
    * All parts of this program not in the public domain are:
    *
    *   Copyright (C) 1987, 2012 David N. Williams
    *
    * This work is licensed under the Creative Commons Attribution-
    * Share Alike 2.5 License.  To view a copy of this license, visit
    * http://creativecommons.org/licenses/by-sa/2.5/ or send a letter
    * to Creative Commons, 543 Howard Street, 5th Floor, San
    * Francisco, California, 94105, USA.
    *
    *   Input:  RdotU from rdotu.sav
    *  Output:  none
    *
    * This program loads the output rdotu.sav of rdotu.frm and
    * checks it against Schoonschip results pasted from rdotu.out.
    
    *** DECLARATIONS
    
    off statistics;
    
    S eps;
    S c;
    S snorm;
    I mu,nu;
    V u,a,j,n;
    T d1y,d2y,d3y,d1u,d2u,d1a;
    Load rdotu.sav;
 RdotU loaded
    
    S norm, aDa, p;
    I Naa, Nab;
    
    *** MODULES
    
    * NeXT Schoonschip output, May 12, 1993
    * D1..., D2..., D3... replaced by d1..., d2..., d3... 28-Jul-12
    
    L RdotUnext =
     + eps
      * ( c*norm )
     
     + eps^3
      * ( - 1/8*c^-3*norm^3*aDa )
     
     + d1y(Naa,n)*d1y(Nab,n)*d1y(u,n)*d1y(u,n)*d1u(Naa,n)*d1u(Nab,n)*eps^3
      * ( - 1/2*c^-3*norm^-3 )
     
     + d1y(Naa,n)*d1y(Nab,n)*d1y(u,n)*d2y(Naa,n,n)*d1u(Nab,n)*eps^3
      * ( 1/4*c^-1*norm^-3 )
     
     + d1y(Naa,n)*d1y(Nab,n)*d1y(u,n)*d2y(Nab,n,n)*d1u(Naa,n)*eps^3
      * ( 1/4*c^-1*norm^-3 )
     
     + d1y(Naa,n)*d1y(Nab,n)*d2y(Naa,n,n)*d2y(Nab,n,n)*eps^3
      * ( - 1/8*c*norm^-3 )
     
     + d1y(Naa,n)*d1y(Nab,n)*d2y(Naa,n,n)*d1u(Nab,n)*eps^3
      * ( 1/4*norm^-2 )
     
     + d1y(Naa,n)*d1y(Nab,n)*d2y(Nab,n,n)*d1u(Naa,n)*eps^3
      * ( - 1/4*norm^-2 )
     
     + d1y(Naa,n)*d1y(Nab,n)*d1u(Naa,n)*d1u(Nab,n)*eps^3
      * ( 1/2*c^-1*norm^-1 )
     
     + d1y(Naa,n)*d1y(u,n)*d1y(u,n)*d1y(u,n)*d1y(a,n)*d1u(Naa,n)*eps^3
      * ( 1/2*c^-5*norm^-3 )
     
     + d1y(Naa,n)*d1y(u,n)*d1y(u,n)*d1y(a,n)*d2y(Naa,n,n)*eps^3
      * ( - 1/4*c^-3*norm^-3 )
     
     + d1y(Naa,n)*d1y(u,n)*d1y(u,n)*d2y(u,n,n)*d1u(Naa,n)*eps^3
      * ( - 1/2*c^-3*norm^-3 )
     
     + d1y(Naa,n)*d1y(u,n)*d1y(u,n)*d1a(Naa,n)*eps^3
      * ( - 1/2*c^-3*norm^-1 )
     
     + d1y(Naa,n)*d1y(u,n)*d1y(a,n)*d1u(Naa,n)*eps^3
      * ( - 1/2*c^-3*norm^-1 )
     
     + d1y(Naa,n)*d1y(u,n)*d2y(Naa,n,n)*d2y(u,n,n)*eps^3
      * ( 1/4*c^-1*norm^-3 )
     
     + d1y(Naa,n)*d1y(u,n)*d1u(Naa,n)*eps^2
      * ( c^-1*norm^-1 )
     
     + d1y(Naa,n)*d1y(u,n)*d2u(Naa,n,n)*eps^3
      * ( 1/2*c^-1*norm^-1 )
     
     + d1y(Naa,n)*d1y(a,n)*d2y(Naa,n,n)*eps^3
      * ( - 1/4*c^-1*norm^-1 )
     
     + d1y(Naa,n)*d2y(Naa,n,n)*eps^2
      * ( - 1/2*c*norm^-1 )
     
     + d1y(Naa,n)*d2y(u,n,n)*d1u(Naa,n)*eps^3
      * ( 1/2*c^-1*norm^-1 )
     
     + d1y(Naa,n)*d3y(Naa,n,n,n)*eps^3
      * ( - 1/6*c*norm^-1 )
     
     + d1y(Naa,n)*d1a(Naa,n)*eps^3
      * ( 1/2*c^-1*norm )
     
     + d1y(u,n)*eps^3
      * ( - 1/3*c^-4*norm^2*aDa )
     
     + d1y(u,n)*d1y(u,n)*eps^3
      * ( - 1/4*c^-5*norm*aDa )
     
     + d1y(u,n)*d1y(u,n)*d1y(u,n)*d1y(u,n)*eps^3
      * ( 1/24*c^-7*norm^-1*aDa )
     
     + d1y(u,n)*d1y(u,n)*d1y(u,n)*d1y(u,n)*d1y(a,n)*d1y(a,n)*eps^3
      * ( - 1/8*c^-7*norm^-3 )
     
     + d1y(u,n)*d1y(u,n)*d1y(u,n)*d1y(a,n)*d2y(u,n,n)*eps^3
      * ( 1/4*c^-5*norm^-3 )
     
     + d1y(u,n)*d1y(u,n)*d1y(u,n)*d1y(j,n)*eps^3
      * ( 1/6*c^-5*norm^-1 )
     
     + d1y(u,n)*d1y(u,n)*d1y(a,n)*eps^2
      * ( - 1/2*c^-3*norm^-1 )
     
     + d1y(u,n)*d1y(u,n)*d1y(a,n)*d1y(a,n)*eps^3
      * ( 1/4*c^-5*norm^-1 )
     
     + d1y(u,n)*d1y(u,n)*d2y(u,n,n)*d2y(u,n,n)*eps^3
      * ( - 1/8*c^-3*norm^-3 )
     
     + d1y(u,n)*d1y(u,n)*d2y(a,n,n)*eps^3
      * ( - 1/4*c^-3*norm^-1 )
     
     + d1y(u,n)*d1y(a,n)*d2y(u,n,n)*eps^3
      * ( - 1/4*c^-3*norm^-1 )
     
     + d1y(u,n)*d1y(j,n)*eps^3
      * ( - 1/2*c^-3*norm )
     
     + d1y(u,n)*d2y(Naa,n,n)*d1u(Naa,n)*eps^3
      * ( 1/2*c^-1*norm^-1 )
     
     + d1y(u,n)*d2y(u,n,n)*eps^2
      * ( 1/2*c^-1*norm^-1 )
     
     + d1y(u,n)*d3y(u,n,n,n)*eps^3
      * ( 1/6*c^-1*norm^-1 )
     
     + d1y(a,n)*eps^2
      * ( 1/2*c^-1*norm )
     
     + d1y(a,n)*d1y(a,n)*eps^3
      * ( - 1/8*c^-3*norm )
     
     + d1y(j,n)*eps^3
      * ( - 1/3*c^-2*norm^2 )
     
     + d2y(Naa,n,n)*d2y(Naa,n,n)*eps^3
      * ( - 1/8*c*norm^-1 )
     
     + d2y(u,n,n)*d2y(u,n,n)*eps^3
      * ( 1/8*c^-1*norm^-1 )
     
     + d2y(a,n,n)*eps^3
      * ( 1/4*c^-1*norm ) + 0;
    
    id aDa = a.a;
    id norm = snorm;
    id norm^p? = snorm^p;
    sum Naa,Nab;
    .sort
    
    hide RdotUnext;
    
    G diff = RdotU - RdotUnext;
    b eps,snorm;
    print;
    .end

   diff = 0;

  0.00 sec out of 0.00 sec