>

restart;

>	n:=sqrt(x^2+y^2+z^2);

>	with(plots);

Warning, the name changecoords has been redefined

>

z:=1;

>	p1:=plot3d([x/n,y/n,z/n],x=-1..1,y=-1..1,scaling=constrained,color=red):

>	z:=-1;p2:=plot3d([x/n,y/n,z/n],x=-1..1,y=-1..1,scaling=constrained,color=red):

>	unassign('z');

>

x:=1;

>	p3:=plot3d([x/n,y/n,z/n],z=-1..1,y=-1..1,scaling=constrained,color=blue):

>	x:=-1;p4:=plot3d([x/n,y/n,z/n],z=-1..1,y=-1..1,scaling=constrained,color=blue):

>	unassign('x');

>

y:=1;

>	p5:=plot3d([x/n,y/n,z/n],x=-1..1,z=-1..1,scaling=constrained,color=green):

>	y:=-1;p6:=plot3d([x/n,y/n,z/n],x=-1..1,z=-1..1,scaling=constrained,color=green):

>	display(p1,p2,p3,p4,p5,p6);

>

>