Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://num-meth.srcc.msu.ru/english/zhurnal/tom_2001/pdf/art3_6.pdf
Äàòà èçìåíåíèÿ: Mon Dec 16 17:35:59 2002
Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 23:24:00 2012
Êîäèðîâêà: koi8-r
¸ e ë F8 8 "ù 8 e 0 C÷ Cù d Ph 2n " ø |8 r 9 6 7 z z 2ò r ñ 24 "r g " i n "f d " Fd h3 "8 2 r n i f á 2n 8 h e ¨ q V l V " ¸ e R " " Vd " h

ç

2 u v 8 Pu i

æ " V

H n i û ü â df "

l "

" Pu á n i P

æ k "ó i g

f

¸ f "

k Pö n " F d

F m "

F m d "

xý "8 ð i

"¤ ÿ ¥ ¥ ¶ ï ¥ õ ½ l d " g Cë ôå 7¤ ä wq P© d P " d R " d "° · ¶ 5° » «

"¤ ÿ ¥ ï "î í ½ i " d Hj n " g " " i V " r ¸ e " d V " l " i w "ì 8 g Hë " l " i ¤ õ f ê ´ ¤ ä 7¤ "r d H ä ° ¤ ¶° " R V g " ° "¼ ¶ 5° » ° g " " d " i V " Vf g F m " d q " á e f 2 "à 8 `9 "þ R° ¿ ¨ ¤ ¶ ®° e e i q R k © " d " H P ¶°

i " Pi e "k P´ " R f w ö C8 ¾ Rl k d " i " d Vj P e Rk H¾ " H " d F " « ° ½

4 e F m " 2d s 2e " d R " " Hf g g o g " d j e g " d " á 2e l " Hd R « ° ¤ ¶ V° é ¤ èõ ã f ¶ ®° e " fj «« » " Pd R° ¿ ¨ ¤ Cd ¶ 5° » «

" " " " i «« »

q P j R g " " d l F i f g d V V 2 g l r ¬ " i f g d V " h d h He g " f H¹ º f ¡ ¢ ²² "f n " "r · C¶ n " Hl 2j q F ¨ q H 2n l R " r " " l 9 " d " d R l P d ¸ e ¦ ¥ ¤ ¢ l " µ ²² ¡ ¢ ²P² ¸ 0 l " d ® ¬ f

n " Rl 2j o k q s e F m ¸ "8 g " d F " Rd " l " i H f " d P 7 ²² w¹ º ¶ c¸

g " d R i d " i " d Fj l P 2j w8 ´ ¸ e " d " ¸ e ¸ e q o R g Rf g g " d j s e 8 7 ± ° ¡ ¢ ° " r 2e f m " 2d H g g g F m " d q " ¯ 8 § ¨ ¥ ¨ ¤ § ¬ " « « « 9 ª ©

k " o F p ¸ e R o 2 P 2 q e d

i " Ri g g f " d " d " ¸i ¡ ¢ " h d h Fe q i g 2 m " d q " P Ri r w| d " " d " f l

R

g

n P " H w Cu o v Cq



k k k

e i h "r 2j n " " i l P F m ¸ Ve d " 8 e " " j P " d R ¸e d " f "f d ¸ h " " e " Vd f l d n " " ¸ q d R " Cd h e ¸e e " d e F n 2 q n ¸ ¸e g F

2j j h d R " e " e " d m

F H " d d P

Vi R C j d "

e 2 H 2j ¸ d i n " j " d

r ¸ e " 5 i

e " e R R " d

n h

hf k d R g

h 2j q C

8 " n

e " ¸ Hi Rd " Re " d

¸ lg " h ¸ e 2

d 2 g q ¸ C n g

d i F o F 8 ¸

" n g p 2 F m f

¸ ¸j " "u } l d

e i Cf e

" "f l " } d V

l e d ~ "

C Cf i " R i H{ P g

¸ ¸ e l j d R 77 j "f

e " P e Fq F l f g 7 z y u d "

n " p d x

i ¸ d 2 Ht j f

e 8 " F " F p

¸ e " r 9 " Vd H n " d ¸ i " o 9 T

" e

" d d e 8 wy w

Pf R j r e f V " " Hd x

s g ni ¸ Vi 8 y "y

F ¸ e F Cf 6 z H " " " g

g m g g r l i X

F 1 f e g " d C

h d V n

8 "u } ~ } w| V{ 7 7 7 Cz y u x w v Cu Rt " i " s e l F V g " f l "r g h l " d g " ¸ ¸e e " d " d j 2 j "f d " q j F p o " 2d F i " n " d " h e g F m f h l 2j h k " d V i d " i " d Vj 2j Vi i h F g f ¸ e d 2 2 " 2 2 " wy Cx wv "u Ht ¸s Rg X Vg D q rq p i PD CY hA Rg RX PQ FU HW PY HB RE HE VB fe Hd `B cb Ha VY PX B G VD PD Y `W PG PB HI VY RX VQ PW HB 2E HG VU 5Q TS RG CQ PI I HG FE CD CB A @ 9 8 7 6 54 3 21 0 )( $ ¸# ' # & % % $ # ¨ ¤ ! "¤ ! ¤ ¡ § ¨ ¸¨ ¡ © ¦ ¨ § ¤ ¦ ¥ ¤ ¸¢ ¡


ó k


$

He

d



V i e H "





i n ¸ e " d ¸ e " d 2 " d l " r h d h ¸ e z 8 d " F g 8 & B · 2 î 2 r ' P

nd " ¸l hq h i P U9 » ª IHG4 &

¸ " i R l H V h "9 r
1

e " 8 Hz " i B·

" 2 e C 2 h 2j

2 " h e

"f d " d

d Ci F p h Pe g d ¸ e

g " " "

l Vd

8 r

i P 2 g " i f q F g " ¸ e g 1 Vf "

H d " ¨

n ¸ e ¸ r9

Pf h

W l f

" d R el V W n ¸ e " d g V R R R

" ¸ g

l h d ¸e

d "

" ¸ @

" f

$ F

9



p

Hd " R " d " i k f " d `Y ó H æ 2 SX n h i g p g " " r " h d h 8 d " F " " f " d d H j " V f n n Pf î (1 e ¸ e d R î S@ r · "$ 2¸ ¥ r 9 & 7 ¨7

" d " i H ' ' n h f V ¸ "

« ¥ ¸ ¥ P I H R 4 e ¸ e d R r · F$ F» ¥ g 2f ¥ E ¥ D V¥ @ " d " i k f " Hd " "f $ 9 « « 8 î 2 î 1 î @ î 0@ d " 89 « 7 î$

" "f d ¸ e " d " V " i l " d · F$ ¥ P I H G i r d H f m " Vi g i " i " d e " n P o ¸ e

i V V " d F» ¥ R 4 &î j " V Rf ¸ i

" Fl " l P i h ¥ ¥ P

C D ¸ e f

F

d

" P

r 2 F " d 8 2 2 d

d n p o d ¸ T '

PIHR4

" d V " ¸ e ¸ e f m " i R g 2 h ¸ ñ " "f R " " d i d l Hl dl V " " h 2 f h h ¸ ¡ ¡ ò ê 4 ê 6 z

C n H $ 09 7 V

¸ d H n ' C

V " ¸ e

" e "

l R " d d

g R l R

" R

& î 2 h ¥ QE ¥ V 2 V d Ri g " l ª « ª $ 9 ¥ CB A $ « 9 7 ª R ¥ P i i V P d " i " d 8 8 87 n h e R " d "° 4 ¸ ¤ 2 & 61

¸ ¤ 1

¤$ ¸ e ¤

i

i

di ¸ " d

d " « l " "

"° 4 ° " d i " 2l " "f d i F p m d f r " d " i " d V i l " d h Hi n " " R P i " H n " d e ® h " R " d " i k f " d n k 2 r h " F " d " i k h e " Pd P l á e f " q " e

o q ¸ d " "f o n g " k Hq Vd " P g

i " h d

P V 5e " f "

" d Pl g Pj

Ri f F

¼¤2 d " ¸ e g mj p o " i

i i " f l ¸

" F " Vl g F F 2 H r p g n d i d
$

d V e "

l " l o "1 F f i

il ° ¼¤ d 8i l m "

d

e k © r V n q " R d "

" P 3 2j i "

f d " r¤ F m ¸ d 8

h

2

d ¤ 1 Hi

p d Pd Cf " P F wy

d

g1

" d 2 d o Rq " d "¼ % f dq e " j "

2n i k g g " 2k

d " " % V

dl ¸ g

i r 2j "d r 2 " d "f d f

z V " " w8

h9 8 " h g i #

l f k ¤ wx

Pd " ¸ d " h " 2 ¸ e f q F PD

9 á e g 2 l Fd Fd " g y h " Fi "

V F u ¤

F ¸ ¸e j " m l r "

" m e i ¨¥ ¥

sf " d r " f " "!

" d ¸ e V V d

d d C g § t

" j d " i

n i f " d

" " 8 "r j 7 t

wd ö g " µ 8 " "t ¨

9 s " ¨§ wy

C ` "f d " i ¸ d h e ¤¦ F

2 e ¤ wx

" "

w Hq i i " H ¸

d 5d

" r "u P ef ¤

" " f e r

d ¸ w| P

"
¤¸

f } P ¸ ¸q w © §

2 5 " r P9 i 0$ ¸ Ce e " d " j 2j n d ~ } H n o ¸ Ve 2 2 w u wy

d ¸ C i " e C{ e f ©

F " 3 «« " i e " d ¸ e " n Hf i " g 77 " r° d " rl o ¤ u

d d ¸ " d 7 z d i i

xw Pe ) ° " sd wx

e R

j He ¸ n 2 g ql Pt F ' H° 2j " l H

g k i " i g k " kf

g P ¸ F Cd i

H m sf h f

Rf

j r F ¸e " r f H e



e

q

"

p



" j

j

f e h m q

F h 2 i ¸ e

n " e " f " g "

n Fe " d n g f

q " i ¸ Ve q H " d " f " C

d R d

F d " h d " n F p

l n " e ¸ q H

Vl j F p g f 2j V i f "r f " ¸ n " f d 2j ¸ e H g " d " w "r " C ¸ e " d hi ó

" V Vi F p o f 2 C " 2n r f r 6 r z P F m f " d Ch 2j q C g wf "r g " d l i " F m ¸ e l F Hf d d ¸ e Hf i " R " F m g F m f g " d f ¸ Ce g f e " " d "f " ¸ 2e " l " d q H " ¸ Fe q C h e Vl " "8 Re j i h e " h Ve g l f " d w8 ´ ¸ e Cl i " á n h i 2 g o g g ¤

f n h ¸ ¸e Vl l g õ

2n f " e " " g h " o

"f ne 8 d e " dd d e d F zn " ¸ r ¸

e Rl " d e 6 g

¸ d Ce g ë

l e o f " h " d

s " d d P e " "

d Pj C " R i Hq

"

i 9 7 d e d " ¸ yu h %° r i F Vy ¤ 8 V h e " i d " o Hd " i 2j g l " i

j i " 8 ¸ H C

© 2 k H© e ¸ ¸ e " d F Pi 2e v Cu g % " V i x "y ¸ e H h " 5d R e e "f P "r h "



" R Fl

o ¸ i d " Pi " m "

" i e " " " m f f © ¨§ " " " V h "

¸ 2 " F " f " i r jf ¤¦ n q Pf n l Re

" h " jf " d " 2e r " j¡

f d " h ¸ e

Hi ¤ H n 2 g " i Rl n d ¸ " % (' d ¥ i 2 " Hd " e q " l H g " g e F p

H 1 h Hf ¸ H e H d & " ¸ Ht r Fi d F dl w f d F q



F r¤ 9 " V e il 2 i % h Ve g¢ ¸ e H n f ¸ e d o " d e Hi q " mg `

¸i e " d d " d
(

$ ¸# ' # & % % $ # ¨ ¤ ! "¤ ! F ¤ ¡ § ¨ F ¸¨ ¡ © ¦ ¨ § ¤ ¦ ¥ ¤ ¸¢ ¡


$ k

P ¸¡ jil " Pd ¨V 8 l

¸ e q r n $ P ¢¡ P

d ¸

l " e F m

r i " d ¸

¸ " d R e

e C n

" d 7 d "

i l " d $ i "

¸ P d P

e ¸ ¸ e " g

2e h d

F H

p g " i ¸ e ql

" ¸ " 5d

"f e H

d ¸

Ci e

" Vd g i

l

q

" d r i

" ¸ e

i 2 P gd d

k g " "

f i i

f " " d " d

Pd R j n

n f

m

" " ¸ e

l " g

lr f

7 " g

$ d r "

" " Ri

¸ Vd

e " ¸ e

d R i ¸

¸ e ¸ e f ¸ e el " i

¸ e " d á ¸ e ¸ e F p g " " Rd H i " n " d " h Ve g F m R F p i " F " d " i k f " Hd " "f d Ri g " l g p o p " d l "r g h i P V 8 z 8 w
s si , 2 ¬ R01 ,
0

9 P "W n h i 2 V g h e p " d l r " h d h e h i 2 P ¸ e d µ 8 d " i " d R i l " d R j " i " H 2j i H i " o F m H j F p o r ¸ e " d " ¸ e ¸ He g " Vl g p o p " d l "r g h i P V 8 9 8 w
s si+ 1 , 0 ¬ ti+
1


(

$ ¸# ' # & % % $ # ¨ ¤ ! "¤ ! ¤ ¡ § ¨ ¸¨ ¡ © ¦ ¨ § ¤ ¦ ¥ ¤ ¸¢ ¡

Lt ( L , Tf , Nt) :=

Times ( L) :=

Ô s0 , 1 Æ R0 ¬ Ã Â Õ s0 , 2 Ü

s0 , 2 ¬ L0 ,

s0 , 1 ¬ L0 ,

s0 , 0 ¬ 0

k¬1

DT ¬ 0

Dt ¬

for i í 1 .. Nt

otherwise k¬k+1

si , 0 ¬ Dt â i

if si , 0 ¸ Lk ,

Tf Nt

si , 1 ¬ R00 , 0 + Lk-1 , 2 â si , 0 - si-1 , 0 - DT

si , 1 ¬ R00 , R0 ¬
1 0

si , 2 ¬ R01 , 0 + Lk-1 , 3 â

DT ¬ si , 0 - Lk-1 ,

v0 ¬ Ã

s0 , 0 ¬ 0

t0 ¬ 0

for i í 0 .. rows ( L) - 2

ti+ 1 ¬ ti +

r¬Ã

ÔL Ã ÃL Õ

5

ÕL

ÔL

k-1 , 0 + Lk-1 , 2 â DT

k - 1 , 1 + L k- 1 , 3

ÕL

ÔL

0

0,3

0,2

i+ 1 , 0

i+ 1 , 1

Æ Â Ü

5

r v0

Æ ÔL Â-Ã Ü ÕL

( (si

Æ Â â DT Â Ü
, 0 - si-1 , 0 - DT

i, 1

i, 0

Æ Â Ü

) )


j d

" "f d " P 8 ¸ e " d ¸ e l " k " R f " k " Rd 2n " ¸ e k V n ¸ he 2j

si g ¸ ¸ Ve g

2e l

¸ R

e 2j f ¸

Rj " ® " Hl "r " " e " d f $ P
35 0

wl r ¸ g

d e " "

¸ e " H 2j e ¸

d V

" n "f F

á ¸ i l " Pd d " m ¸ He n « @ 7 g " il

e " ¸ q

e " d H g

2j f

f V

j

R 2j F 1 8 e r7 " "

i S$ V 2j

H 2j j g Pl

¸ d

Rf g e "

f " j r

2

7 " d q r n l "

P i

k "f d d " í

¡



F p g "î í g 2j 7 2¸ i " i r68

o d " "f n r "î " " " d

í Rl R 2j

i " d j "r f g d 2j ¸ e " d " d i H l " i d " d P d " d i 2j l V l Ff ¯ n " ál e f " d V " V 2 Hq F m ¸ e
q " i h ¯ "8 " " d f i j F ¸ e

Rf g Cf " d k í i h á o " F 2 à
©¨ "î í

9 P R i g P n

l f " " 9 7 h V Vq d "

"

¡

2j l

q

R " d

" "

" " f

" d d " d

R d " 2

r

" " h Fe Fq $ P ¸¡



F m " h ¸

Pf f 7 e

g

9 l

l i h l

H 8$

"

" l ¸ e " d j g µ

¡

¸ e Rl g 2 V g Rf " ¸ e i Pf ¸ e " Pd V 8 6 8 w
0 0 50 10 0 15 0

k r k k



e Fe j "



Pd

F



d f

m 9

¸

q p " d d d

"

d "f d p o m P m

7 0$ i ò ò W e " i f r H v d ¸ j " i F

f i " g i " k

m " d " h C k " d H ' "

¸ Re g Re ¸ ¸ d h © Vl d i i P h e k " R ó Hd

g " " j 2 " i F il H " H n

" "f F " p " d n n F

" n d " d "

g Ri i f f F " " R

" Cl m d 2 ò i d " g H " i h Pe n

p ¨q r d H d d "

n f g ò n i k

F p r g W r 4 4 " " l f d P

m h " W ò F " d i

f d h l l ó i " C C X r 4 p d "

8 r d " e " d C d P n g ¨X "ó df ¡ ó i " d sd e "

" g i P V

§¦ v 4 d " d d n " H "6 r " " e " h R 8

¸ " d H ¸ 4 ói r i l " d " R r e " d Vd "

e ¥¤ ó v W Fi " H n i h n 8

e " d ¸ e 2 i ò " " C

e ¸ Re d d " Vd " W r W ó h n " n " F " H 4 He " P d f

P ¸ d d 2j " d i d

¸ ¸ e e " i " "r

i Fe " " V Re " H

e ¸ d l h " H k g

l F n ¸ 2 d n i f

e "f " " "

$ ¸# ' # & % % $ # ¨ ¤ ! "¤ ! F ¤ ¡ § ¨ F ¸¨ ¡ © ¦ ¨ § ¤ ¦ ¥ ¤ ¸¢ ¡

R0 := Ã

r ¸ g " d ¸ e d i " h Pe " H n " f

Tf :=

K := 133

D := Ã

Dt :=

Dti

Ô Õ

Ô -3.3 Æ Â Õ 3.63 Ü

0.01

10

15

55À5ß )rows( L55) (L66À5ß )rows(L66)-

(L

0.0

s

5

for i í 0 .. Nt

2 si ¬ È (Lt55À1ß )i - (Lt66À1ß )iÛ + È (Lt55À2ß )i - (Lt66À2ß )iÛ

Æ Â Ü

L55 := augment( L5 , T5 )

Lt55 := Lt ( L55 , Tf , Nt)

T5 := Times ( L5)

v := Ã

R01 := R0 + D

Ê

Ô 1.5 Æ Â Õ -1.8 Ü

1

1

otherwise

Àß Àß if (L55 5 )rows( L55) -1 ¸ (L66 5 )rows( L66) -

Nt := 2000

L5 := L4 ( R0 , v , Bound , Bound2 , K)

L6 := L4 ( R01 , v , Bound , Bound2 , K)

T6 := Times ( L6)

i := 0 .. Nt

Ø

L66 := augment( L6 , T6 )

Lt66 := Lt ( L66 , Tf , Nt)

Ê

Lt55i , 0
2 00 2 50

Ø

2

1

3 00

e

i "
¸

e " g v Ci n " d " 7 " d

w r d " "

¢¡ ó h " d P n F p o

" R Rd " ó h He





¸ 2e " Hd " l " d ¸ e C n " "

P i Hl " " d n d "

2~ F i P

" i " X e d " ¸ d " "

" d " H ó p " d " P f H

d 2j ¸ m V P d d

" " f il d ee d

h d d " ¸ " i i

V i h " d P ((


¸¡

k " r n P n " H V V " " d "î í r " l " Pd 2n F " r ¸ e 2l g 2 H e F p o g l r e " ¸ e



i h " 2d

g ¨ "

r V " q P d V ¸ e

$ P h "

h $ V F i ¸f







d e R4 i
¸

r h

¸ e $ V $ P m "

¡ " d q " 2 2n

sl e ¸ e ¸ e P 2



$ V

l Hl g l Fl " w "

g 1 8 i " F 2j " d " h He n " " d i ¸ e " d " i V f m ¸ e R " d "r "f V i " i r 2n d ¸ e f m ¸ `e C l P F m ¸ e n d " i " d g q l " d " i d k Pi " i r 7 F$ ¸ e " d ¸ e ¸ e ¸ e i l q r ¸ e g d " i " d f m ¸ Ce g g " "8 l d " g r 2 1 P 7 ¨ q ¸ e Vl 2 V " " d "8 g F m " d q " i " R e Y G F4 R 2 Y R 4 r " 1 Y G U4 ¸ e l " d " i P g F m " d q " s h i " "8 ($ ¸ e " d ¸ e ¸ Ce ¸ e i l "q r ¸ e g d " i " d f m ¸ Pe g g " " Cd ¸ e Hf " r i V V

" " g " P g Hl

d m d n

n

g g R j ó H r 4 k F r n

F Ff " " i

p d " "

og i F 7 n

F p 77

Fg r C9 F "

od "f d F m d " i

d " " d

8 ¸ q q " g e d

"f d e 8 8 F m F

" " m "

q fg h R d q "

e " " d f R

" d d " i p d F

" g g

f o i

R F " " Hd e " Rf "
"W n " d d

Vi 2

g d V g

2$

2

" m i 8 ¸ e "

i



d F

V " 8 e i R h 2







e ¸

" r ¸ g d ¸ F " " " w g n Fq g V F "

"

n " d H " i " 6 p o ¸ Y z d ¸ e 9 g¨ ¸ d k "



d e d

e Ve

2j k "8 2k

2 0X n Hd " i " " d q " h i 2 f f F 8 m ¸ e l H g

æ 8 m F

H g f

h f C 2j 8 l g

i 2 C 2n d " r " d " g l " i " d $ 2 " d " g ¸ e

V ¸

F " d " i h se n g $ m ¸ ¸ e 2 F "

p " " e f 1 Y e n

d " " " d

l "f Rl s R i i

r d F d µ d 4 r " d ¸

i dq $ « 8 $ 2 " d i f

" g 8

æ 2 X " F p o « î$ 8 8 f 8 $ 2 G $ l d ¸

5 « i s 1 e P n ¸

H o 2 $ P Y n e

"ó p ¸ e ¤ 6$ r " d Y F4 s e d h

" g Rf e

l ¸ R ¸ " "

m " 8 h e " 4r e Cd

dq F ¸ d " i "r g F ® $ ¸ F

1 l e m

P R ò $ i Cf g " i P ¸ ¸ F m " Y " d á " d q " d

V 8 F g F V 8 e ¸ dq G l V "

g " Rd " d 1 ¸ e F p " d ¸ 2 " R 4 " d h " d R

" 2 2d d l " ¸ "

R " Y " m " e d

Ri G P q d s " e d H l

"f d 04 " r " l 2 ¸ e P

" r "f " P i hd ¸ f d sf g f " d "

" d " ¸ h e l 5 ¸ r

i i " " Cd d " F " 28 Rq R f " f " d " R e Fe " P V e j " d f j R e " R

q P g

" "

" g e j e ¸ H d dq d ¦ 8 Vl "ò 8 FW pg d o f i " r

F P " d ¸ e f d " n e g F 2 1 F " " ¥© F ¸ e " d F e " d " i C n C 2 j 2e l g 5

Vi " " g 2 h " d d ¸ h " d ` j "

e " df d i i i d ¸ e

q " F h e e i ¸ 2e d V 2 $ P
40 0

¸¡

¸ e Cl d i "8 8 w
4 1 .10



d

¦





g e 2n k " g

g R d ¸ e "k H F o g

Vg d H " ´ 8 ¸ Vg

i r §¦ f d e " " f g Pd

f i 2j

ñ á 2 r " C

¸ e l F P ¸ ¸

" V n i " " "r " 8 Rd ¥¸ ¤¸ ¢¡ y u x 2 l " d d d i

i X 2e e "

" wv C

¸ e C Cu g i ¸

" n Vt

" l ¸ e F

Rf ed m

" i " ¸ R " d e

o e i

f "

j " $ h á d F

d i 5

r 2

" F P ¡ o ¸ e ¸ e q j 8

d " d @

" d ¸ " F g " 8 P

e "f p 8 f

i

2j Hl

d V P j ¸ e

" " d P


(

$ ¸# ' # & % % $ # ¨ ¤ ! "¤ ! ¤ ¡ § ¨ ¸¨ ¡ © ¦ ¨ § ¤ ¦ ¥ ¤ ¸¢ ¡

Dtsi

3 1 .10

0.01

0.1

1

0

Dts := CFFT ( Dt)

2 00

i := 1 ..

i
600

rows ( Dts) - 1

2

800

10 0 0

8 f "

µ

i " " "

l e

i g " Pd " d

l h Hf F H H

l C e l " gl $ P ¡ 8 8 r C

F 8 rfg i 2 i r

2e h 8 h ¸ 2e d

d d n ¸

d i e "

d F i ¸ " d F d H




§¦

¥

ki 8 k g

" d " d C f "î

Pj 2j h ¸ "î í g

s í ¤¡

e "

F ¸ e Hl

m g



Ci i P C Vl " 2

h nf e F

R g áf F m p o g

¸ e d i P ¨¡ 8 ¸ "î

n í l

¢¡

d F p r e

V o "

2 " Rl ¸ e

" d " 2d

" e f " d " d d P " d P i V

f d C j ¸ e e ¸

2 " d r H i

"f d g Rf

" 2 "

C H

" Pq ¸ e i 2n

á " " d ¸ e " 8 ¸ i "

l " d " f " h e F m n F p o F m

k d Ri r $ V



Vi e ¸ i



¥ Ve P

" " d 8 2 k " d n F k 5f o d " i

C@ l " d q " r " k f d " i " d V 1 Y p o " d V

¸ $ P ¡ e g ¸ 8 " " d f g 2 8 R4 r2 R g 2

e l V g o h V n d i F 2 i " f 7 $ ¸ i ¸ e " d C ` o h o R " s e F p o g l r e " ¸ e " "ò 8 Vd " s 2e " f " d 1 Y G U4 H 2 Y R 4 r " C g " "f d " Pq g Pf " d e F p o g l r e " ¸ e " 8 r i h " d V " V " l " d q "

" e 2d r 1 " 2d $ H

" P ¸ " i Y 2 2 $ V

g ¸ ¸ e G4 V ¸ e $ P

h d



" Ve C w "8 Vf e ¸ F m " i f n "

¸¡



lr ¸ e F ¸ e " F m e dq " ¸ e ¸

¸ e "î e r " f 2n d F á " n " i 2n l Cl d í p d R " V e

o ¸ C q " Cl H

H e d " d

"g wy l " Hd 2j F i " l " i " " n i i F

w wy 8¦ f " n n " ¸ i l 8 n

x " H " R j " ¸ " F " "

m ¸ " f 2e " Fd "

i 2 " " h

g 2j " " d 2

¥© g F 21 ñ " H ¤¡ ò " d g h1 4 d Ce g ¸1

o i F " d " i g F m " d q " H g H g " d P R g 2 V g Rf " ¸ e i Pf ¸ e " Pd V 8 ò 8 w
( X , Y , Z) X := P1 À0ß P1 := P ( L6 , T6 , 1000) À1ß Y := P1
P ( L , Times , N) :=
à ÕÊ Ê ÕÞÏÄ ß ÏÞß ß ÅÕÁÞ×ÛÅÈ ÖÏÞÌÕÆÅÈ

" d " i g F m " d q " R g " Rl g p o p " d l "r g h i P ¨V 8 "4 8 w
V Vi , 3 ¬ L j-1 , 3

$ ¸# ' # & % % $ # ¨ ¤ ! "¤ ! F ¤ ¡ § ¨ F ¸¨ ¡ © ¦ ¨ § ¤ ¦ ¥ ¤ ¸¢ ¡

Z := P1

N1 ¬ rows ( Times ) - 1

T1 ¬ Times N1 , 0

DT ¬

j¬0

for i í 0 .. N - 1

t ¬ DT â ( i + 1)

Vi , 2 ¬ L j-1 , 2

Vi , 1 ¬ L j-1 , 1 + L j-1 , 3 â Dti

while t > Times j , 0

Dti ¬ t - Times j-1 , 0

Vi , 0 ¬ L j-1 , 0 + L j-1 , 2 â Dti

j¬j+1

À2ß

T1 N

( X 2 , Y 2 , Z 2)

é ÁÞÄ ÏÞÐÌØÈ ÀÕÊ Ê ÕÞÏÄ


(


9 7 7 z8 ò 78 ñ z o h i f " Rd V

f q q ¸ h D w4 z Y z G 4 9 F R 8 7 @ yw 8 4 w 6 8 a 4 e f e r f f e ¸ r g ¨q q s q 2 ¦n 2m 2m im 5l 2k yj V YV e 8 4 c @ D S 9 v e @ 9 a 8 9 @ "w W6 Q F w @ S @ c f { { r o { ¸ g ¨q q s q 2 ¦n 2m 2m im 5l 2k yj V iV ¥` A 8 S 6 4 6 e @ w 6 8 a ¥4 A F 9 4 YC 4 4 c @ D S 9 v 4 4 S G 9 @ c y yh

w

l { o s p p c e

Y ¸z F 6 @ F Q )4 4 S 6 @ 9 S 9 8 7 t R D t uS @ C F w 4 9 S G ¥F A 5S ¥ t w 4 f ¦4 h I t 9 e 8 9 4 7 8 u9 WE 8 4 9 8 C 8 ¦E X ¸0 p q ¸ h F w t 5F l Y z 2 A 8 S 6 4 w6 Q 4 w 6 8 a ¥4 A F 9 4 wC T S 6 8S D F 4 7 7 4 uR Q "4 5P I F E 4 8 6 6 F wD F E S 6 e @ E y6 8 4 w 6 8 a 4 C @ d ¦D ¨v ¥4 A G 4 9 "8 u F D S @ ¥4 A 4 d t F Q ¥8 iA e @ w 6 8 a 4 S 6 4 S F S h6 ¥` A 8 S 6 4 g6 e @ w 6 8 a ¥4 A F 9 4 HC e @ 9 C Y@ G 7 C ` f 8 S @ c 4 Hd e y x y y V V yx uE F 9 8 C 8 E 4 D c 4 T S F S W6 G 4 6 D 8 uE G F 9 9 @ D S w 8 7 b4 ` C @ S ¥8 BA 28 ` 9 T 7 8 S 4 7 6 4 2a ¦` YX V WV )U T S 6 5F "3 E @ C D F 4 7 7 4 HR Q "4 5P I F E 4 8 6 6 F HD D po @ v ut G4 c "t e w6 s 9F i6 r 8a sr E ¥` o 4 q @D A e ¸ C@ p 4 d be 78

| y ' & & 9a4S6 e r p ¸ 4 W6 8 4 s ¨q p o D I49F hg BA 8 @

¦0 @ q 6 r dv i C ¥@

F w6 q E ff 9

Q@ r 8 2 @ ¥e 98

¨ "( ¦S g a ¥4 ¦n 96 76

m hf A 2m @ 84 54

w0 ¥0 w0 ¥0 {q F 9 24 2m im R @ ¥ 9FE 23 20

5F "0 "0 q 5z "0 2k yj "0 @ ¥4 ¥A

S@ ¥0 ¥0 " wl j ¥0 5l ¥0 D 1 ¨0 )(

y0 y y ¸ h }' }' yx y V A ' & "%

Y F V V FD uV

% 4 % d $#

% % ¥( % ¥( t ~ } ¸| ~ } ¸| @C ¥( 49 ¥(

w ¨ ¨ 5F ¸ ¸ E ¨ w ¨ @D "!

¸ ¥ ¨ ¥ ¸ ¸ ¨ ¥ ¥¤

¥© ¨§ ¦¡ ¥¤ ¸¢ ¡

" "

í i l

kl k
(

F m

wq ¸ "

h r ¸ e Hd 2

H F eg ¸

r m d Cj ¸e

e F m ¸ he " i e " d

" $ P " "

Pd n ¸¡ d i d

i f j

¸ e m ¸ e l r

g " j "f

Pj ¸ d

"8 g " Vl g " "r "

F i 8 " " sd l " q

" d l V di i ¸ j F

F ¸ 8 p o

F e "u

" F }

m

g " i " sl ~ } " d

¸ e r n d F " r $ " w| { 7 7

V " P d 7 áz F p

g Cf i " ¸¡ ¸ e y o "

d ux i

R ¸ e " sd wv

f Hl g " d R d Cu t " l w r

ql r7 i

" d C " Vd $ " d " C j

F d 2 i

m ¸ e

f d ¸ " ¸

" d h H e f " q e "

8 ¸ d " d

e ¸ C l " f i

" d i V ¸g e H fg ¸d r

e " ¸ e d

k d " h "

$ ¸# ' # & % % $ # ¨ ¤ ! "¤ ! ¤ ¡ § ¨ ¸¨ ¡ © ¦ ¨ § ¤ ¦ ¥ ¤ ¸¢ ¡