User:Tom Allen/ta/AddendumX04P08B

From J Wiki
Jump to navigation Jump to search



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 ...


SX04P08B00C.jpg


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'


JX04P08B00C.jpg
... 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'


JX04P08B01D.jpg
... 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'


JX04P08B02C.jpg
... (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'


JX04P08B03D.jpg
... 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'


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