User:Tom Allen/TensorExperimentsB13
Jump to navigation
Jump to search
...
... ADDENDUM ...
...
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
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'
... a standing wave (x0 direction) ...
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'
... a standing wave (x1 direction) ...
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
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'
... energy density ...
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'
... electric field components and charge density ...
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'
... magnetic field components ...
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'
... Poynting vector components ...
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'
...
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