User:Tom Allen/TensorExperimentsZ19

gXsAll=. _1  1  1  1  1  0"_

Initu =.  _2.1095p1 , 0        , 0 , 0.5p1"_
Initdu=. ( 0        , 0.013685 , 0 , 0.1  "_)*%@(1"_)

Initds=.%:@|@((gXS([(0{])@mcv uRef@hkuR),:@Initu)+/^:2@:*(]*/])@Initdu)

vS00xL   =.Initu
vS01dxL  =.0,0,0,0:
vS02d1L  =.Initdu%Initds
vS06ds   =.Initds
vS07count=.(0,96000"_)*1"_

vSgDd2=.gXS([-@([(0{])@E1cn uRef@hkuR),:@(S03xH{]))]

pdDR=:vSgD@vS''

NB. ... pdDR to disk file ...

pdDR arrayPut   'pdDR'
     arrayClose 'pdDR'

NB. ... establish pdDR if necessary ...

                arrayOpen  'pdDR'
pdDR=:96001 8 4 arrayGet   'pdDR'   NB. ... shape required ...
                arrayClose 'pdDR'

load 'plot'

pd ((2.1095p1*cos@]"0);2.1095p1*sin@]"0) pTsteps _ 0p1 2p1 160

pd <"1@(0|:])@:((cu0,cu1)@(S00xL{])"2) pdDR

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

... waves ...

---

Download MoinMoin source: File:TensorExperimentsZ19.ijs

---