User:Tom Allen/TensorExperimentsZ20
Jump to navigation
Jump to search
gXsAll=. _1 1 1 1 1 0"_ Initu =. _2.05p1 , 0 , 0 , 0.5p1"_ Initdu=. ( 0 , 0.01331 , 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,3724"_)*1"_ vSgDd2=.gXS([-@([(0{])@E1cn uRef@hkuR),:@(S03xH{]))] pdDA=:vSgD@vS'' NB. ... pdDA to disk file ... pdDA arrayPut 'pdDA' arrayClose 'pdDA' Initu =. _2.20p1 , 0 , 0 , 0.5p1"_ Initdu=. ( 0 , 0.01375 , 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,3362"_)*1"_ vSgDd2=.gXS([-@([(0{])@E1cn uRef@hkuR),:@(S03xH{]))] pdDB=:vSgD@vS'' NB. ... pdDB to disk file ... pdDB arrayPut 'pdDB' arrayClose 'pdDB'
load 'plot' pd ((2.05p1*cos@]"0);2.05p1*sin@]"0) pTsteps _ 0p1 2p1 160 pd <"1@(0|:])@:((cu0,cu1)@(S00xL{])"2) pdDA pd ((2.20p1*cos@]"0);2.20p1*sin@]"0) pTsteps _ 0p1 2p1 160 pd <"1@(0|:])@:((cu0,cu1)@(S00xL{])"2) pdDB pd 'aspect 1' pd 'pdf' pd 'show'
Download MoinMoin source: File:TensorExperimentsZ20.ijs