(commands used to produce figure above)

>	p1:=plottools[cuboid]([-1/2,-1/2,-1],[1/2,1/2,1],color=blue): p2:=plots[fieldplot3d](<-y,x,1/2/Pi>,x=-1..1,y=-1..1,z=-1.2..1.2,color=red,arrows=THICK,grid=[9,9,9]): plots[display](p1,p2,orientation=[-60,60]);

>