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