User:Tom Allen/TensorExperimentsR20
Jump to navigation
Jump to search
... charge negative (_1) ...
vSgDd2M=.[([(0{])@rfE1cn xRef@hkxR03),:@(S03xH{]) vSgDd2 =.gXR(-@(([([(0{])@ch2k xRef@hkxR03),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)] Initds2=.(gXR([(0{])@mcv xRef@hkxR03),:@Initx)+/^:2@:*(]*/])@Initdx Initds =.%:@|@Initds2 gXRAll=:,:@(0.5 1 2 8 1 1 1 _0.2 0 _1"_) Initx =: 1.1000 , 0.5p1 , 0.25p1 , 0p1 "_ Initdx=: (0 , 0 , 0.001 , 0.01"_)*%@(8"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,3078"_)*8"_ pdED =:vSgD@vS'' pdED arrayPut 'pdED' arrayClose 'pdED'
NB. ... execute (ijx) ... arrayOpen 'pdED' pdED=:24625 8 4 arrayGet 'pdED' NB. ... shape required ... arrayClose 'pdED'
load 'plot' pd ((1.1000*cos@]"0);1.1000*sin@]"0) pTsteps _ 0p1 2p1 160 pd vSgDplot pdED pd 'aspect 1' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... (2p1+0.25p1),2p1*5 7.06858 31.4159 {.{:pdED 1.03452 1.5708 7.06367 31.4589
... charge neutral (0) ...
vSgDd2M=.[([(0{])@rfE1cn xRef@hkxR03),:@(S03xH{]) vSgDd2 =.gXR(-@(([([(0{])@ch2k xRef@hkxR03),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)] Initds2=.(gXR([(0{])@mcv xRef@hkxR03),:@Initx)+/^:2@:*(]*/])@Initdx Initds =.%:@|@Initds2 gXRAll=:,:@(0.5 1 2 8 1 1 1 _0.2 0 0"_) Initx =: 1.1000 , 0.5p1 , 0.25p1 , 0p1 "_ Initdx=: (0 , 0 , 0.001 , 0.01"_)*%@(8"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,3067"_)*8"_ pdFD =:vSgD@vS'' pdFD arrayPut 'pdFD' arrayClose 'pdFD'
NB. ... execute (ijx) ... arrayOpen 'pdFD' pdFD=:24537 8 4 arrayGet 'pdFD' NB. ... shape required ... arrayClose 'pdFD'
load 'plot' pd vSgDplot pdFD pd ((1.1000*cos@]"0);1.1000*sin@]"0) pTsteps _ 0p1 2p1 160 pd 'aspect 1' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... (2p1+0.25p1),2p1*5 7.06858 31.4159 {.{:pdFD 1.07986 1.5708 7.0639 31.424
... charge positive (1) ...
vSgDd2M=.[([(0{])@rfE1cn xRef@hkxR03),:@(S03xH{]) vSgDd2 =.gXR(-@(([([(0{])@ch2k xRef@hkxR03),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)] Initds2=.(gXR([(0{])@mcv xRef@hkxR03),:@Initx)+/^:2@:*(]*/])@Initdx Initds =.%:@|@Initds2 gXRAll=:,:@(0.5 1 2 8 1 1 1 _0.2 0 1"_) Initx =: 1.1000 , 0.5p1 , 0.25p1 , 0p1 "_ Initdx=: (0 , 0 , 0.001 , 0.01"_)*%@(8"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,3066"_)*8"_ pdDW =:vSgD@vS'' pdDW arrayPut 'pdDW' arrayClose 'pdDW'
NB. ... execute (ijx) ... arrayOpen 'pdDW' pdDW=:24529 8 4 arrayGet 'pdDW' NB. ... shape required ... arrayClose 'pdDW'
load 'plot' pd ((1.1000*cos@]"0);1.1000*sin@]"0) pTsteps _ 0p1 2p1 160 pd ((1.1000*cos@]"0);1.1000*sin@]"0) pTsteps _ 0p1 2p1 160 pd vSgDplot pdDW pd 'aspect 1' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... (2p1+0.25p1),2p1*5 7.06858 31.4159 {.{:pdDW 1.08333 1.5708 7.06849 31.4057