フーリエ変換サブルーチン
V2 Fourier Transformation Subroutine
Closed; even; 50 storage locations + 34 storage locations for S2 
and 16 storage locations for working spaces.
Working spaces: 0w 〜 15w.
Puts Cm in 12wL and Sm in 14wL, where 2m = C(3)
and

Cm =  2-15Σr=0n-1arcos(2mrπ/n),
Sm =  2-15Σr=0n-1arsin(2mrπ/n),

where ar=C(A+r).

Presets parameters:
  0n = number of points,
  0a = A location of a0
  0h = 15,       ┐
  0s = 50f,      │ set by the subroutine tape.
  0w = 34s = 84f,┘

Notes:
1. Scale factor is set to 2-15, but it can be changes to 2-h
by re-setting the h-parameter.
2. The subroutine places -2m/n ,2-3cos(2mπ/n) and
 2-3sin(2mπ/n) in 6wL, 8wL, 10wL respectively.


 0f: a 40   ┐plant link
 1   x 46r  ┘
 2   q  2   ┐
 3   n  3   │-2m/n to 6wL
 4   d 62   │
 5   tl 6w  ┘
 6   tlo    ┐
 7   it     │2-3cos(2mπ/n)
 8   jls    │to 8wL
 9   tl 8w  ┘
10   p  2   ┐
11   tlw    │clear working spaces
12   tl 4w  ┘
13   p 47r  ┐
14   x 25r  │
15   s 48r  │
16   z 32r  │
17   q 49r  │
18   pl 4w  │
19   wl 8w  │
20   ll 4   │
21   slw    │
22   tl 2w  │iterative cycle
23   pl 4w  │
24   tlw    │to compute



25   p (  ) │Si+1=
26   rh     │2Sicos(2mπ/n)
27   al 2w  │-Si-1+2-hai
28   tl 4w  │
29   p 25r  │
30   a 49   │
31   jl 14r ┘
32   pl 4w  ┐
33   vl 8w  │Cm to 12wL
34   ll 3   │
35   slw    │
36   tl 12w ┘
37   nl 6w  ┐
38   a 21s  │
39   tlo    │2-3sin(2mπ/n)
40   it     │to 10wL
41   jls    │
42   tl 10w ┘
43   vl 4w  ┐
44   ll 3   │Sm to 14wL
45   tl14w  ┘
46   jl (  )  link
47 | 0a
48 | 0na
49 | t        round off constant

                                        Date   60-01-11
                                        Author H. Takahasi