!AWJan92/Goodies/ProgTheARM/FFT5

Home » Recent acquisitions » Acorn ADFS disks » adfs_ArchimedesWorld_199201.adf » January92 » !AWJan92/Goodies/ProgTheARM/FFT5

This website contains an archive of files for the Acorn Electron, BBC Micro, Acorn Archimedes, Commodore 16 and Commodore 64 computers, which Dominic Ford has rescued from his private collection of floppy disks and cassettes.

Some of these files were originally commercial releases in the 1980s and 1990s, but they are now widely available online. I assume that copyright over them is no longer being asserted. If you own the copyright and would like files to be removed, please contact me.

Tape/disk:	Home » Recent acquisitions » Acorn ADFS disks » adfs_ArchimedesWorld_199201.adf » January92
Filename:	!AWJan92/Goodies/ProgTheARM/FFT5
Read OK:	✔
File size:	094B bytes
Load address:	0000
Exec address:	0000

File contents

   10REM > FFT5
   20
   30MODE 12
   40
   50N%=128
   60
   70n2%=2*N%
   80
   90DIM data(n2%)
   91REPEAT
  100K=RND(15)
  110J=RND(96)+16
  120CLS
  130PRINT"Frequencies: "K,J
  140
  150GCOL7
  160
  170FOR I%=0TO N%-1
  180  ii%=2*I%+1
  190  data(ii%)   = SINRAD(I%*K) + SINRAD(I%*J)
  200  data(ii%+1) = 0
  210  MOVE 4*I%,512
  220  DRAW 4*I%,512+256*data(ii%)
  230NEXT
  240
  250PROCfour1(N%,1)
  260
  270GCOL1
  280
  290FOR I%=0 TO (N%-1) DIV 2
  300  ii%=2*I%+1
  310  MOVE 8*I%,0
  320  DRAW 8*I%,8*ABS(data(ii%))
  330NEXT
  340
  350REM How about chucking out everything above the ?th bucket
  360
  370FOR I% = (N% DIV 8) TO (7 * N%) DIV 8
  380  ii%=2*I%+1
  390  data(ii%)  =0
  400  data(ii%+1)=0
  420NEXT
  430
  440GCOL2
  450
  460FOR I%=0 TO (N%-1) DIV 2
  470  ii%=2*I%+1
  480  MOVE 640+8*I%,0
  490  DRAW 640+8*I%,8*ABS(data(ii%))
  500NEXT
  510
  520
  530REM Do the inverse and draw it again
  540
  550PROCfour1(N%,-1)
  560
  570GCOL3
  580
  590FOR I%=0TO N%-1
  600  ii%=2*I%+1
  610  MOVE 640+4*I%,512
  620  DRAW 640+4*I%,512+256*data(ii%)/N%
  630NEXT
  640
  641UNTILGET=13
  650END
  660
  670
  680REM data() accessed by refernece
  690DEFPROCfour1(nn%,isign%)
  700
  710LOCALii%,jj%,n%,mmax%,m%,j%,istep%,i%
  720LOCALwtemp,wr,wpr,wpi,wi,theta
  730LOCALtempr,tempi,wrs,wis
  740
  750LOCALs
  760
  770n%=2*nn%
  780j%=1
  790REM This is the bit-reversal section of the routine
  800FOR ii%=1 TO nn%
  810  i%=2*ii%-1
  820  IF j% > i% THEN
  830    REM Swap the real and imarinary parts of the two
  840    tempr      = data(j%) : tempi      = data(j%+1)
  850    data(j%)   = data(i%) : data(j%+1) = data(i%+1)
  860    data(i%)   = tempr    : data(i%+1) = tempi
  870  ENDIF
  880  m% = n% DIV 2
  890  WHILE (m%>=2) AND (j%>m%)
  900    j% = j% - m%
  910    m% = m% DIV 2
  920  ENDWHILE
  930  j% = j% + m%
  940NEXT
  950
  960mmax% = 2
  970
  980REM Now the Danielson-Lanczos recurrence
  990
 1000WHILE n% > mmax%
 1010  istep% = 2*mmax%
 1020  theta  = 2 * PI / (isign%*mmax%)
 1030
 1040        s=SIN(0.5*theta)
 1050  wpr    = -2.0 * s*s
 1060  wpi    = SIN(theta)
 1070
 1080  wr = 1.0
 1090  wi = 0.0
 1100  FOR ii% = 1 TO mmax% DIV 2
 1110    m% = 2*ii%-1
 1120    wrs = wr
 1130    wis = wi
 1140    FOR jj%=0 TO (n%-m%) DIV istep%
 1150      i% = m% + jj%*istep%
 1160      j% = i%+mmax%
 1170      tempr = wrs*data(j%)   - wis*data(j%+1)
 1180      tempi = wrs*data(j%+1) + wis*data(j%)
 1190      data(j%)   = data(i%)   - tempr
 1200      data(j%+1) = data(i%+1) - tempi
 1210      data(i%)   = data(i%)   + tempr
 1220      data(i%+1) = data(i%+1) + tempi
 1230    NEXT
 1240    wtemp = wr
 1250    wr = wr*wpr -    wi*wpi + wr
 1260    wi = wi*wpr + wtemp*wpi  +wi
 1270  NEXT
 1280  mmax% = istep%
 1290ENDWHILE
 1300ENDPROC
 1310
 1320
 1330
 1340
 1350
 1360
 1370
 1380
 1390
 1400
 1410
 1420
 1430
 1440
 1450
 1460


� > FFT5

� 12
(
2
N%=128
<
Fn2%=2*N%
P
Z� data(n2%)
[�
dK=�(15)
nJ=�(96)+16
x�
��"Frequencies: "K,J
�
��7
�
�� I%=0� N%-1
�  ii%=2*I%+1
�'  data(ii%)   = ��(I%*K) + ��(I%*J)
�  data(ii%+1) = 0
�  � 4*I%,512
�  � 4*I%,512+256*data(ii%)
��
�
��four1(N%,1)

�1

"� I%=0 � (N%-1) � 2
,  ii%=2*I%+1
6  � 8*I%,0
@  � 8*I%,8*�(data(ii%))
J�
T
^<� How about chucking out everything above the ?th bucket
h
r"� I% = (N% � 8) � (7 * N%) � 8
|  ii%=2*I%+1
�  data(ii%)  =0
�  data(ii%+1)=0
��
�
��2
�
�� I%=0 � (N%-1) � 2
�  ii%=2*I%+1
�  � 640+8*I%,0
�  � 640+8*I%,8*�(data(ii%))
��
�

&� Do the inverse and draw it again

&�four1(N%,-1)
0
:�3
D
N� I%=0� N%-1
X  ii%=2*I%+1
b  � 640+4*I%,512
l%  � 640+4*I%,512+256*data(ii%)/N%
v�
�
�	��=13
��
�
�
�"� data() accessed by refernece
���four1(nn%,isign%)
�
�%�ii%,jj%,n%,mmax%,m%,j%,istep%,i%
��wtemp,wr,wpr,wpi,wi,theta
��tempr,tempi,wrs,wis
�
��s
�
n%=2*nn%
j%=1
5� This is the bit-reversal section of the routine
 � ii%=1 � nn%
*  i%=2*ii%-1
4  � j% > i% �
>6    � Swap the real and imarinary parts of the two
H7    tempr      = data(j%) : tempi      = data(j%+1)
R7    data(j%)   = data(i%) : data(j%+1) = data(i%+1)
\2    data(i%)   = tempr    : data(i%+1) = tempi
f  �
p  m% = n% � 2
z  ȕ (m%>=2) � (j%>m%)
�    j% = j% - m%
�    m% = m% � 2
�  �
�  j% = j% + m%
��
�
�
mmax% = 2
�
�*� Now the Danielson-Lanczos recurrence
�
�ȕ n% > mmax%
�  istep% = 2*mmax%
�%  theta  = 2 * � / (isign%*mmax%)

        s=�(0.5*theta)
  wpr    = -2.0 * s*s
$  wpi    = �(theta)
.
8  wr = 1.0
B  wi = 0.0
L  � ii% = 1 � mmax% � 2
V    m% = 2*ii%-1
`    wrs = wr
j    wis = wi
t"    � jj%=0 � (n%-m%) � istep%
~      i% = m% + jj%*istep%
�      j% = i%+mmax%
�1      tempr = wrs*data(j%)   - wis*data(j%+1)
�/      tempi = wrs*data(j%+1) + wis*data(j%)
�)      data(j%)   = data(i%)   - tempr
�)      data(j%+1) = data(i%+1) - tempi
�)      data(i%)   = data(i%)   + tempr
�)      data(i%+1) = data(i%+1) + tempi
�	    �
�    wtemp = wr
�$    wr = wr*wpr -    wi*wpi + wr
�$    wi = wi*wpr + wtemp*wpi  +wi
�  �
  mmax% = istep%

�
�

(
2
<
F
P
Z
d
n
x
�
�
�
�
�
�
�

00000000  0d 00 0a 0c f4 20 3e 20  46 46 54 35 0d 00 14 04  |..... > FFT5....|
00000010  0d 00 1e 08 eb 20 31 32  0d 00 28 04 0d 00 32 0a  |..... 12..(...2.|
00000020  4e 25 3d 31 32 38 0d 00  3c 04 0d 00 46 0c 6e 32  |N%=128..<...F.n2|
00000030  25 3d 32 2a 4e 25 0d 00  50 04 0d 00 5a 0f de 20  |%=2*N%..P...Z.. |
00000040  64 61 74 61 28 6e 32 25  29 0d 00 5b 05 f5 0d 00  |data(n2%)..[....|
00000050  64 0b 4b 3d b3 28 31 35  29 0d 00 6e 0e 4a 3d b3  |d.K=.(15)..n.J=.|
00000060  28 39 36 29 2b 31 36 0d  00 78 05 db 0d 00 82 17  |(96)+16..x......|
00000070  f1 22 46 72 65 71 75 65  6e 63 69 65 73 3a 20 22  |."Frequencies: "|
00000080  4b 2c 4a 0d 00 8c 04 0d  00 96 06 e6 37 0d 00 a0  |K,J.........7...|
00000090  04 0d 00 aa 10 e3 20 49  25 3d 30 b8 20 4e 25 2d  |...... I%=0. N%-|
000000a0  31 0d 00 b4 10 20 20 69  69 25 3d 32 2a 49 25 2b  |1....  ii%=2*I%+|
000000b0  31 0d 00 be 27 20 20 64  61 74 61 28 69 69 25 29  |1...'  data(ii%)|
000000c0  20 20 20 3d 20 b5 b2 28  49 25 2a 4b 29 20 2b 20  |   = ..(I%*K) + |
000000d0  b5 b2 28 49 25 2a 4a 29  0d 00 c8 15 20 20 64 61  |..(I%*J)....  da|
000000e0  74 61 28 69 69 25 2b 31  29 20 3d 20 30 0d 00 d2  |ta(ii%+1) = 0...|
000000f0  10 20 20 ec 20 34 2a 49  25 2c 35 31 32 0d 00 dc  |.  . 4*I%,512...|
00000100  1e 20 20 df 20 34 2a 49  25 2c 35 31 32 2b 32 35  |.  . 4*I%,512+25|
00000110  36 2a 64 61 74 61 28 69  69 25 29 0d 00 e6 05 ed  |6*data(ii%).....|
00000120  0d 00 f0 04 0d 00 fa 10  f2 66 6f 75 72 31 28 4e  |.........four1(N|
00000130  25 2c 31 29 0d 01 04 04  0d 01 0e 06 e6 31 0d 01  |%,1).........1..|
00000140  18 04 0d 01 22 17 e3 20  49 25 3d 30 20 b8 20 28  |....".. I%=0 . (|
00000150  4e 25 2d 31 29 20 81 20  32 0d 01 2c 10 20 20 69  |N%-1) . 2..,.  i|
00000160  69 25 3d 32 2a 49 25 2b  31 0d 01 36 0e 20 20 ec  |i%=2*I%+1..6.  .|
00000170  20 38 2a 49 25 2c 30 0d  01 40 1b 20 20 df 20 38  | 8*I%,0..@.  . 8|
00000180  2a 49 25 2c 38 2a 94 28  64 61 74 61 28 69 69 25  |*I%,8*.(data(ii%|
00000190  29 29 0d 01 4a 05 ed 0d  01 54 04 0d 01 5e 3c f4  |))..J....T...^<.|
000001a0  20 48 6f 77 20 61 62 6f  75 74 20 63 68 75 63 6b  | How about chuck|
000001b0  69 6e 67 20 6f 75 74 20  65 76 65 72 79 74 68 69  |ing out everythi|
000001c0  6e 67 20 61 62 6f 76 65  20 74 68 65 20 3f 74 68  |ng above the ?th|
000001d0  20 62 75 63 6b 65 74 0d  01 68 04 0d 01 72 22 e3  | bucket..h...r".|
000001e0  20 49 25 20 3d 20 28 4e  25 20 81 20 38 29 20 b8  | I% = (N% . 8) .|
000001f0  20 28 37 20 2a 20 4e 25  29 20 81 20 38 0d 01 7c  | (7 * N%) . 8..||
00000200  10 20 20 69 69 25 3d 32  2a 49 25 2b 31 0d 01 86  |.  ii%=2*I%+1...|
00000210  13 20 20 64 61 74 61 28  69 69 25 29 20 20 3d 30  |.  data(ii%)  =0|
00000220  0d 01 90 13 20 20 64 61  74 61 28 69 69 25 2b 31  |....  data(ii%+1|
00000230  29 3d 30 0d 01 a4 05 ed  0d 01 ae 04 0d 01 b8 06  |)=0.............|
00000240  e6 32 0d 01 c2 04 0d 01  cc 17 e3 20 49 25 3d 30  |.2......... I%=0|
00000250  20 b8 20 28 4e 25 2d 31  29 20 81 20 32 0d 01 d6  | . (N%-1) . 2...|
00000260  10 20 20 69 69 25 3d 32  2a 49 25 2b 31 0d 01 e0  |.  ii%=2*I%+1...|
00000270  12 20 20 ec 20 36 34 30  2b 38 2a 49 25 2c 30 0d  |.  . 640+8*I%,0.|
00000280  01 ea 1f 20 20 df 20 36  34 30 2b 38 2a 49 25 2c  |...  . 640+8*I%,|
00000290  38 2a 94 28 64 61 74 61  28 69 69 25 29 29 0d 01  |8*.(data(ii%))..|
000002a0  f4 05 ed 0d 01 fe 04 0d  02 08 04 0d 02 12 26 f4  |..............&.|
000002b0  20 44 6f 20 74 68 65 20  69 6e 76 65 72 73 65 20  | Do the inverse |
000002c0  61 6e 64 20 64 72 61 77  20 69 74 20 61 67 61 69  |and draw it agai|
000002d0  6e 0d 02 1c 04 0d 02 26  11 f2 66 6f 75 72 31 28  |n......&..four1(|
000002e0  4e 25 2c 2d 31 29 0d 02  30 04 0d 02 3a 06 e6 33  |N%,-1)..0...:..3|
000002f0  0d 02 44 04 0d 02 4e 10  e3 20 49 25 3d 30 b8 20  |..D...N.. I%=0. |
00000300  4e 25 2d 31 0d 02 58 10  20 20 69 69 25 3d 32 2a  |N%-1..X.  ii%=2*|
00000310  49 25 2b 31 0d 02 62 14  20 20 ec 20 36 34 30 2b  |I%+1..b.  . 640+|
00000320  34 2a 49 25 2c 35 31 32  0d 02 6c 25 20 20 df 20  |4*I%,512..l%  . |
00000330  36 34 30 2b 34 2a 49 25  2c 35 31 32 2b 32 35 36  |640+4*I%,512+256|
00000340  2a 64 61 74 61 28 69 69  25 29 2f 4e 25 0d 02 76  |*data(ii%)/N%..v|
00000350  05 ed 0d 02 80 04 0d 02  81 09 fd a5 3d 31 33 0d  |............=13.|
00000360  02 8a 05 e0 0d 02 94 04  0d 02 9e 04 0d 02 a8 22  |..............."|
00000370  f4 20 64 61 74 61 28 29  20 61 63 63 65 73 73 65  |. data() accesse|
00000380  64 20 62 79 20 72 65 66  65 72 6e 65 63 65 0d 02  |d by refernece..|
00000390  b2 17 dd f2 66 6f 75 72  31 28 6e 6e 25 2c 69 73  |....four1(nn%,is|
000003a0  69 67 6e 25 29 0d 02 bc  04 0d 02 c6 25 ea 69 69  |ign%).......%.ii|
000003b0  25 2c 6a 6a 25 2c 6e 25  2c 6d 6d 61 78 25 2c 6d  |%,jj%,n%,mmax%,m|
000003c0  25 2c 6a 25 2c 69 73 74  65 70 25 2c 69 25 0d 02  |%,j%,istep%,i%..|
000003d0  d0 1e ea 77 74 65 6d 70  2c 77 72 2c 77 70 72 2c  |...wtemp,wr,wpr,|
000003e0  77 70 69 2c 77 69 2c 74  68 65 74 61 0d 02 da 18  |wpi,wi,theta....|
000003f0  ea 74 65 6d 70 72 2c 74  65 6d 70 69 2c 77 72 73  |.tempr,tempi,wrs|
00000400  2c 77 69 73 0d 02 e4 04  0d 02 ee 06 ea 73 0d 02  |,wis.........s..|
00000410  f8 04 0d 03 02 0c 6e 25  3d 32 2a 6e 6e 25 0d 03  |......n%=2*nn%..|
00000420  0c 08 6a 25 3d 31 0d 03  16 35 f4 20 54 68 69 73  |..j%=1...5. This|
00000430  20 69 73 20 74 68 65 20  62 69 74 2d 72 65 76 65  | is the bit-reve|
00000440  72 73 61 6c 20 73 65 63  74 69 6f 6e 20 6f 66 20  |rsal section of |
00000450  74 68 65 20 72 6f 75 74  69 6e 65 0d 03 20 11 e3  |the routine.. ..|
00000460  20 69 69 25 3d 31 20 b8  20 6e 6e 25 0d 03 2a 10  | ii%=1 . nn%..*.|
00000470  20 20 69 25 3d 32 2a 69  69 25 2d 31 0d 03 34 11  |  i%=2*ii%-1..4.|
00000480  20 20 e7 20 6a 25 20 3e  20 69 25 20 8c 0d 03 3e  |  . j% > i% ...>|
00000490  36 20 20 20 20 f4 20 53  77 61 70 20 74 68 65 20  |6    . Swap the |
000004a0  72 65 61 6c 20 61 6e 64  20 69 6d 61 72 69 6e 61  |real and imarina|
000004b0  72 79 20 70 61 72 74 73  20 6f 66 20 74 68 65 20  |ry parts of the |
000004c0  74 77 6f 0d 03 48 37 20  20 20 20 74 65 6d 70 72  |two..H7    tempr|
000004d0  20 20 20 20 20 20 3d 20  64 61 74 61 28 6a 25 29  |      = data(j%)|
000004e0  20 3a 20 74 65 6d 70 69  20 20 20 20 20 20 3d 20  | : tempi      = |
000004f0  64 61 74 61 28 6a 25 2b  31 29 0d 03 52 37 20 20  |data(j%+1)..R7  |
00000500  20 20 64 61 74 61 28 6a  25 29 20 20 20 3d 20 64  |  data(j%)   = d|
00000510  61 74 61 28 69 25 29 20  3a 20 64 61 74 61 28 6a  |ata(i%) : data(j|
00000520  25 2b 31 29 20 3d 20 64  61 74 61 28 69 25 2b 31  |%+1) = data(i%+1|
00000530  29 0d 03 5c 32 20 20 20  20 64 61 74 61 28 69 25  |)..\2    data(i%|
00000540  29 20 20 20 3d 20 74 65  6d 70 72 20 20 20 20 3a  |)   = tempr    :|
00000550  20 64 61 74 61 28 69 25  2b 31 29 20 3d 20 74 65  | data(i%+1) = te|
00000560  6d 70 69 0d 03 66 07 20  20 cd 0d 03 70 11 20 20  |mpi..f.  ...p.  |
00000570  6d 25 20 3d 20 6e 25 20  81 20 32 0d 03 7a 1a 20  |m% = n% . 2..z. |
00000580  20 c8 95 20 28 6d 25 3e  3d 32 29 20 80 20 28 6a  | .. (m%>=2) . (j|
00000590  25 3e 6d 25 29 0d 03 84  14 20 20 20 20 6a 25 20  |%>m%)....    j% |
000005a0  3d 20 6a 25 20 2d 20 6d  25 0d 03 8e 13 20 20 20  |= j% - m%....   |
000005b0  20 6d 25 20 3d 20 6d 25  20 81 20 32 0d 03 98 07  | m% = m% . 2....|
000005c0  20 20 ce 0d 03 a2 12 20  20 6a 25 20 3d 20 6a 25  |  .....  j% = j%|
000005d0  20 2b 20 6d 25 0d 03 ac  05 ed 0d 03 b6 04 0d 03  | + m%...........|
000005e0  c0 0d 6d 6d 61 78 25 20  3d 20 32 0d 03 ca 04 0d  |..mmax% = 2.....|
000005f0  03 d4 2a f4 20 4e 6f 77  20 74 68 65 20 44 61 6e  |..*. Now the Dan|
00000600  69 65 6c 73 6f 6e 2d 4c  61 6e 63 7a 6f 73 20 72  |ielson-Lanczos r|
00000610  65 63 75 72 72 65 6e 63  65 0d 03 de 04 0d 03 e8  |ecurrence.......|
00000620  11 c8 95 20 6e 25 20 3e  20 6d 6d 61 78 25 0d 03  |... n% > mmax%..|
00000630  f2 16 20 20 69 73 74 65  70 25 20 3d 20 32 2a 6d  |..  istep% = 2*m|
00000640  6d 61 78 25 0d 03 fc 25  20 20 74 68 65 74 61 20  |max%...%  theta |
00000650  20 3d 20 32 20 2a 20 af  20 2f 20 28 69 73 69 67  | = 2 * . / (isig|
00000660  6e 25 2a 6d 6d 61 78 25  29 0d 04 06 04 0d 04 10  |n%*mmax%).......|
00000670  1a 20 20 20 20 20 20 20  20 73 3d b5 28 30 2e 35  |.        s=.(0.5|
00000680  2a 74 68 65 74 61 29 0d  04 1a 19 20 20 77 70 72  |*theta)....  wpr|
00000690  20 20 20 20 3d 20 2d 32  2e 30 20 2a 20 73 2a 73  |    = -2.0 * s*s|
000006a0  0d 04 24 17 20 20 77 70  69 20 20 20 20 3d 20 b5  |..$.  wpi    = .|
000006b0  28 74 68 65 74 61 29 0d  04 2e 04 0d 04 38 0e 20  |(theta)......8. |
000006c0  20 77 72 20 3d 20 31 2e  30 0d 04 42 0e 20 20 77  | wr = 1.0..B.  w|
000006d0  69 20 3d 20 30 2e 30 0d  04 4c 1b 20 20 e3 20 69  |i = 0.0..L.  . i|
000006e0  69 25 20 3d 20 31 20 b8  20 6d 6d 61 78 25 20 81  |i% = 1 . mmax% .|
000006f0  20 32 0d 04 56 14 20 20  20 20 6d 25 20 3d 20 32  | 2..V.    m% = 2|
00000700  2a 69 69 25 2d 31 0d 04  60 10 20 20 20 20 77 72  |*ii%-1..`.    wr|
00000710  73 20 3d 20 77 72 0d 04  6a 10 20 20 20 20 77 69  |s = wr..j.    wi|
00000720  73 20 3d 20 77 69 0d 04  74 22 20 20 20 20 e3 20  |s = wi..t"    . |
00000730  6a 6a 25 3d 30 20 b8 20  28 6e 25 2d 6d 25 29 20  |jj%=0 . (n%-m%) |
00000740  81 20 69 73 74 65 70 25  0d 04 7e 1e 20 20 20 20  |. istep%..~.    |
00000750  20 20 69 25 20 3d 20 6d  25 20 2b 20 6a 6a 25 2a  |  i% = m% + jj%*|
00000760  69 73 74 65 70 25 0d 04  88 17 20 20 20 20 20 20  |istep%....      |
00000770  6a 25 20 3d 20 69 25 2b  6d 6d 61 78 25 0d 04 92  |j% = i%+mmax%...|
00000780  31 20 20 20 20 20 20 74  65 6d 70 72 20 3d 20 77  |1      tempr = w|
00000790  72 73 2a 64 61 74 61 28  6a 25 29 20 20 20 2d 20  |rs*data(j%)   - |
000007a0  77 69 73 2a 64 61 74 61  28 6a 25 2b 31 29 0d 04  |wis*data(j%+1)..|
000007b0  9c 2f 20 20 20 20 20 20  74 65 6d 70 69 20 3d 20  |./      tempi = |
000007c0  77 72 73 2a 64 61 74 61  28 6a 25 2b 31 29 20 2b  |wrs*data(j%+1) +|
000007d0  20 77 69 73 2a 64 61 74  61 28 6a 25 29 0d 04 a6  | wis*data(j%)...|
000007e0  29 20 20 20 20 20 20 64  61 74 61 28 6a 25 29 20  |)      data(j%) |
000007f0  20 20 3d 20 64 61 74 61  28 69 25 29 20 20 20 2d  |  = data(i%)   -|
00000800  20 74 65 6d 70 72 0d 04  b0 29 20 20 20 20 20 20  | tempr...)      |
00000810  64 61 74 61 28 6a 25 2b  31 29 20 3d 20 64 61 74  |data(j%+1) = dat|
00000820  61 28 69 25 2b 31 29 20  2d 20 74 65 6d 70 69 0d  |a(i%+1) - tempi.|
00000830  04 ba 29 20 20 20 20 20  20 64 61 74 61 28 69 25  |..)      data(i%|
00000840  29 20 20 20 3d 20 64 61  74 61 28 69 25 29 20 20  |)   = data(i%)  |
00000850  20 2b 20 74 65 6d 70 72  0d 04 c4 29 20 20 20 20  | + tempr...)    |
00000860  20 20 64 61 74 61 28 69  25 2b 31 29 20 3d 20 64  |  data(i%+1) = d|
00000870  61 74 61 28 69 25 2b 31  29 20 2b 20 74 65 6d 70  |ata(i%+1) + temp|
00000880  69 0d 04 ce 09 20 20 20  20 ed 0d 04 d8 12 20 20  |i....    .....  |
00000890  20 20 77 74 65 6d 70 20  3d 20 77 72 0d 04 e2 24  |  wtemp = wr...$|
000008a0  20 20 20 20 77 72 20 3d  20 77 72 2a 77 70 72 20  |    wr = wr*wpr |
000008b0  2d 20 20 20 20 77 69 2a  77 70 69 20 2b 20 77 72  |-    wi*wpi + wr|
000008c0  0d 04 ec 24 20 20 20 20  77 69 20 3d 20 77 69 2a  |...$    wi = wi*|
000008d0  77 70 72 20 2b 20 77 74  65 6d 70 2a 77 70 69 20  |wpr + wtemp*wpi |
000008e0  20 2b 77 69 0d 04 f6 07  20 20 ed 0d 05 00 14 20  | +wi....  ..... |
000008f0  20 6d 6d 61 78 25 20 3d  20 69 73 74 65 70 25 0d  | mmax% = istep%.|
00000900  05 0a 05 ce 0d 05 14 05  e1 0d 05 1e 04 0d 05 28  |...............(|
00000910  04 0d 05 32 04 0d 05 3c  04 0d 05 46 04 0d 05 50  |...2...<...F...P|
00000920  04 0d 05 5a 04 0d 05 64  04 0d 05 6e 04 0d 05 78  |...Z...d...n...x|
00000930  04 0d 05 82 04 0d 05 8c  04 0d 05 96 04 0d 05 a0  |................|
00000940  04 0d 05 aa 04 0d 05 b4  04 0d ff                 |...........|
0000094b

Our privacy policy is here.
Last updated: 10 Jun 2024, 13:20 UTC
Website designed by .