User:Tom Allen/ta/AddendumX04P06D
< User:Tom Allen | ta
Jump to navigation
Jump to search
FOR THIS MODEL
... conservation of total energy (+/-) ...
load 'plot' pd 'reset' uQ0=:,:@(1 1 2 8 1 1 1 _0.5 0"_) x46DPsel=:(S00yL{])"2@(i.@(2516"_){])@(,@(80000$1 16$(15$0:),1:)#}:@]) x46DEtot=:(uQ0([Ocndc yRef@hkyX04A )]) x46DPsel pdDB x46Ddy0=:((<0 0){])"2@(x46DEtot"_) x46Ddy1=:((<1 1){])"2@(x46DEtot"_) x46Ddy2=:((<2 2){])"2@(x46DEtot"_) x46Ddy3=:((<3 3){])"2@(x46DEtot"_) pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2 )]) x46DPsel pdDB NB. ... rate of change (y0 y1 y2) ... pd (uQ0(cy3"_ 1;x46Ddy3 )]) x46DPsel pdDB NB. ... rate of change (energy density) ... pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2+x46Ddy3)]) x46DPsel pdDB NB. ... conservation of energy ... pd 'pdf'
... conservation of charge (+/-) ...
load 'plot' pd 'reset' uQ0=:,:@(1 1 2 8 1 1 1 _0.5 0"_) x46DPsel=:(S00yL{])"2@(i.@(1246"_){])@(,@(80000$1 16$(15$0:),1:)#}:@]) NB. ... at beginning ... x46DEtot0=:(uQ0([Ocndc yRef@hkyX04A )]) x46DPsel pdDB x46Ddy0=:((<0 0){])"2@(x46DEtot0"_) x46Ddy1=:((<1 1){])"2@(x46DEtot0"_) x46Ddy2=:((<2 2){])"2@(x46DEtot0"_) x46Ddy3=:((<3 3){])"2@(x46DEtot0"_) pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2 )]) x46DPsel pdDB NB. ... rate of change (y0 y1 y2) ... pd (uQ0(cy3"_ 1;x46Ddy3 )]) x46DPsel pdDB NB. ... rate of change (energy density) ... pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2+x46Ddy3)]) x46DPsel pdDB NB. ... conservation of energy ... x46DEtot1=:(uQ0([JEMcndn yRef@hkyX04A )]) x46DPsel pdDB x46Ddy0=:((<0 0){])"2@(x46DEtot1"_) x46Ddy1=:((<1 1){])"2@(x46DEtot1"_) x46Ddy2=:((<2 2){])"2@(x46DEtot1"_) x46Ddy3=:((<3 3){])"2@(x46DEtot1"_) pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2 )]) x46DPsel pdDB NB. ... rate of change (y0 y1 y2) ... pd (uQ0(cy3"_ 1;x46Ddy3 )]) x46DPsel pdDB NB. ... rate of change (charge density) ... pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2+x46Ddy3)]) x46DPsel pdDB NB. ... conservation of charge ... pd 'pdf'
load 'plot' pd 'reset' uQ0=:,:@(1 1 2 8 1 1 1 _0.5 0"_) x46DPsel=:(S00yL{])"2@((78112+i.@(1246"_)){])@(,@(80000$1 16$(15$0:),1:)#}:@]) NB. ... after beginning ... x46DEtot2=:(uQ0([Ocndc yRef@hkyX04A )]) x46DPsel pdDB x46Ddy0=:((<0 0){])"2@(x46DEtot2"_) x46Ddy1=:((<1 1){])"2@(x46DEtot2"_) x46Ddy2=:((<2 2){])"2@(x46DEtot2"_) x46Ddy3=:((<3 3){])"2@(x46DEtot2"_) pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2 )]) x46DPsel pdDB NB. ... rate of change (y0 y1 y2) ... pd (uQ0(cy3"_ 1;x46Ddy3 )]) x46DPsel pdDB NB. ... rate of change (energy density) ... pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2+x46Ddy3)]) x46DPsel pdDB NB. ... conservation of energy ... x46DEtot3=:(uQ0([JEMcndn yRef@hkyX04A )]) x46DPsel pdDB x46Ddy0=:((<0 0){])"2@(x46DEtot3"_) x46Ddy1=:((<1 1){])"2@(x46DEtot3"_) x46Ddy2=:((<2 2){])"2@(x46DEtot3"_) x46Ddy3=:((<3 3){])"2@(x46DEtot3"_) pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2 )]) x46DPsel pdDB NB. ... rate of change (y0 y1 y2) ... pd (uQ0(cy3"_ 1;x46Ddy3 )]) x46DPsel pdDB NB. ... rate of change (charge density) ... pd (uQ0(cy3"_ 1;x46Ddy0+x46Ddy1+x46Ddy2+x46Ddy3)]) x46DPsel pdDB NB. ... conservation of charge ... pd 'pdf'
Download: File:LX04P06D00D.txt
Download: File:LX04P06D01D.txt
Download: File:LX04P06D02C.txt