User:Tom Allen/ta/AddendumX06P03

From J Wiki
Jump to navigation Jump to search 


load 'plot'

pd 'reset'

pU0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_)

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.00p1 0.00p1 0"_)*10"_
pd 'itemcolor red'
pd (pU0(cx0"_ 1;[(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.15p1 0.15p1 0"_)*10"_
pd 'itemcolor green'
pd (pU0(cx0"_ 1;[(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.35p1 0.35p1 0"_)*10"_
pd 'itemcolor purple'
pd (pU0(cx0"_ 1;[(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.50p1 0.50p1 0"_)*10"_
pd 'itemcolor pink'
pd (pU0(cx0"_ 1;[(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.65p1 0.65p1 0"_)*10"_
pd 'itemcolor purple'
pd (pU0(cx0"_ 1;[(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.85p1 0.85p1 0"_)*10"_
pd 'itemcolor green'
pd (pU0(cx0"_ 1;[(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 1.00p1 1.00p1 0"_)*10"_
pd 'itemcolor red'
pd (pU0(cx0"_ 1;[(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

pd 'itemcolor blue'
pd (((0.8704"_)@]"0);]"0) pTsteps _ _4 4 160
pd (((1.3235"_)@]"0);]"0) pTsteps _ _4 4 160

pd 'itemcolor black'
pd (((0.9675"_)@]"0);]"0) pTsteps _ _9 7 160
pd (((1.2832"_)@]"0);]"0) pTsteps _ _9 7 160

pd 'canvas plot.html 666 500'

'finished' 1!:2[2


JX06P0300C.png 


pU0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_)

vSgDd2=:pU0([([(0{])@E1cn yRef@hkyX06B),:@(S03yH{]))]


Initx =:  0.8700 , 0.5p1 , 0.25p1 , 0p1"_
Initdx=: (0      , 0     , 0      , 0.1"_)*%@InitInt=:32"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,1200"_)*InitInt
pdDR     =:vSgD@vS''

pdDR arrayPut   'pdDR'
     arrayClose 'pdDR'


Initx =:  0.5500 , 0.5p1 , 0.25p1 , 0p1"_
Initdx=: (0      , 0     , 0      , 0.1"_)*%@InitInt=:32"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,1200"_)*InitInt
pdDRlow     =:vSgD@vS''

pdDRlow arrayPut   'pdDRlow'
        arrayClose 'pdDRlow'


Initx =:  0.9675 , 0.5p1 , 1.25p1 , 0p1"_
Initdx=: (0      , 0     , 0      , 0.1"_)*%@InitInt=:32"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,1200"_)*InitInt
pdDRhigh     =:vSgD@vS''

pdDRhigh arrayPut   'pdDRhigh'
         arrayClose 'pdDRhigh'


Initx =:  1.3200 , 0.5p1 , 0.25p1 , 0p1"_
Initdx=: (0      , 0     , 0      , 0.1"_)*%@InitInt=:32"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,1200"_)*InitInt
pdDS     =:vSgD@vS''

pdDS arrayPut   'pdDS'
     arrayClose 'pdDS'


Initx =:  1.2832 , 0.5p1 , 0.25p1 , 0p1"_
Initdx=: (0      , 0     , 0      , 0.1"_)*%@InitInt=:32"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,1200"_)*InitInt
pdDSlow     =:vSgD@vS''

pdDSlow arrayPut   'pdDSlow'
        arrayClose 'pdDSlow'


Initx =:  1.4500 , 0.5p1 , 1.25p1 , 0p1"_
Initdx=: (0      , 0     , 0      , 0.1"_)*%@InitInt=:32"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,1200"_)*InitInt
pdDShigh     =:vSgD@vS''

pdDShigh arrayPut   'pdDShigh'
         arrayClose 'pdDShigh'

'finished' 1!:2[2



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

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

                    arrayOpen  'pdDRlow'
pdDRlow =:38401 8 4 arrayGet   'pdDRlow'    NB. ... shape required ...
                    arrayClose 'pdDRlow'

                    arrayOpen  'pdDRhigh'
pdDRhigh=:38401 8 4 arrayGet   'pdDRhigh'   NB. ... shape required ...
                    arrayClose 'pdDRhigh'

                    arrayOpen  'pdDS'
pdDS    =:38401 8 4 arrayGet   'pdDS'       NB. ... shape required ...
                    arrayClose 'pdDS'

                    arrayOpen  'pdDSlow'
pdDSlow =:38401 8 4 arrayGet   'pdDSlow'    NB. ... shape required ...
                    arrayClose 'pdDSlow'

                    arrayOpen  'pdDShigh'
pdDShigh=:38401 8 4 arrayGet   'pdDShigh'   NB. ... shape required ...
                    arrayClose 'pdDShigh'

'finished' 1!:2[2



load 'plot'

pd 'reset'

pd 'itemcolor blue'
pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDR
pd 'itemcolor red'
pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDRlow
pd 'itemcolor green'
pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDRhigh
pd 'itemcolor blue'
pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDS
pd 'itemcolor purple'
pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDSlow
pd 'itemcolor pink'
pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDShigh

pd 'itemcolor Black'
pd ((0.9675*cos@]"0);0.9675*sin@]"0) pTsteps _ 0p1 2p1 160
pd ((1.2832*cos@]"0);1.2832*sin@]"0) pTsteps _ 0p1 2p1 160

pd 'aspect 1'
pd 'canvas plot.html 666 500'

'finished' 1!:2[2


JX06P0301C.png


GX06P0300C.png 


load 'plot'

pd 'reset'

incy0=:(2-~/\])@:(cy0"1)
incy1=:(2-~/\])@:(cy1"1)
incy2=:(2-~/\])@:(cy2"1)
incy3=:(2-~/\])@:(cy3"1)
veloc=:(((incy0^2:)+incy1^2:)^0.5"_)%incy3

pd (}:@((0.1%32"_)*i.@(38401"_));veloc)@:((S00yL{])"2) (i.@(38401"_){]) pdDR
pd (}:@((0.1%32"_)*i.@(38401"_));veloc)@:((S00yL{])"2)                  pdDRlow
pd (}:@((0.1%32"_)*i.@(38401"_));veloc)@:((S00yL{])"2)                  pdDRhigh

pd 'canvas plot.html 666 500'

'finished' 1!:2[2


JX06P0302C.png 


load 'plot'

pd 'reset'

incy0=:(2-~/\])@:(cy0"1)
incy1=:(2-~/\])@:(cy1"1)
incy2=:(2-~/\])@:(cy2"1)
incy3=:(2-~/\])@:(cy3"1)
veloc=:(((incy0^2:)+incy1^2:)^0.5"_)%incy3

pd (}:@((0.1%32"_)*i.@(38401"_));veloc)@:((S00yL{])"2) (i.@(38401"_){]) pdDS
pd (}:@((0.1%32"_)*i.@(38401"_));veloc)@:((S00yL{])"2)                  pdDSlow
pd (}:@((0.1%32"_)*i.@(38401"_));veloc)@:((S00yL{])"2)                  pdDShigh

pd 'canvas plot.html 666 500'

'finished' 1!:2[2


JX06P0303C.png 



load 'plot'

pd 'reset'

pU0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_)

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.00p1 0.00p1 0"_)*10"_
pd 'itemcolor red'
pd (pU0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.15p1 0.15p1 0"_)*10"_
pd 'itemcolor green'
pd (pU0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.35p1 0.35p1 0"_)*10"_
pd 'itemcolor purple'
pd (pU0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.50p1 0.50p1 0"_)*10"_
pd 'itemcolor pink'
pd (pU0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.65p1 0.65p1 0"_)*10"_
pd 'itemcolor purple'
pd (pU0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.85p1 0.85p1 0"_)*10"_
pd 'itemcolor green'
pd (pU0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

xCpts=:_ 0.4 1.7 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 1.00p1 1.00p1 0"_)*10"_
pd 'itemcolor red'
pd (pU0(cx0"_ 1;[-@(0{])"1@E1cn xRef@([hkyX06B zfyc"_ 1))]) pTtotal xCpts''

pd 'itemcolor blue'
pd (((1.0909"_)@]"0);]"0) pTsteps _ _4 4 160

pd 'itemcolor black'
pd (((0.9675"_)@]"0);]"0) pTsteps _ _9 7 160
pd (((1.2832"_)@]"0);]"0) pTsteps _ _9 7 160

pd 'canvas plot.html 666 500'

'finished' 1!:2[2


JX06P0304C.png 


pU0=:,:@(0.05 1 2 8 1 1 1 _0.1 0"_)

vSgDd2=:pU0([([-@(0{])@E1cn yRef@hkyX06B),:@(S03yH{]))]


Initx =:  1.0911 , 0.5p1 , 0.25p1 , 0p1"_
Initdx=: (0      , 0     , 0      , 0.1"_)*%@InitInt=:32"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,1200"_)*InitInt
pdDT     =:vSgD@vS''

pdDT arrayPut   'pdDT'
     arrayClose 'pdDT'


Initx =:  0.9675 , 0.5p1 , 0.25p1 , 0p1"_
Initdx=: (0      , 0     , 0      , 0.1"_)*%@InitInt=:32"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,1200"_)*InitInt
pdDTlow     =:vSgD@vS''

pdDTlow arrayPut   'pdDTlow'
        arrayClose 'pdDTlow'


Initx =:  1.2832 , 0.5p1 , 1.25p1 , 0p1"_
Initdx=: (0      , 0     , 0      , 0.1"_)*%@InitInt=:32"_

Inity =:zfyc@Initx
Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx

Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy

vS00yL   =:Inity
vS01dyL  =:0,0,0,0:
vS02d1L  =:Initdy%Initds
vS06ds   =:Initds
vS07count=:(0,1200"_)*InitInt
pdDThigh     =:vSgD@vS''

pdDThigh arrayPut   'pdDThigh'
         arrayClose 'pdDThigh'

'finished' 1!:2[2



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

                    arrayOpen  'pdDT'
pdDT    =:38401 8 4 arrayGet   'pdDT'       NB. ... shape required ...
                    arrayClose 'pdDT'

                    arrayOpen  'pdDTlow'
pdDTlow =:38401 8 4 arrayGet   'pdDTlow'    NB. ... shape required ...
                    arrayClose 'pdDTlow'

                    arrayOpen  'pdDThigh'
pdDThigh=:38401 8 4 arrayGet   'pdDThigh'   NB. ... shape required ...
                    arrayClose 'pdDThigh'

'finished' 1!:2[2



load 'plot'

pd 'reset'

pd 'itemcolor blue'
pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDT
pd 'itemcolor red'
pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDTlow
pd 'itemcolor green'
pd <"1@(0|:])@:((cy0,cy1)@(S00yL{])"2) pdDThigh

pd 'itemcolor Black'
pd ((0.9675*cos@]"0);0.9675*sin@]"0) pTsteps _ 0p1 2p1 160
pd ((1.2832*cos@]"0);1.2832*sin@]"0) pTsteps _ 0p1 2p1 160

pd 'aspect 1'
pd 'canvas plot.html 666 500'

'finished' 1!:2[2


JX06P0305.png


GX06P0301C.png 


load 'plot'

pd 'reset'

incy0=:(2-~/\])@:(cy0"1)
incy1=:(2-~/\])@:(cy1"1)
incy2=:(2-~/\])@:(cy2"1)
incy3=:(2-~/\])@:(cy3"1)
veloc=:(((incy0^2:)+incy1^2:)^0.5"_)%incy3

pd (}:@((0.1%32"_)*i.@(38401"_));veloc)@:((S00yL{])"2) (i.@(38401"_){]) pdDT
pd (}:@((0.1%32"_)*i.@(38401"_));veloc)@:((S00yL{])"2)                  pdDTlow
pd (}:@((0.1%32"_)*i.@(38401"_));veloc)@:((S00yL{])"2)                  pdDThigh

pd 'canvas plot.html 666 500'

'finished' 1!:2[2


JX06P0306C.png


Download: File:LX06P0300C.txt
Download: File:LX06P0301C.txt


Next Page: User:Tom Allen/ta/AddendumX06P03A
Prev Page: User:Tom Allen/ta/AddendumX06P02A

Retrieved from "https://code.jsoftware.com/mediawiki/index.php?title=User:Tom_Allen/ta/AddendumX06P03&oldid=49125"