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