User:Tom Allen/ta/AddendumX04P08B

---

FOR THIS MODEL

---

---

... geodesics calculated using -@E1cn do not exist in the real world ...
... where geodesics are calculated using -@E1cn the containment of neutral particles is not supported ...
... only fully ionized plasma can exist where the containment of neutral particles is not supported ...

---

load 'plot'

pd 'reset'

uQ0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_)

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.00p1 0.00p1 0"_)*10"_

pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04A zfyc"_ 1))]) pTtotal xCpts''
pd (uQ0(cx0"_ 1;[-@(1{])"1@E1cn xRef@([hkyX04A zfyc"_ 1))]) pTtotal xCpts''
pd (uQ0(cx0"_ 1;[-@(2{])"1@E1cn xRef@([hkyX04A zfyc"_ 1))]) pTtotal xCpts''
pd (uQ0(cx0"_ 1;[-@(3{])"1@E1cn xRef@([hkyX04A zfyc"_ 1))]) pTtotal xCpts''

pd (((0.9675"_)@]"0);]"0) pTsteps _ _7 9 160
pd (((1.2832"_)@]"0);]"0) pTsteps _ _7 9 160

pd 'pdf'

... divergence of spacetime in x-coordinates (synchronized in 3D space at initial condition time=0) ...

uQ0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_)

vSgDd2=:uQ0([([-@(0{])@E1cn yRef@hkyX04A),:@(S03yH{]))]

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

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

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

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,620"_)*InitInt
pdDL     =:vSgD@vS''

pdDL arrayPut   'pdDL'
     arrayClose 'pdDL'

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

                arrayOpen  'pdDL'
pdDL=:79361 8 4 arrayGet   'pdDL'   NB. ... shape required ...
                arrayClose 'pdDL'

load 'plot'

pd 'reset'

pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDL

pd ((0.9675*cos@]"0);0.9675*sin@]"0) pTsteps _ 0p1 2p1 160
pd ((1.2832*cos@]"0);1.2832*sin@]"0) pTsteps _ 0p1 2p1 160

pd 'aspect 1'
pd 'pdf'

... geodesic (does not support the containment of neutral particles) ...

load 'plot'

pd 'reset'

incy0=:(2-~/\])@:(cy0"1)
incy1=:(2-~/\])@:(cy1"1)
incy2=:(2-~/\])@:(cy2"1)
incy3=:(2-~/\])@:(cy3"1)
veloc=:(((incy0^2:)+incy1^2:)^0.5"_)%incy3

pd (}:@:(cy3"1);veloc)@:((S00yL{])"2) pdDL

pd 'pdf'

... (velocity along geodesic)>1 and not supported ...

uQ0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_)

vSgDd2=:uQ0([([-@(0{])@E1cn yRef@hkyX04A),:@(S03yH{]))]

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

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

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

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,600"_)*InitInt
pdDM     =:vSgD@vS''

pdDM arrayPut   'pdDM'
     arrayClose 'pdDM'

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

                arrayOpen  'pdDM'
pdDM=:76801 8 4 arrayGet   'pdDM'   NB. ... shape required ...
                arrayClose 'pdDM'

load 'plot'

pd 'reset'

pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDM

pd ((0.9675*cos@]"0);0.9675*sin@]"0) pTsteps _ 0p1 2p1 160
pd ((1.2832*cos@]"0);1.2832*sin@]"0) pTsteps _ 0p1 2p1 160

pd 'aspect 1'
pd 'pdf'

... geodesic (does not support the containment of neutral particles) ...

load 'plot'

pd 'reset'

incy0=:(2-~/\])@:(cy0"1)
incy1=:(2-~/\])@:(cy1"1)
incy2=:(2-~/\])@:(cy2"1)
incy3=:(2-~/\])@:(cy3"1)
veloc=:(((incy0^2:)+incy1^2:)^0.5"_)%incy3

pd (}:@:(cy3"1);veloc)@:((S00yL{])"2) pdDM

pd 'pdf'

... (velocity along geodesic)>1 and not supported ...