User:Tom Allen/TensorExperimentsZ09

...

... E1cn (geodesics) ...

...

...

... one sign ...

...

load 'plot'

vSgDd2M=.[([(0{])@([(xdz smx"2 1 E1cn)zRef)hkxM),:@(S03xH{])
vSgDd2 =.gXZ(-@(([([(0{])@ch2k xRef@hkxM),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)]

Initds2=.(gXZ([(0{])@mcv xRef@hkxM),:@Initx)+/^:2@:*(]*/])@Initdx

gXcw  =:  1  0  _  _  _  _"_
gXbAll=: _1  1  1  1  1  0"_


Initx =: 2p1 , 0.5p1 ,  1p1  , 0.5p1"_
Initdx=: 0   , 0     , _0.01 , 0    "_

vS00xL   =:Initx
vS01dxL  =:0,0,0,0:
vS02d1L  =:Initd1
vS06ds   =:Initds
vS07count=:0,628"_
pd vSgDplot pdDA=:vSgD@vS''


Initx =: 3p1 , 0.5p1 ,  1p1  , 0.5p1"_
Initdx=: 0   , 0     , _0.01 , 0    "_

vS00xL   =:Initx
vS01dxL  =:0,0,0,0:
vS02d1L  =:Initd1
vS06ds   =:Initds
vS07count=:0,628"_
pd vSgDplot pdDB=:vSgD@vS''


Initx =: 4p1 , 0.5p1 ,  1p1  , 0.5p1"_
Initdx=: 0   , 0     , _0.01 , 0    "_

vS00xL   =:Initx
vS01dxL  =:0,0,0,0:
vS02d1L  =:Initd1
vS06ds   =:Initds
vS07count=:0,628"_
pd vSgDplot pdDC=:vSgD@vS''


pd 'aspect 1'
pd 'pdf'
pd 'show'

...

NB. ... execute (ijx) ...

NB. ... for pdDA ...

   $pdDA
629 8 4

NB. ... x0 ...

   *./ (0=]) (2p1-])    cx0"1@:((S00xL{])"2) pdDA
1

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=]) (0.50p1-]) cx3"1@:((S00xL{])"2) pdDA
1

NB. ... for pdDB ...

   $pdDB
629 8 4

NB. ... etc ...

---

Download MoinMoin source: File:TensorExperimentsZ09.ijs

---