User:Tom Allen/TensorExperimentsR21

From J Wiki
Jump to navigation Jump to search


GTER2100E.jpg


load 'plot'

gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2100G.jpg
... divergence of x0 (rfE1cn) ...


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

   gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

   xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 0,_ 0 2p1 26,:_ 0 2p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(1{])"1@rfE1cn xRef@hkxR03)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


load 'plot'

gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_

pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_

pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_

pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_

pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_

pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_

pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_

pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2101D.jpg
... divergence of x2 (rfE1cn) ...


load 'plot'

gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2102D.jpg
... divergence of x3 (rfE1cn) ...


load 'plot'

gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2103D.jpg
... electric field (x0) ...


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

   gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

   xCpts=:_ 0.4 2.3 26,_ 0.5p1 0.5p1 0,_ 0 2p1 26,:_ 0 2p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(1{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


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

   gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

   xCpts=:_ 0.4 2.3 26,_ 0.4p1 0.6p1 0,_ 0 2p1 26,:_ 0 2p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


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

   gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

   xCpts=:_ 0.4 2.0 26,_ 0.5p1 0.5p1 0,_ 0 2p1 26,:_ 0 2p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(0{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


load 'plot'

gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_

pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_

pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_

pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_

pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_

pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_

pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_

pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2104E.jpg
... magnetic field (x1) ...


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

   gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

   xCpts=:_ 0.4 2.0 26,_ 0.5p1 0.5p1 0,_ 0 2p1 26,:_ 0 2p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(2{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


load 'plot'

gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2105G.jpg
... charge density ...


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=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

Initx =:  1.03 , 0.5p1 , 0.25p1 , 0p1   "_
Initdx=: (0    , 0     , 0.01   , 0.0100"_)*%@(1"_)

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:(0,628"_)*1"_
pdDQ      =:vSgD@vS''

pdDQ arrayPut   'pdDQ'
     arrayClose 'pdDQ'

pdDQds=:Initds2''

Initx =:  1.03 , 0.5p1 , 0.25p1 , 0p1   "_
Initdx=: (0    , 0     , 0.01   , 0.0106"_)*%@(1"_)

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:(0,628"_)*1"_
pdDF      =:vSgD@vS''

pdDF arrayPut   'pdDF'
     arrayClose 'pdDF'

pdDFds=:Initds2''


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

              arrayOpen  'pdDQ'
pdDQ=:629 8 4 arrayGet   'pdDQ'   NB. ... shape required ...
              arrayClose 'pdDQ'

              arrayOpen  'pdDF'
pdDF=:629 8 4 arrayGet   'pdDF'   NB. ... shape required ...
              arrayClose 'pdDF'


load 'plot'

pd ((1.03*cos@]"0);1.03*sin@]"0) pTsteps _ 0p1 2p1 160   NB. ... ds=0 (EM wave) ...
pd vSgDplot pdDQ                                         NB. ... >speed limit   ...
pd vSgDplot pdDF                                         NB. ... <speed limit   ...

pd 'aspect 1'
pd 'pdf'
pd 'show'


JTER2106K.jpg
... geodesics ...


load 'plot'

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);incy2)@:(rfyc"1)@:((S00xL{])"2) (i.@(629"_){]) pdDQ    NB. ... incy2 is zero ...

pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) (i.@(629"_){]) pdDQ

pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) (i.@(629"_){]) pdDF

pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) (i.@(629"_){]) pdDF    NB. ... incy2 is zero ...

pd 'pdf'
pd 'show'


JTER2107I.jpg
... velocity with elapsed time along geodesics ...


GTER2101D.jpg


...
... try to contain the geodesic in a standing wave ...
...


load 'plot'

gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.00p1 0.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.25p1 0.25p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.40p1 0.40p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.50p1 0.50p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.60p1 0.60p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.75p1 0.75p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 1.00p1 1.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2108G.jpg
... a standing wave (8 cycles x2) ...


load 'plot'

gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.00p1 0.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.25p1 0.25p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.40p1 0.40p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.50p1 0.50p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.60p1 0.60p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.75p1 0.75p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 1.00p1 1.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2109H.jpg
... divergence of x0 (rfE1cn) ...


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

   gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
   gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

   xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(1{])"1@rfE1cn xRef@hkxR03S)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


load 'plot'

gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.00p1 0.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.25p1 0.25p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.40p1 0.40p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.50p1 0.50p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.60p1 0.60p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.75p1 0.75p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 1.00p1 1.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2110D.jpg
... divergence of x2 (rfE1cn) ...


load 'plot'

gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.00p1 0.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.25p1 0.25p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.40p1 0.40p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.50p1 0.50p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.60p1 0.60p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.75p1 0.75p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 1.00p1 1.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2111E.jpg
... divergence of x3 (rfE1cn) ...


load 'plot'

gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.00p1 0.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.25p1 0.25p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.40p1 0.40p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.50p1 0.50p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.60p1 0.60p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.75p1 0.75p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 1.00p1 1.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2112C.jpg
... electric field x0 (standing wave) ...


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

   gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
   gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

   xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(1{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


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

   gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
   gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

   xCpts=:_ 0.4 1.7 26,_ 0.4p1 0.6p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


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

   gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
   gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

   xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(0{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


load 'plot'

gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.00p1 0.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.25p1 0.25p1 0*0.125"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.40p1 0.40p1 0*0.125"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.50p1 0.50p1 0*0.125"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.60p1 0.60p1 0*0.125"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.75p1 0.75p1 0*0.125"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 1.00p1 1.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2113C.jpg
... magnetic field x1 (standing wave) ...


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

   gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
   gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

   xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_

   (0=]) ((2^_44)qtsz]) (gXR([(2{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


load 'plot'

gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.00p1 0.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.25p1 0.25p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.40p1 0.40p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.50p1 0.50p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.60p1 0.60p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.75p1 0.75p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 1.00p1 1.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2114C.jpg
... energy flow EM x2 (standing wave) ...


load 'plot'

gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.00p1 0.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.25p1 0.25p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.40p1 0.40p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.50p1 0.50p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.60p1 0.60p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.75p1 0.75p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts''

xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 1.00p1 1.00p1 0*0.125"_

pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER2115C.jpg
... charge density (standing wave) ...


vSgDd2M=.[([(0{])@rfE1cn xRef@hkxR03S),:@(S03xH{])
vSgDd2 =.gXR(-@(([([(0{])@ch2k xRef@hkxR03S),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)]

Initds2=.(gXR([(0{])@mcv xRef@hkxR03S),:@Initx)+/^:2@:*(]*/])@Initdx
Initds =.%:@|@Initds2

gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_)
gXSAll=:,:@(0.15 1 2 8 1 1 8  8 0 0"_)

Initx =:  0.94 , 0.5p1 , 0.03125p1 , 0p1 "_
Initdx=: (0    , 0     , 0.001     , 0.04"_)*%@(1"_)

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:(0,1000"_)*1"_
pdDT      =:vSgD@vS''

pdDT arrayPut   'pdDT'
     arrayClose 'pdDT'

Initx =:  0.94 , 0.5p1 , 0.03125p1 , 0p1 "_
Initdx=: (0    , 0     , 0.001     , 0.04"_)*%@(2"_)

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:(0,1000"_)*2"_
pdDU      =:vSgD@vS''

pdDU arrayPut   'pdDU'
     arrayClose 'pdDU'

Initx =:  0.94 , 0.5p1 , 0.03125p1 , 0p1 "_
Initdx=: (0    , 0     , 0.001     , 0.04"_)*%@(4"_)

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:(0,1000"_)*4"_
pdDV      =:vSgD@vS''

pdDV arrayPut   'pdDV'
     arrayClose 'pdDV'

Initx =:  0.94 , 0.5p1 , 0.03125p1 , 0p1 "_
Initdx=: (0    , 0     , 0.001     , 0.04"_)*%@(80"_)

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:(0,1000"_)*80"_
pdDW      =:vSgD@vS''

pdDW arrayPut   'pdDW'
     arrayClose 'pdDW'


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

                arrayOpen  'pdDT'
pdDT=:01001 8 4 arrayGet   'pdDT'   NB. ... shape required ...
                arrayClose 'pdDT'

                arrayOpen  'pdDU'
pdDU=:02001 8 4 arrayGet   'pdDU'   NB. ... shape required ...
                arrayClose 'pdDU'

                arrayOpen  'pdDV'
pdDV=:04001 8 4 arrayGet   'pdDV'   NB. ... shape required ...
                arrayClose 'pdDV'

                arrayOpen  'pdDW'
pdDW=:80001 8 4 arrayGet   'pdDW'   NB. ... shape required ...
                arrayClose 'pdDW'


load 'plot'

pd ((0.94*cos@]"0);0.94*sin@]"0) pTsteps _ 0p1 0.0625p1 160
pd vSgDplot pdDT
pd vSgDplot pdDU
pd vSgDplot pdDV
pd vSgDplot pdDW

pd 'aspect 1'
pd 'pdf'
pd 'show'


JTER2116C.jpg
... a geodesic ...


load 'plot'

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)@:(rfyc"1)@:((S00xL{])"2) pdDT
pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) pdDT    NB. ... incy2 is zero ...

pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) pdDU
pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) pdDU    NB. ... incy2 is zero ...

pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) pdDV
pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) pdDV    NB. ... incy2 is zero ...

pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) pdDW
pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) pdDW    NB. ... incy2 is zero ...

pd 'pdf'
pd 'show'


JTER2117C.jpg
... as interval of integration -> 0 maybe confined ...



Download: File:LTER2100E.txt
Download: File:LTER2101D.txt