User:Tom Allen/ta/AddendumX04P08B
< User:Tom Allen | ta
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 ...
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'
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'
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'
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'