User:Tom Allen/TensorExperimentsB09
Jump to navigation
Jump to search
...
...
...
... 5.75/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 =: 5.75p1 , 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'
... geodesic ...
NB. ... execute (ijx) ...
NB. ... using pdDA from the plot above ...
$pdDA
32001 8 4
NB. ... x0 ...
A=:cx0"1@:((S00xL{])"2) pdDA
(5.75p1-<./) A
0.0469877
(>./-5.75p1"_) A
0.252318
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^_45)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'
... x0 ...
...
...
...
... 5.875/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 =: 5.875p1 , 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'
... geodesic ...
NB. ... execute (ijx) ...
NB. ... using pdDA from the plot above ...
$pdDA
32001 8 4
NB. ... x0 ...
A=:cx0"1@:((S00xL{])"2) pdDA
(5.75p1-<./) A
0.301802
(>./-5.75p1"_) A
0.557733
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.@(27600"_){]) pdDA
pd 'pdf'
pd 'show'
... x0 ...
...
...
...
... 5.75/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 =: 5.75p1 , 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'
... geodesic ...
NB. ... execute (ijx) ...
NB. ... using pdDA from the plot above ...
$pdDA
32001 8 4
NB. ... x0 ...
A=:cx0"1@:((S00xL{])"2) pdDA
(5.75p1-<./) A
0.145307
(>./-5.75p1"_) A
0.216068
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'
... x0 ...
load 'plot'
NB. ... using pdDA from the plot above ...
pd (-@cx2"1;cx3"1)@:((S00xL{])"2) (i.@(32000"_){]) pdDA
pd 'pdf'
pd 'show'
... x3 ...
Download MoinMoin source: File:TensorExperimentsB09.ijs