User:Tom Allen/ta/AddendumX04P01
< User:Tom Allen | ta
Jump to navigation
Jump to search
load 'plot' pd 'reset' uQ0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_) xCpts=:_ 0.2 1.4 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@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.2 1.4 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.25p1 0.25p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.2 1.4 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.40p1 0.40p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.2 1.4 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.50p1 0.50p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.2 1.4 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.60p1 0.60p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.2 1.4 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.75p1 0.75p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.2 1.4 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 1.00p1 1.00p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' pd 'pdf'
... an appropriate scale (or gauge) must be assigned ...
load 'plot' pd 'reset' uQ0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_) xCpts=:_ 0.7 2.8 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@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.7 2.8 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.25p1 0.25p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.7 2.8 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.40p1 0.40p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.7 2.8 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.50p1 0.50p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.7 2.8 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.60p1 0.60p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.7 2.8 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.75p1 0.75p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' xCpts=:_ 0.7 2.8 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 1.00p1 1.00p1 0"_)*10"_ pd (uQ0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX04B zfyc"_ 1))]) pTtotal xCpts'' pd 'pdf'
... an appropriate scale (or gauge) must be assigned ...
uQ0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_) vSgDd2=:uQ0([-@([(0{])@E1cn yRef@hkyX04B),:@(S03yH{]))] Initx =: 1.01 , 0.5p1 , 0.25p1 , 0p1"_ Initdx=: (0 , 0 , 0.01 , 0.1"_)*%@InitInt=:32"_ 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,6280"_)*InitInt pdDK =:vSgD@vS'' pdDK arrayPut 'pdDK' arrayClose 'pdDK' Initx =: 0.53 , 0.5p1 , 0.25p1 , 0p1"_ Initdx=: (0 , 0 , 0.01 , 0.1"_)*%@InitInt=:32"_ 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,6280"_)*InitInt pdDL =:vSgD@vS'' pdDL arrayPut 'pdDL' arrayClose 'pdDL' Initx =: 2.400 , 0.5p1 , 0.25p1 , 0p1"_ Initdx=: (0 , 0 , 0.01 , 0.1"_)*%@InitInt=:4"_ 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,72000"_)*InitInt pdDM =:vSgD@vS'' pdDM arrayPut 'pdDM' arrayClose 'pdDM'
NB. ... execute (ijx) ... arrayOpen 'pdDK' pdDK=:200961 8 4 arrayGet 'pdDK' NB. ... shape required ... arrayClose 'pdDK' arrayOpen 'pdDL' pdDL=:200961 8 4 arrayGet 'pdDL' NB. ... shape required ... arrayClose 'pdDL' arrayOpen 'pdDM' pdDM=:288001 8 4 arrayGet 'pdDM' NB. ... shape required ... arrayClose 'pdDM'
load 'plot' pd 'reset' pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDK pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDL pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDM pd 'aspect 1' pd 'pdf'
... three geodesics (blue red and green) in a single rotating energy distribution (y2=0) ...
load 'plot' pd 'reset' pd (cy3"1;(((cy0^2:)+(cy1^2:)+cy2^2:)^0.5"_)"1)@:((S00yL{])"2) ((0+i.@(20000"_)){]) pdDK NB. ... colour blue ... pd 'show'
... a pattern (based on the geodesic) ...
load 'plot' pd 'reset' pd (cy3"1;(((cy0^2:)+(cy1^2:)+cy2^2:)^0.5"_)"1)@:((S00yL{])"2) ((0+i.@(00001"_)){]) pdDL NB. ... change colour ... pd (cy3"1;(((cy0^2:)+(cy1^2:)+cy2^2:)^0.5"_)"1)@:((S00yL{])"2) ((0+i.@(05000"_)){]) pdDL NB. ... colour red ... pd 'show'
... a pattern (based on the geodesic) ...
load 'plot' pd 'reset' pd (cy3"1;(((cy0^2:)+(cy1^2:)+cy2^2:)^0.5"_)"1)@:((S00yL{])"2) ((0+i.@(00001"_)){]) pdDM NB. ... change colour ... pd (cy3"1;(((cy0^2:)+(cy1^2:)+cy2^2:)^0.5"_)"1)@:((S00yL{])"2) ((0+i.@(00001"_)){]) pdDM NB. ... change colour ... pd (cy3"1;(((cy0^2:)+(cy1^2:)+cy2^2:)^0.5"_)"1)@:((S00yL{])"2) ((0+i.@(14000"_)){]) pdDM NB. ... colour green ... pd 'show'
... a pattern (based on the geodesic) ...
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 (((S06ds{])*i.@(S0700nmbr{]))@{:;veloc@:((S00yL{])"2)) (i.@(6280*32"_){]) pdDK NB. ... velocity along geodesic (blue) ... pd (((S06ds{])*i.@(S0700nmbr{]))@{:;veloc@:((S00yL{])"2)) (i.@(6000*32"_){]) pdDL NB. ... velocity along geodesic (red) ... pd (((S06ds{])*i.@(S0700nmbr{]))@{:;veloc@:((S00yL{])"2)) (i.@(6000*04"_){]) pdDM NB. ... velocity along geodesic (green) ... pd (((S06ds{])*i.@(S0700nmbr{]))@{:;incy2@:((S00yL{])"2)) (i.@(6280*32"_){]) pdDK NB. ... y2=0 ... pd (((S06ds{])*i.@(S0700nmbr{]))@{:;incy2@:((S00yL{])"2)) (i.@(6000*32"_){]) pdDL NB. ... y2=0 ... pd (((S06ds{])*i.@(S0700nmbr{]))@{:;incy2@:((S00yL{])"2)) (i.@(6000*04"_){]) pdDM NB. ... y2=0 ... pd 'pdf'
... velocity along geodesics as a fraction of the velocity of electromagnetic waves (temperature) ...