User:Tom Allen/ta/AddendumX04P04

... these plots show the charge density at points in 4-dimensional spacetime ...

... red is high ...
... blue is low ...

load 'plot'

pd 'reset'

uQ0=:,:@(1 1 2 8 1 1 1 _0.5 0"_)

xCpts=:_ 0.4 1.7 30,_ 0.5p1 0.5p1 0,_ 0 2p1 30,:_ 0.00p1 0.00p1 0"_

pd 'viewpoint 0.001 0 100'
pd 'type surface'

pd (uQ0((cx0*cos@cx2)"_ 1;(cx0*sin@cx2)"_ 1;([-@(3{])"1@JEMcn yRef@([hkyX04A zfyc"_ 1))"_ 2)]) pTtotal3D xCpts''

pd 'aspect 1'
pd 'pdf'

... charge density at time y3=0 ...

load 'plot'

pd 'reset'

uQ0=:,:@(1 1 2 8 1 1 1 _0.5 0"_)

xCpts=:_ 0.4 1.7 30,_ 0.5p1 0.5p1 0,_ 0 2p1 30,:_ 0.50p1 0.50p1 0"_

pd 'viewpoint 0.001 0 100'
pd 'type surface'

pd (uQ0((cx0*cos@cx2)"_ 1;(cx0*sin@cx2)"_ 1;([-@(3{])"1@JEMcn yRef@([hkyX04A zfyc"_ 1))"_ 2)]) pTtotal3D xCpts''

pd 'aspect 1'
pd 'pdf'

... charge density at time y3=0.50p1 ...

load 'plot'

pd 'reset'

uQ0=:,:@(1 1 2 8 1 1 1 _0.5 0"_)

xCpts=:_ 0.4 1.7 30,_ 0.5p1 0.5p1 0,_ 0 2p1 30,:_ 1.00p1 1.00p1 0"_

pd 'viewpoint 0.001 0 100'
pd 'type surface'

pd (uQ0((cx0*cos@cx2)"_ 1;(cx0*sin@cx2)"_ 1;([-@(3{])"1@JEMcn yRef@([hkyX04A zfyc"_ 1))"_ 2)]) pTtotal3D xCpts''

pd 'aspect 1'
pd 'pdf'

... charge density at time y3=1.00p1 ...

load 'plot'

pd 'reset'

uQ0=:,:@(1 1 2 8 1 1 1 _0.5 0"_)

xCpts=:_ 0.4 1.7 30,_ 0.5p1 0.5p1 0,_ 0 2p1 30,:_ 1.50p1 1.50p1 0"_

pd 'viewpoint 0.001 0 100'
pd 'type surface'

pd (uQ0((cx0*cos@cx2)"_ 1;(cx0*sin@cx2)"_ 1;([-@(3{])"1@JEMcn yRef@([hkyX04A zfyc"_ 1))"_ 2)]) pTtotal3D xCpts''

pd 'aspect 1'
pd 'pdf'

... charge density at time y3=1.50p1 ...

load 'plot'

pd 'reset'

uQ0=:,:@(1 1 2 8 1 1 1 _0.5 0"_)

xCpts=:_ 0.4 1.7 30,_ 0.5p1 0.5p1 0,_ 0 2p1 30,:_ 2.00p1 2.00p1 0"_

pd 'viewpoint 0.001 0 100'
pd 'type surface'

pd (uQ0((cx0*cos@cx2)"_ 1;(cx0*sin@cx2)"_ 1;([-@(3{])"1@JEMcn yRef@([hkyX04A zfyc"_ 1))"_ 2)]) pTtotal3D xCpts''

pd 'aspect 1'
pd 'pdf'

... charge density at time y3=2.00p1 ...