User:Tom Allen/ta/AddendumX02P02

vZ0=:0.05 1 2 8 1 1 1 _1 0,:0"_

vSgDd2=:vZ0([([(0{])@E1cn zRef@([hkyX02C vZyc"_ 1)),:@(S03yH{]))]

Initx =:  0.7500 , 0.5p1 ,  0.25p1 , 0p1"_
Initdx=: (0      , 0     ,  0      , 0.1"_)*%@InitInt=:8"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initz =:vZ0 vZzc Inity
Initdz=:vZ0((vZgamma*(0{])+-@vZRv00*(3{])),(1{]),(2{]),vZgamma*(3{])+-@vZRv00*(0{]))Initdy

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Initz
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdz%Initds
vS06ds   =:Initds
vS07count=:(0,10000"_)*InitInt
pdDS1A  =:vSgD@vS''

pdDS1A arrayPut   'pdDS1A'
       arrayClose 'pdDS1A'

vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.01"_

vSgDd2=:vZ0([([(0{])@E1cn zRef@([hkyX02C vZyc"_ 1)),:@(S03yH{]))]

Initx =:  0.7500 , 0.5p1 ,  0.25p1 , 0p1"_
Initdx=: (0      , 0     ,  0      , 0.1"_)*%@InitInt=:8"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initz =:vZ0 vZzc Inity
Initdz=:vZ0((vZgamma*(0{])+-@vZRv00*(3{])),(1{]),(2{]),vZgamma*(3{])+-@vZRv00*(0{]))Initdy

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Initz
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdz%Initds
vS06ds   =:Initds
vS07count=:(0,10000"_)*InitInt
pdDS1B  =:vSgD@vS''

pdDS1B arrayPut   'pdDS1B'
       arrayClose 'pdDS1B'

vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.90"_

vSgDd2=:vZ0([([(0{])@E1cn zRef@([hkyX02C vZyc"_ 1)),:@(S03yH{]))]

Initx =:  0.7500 , 0.5p1 ,  0.25p1 , 0p1"_
Initdx=: (0      , 0     ,  0      , 0.1"_)*%@InitInt=:8"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initz =:vZ0 vZzc Inity
Initdz=:vZ0((vZgamma*(0{])+-@vZRv00*(3{])),(1{]),(2{]),vZgamma*(3{])+-@vZRv00*(0{]))Initdy

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Initz
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdz%Initds
vS06ds   =:Initds
vS07count=:(0,10000"_)*InitInt
pdDS1C  =:vSgD@vS''

pdDS1C arrayPut   'pdDS1C'
       arrayClose 'pdDS1C'

vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.95"_

vSgDd2=:vZ0([([(0{])@E1cn zRef@([hkyX02C vZyc"_ 1)),:@(S03yH{]))]

Initx =:  0.7500 , 0.5p1 ,  0.25p1 , 0p1"_
Initdx=: (0      , 0     ,  0      , 0.1"_)*%@InitInt=:8"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initz =:vZ0 vZzc Inity
Initdz=:vZ0((vZgamma*(0{])+-@vZRv00*(3{])),(1{]),(2{]),vZgamma*(3{])+-@vZRv00*(0{]))Initdy

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Initz
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdz%Initds
vS06ds   =:Initds
vS07count=:(0,10000"_)*InitInt
pdDS1D  =:vSgD@vS''

pdDS1D arrayPut   'pdDS1D'
       arrayClose 'pdDS1D'

NB. ... execute (ijx) ...

                  arrayOpen  'pdDS1A'
pdDS1A=:80001 8 4 arrayGet   'pdDS1A'   NB. ... shape required ...
                  arrayClose 'pdDS1A'

                  arrayOpen  'pdDS1B'
pdDS1B=:80001 8 4 arrayGet   'pdDS1B'   NB. ... shape required ...
                  arrayClose 'pdDS1B'

                  arrayOpen  'pdDS1C'
pdDS1C=:80001 8 4 arrayGet   'pdDS1C'   NB. ... shape required ...
                  arrayClose 'pdDS1C'

                  arrayOpen  'pdDS1D'
pdDS1D=:80001 8 4 arrayGet   'pdDS1D'   NB. ... shape required ...
                  arrayClose 'pdDS1D'

load 'plot'

pd 'reset'

vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.01"_

pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2))                 (i.@(8001"_){]) pdDS1A    NB. ... blue ...

pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2))                 (i.@(8001"_){]) pdDS1B    NB. ... red ...

pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@(vZ0 vZyc"_ 1])@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1B    NB. ... green overwrites blue ...

pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@(vZ0 vZzc"_ 1])@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1A    NB. ... purple overwrites red ...

pd 'pdf'

load 'plot'

pd 'reset'

incz0=:(2-~/\])@:(cz0"1)
incz1=:(2-~/\])@:(cz1"1)
incz2=:(2-~/\])@:(cz2"1)
incz3=:(2-~/\])@:(cz3"1)
veloc=:(((incz0^2:)+incz1^2:)^0.5"_)%incz3

pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1A
pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1B

pd 'pdf'

load 'plot'

pd 'reset'

pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2)) pdDS1A    NB. ... blue ...

pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2)) pdDS1C    NB. ... red ...

pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2)) pdDS1D    NB. ... green ...

vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.90"_

pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@(vZ0 vZzc"_ 1])@:((S00yL{])"2)) pdDS1A    NB. ... purple overwrites red ...

vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.95"_

pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@(vZ0 vZzc"_ 1])@:((S00yL{])"2)) pdDS1A    NB. ... pink overwrites green ...

pd 'pdf'

load 'plot'

pd 'reset'

incz0=:(2-~/\])@:(cz0"1)
incz1=:(2-~/\])@:(cz1"1)
incz2=:(2-~/\])@:(cz2"1)
incz3=:(2-~/\])@:(cz3"1)
veloc=:(((incz0^2:)+incz1^2:)^0.5"_)%incz3

pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) pdDS1A
pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) pdDS1C
pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) pdDS1D

pd 'pdf'

---

Next Page: User:Tom Allen/ta/AddendumX02P03
Prev Page: User:Tom Allen/ta/AddendumX02P01

---