User:Tom Allen/TensorExperimentsB13

From J Wiki
Jump to navigation Jump to search


...

... ADDENDUM ...

...


GTEB1300C.jpg


NB. ... script experimentalZ.ijs ...

NB. ... these verbs are valid only in Cartesian coordinates ...

FEMcvcv    =:+/"1@((<0 1)|:])"4@B2k
FEMcvcvdn  =:+/"1@((<0 1)|:])"5@B2kdc
FEMcvcvdndn=:+/"1@((<0 1)|:])"6@B2kdcdc



NB. ... script experimentalZ.ijs ...

NB. ... these verbs are valid only in Cartesian coordinates ...

NB. ... electric and magnetic fields ...

FEMcvcv01=:((<0 1){])"2@FEMcvcv
FEMcvcv02=:((<0 2){])"2@FEMcvcv
FEMcvcv03=:((<0 3){])"2@FEMcvcv
FEMcvcv10=:((<1 0){])"2@FEMcvcv
FEMcvcv12=:((<1 2){])"2@FEMcvcv
FEMcvcv13=:((<1 3){])"2@FEMcvcv
FEMcvcv20=:((<2 0){])"2@FEMcvcv
FEMcvcv21=:((<2 1){])"2@FEMcvcv
FEMcvcv23=:((<2 3){])"2@FEMcvcv
FEMcvcv30=:((<3 0){])"2@FEMcvcv
FEMcvcv31=:((<3 1){])"2@FEMcvcv
FEMcvcv32=:((<3 2){])"2@FEMcvcv

NB. ... used for Maxwell's equation (2) ...

FEMcvcvdn231=:((<2 3 1){])"3@FEMcvcvdn
FEMcvcvdn312=:((<3 1 2){])"3@FEMcvcvdn
FEMcvcvdn123=:((<1 2 3){])"3@FEMcvcvdn

FEMcvcvdn302=:((<3 0 2){])"3@FEMcvcvdn
FEMcvcvdn023=:((<0 2 3){])"3@FEMcvcvdn
FEMcvcvdn230=:((<2 3 0){])"3@FEMcvcvdn

FEMcvcvdn013=:((<0 1 3){])"3@FEMcvcvdn
FEMcvcvdn130=:((<1 3 0){])"3@FEMcvcvdn
FEMcvcvdn301=:((<3 0 1){])"3@FEMcvcvdn

FEMcvcvdn120=:((<1 2 0){])"3@FEMcvcvdn
FEMcvcvdn201=:((<2 0 1){])"3@FEMcvcvdn
FEMcvcvdn012=:((<0 1 2){])"3@FEMcvcvdn

NB. ... used for Maxwell's equation (1) ...

FEMcvcvdn011=:((<0 1 1){])"3@FEMcvcvdn
FEMcvcvdn022=:((<0 2 2){])"3@FEMcvcvdn
FEMcvcvdn033=:((<0 3 3){])"3@FEMcvcvdn

FEMcvcvdn100=:((<1 0 0){])"3@FEMcvcvdn
FEMcvcvdn122=:((<1 2 2){])"3@FEMcvcvdn
FEMcvcvdn133=:((<1 3 3){])"3@FEMcvcvdn

FEMcvcvdn200=:((<2 0 0){])"3@FEMcvcvdn
FEMcvcvdn211=:((<2 1 1){])"3@FEMcvcvdn
FEMcvcvdn233=:((<2 3 3){])"3@FEMcvcvdn

FEMcvcvdn300=:((<3 0 0){])"3@FEMcvcvdn
FEMcvcvdn311=:((<3 1 1){])"3@FEMcvcvdn
FEMcvcvdn322=:((<3 2 2){])"3@FEMcvcvdn



NB. ... script experimentalZ.ijs ...

NB. ... these verbs are valid only in Cartesian coordinates ...

EEMcn0=:FEMcvcv30
EEMcn1=:FEMcvcv31
EEMcn2=:FEMcvcv32

EEMcn=:(0|:])@(EEMcn0,EEMcn1,:EEMcn2)

HEMcn0=:FEMcvcv12
HEMcn1=:FEMcvcv20
HEMcn2=:FEMcvcv01

HEMcn=:(0|:])@(HEMcn0,HEMcn1,:HEMcn2)



NB. ... script experimentalZ.ijs ...

NB. ... these verbs are valid only in Cartesian coordinates ...

PEMcn0=:(EEMcn1*HEMcn2)+-@(EEMcn2*HEMcn1)
PEMcn1=:(EEMcn2*HEMcn0)+-@(EEMcn0*HEMcn2)
PEMcn2=:(EEMcn0*HEMcn1)+-@(EEMcn1*HEMcn0)

PEMcn=:(0|:])@(PEMcn0,PEMcn1,:PEMcn2)



NB. ... script experimentalZ.ijs ...

NB. ... these verbs are valid only in Cartesian coordinates ...

JEMcn  =:+/"1@((<1 2)|:])"3@FEMcvcvdn
JEMcndn=:+/"1@((<1 2)|:])"4@FEMcvcvdndn


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

   gXcw   =:  1  0  _    _   _"_
   gXpAll =:  1  1  2.2  20  _1"_
   gXq2All=:  1  1  2.2  20   _"_
   gXq3All=:  0  0  1    1    0"_
   gXq4All=:  1  1  2.2  1    _"_
   gXq5All=:  2  0  1    1   _2"_
   gXq6All=:  1  1  2.2  0    _"_
   gXq7All=:  1  1  3    1   _2"_

   xCpts=:_ 2.78p1 5.12p1 50,_ 0.000001 1p1 50,_ 0.000001 2p1 50,:_ 0 4p1 50"_

NB. ... Maxwell's equation (2) ...

   *./^:3"3 (0=]) ((2^_44)qtsz]) (]+(0|:])+0 1|:])"3 (gXB([+/"1@((<0 1)|:])"5@B2kdc yRef@hkyI)]) ((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. ... derivatives ...

   A=:(gXB([FEMcvcvdn yRef@hkyI)])                               ((0;50"_)pTRandom xCpts)''
   B=:(gXB(0|:[:(gXB([(0{])@FEMcvcv yRef@hkyI),:@])D.1])"_ 1])   ((0;50"_)pTRandom xCpts)''

   A (((2^_20);2^_07)qteq[;])"3 B
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

   A=:(gXB([FEMcvcvdndn yRef@hkyI)])                             ((0;50"_)pTRandom xCpts)''
   B=:(gXB(0|:[:(gXB([(0{])@FEMcvcvdn yRef@hkyI),:@])D.1])"_ 1]) ((0;50"_)pTRandom xCpts)''

   A (((2^_17);2^_07)qteq[;])"4 B
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 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. ... Maxwell's equation (2) ...

   *./^:3"3 (0=]) ((2^_44)qtsz]) (]+(0|:])+0 1|:])"3 (gXB([FEMcvcvdn yRef@hkyI)]) ((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. ... electric and magnetic fields ...

   $ (gXB([EEMcn xRef@hkyI)]) ((0;50"_)pTRandom xCpts)''
50 3
   $ (gXB([HEMcn xRef@hkyI)]) ((0;50"_)pTRandom xCpts)''
50 3

NB. ... Maxwell's equation (2) ...

   (0=]) ((2^_44)qtsz]) (gXB([(FEMcvcvdn231+FEMcvcvdn312+FEMcvcvdn123)yRef@hkyI)]) ((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

   (0=]) ((2^_44)qtsz]) (gXB([(FEMcvcvdn302+FEMcvcvdn023+FEMcvcvdn230)yRef@hkyI)]) ((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

   (0=]) ((2^_44)qtsz]) (gXB([(FEMcvcvdn013+FEMcvcvdn130+FEMcvcvdn301)yRef@hkyI)]) ((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

   (0=]) ((2^_44)qtsz]) (gXB([(FEMcvcvdn120+FEMcvcvdn201+FEMcvcvdn012)yRef@hkyI)]) ((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. ... Poynting vector ...

   $ (gXB([PEMcn yRef@hkyI)]) ((0;50"_)pTRandom xCpts)''
50 3

NB. ... Maxwell's equation (1) ...

NB. ... derivatives ...

   A=:(gXB([JEMcndn yRef@hkyI)])                             ((0;50"_)pTRandom xCpts)''
   B=:(gXB(0|:[:(gXB([(0{])@JEMcn yRef@hkyI),:@])D.1])"_ 1]) ((0;50"_)pTRandom xCpts)''

   A (((2^_17);2^_07)qteq[;])"2 B
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 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. ... conservation of charge ...

NB. ... divergence is 0 ...

   (0=]) ((2^_44)qtsz]) (gXB([+/@((<0 1)|:])"2@JEMcndn yRef@hkyI)]) ((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


GTEB1302C.jpg


load 'plot'

gXcw   =:  1  0  _    _   _"_
gXpAll =:  1  1  2.2  20  _1"_
gXq2All=:  1  1  2.2  20   _"_
gXq3All=:  0  0  1    1    0"_
gXq4All=:  1  1  2.2  1    _"_
gXq5All=:  2  0  1    1   _2"_
gXq6All=:  1  1  2.2  0    _"_
gXq7All=:  1  1  3    1   _2"_

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 1.75p1 1.75p1 0"_

pd (gXB([((0{])"1@xpts;(0{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.00p1 2.00p1 0"_

pd (gXB([((0{])"1@xpts;(0{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.15p1 2.15p1 0"_

pd (gXB([((0{])"1@xpts;(0{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.25p1 2.25p1 0"_

pd (gXB([((0{])"1@xpts;(0{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.35p1 2.35p1 0"_

pd (gXB([((0{])"1@xpts;(0{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.50p1 2.50p1 0"_

pd (gXB([((0{])"1@xpts;(0{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.75p1 2.75p1 0"_

pd (gXB([((0{])"1@xpts;(0{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

pd 'pdf'
pd 'show'


Egraphic1300C.jpg
... a standing wave (x0 direction) ...


GTEB1303C.jpg


load 'plot'

gXcw   =:  1  0  _    _   _"_
gXpAll =:  1  1  2.2  20  _1"_
gXq2All=:  1  1  2.2  20   _"_
gXq3All=:  0  0  1    1    0"_
gXq4All=:  1  1  2.2  1    _"_
gXq5All=:  2  0  1    1   _2"_
gXq6All=:  1  1  2.2  0    _"_
gXq7All=:  1  1  3    1   _2"_

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 1.75p1 1.75p1 0"_

pd (gXB([((0{])"1@xpts;(1{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.00p1 2.00p1 0"_

pd (gXB([((0{])"1@xpts;(1{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.15p1 2.15p1 0"_

pd (gXB([((0{])"1@xpts;(1{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.25p1 2.25p1 0"_

pd (gXB([((0{])"1@xpts;(1{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.35p1 2.35p1 0"_

pd (gXB([((0{])"1@xpts;(1{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.50p1 2.50p1 0"_

pd (gXB([((0{])"1@xpts;(1{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.75p1 2.75p1 0"_

pd (gXB([((0{])"1@xpts;(1{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

pd 'pdf'
pd 'show'


Egraphic1301C.jpg
... a standing wave (x1 direction) ...


GTEB1304C.jpg


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

   gXcw   =:  1  0  _    _   _"_
   gXpAll =:  1  1  2.2  20  _1"_
   gXq2All=:  1  1  2.2  20   _"_
   gXq3All=:  0  0  1    1    0"_
   gXq4All=:  1  1  2.2  1    _"_
   gXq5All=:  2  0  1    1   _2"_
   gXq6All=:  1  1  2.2  0    _"_
   gXq7All=:  1  1  3    1   _2"_

   xCpts=:_ 2.78p1 5.12p1 50,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 1.75p1 3.75p1 50"_

   (0=]) ((2^_44)qtsz]) (gXB([(2{])"1@Ocn xRef@hkyI)]) (gXB yc]) ((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


GTEB1305C.jpg


load 'plot'

gXcw   =:  1  0  _    _   _"_
gXpAll =:  1  1  2.2  20  _1"_
gXq2All=:  1  1  2.2  20   _"_
gXq3All=:  0  0  1    1    0"_
gXq4All=:  1  1  2.2  1    _"_
gXq5All=:  2  0  1    1   _2"_
gXq6All=:  1  1  2.2  0    _"_
gXq7All=:  1  1  3    1   _2"_

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 1.75p1 1.75p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.00p1 2.00p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.15p1 2.15p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.25p1 2.25p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.35p1 2.35p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.50p1 2.50p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.75p1 2.75p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 3.00p1 3.00p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 3.15p1 3.15p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 3.25p1 3.25p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 3.35p1 3.35p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 3.50p1 3.50p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

xCpts=:_ 2.78p1 5.12p1 200,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 3.75p1 3.75p1 0"_

pd (gXB([((0{])"1@xpts;(3{])"1@Ocn)xRef@hkyI)]) (gXB yc]) pTtotal xCpts''

pd 'pdf'
pd 'show'


Egraphic1302C.jpg
... energy density ...


GTEB1306C.jpg


load 'plot'

gXcw   =:  1  0  _    _   _"_
gXpAll =:  1  1  2.2  20  _1"_
gXq2All=:  1  1  2.2  20   _"_
gXq3All=:  0  0  1    1    0"_
gXq4All=:  1  1  2.2  1    _"_
gXq5All=:  2  0  1    1   _2"_
gXq6All=:  1  1  2.2  0    _"_
gXq7All=:  1  1  3    1   _2"_

xCpts=:_ 3.75p1 3.75p1 0,_ 0.000001p1 1.0p1 2000,_ 0.25p1 0.25p1 0,:_ 2.00p1 2.00p1 0"_

pd (gXB([(cx1"1@xpts;(0{])"1@(}:"1@}:"2@xdy smx"2 1 EEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... blue    electric field  x0 ...

pd (gXB([(cx1"1@xpts;(1{])"1@(}:"1@}:"2@xdy smx"2 1 EEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... red     electric field  x1 ...

pd (gXB([(cx1"1@xpts;(2{])"1@(}:"1@}:"2@xdy smx"2 1 EEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... green   electric field  x2 ...

pd (gXB([(cx1"1@xpts;(3{])"1@(                      JEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... purple  charge density     ...

pd 'pdf'
pd 'show'


Egraphic1303C.jpg
... electric field components and charge density ...


GTEB1307C.jpg


load 'plot'

gXcw   =:  1  0  _    _   _"_
gXpAll =:  1  1  2.2  20  _1"_
gXq2All=:  1  1  2.2  20   _"_
gXq3All=:  0  0  1    1    0"_
gXq4All=:  1  1  2.2  1    _"_
gXq5All=:  2  0  1    1   _2"_
gXq6All=:  1  1  2.2  0    _"_
gXq7All=:  1  1  3    1   _2"_

xCpts=:_ 3.75p1 3.75p1 0,_ 0.000001p1 1.0p1 2000,_ 0.25p1 0.25p1 0,:_ 2.00p1 2.00p1 0"_

pd (gXB([(cx1"1@xpts;(0{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... blue    magnetic field  x0 ...

pd (gXB([(cx1"1@xpts;(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... red     magnetic field  x1 ...

pd (gXB([(cx1"1@xpts;(2{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... green   magnetic field  x2 ...

pd 'pdf'
pd 'show'


Egraphic1304C.jpg
... magnetic field components ...


GTEB1308C.jpg


load 'plot'

gXcw   =:  1  0  _    _   _"_
gXpAll =:  1  1  2.2  20  _1"_
gXq2All=:  1  1  2.2  20   _"_
gXq3All=:  0  0  1    1    0"_
gXq4All=:  1  1  2.2  1    _"_
gXq5All=:  2  0  1    1   _2"_
gXq6All=:  1  1  2.2  0    _"_
gXq7All=:  1  1  3    1   _2"_

xCpts=:_ 3.75p1 3.75p1 0,_ 0.000001p1 1.0p1 2000,_ 0.25p1 0.25p1 0,:_ 2.00p1 2.00p1 0"_

pd (gXB([(cx1"1@xpts;(0{])"1@(}:"1@}:"2@xdy smx"2 1 PEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... blue    Poynting vector  x0 ...

pd (gXB([(cx1"1@xpts;(1{])"1@(}:"1@}:"2@xdy smx"2 1 PEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... red     Poynting vector  x1 ...

pd (gXB([(cx1"1@xpts;(2{])"1@(}:"1@}:"2@xdy smx"2 1 PEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... green   Poynting vector  x2 ...

pd 'pdf'
pd 'show'


Egraphic1305C.jpg
... Poynting vector components ...


GTEB1309C.jpg


load 'plot'

gXcw   =:  1  0  _    _   _"_
gXpAll =:  1  1  2.2  20  _1"_
gXq2All=:  1  1  2.2  20   _"_
gXq3All=:  0  0  1    1    0"_
gXq4All=:  1  1  2.2  1    _"_
gXq5All=:  2  0  1    1   _2"_
gXq6All=:  1  1  2.2  0    _"_
gXq7All=:  1  1  3    1   _2"_

xCpts=:_ 3.75p1 3.75p1 0,_ 0.000001p1 1.0p1 2000,_ 0.25p1 0.25p1 0,:_ 2.00p1 2.00p1 0"_

pd (gXB([(cx1"1@xpts;(3{])"1@(                      Ocn  ))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... blue    energy density     ...

pd (gXB([(cx1"1@xpts;(3{])"1@(                      JEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... red     charge density     ...

pd (gXB([(cx1"1@xpts;(0{])"1@(}:"1@}:"2@xdy smx"2 1 EEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... green   electric field  x0 ...

pd (gXB([(cx1"1@xpts;(1{])"1@(}:"1@}:"2@xdy smx"2 1 EEMcn))yRef@hkyI)]) (gXB yc]) pTtotal xCpts''    NB. ... purple  electric field  x1 ...

pd 'pdf'
pd 'show'


Egraphic1306C.jpg
...



Download: File:LTEB1300C.txt
Download: File:LTEB1302C.txt
Download: File:LTEB1303C.txt
Download: File:LTEB1304C.txt
Download: File:LTEB1305C.txt
Download: File:LTEB1306C.txt
Download: File:LTEB1307C.txt
Download: File:LTEB1308C.txt
Download: File:LTEB1309C.txt




Download MoinMoin source: File:TensorExperimentsB13.ijs