User:Tom Allen/TensorExperimentsR20

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

... geodesic ...

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'

... geodesic ...

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'

... geodesic ...

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

   (2p1+0.25p1),2p1*5
7.06858 31.4159
   {.{:pdDW
1.08333 1.5708 7.06849 31.4057