User:Tom Allen/TensorExperimentsR11
Jump to navigation
Jump to search
load 'plot'
gXRAll=:,:@(0 1 2 8 0 0 0 0 0 0"_)
xCpts=:_ 0.01 0.55 1000,_ 0.5p1 0.5p1 0,_ 0 0 0,:_ 0 0 0"_
pd (gXR({.@(0|:]);((< 0 0 0 0){])@cfXX4@rfaV05XX"_ 1)]) pTtotal xCpts''
pd (gXR({.@(0|:]);((<0 0 0 0 0){])@cfXX5@rfaV05XX"_ 1)]) pTtotal xCpts''
pd 'pdf'
pd 'show'
... the derivative of blue is 0 at 0 ...
NB. ... script experimentalR.ijs ... rfaV21=:rfaV02*rfaV04 rfaV21XX=:[cfMP4 rfaV02XX,:rfaV04XX
NB. ... script experimentalR.ijs ... rfaV22=:gXRa04*cx0*cos@cx1 rfaV22XX=:[cfMP4 gXRa04XX,:[cfMP4 cx0XX,:[cfFN4(cos@cx1;-@sin@cx1;-@cos@cx1;sin@cx1;cos@cx1;-@sin@cx1),:cx1XX
NB. ... script experimentalR.ijs ... rfaV23=:rfaV22^gXRa05 rfaV22U0=: rfaV22^gXRa05 rfaV22U1=: gXRa05*rfaV22^gXRa05+_1: rfaV22U2=: (gXRa05+_1:)*gXRa05*rfaV22^gXRa05+_2: rfaV22U3=: (gXRa05+_2:)*(gXRa05+_1:)*gXRa05*rfaV22^gXRa05+_3: rfaV22U4=: (gXRa05+_3:)*(gXRa05+_2:)*(gXRa05+_1:)*gXRa05*rfaV22^gXRa05+_4: rfaV22U5=:(gXRa05+_4:)*(gXRa05+_3:)*(gXRa05+_2:)*(gXRa05+_1:)*gXRa05*rfaV22^gXRa05+_5: rfaV23XX=:[cfFN4(rfaV22U0;rfaV22U1;rfaV22U2;rfaV22U3;rfaV22U4;rfaV22U5),:rfaV22XX
NB. ... script experimentalR.ijs ... rfaV24=:^@rfaV23 rfaV24XX=:[cfFN4(^@rfaV23;^@rfaV23;^@rfaV23;^@rfaV23;^@rfaV23;^@rfaV23),:rfaV23XX
NB. ... script experimentalR.ijs ... rfaV25=:^@-@rfaV23 rfaV25XX=:[cfFN4(^@-@rfaV23;-@^@-@rfaV23;^@-@rfaV23;-@^@-@rfaV23;^@-@rfaV23;-@^@-@rfaV23),:rfaV23XX
NB. ... script experimentalR.ijs ... rfaV26=:rfaV24+rfaV25 rfaV26XX=:rfaV24XX+&.>rfaV25XX
NB. ... script experimentalR.ijs ... rfaV27=:rfaV26^_1: rfaV26U0=: rfaV26^_1: rfaV26U1=: _1*rfaV26^_2: rfaV26U2=: _2*_1*rfaV26^_3: rfaV26U3=: _3*_2*_1*rfaV26^_4: rfaV26U4=: _4*_3*_2*_1*rfaV26^_5: rfaV26U5=:_5*_4*_3*_2*_1*rfaV26^_6: rfaV27XX=:[cfFN4(rfaV26U0;rfaV26U1;rfaV26U2;rfaV26U3;rfaV26U4;rfaV26U5),:rfaV26XX
load 'plot'
pd (];>@(0{])@rfzM1XX"0) pTsteps _ _5 5 1000
pd (];>@(0{])@rfzM2XX"0) pTsteps _ _5 5 1000
pd (];>@(0{])@rfzM3XX"0) pTsteps _ _5 5 1000
pd (];>@(0{])@rfzM4XX"0) pTsteps _ _5 5 1000
pd (];>@(0{])@rfzM5XX"0) pTsteps _ _5 5 1000
gXRAll=:,:@(0 0 0 0 1 1 0 0 0 0"_)
pd ((cx0*cos@cx1)"1;gXR cfXX0@rfaV27XX"_ 1]) pTtotal@(_ 5 5 0,_ 0.99p1 0.01p1 1000,_ 0 0 0,:_ 0 0 0"_)'' NB. ... overwrite ...
gXRAll=:,:@(0 0 0 0 1 2 0 0 0 0"_)
pd ((cx0*cos@cx1)"1;gXR cfXX0@rfaV27XX"_ 1]) pTtotal@(_ 5 5 0,_ 0.99p1 0.01p1 1000,_ 0 0 0,:_ 0 0 0"_)'' NB. ... overwrite ...
gXRAll=:,:@(0 0 0 0 1 3 0 0 0 0"_)
pd ((cx0*cos@cx1)"1;gXR cfXX0@rfaV27XX"_ 1]) pTtotal@(_ 5 5 0,_ 0.99p1 0.01p1 1000,_ 0 0 0,:_ 0 0 0"_)'' NB. ... overwrite ...
gXRAll=:,:@(0 0 0 0 1 4 0 0 0 0"_)
pd ((cx0*cos@cx1)"1;gXR cfXX0@rfaV27XX"_ 1]) pTtotal@(_ 5 5 0,_ 0.99p1 0.01p1 1000,_ 0 0 0,:_ 0 0 0"_)'' NB. ... overwrite ...
gXRAll=:,:@(0 0 0 0 1 5 0 0 0 0"_)
pd ((cx0*cos@cx1)"1;gXR cfXX0@rfaV27XX"_ 1]) pTtotal@(_ 2 2 0,_ 0.99p1 0.01p1 1000,_ 0 0 0,:_ 0 0 0"_)'' NB. ... overwrite ...
pd 'pdf'
pd 'show'
NB. ... execute (ijx) ... gXRAll=:,:@(0 0 0 0 1 1 0 0 0 0"_) xCpts=:_ 5 5 0,_ 0.99p1 0.01p1 1000,_ 0 0 0,:_ 0 0 0"_ (gXR((rfaV27-:cfXX0@rfaV27XX)"_ 1)]) ((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 (gXR(cfXX1@rfaV27XX"_ 1(((2^_44);2^_18)qteq[;])"1(0|:[:(gXR cfXX0@rfaV27XX])D.1])"_ 1)]) ((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 (gXR(cfXX2@rfaV27XX"_ 1(((2^_44);2^_17)qteq[;])"2(0|:[:(gXR cfXX1@rfaV27XX])D.1])"_ 1)]) ((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 (gXR(cfXX3@rfaV27XX"_ 1(((2^_44);2^_16)qteq[;])"3(0|:[:(gXR cfXX2@rfaV27XX])D.1])"_ 1)]) ((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 (gXR(cfXX4@rfaV27XX"_ 1(((2^_44);2^_13)qteq[;])"4(0|:[:(gXR cfXX3@rfaV27XX])D.1])"_ 1)]) ((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 0 0 0 1 1 0 0 0 0"_) xCpts=:_ 5 5 0,_ 0.99p1 0.01p1 1000,_ 0 0 0,:_ 0 0 0"_ A=:(gXR(cfXX1@([rfTT4 rfxXXS,:[rfaV27XX rfxc)"_ 1)]) ((0;50"_)rfyc"1@pTRandom xCpts)'' B=:(gXR((0|:[:(gXR cfXX0@([rfTT4 rfxXXS,:[rfaV27XX rfxc)])D.1])"_ 1)]) ((0;50"_)rfyc"1@pTRandom xCpts)'' A (((2^_28);2^_21)qteq[;])"1 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=:(gXR(cfXX2@([rfTT4 rfxXXS,:[rfaV27XX rfxc)"_ 1)]) ((0;50"_)rfyc"1@pTRandom xCpts)'' B=:(gXR((0|:[:(gXR cfXX1@([rfTT4 rfxXXS,:[rfaV27XX rfxc)])D.1])"_ 1)]) ((0;50"_)rfyc"1@pTRandom xCpts)'' A (((2^_29);2^_20)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 A=:(gXR(cfXX3@([rfTT4 rfxXXS,:[rfaV27XX rfxc)"_ 1)]) ((0;50"_)rfyc"1@pTRandom xCpts)'' B=:(gXR((0|:[:(gXR cfXX2@([rfTT4 rfxXXS,:[rfaV27XX rfxc)])D.1])"_ 1)]) ((0;50"_)rfyc"1@pTRandom xCpts)'' A (((2^_27);2^_17)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=:(gXR(cfXX4@([rfTT4 rfxXXS,:[rfaV27XX rfxc)"_ 1)]) ((0;50"_)rfyc"1@pTRandom xCpts)'' B=:(gXR((0|:[:(gXR cfXX3@([rfTT4 rfxXXS,:[rfaV27XX rfxc)])D.1])"_ 1)]) ((0;50"_)rfyc"1@pTRandom xCpts)'' A (((2^_26);2^_18)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
Download: File:LTER1100C.txt
Download: File:LTER1101C.txt
Download: File:LTER1102C.txt
Download: File:LTER1103C.txt
Download: File:LTER1104C.txt
Download: File:LTER1105C.txt
Download: File:LTER1106C.txt
Download: File:LTER1107D.txt