User:Tom Allen/TensorExperimentsR21
Jump to navigation
Jump to search
load 'plot' gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 0,_ 0 2p1 26,:_ 0 2p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(1{])"1@rfE1cn xRef@hkxR03)]) ((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' gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_ pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_ pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_ pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_ pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_ pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_ pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_ pd (gXR(cx0"_ 1;[(2{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@rfE1cn xRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.3 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 2.3 26,_ 0.5p1 0.5p1 0,_ 0 2p1 26,:_ 0 2p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(1{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) ((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. ... execute (ijx) ... gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 2.3 26,_ 0.4p1 0.6p1 0,_ 0 2p1 26,:_ 0 2p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) ((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. ... execute (ijx) ... gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 2.0 26,_ 0.5p1 0.5p1 0,_ 0 2p1 26,:_ 0 2p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(0{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) ((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' gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_ pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_ pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_ pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_ pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_ pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_ pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_ pd (gXR(cx0"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 2.0 26,_ 0.5p1 0.5p1 0,_ 0 2p1 26,:_ 0 2p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(2{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) ((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' gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.00p1 0.00p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.25p1 0.25p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.40p1 0.40p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.50p1 0.50p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.60p1 0.60p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0.75p1 0.75p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 1.8 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 1.00p1 1.00p1 0"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
vSgDd2M=.[([(0{])@rfE1cn xRef@hkxR03),:@(S03xH{]) vSgDd2 =.gXR(-@(([([(0{])@ch2k xRef@hkxR03),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)] Initds2=.(gXR([(0{])@mcv xRef@hkxR03),:@Initx)+/^:2@:*(]*/])@Initdx Initds =.%:@|@Initds2 gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) Initx =: 1.03 , 0.5p1 , 0.25p1 , 0p1 "_ Initdx=: (0 , 0 , 0.01 , 0.0100"_)*%@(1"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,628"_)*1"_ pdDQ =:vSgD@vS'' pdDQ arrayPut 'pdDQ' arrayClose 'pdDQ' pdDQds=:Initds2'' Initx =: 1.03 , 0.5p1 , 0.25p1 , 0p1 "_ Initdx=: (0 , 0 , 0.01 , 0.0106"_)*%@(1"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,628"_)*1"_ pdDF =:vSgD@vS'' pdDF arrayPut 'pdDF' arrayClose 'pdDF' pdDFds=:Initds2''
NB. ... execute (ijx) ... arrayOpen 'pdDQ' pdDQ=:629 8 4 arrayGet 'pdDQ' NB. ... shape required ... arrayClose 'pdDQ' arrayOpen 'pdDF' pdDF=:629 8 4 arrayGet 'pdDF' NB. ... shape required ... arrayClose 'pdDF'
load 'plot' pd ((1.03*cos@]"0);1.03*sin@]"0) pTsteps _ 0p1 2p1 160 NB. ... ds=0 (EM wave) ... pd vSgDplot pdDQ NB. ... >speed limit ... pd vSgDplot pdDF NB. ... <speed limit ... pd 'aspect 1' pd 'pdf' pd 'show'
load 'plot' incy0=.(2-~/\])@:(cy0"1) incy1=.(2-~/\])@:(cy1"1) incy2=.(2-~/\])@:(cy2"1) incy3=.(2-~/\])@:(cy3"1) veloc=.(((incy0^2:)+incy1^2:)^0.5"_)%incy3 pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) (i.@(629"_){]) pdDQ NB. ... incy2 is zero ... pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) (i.@(629"_){]) pdDQ pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) (i.@(629"_){]) pdDF pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) (i.@(629"_){]) pdDF NB. ... incy2 is zero ... pd 'pdf' pd 'show'
...
... try to contain the geodesic in a standing wave ...
...
load 'plot' gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.00p1 0.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.25p1 0.25p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.40p1 0.40p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.50p1 0.50p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.60p1 0.60p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 0.75p1 0.75p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 2.03125p1 400,:_ 1.00p1 1.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@Ocn xRef@hkxR03S)]) pTtotal xCpts'' pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.00p1 0.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.25p1 0.25p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.40p1 0.40p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.50p1 0.50p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.60p1 0.60p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.75p1 0.75p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 1.00p1 1.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(1{])"1@rfE1cn xRef@hkxR03S)]) ((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' gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.00p1 0.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.25p1 0.25p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.40p1 0.40p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.50p1 0.50p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.60p1 0.60p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.75p1 0.75p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 1.00p1 1.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.00p1 0.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.25p1 0.25p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.40p1 0.40p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.50p1 0.50p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.60p1 0.60p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 0.75p1 0.75p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 400,:_ 1.00p1 1.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfE1cn xRef@hkxR03S)]) pTtotal xCpts'' pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.00p1 0.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.25p1 0.25p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.40p1 0.40p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.50p1 0.50p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.60p1 0.60p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.75p1 0.75p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 1.00p1 1.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) pTtotal xCpts'' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(1{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) ((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. ... execute (ijx) ... gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.4 1.7 26,_ 0.4p1 0.6p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03S)]) ((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. ... execute (ijx) ... gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(0{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) ((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' gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.00p1 0.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.25p1 0.25p1 0*0.125"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.40p1 0.40p1 0*0.125"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.50p1 0.50p1 0*0.125"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.60p1 0.60p1 0*0.125"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.75p1 0.75p1 0*0.125"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 1.00p1 1.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) pTtotal xCpts'' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.4 1.7 26,_ 0.5p1 0.5p1 26,_ 0.03125p1 0.28125p1 26,:_ 0 0.25p1 26"_ (0=]) ((2^_44)qtsz]) (gXR([(2{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03S)]) ((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' gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.00p1 0.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.25p1 0.25p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.40p1 0.40p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.50p1 0.50p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.60p1 0.60p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.75p1 0.75p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 1.00p1 1.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03S)]) pTtotal xCpts'' pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.00p1 0.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.25p1 0.25p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.40p1 0.40p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.50p1 0.50p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.60p1 0.60p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 0.75p1 0.75p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts'' xCpts=:_ 0.94 0.94 0,_ 0.5p1 0.5p1 0,_ 0.03125p1 0.28125p1 200,:_ 1.00p1 1.00p1 0*0.125"_ pd (gXR(cx2"_ 1;[(3{])"1@rfJEMcn yRef@hkxR03S)]) pTtotal xCpts'' pd 'pdf' pd 'show'
vSgDd2M=.[([(0{])@rfE1cn xRef@hkxR03S),:@(S03xH{]) vSgDd2 =.gXR(-@(([([(0{])@ch2k xRef@hkxR03S),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)] Initds2=.(gXR([(0{])@mcv xRef@hkxR03S),:@Initx)+/^:2@:*(]*/])@Initdx Initds =.%:@|@Initds2 gXRAll=:,:@(0.15 1 2 8 1 1 8 _8 0 0"_) gXSAll=:,:@(0.15 1 2 8 1 1 8 8 0 0"_) Initx =: 0.94 , 0.5p1 , 0.03125p1 , 0p1 "_ Initdx=: (0 , 0 , 0.001 , 0.04"_)*%@(1"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,1000"_)*1"_ pdDT =:vSgD@vS'' pdDT arrayPut 'pdDT' arrayClose 'pdDT' Initx =: 0.94 , 0.5p1 , 0.03125p1 , 0p1 "_ Initdx=: (0 , 0 , 0.001 , 0.04"_)*%@(2"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,1000"_)*2"_ pdDU =:vSgD@vS'' pdDU arrayPut 'pdDU' arrayClose 'pdDU' Initx =: 0.94 , 0.5p1 , 0.03125p1 , 0p1 "_ Initdx=: (0 , 0 , 0.001 , 0.04"_)*%@(4"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,1000"_)*4"_ pdDV =:vSgD@vS'' pdDV arrayPut 'pdDV' arrayClose 'pdDV' Initx =: 0.94 , 0.5p1 , 0.03125p1 , 0p1 "_ Initdx=: (0 , 0 , 0.001 , 0.04"_)*%@(80"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count =:(0,1000"_)*80"_ pdDW =:vSgD@vS'' pdDW arrayPut 'pdDW' arrayClose 'pdDW'
NB. ... execute (ijx) ... arrayOpen 'pdDT' pdDT=:01001 8 4 arrayGet 'pdDT' NB. ... shape required ... arrayClose 'pdDT' arrayOpen 'pdDU' pdDU=:02001 8 4 arrayGet 'pdDU' NB. ... shape required ... arrayClose 'pdDU' arrayOpen 'pdDV' pdDV=:04001 8 4 arrayGet 'pdDV' NB. ... shape required ... arrayClose 'pdDV' arrayOpen 'pdDW' pdDW=:80001 8 4 arrayGet 'pdDW' NB. ... shape required ... arrayClose 'pdDW'
load 'plot' pd ((0.94*cos@]"0);0.94*sin@]"0) pTsteps _ 0p1 0.0625p1 160 pd vSgDplot pdDT pd vSgDplot pdDU pd vSgDplot pdDV pd vSgDplot pdDW pd 'aspect 1' pd 'pdf' pd 'show'
load 'plot' incy0=.(2-~/\])@:(cy0"1) incy1=.(2-~/\])@:(cy1"1) incy2=.(2-~/\])@:(cy2"1) incy3=.(2-~/\])@:(cy3"1) veloc=.(((incy0^2:)+incy1^2:)^0.5"_)%incy3 pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) pdDT pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) pdDT NB. ... incy2 is zero ... pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) pdDU pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) pdDU NB. ... incy2 is zero ... pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) pdDV pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) pdDV NB. ... incy2 is zero ... pd (}:@:(cy3"1);veloc)@:(rfyc"1)@:((S00xL{])"2) pdDW pd (}:@:(cy3"1);incy2)@:(rfyc"1)@:((S00xL{])"2) pdDW NB. ... incy2 is zero ... pd 'pdf' pd 'show'
Download: File:LTER2100E.txt
Download: File:LTER2101D.txt