
Codemist Standard Lisp 3.54 for DEC Alpha: May 23 1994
Dump file created: Mon May 23 10:39:11 1994
REDUCE 3.5, 15-Oct-93 ...
Memory allocation: 6023424 bytes

+++ About to read file tstlib.red


n := 4;


n := 4


on rational, rat;


off allfac;



array p(n/2+2);



harmonic u,v,w,x,y,z;



weight e=1, b=1, d=1, a=1;



%% Step1: Solve Kepler equation
bige := fourier 0;


bige := 0

for k:=1:n do <<
  wtlevel k;
  bige:=fourier e * hsub(fourier(sin u), u, u, bige, k);
>>;


write "Kepler Eqn solution:", bige$


                             1   4                3   3
Kepler Eqn solution: - [( - ---*e )sin[4u] + ( - ---*e )sin[3u] + (
                             3                    8

 1   4    1   2             1   3
---*e  - ---*e )sin[2u] + (---*e  - e)sin[u]]
 6        2                 8


%% Ensure we do not calculate things of too high an order
wtlevel n;



%% Step 2: Calculate r/a in terms of e and l
dd:=-e*e;


          2
dd :=  - e
 hh:=3/2;


       3
hh := ---
       2
 j:=1;


j := 1
 cc := 1;


cc := 1

for i:=1:n/2 do <<
  j:=i*j; hh:=hh-1; cc:=cc+hh*(dd^i)/j
>>;



bb:=hsub(fourier(1-e*cos u), u, u, bige, n);


            1   4                3   3             1   4    1   2
bb := [( - ---*e )cos[4u] + ( - ---*e )cos[3u] + (---*e  - ---*e )cos
            3                    8                 3        2

               3   3                1   2
      [2u] + (---*e  - e)cos[u] + (---*e  + 1)]
               8                    2

aa:=fourier 1+hdiff(bige,u);


         4   4             9   3                1   4    2
aa := [(---*e )cos[4u] + (---*e )cos[3u] + ( - ---*e  + e )cos[2u] + 
         3                 8                    3

           1   3
      ( - ---*e  + e)cos[u] + 1]
           8
 ff:=hint(aa*aa*fourier cc,u);


               103   4                13   3             11   4
ff :=  - [( - -----*e )sin[4u] + ( - ----*e )sin[3u] + (----*e
               96                     12                 24

          5   2             1   3                  1   4
       - ---*e )sin[2u] + (---*e  - 2*e)sin[u] + (---*e  - 1)]
          4                 4                      8



%% Step 3: a/r and f
uu := hsub(bb,u,v);


            1   4                3   3             1   4    1   2
uu := [( - ---*e )cos[4v] + ( - ---*e )cos[3v] + (---*e  - ---*e )cos
            3                    8                 3        2

               3   3                1   2
      [2v] + (---*e  - e)cos[v] + (---*e  + 1)]
               8                    2
 uu:=hsub(uu,e,b);


            1   4                3   3             1   4    1   2
uu := [( - ---*b )cos[4v] + ( - ---*b )cos[3v] + (---*b  - ---*b )cos
            3                    8                 3        2

               3   3                1   2
      [2v] + (---*b  - b)cos[v] + (---*b  + 1)]
               8                    2

vv := hsub(aa,u,v);


         4   4             9   3                1   4    2
vv := [(---*e )cos[4v] + (---*e )cos[3v] + ( - ---*e  + e )cos[2v] + 
         3                 8                    3

           1   3
      ( - ---*e  + e)cos[v] + 1]
           8
 vv:=hsub(vv,e,b);


         4   4             9   3                1   4    2
vv := [(---*b )cos[4v] + (---*b )cos[3v] + ( - ---*b  + b )cos[2v] + 
         3                 8                    3

           1   3
      ( - ---*b  + b)cos[v] + 1]
           8

ww := hsub(ff,u,v);


               103   4                13   3             11   4
ww :=  - [( - -----*e )sin[4v] + ( - ----*e )sin[3v] + (----*e
               96                     12                 24

          5   2             1   3                  1   4
       - ---*e )sin[2v] + (---*e  - 2*e)sin[v] + (---*e  - 1)]
          4                 4                      8
 ww:=hsub(ww,e,b);


               103   4                13   3             11   4
ww :=  - [( - -----*b )sin[4v] + ( - ----*b )sin[3v] + (----*b
               96                     12                 24

          5   2             1   3                  1   4
       - ---*b )sin[2v] + (---*b  - 2*b)sin[v] + (---*b  - 1)]
          4                 4                      8


%% Step 4: Substitute f and f' into S
yy:=ff-ww;


         103   4             13   3                11   4    5   2
yy := [(-----*e )sin[4u] + (----*e )sin[3u] + ( - ----*e  + ---*e )
         96                  12                    24        4

                     1   3                     103   4
      sin[2u] + ( - ---*e  + 2*e)sin[u] + ( - -----*b )sin[4v] + (
                     4                         96

          13   3             11   4    5   2             1   3
       - ----*b )sin[3v] + (----*b  - ---*b )sin[2v] + (---*b  - 2*b)
          12                 24        4                 4

                 1   4    1   4
      sin[v] + (---*b  - ---*e )]
                 8        8
 zz:=ff+ww;


               103   4                13   3             11   4
zz :=  - [( - -----*e )sin[4u] + ( - ----*e )sin[3u] + (----*e
               96                     12                 24

          5   2             1   3                     103   4
       - ---*e )sin[2u] + (---*e  - 2*e)sin[u] + ( - -----*b )sin[4v]
          4                 4                         96

              13   3             11   4    5   2             1   3
       + ( - ----*b )sin[3v] + (----*b  - ---*b )sin[2v] + (---*b
              12                 24        4                 4

                        1   4    1   4
       - 2*b)sin[v] + (---*b  + ---*e  - 2)]
                        8        8

xx:=hsub(fourier((1-d*d)*cos(u)),u,u-v+w-x-y+z,yy,n)+
    hsub(fourier(d*d*cos(v)),v,u+v+w+x+y-z,zz,n);


               625   4                       4     3
xx :=  - [( - -----*e )cos[5u-v+w-x-y+z] + (---*b*e )cos[4u+w-x-y+z]
               384                           3

              4   3                          4     3
       + ( - ---*e )cos[4u-v+w-x-y+z] + ( - ---*b*e )cos[4u-2v+w-x-y+
              3                              3

             9   2  2                       17   2  2
      z] + (---*d *e )cos[3u+v+w+x+y-z] + (----*d *e )sin[3u+v+w+x+y-
             8                              12

             9    2  2                        9    4
      z] + (----*b *e )cos[3u+v+w-x-y+z] + (-----*e )cos[3u+v-w+x+y-z
             64                              128

            9     2                     9   2  2    9   2  2
      ] + (---*b*e )cos[3u+w-x-y+z] + (---*b *e  + ---*d *e
            8                           8           8

          27   4    9   2                          9     2
       + ----*e  - ---*e )cos[3u-v+w-x-y+z] + ( - ---*b*e )cos[3u-2v+
          16        8                              8

                      81   2  2                          2
      w-x-y+z] + ( - ----*b *e )cos[3u-3v+w-x-y+z] + (b*d *e)cos[2u+2
                      64

                         2                          1    3
      v+w+x+y-z] + (2*b*d *e)sin[2u+2v+w+x+y-z] + (----*b *e)cos[2u+2
                                                    12

                     1      3                        2
      v+w-x-y+z] + (----*b*e )cos[2u+2v-w+x+y-z] + (d *e)cos[2u+v+w+x
                     12

                2   2                         1   2
      +y-z] + (---*d *e)sin[2u+v+w+x+y-z] + (---*b *e)cos[2u+v+w-x-y+
                3                             8

             1    3                            2
      z] + (----*e )cos[2u+v-w+x+y-z] + ( - b*d *e)cos[2u+w+x+y-z] + 
             12

               2                            2      5     3
      ( - 2*b*d *e)sin[2u+w+x+y-z] + ( - b*d *e - ---*b*e  + b*e)cos[
                                                   4

                         1      3                     2      2
      2u+w-x-y+z] + ( - ----*b*e )cos[2u-w+x+y-z] + (b *e + d *e
                         12

          5   3                           5   3        2      5     3
       + ---*e  - e)cos[2u-v+w-x-y+z] + (---*b *e + b*d *e + ---*b*e
          4                               4                   4

                                       9   2
       - b*e)cos[2u-2v+w-x-y+z] + ( - ---*b *e)cos[2u-3v+w-x-y+z] + (
                                       8

          4   3                          9   2  2
       - ---*b *e)cos[2u-4v+w-x-y+z] + (---*b *d )cos[u+3v+w+x+y-z]
          3                              8

           17   2  2                        9    4
       + (----*b *d )sin[u+3v+w+x+y-z] + (-----*b )cos[u+3v+w-x-y+z]
           12                              128

           9    2  2                         2
       + (----*b *e )cos[u+3v-w+x+y-z] + (b*d )cos[u+2v+w+x+y-z] + (
           64

       2     2                       1    3                       1
      ---*b*d )sin[u+2v+w+x+y-z] + (----*b )cos[u+2v+w-x-y+z] + (---
       3                             12                           8

          2                          2  2    2  2    1   2
      *b*e )cos[u+2v-w+x+y-z] + ( - b *d  - d *e  + ---*d )cos[u+v+w+
                                                     3

                      2  2      2  2    2   2
      x+y-z] + ( - 2*b *d  - 2*d *e  + ---*d )sin[u+v+w+x+y-z] + (
                                        3

          1    4    1   2  2    1   2  2    1   2
       - ----*b  - ---*b *d  - ---*b *e  + ---*b )cos[u+v+w-x-y+z] + 
          48        8           8           8

           1   2  2    1   2  2    1    4    1   2
      ( - ---*b *e  - ---*d *e  - ----*e  + ---*e )cos[u+v-w+x+y-z]
           8           8           48        8

                2                       2     2
       + ( - b*d )cos[u+w+x+y-z] + ( - ---*b*d )sin[u+w+x+y-z] + (
                                        3

            2      2                           1     2
       - b*d  - b*e  + b)cos[u+w-x-y+z] + ( - ---*b*e )cos[u-w+x+y-z]
                                               8

              1   2  2                      7    2  2
       + ( - ---*b *d )cos[u-v+w+x+y-z] + (----*b *d )sin[u-v+w+x+y-z
              8                             12

               7    4    2  2    2  2    2    2  2    2    7    4
      ] + ( - ----*b  - b *d  - b *e  + b  - d *e  + d  - ----*e
               64                                          64

               2                          1   4    1   4
            + e  - 1)cos[u-v+w-x-y+z] + (---*b  - ---*e )sin[u-v+w-x-
                                          8        8

                  1    2  2                         1   2  2
      y+z] + ( - ----*b *e )cos[u-v-w+x+y-z] + ( - ---*d *e )cos[u-v-
                  64                                8

                      7    2  2                      5   3      2
      w-x-y+z] + ( - ----*d *e )sin[u-v-w-x-y+z] + (---*b  + b*d
                      12                             4

            2                           27   4    9   2  2
       + b*e  - b)cos[u-2v+w-x-y+z] + (----*b  + ---*b *d
                                        16        8

          9   2  2    9   2                          4   3
       + ---*b *e  - ---*b )cos[u-3v+w-x-y+z] + ( - ---*b )cos[u-4v+w
          8           8                              3

                     625   4                       4   3
      -x-y+z] + ( - -----*b )cos[u-5v+w-x-y+z] + (---*b *e)cos[4v-w+x
                     384                           3

                9   2                            2
      +y-z] + (---*b *e)cos[3v-w+x+y-z] + ( - b*d *e)cos[2v+w+x+y-z]
                8

                  2                          1    3
       + ( - 2*b*d *e)sin[2v+w+x+y-z] + ( - ----*b *e)cos[2v+w-x-y+z]
                                             12

              5   3        2                                2
       + ( - ---*b *e - b*d *e + b*e)cos[2v-w+x+y-z] + ( - d *e)cos[v
              4

                       2   2                         1   2
      +w+x+y-z] + ( - ---*d *e)sin[v+w+x+y-z] + ( - ---*b *e)cos[v+w-
                       3                             8

                    2      2                            2
      x-y+z] + ( - b *e - d *e + e)cos[v-w+x+y-z] + (b*d *e)cos[w+x+y

                  2                      2
      -z] + (2*b*d *e)sin[w+x+y-z] + (b*d *e - b*e)cos[w-x-y+z]]


%% Step 5: Calculate R
zz:=bb*vv;


            1   4                3      3                  3   3
zz := [( - ---*e )cos[4u] + ( - ----*b*e )cos[3u+v] + ( - ---*e )cos[
            3                    16                        8

                 3      3                  1   2  2
      3u] + ( - ----*b*e )cos[3u-v] + ( - ---*b *e )cos[2u+2v] + (
                 16                        4

          1     2               1   4    1   2                1     2
       - ---*b*e )cos[2u+v] + (---*e  - ---*e )cos[2u] + ( - ---*b*e
          4                     3        2                    4

                        1   2  2                   9    3
      )cos[2u-v] + ( - ---*b *e )cos[2u-2v] + ( - ----*b *e)cos[u+3v]
                        4                          16

              1   2                 1    3      3      3    1
       + ( - ---*b *e)cos[u+2v] + (----*b *e + ----*b*e  - ---*b*e)
              2                     16          16          2

                   3   3                1    3      3      3
      cos[u+v] + (---*e  - e)cos[u] + (----*b *e + ----*b*e
                   8                    16          16

          1                      1   2                    9    3
       - ---*b*e)cos[u-v] + ( - ---*b *e)cos[u-2v] + ( - ----*b *e)
          2                      2                        16

                    4   4             9   3                1   4
      cos[u-3v] + (---*b )cos[4v] + (---*b )cos[3v] + ( - ---*b
                    3                 8                    3

          1   2  2    2                1   3    1     2
       + ---*b *e  + b )cos[2v] + ( - ---*b  + ---*b*e  + b)cos[v] + 
          2                            8        2

        1   2
      (---*e  + 1)]
        2
 yy:=zz*zz*vv;


            1   4                3     3                  1   3
yy := [( - ---*e )cos[4u] + ( - ---*b*e )cos[3u+v] + ( - ---*e )cos[3
            6                    8                        4

                3     3                  9   2  2
      u] + ( - ---*b*e )cos[3u-v] + ( - ---*b *e )cos[2u+2v] + (
                8                        8

          3     2                  3   2  2    1   4    1   2
       - ---*b*e )cos[2u+v] + ( - ---*b *e  + ---*e  - ---*e )cos[2u]
          4                        4           6        2

              3     2                  9   2  2
       + ( - ---*b*e )cos[2u-v] + ( - ---*b *e )cos[2u-2v] + (
              4                        8

          53   3                    9   2                    27   3
       - ----*b *e)cos[u+3v] + ( - ---*b *e)cos[u+2v] + ( - ----*b *e
          8                         2                        8

          3     3                           2      1   3
       + ---*b*e  - 3*b*e)cos[u+v] + ( - 3*b *e + ---*e  - 2*e)cos[u]
          8                                        4

              27   3      3     3                         9   2
       + ( - ----*b *e + ---*b*e  - 3*b*e)cos[u-v] + ( - ---*b *e)cos
              8           8                               2

                    53   3                 77   4             53   3
      [u-2v] + ( - ----*b *e)cos[u-3v] + (----*b )cos[4v] + (----*b )
                    8                      8                  8

                  7   4    27   2  2    9   2             27   3
      cos[3v] + (---*b  + ----*b *e  + ---*b )cos[2v] + (----*b
                  2        4            2                 8

          9     2                  15   4    9   2  2    3   2
       + ---*b*e  + 3*b)cos[v] + (----*b  + ---*b *e  + ---*b
          2                        8         4           2

          3   2
       + ---*e  + 1)]
          2


on fourier;


*** Domain mode rational changed to fourier 


p(0):= fourier 1;


p(0) := [1]
 p(1) := xx;


                 625   4                       4     3
p(1) :=  - [( - -----*e )cos[5u-v+w-x-y+z] + (---*b*e )cos[4u+w-x-y+z
                 384                           3

                 4   3                          4     3
        ] + ( - ---*e )cos[4u-v+w-x-y+z] + ( - ---*b*e )cos[4u-2v+w-x
                 3                              3

                  9   2  2                       17   2  2
        -y+z] + (---*d *e )cos[3u+v+w+x+y-z] + (----*d *e )sin[3u+v+w
                  8                              12

                    9    2  2                        9    4
        +x+y-z] + (----*b *e )cos[3u+v+w-x-y+z] + (-----*e )cos[3u+v-
                    64                              128

                     9     2                     9   2  2    9   2  2
        w+x+y-z] + (---*b*e )cos[3u+w-x-y+z] + (---*b *e  + ---*d *e
                     8                           8           8

            27   4    9   2                          9     2
         + ----*e  - ---*e )cos[3u-v+w-x-y+z] + ( - ---*b*e )cos[3u-2
            16        8                              8

                          81   2  2                          2
        v+w-x-y+z] + ( - ----*b *e )cos[3u-3v+w-x-y+z] + (b*d *e)cos[
                          64

                               2                          1    3
        2u+2v+w+x+y-z] + (2*b*d *e)sin[2u+2v+w+x+y-z] + (----*b *e)
                                                          12

                               1      3                        2
        cos[2u+2v+w-x-y+z] + (----*b*e )cos[2u+2v-w+x+y-z] + (d *e)
                               12

                              2   2                         1   2
        cos[2u+v+w+x+y-z] + (---*d *e)sin[2u+v+w+x+y-z] + (---*b *e)
                              3                             8

                              1    3                            2
        cos[2u+v+w-x-y+z] + (----*e )cos[2u+v-w+x+y-z] + ( - b*d *e)
                              12

                                   2                            2
        cos[2u+w+x+y-z] + ( - 2*b*d *e)sin[2u+w+x+y-z] + ( - b*d *e

            5     3                              1      3
         - ---*b*e  + b*e)cos[2u+w-x-y+z] + ( - ----*b*e )cos[2u-w+x+
            4                                    12

                 2      2      5   3
        y-z] + (b *e + d *e + ---*e  - e)cos[2u-v+w-x-y+z] + (
                               4

         5   3        2      5     3
        ---*b *e + b*d *e + ---*b*e  - b*e)cos[2u-2v+w-x-y+z] + (
         4                   4

            9   2                             4   3
         - ---*b *e)cos[2u-3v+w-x-y+z] + ( - ---*b *e)cos[2u-4v+w-x-y
            8                                 3

                9   2  2                       17   2  2
        +z] + (---*b *d )cos[u+3v+w+x+y-z] + (----*b *d )sin[u+3v+w+x
                8                              12

                   9    4                       9    2  2
        +y-z] + (-----*b )cos[u+3v+w-x-y+z] + (----*b *e )cos[u+3v-w+
                  128                           64

                     2                       2     2
        x+y-z] + (b*d )cos[u+2v+w+x+y-z] + (---*b*d )sin[u+2v+w+x+y-z
                                             3

              1    3                       1     2
        ] + (----*b )cos[u+2v+w-x-y+z] + (---*b*e )cos[u+2v-w+x+y-z]
              12                           8

                2  2    2  2    1   2                           2  2
         + ( - b *d  - d *e  + ---*d )cos[u+v+w+x+y-z] + ( - 2*b *d
                                3

              2  2    2   2                         1    4
         - 2*d *e  + ---*d )sin[u+v+w+x+y-z] + ( - ----*b
                      3                             48

            1   2  2    1   2  2    1   2
         - ---*b *d  - ---*b *e  + ---*b )cos[u+v+w-x-y+z] + (
            8           8           8

            1   2  2    1   2  2    1    4    1   2
         - ---*b *e  - ---*d *e  - ----*e  + ---*e )cos[u+v-w+x+y-z]
            8           8           48        8

                  2                       2     2
         + ( - b*d )cos[u+w+x+y-z] + ( - ---*b*d )sin[u+w+x+y-z] + (
                                          3

              2      2                           1     2
         - b*d  - b*e  + b)cos[u+w-x-y+z] + ( - ---*b*e )cos[u-w+x+y-
                                                 8

                  1   2  2                      7    2  2
        z] + ( - ---*b *d )cos[u-v+w+x+y-z] + (----*b *d )sin[u-v+w+x
                  8                             12

                     7    4    2  2    2  2    2    2  2    2
        +y-z] + ( - ----*b  - b *d  - b *e  + b  - d *e  + d
                     64

                     7    4    2                          1   4
                  - ----*e  + e  - 1)cos[u-v+w-x-y+z] + (---*b
                     64                                   8

            1   4                         1    2  2
         - ---*e )sin[u-v+w-x-y+z] + ( - ----*b *e )cos[u-v-w+x+y-z]
            8                             64

                1   2  2                         7    2  2
         + ( - ---*d *e )cos[u-v-w-x-y+z] + ( - ----*d *e )sin[u-v-w-
                8                                12

                   5   3      2      2
        x-y+z] + (---*b  + b*d  + b*e  - b)cos[u-2v+w-x-y+z] + (
                   4

         27   4    9   2  2    9   2  2    9   2
        ----*b  + ---*b *d  + ---*b *e  - ---*b )cos[u-3v+w-x-y+z] + 
         16        8           8           8

             4   3                          625   4
        ( - ---*b )cos[u-4v+w-x-y+z] + ( - -----*b )cos[u-5v+w-x-y+z]
             3                              384

             4   3                       9   2
         + (---*b *e)cos[4v-w+x+y-z] + (---*b *e)cos[3v-w+x+y-z] + (
             3                           8

              2                              2
         - b*d *e)cos[2v+w+x+y-z] + ( - 2*b*d *e)sin[2v+w+x+y-z] + (

            1    3                          5   3        2
         - ----*b *e)cos[2v+w-x-y+z] + ( - ---*b *e - b*d *e + b*e)
            12                              4

                               2                         2   2
        cos[2v-w+x+y-z] + ( - d *e)cos[v+w+x+y-z] + ( - ---*d *e)sin[
                                                         3

                          1   2                         2      2
        v+w+x+y-z] + ( - ---*b *e)cos[v+w-x-y+z] + ( - b *e - d *e
                          8

                                  2                        2
         + e)cos[v-w+x+y-z] + (b*d *e)cos[w+x+y-z] + (2*b*d *e)sin[w+

                     2
        x+y-z] + (b*d *e - b*e)cos[w-x-y+z]]

for i := 2:n/2+2 do <<
  wtlevel n+4-2i;
  p(i) := fourier ((2*i-1)/i)*xx*p(i-1) - fourier ((i-1)/i)*p(i-2);
>>;



wtlevel n;


for i:=n/2+2 step -1 until 3 do p(n/2+2):=fourier(a*a)*zz*p(n/2+2)+p(i-1);



yy*p(n/2+2);


   27   4                               25     3
[(----*e )cos[6u-2v+2w-2x-2y+2z] + ( - ----*b*e )cos[5u-v+2w-2x-2y+2z
   32                                   64

      25   3                            75   2  2
] + (----*e )cos[5u-2v+2w-2x-2y+2z] + (----*a *e )cos[5u-3v+3w-3x-3y+
      32                                64

        175     3                               13   2  2
3z] + (-----*b*e )cos[5u-3v+2w-2x-2y+2z] + ( - ----*d *e )cos[4u+2w]
        64                                      8

          2  2                   1    4                3     2
 + ( - 2*d *e )sin[4u+2w] + ( - ----*e )cos[4u] + ( - ---*b*e )cos[4u
                                 24                    8

                       15   2
-v+2w-2x-2y+2z] + ( - ----*a *b*e)cos[4u-2v+3w-3x-3y+3z] + (
                       16

    15   2  2    3   2  2    15   4    3   2
 - ----*b *e  - ---*d *e  - ----*e  + ---*e )cos[4u-2v+2w-2x-2y+2z]
    8            2           8         4

     15   2                              21     2
 + (----*a *e)cos[4u-3v+3w-3x-3y+3z] + (----*b*e )cos[4u-3v+2w-2x-2y+
     16                                  8

        35   4                            75   2
2z] + (----*a )cos[4u-4v+4w-4x-4y+4z] + (----*a *b*e)cos[4u-4v+3w-3x-
        64                                16

           51   2  2                               9     2
3y+3z] + (----*b *e )cos[4u-4v+2w-2x-2y+2z] + ( - ---*b*d *e)cos[3u+v
           8                                       4

            7     2                    1    3
+2w] + ( - ---*b*d *e)sin[3u+v+2w] + (----*b *e)cos[3u+v+2w-2x-2y+2z]
            2                          64

        3      3                  3   2                     2
 + ( - ----*b*e )cos[3u+v] + ( - ---*d *e)cos[3u+2w] + ( - d *e)sin[3
        32                        2

             1    3                5   2  2
u+2w] + ( - ----*e )cos[3u] + ( - ---*a *d )cos[3u-v+3w-x-y+z] + (
             16                    8

    5   2  2                        5    2  2
 - ---*a *d )sin[3u-v+3w-x-y+z] + (----*a *b )cos[3u-v+3w-3x-3y+3z]
    4                               64

        9     2                    1     2
 + ( - ---*b*d *e)cos[3u-v+2w] + (---*b*d *e)sin[3u-v+2w] + (
        4                          2

 3    3      3     2      57     3    3
----*b *e + ---*b*d *e + ----*b*e  - ---*b*e)cos[3u-v+2w-2x-2y+2z] + 
 64          4            64          8

     9    2  2                          3      3
( - ----*a *e )cos[3u-v+w-x-y+z] + ( - ----*b*e )cos[3u-v] + (
     64                                 32

    5   2                                 15   2      3   2
 - ---*a *b)cos[3u-2v+3w-3x-3y+3z] + ( - ----*b *e - ---*d *e
    8                                     8           2

    57   3    3                                  15   2  2
 - ----*e  + ---*e)cos[3u-2v+2w-2x-2y+2z] + ( - ----*a *b
    32        4                                  4

    15   2  2    15   2  2    5   2
 - ----*a *d  - ----*a *e  + ---*a )cos[3u-3v+3w-3x-3y+3z] + (
    8            4            8

    369   3      21     2      399     3    21
 - -----*b *e - ----*b*d *e - -----*b*e  + ----*b*e)cos[3u-3v+2w-2x-2
    64           4             64           8

          25   2                              51   2
y+2z] + (----*a *b)cos[3u-4v+3w-3x-3y+3z] + (----*b *e)cos[3u-4v+2w-2
          8                                   8

             635   2  2                            845   3
x-2y+2z] + (-----*a *b )cos[3u-5v+3w-3x-3y+3z] + (-----*b *e)cos[3u-5
             64                                    64

                      1   4                            1   4
v+2w-2x-2y+2z] + ( - ---*d )cos[2u+2v+2w+2x+2y-2z] + (---*d )sin[2u+2
                      4                                3

                      11   2  2                      13   2  2
v+2w+2x+2y-2z] + ( - ----*b *d )cos[2u+2v+2w] + ( - ----*b *d )sin[2u
                      4                              4

            1    4                            2  2
+2v+2w] + (----*b )cos[2u+2v+2w-2x-2y+2z] + (d *e )cos[2u+2v+2x+2y-2z
            32

         3   2  2                            9    2  2
] + ( - ---*d *e )sin[2u+2v+2x+2y-2z] + ( - ----*b *e )cos[2u+2v] + (
         4                                   32

    3    4                               7     2
 - ----*e )cos[2u+2v-2w+2x+2y-2z] + ( - ---*b*d )cos[2u+v+2w] + (
    64                                   4

    3     2                  1    3
 - ---*b*d )sin[2u+v+2w] + (----*b )cos[2u+v+2w-2x-2y+2z] + (
    2                        64

    3      2                  7   2  2    1   4    17   2  2    1   2
 - ----*b*e )cos[2u+v] + ( - ---*b *d  + ---*d  + ----*d *e  - ---*d
    16                        4           2        4            2

                1   2  2    4    9   2  2    2
)cos[2u+2w] + (---*b *d  + d  + ---*d *e  - d )sin[2u+2w] + (
                2                2

    3    2                            3    2  2    3   2  2    1    4
 - ----*a *b*e)cos[2u+w-x-y+z] + ( - ----*b *e  + ---*d *e  + ----*e
    16                                16           4           24

    1   2             1     2                     3     2
 - ---*e )cos[2u] + (---*b*d )cos[2u-v+2w] + ( - ---*b*d )sin[2u-v+2w
    8                 4                           2

      3    3    3     2    15     2    3
] + (----*b  + ---*b*d  + ----*b*e  - ---*b)cos[2u-v+2w-2x-2y+2z] + (
      64        4          16          8

    3    2                            3      2               45   2
 - ----*a *e)cos[2u-v+w-x-y+z] + ( - ----*b*e )cos[2u-v] + (----*a *b
    16                                16                     16

                              3   2  2                      13   2  2
*e)cos[2u-2v+3w-3x-3y+3z] + (---*b *d )cos[2u-2v+2w] + ( - ----*b *d
                              2                             4

                   5    4    39   4    15   2  2    75   2  2
)sin[2u-2v+2w] + (----*a  + ----*b  + ----*b *d  + ----*b *e
                   16        64        4            16

                      15   2    3   4    15   2  2    3   2    69   4
                   - ----*b  + ---*d  + ----*d *e  - ---*d  + ----*e
                      8         4        4            2        64

                      15   2    3
                   - ----*e  + ---)cos[2u-2v+2w-2x-2y+2z] + (
                      8         4

    3    4    3    4                               9    2
 - ----*b  + ----*e )sin[2u-2v+2w-2x-2y+2z] + ( - ----*a *b*e)cos[2u-
    16        16                                   16

                   9    2  2                1   2  2
2v+w-x-y+z] + ( - ----*b *e )cos[2u-2v] + (---*d *e )cos[2u-2v-2x-2y+
                   32                       4

        3   2  2                            45   2
2z] + (---*d *e )sin[2u-2v-2x-2y+2z] + ( - ----*a *e)cos[2u-3v+3w-3x-
        4                                   16

              369   3    21     2    105     2    21
3y+3z] + ( - -----*b  - ----*b*d  - -----*b*e  + ----*b)cos[2u-3v+2w-
              64         4           16           8

                 225   2                                   115   4
2x-2y+2z] + ( - -----*a *b*e)cos[2u-4v+3w-3x-3y+3z] + ( - -----*b
                 16                                         8

    51   2  2    255   2  2    51   2                            845
 - ----*b *d  - -----*b *e  + ----*b )cos[2u-4v+2w-2x-2y+2z] + (-----
    4            16            8                                 64

  3                            1599   4                            1
*b )cos[2u-5v+2w-2x-2y+2z] + (------*b )cos[2u-6v+2w-2x-2y+2z] + (---
                                64                                 4

    2                          3     2
*b*d *e)cos[u+3v+2x+2y-2z] + (---*b*d *e)sin[u+3v+2x+2y-2z] + (
                               2

    53   3                    49     3
 - ----*b *e)cos[u+3v] + ( - ----*b*e )cos[u+3v-2w+2x+2y-2z] + (
    32                        64

    1   2                          2
 - ---*d *e)cos[u+2v+2x+2y-2z] + (d *e)sin[u+2v+2x+2y-2z] + (
    2

    9   2                    7    3                           23
 - ---*b *e)cos[u+2v] + ( - ----*e )cos[u+2v-2w+2x+2y-2z] + (----*b
    8                        32                               4

  2                   13     2                      3    3
*d *e)cos[u+v+2w] + (----*b*d *e)sin[u+v+2w] + ( - ----*b *e)cos[u+v+
                      2                             64

                    3   2  2                         3   2  2
2w-2x-2y+2z] + ( - ---*a *d )cos[u+v+w+x+y-z] + ( - ---*a *d )sin[u+v
                    4                                2

              33   2  2                         7     2
+w+x+y-z] + (----*a *b )cos[u+v+w-x-y+z] + ( - ---*b*d *e)cos[u+v+2x+
              64                                4

           3     2                            27   3      9     2
2y-2z] + (---*b*d *e)sin[u+v+2x+2y-2z] + ( - ----*b *e + ---*b*d *e
           2                                  32          2

    3      3    3                   33   2  2
 + ----*b*e  - ---*b*e)cos[u+v] + (----*a *e )cos[u+v-w+x+y-z] + (
    32          4                   64

 7      3                          5   2                   2
----*b*e )cos[u+v-2w+2x+2y-2z] + (---*d *e)cos[u+2w] + (3*d *e)sin[u+
 64                                2

        3   2                         3   2        2      1    3
2w] + (---*a *b)cos[u+w-x-y+z] + ( - ---*b *e + 3*d *e + ----*e
        8                             4                   16

    1               7     2                   5     2
 - ---*e)cos[u] + (---*b*d *e)cos[u-v+2w] + (---*b*d *e)sin[u-v+2w]
    2               4                         2

        9    3      9     2      39     3    9
 + ( - ----*b *e - ---*b*d *e - ----*b*e  + ---*b*e)cos[u-v+2w-2x-2y+
        64          4            64          8

        3   2  2    33   2  2    3   2  2    3   2
2z] + (---*a *b  - ----*a *d  + ---*a *e  + ---*a )cos[u-v+w-x-y+z]
        4           8            4           8

        27   3      9     2      3      3    3
 + ( - ----*b *e + ---*b*d *e + ----*b*e  - ---*b*e)cos[u-v] + (
        32          2            32          4

    3     2                         5     2
 - ---*b*d *e)cos[u-v-2x-2y+2z] + (---*b*d *e)sin[u-v-2x-2y+2z] + (
    4                               2

 45   2      9   2      39   3    9                              9
----*b *e + ---*d *e + ----*e  - ---*e)cos[u-2v+2w-2x-2y+2z] + (---
 8           2          32        4                              8

  2                            9   2                 3   2
*a *b)cos[u-2v+w-x-y+z] + ( - ---*b *e)cos[u-2v] + (---*d *e)cos[u-2v
                               8                     2

                  2                          285   2  2
-2x-2y+2z] + ( - d *e)sin[u-2v-2x-2y+2z] + (-----*a *e )cos[u-3v+3w-3
                                             64

             1107   3      63     2      273     3    63
x-3y+3z] + (------*b *e + ----*b*d *e + -----*b*e  - ----*b*e)cos[u-3
              64           4             64           8

                   159   2  2                          5   2  2
v+2w-2x-2y+2z] + (-----*a *b )cos[u-3v+w-x-y+z] + ( - ---*a *d )cos[u
                   64                                  8

                    5   2  2                             53   3
-3v+w-3x-3y+3z] + (---*a *d )sin[u-3v+w-3x-3y+3z] + ( - ----*b *e)cos
                    4                                    32

           21     2                             11     2
[u-3v] + (----*b*d *e)cos[u-3v-2x-2y+2z] + ( - ----*b*d *e)sin[u-3v-2
           4                                    2

                153   2                                2535   3
x-2y+2z] + ( - -----*b *e)cos[u-4v+2w-2x-2y+2z] + ( - ------*b *e)cos
                 8                                      64

                          63   2  2                         19   2  2
[u-5v+2w-2x-2y+2z] + ( - ----*b *d )cos[4v+2x+2y-2z] + ( - ----*b *d
                          8                                 2

                      77   4             255   2  2
)sin[4v+2x+2y-2z] + (----*b )cos[4v] + (-----*b *e )cos[4v-2w+2x+2y-2
                      32                 16

          11     2                         7     2
z] + ( - ----*b*d )cos[3v+2x+2y-2z] + ( - ---*b*d )sin[3v+2x+2y-2z]
          4                                2

     53   3             105     2                         17   2  2
 + (----*b )cos[3v] + (-----*b*e )cos[3v-2w+2x+2y-2z] + (----*b *d
     32                 16                                4

    1   4    7   2  2    1   2                      9   2  2    4
 + ---*d  - ---*d *e  - ---*d )cos[2v+2x+2y-2z] + (---*b *d  + d
    2        4           2                          2

    1   2  2    2                      7   4    27   2  2
 + ---*d *e  - d )sin[2v+2x+2y-2z] + (---*b  - ----*b *d
    2                                  8        4

    27   2  2    9   2                45   2
 + ----*b *e  + ---*b )cos[2v] + ( - ----*a *b*e)cos[2v-w+x+y-z] + (
    16           8                    16

    75   2  2    15   2  2    15   2                         5     2
 - ----*b *e  - ----*d *e  + ----*e )cos[2v-2w+2x+2y-2z] + (---*b*d )
    16           4            8                              4

                    1     2                     27   3    9     2
cos[v+2x+2y-2z] + (---*b*d )sin[v+2x+2y-2z] + (----*b  - ---*b*d
                    2                           32        2

    9     2    3                  15   2
 + ---*b*e  + ---*b)cos[v] + ( - ----*a *e)cos[v-w+x+y-z] + (
    8          4                  16

    15     2                           25   2  2
 - ----*b*e )cos[v-2w+2x+2y-2z] + ( - ----*d *e )cos[2w] + (
    16                                 8

    7   2  2                15   2                      5   2  2
 - ---*d *e )sin[2w] + ( - ----*a *b*e)cos[w-x-y+z] + (---*b *d )cos[
    2                       16                          8

                 2  2                   9    4    15   4    9   2  2
2x+2y-2z] + ( - b *d )sin[2x+2y-2z] + (----*a  + ----*b  - ---*b *d
                                        64        32        4

    9    2  2    3   2    7   4    9   2  2    3   2    3   2    1
 + ----*b *e  + ---*b  + ---*d  - ---*d *e  - ---*d  + ---*e  + ---)]
    16           8        6        4           2        8        4


showtime;


Time: 28616 ms  plus GC time: 500 ms
end;
(camal 28616 500)


End of Lisp run after 28.63+1.16 seconds
