User:Tom Allen/TensorExperimentsR07
Jump to navigation
Jump to search
NB. ... script experimentalR.ijs ... gXR=:gXRAll
NB. ... script experimentalR.ijs ...
gXRa00=:(<0 0){[
gXRa01=:(<0 1){[
gXRa02=:(<0 2){[
gXRa03=:(<0 3){[
gXRa04=:(<0 4){[
gXRa05=:(<0 5){[
gXRa06=:(<0 6){[
gXRa07=:(<0 7){[
gXRa08=:(<0 8){[
gXRa09=:(<0 9){[
gXRa00XX=:gXRa00;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXRa01XX=:gXRa01;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXRa02XX=:gXRa02;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXRa03XX=:gXRa03;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXRa04XX=:gXRa04;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXRa05XX=:gXRa05;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXRa06XX=:gXRa06;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXRa07XX=:gXRa07;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXRa08XX=:gXRa08;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXRa09XX=:gXRa09;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
gXR400XX=:0;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);_"_
gXR500XX=:0;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
NB. ... execute (ijx) ... gXRAll=:,:@(0 1 2 3 4 5 6 7 8 9"_) (gXR(gXRa00,gXRa01,gXRa02,gXRa03,gXRa04,gXRa05,gXRa06,gXRa07,gXRa08,gXRa09)])'' 0 1 2 3 4 5 6 7 8 9
NB. ... script experimentalR.ijs ... cx0XX=:cx0;1 0 0 0;((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0: cx1XX=:cx1;0 1 0 0;((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0: cx2XX=:cx2;0 0 1 0;((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0: cx3XX=:cx3;0 0 0 1;((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
NB. ... script experimentalR.ijs ... rfaV00=:gXRa01*cx0*sin@cx1 rfaV00XX=:[cfMP5 gXRa01XX,:[cfMP5 cx0XX,:[cfFN5(sin@cx1;cos@cx1;-@sin@cx1;-@cos@cx1;sin@cx1;cos@cx1),:cx1XX
NB. ... script experimentalR.ijs ... rfaV01=:rfaV00^gXRa03 rfaV00U0=: rfaV00^gXRa03 rfaV00U1=: gXRa03*rfaV00^gXRa03+_1: rfaV00U2=: (gXRa03+_1:)*gXRa03*rfaV00^gXRa03+_2: rfaV00U3=: (gXRa03+_2:)*(gXRa03+_1:)*gXRa03*rfaV00^gXRa03+_3: rfaV00U4=: (gXRa03+_3:)*(gXRa03+_2:)*(gXRa03+_1:)*gXRa03*rfaV00^gXRa03+_4: rfaV00U5=:(gXRa03+_4:)*(gXRa03+_3:)*(gXRa03+_2:)*(gXRa03+_1:)*gXRa03*rfaV00^gXRa03+_5: rfaV01XX=:[cfFN5(rfaV00U0;rfaV00U1;rfaV00U2;rfaV00U3;rfaV00U4;rfaV00U5),:rfaV00XX
NB. ... script experimentalR.ijs ... rfaV02=:1-^@-@rfaV01 rfaV02XX=:[cfFN5((1+-@^@-@rfaV01);^@-@rfaV01;-@^@-@rfaV01;^@-@rfaV01;-@^@-@rfaV01;^@-@rfaV01),:rfaV01XX
NB. ... script experimentalR.ijs ... rfaV03=:rfaV00^gXRa02 rfaV00xU0=: rfaV00^gXRa02 rfaV00xU1=: gXRa02*rfaV00^gXRa02+_1: rfaV00xU2=: (gXRa02+_1:)*gXRa02*rfaV00^gXRa02+_2: rfaV00xU3=: (gXRa02+_2:)*(gXRa02+_1:)*gXRa02*rfaV00^gXRa02+_3: rfaV00xU4=: (gXRa02+_3:)*(gXRa02+_2:)*(gXRa02+_1:)*gXRa02*rfaV00^gXRa02+_4: rfaV00xU5=:(gXRa02+_4:)*(gXRa02+_3:)*(gXRa02+_2:)*(gXRa02+_1:)*gXRa02*rfaV00^gXRa02+_5: rfaV03XX=:[cfFN5(rfaV00xU0;rfaV00xU1;rfaV00xU2;rfaV00xU3;rfaV00xU4;rfaV00xU5),:rfaV00XX
NB. ... script experimentalR.ijs ... rfaV04=:rfaV03^_1: rfaV03U0=: rfaV03^_1: rfaV03U1=: _1*rfaV03^_2: rfaV03U2=: _2*_1*rfaV03^_3: rfaV03U3=: _3*_2*_1*rfaV03^_4: rfaV03U4=: _4*_3*_2*_1*rfaV03^_5: rfaV03U5=:_5*_4*_3*_2*_1*rfaV03^_6: rfaV04XX=:[cfFN5(rfaV03U0;rfaV03U1;rfaV03U2;rfaV03U3;rfaV03U4;rfaV03U5),:rfaV03XX
NB. ... script experimentalR.ijs ... rfaV05=:rfaV02*rfaV04 rfaV05XX=:[cfMP5 rfaV02XX,:rfaV04XX
NB. ... execute (ijx) ... gXRAll=:,:@(0 1 2 9 0 0 0 0 0 0"_) xCpts=:_ 0.7 1.3 26,_ 0.4p1 0.6p1 26,_ 0 2p1 26,:_ 0 20 26"_ (gXR((rfaV05-:cfXX0@rfaV05XX)"_ 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@rfaV05XX"_ 1(((2^_44);2^_17)qteq[;])"1(0|:[:(gXR cfXX0@rfaV05XX])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@rfaV05XX"_ 1(((2^_44);2^_15)qteq[;])"2(0|:[:(gXR cfXX1@rfaV05XX])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@rfaV05XX"_ 1(((2^_44);2^_14)qteq[;])"3(0|:[:(gXR cfXX2@rfaV05XX])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@rfaV05XX"_ 1(((2^_44);2^_14)qteq[;])"4(0|:[:(gXR cfXX3@rfaV05XX])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(cfXX5@rfaV05XX"_ 1(((2^_44);2^_10)qteq[;])"5(0|:[:(gXR cfXX4@rfaV05XX])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
load 'plot'
gXRAll=:,:@(0 1 2 9 0 0 0 0 0 0"_)
xCpts=:_ 1e_16 2.5 1000,_ 0.5p1 0.5p1 0,_ 0 0 0,:_ 0 0 0"_
pd (gXR({.@(0|:]); cfXX0@rfaV05XX"_ 1)]) pTtotal xCpts''
pd (gXR({.@(0|:]);(0{])@cfXX1@rfaV05XX"_ 1)]) pTtotal xCpts'' NB. ... 1st derivative ...
NB. ... overwrite ...
pd (];>@(0{])@rffm9XX"0) pTsteps _ 1e_16 2.5 1000
pd (];>@(1{])@rffm9XX"0) pTsteps _ 1e_16 2.5 1000 NB. ... 1st derivative ...
pd 'pdf'
pd 'show'
load 'plot'
gXRAll=:,:@(0 1 2 9 0 0 0 0 0 0"_)
xCpts=:_ 0.7 2 1000,_ 0.5p1 0.5p1 0,_ 0 0 0,:_ 0 0 0"_
pd (gXR({.@(0|:]);{.@(0|:])^_2: )]) pTtotal xCpts''
pd (gXR({.@(0|:]);cfXX0@rfaV05XX"_ 1)]) pTtotal xCpts''
pd 'pdf'
pd 'show'
Download: File:LTER0700D.txt
Download: File:LTER0701C.txt
Download: File:LTER0702C.txt
Download: File:LTER0703C.txt
Download: File:LTER0704C.txt
Download: File:LTER0705D.txt