Artifact a20301fe508a879e019b1bd14b2b5ee730ca25a07eacb0cc2df8cbff1e88b0cb:
- File
r35/xlog/camal.log
— part of check-in
[f2fda60abd]
at
2011-09-02 18:13:33
on branch master
— Some historical releases purely for archival purposes
git-svn-id: https://svn.code.sf.net/p/reduce-algebra/code/trunk/historical@1375 2bfe0521-f11c-4a00-b80e-6202646ff360 (user: arthurcnorman@users.sourceforge.net, size: 37514) [annotate] [blame] [check-ins using] [more...]
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