>with(Mambo):

>DeclareObservers(W):
>DeclarePoints(W,seq(cat(S,k),k=1..5)):
>DeclareTriads(w):
>DefineObservers([W,W,w]):
>h:=sqrt(2)*R:
>DefinePoints([W,S1,w,R,R,-h],[W,S2,w,R,-R,-h],[W,S3,w,-R,-R,-h],[W,S4,w,-R,R,-h],[W,S5,NullVector()]):

>DefineObjects(seq([W,'Sphere',point=cat(S,k),radius=R1,color=cat("{1,",convert(1-(k-1)/4,string),",0}")],k=1..5)):

>GeometryOutput(main=W,parameters=[R=1,R1=1],checktree,checkargs,filename="pyramid.geo");
>MotionOutput(ode={},parameters=[R=1,R1=1],filename="pyramid.dyn");