User:Tom Allen/TensorExperimentsB10
Jump to navigation
Jump to search
...
...
...
... 6.25/1.75/8/4 ...
load 'plot' vSgDd2M=.[([(0{])@(([mcn xRef)smx"2 1[(ydx((0|:[)smx])"2 1 B2kCEr1+SBkCEr1)yRef)hkxI),:@(S03xH{]) vSgDd2 =.gXB(-@(([([(0{])@ch2k xRef@hkxI),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)] Initds2=.(gXB([(0{])@mcv xRef@hkxI),:@Initx)+/^:2@:*(]*/])@Initdx gXcw =: 1.5 0 _ _ _"_ gXpAll =: 1 1 1 9 _1"_ gXq2All=: 1 1 1 9 _"_ gXq3All=: 0 0 1 1 0"_ gXq4All=: 1 1 1 1 _"_ gXq5All=: 2 0 1 1 _2"_ gXq6All=: 1 1 1 0 _"_ gXq7All=: 1 1 3 1 _2"_ Initx =: 6.25p1 , 0.5p1 , 1p1 , 1.75p1"_ Initdx=: (0 , 0 , _0.01 , 0 "_)*%@(4"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,8000"_)*4"_ pd vSgDplot pdDA=:vSgD@vS'' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... NB. ... using pdDA from the plot above ... $pdDA 32001 8 4 NB. ... x0 ... A=:cx0"1@:((S00xL{])"2) pdDA (6.25p1-<./) A 0.0274124 (>./-6.25p1"_) A 0.142113 NB. ... x1 ... *./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDA 1 NB. ... x2 decreases continuously ... *./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDA 1 NB. ... x3 ... *./ (0=]) ((2^_47)qtsz]) (1.75p1-]) cx3"1@:((S00xL{])"2) pdDA 1
load 'plot' NB. ... using pdDA from the plot above ... pd (-@cx2"1;cx0"1)@:((S00xL{])"2) (i.@(20000"_){]) pdDA pd 'pdf' pd 'show'
...
...
...
... 6.375/1.75/8/4 ...
load 'plot' vSgDd2M=.[([(0{])@(([mcn xRef)smx"2 1[(ydx((0|:[)smx])"2 1 B2kCEr1+SBkCEr1)yRef)hkxI),:@(S03xH{]) vSgDd2 =.gXB(-@(([([(0{])@ch2k xRef@hkxI),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)] Initds2=.(gXB([(0{])@mcv xRef@hkxI),:@Initx)+/^:2@:*(]*/])@Initdx gXcw =: 1.5 0 _ _ _"_ gXpAll =: 1 1 1 9 _1"_ gXq2All=: 1 1 1 9 _"_ gXq3All=: 0 0 1 1 0"_ gXq4All=: 1 1 1 1 _"_ gXq5All=: 2 0 1 1 _2"_ gXq6All=: 1 1 1 0 _"_ gXq7All=: 1 1 3 1 _2"_ Initx =: 6.375p1 , 0.5p1 , 1p1 , 1.75p1"_ Initdx=: (0 , 0 , _0.01 , 0 "_)*%@(4"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,8000"_)*4"_ pd vSgDplot pdDA=:vSgD@vS'' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... NB. ... using pdDA from the plot above ... $pdDA 32001 8 4 NB. ... x0 ... A=:cx0"1@:((S00xL{])"2) pdDA (6.25p1-<./) A 0.414795 (>./-6.25p1"_) A 0.583729 NB. ... x1 ... *./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDA 1 NB. ... x2 decreases continuously ... *./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDA 1 NB. ... x3 ... *./ (0=]) ((2^_48)qtsz]) (1.75p1-]) cx3"1@:((S00xL{])"2) pdDA 1
load 'plot' NB. ... using pdDA from the plot above ... pd (-@cx2"1;cx0"1)@:((S00xL{])"2) (i.@(27600"_){]) pdDA pd 'pdf' pd 'show'
...
...
...
... 6.25/2.00/8/4 ...
load 'plot' vSgDd2M=.[([(0{])@(([mcn xRef)smx"2 1[(ydx((0|:[)smx])"2 1 B2kCEr1+SBkCEr1)yRef)hkxI),:@(S03xH{]) vSgDd2 =.gXB(-@(([([(0{])@ch2k xRef@hkxI),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)] Initds2=.(gXB([(0{])@mcv xRef@hkxI),:@Initx)+/^:2@:*(]*/])@Initdx gXcw =: 1.5 0 _ _ _"_ gXpAll =: 1 1 1 9 _1"_ gXq2All=: 1 1 1 9 _"_ gXq3All=: 0 0 1 1 0"_ gXq4All=: 1 1 1 1 _"_ gXq5All=: 2 0 1 1 _2"_ gXq6All=: 1 1 1 0 _"_ gXq7All=: 1 1 3 1 _2"_ Initx =: 6.25p1 , 0.5p1 , 1p1 , 2.00p1"_ Initdx=: (0 , 0 , _0.01 , 0 "_)*%@(4"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,8000"_)*4"_ pd vSgDplot pdDA=:vSgD@vS'' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... NB. ... using pdDA from the plot above ... $pdDA 32001 8 4 NB. ... x0 ... A=:cx0"1@:((S00xL{])"2) pdDA (6.25p1-<./) A 0.0108606 (>./-6.25p1"_) A 0.0789094 NB. ... x1 ... *./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDA 1 NB. ... x2 decreases continuously ... *./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDA 1
load 'plot' NB. ... using pdDA from the plot above ... pd (-@cx2"1;cx0"1)@:((S00xL{])"2) (i.@(18000"_){]) pdDA pd 'pdf' pd 'show'
load 'plot' NB. ... using pdDA from the plot above ... pd (-@cx2"1;cx3"1)@:((S00xL{])"2) (i.@(32000"_){]) pdDA pd 'pdf' pd 'show'
Download MoinMoin source: File:TensorExperimentsB10.ijs