User:Tom Allen/ta/AddendumX02P02
< User:Tom Allen | ta
Jump to navigation
Jump to search
vZ0=:0.05 1 2 8 1 1 1 _1 0,:0"_ vSgDd2=:vZ0([([(0{])@E1cn zRef@([hkyX02C vZyc"_ 1)),:@(S03yH{]))] Initx =: 0.7500 , 0.5p1 , 0.25p1 , 0p1"_ Initdx=: (0 , 0 , 0 , 0.1"_)*%@InitInt=:8"_ Inity =:zfyc@Initx Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx Initz =:vZ0 vZzc Inity Initdz=:vZ0((vZgamma*(0{])+-@vZRv00*(3{])),(1{]),(2{]),vZgamma*(3{])+-@vZRv00*(0{]))Initdy Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy vS00yL =:Initz vS01dyL =:0,0,0,0: vS02d1L =:Initdz%Initds vS06ds =:Initds vS07count=:(0,10000"_)*InitInt pdDS1A =:vSgD@vS'' pdDS1A arrayPut 'pdDS1A' arrayClose 'pdDS1A'
vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.01"_ vSgDd2=:vZ0([([(0{])@E1cn zRef@([hkyX02C vZyc"_ 1)),:@(S03yH{]))] Initx =: 0.7500 , 0.5p1 , 0.25p1 , 0p1"_ Initdx=: (0 , 0 , 0 , 0.1"_)*%@InitInt=:8"_ Inity =:zfyc@Initx Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx Initz =:vZ0 vZzc Inity Initdz=:vZ0((vZgamma*(0{])+-@vZRv00*(3{])),(1{]),(2{]),vZgamma*(3{])+-@vZRv00*(0{]))Initdy Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy vS00yL =:Initz vS01dyL =:0,0,0,0: vS02d1L =:Initdz%Initds vS06ds =:Initds vS07count=:(0,10000"_)*InitInt pdDS1B =:vSgD@vS'' pdDS1B arrayPut 'pdDS1B' arrayClose 'pdDS1B'
vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.90"_ vSgDd2=:vZ0([([(0{])@E1cn zRef@([hkyX02C vZyc"_ 1)),:@(S03yH{]))] Initx =: 0.7500 , 0.5p1 , 0.25p1 , 0p1"_ Initdx=: (0 , 0 , 0 , 0.1"_)*%@InitInt=:8"_ Inity =:zfyc@Initx Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx Initz =:vZ0 vZzc Inity Initdz=:vZ0((vZgamma*(0{])+-@vZRv00*(3{])),(1{]),(2{]),vZgamma*(3{])+-@vZRv00*(0{]))Initdy Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy vS00yL =:Initz vS01dyL =:0,0,0,0: vS02d1L =:Initdz%Initds vS06ds =:Initds vS07count=:(0,10000"_)*InitInt pdDS1C =:vSgD@vS'' pdDS1C arrayPut 'pdDS1C' arrayClose 'pdDS1C'
vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.95"_ vSgDd2=:vZ0([([(0{])@E1cn zRef@([hkyX02C vZyc"_ 1)),:@(S03yH{]))] Initx =: 0.7500 , 0.5p1 , 0.25p1 , 0p1"_ Initdx=: (0 , 0 , 0 , 0.1"_)*%@InitInt=:8"_ Inity =:zfyc@Initx Initdy=:zfyc@(Initx+Initdx)+-@zfyc@Initx Initz =:vZ0 vZzc Inity Initdz=:vZ0((vZgamma*(0{])+-@vZRv00*(3{])),(1{]),(2{]),vZgamma*(3{])+-@vZRv00*(0{]))Initdy Initds=:|@((((0{])^2:)+((1{])^2:)+((2{])^2:)+-@((3{])^2:))^0.5"_)@Initdy vS00yL =:Initz vS01dyL =:0,0,0,0: vS02d1L =:Initdz%Initds vS06ds =:Initds vS07count=:(0,10000"_)*InitInt pdDS1D =:vSgD@vS'' pdDS1D arrayPut 'pdDS1D' arrayClose 'pdDS1D'
NB. ... execute (ijx) ... arrayOpen 'pdDS1A' pdDS1A=:80001 8 4 arrayGet 'pdDS1A' NB. ... shape required ... arrayClose 'pdDS1A' arrayOpen 'pdDS1B' pdDS1B=:80001 8 4 arrayGet 'pdDS1B' NB. ... shape required ... arrayClose 'pdDS1B' arrayOpen 'pdDS1C' pdDS1C=:80001 8 4 arrayGet 'pdDS1C' NB. ... shape required ... arrayClose 'pdDS1C' arrayOpen 'pdDS1D' pdDS1D=:80001 8 4 arrayGet 'pdDS1D' NB. ... shape required ... arrayClose 'pdDS1D'
load 'plot' pd 'reset' vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.01"_ pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1A NB. ... blue ... pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1B NB. ... red ... pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@(vZ0 vZyc"_ 1])@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1B NB. ... green overwrites blue ... pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@(vZ0 vZzc"_ 1])@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1A NB. ... purple overwrites red ... pd 'pdf'
load 'plot' pd 'reset' incz0=:(2-~/\])@:(cz0"1) incz1=:(2-~/\])@:(cz1"1) incz2=:(2-~/\])@:(cz2"1) incz3=:(2-~/\])@:(cz3"1) veloc=:(((incz0^2:)+incz1^2:)^0.5"_)%incz3 pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1A pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) (i.@(8001"_){]) pdDS1B pd 'pdf'
load 'plot' pd 'reset' pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2)) pdDS1A NB. ... blue ... pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2)) pdDS1C NB. ... red ... pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@:((S00yL{])"2)) pdDS1D NB. ... green ... vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.90"_ pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@(vZ0 vZzc"_ 1])@:((S00yL{])"2)) pdDS1A NB. ... purple overwrites red ... vZ0=:0.05 1 2 8 1 1 1 _1 0,:_0.95"_ pd (((S06ds{])@{:*i.@{.@$)@];cz0"1@(vZ0 vZzc"_ 1])@:((S00yL{])"2)) pdDS1A NB. ... pink overwrites green ... pd 'pdf'
load 'plot' pd 'reset' incz0=:(2-~/\])@:(cz0"1) incz1=:(2-~/\])@:(cz1"1) incz2=:(2-~/\])@:(cz2"1) incz3=:(2-~/\])@:(cz3"1) veloc=:(((incz0^2:)+incz1^2:)^0.5"_)%incz3 pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) pdDS1A pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) pdDS1C pd (}:@((S06ds{])@{:*i.@{.@$)@];veloc@:((S00yL{])"2)) pdDS1D pd 'pdf'
Next Page: User:Tom Allen/ta/AddendumX02P03
Prev Page: User:Tom Allen/ta/AddendumX02P01